arxiv:math/ v1 [math.ap] 10 Jun 2003

Transcription

1 arxiv:math/ v1 [math.ap] 10 Jun 003 STAILITY OF SOLUTIONS OF QUASILINEAR PARAOLIC EQUATIONS GIUSEPPE MARIA COCLITE AND HELGE HOLDEN Abstract. We bound the difference between solutions u and v of u t = a u+ div xf +h and v t = b v +div xg +k with initial data ϕ and ψ, respectively, by ut, vt, L p E A E t ϕ ψ ρp L R n +t a b + x f x g + f u g u + h k ρp E ηp. Here all functions a, f, and h are smooth and bounded, and may depend on u, x R n, and t. The functions a and h may in addition depend on u. Identical assumptions hold for the functions that determine the solutions v. Furthermore, E R n is assumed to be a bounded set, and ρ p and η p are fractions that depend on n and p. The diffusion coefficients a and b are assumed to be strictly positive and the initial data are smooth. 1. Introduction We show that one can bound the difference between solutions u and v of u t = at,x,u, u u+div x ft,x,u +ht,x,u, u, x R n, 0 < t < T, u0,x = ϕx, x R n, 1.1 and v t = bt,x,v, v v +div x gt,x,v +kt,x,v, v, x R n, 0 < t < T, v0,x = ψx, x R n, 1. respectively. The assumptions are that the diffusion coefficients a and b are bounded from below by a strictly positive constant. All functions a, f, h, etc, as well as the initial data ϕ, etc, are assumed to be smooth and bounded. We are interested in estimating the local L p -norm of ut, vt, over any bounded subset E R n in terms of norm differences of the initial data as well as a and b, etc. In the hyperbolic case, that is, a = b = 0, the classical result of Kuznetsov [1] and Lucier [14] see also [9, Ch. ] reads ut, vt, L1 R ϕ ψ L 1 R +t min{t.v.ϕ,t.v.ψ} f g Lip in the one-dimensional case n = 1 where f = fu, g = gu and h = k = 0. Here T.V.φ denotes the total variation of the function φ and f Lip denotes the Date: November 9, Mathematics Subject Classification. Primary: 35K15, 3530; Secondary: 35K10, 3535, Key words and phrases. Stability, quasilinear parabolic partial differential equation, diffusion. Partially supported by the emata program of the Research Council of Norway and the European network HYKE, contract HPRN-CT The first author gratefully acknowledges the support and hospitality of the Centre of Mathematics for Applications. 1

2 G. M. COCLITE AND H. HOLDEN Lipschitz semi-norm. Recently, ianchini and Colombo [] showed flux stability in the case of hyperbolic systems on the line. Indeed, they established the estimate ut, vt, L1 R Ct Df Dg C 0 Ω for solutions u and v of u t +fu x = 0, v t +gu x = 0, respectively with u t=0 = v t=0. The usual assumptions on the flux functions and the initial conditions apply, see []. The dependence in a of the solution u of the equation u t au = 0 is treated in [1], assuming only that a is nondecreasing, and thereby allowing degenerate diffusion. However, no explicit stability estimate is provided. Otto [15] studied the equation u t div x a u,u+hu = 0 with a continuous and monotone nondecreasing. Under certain assumptions he proved that u 1 t u t 1 explt u 1 0 u 0 1. y extending Kružkov s famous doubling of variables method, ouchut and Perthame [3] showed that u 1 t, u t, L 1 R n u 0 1 u 0 L 1 R n +CT.V.u0 1 tlipa when u j satisfies u t + div x f = au j with initial data u 0 j, j = 1,. Here a is assumed to be Lipschitz and nondecreasing. Closer to the approach of this paper, Cockburn and Gripenberg [6] established the estimate u 1 t, u t, L1 R n T.V.ϕ t f 1 f +4 tn a 1 a for solutions u j, j = 1, of u j,t = div x f j + a j u j, u j t=0 = ϕ. Allowing for explicit spatial dependence in the flux function, Evje, Karlsen, and Risebro [8, 11] showed stability for solutions of u j,t +div x k j xf j u = A j u, u j t=0 = u 0 j, in the sense that u 1 t, u t, L 1 R n u 0 1 u 0 L 1 R n +tc k 1 k,bv + f 1 f,lip + tc A 1 A where,bv and,lip is the sum of the sup-norm and the V-norm and the sum of sup-norm and the Lipschitz norm, respectively. Here A j is allowed to be degenerate. Karlsen and Ohlberger [10] established L 1 contractivity of solutions of u t +div x Vt,xfu = Kt,x Au+qt,x,u. Recently, Chen and Karlsen [5] established the estimate u 1 t, u t, L 1 R n u 0 1 u 0 L1 R n +tc f 1 f + t A 1 A A 1 A 1/

