arxiv: v1 [math.pr] 28 Dec 2013

Transcription

1 One-dimensional heat equation with discontinuous conductance Zhen-Qing Chen and Mounir Zili arxiv:3.7396v [math.pr] 8 Dec 3 December 3, 3 Abstract We study a second-order parabolic equation with divergence form elliptic operator, having piecewise constant diffusion coefficients with two points of discontinuity. Such partial differential equations appear in the modelization of diffusion phenomena in medium consisting of three kind of materials. Using probabilistic methods, we present an explicit expression of the fundamental solution under certain conditions. We also derive small-time asymptotic expansion of the PDE s solutions in the general case. The obtained results are directly usable in applications. AMS Mathematics Subject Classification: Primary 6H, 6J55; Secondary 35K, 35K8 Keywords and phrases: Stochastic differential equations, semimartingale local time, strong solution, heat kernel, asymptotic expansion Introduction In many practical applications, one is encountered with heat propagations in heterogeneous media consisting of different kind of materials. Mathematically, such kind of heat propagation is modeled by heat equations with divergence form elliptic operators having discontinuous diffusion coefficients. In this paper we exam one-dimensional model where the medium consisting of three different kind of materials. Let a, a i,ρ i, i =,,3 be positive constants and define Set Ax = a {x } +a {<x a} +a 3 {a<x} and ρx = ρ {x } +ρ {<x a} +ρ 3 {a<x}. L = d ρx dx which is a self-adjoint operator in L R;ρxdx. The parabolic equation ρxax d,. dx ut, x t = Lut,x. Research partially supported by NSF Grant DMS-676 and NNSFC Grant 8.

2 describes the heat propagation in the line that consists of three different kind of material. It is wellknown that there is an m-symmetric diffusion process X associated with L, where mdx := ρxdx. Since both Ax and ρx are bounded between two positive constants, it is known see [7] that X possesses a jointly continuous transition density function qt, x, y with respect to the measure mdx on R that is symmetric inxandy and enjoys the following Aronsontype Gaussian estimates: there exists positive constants C,C so that for every t > and x,y R, C t / exp C x y /t qt,x,y C t / exp x y /C t..3 Using Dirichlet form theory, the following semimartingale representation of X is derived in [, Proposition 5]: dx t = AX t db t + ρ a ρ a d L t + ρ 3a 3 ρ a d L a ρ a +ρ a ρ 3 a 3 +ρ a t,.4 where B is one-dimensional Brownian motion, L and L a are symmetric semimartingale local times of X at and a, respectively. Here for a semimartingale Y, its symmetric semimartingale local time L w t Y at w R is defined to be L w t Y := lim ε ε [w ε,w+ε] Y s d Y s,.5 where Y is the predictable quadratic variation process of Y. We define the semimartingale local time for Y at level w R to be L w t Y = lim ε ε [w,w+ε] Y s d Y s..6 It is known see [] that equation.4 admits a unique strong solution. Moreover, X is a strong solution to.4 if and only if it is a strong solution to dx t = AX t db t + ρ a dl ρ a tx+ ρ a dl a ρ 3 a tx..7 3 The purpose of this paper is to study the parabolic equation. with initial value u,x = hx and solve it as explicitly as possible, where h is a piecewise C N+ function on R of the form hx = h x {x } +h x {<x a} +h 3 x {x>a},.8 with h i being bounded with continuously differentiable up to order N +, i =,,3, for some integer N. This kind of PDEs appears in the modelization of diffusion phenomena in many areas, for example, in geophysics [9], ecology [4], biology [] and so on. The non-smoothness of the coefficients reflects the multilayered medium. The solution ut, x of. with initial value hx admits a probabilistic representation: ut,x = E x [hx t ], t..9 Here the subscript x in E x means the expectation is taken with respect to the law of the diffusion process X that starts from x at time t =. In the first part of this paper, we show that, when either ρ a = ρ a or ρ a = ρ 3 a3, we can reduce the diffusion process of.3 to skew Brownian motion through transformations by scale functions. Thus in this case, we are able to derive an explicit formula for the transition density function qt, x, y of the diffusion process X with respect to the measure mdx. The kernel qt, x, y is the fundamental solution also called

