Continuity of the density of a gas flow in a porous medium
Continuity of the density of a gas flow in a porous medium
复制标题
多孔介质中气流密度的连续性
DOI:
10.1090/s0002-9947-1979-0534112-2
复制
发表时间:
1979
影响因子:
1.3
通讯作者:
A. Friedman
中科院分区:
文献类型:
--
作者:
L. Caffarelli;A. Friedman
The equation of gas in a porous medium is a degenerate nonlinear parabolic equation. It is known that a unique generalized solution exists. In this paper it is proved that the generalized solution is continuous. 0. Introduction. The density u(x, t) of gas in a porous medium satisfies the equation du/dt = Aum (m > 1) (0.1) for x G R", t > 0, and an initial condition m(x, 0) = m0(x). (0.2) Here u0(x) > 0 and u(x, t) > 0. The equation (0.1) is a nonlinear parabolic equation, degenerating at the points where u = 0. The concept of a solution of (0.1), (0.2) is taken in some weak sense (to be defined precisely in §1). The purpose of this paper is to prove that m(x, t) is continuous. (0.3) This result is known for n = 1; see [9], [10], [1] and [5]. In § 1 we state this result more precisely, giving also a uniform modulus of continuity. In §§2 and 3 we establish preliminary estimates. The proof of (0.3) for t > 0 is given in §4 and, for / = 0, in §5. 1. The main results. Let m0(x) be a function defined in R" and satisfying: 0 < m0(x) < TV (TV < oo), (1.1) f (u0(x))2 dx < oo, (1.2) u0(x) is continuous in R", and uniformly Holder continuous in every compact set where m0 > 0. Received by the editors August 7, 1978. AMS (MOS) subject classifications (1970). Primary 35K55. 'This work is partially supported by National Science Foundation Grants 74 06375 A01 and MC 575-21416 AOL © 1979 American Mathematical Society 0002-9947/79/0000-035 3/$04.75 99 License or copyright restrictions may apply to redistribution; see http://www.ams.org/journal-terms-of-use 100 L. A. CAFFARELLI AND AVNER FRIEDMAN We consider the Cauchy problem du/dt = Aum in R" X (0, oo), (1.4) m(x, 0) = Mq(x) inR", (1.5) where m is a fixed number, m > 1. By a solution of (1.4), (1.5) we mean a function u(x, t) such that, for any T < oo, fT f \(u(x, t))2 + \Vxum(x, t)\2] dxdt<<x> (1.6) J0 JR»L and /oX("í " V*"m ' Vx/) dX dt + hUo(x)Áx) ** = ° (1-7) for any continuously differentiable function / with compact support in R" X [0, T). We recall [11] that under the conditions (1.1), (1.2), there exists a unique solution. Other concepts of a solution can be given which allow for a different decay condition at x = oo than in (1.2). The results of this paper are not affected by working with these other concepts of a solution. The solution u(x, t) can be obtained as a limit of solutions m,(x, t) (tj|0) of the equation (1.4) with the initial condition u(x, 0) = m0(x) + 7] in R"; (1.8) see [11]. Notice that the solution uv of (1.4), (1.8) is taken in the classical sense, uv < u^ if tj < 17', and tj < uv(x, t) < TV + 7} in R" X (0, 00). (1.9) We define the parabolic distance between two points (x1, f1), (x2, t2) by d((x\t>),(x2,t2))=\x>-x2\+\tl-t2\X/2. When we shall speak of a modulus of continuity of a function v(x, t), we shall always mean the distance between two points to be the parabolic distance. We now introduce two moduli of continuity: <oe(/-) = qiogr|"e (0<6<2/h), (1.10) <b(r) = C2-cllos'l'/2 (1.11) where C > 0, c > 0. The main result of the paper is stated in the following theorem. Theorem 1.1 (i) The solution u of (1.4), (1.5) is continuous in R" X [0, 00); (ii) For any 80 > 0, um has a modulus of continuity ae(r) in R" X [8Q, 00) for any 0 < e < 2/n, if n > 3, and a modulus of continuity <b(r) in R" X [50, 00) if « = 2. License or copyright restrictions may apply to redistribution; see http://www.ams.org/journal-terms-of-use CONTINUITY OF THE DENSITY OF A GAS FLOW 101 The proof of continuity for t > 0 and the proof of (ii) are given in §4. It will become obvious from the proof that the condition (1.3) is not required for this part of the theorem. The proof of continuity for t = 0 is given in §5. §§2 and 3 develop some estimates needed in §4. 2. Preliminary lemmas. In this section and in §3 we obtain various auxiliary results for the solution ^(x, t) of (1.4), (1.8). For simplicity we shall denote this solution by u(x, t); we also take 0 < tj < 1. All the estimates which we shall obtain, and all the constants will be independent of tj. We set M = TV + 1, so that, by (1.9), 0 < m(x, t) < M. (2.1) Lemma 2.1. The following inequalities hold: '£■>-"-7, (2-2) at m — I dum m _ ,„ _,, t-r>-r Um. (2.3) at m — I This result is due to Aronson and Benilan [4]. Since the proof is short, we briefly give it here. The function w = t(du/dt) satisfies dw/dt = mA(um~lw) + AmTM. The function z = — u/(m — 1) satisfies the same equation, and z(x, 0) < 0 = w(x, 0) (w is continuous at t = 0 if m0(x) is smooth). By comparison, then, (2.2) follows if u0(x) is smooth; for general m0, one uses approximation. The inequality (2.3) follows immediately from (2.2). Let ô0 be a fixed positive number. For any x° G R", t° > 250, we introduce the sets Cr,h{x°, t°) = {(x, t); \x-x°\<r,t°-h<t< t°}, Br(x°)={x;\x-x°\<r}, where h < 80. Denote the volume of Br(x°) by \Br(x°)\. Set C0 = m/ (m l)80 (2.4) and let h0 = h0(S0) be any positive number satisfying: e2C°*°<f, h0<80. (2.5) Lemma 2.2. For any x° G Rn, t° > 250, X>0, r > 0, 0 <h < h0(80), the following is true: if -i-— f um(x, t°h)dx>X (2.6) \Br(x°)\ Jb,^ and if h > Mynr2/X, (2.7) License or copyright restrictions may apply to redistribution; see http://www.ams.org/journal-terms-of-use 102 L. A. CAFFARELLI AND AVNER FRIEDMAN then um(x°, t°) > i X; (2.8) here yn is a positive number depending only on the dimension n. Proof. By (2.3), (2.4), dum/dt > (m/(m \)t)um > C0um if t° h0 <t<t° (since t > 280 h0 > 80). Hence Mm(x°, t°) > e-co<'-'°+A>M'"(x0, t) > e-c^">um(x0, t). (2.9) Again, by (2.3), (2.4), the function 1 w(t) =-— f um(x, t) dx |Rr(x°)|V°) satisfies <p'(t) > (m/(m l)t)tp(t) > -C0tp(0 so that <p(t) > e-c»('-'°+*V('o h)> e-c^(t0 h)> e~c^X, (2.10) where (2.6) was used in the last inequality. Suppose n > 3 and let C7(P) = p2'" r2'" ((n 2)/2)r-"(r2 p2), p = |x x°|. (2.11) Notice that G(r) = 0, G'(r) = 0 Since G'(p) < 0 if p < r, G(p) is positive in Br(x°). By Green's formula, um(x°, t)=\y( GAum dx + —i— f um(x, t) dx (2.12) JBr{x°) \Br(x°)\ JBr(X») where y„ is a positive constant depending only on n. Suppose now that the assertion (2.8) is not true. Then (2.9) gives um(x°, t) < i Xec^o. Substituting this and (2.10) into (2.12) and using the first inequality of (2.5), we obtain Ï A < i y„ f GAum dx. JBr{x°) Integrating this inequality with respect to t, t° — h < t < t°, we get Xh<y„[' [ G(p)Aum(x, t) dx dt Jt°-h JBr(x°)