We study the Cauchy problem for a damped wave equation. Our main result is the th order approximation of solutions by the solution of the corresponding parabolic equation for arbitrary in the framework. This means that our results can be applied to nonlinear damped wave equations.
In this paper we consider the Cauchy problem for a damped wave equation
where and are given functions.
Our purpose is to study the large time behavior of the solution of (1.1). More precisely we establish the approximation of the solution u by the solution of the corresponding parabolic equation
in the framework. We use the Fourier multiplier representation of the solution and apply the Fourier multiplier theorem. As a corollary, we obtain the approximation of the solution u by the Gauss kernel
and we observe the optimality of our results with respect to the order of the decay rate.
Let us describe some previous results on decay estimates and the asymptotic behavior of the solutions for the damped wave equation (1.1). In a pioneering work in this direction, Matsumura [24] obtained classical decay estimates in the base spaces. After that, it has been asserted that the damped wave equation has a diffusive structure as ([2,21] and [5]). This means that the solution of the damped wave equation is approximated by the Gauss kernel when t is large. Subsequently, many authors showed the diffusion phenomena of the solution of (1.1). The more general case is also well investigated (including nonlinear perturbations [8,10,11,13–15,18–28,30–34] and see also references therein). The previous works show that, up to the second order, the asymptotic behavior of the solution of (1.1) is same as that of solutions to (1.2). The main goal of this paper is a description of the asymptotic profiles of solutions to (1.1) for . Regarding the above results for (1.1) and (1.2), one may expect that the asymptotic behavior of the solutions of (1.1) is as same as that of (1.2) for arbitrary order. However from the third-order expansion by the Gauss kernel, the asymptotic behavior of the solution of (1.1) is different from that of the solution of (1.2). Formally speaking, this implies that the effect of in (1.1) is not negligible for the third-order expansion. In fact, Gallay and Raugel [8] considered the variable coefficient case with nonlinear perturbations for and they suggested that this property does not hold for higher-order approximation [8, Remark 4, p. 51]. Furthermore, Orive, Zuazua and Pazoto [29] considered the variable coefficient case and they proved the higher-order expansion of the arbitrary order by the Bloch wave decomposition in for all . In contrast, their results seem to be difficult to simply apply to the arbitrary-order expansion of the solution of the nonlinear problem. Thus we only consider the constant coefficient case and we obtain the sharp higher-order asymptotic expansion in the framework under the minimal regularity condition for the initial data. Especially, our results can deal with the asymptotic behavior of solution of (1.1) in , . One can easily guess that our results are useful to show the higher-order expansion of the solution to the nonlinear problem.
To state our main results, we introduce some notation. We denote by the integer satisfying , e.g. for ,
We denote the weighted space of Nth order and the corresponding weighted Sobolev space by
where we identify . For j, k and , we put
and
In particular, the norm is denoted by . We also put
for , which appear as the coefficients of the asymptotic expansion of the solution of (1.1) by the Gauss kernel depending on the initial data.
We can now formulate our main results.
Let,and. Assume that, then there exists a constant such that the unique solution of (1.1) satisfies the estimatewhere,are defined by (1.5) andis defined by (1.3).
Combining Theorem 1.1 and the weighted estimate (1.16), we have the approximation of the solution by the Gauss kernel. This implies the lower bound of the estimate (1.7).
Let,and. Assume that, then the unique solution of (1.1) is in the classand satisfies the estimateas, whereandare defined by (1.6),,are defined by (1.5) andis defined by (1.3).
Since the estimate (1.7) in Theorem 1.1 requires only minimal regularity in the right-hand side, we can apply it to obtain sharp estimates of the solutions of the nonlinear damped wave equation such as
For the nonlinear problem (1.9), the results will be given in the separated paper [31].
For the Cauchy problem (1.2), the decay estimate and the asymptotic expansion of the solution by the Gauss kernel are well known.
Letand. Assume that. Then there exists a unique solutionof (1.2) in the classsatisfying the estimatesandasfor, whereare defined byfor.
Especially, the third-order expansion of the solution of (1.1) by the Gauss kernel is given by
as for . On the other hand, the third-order expansion of the solution of (1.2) with the data ( in the estimate (1.10)) is given by
as for . Then we can observe the difference between the asymptotic behavior of the solution of (1.1) and that of (1.2).
Here we recall that the solution formula of (1.1) which is given by the following Fourier multiplier representation (see e.g. [13])
where the evolution operators and are given by
Here we denote the Fourier and Fourier inverse transform by and . We also denote the Fourier transform of g by .
The crucial part of the demonstrating the estimate (1.7) is the construction of the higher-order approximations of the operators and in the Fourier space. Namely, the approximations
for small and the Fourier multiplier theory enable us to obtain the higher-order expansion of the solutions for the linear problem (1.1) (see also the results [7,12] and [30] for the study of the asymptotic behavior of the solution to the linear problem (1.1) which has the intersection with the main topics in this paper).
The second important part is to establish the weighted estimate
where u is the solution of (1.1) and . The estimate (1.16) shows some consistency criterion. That is, if the data is in then is also in for all .
Here we remark that Kawakami and Ueda [19] have obtained the estimate (1.16) by the different procedure. We will give another proof. We also note that in [19] they proved that the asymptotic profile of the solution to the nonlinear problem (1.9) coincides with that of the solution to the corresponding nonlinear parabolic problem
up to the second-order expansion, using the estimate (1.16) in combination with the well-known estimates
for , , which are shown by [13,23,27] and [26]. By contrast, we can regard our approximation (1.15) as an improvement of the estimates (1.17). As a consequence, Theorem 1.1 and Corollary 1.2 imply the difference between the solution of (1.1) and that of the corresponding parabolic problem (1.2).
We conclude this section by giving some notation used in this paper. Let denote the Fourier transform of f defined by
with . Also, let or denote the inverse Fourier transform. For and , let be the Sobolev space:
where is the Lebesgue space as usual and . is the homogeneous Sobolev space defined by
For the notation of the function spaces, the domain is often abbreviated. In the following C denotes a positive constant. It may change from line to line.
The paper is organized as follows. Section 2 presents some preliminaries. In Section 3, we show an appropriate decomposition of the propagators and in the Fourier space. Section 4 is devoted to the proof of Theorem 1.1, which is based on the Fourier multiplier theory. In Section 5, we show the weighted estimate of the solution of (1.1). We also prove Corollary 1.2 by applying Theorem 1.1 and the weighted estimate (1.16).
Preliminaries
In this section, we prepare several lemmas for the proofs of our results.
Estimates for the evolution operators
We begin with mentioning the – estimates for the evolution operators and .
The next lemma is useful to obtain decay estimates.
(i) Let ,,and. There exists a constant (independent of t) such that
(ii) Let ,,and. There exists a constant (independent of t) such that
The estimate (2.4) is well known (see e.g. Takeda and Yoshikawa [32]). Here we only prove the estimate (2.5). Observing that the integrand is radial, we see that
Changing the integral variable and using for , we see that
which is the desired estimate and the proof is complete. □
Decomposition of the propagators in the Fourier space
In this section, we use the Fourier splitting method to decompose the evolution operators and into low, middle and high frequency parts. In what follows, we introduce the cut-off functions which will be used in the proofs to aligned to each frequency parts, respectively. Let , and be radial cut-off functions satisfying that
Then we have the following decomposition formulas of and .
Here we only show the proof of Lemma 3.1. Lemma 3.2 is proved by the same argument as in the proof of Lemma 3.1.
At first, the direct calculation yields that
Therefore using , we see that
On the other hand, we can rewrite the second term in the left-hand side of (3.1) as
Thus, taking a difference between (3.5) and (3.6), we obtain (3.1) and Lemma 3.1 follows. □
Proof of main theorem
Our purpose in this section is to prove Theorem 1.1. According to the definitions of and (1.13) and (1.14), it is sufficient to show the following proposition.
Let,,and. Then there exists a constant such that
Once we obtain Proposition 4.1, we immediately have Theorem 1.1.
We apply Proposition 4.1 to the representation (1.12) and we obtain the estimate (1.7). The proof of Theorem 1.1 is complete. □
When , the linear estimates (2.1) and (2.2) show that
Then we have (4.1) and (4.2) for .
In what follows, we assume . Here we remark that the norms of the latter terms in (3.1) decay exponentially. Namely it is well known that
for , (cf. [13,23,26,27]). In the same manner, we can see that in the right-hand side of (3.3), the norms of the latter three terms decay exponentially (cf. [13,23,26,27]). Namely, it is also well known that
for , .
By this observation, the proof of Proposition 4.1 is reduced to the following lemma.
Let,,and. Then there exists a constant such thatwhereandare defined by (3.2) and (3.4) respectively,andare defined by (1.5).
For our purpose, we claim that
for . Indeed, once we have the estimates (4.9)–(4.12), the estimates (4.9), (4.10) and the Riesz–Thorin interpolation theorem (cf. [3]) show the estimate (4.7). Similarly, the estimates (4.11), (4.12) and the Riesz–Thorin interpolation theorem directly yield the estimate (4.8).
At first, we show the estimate (4.9). By Taylor’s expansion, the integrand is divided into two parts:
Here we define () by
where is defined by (1.4) and satisfies .
Thus we see that there exist some such that
and then we can rewrite () as
Then, for , it follows from the Hausdorff–Young inequality and (2.4) that
where we have used the fact that there exists a constant (independent of ξ) such that
To show the estimate of , we also remark that there exists a constant (independent of ξ) such that
Applying the Hausdorff–Young inequality, (2.4) and (4.18), we see that
Combining the estimates (4.16) and (4.19), we get the estimate (4.9).
Secondly, we prove the estimate (4.11). To do so, for , we let
Then the integrand ( is defined by (3.4)) is divided into three parts:
where we use (4.14) and (4.15). By the estimate and the similar fashion of the derivation of the estimates (4.16) and (4.19), we reached at the estimate
The third term is also estimated by the similar way. Indeed, we again use the Taylor expansion of ψ up to th order and
for some . Then we apply the Hausdorff–Young inequality and the estimate (2.4) to have
where we have used the estimate
for , which is independent of ξ. Combining the estimates (4.21) and (4.22), we obtain the estimate (4.11).
Next we show the estimates (4.10) and (4.12) for . In what follows, we denote Fourier multipliers in () (defined by (4.13)), () (defined by (4.20)) by
Namely we rewrite () and () as
Our strategy for the proof of (4.9) and (4.12) is the use of the Carleson–Beurling inequality (Lemma 2.4). Then we need to calculate the norm, norm (for ) and norm (for ) of multipliers. That is, we claim that
for , where we identify . We note that once we have the claim for and (4.25) and (4.26), the estimates for , and (4.27), (4.28) and (4.29) immediately follow since satisfies that there exists a constant (independent of ξ) such that
on the support of . Thus we concentrate on the proof of estimates (4.25) and (4.26).
For the proof of (4.25), when , the estimate follows by the same methods as in the proofs of (4.16) and (4.19) since we have already shown the estimates of the multipliers and in the procedure of (4.16) and (4.19). For , we recall that for and when . The multipliers in and in are radial and then it suffices to show the estimates
We denote briefly by τ. Noting that and , we have
To obtain (4.31), we prepare point-wise estimates for , and . By the definition of (4.23) and the estimate (4.17), we see
For , observing that
we arrive at the estimate
where we have used the estimates
on the support of to estimate the first term and the second term in the right-hand side of (4.35).
To obtain the point-wise estimate for , we see that the direct calculation yields that
For the first term in the right-hand side of (4.38), we have
on the support of , since
Thus by the estimate (4.39) and the similar method as in (4.37), we have
on the support of . Combining (4.33) and the estimates (4.34), (4.36) and (4.40), we obtain
on the support of . Therefore, by the estimate (2.4), it is easy to verify that
where we remark that the supports of and do not include the neighborhood of the origin . The estimate (4.25) for is easily shown by the interpolation inequality
Then we have proved the estimate (4.25) for .
For the proof of (4.32), by the similar method as in (4.31), we also have
where is defined by (4.23). In addition, the direct calculation yields that
Then as in the proof of the estimates (4.34), (4.36) and (4.40), we can assert that
Combining these inequalities gives
Therefore we see that
where we have again used the fact that the supports of and do not include the neighborhood of the origin . The estimate (4.26) for is derived by the interpolation inequality (4.41) again. Thus we have proved (4.26). Due to the previous observation and (4.30), we also have the estimates (4.27)–(4.29). Here we recall the definition of () ((3.2) and (3.4)) and the representation formulas (4.24):
Then by the estimates (4.25) and (4.26), we have
Similarly, we use the estimates (4.27)–(4.29) to obtain
Then the estimates (4.42), (4.43) and Lemma 2.4 show that
Therefore we have
by Lemma 2.3. This means the estimates (4.10) and (4.12) for , and the proof of Lemma 4.2 is complete. □
By Lemma 4.2, the estimates (4.6), (4.5), (4.3) and (4.4), we have proved the estimates (4.1) and (4.2), and the proof of Proposition 4.1 is complete. □
Weighted estimates
When , the representation formula of is well known (see Courant and Hilbert [6])
where is the modified Bessel function of order ν defined by
By changing the variable , we see that
For , such decomposition firstly pointed out by Nishihara [27].
Here, we prepare the point-wise estimate for the modified Bessel function for the proof of .
There exists a constantsuch thatfor.
The proof of Lemma 5.1 is immediately seen the representation formulas of the modified Bessel functions. See [1].
Now, we are in a position to show the weighted estimates for the evolution operators and .
Letand. Then there exists such that
The estimates (5.3) and (5.4) for are obtained in [19] by another method.
When , the proof is easy. So we omit the proof here. We only consider the case . At first, we show the estimate (5.3). For , we recall the representation formula of (5.1). Then we use the estimate (5.2) and
to have
Therefore the Fubini theorem and the estimate (2.5) show that
Thus, we have the estimate (5.3) for .
When , by (5.1) and (5.5), we see that
Therefore the Fubini theorem and (2.5) again yield
where we have used the fact that
for . Namely, we proved the estimate (5.3) for .
For the proof of , we recall the decomposition of (5.1). According to [27], we define and by
By the estimate (5.2), we see that
Again, we apply the Fubini theorem and (2.5) to have
For the term , the Fubini theorem yields
Then, combining the estimates (5.8) and (5.9) gives
and we have proved the estimate (5.3) for . Thus we have (5.3) for . Finally we prove (5.4). For this purpose we note that
When , the weight estimate (5.4) is shown by the same method as in (5.3). Indeed, for , Marcati and Nishihara [23] proved that
where . Observing that
we use the estimates (5.3), (5.2) and (5.5) to see that
Therefore we apply the Fubini theorem, the estimates (2.5) and (5.3) to have
For , applying the estimates (5.3) and (5.11) directly yields the desired estimate (5.4) for .
When , we have
by (5.1). Therefore, as in the proof of (5.11), we arrive at the estimate
For , the estimate (5.12) means the desired estimate (5.4).
To show the case , we observe that
by (5.1), (5.6) and (5.7). For the second term in the right-hand side, we apply the well-known formula (see Nishihara [27]) to see that
and
Therefore we obtain the point-wise estimate for
In (5.13), we only need to estimate the second term in the right-hand side
Indeed, by the estimates (5.3) and (5.9), we have
As in the proof of (2.5), we see that
for and . Then we obtain the estimate:
For , the combination of the estimates (5.8) and (5.14) implies the estimate (5.4) for . Thus, we complete the proof of the estimate (5.4) for . □
Summing up Theorem 1.1 and Proposition 5.2, we have the approximation of the solution to (1.1) by the Gauss kernel.
By Lemma 2.1 and Proposition 5.2, we easily see that the solution is in the class . Moreover, the estimates (1.7) and (2.3) yield that
as . Then we have the estimate (1.8). The proof of Corollary 1.2 is complete. □
Footnotes
Acknowledgements
The author would like to express sincere gratitude to Professor Takayoshi Ogawa for valuable discussion. The author is grateful to Professor Shuji Yoshikawa for his many encouragement and helpful comment. The author is also grateful to Professor Tatsuki Kawakami for his valuable conversation. The author would like to thank the referee for several helpful suggestions.
References
1.
M.Abramowitz and I.A.Stegun (eds), Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables, Dover Publications Inc., New York, 1992. (Reprint of the 1972 edition.)
2.
H.Bellout and A.Friedman, Blow-up estimates for a nonlinear hyperbolic heat equation, SIAM J. Math. Anal.20 (1989), 354–366.
3.
J.Bergh and J.Löfström, Interpolation Spaces, Grundlehren der mathematischen Wissenschaften, Vol. 223, Springer, Berlin, New York, 1976.
4.
P.Brenner, V.Thomée and L.Wahlbin, Besov Spaces and Applications to Difference Methods for Initial Value Problems, Lecture Notes in Mathematics, Vol. 434, Springer, Berlin, New York, 1975.
5.
R.Chill and A.Haraux, An optimal estimate for the difference of solutions of two abstract evolution equations, J. Differential Equations193 (2003), 385–395.
6.
R.Courant and D.Hilbert, Methods of Mathematical Physics, Vol. II: Partial Differential Equations, Interscience Publishers(a division of John Wiley & Sons), New York, London, 1962.
7.
W.Dan and Y.Shibata, On a local energy decay of solutions of a dissipative wave equation, Funkcial. Ekvac.38 (1995), 545–568.
8.
T.Gallay and G.Raugel, Scaling variables and asymptotic expansions in damped wave equations, J. Differential Equations150 (1998), 42–97.
9.
M.-H.Giga, Y.Giga and J.Saal, Nonlinear Partial Differential Equations. Asymptotic Behavior of Solutions and Self-Similar Solutions, Progress in Nonlinear Differential Equations and Their Applications, Vol. 79, Birkhäuser, Boston, MA, 2010.
10.
N.Hayashi, E.Kaikina and P.I.Naumkin, Damped wave equation with super critical nonlinearities, Differential Integral Equations17 (2004), 637–652.
11.
N.Hayashi, E.Kaikina and P.I.Naumkin, Damped wave equation with a critical nonlinearity, Trans. Amer. Math. Soc.358 (2006), 1165–1185.
12.
F.Hirosawa and J.Wirth, -theory of damped wave equations with stabilisation, J. Math. Anal. Appl.343 (2008), 1022–1035.
13.
T.Hosono and T.Ogawa, Large time behavior and – estimate of 2-dimensional nonlinear damped wave equations, J. Differential Equations203 (2004), 82–118.
14.
R.Ikehata, T.Miyaoka and T.Nakatake, Decay estimates of solutions for dissipative wave equations in with lower power nonlinearities, J. Math. Soc. Japan56 (2004), 365–373.
15.
R.Ikehata and K.Tanizawa, Global existence of solutions for semilinear wave equations in with non-compactly supported initial data, Nonlinear Anal. T.M.A.61 (2005), 1189–1208.
16.
K.Ishige, M.Ishiwata and T.Kawakami, The decay of the solutions for the heat equation with a potential, Indiana Univ. Math. J.58 (2009), 2673–2707.
17.
K.Ishige and T.Kawakami, Refined asymptotic profiles for a semilinear heat equation, Math. Ann.353 (2012), 161–192.
18.
G.Karch, Selfsimilar profiles in large time asymptotics of solutions to damped wave equations, Studia Math.143 (2000), 175–197.
19.
T.Kawakami and Y.Ueda, Asymptotic profiles to the solutions for a nonlinear damped wave equations, Differential Integral Equations26 (2013), 781–814.
20.
S.Kawashima, M.Nakao and K.Ono, On the decay property of solutions to the Cauchy problem of the semilinear wave equation with a dissipative term, J. Math. Soc. Japan47 (1995), 617–653.
21.
T.-T.Li, Nonlinear heat conduction with finite speed of propagation, in: China–Japan Symposium on Reaction–Diffusion Equations and Their Applications and Computational Aspects, Shanghai, 1994, World Sci. Publ., River Edge, NJ, 1997, pp. 81–91.
22.
T.-T.Li and Y.Zhou, Breakdown of the solutions to , Discrete Contin. Dynam. Systems1 (1995), 503–520.
23.
P.Marcati and K.Nishihara, The – estimates of solutions to one-dimensional damped wave equations and their application to compressible flow through porous media, J. Differential Equations191 (2003), 445–469.
24.
A.Matsumura, On the asymptotic behavior of solutions of semilinear wave equations, Publ. Res. Inst. Sci. Kyoto Univ.12 (1976), 169–189.
25.
M.Nakao and K.Ono, Existence of global solutions to the Cauchy problem for the semilinear dissipative wave equations, Math. Z.214 (1993), 325–342.
26.
T.Narazaki, – estimates for damped wave equations and their applications to semi-linear problem, J. Math. Soc. Japan56 (2004), 585–626.
27.
K.Nishihara, – estimates of solutions to the damped wave equation in 3-dimensional space and their application, Math. Z.244 (2003), 631–649.
28.
K.Ono, Global existence and asymptotic behavior of small solutions for semilinear dissipative wave equations, Discrete Contin. Dynam. Systems9 (2003), 651–662.
29.
R.Orive, E.Zuazua and A.F.Pazoto, Asymptotic expansion for damped wave equations with periodic coefficients, Math. Methods Appl. Sci.11 (2001), 1285–1310.
30.
P.Radu, G.Todorova and B.Yordanov, Higher order energy decay rates for damped wave equations with variable coefficients, Discrete Contin. Dyn. Syst. Ser. S2 (2009), 609–629.
31.
H.Takeda, Large time behavior of solutions for a nonlinear damped wave equation, submitted for publication.
32.
H.Takeda and S.Yoshikawa, Asymptotic profiles of solutions for the isothermal Falk–Konopka system of shape memory alloys with weak damping, Asymptotic Anal.81 (2013), 331–372.
33.
G.Todorova and B.Yordanov, Critical exponent for a nonlinear wave equation with damping, J. Differential Equations174 (2001), 464–489.
34.
Q.Zhang, A blow-up result for a nonlinear wave equation with damping: The critical case, C. R. Acad. Sci. Paris333 (2001), 109–114.