3 heat kernel qt,x,y of L. Skew Brownian motions are a subclass of diffusion processes X of.3 with a = a = a 3 = and ρ = ρ 3. It is first introduced by Ito and McKean in 963, and has since been studied extensively by many authors; see, for example, [,, 3, 5] and the references therein. However we do not know the explicit formula for qt,x,y in the general case. In [6], formulas for the fundamental solutions are derived in terms of inverse Fourier transform and Green functions, but they are not explicit. Though we do not have explicit formula for its heat kernel, in the second part of this paper, we are able to derive asymptotic expansion in small time of solutions of. in the general case with piecewise C N+ initial data, by making use of the explicit heat kernel obtained under the special case a,ρ = a 3,ρ 3. The small-time asymptotic expansion we get in this paper is explicit and can be directly used in applications. This extends the work of [7] and of [8], where only one discontinuity is studied. In [5] and [], numerical methods are proposed to study.4 with discontinuous coefficients. The rest of the paper is organized as follows. In the section, we derive explicit formula for the fundamental solution of L in the case when either ρ a = ρ a or ρ a = ρ 3 a3. We then study the small-time asymptotic expansion of the solution of. in the general case for any positive constants a, a i andρ i. Extensions to themore general piecewise constant coefficient case with more than two discontinuities are mentioned at the end of the paper. Fundamental solution In view of the probabilistic representation.9, properties of solution to the heat equation.9 can be deduced from the properties of the diffusion process X of.7. For notational simplicity, let s rewrite SDE.7 as { dy x t = Y x = x R, p {Y x t } +q {<Y x t a} +r {a<y x t } db t + α dl t Y x + β dla t Y x, where p, q, r are positive constants, α, β, and B a one-dimensional Brownian motion. The superscript x indicates the process Y x starts from x at t =. In the particular case when p = q = r = and β =, Yt x is just a Skew Brownian motion of parameter α cf. [5], and in the case when p = q = r =, Yt x is the double-skewed Brownian motion cf. []. The next two theorems are known; see [, Theorems. and.] or [8]. We give a proof here not only for reader s convenience but also some formulas in the proof will be used later to derive the transition density function of Y. Theorem. For every α < and β <, SDE. has a unique strong solution for every x R.. Proof. Define x/p for x <, sx = x/q for x [,a], x a/r +a/q for x > a. Clearly, s is strictly increasing and one-to-one. Let σ denote the inverse of s: px for x <, σx = qx for x [,sa], rx sa+saq for x > sa...3 3

4 Let s l and σ l denote the left hand derivative of s and σ, respectively; that is, /p for x, p for x, s l x = /q for x,a], and σ l x = q for x,sa], /r for x > a, r for x > sa, The second derivative σ dx = q pδ {} dx + r qδ {sa} dx. Here δ {z} denotes the Dirac measure concentrated at point z. Since both qα p + and rβ q + are strictly less than, by Theorem. of [], the following SDE has a unique strong solution Z = {Z t, t } with Z = sx: dz t = db t + qα + dl p tz+ rβ + dl sa t Z..4 q Let Y = σz. We will show that Y is a solution to.. Clearly Y = x. Function σ is a piecewise linear function and can be expressed as a difference of two convex functions. By Ito-Tanaka s formula see [3, Theorem VI..5], Y is a semimartingale and in fact, for t, Y t = σz + L w t Zσ dw = x+ By Corollary VI..9 of Revuz and Yor [3], for t, L sa t Z = lim ε ε = r lim ε σ l Z sdz s + p {Yt } +q {<Yt a} +r {Yt>a} db t + αq L tz+ βr Lsa t Z..5 [sa,sa+ε Z s d Z s = lim ε rε t ε [a,a+rε Y s ds [a,a+rε Y s d Y s = L a ty/r..6 A similar calculation shows that L tz = L ty/q. Thus we have from.5 that Y = σz is a strong solution for. with Y = x. We now examine pathwise uniqueness. Suppose that Y x is a solution of.. Since s is the difference of two convex functions, we can apply Ito-Tanaka s formula to get sy x t = sx+ s l Y x s dy x s + R L w t Y x s dw..7 Let Z = sy x. A similar calculation as above shows that Z satisfies SDE.4 with initial value Z = sx. Since by Theorem. of [] solutions to.4 is unique, hence so is solution Y x to.. Theorem. The process Y x is a strong solution of the SDE., if and only if Y x is a strong solution to dyt x = p {Y x t } +q {<Y x t a} +r {a<y x t } db t + α α d L ty x + β β d L a ty x.8 with Y x = x. 4

5 Proof. This is due to the fact that, for solution Y x of. see the proof of Theorem. in [], L t Y x = α L t Y x, L a t Y x = β L a t Y x,.9 and L w t Y x = L w Y x for w / {,a}. Corollary.3 For every p,q >,α < and x R, the following SDE has a unique strong solution for every x R: { dx x t = p {X x t } +q {<X x t } db t + α dl t Xx, X x = x R.. Moreover, the transition probability density of the diffusion X x is given by p X {y } t,x,y = + { {y>} exp fx fy πt p q t + p+qα }. fx + fy signy exp p qα t where fy = y p {y } + y q {y>}. Proof. Taking q = r and β =, we have by Theorem. that for every p,q >, α < and x R, dxt x = p {X x t } +q {<X x t } db t + α dl tx x,. has a unique strong solution with X x = x. Let s be the function defined by.. We see from.4 that Z := sx x satisfies SDE.4 with q = r and β = ; that is, dz t = db t + qα + dl p tz..3 Thus Z is a skew Brownian, whose transition density function is explicitly known. Using Theorem., we can rewrite Z in terms of its symmetric semimartingale local time L Z: dz t = db t + p+qα p qα d L tz..4 Thus by [5] see also Exercise III..6 on page 8 of [3], Z has a transition density function p Z t,x,y with respect to the Lebesgue measure on R given by p Z t,x,y = πt exp x y The desired formula. now follows from the fact that t p X t,x,y = + p+qα signy exp p qα σ l sypz t,sx,sy x + y t. 5

6 Remark.4 iwhenα = p /q, Corollary.3inparticularrecovers theformulagiven in[, Theorem 5.] see the proof of Proposition 5. of [] for a more explicit one for the transition density function of X x in.. Note that in this case, the more compact formula. in this paper is equivalent to the form given in []. Using the connection to skew Brownian motion established in the proof of Corollary.3, one can derive the joint distribution of X x given by.8, its semimartingale respectively, symmetric semimartingale local time L X x respectively, L X x at level and its occupation time on [, from that of the corresponding skew Brownian motion. The latter can be found in [, Theorem.]. Same remark applies to the diffusion process Y x satisfying SDE.6 below. ii The diffusion process X of. can be called oscillating skew Brownian motion, or skewed oscillating Brownian motion assuggested bywellner [6]. Whenα = p /q, themarginal distribution of X x of. characterized by its transition density function p X t,x,y is the Fechner distribution in statistics; see [6]. This corresponds exactly to the case when the density function p X t,,y is continuous in y. For generator L in. or its associated SDE.7, the condition α = p /q corresponds precisely to the case of ρ = ρ 3. The next theorem allows us to obtain the transition density function for solution Y of. when β = q r. Theorem.5 Let p,q,r,a > and α <. For x R, let X x denote the unique strong solution of. with initial value X x = x. The process Y defined by: a+ r q Xx t a+ Xt x a if Y = x a Y t =.5 a+ r q X q r x a+a t a + X q r x a+a t a if Y = x > a. is the unique strong solution of the stochastic differential equation: dyt x = p {Y x t } +q {<Y x t a} +r {a<y x t } db t + α dl t Y x + q dl a t r Y x.6 with initial value Y x = x. Proof. We will present the proof just in the case where x a. The other case is similar. Consider the bijective function defined by: ϕx = a+r/qx a + x a for x R. Then Y = ϕx x. Note that ϕx = x for x a, and ϕx > a if and only if x > a. Since ϕ is the difference of two convex functions, applying the Ito-Tanaka formula to Y = ϕx x, we obtain: Y t = Y + = x+ = x+ ϕ l Xx sdxs x + R {X xs a} + r q {X xs >a} dxs x + p {X x s } +q {<X x s a} +r {X x s >a} L w t X x ϕ dx r q L a tx x db t 6

7 + α dl tx x + r q L a tx x = x+ p {Y x s } +q {<Y x s a} +r {Y x s >a} db t + α dl t Xx + r q L a t Xx..7 By the same reasoning as that for.6, we have L tx x = L ty and L a tx x = q r La ty..8 This together with.7 shows that Y is a strong solution to SDE.6 with Y = x. The uniqueness follows from Theorem.. Corollary.6 When β = q r, the solution Y to SDE. has a transition probability density p Y t,x,y with respect to the Lebesgue measure on R given by p Y t x,y = {y } + {<y a} + { {y>a} exp sx sy πt p q r t + p+qα signy exp p qα sx + sy t }.9 Proof. This follows immediately from Theorem.5 and Corollary.3. Remark.7 SDE. with β = q r corresponds exactly to SDE.7 with ρ a = ρ 3 a3. If we use qt,x,y to denote the transition density function of Y with respect to its symmetrizing measure mdx = ρxdx, then clearly we have qt,x,y = p Y t,x,y/ρy. 3 Small-time asymptotic expansion Inthis section, wewill give an explicit asymptotic expansion, in small time, of thefunction ux,t = EhYt x, in the general case where p,q,r > and α,β,. Definition 3. For every k N, we denote by erfc k the function defined on R by: erfc k z = + u k e u du, z R. π z Observe that erfc k = π Γ k+, where Γ is the Euler Gamma function. The following elementary properties are straightforward from its definition. Lemma 3. For every given integer n, 7

8 i As z +, erfc k z = oz n. Γ ii As z : erfc k z = π + k k+ +o z n. Here the notation o z n o z n as z + respectively, z means that lim z + z n = o z n respectively, lim z z n =. Similar notation will also be used for ot n as t. Note that by Corollary.3, SDE { ds x t = q {S x t a} +r {a<s x t } db t + β dla ts x S x = x R 3. has a unique strong solution for every x R and it has a transition probability density with respect to the Lebesgue measure on R given by {y a} pt,x,y = + { {y>a} exp gx gy πt q r t q +rβ gx + gy + signy a exp q rβ t }, 3. where gx = x q {x a} + x r {x>a}. The following estimation is a key in small-time asymptotic expansion for the solution of heat equation.. It reduces the case to diffusions of the form. for which we have explicit knowledge about their heat kernels. Lemma 3.3 Define X x, solution of sde. if x a Ξ x = S x, solution of sde3. if x > a 3.3. Then there exist positive constants c and c such that for every t >, P Ξ s Ys x for some s [,t] c exp c a /t, where Y x is the diffusion process of.. Proof. We will prove the lemma in the case where x a. The other case is similar. For t >, let us denote by τ t = inf{t ;Ξ s = a} = inf{t ;Xs x = a}, and A t = { } { } ω Ω : Ξ s ω Ys x ω for some s [,t] = ω Ω : Xsω x Ys x ω for some s [,t]

9 It is clear that, A t {τ t t} and consequently, PA t Pτ t t P x sup Xs x x a x P x sup Xs x x a/. s t s t Here P x is the probability law of X x. Since X is a diffusion whose transition density function qt,x,y with respect to the measure mdx = ρxdx enjoys the Aronson type Gaussian estimate.3, we have by a similar argument as that for [4, Lemma II..] that P x sup Xs x x a/ c e c a /t. s t This proves the lemma. Corollary 3.4 Let Ξ x be defined by 3.3. Then for every integer n and time t > small, [ ] ] E hyt x = E [hξ x t +ot n. Proof. Let us consider the case where x a ; the other case is similar. For all t, [ ] [ ] [ ] E hyt x = E hyt x ;A t +E hyt x ;Ac t, where A t is the set defined in 3.4, and A c is the complementary of A. Since Xs xw = Y s x w for all s [,t] on A c t, we have [ ] [ ] [ ] E hyt x = E hxt x ;A c t +E hyt x ;A t [ ] [ = E hxt x +E hyt x hxt x ;A t ]. The desired result now follows from Cauchy-Schwarz inequality, the fact that h is bounded and Lemma 3.3. Let h be a piecewise C N+ function of the form.8. By Taylor expansion, for every x R\{,a}, N h j x hy = y x j +O y x N+, j! and for x {,a} and i =,,3, h i y = j= N j= h j i x y x j +O y x N+. j! Here the notation O y x N+ means that there are constants C,δ > so that the term O y x N+ is no larger than C y x N+ for any y R with y x < δ. Thus we have the following. 9

10 Proposition 3.5 If x R\{,a}, and if x {,a}, E[hΞ x t] = N j= [ j! hj xe Ξ x t x j] +O Ξ x t x N+, 3.5 E[hΞ x t] = + N j= N j= ] j! hj [Ξ xe x t x j {Ξ x t } + N j= ] j! hj [Ξ xe x t x j {<Ξ x t a} ] j! hj 3 [Ξ xe x t xj {Ξ x t >a} +O Ξ x t x N Using the transition probability densities of the diffusions X x and S x given in corollary.3 and in 3., we can compute the expectations in Proposition 3.6 Let k be an integer. i For x <, [ E Ξ x t x k] x = k p k k/ t k/ erfc k p + t k j= k x j/ t j/ A j xerfc j j p, t where A j x = x k j p+qα p j k+ k j +pq j q p qα p k j and k = j k! j!k j!. ii For < x < a, where [ E Ξ x t x k] x = k/ q k t k/ erfc k q + t B j x = iii For x > a, [ E k j= k x j/ t j/ B j xerfc j j q, t x k j qα p j p k j p qα q +p+qα q j k j. Ξ x t xk] = k/ r k t k/ erfc k a x r t + k j= k j j/ t j/ C j xerfc j x a r t where C j x = x ak j rβ j q j q k j q rβ r +q +rβ r j k j k j. As [ a consequence of the above Proposition, and by Lemma 3., we obtain the following expansion of E Ξ x t x k] for x R\{,a} and small time t >.,

11 Corollary 3.7 For any positive integers n > k, x R\{,a}, as t +, [ E Ξ x t x k] = k/ D k xt k/ +ot n, where D k x = π + k Γ k + [ p k {x<} +q k {<x<a} +r k {x>a}]. 3.7 Using again the transition probability densities of the diffusions X x and S x, we get Proposition 3.8 For x {,a} and integer k, q α k+ [ ] E Ξ x t x k {Ξ x t } = p qα k p ktk/ Γ when x = π r β a q rβ k q k t k/ erfc k q when x = a, t [ ] E Ξ x t x k k + a {<Ξ x t a} = A k xq k t k/ π Γ erfc k q, t p a [ ] E Ξ x p qα qk k/ t k/ erfc k q when x = t t xk {Ξ x t >a} = q rkk/ t k/ k + Γ when x = a. q rβ π Here p p qα A k x = r β k q rβ when x = when x = a. 3.8 From Proposition 3.8 and Lemma 3. we deduce: Corollary 3.9 For x {,a} and integers n > k, [ ] q α E Ξ x t x k {Ξ x t } = p qα k p kk/ t k/ k+ Γ when x = π ot n when x = a; [ ] E Ξ x t x k <Ξ x t a = A k xq k t k/ k+ Γ +ot n, π [ ] ot n when x = E Ξ x t xk Ξ x t >a = q rkk/ t k/ k+ Γ when x = a, q rβ π where A k x is defined in 3.8. Combining Proposition 3.5 with Corollaries 3.7 and 3.9, we have Theorem 3. For small time and x R, N E[hΞ x t ] = b k xt k/ +Ot N+/, k=

12 where k/ D k x k hx k! x k k/ b k x = k!p qα π Γk+ k/ k!q rβ π Γk + with the function D k x given by 3.7. k h x k q α k p k + k h x k pq k k h a x k r β k q k + i h 3 a x k qr k when x R\{,a} when x = when x = a 3.9 By Theorem 3. and Corollary 3.4 we get: Corollary 3. For every x R, when t > is small, N E[hYt x ] = b k xt k/ +Ot N+/, k= where the function b k x is defined by 3.9. Remark 3. i When a = and r = q, we recover the expansion obtained in [7]. ii Employing the same approach used in this section, one can obtain without any difficulty a similar small-time expansion for the solution of the heat equation. where Ax and ρx are piecewise constant with n 3 points of discontinuity; that is, in case of the following SDE: dy x t = Y = x R, AY x t db t + n i= α i dla i t Y x, where a = < a <... < a n, p i,+, α i, for i =,,...,n and AY x t = p {Y x n t a } + i= p i {ai <Y x t a i+} +p N {an<y x t }. See [] for the existence and uniqueness of strong solution of this SDE. References [] T. Appuhamillage, V. Bokil, E. Thomann, E. Waymire and B. Wood, Occupation and local times for skew Brownian motion with applications to dispersion across an interface. Ann. Appl. Probab., Corrections: Ann. Appl. Probab., 5-5. [] R. F. Bass and Z.-Q. Chen, One-dimensional stochastic differential equations with singular and degenerate coefficients. Sankhya 67 5, [3] K. Burdzy and Z.-Q. Chen, Local time flow related to skew Brownian motion. Ann. Probab. 9,

13 [4] R. Cantrell, C. Cosner, Diffusion models for population dynamics incorporating individual behavior at boundaries: Applications to refuge design. Theor. Population Biology, 55, 999, [5] P. Étoré, On random walk simulation of one-dimonsional diffusion processes with discontinuous coefficients. Electron. J. Probab. 6, [6] B. Gaveau, M. Okada, T. Okada, Second order differential operators and Dirichlet integrals with singular coefficients. Tohoku Math. J , [7] P. Gyrya and L. Saloff-Coste. Neumann and Dirichlet Heat Kernels in Inner Uniform Domains. Paris: Socit mathmatique de France,. [8] J. G. Le Gall, One-dimensional stochastic differential equations involving the local times of the unknown process. In Stochastic Analysis and Applications Swansea, 983. Lecture Notes in Math , 58. Springer, Berlin. [9] A. Lejay, Monte Carlo methods for fissured porous media: a gridless approach. Monte Carlo Methods Appl., 4, [] A. Lejay and M. Martinez, A scheme for simulating one-dimensional diffusion processes with discontinuous coefficients. Annals Appl. Probab. 6 6, [] S. Nicas, Some results on spectral theory over networks, applied to nerve impulse transmission. Orthoginal Polynomials and Applications Bar-le-Duc, 984. Lect. notes Math. 7. Springer, [] J. M. Ramirez, Multi-skewed Brownian motion and diffusion in layered Media. Proc. Amri. Math. Soc. 39, [3] D. Revuz and M. Yor, Continuous Martingales and Brownian Motion. Springer-Verlag, Berlin, 99. [4] D. W. Stroock, Diffusion semigroups corresponding to uniformly elliptic divergence form operators. Séminaire de Probabilités, XXII, , Lecture Notes in Math., 3, Springer, Berlin, 988. [5] J. Walsh, Diffusion with discontinuous local time. Astérisques , [6] J. Wellner, Fechner s distribution and connections to skew Brownian motion. Preprint, 3. [7] M. Zili, Développement asymptotique en temps petits de la solution d une équation aux dérivées partielles de type parabolique généralisée au sens des distributions-mesures. Note des Comptes Rendues de l Académie des Sciences de Paris, t. 3, Série I, p. 49-5, 995. [8] M. Zili, Construction d une solution fondamentale d une équation aux dérivées partielles à coefficients constants par morceaux. Bull. Sci. Math , Zhen-Qing Chen Department of Mathematics, University of Washington, Seattle, WA 9895, USA zqchen@uw.edu Mounir Zili Department of Mathematics, Faculty of sciences of Monastir, Tunisia Mounir.Zili@fsm.rnu.tn 3