ASYMPTOTIC LIMIT AND DECAY ESTIMATES FOR A CLASS OF
DISSIPATIVE LINEAR HYPERBOLIC SYSTEMS
IN SEVERAL DIMENSIONS
August 4, 2019Thinh Tien NGUYEN111Department of Mathematics, Gran Sasso Science Institute – Istituto Nazionale di Fisica Nucleare, Viale Francesco Crispi, 7 - 67100 L’Aquila (ITALY), firstname.lastname@example.org
Abstract. In this paper, we study the large-time behavior of solutions to a class of partially dissipative linear hyperbolic systems with applications in velocity-jump processes in several dimensions. Given integers , let be a matrix-vector, where , and let be not required to be symmetric but have one single eigenvalue zero, we consider the Cauchy problem for linear systems having the form
Under appropriate assumptions, we show that the solution is decomposed into , where has the asymptotic profile which is the solution, denoted by , of a parabolic equation and decays at the rate as in any -norm, and decays exponentially in -norm, provided for . Moreover, decays at the optimal rate as if the system satisfies a symmetry property. The main proofs are based on asymptotic expansions of the solution in the frequency space and the Fourier analysis.
Keywords: Large-time behavior, Dissipative linear hyperbolic systems, Asymptotic expansions.
MSC2010: 35L45, 35C20.
Consider the Cauchy problem for partially dissipative linear hyperbolic systems
where and , not required to be symmetric. The system (1.1) can be regarded as discrete-velocity models where determines the velocities of moving particles and gives the transition rates of the velocities after collisions among the particles in the system. For instance, this type of dissipative linear systems arises in the Goldstein–Kac model [6, 8] and the model of neurofilament transport in axons .
The large-time behavior of the solution to (1.1) in terms of decay estimates has been established for years. It follows from  that under appropriate assumptions on and , if is the solution to (1.1) with the initial data , then one has
for some positive constants and . Moreover, the estimate (1.2) was generalized in , where can be written in the conservative-dissipative form with , a positive definite matrix not required to be symmetric. The authors in  also showed that the conservative part , describing the low dynamics of to the equilibrium manifold , satisfies the estimate
while the dissipative part of decays exponentially in , where is the solution to a parabolic system given by applying the Chapman–Enskog expansion method. A more exact asymptotic parabolic-limit of was established for a class of the generalized Goldstein–Kac system in , where the matrix is symmetric. The solution to (1.1) then satisfies the estimate
where is obtained based on an exhaustive analysis of the dispersion relation and on the application of a variant of the Kirchoff’s matrix tree theorem from the graph theory. Recently,  has optimized the above decay estimates for the solution to (1.1) in one dimension , with an explicit parabolic-limit and a corrector , namely, one has
where , solving a parabolic system arising in the low-frequency analysis decays diffusively, and solving a hyperbolic system arising in the high-frequency analysis decays exponentially. The decay estimate (1.5) is remarkable since it holds for general and ranging over . Such kind of decay estimates is very well-known e.g. the - decay estimate for the linear damped wave equation as in [11, 15, 7, 14].
To obtain (1.5) in one dimension , one primarily considers the asymptotic expansions of the fundamental solution to the system (1.1) in the Fourier space, divided into the low frequency, the intermediate frequency and the high frequency which naturally produce the time-asymptotic profile. Then, by an interpolation argument once the - estimate and the - estimate for are accomplished, one obtains the desired - estimate for any . The same strategy will be applied to the system (1.1) in several dimensions in this paper. Nevertheless, difficulties occur as the dimension increases. For instance, as mentioned in , one cannot expect the estimate
hold in general since for large time, , where is the left eigenvector associated with the eigenvalue of , behaves as the solution to the reduced system
where is the right eigenvector associated with the eigenvalue of , and thus, it is known in  that (1.6) is not true in general. The estimate (1.6) in fact depends strongly on a uniform parabolic operator. Nonetheless, this obstacle can be defeated if as in  or if is a simple eigenvalue of since the system (1.7) then becomes scalar and it allows us to obtain (1.6) as we will see in this paper. Another difficulty arises in the high-frequency analysis due to the loss of integrability and the fact that one cannot perform a uniform expansion of the fundamental solution as the dimension increases. Hence, the corrector as in (1.5) cannot be obtained trivially.
The aim of this paper is to study the - decay estimate for the conservative part of the solution to the system (1.1) in several dimensions for general and in in order to generalize (1.2), (1.3) and (1.4), where is not required to be symmetric but has one single eigenvalue zero. The - estimate as in (1.5) for the multi-dimensional case is still a challenge for the author.
For , consider the operators
where and . We start with the following reasonable assumptions. Condition . [Hyperbolicity] for is uniformly diagonalizable with real linear eigenvalues i.e. there is an invertible matrix for satisfying
for any matrix norm, such that is a diagonal matrix whose nonzero entries are real linear in .
Condition . [Diagonalizing matrix] There is a matrix uniformly diagonalizing such that is a constant matrix.
Condition . [Partial dissipation] The spectrum of is decomposed into where is simple and .
Moreover, the requisite condition for the decay of the solution to (1.1), strictly related to the Shizuta–Kawashima condition: the eigenvectors of do not belong to the kernel of for any (see [10, 17, 19] and therein), is given by Condition . [Uniform dissipation] There is a constant such that for any eigenvalue of in (1.8) for , one has
Remark 1.1 (Relaxing the conditions and ).
The requirement of the linearity of the eigenvalues of the matrix satisfying the condition and the existence of the matrix satisfying the condition can be omitted by considering the dissipative structures proposed in [18, 3]. Nonetheless, the structures in [18, 3] require that the system (1.1) is Friedrich symmetrizable while in our case, the matrix is only uniformly diagonalizable. The advantage of the linearity of the eigenvalues of the matrix and the existence of the matrix is that one can construct the high-frequency asymptotic expansion of in (1.8) after subtracting a suitable Lebesgue measure zero set.
We now construct the asymptotic parabolic-limit of the solution to (1.1). Let be an oriented closed curve in the resolvent set of such that it encloses zero except for the other eigenvalues of . One sets
which are the eigenprojection and the reduced resolvent coefficient associated with the eigenvalue zero of . We consider the Cauchy problem
where and is positive definite with scalar entries
Theorem 1.2 (- decay estimates).
Let be the solution to the Cauchy problem (1.1) with the initial data for . Under the assumptions , , and , the solution is decomposed into
and is the remainder, where is the eigenprojection associated with the eigenvalue of in (1.8) converging to as and is a cut-off function with support contained in the ball , valued in , for small . Moreover, for any and , one has
where is the solution to (1.10) with the initial data , and one has
for some constant and for all .
Remark 1.3 (Finite speed of propagation).
In the case where the solution to the system (1.1) has finite speed of propagation, since the fundamental solution associated with has compact support contained in the wave cone for some constant , one can decompose into , where
and is the remainder, where is a cut-off function with support contained in the ball , valued in , for any , and the estimates (1.13) and (1.14) still hold for . This fact will be proved in the subsequent sections. For instance, it is the case where the system (1.1) is Friedrich symmetrizable. Nonetheless, in one dimension , the case can be treated since the Cauchy integral theorem holds for the whole complex plane, and thus, one can use the estimates for the asymptotic expansion of the fundamental solution in the high frequency after changing paths of integrals of holomorphic functions (see ).
Moreover, consider the one-dimensional linear Goldstein–Kac system
It can be checked easily that satisfies the linear damped wave equation
where and are appropriate initial data. It then follows from  that
for any and , where is the solution to the heat equation
Without regarding the exponentially decaying term in (1.15), there is a difference of a quantity of between the decay rates (1.15) and (1.13). The difference can be explained by a symmetry property that the one-dimensional linear Goldstein–Kac system possesses. Such kind of symmetry properties is already studied in  based on the existence of an invertible matrix commuting with and anti-commuting with the matrix of one-dimensional dissipative linear hyperbolic systems. More general, in several dimensions , the symmetry property is given by Condition . [Symmetry] There is an invertible matrix for such that
where is given by (1.8) for .
In practice, it is easy to check the condition if there is a constant invertible matrix satisfying and for all since by definition.
We will show that under the conditions , and , the decay rate in the estimate (1.13) increases. We primarily refine the asymptotic profile .
where with matrix entries
Theorem 1.5 (Optimal decay rate).
The paper is organized as follows. Section 2 is devoted to proofs and examples of Theorem 1.2 and Theorem 1.5, where the proofs are based on the estimates obtained in Section 5. In order to prove these estimates in Section 5, we primarily invoke some useful tools of the Fourier analysis and the perturbation analysis in Section 3. With these tools, we construct the asymptotic expansions of the operator in (1.8) in Section 4 in order to obtain the asymptotic expansions of the fundamental solution to the system (1.1) to be able to prove the estimates in Section 5.
Notations and Definitions
Let be a function from to a Banach space equipped with norm , we define the Lebesgue spaces for consisting of functions satisfying
Let be the multi-index with . One denotes by
where , the partial derivatives of a smooth function on . Then, for smooths functions and on , we have the Leibniz rule
is the multi-index binomial coefficient, means that for all and the difference is defined by
The Schwartz space is the set of smooth functions on such that for any , we have
One denotes by the dual space of and is called a tempered distribution. For , the Fourier transform is defined by
where is the usual scalar product on , and the inverse Fourier transform of also denoted by is given by
On the other hand, we can define the Fourier transform of tempered distributions by the inner product on , namely
Let , the Sobolev space consists of tempered distributions such that and
Let , is called a Fourier multiplier on for if the convolution for all and if
The linear space of all such is denoted by equipped with norm .
Consider also the kernel associated with the system (1.16). One has
Proof of Theorem 1.2.
Moreover, by the relation (2.1), one has
where is given by (2.2).
On the other hand, we decompose
and is the remainder, where is the eigenprojection associated with the eigenvalue of in (1.8) converging to as and is a cut-off function with support contained in the ball , valued in , for small .
where and is a cut-off function with support contained in , valued in , for large .
Consider the three-dimensional linear Goldstein–Kac system (1.1) where
where for and , and the initial data is . Let for , we also consider the system (1.10) where and if , the matrix is given by
Moreover, the initial data is chosen as
Theorem 1.2 then implies that the solution to the three-dimensional Goldstein–Kac system can be decomposed into such that the difference decays in at the rate with respect to in as for any , where is the solution to the above system (1.10). The formulas of and in fact coincide the formulas obtained by using the graph theory as in Example 3.3 p. 412 in .
We give the proof of Theorem 1.5.
Proof of Theorem 1.5.
Consider the two-dimensional linearized isentropic Euler equations with damping
which can be written in the vectorial form
where with the initial data and
Moreover, the matrix satisfying the condition and the matrix satisfying the condition are given by
Then, Theorem 1.5 implies that , where has the asymptotic profile, which is the solution to the Cauchy problem
The proofs of Theorem 2.1 in  then implies that decays in at the rate with respect to in as for any , where is the solution to the heat equation
Remark 2.3 (Proof of the case of finite speed of propagation).
where is a cut-off function with support contained in the ball , valued in , for any . The proof is then similar to the above proofs. Moreover, this property holds for the above two examples since they are in fact symmetric hyperbolic systems.
3. Useful lemmas
3.1. Fourier analysis
We introduce here the two well-known inequalities which are the Young inequality and the complex interpolation inequality. On the other hand, we also introduce a powerful Fourier multiplier estimate which is the estimate (3.1) given by Lemma 3.3. The multiplier estimates are very helpful to study the - estimate for .
Lemma 3.1 (Young’s inequality).
For satisfying and any and , one has and the inequality
See the proof of Lemma 1.4 p. 5 in . ∎
Lemma 3.2 (Complex interpolation inequality).
Consider a linear operator which continuously maps into for with . Let be such that
then continuously maps into and one has
See the proof of Corollary 1.12 p. 12 in . ∎
Lemma 3.3 (Carlson–Beurling).