3 STAILITY OF PARAOLIC EQUATIONS 3 for solutions of u j,t +div x f j u j = A j u j u j with initial data u j t=0 = u 0 j. We consider here the strictly parabolic case where the diffusion constant is not allowed to decrease to zero. However, we allow full explicit spatial and temporal dependence in all parameters. In addition, we let the diffusion and source depend explicitly on the gradient of the unknown u. All parameters, including the initial data are assumed to be smooth. Existence of regular bounded solutions is secured by classical results, see [13]. The question is to obtain explicit stability estimates. Our main result reads as follows. Let u and v denote solutions of 1.1 and 1., respectively. Then ut, vt, L p E A E t ϕ ψ ρp L R n +t a b L R 0 + x f x g L R + f u g u L R + h k L R 0 ρp E η p, where { p { 1, if 1 p, p ρ p := 1 p, if < p <, η p := + 1 n, if 1 p, 1 np, if < p <, { E p/p+1/n + E 1/p, if 1 p, A E t := C 1+t p /p E 1/np + E 1/p, if < p <, { t, if 1 p, t := C t+t /p, if < p <, for any bounded connected set E R n with Lipschitz boundary. Here R 0 = [0,T] E [ K 1,K 1 ] [ K,K ] and R = [0,T] E [ K 1,K 1 ]. As a particular example we note that for solutions u and v of u t = at,x,u, u u, v t = bt,x,v, v v with initial conditions u t=0 = ϕ and v t=0 = ψ, we find ut, vt, L E C E 1/n + E 1/ ϕ ψ L R n +Ct E 1/n a b 1/ L R 0. Our proof is based on a homotopy argument, inspired by [4]. Introducing u θ,t = θa+1 θb u θ +div x θf +1 θg +θh+1 θk, u θ t=0 = θϕ+1 θψ, we see that u 0 = u and u 1 = v. Thus u θ interpolates between u for θ = 0 and v for θ = 1. The key estimate establishes that ut, vt, L p E dist L p E ut,,vt, length L p E uθ t, 1 = u θ θ t, L dθ, p E and we establish θ-independent estimates for u θ / θ. 0

4 4 G. M. COCLITE AND H. HOLDEN. Fundamental assumptions Fix T > 0. Let u = ut,x and v = vt,x be the bounded solution of the quasilinear initial value problem see [13] u t = at,x,u, u u+div x ft,x,u +ht,x,u, u, x R n, 0 < t < T, u0,x = ϕx, x R n,.1 and v t = bt,x,v, v v +div x gt,x,v +kt,x,v, v, x R n, 0 < t < T, v0,x = ψx, x R n,. respectively. Here and f = f 1,...,f n,g = g 1,...,g n : R R n R R n, n div x ft,x,u = fj t,x,u, u x j = j=1 n j=1 f j,xj + f j u u x j = x f +f u u. Observe that x f is a scalar. The divergence operator div x always acts on the spatial variables only. y q a similarly for b, h, and k we denote the gradient of a with respect to the final n variables where u usually sits. Our fundamental assumptions are H 1 the viscous coefficients a and b are of class C 3 [0,T] R n R R n such that 0 < a a,,, a <, a C3 [0,T] R n R R n k 1, 0 < b b,,, b <, b C 3 [0,T] R n R R n k 1,.3 for some positive constants a, a, b, b, k 1 ; H the convective terms f and g are of class C 3 [0,T] R n R and the source terms h and k are of class C 3 [0,T] R n R R n such that for all i, j, l {1,...,n} and any Φ {f 1,...,f n,g 1,...,g n,h,k} the following quantities Φ L, Φ L 3 Φ L,,.4 x i x i x j x i x j x l Φ u L, Φ L u, Φ x i u L, q h L, q k L, are all bounded by a positive constant k ; H 3 the initial data ϕ and ψ are of class C R n such that for a positive constant k 3. ϕ C R n, ψ C R n k 3,.5

5 STAILITY OF PARAOLIC EQUATIONS 5 Lemma.1 L -bounds on u and v. Fix T > 0. y [13] there exist positive constants K 1, K, K 3 such that u L [0,T] R n, v L [0,T] R n K 1, u L x, v L i [0,T] R n x K,.6 i [0,T] R n u L x i x, v L j [0,T] R n x i x K 3, j [0,T] R n for all i, j {1,...,n} where K 1, K, K 3 depend only on T, n, a, a, b, b, k 1, k, and k The homotopy argument Our approach is based on the following homotopy argument. Let 0 θ 1. The function u θ interpolates between the functions u and v. More precisely, denote by u θ the solution of the quasilinear initial value problem u θ,t = θat,x,u θ, u θ +1 θbt,x,u θ, u θ u θ +div x θft,x,uθ +1 θgt,x,u θ +θht,x,u θ, u θ +1 θkt,x,u θ, u θ, x R n, 0 < t < T, u θ 0,x = θϕx+1 θψx, x R n. 3.1 Clearly u 0 = v, u 1 = u. Indeed θ u θ t, is a curve joining vt, and ut,, and ut, vt, L p E dist L p E ut,,vt, lengthl p E uθ t,, 3. for each 0 t T, E R n measurable set and 1 p. Lemma 3.1 L -bounds on u θ. y [13, Theorem V 8.1], there exist positive constants K 1, K, K 3 depending only on T, n, a, a, b, b, k 1, k and k 3 such that for each 0 θ 1 and i, j {1,...,n}. u θ L [0,T] R n K 4, u θ L x K 5, 3.3 i [0,T] R n u θ L x i x K 6, j [0,T] R n Lemma 3. Smoothness of θ u θ. Assume H 1, H, and H 3. The curve is of class C 1. In particular, we infer θ [0,1] u θ t, C R n length L p E uθ t, = 1 for each 0 t T and E R n measurable set. 0 u θ θ t, L dθ, 3.4 p E

6 6 G. M. COCLITE AND H. HOLDEN Proof. Consider the map F: D C ]0,T[ R n C [0,T] R n, Fθ,ωt,x := ω t t,x θa t,x,ωt,x, ωt,x +1 θb t,x,ωt,x, ωt,x ωt,x div x θft,x,ωt,x+1 θgt,x,ωt,x θht,x,ωt,x, ωt,x+1 θkt,x,ωt,x, ωt,x, where { D := θ,ω [0,1] C ]0,T[ R n C [0,T] R n } ω0, = θϕ+1 θψ. From the definition of u θ, Observe that F is of class C 1 and F θ θ,ω = bt,x,ω, ω at,x,ω, ω ω To compute we find Fθ,u θ 0, 0 θ div x gt,x,ω ft,x,ω +kt,x,ω, ω ht,x,ω, ω. F ω θ,ω[ θ,z ] = F ε θ,ω +εz ε=0 Fθ,ω +εz = ω t +ε z t θa t,x,ω +εz, ω+ε z +1 θbt,x,ω +εz, ω+ε z ω +ε z div x θft,x,ω +εz+1 θgt,x,ω +εz θht,x,ω +εz, ω+ε z+1 θkt,x,ω +εz, ω+ε z, F z θ,ω +εz = ε t θat,x,ω +εz, ω+ε z +1 θbt,x,ω +εz, ω+ε z z t,x,ω +εz, ω+ε z θ a ω +1 θ b ω t,x,ω +εz, ω+ε z z ω +ε z θ q at,x,ω +εz, ω +ε z +1 θ q bt,x,ω +εz, ω+ε z z ω +ε z div x θ f t,x,ω +εz+1 θ g ω ω t,x,ω +εz z θ h t,x,ω +εz, ω+ε z ω +1 θ k ω t,x,ω +εz, ω+ε z z

7 STAILITY OF PARAOLIC EQUATIONS 7 θ q ht,x,ω +εz, ω+ε z +1 θ q kt,x,ω +εz, ω+ε z z. Thus F ω θ,ω[ θ,z ] = z t θat,x,ω, ω+1 θbt,x,ω, ω z θ a ω t,x,ω, ω+1 θ b ω t,x,ω, ω z ω θ q at,x,ω, ω+1 θ q bt,x,ω, ω z ω div x θ f ω t,x,ω+1 θ g ω t,x,ω z θ h ω t,x,ω, ω+1 θ k ω t,x,ω, ω z θ q ht,x,ω, ω+1 θ q kt,x,ω, ω z, θ,ω,θ,z D. Observe that θ,z D satisfies the equation F ω θ,ω[ θ,z ] = ζ if and only if z is solution of the linear initial value problem z t = θat,x,ω, ω+1 θbt,x,ω, ω z + θa ω t,x,ω, ω+1 θb ω t,x,ω, ω ωz + θ q at,x,ω, ω+1 θ q bt,x,ω, ω z ω θfω +div x t,x,ω+1 θg ω t,x,ω z + θh ω t,x,ω, ω+1 θk ω t,x,ω, ω z + θ q ht,x,ω, ω+1 θ q kt,x,ω, ω z +ζt,x, x R n, 0 < t < T, z0,x = θ ϕx+1 θ ψx, x R n. Since this problem is well-posed see [13, Theorem IV 5.1], F ω θ,ω is invertible. y the implicit function theorem, the curve θ u θ is of class C 1 and clearly 3.4 holds. This concludes the proof. Differentiating equation 3.1 with respect to θ, we have u θ t θ = θat,x,u θ, u θ +1 θbt,x,u θ, u θ uθ θ + θ a u t,x,u θ, u θ +1 θ b u t,x,u u θ θ, u θ u θ θ uθ +θ q at,x,u θ, u θ +1 θ q bt,x,u θ, u θ u θ θ + at,x,u θ, u θ bt,x,u θ, u θ u θ 3.6

8 8 G. M. COCLITE AND H. HOLDEN Denoting f + u t,x,u θ g u t,x,u θ u θ + x ft,x,u θ x gt,x,u θ + θ f u t,x,u θ+1 θ g u t,x,u uθ θ θ + θ x f u t,x,u θ+1 θ x g u t,x,u uθ θ θ + θ f u t,x,u θ+1 θ g u θ u t,x,u θ u θ θ +ht,x,u θ, u θ kt,x,u θ, u θ + θ h u t,x,u θ, u θ +1 θ k u t,x,u uθ θ, u θ θ +θ q ht,x,u θ, u θ +1 θ q kt,x,u θ, u θ z θ t,x := u θ θ, αθ,t,x := θa+1 θb, uθ θ βθ,t,x := θ q a+1 θ q b u θ +θ f u +1 θ g u + θ q h+1 θ q k, γθ,t,x := θ a u +1 θ b u θ u +θ x f u +1 θ x g u +θ h u +1 θ k u + θ f g u +1 θ u u θ, σθ,t,x := a b f u θ + u g u θ u + x f x g +h k, for each 0 θ 1, t 0, x R n, there results z θ t = αθ,t,x z θ +βθ,t,x z θ +γθ,t,xz θ +σθ,t,x, Moreover, observe that. 0 θ 1, 0 < t < T, x R n. 3.7 z θ 0,x = ϕx ψx, 0 θ 1, x R n. 3.8 Lemma 3.3 L -bounds on α, β, γ. From the definition of α,.3 and 3.3, we have where 0 < α α,, α, α L k 1 1+K 5 +nk α := min{a,b }, α := max{a,b }.

9 STAILITY OF PARAOLIC EQUATIONS 9 Moreover, from the definition of β and.4, we infer where β L = sup β j L K 7, 3.10 j=1,...,n K 7 := nk 1 K 6 +k. Finally, from the definition of γ,.3,.4 and 3.3, we find γ L K 8, 3.11 where K 8 := nk 1 K 6 +n+1+nk 5 k. Lemma 3.4 L -bounds on z θ. Assume H 1, H, and H 3. There exists a positive constant C 1 depending only on T, n, a, a, b, b, k 1, k, and k 3 such that for each 0 t T and 0 θ 1. z θ t, L R n C 1 t+ ϕ ψ L R n, 3.1 Proof. To simplify the notation we let w denote the solution of 3.7, that is, Linearity implies that w t = α w+β w+γw+σ, w t=0 = w where w 1 and w solve w = w 1 +w w 1,t = α w 1 +β w 1 +γw 1, w 1 t=0 = w 0, w,t = α w +β w +γw +σ, w t=0 = 0, respectively. We infer from [13, p. 389] that w 1 t,x = Gt,0,x,ξw 0 ξdξ, R n t w t,x = Gt,τ,x,ξστ,ξdξdτ, 0 R n where G is the Green s function. For t [0,T] for some fixed T positive we find w t,x Ct σ. Introduce z = w 1 w 0 which satisfies the equation for w with σ = α w 0 + βdiv x w 0 +γw 0. Thus w 1 t,x zt,x + w 0 x w 0 +Ct α w 0 +β w 0 +γw 0. Observe that in the previous lemma, the smoothness of the initial condition enters in a crucial way. With less regularity we get the familiar Ot 1/ behavior near t = 0 see, e.g., [9, Sec. 4.4].

10 10 G. M. COCLITE AND H. HOLDEN 4. Stability of quasilinear parabolic equations We begin with the following lemma. Lemma 4.1 Poincaré-type inequality. There exists a positive constant Λ 0, depending only on n, such that f dx Λ 0 /n f dx+λ 0 1/n f dx, 4.1 for each f C R n and R n bounded connected set with Lipschitz boundary. In the case n = 1 we mean f dx = fx 0, for some x 0. The proof of this lemma is more or less classical see [16, Theorem A.9] and [7, Lemma A.] and the dependence of the coefficients on the measure of the domain is consequence of a standard rescaling argument. Now we prove the key estimate in the L -norm for the map z θ. Lemma 4. Case p = : Energy estimate. Assume H 1, H, and H 3. Then there exists a positive constant C depending only on T,n,a,a,b,b,k 1,k, and k 3 such that z θ t, L E C E 1/n + E 1/ ϕ ψ L R n +C t E 1/n a b L R 0 + x f x g L R 4. + f u g u L R + h k L R 0 1/, for each 0 t T, 0 θ 1 and E R n bounded connected set with Lipschitz boundary. Here R 0 = [0,T] E [ K 1,K 1 ] [ K,K ] and R = [0,T] E [ K 1,K 1 ]. Proof. Let R n be a ball and 0 < t < T. Then by 3.7 we find d 1 dt z θt,xdx = z θ z θ,t dx 4.3 = αz θ z θ dx+ z θ β z θ dx+ γzθ dx+ σz θ dx. Observe that, by 3.11, and, by 3.10, y Lemma 3.4 and.5, σz θ dx γz θ dx K 8 z θ dx, 4.4 z θ β z θ dx 1 β zθ α dx+ α z θ dx K 7 zθ dx+ α z θ dx. α 4 σ z θ dx = σ σ zθ dx

11 where STAILITY OF PARAOLIC EQUATIONS 11 1 σ dx+ 1 1 σ L R +C 1 t σ L R σ zθ dx ϕ ψ L R n σ L R K 9 σ L R +C 1 t σ L R, K 9 := 1+8k 3. Moreover, by the divergence theorem we have αz θ z θ dx = αz θ z θ νdx α z θ z θ dx α z θ dx αz θ z θ νdx+ 1 α zθ dx 4.7 α + α z θ dx α z θ dx = αz θ z θ νdx+ 1 α zθ α dx α z θ dx αz θ z θ νdx+ α L zθ dx α z θ dx, α whereν isthe externalnormalto andin the casen = 1, = {x 1,x },x 1 < x, we mean αz θ z θ νdx = αθ,t,x z θ t,x z θ,x t,x αθ,t,x 1 z θ t,x 1 z θ,x t,x 1. Substituting 4.4, 4.5, 4.6, 4.7 in 4.3 we obtain d 1 dt z θt,xdx α z θ dx+ αz θ z θ νdx 4 + K 8 + K 7 + α L zθ dx 4.8 α α y Lemma 4.1 and the assumptions on, z θ 1 dx Λ 0 /n so by Lemma 3.4, 4.8 and 4.9, d 1 dt z θ t,xdx + K 9 σ L R +C 1 t σ L R. K 8 + K 7 + α L α α zθ dx+ 1 z 1/n θ dx, 4.9 α z 4Λ 0 /n θ dx + K 9 σ L R +C1 t σ L R + αz θ z θ νdx+ α z 4 1/n θ dx 4.10

12 1 G. M. COCLITE AND H. HOLDEN K 8 + K 7 + α L α z α α 4Λ 0 /n θ dx + K 9 σ L R +C1 t σ L R + αz θ z θ νdx+ α C1 t α ϕ ψ L + R n, 1/n 1/n for some constant α > 0 assuming that, say, e.g., 1. We will eventually choose < δ < 1 sufficiently small maybe dependent on σ L R and Λ sufficiently large independent of σ L R so that Furthermore, C1 t σ L R + α C1 t Λt 1/n 1/n σ L R K 9 σ L R + αz θ z θ νdx Λ σ L R. There exists ω > 0 independent of σ L R such that α 4Λ 0 K /n 8 K 7 α L α α Substituting 4.11 and 4.1 in 4.10, we have d zθ dt t,xdx ω z /n θ t,xdx+λ σ L R ω /n Λt 1/n σ L R + α 1/n ϕ ψ L R n y the Gronwall inequality and 3.8, we have zθt,xdx exp ωt zθ0,xdx /n +Λexp ωt t exp ωτ /n 0 /n σ L Rdτ +exp ωt /n +exp ωt /n exp ωt /n t 0 t +Λ 1+/n σ L ω R + Λt ω σ L R 1/n exp ωτ /n Λτ 1/n σ L Rdτ exp ωτ α 0 /n 1/n ϕ ψ L R n dτ dx ϕx ψx exp ωt /n 1 exp ωt /n

13 Observe that, /n ω STAILITY OF PARAOLIC EQUATIONS 13 + α ω ϕ ψ L R n 1/n 1 exp ωt. /n 1 exp ωt /n and, by.3,.4 and Remark 3.1, exp ωt /n 1, t 0, t, 1 exp ωt /n 1, σ L R K 0 a b L R + x f x g L R + f u g u L R + h k L R, 0 t T, for some positive constant K 0, then, from 4.14 and since < 1, zθ t,xdx dx ϕx ψx +K 0 Λ /n t a b L R + x f x g L R ω α 1/n + f u g u L R + h k L R + ϕ ψ L ω R n + α ω 1/n ϕ ψ L R n +K 0 Λ /n t a b L ω R + x f x g L R + f u g u L R + h k L R Let now Ẽ E be a connected set such that interior of Ẽ contains the closure of E, dist Ẽ, E > 0, and Ẽ = E. Since the closure of E is compact, we can cover it with finitely many balls 1,..., m R n, that is, E j j. We may choose the balls such that j j is contained in the interior of Ẽ, and thus j j k j Ẽ = E. j=1 We assume that both j 1 and j δ < 1. Thus the result 4.15 holds and we may sum the inequality over all balls 1,..., m R n, which yields zθ E + E t,xdx C 1/n ϕ ψ L R n E which proves 4.. +C E 1/n t a b L R + x f x g L R + f u g u L R + h k L R 4.16 This proves the following result. Theorem 4.3. Fix T > 0. Let u = ut,x and v = vt,x be the classical solution of.1 and., respectively, with a = at,x,y,q and b = bt,x,y,q satisfying

14 14 G. M. COCLITE AND H. HOLDEN H 1, f = ft,x,y, g = gt,x,y, h = ht,x,y,q, and k = kt,x,y,q satisfying H, and ϕ and ψ satisfying H 3. Then there exists a positive constant C depending only on T, n, a, a, b, b, k 1, k, and k 3 such that ut, vt, L E C E 1/n + E 1/ ϕ ψ L R n +Ct a b L R 0 + x f x g L R f u g u L R + h k L R 0 1/ E 1/n, for all 0 t T with R := [0,T] E [ K 1,K 1 ], R 0 := [0,T] E [ K 1,K 1 ] [ K,K ] where E R n is bounded connected set with Lipschitz boundary. Proof. Direct consequence of 3., 3.4 and Lemmas 4., 5.1, and Estimates in L p E We want to extend the estimate of Theorem 4.3 to general p. Lemma 5.1 Case 1 p <. Assume H 1, H, and H 3. There exists a positive constant C 3 depending only on T, n, a, a, b, b, k 1, k and k 3 such that z θ t, L p E C 3 E p/p+1/n + E 1/p ϕ ψ L R n +C 3 t a b L R 0 + x f x g L R f u g u L R + h k L R 0 1/ E p/p+1/n, for each 0 t T, E R n bounded connected set with Lipschitz boundary, 0 θ 1 and 1 p <. Proof. y the Hölder inequality, z θ t, p L p E = with So, by 5., then, by Lemma 4., E z p θ t,xdx E 1/q z p θ t, L q E 5. 1/q = E 1/q z pq θ t,xdx, E q := p, q := p. z θ t, p L p E E p/ E p/ zθt,xdx E p/ z θ t, p L E, z θ t, Lp E C E p/p+1/n + E 1/p ϕ ψ L R n +C 3 t a b L R + f g L R

15 STAILITY OF PARAOLIC EQUATIONS 15 + f u g u L R + h k L R 1/ E p/p+1/n, This concludes the proof. Lemma 5. Case p >. Assume H 1, H, and H 3. There exists a positive constant C 4 depending only on T, n, a, a, b, b, k 1, k and k 3 such that z θ t, L p E C 4 1+t p /p E 1/np + E 1/p ϕ ψ /p L R n +C 4 t+t /p a b L R + x f x g L R f u g u L R + h k L R 1/p E 1/np, for each 0 t T, E R n bounded connected set with Lipschitz boundary, 0 θ 1 and < p <. Proof. Observe that z θ t, p L p E = E z p θ t,xdx z θ t, p L R n E z θt,xdx = z θ t, p L R n z θt, L E. Since /p, p /p < 1, by Lemmas 3.4 and 4., we have z θ t, L p E z θ t, p /p L R n z θt, /p L E C p /p 1 t p /p + ϕ ψ p /p L R n [ C /p E 1/np + E 1/p ϕ ψ /p L R n +C /p t /p E 1/np a b L R + x f x g L R 1/p ] + f u g u L R + h k L R C p /p 1 t p /p +k 4 [ C /p E 1/np + E 1/p ϕ ψ /p L R n +C /p t /p E 1/np a b L R + x f x g L R 1/p ] + f u g u L R + h k L R = C p /p 1 t p /p +k 4 C /p E 1/np + E 1/p C p /p 1 t+k 4 t /p E 1/np ϕ ψ /p L R n +C/p a b L R + x f x g L R + f u g u L R + h k L R 1/p,

16 16 G. M. COCLITE AND H. HOLDEN where k 4 is a positive constant such that k 3 p /p k 4, < p <. Since the maps < p < C p /p 1, C /p are bounded the proof is done. The following theorem summarizes the result in Theorem 4.3 with the extension to general p. Theorem 5.3. Fix T > 0. Let u = ut,x and v = vt,x be the classical solution of.1 and., respectively, with a = at,x,y,q and b = bt,x,y,q satisfying H 1, f = ft,x,y, g = gt,x,y, h = ht,x,y,q, and k = kt,x,y,q satisfying H, and ϕ and ψ satisfying H 3. Then there exists a positive constant C depending only on T, n, a, a, b, b, k 1, k, and k 3 such that ut, vt, Lp E A E t ϕ ψ ρp L R n +t a b L R + x f x g L R f u g u L R + h k L R ρp E η p, with R := [0,T] E [ K 1,K 1 ], R 0 := [0,T] E [ K 1,K 1 ] [ K,K ]. Here { { 1, if 1 p, p p ρ p := 1 p, if < p <, η p := + 1 n, if 1 p,, if < p <, 1 np { E p/p+1/n + E 1/p, if 1 p, A E t := C 1+t p /p E 1/np + E 1/p, if < p <, { t, if 1 p, t := C t+t /p, if < p <, for all 0 t T, where E R n is bounded connected set with Lipschitz boundary and 1 p <. Proof. Direct consequence of 3., 3.4 and Lemmas 4., 5.1, 5.. Acknowledgments. The authors would like to thank Prof. Alberto ressan many useful discussions. References [1] P. enilan and M. G. Crandall. The continuous dependence on φ of solutions of u t φu = 0. Indiana Univ. Math. J [] S. ianchini and R. M. Colombo. On the stability of the standard Riemann semigroup. Proc. Am. Math. Soc [3] F. ouchut and P. Perthame. Kružkov s estimates for scalar conservation laws revisited. Trans. Amer. Math. Soc [4] A. ressan. Contractive metrics for nonlinear hyperbolic systems. Indiana Univ. Math. J [5] G. Q. Chen and K. H. Karlsen. L 1 -framework for continuous dependence and error estimates for quasilinear anisotropic degenerate parabolic equations. Preprint, 003. [6]. Cockburn and G. Gripenberg. Continuous dependence on the nonlinearities of solutions of degenerate parabolic equations. J. Differential Equations

17 STAILITY OF PARAOLIC EQUATIONS 17 [7] G. M. Coclite. An interior estimate for a nonlinear parabolic equation. J. Math. Anal. Appl., to appear. [8] S. Evje, K. H. Karlsen, and N. H. Risebro. A continuous dependence result for nonlinear degenerate parabolic equations with spatially dependent flux function. In: Hyperbolic problems: theory, numerics, applications, Vol. I Magdeburg, 000; Internat. Ser. Numer. Math. 140, irkhäuser, asel, 001, [9] H. Holden and N. H. Risebro. Front Tracking for Hyperbolic Conservation Laws. Springer, New York, 00. [10] K. H. Karlsen and M. Ohlberger. A note on the uniqueness of entropy solutions of nonlinear degenerate parabolic equations. J. Math. Anal. Appl [11] K. H. Karlsen and N. H. Risebro. On the uniqueness and stability of entropy solutions of nonlinear degenerate parabolic equations with rough coefficients. Discrete Contin. Dyn. Syst [1] N. N. Kuznetsov. Accuracy of some approximative methods for computing the weak solutions of a first-order quasi-linear equation. USSR Comput. Math. and Math. Phys. Dokl [13] O. A. Ladyzenskaja, V. A. Solonnikov, and N. N. Ural ceva. Linear and Quasilinear Equations of Parabolic Type. Translations of Mathematical Monographs, vol. 3, American Mathematical Society, Providence, [14]. J. Lucier. A moving mesh numerical method for hyperbolic conservation laws. Math. Comp [15] F. Otto. L 1 -contraction and uniqueness for quasilinear elliptic-parabolic equations. J. Differential Equations [16] M. Struwe. Variational Methods. Applications to Nonlinear Partial Differential Equations and Hamiltonian Systems, Springer, erlin, third edition, 000. Coclite SISSA-ISAS, via eirut -4, I Trieste, Italy address: coclite@sissa.it Holden Department of Mathematical Sciences, Norwegian University of Science and Technology, Alfred Getz vei 1, NO 7491 Trondheim, Norway, and Centre of Mathematics for Applications, P.O. ox 1053, lindern, N 0316 Oslo, Norway address: holden@math.ntnu.no URL: