Global Existence of Solutions to the Timoshenko-Fourier System in Critical Besov Spaces

Abstract

In this work, we are concerned with the Timoshenko-Fourier system in both equal and non-equal wave speeds, which admits a non-symmetric dissipation. Furthermore, the dissipative mechanism of regularity-loss type will occur in the case of non-equal wave speeds. We establish the global-in-time existence of solutions to the Timoshenko-Fourier system in critical Besov spaces with the regularity s= 3 2 .

Share and Cite:

Cao, H. and Fu, X. (2025) Global Existence of Solutions to the Timoshenko-Fourier System in Critical Besov Spaces. Journal of Applied Mathematics and Physics, 13, 3112-3128. doi: 10.4236/jamp.2025.139175.

1. Introduction

We are devoted to the following coupled system of the Timoshenko theory of vibrating beams [1] [2] with additional heat conduction effects according to Fourier’s law:

{ φ tt ( φ x ψ ) x =0, ψ tt [ σ( ψ x ) ] x ( φ x ψ )+γ ψ t +β θ x =0, θ t k θ xx +β ψ tx =0, (1.1)

where t( 0,+ ) is the time variable and x is the spatial variable. The function φ and ψ denote the transverse displacement of the beam from equilibrium and the angle of turn of the beam filaments, θ is the temperature difference of the elastic material, and γ , β and k are physical constants dependent on the elastic and thermal properties of the material. For any η>0 , the smooth function σ( η ) satisfies σ ( η )>0 , we are devoted to the Cauchy problem for the system (1.1), with initial conditions of

( φ, φ t ,ψ, ψ t ,θ )| t=0 =( φ 0 , φ 1 , ψ 0 , ψ 1 , θ 0 )( x ),x.

The linearized system of (1.1) is accordingly written as

{ φ tt ( φ x ψ ) x =0, ψ tt a 2 ψ xx ( φ x ψ )+γ ψ t +β θ x =0, θ t k θ xx +β ψ tx =0, (1.2)

where a>0 represents the speed of sound defined by a 2 = σ ( 0 ) . a=1 represents the same wave speed and a1 implies different wave speeds for the first two equations corresponding to the Timoshenko-Fourier system.

We investigate the global existence of solutions to (1.3) in critical Besov spaces. Xu and Kawashima have developed a theory of local existence of general symmetric hyperbolic systems in critical Besov spaces, which is seen as a generalization of the basic theory of Kato and Majda [4] [5], and since that the non-symmetric dissipation L has no effect on the local-in-time existence, their results can be directly applied to the system (1.3). Based on the observation for the system (1.1), we argue that the viscous term k θ xx in (1.1) does not affect the mathematical entropy, allowing us to obtain a global a priori estimate by considering it from the perspective of the hyperbolic energy methods. The specific ideas are detailed in the following proof of Theorem 1.1. According to the references reviewed, very little has been done for the study of the Timoshenko-Fourier system, although many have achieved modest results for studying the fluid dynamic equations in critical Besov spaces, such as the Navier-Stokes equations in [6]-[9], and [10]-[13] for the Euler equation and related models. Due to non-symmetry of Timoshenko-Fourier system, although the first author and third author in [14] have studied general dissipative systems under the assumptions of dissipative entropy and the Shizuta-Kawashima condition in [15], their result cannot be directly applied to the Timoshenko-Fourier system. The conclusions related to the non-symmetry of relaxation matrices have been studied in [11], which provides the theoretical support for our research. Therefore, we first consider the case that (1.3) has the equal wave speed ( a=1 ), and then the non-equal wave speed ( a1 ), of which we can derive the desired prior estimate by using a fundamental fact that the relation between homogeneous and non-homogeneous Chemin-Lerner spaces in Proposition 2.1. Finally, we construct global solutions pertaining to data in the Besov space B 2,1 3/2 ( ) .

In the case of non-equal waves ( a1 ), there will be regularity loss (a loss of derivatives in the energy estimates), and the system (1.3) admits a weaker dissipation mechanism. To obtain similar result for (1.3) with a=1 , we will establish the global-in-time existence: s c = 3 2 for the Timoshenko-Fourier system with loss of regularity, which will lead to a significant reduction in the regularity requirement for the initial values compared to [16]. Moreover, compared to the case a=1 , when estimating the dissipation with respect to v , we must utilize Proposition 2.1 for the topological relation between u x L ˜ T 2 ( B ˙ 2,1 1/2 ) and u x L ˜ T 2 ( B 2,1 1/2 ) to obtain the dissipation estimate on v .

This paper is arranged as follows. In Section 2, the related lemmas and propositions are shown. In Section 3, we establish several lemmas of estimating Proposition (3.2), which leads to the proof of Theorem 1.1.

Notations. Throughout, we provide some notations. fg denotes fCg , where C>0 is a general constant, the number represented by C in different rows may differ. C( [ 0,T ],X ) ( C 1 ( [ 0,T ],X ) ) denotes the space of continuous (continuously differentiable) functions that take values in [ 0,T ] of the Banach space X . ( f,g,h ) X means f X + g X + h X , where f,g,hX .

Main Results

Rewrite Equation (1.1) as the following Cauchy problem [3]:

{ U t +A U x +LU=B U xx +G ( U ) x , U( x,0 )= U 0 ( x ), (1.3)

where the coefficient matrices A , L and B are given by the following formula:

A=( 0 1 0 0 0 1 0 0 0 0 0 0 0 a 0 0 0 a 0 β 0 0 0 β 0 ),L=( 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 1 0 0 γ 0 0 0 0 0 0 ),

and

B=( 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 k ),

where G( U )=( 0,0,0,g( z ),0 ) with g( z )=σ( z/a )σ( 0 ) σ ( 0 )z/a :=O( z 2 ) near z=0 .

The main result of our work is stated as follows:

Theorem 1.1. Suppose U 0 B 2,1 3/2 ( ) . There exists a positive constant δ 0 such that if

U 0 B 2,1 3/2 δ 0 ,

then the Cauchy problem (1.3) has a unique global classical solution U C ˜ 1 ( + × ) satisfying

U C ˜ ( B 2,1 3/2 ( ) ) C ˜ 1 ( B 2,1 1/2 ( ) ).

In addition, the following energy inequality holds true

U L ˜ T ( B 2,1 3/2 ( ) ) +( y L ˜ T 2 ( B 2,1 3/2 ) + ( v, z x ) L ˜ T 2 ( B 2,1 1/2 ) + u x L ˜ T 2 ( B 2,1 1/2 ) + θ x L ˜ T 2 ( B 2,1 3/2 ) ) C 0 U 0 B 2,1 3/2 ( ) ,

where C 0 is some positive constant.

Remark 1.1. Theorem 1.1 exhibits the global-in-time existence of solutions to the Timoshenko-Fourier system in critical Besov spaces with the regularity s= 3 2 , which is the optimal result for the Timoshenko-Fourier system studied so far based on the dimension of the space and the conditions for system stability.

2. Preliminary

In this section, we mainly show the lemmas and propositions that will be used in the following; for proofs of the propositions, see the references [3] [17].

Lemma 2.1. Let 0< R 1 < R 2 and 1ab .

1) If Suppose f{ ξ n :| ξ | R 1 λ } , for any α0 , then

Λ α f L b λ α+n( 1 a 1 b ) f L a ,

2) If Suppose f{ ξ n : R 1 λ| ξ | R 2 λ } , for any α , then

Λ α f L b λ α f L a .

Lemma 2.2 Let s and 1p,r . then

1) If s>0 , then B p,r s = L p B ˙ p,r s .

2) If s ˜ s , then B p,r s B p,r s ˜ . This inclusion relation is false for the homogeneous Besov spaces.

3) If 1r r ˜ , then B ˙ p,r s B ˙ p, r ˜ s and B p,r s B p, r ˜ s .

4) If 1p p ˜ , then B ˙ p,r s B ˙ p ˜ ,r sn( 1 p 1 p ˜ ) and B p,r s B p ˜ ,r sn( 1 p 1 p ˜ ) .

5) If rθ , then f L ˜ T θ ( B p,r s ) f L T θ ( B p,r s ) ; If rθ , then f L ˜ T θ ( B p,r s ) f L T θ ( B p,r s ) .

6) B p,1 n/p C 0 , B ˙ p,1 n/p C 0 ( 1p< ) .

where C 0 is the spaces of continuous bounded functions that decay at infinity.

Lemma 2.3 Suppose that ϱ>0 and 1p<2 . It holds that

f B ˙ r, ϱ f L p .

with 1 p 1 r = ϱ n . In particular, this holds with ϱ= n 2 ,r=2,p=1 .

Integral existentiality depends on the connection between homogeneous Chemin-Lerner spaces and non-homogeneous Chemin-Lerner spaces, which we will only briefly illustrate here; see [3] for more detailed proofs.

Proposition 2.1. Let s and 1θ,p,r , for any T>0

1) It holds that

L T θ ( L p ) L ˜ T θ ( B ˙ p,r s ) L ˜ T θ ( B p,r s ).

2) Moreover, as s>0 and θr , It holds that

L T θ ( L p ) L ˜ T θ ( B ˙ p,r s )= L ˜ T θ ( B p,r s ).

Proposition 2.2 Let s>0 and 1p,r , then B ˙ p,r s L is an algebra and

fg B ˙ p,r s f L g B ˙ p,r s + g L f B ˙ p,r s .

Let s 1 , s 2 n/p such that s 1 + s 2 >nmax{ 0, 2 p 1 } . Then one has

fg B ˙ p,1 s 1 + s 2 n/p f B ˙ p,1 s 1 g B ˙ p,1 s 2 .

Proposition 2.3 Let 1<p<,1θ and s( n p 1, n p ] . Then there exists a constant C>0 that depends only on s,n such that

{ [ f, Δ ˙ q ]g L p C c q 2 q( s+1 ) f B ˙ p,1 n p +1 g B ˙ p,1 s , [ f, Δ ˙ q ]g L T θ ( L p ) C c q 2 q( s+1 ) f L ˜ T θ 1 ( B ˙ p,1 n p +1 ) g L ˜ T θ 2 ( B ˙ p,1 s ) .

with 1 θ = 1 θ 1 + 1 θ 2 , where the commutators [ , ] is defined by [ f,g ]=fggf and { c q } denotes a sequence such that c q l 1 1 .

Proposition 2.4. Let s>0,1p,r,θ,F W loc [ s ]+3, ( I; ) with F( 0 )=0,T( 0, ] and f L ˜ T θ ( B p,r s ) L T ( L ) , Then there exists a function C that depends only on s,p,r,n and F such that

{ F( f ) F ( 0 )f B ˙ p,r s C( f L ) f B ˙ p,r s 2 , F( f ) F ( 0 )f L ˜ T θ ( B ˙ p,r s ) C( f L T ( L ) ) f L ˜ T θ ( B ˙ p,r s ) 2 .

Proposition 2.5. The following inequality holds

fg L ˜ T θ ( B p,r s ) ( f L T θ 1 ( L ) g L ˜ T θ 2 ( B p,r s ) + g L T θ 3 ( L ) f L ˜ T θ 4 ( B p,r s ) ).

with s>0,1p,1θ, θ 1 , θ 2 , θ 3 , θ 4 and

1 θ = 1 θ 1 + 1 θ 2 = 1 θ 3 + 1 θ 4 .

A direct corollary is that

fg L ˜ T θ ( B p,r s ) f L ˜ T θ 1 ( B p,r s ) g L ˜ T θ 2 ( B p,r s )

with sn/p , 1 θ = 1 θ 1 + 1 θ 2 .

Proposition 2.6. Let s>0,1p,r,ρ,F W loc [ s ]+1, ( I; ) with F( 0 )=0,T( 0, ] and v L ˜ T ρ ( B p,r s ) L T ( L ) , then

F( v ) L ˜ T ρ ( B p,r s ) ( 1+ v L T ( L ) ) [ s+1 ] v L ˜ T ρ ( B p,r s ) .

3. The Proof of Theorem 1.1

Xu and Kawashimsa [3] established a local existence theory for general symmetric hyperbolic systems in critical Besov spaces, which is seen as a generalization of the basic theory established by Kato and Majda [4] [5], and it can be applied to Timoshenko-Fourier system (1.3). Precisely,

Proposition 3.1. Assuming U 0 B 2,1 3/2 , there exists a time T 0 >0 (dependent only on initial data) such that

1) (Existence) The system (1.3) has a unique solution U( t,x ) C 1 ( [ 0, T 0 ]× ) that satisfying U C ˜ T 0 ( B 2,1 3/2 ) C ˜ T 0 1 ( B 2,1 1/2 ) ;

2) (Blow-up Criterion) If the maximal time T * ( > T 0 ) existing of such a solution in finite, then

limsup n T * U( t, ) B 2,1 3/2 =,

if and only if

0 T * ( U( t, ),U( t, ) ) L dt =.

Proof. In order to prove Proposition (3.1) (i), it suffices to show that U C ˜ T 0 ( B 2,1 3/2 ) and U t C ˜ T 0 ( B 2,1 1/2 ) .

We used the stationary cases of estimates of commutator in [18] (Lemma 2.100, p. 112) to (3.26), we arrive at

1 2 d dt N 0 [ Δ ˙ q U ]+γ Δ ˙ q y L 2 2 +k Δ ˙ q θ x L 2 2 y x L Δ ˙ q z L 2 2 + z x L Δ ˙ q z L 2 Δ ˙ q y L 2 + 2 q 2 ( z xx L z B 2,1 1 2 + z x L z xx B 2,1 1 2 ) Δ q y L 2 . (3.4)

Dividing (3.4) by ( Δ q U L 2 2 +ϵ ) 1 2 , we can obtain

d dt ( Δ q U L 2 2 +ϵ ) 1 2 + Δ ˙ q y L 2 Δ ˙ q θ x L 2 + y x L Δ ˙ q z L 2 + z x L Δ ˙ q y L 2 + 2 q 2 ( z xx L z B 2,1 1 2 + z x L z xx B 2,1 1 2 ). (3.5)

Integrating (3.5) on the variable t[ 0, T 0 ] , then taking ϵ0 , and using the estimates of commutators and continuity for the composition in the stationary case (see, e.g., [19]), we have

2 1 2 q Δ q U L t ( L 2 ) 2 1 2 q Δ q U 0 L 2 + 0 t c q ( τ ) ( y, θ x ) B 2,1 1 2 2 dτ + 0 t ( z xx , z x ) B 2,1 1 2 z B 2,1 3 2 dτ , (3.6)

where we used the embedding inequality ( z L z B 2,1 1 2 ) , Lemma 2.2 and Young’s inequality. Summing up (3.6) on q1 gives

U L ˜ t ( B 2,1 3/2 ) U 0 B 2,1 3/2 +ϵ 0 t U,ΔU B 2,1 1/2 2 dτ +C( ϵ ) 0 t U B 2,1 3/2 2 dτ . (3.7)

Let ϵ>0 be a small number, using Gronwall’s inequality, we get

U C ˜ T 0 ( B 2,1 3/2 ). (3.8)

Furthermore, with proposition (2.5) we can deduce that

U t C ˜ T 0 ( B 2,1 1/2 ). (3.9)

This completes Proposition (3.1) (i).

From [3], for a symmetric hyperbolic system it is sufficient to establish the blow-up criterion. Such a criterion can be established for the Timoshenko-Fourier system. We consider the symmetric system (1.3) with LU0 for simplicity, since it is only responsible for the global well-posedness and large time behavior of solutions

t U+A x U=B U xx +G ( U ) x . (3.10)

Applying the homogeneous operator Δ ˙ q to (1.3), we infer that Δ ˙ q U satisfies

t Δ ˙ q U+A Δ ˙ q x U+B Δ ˙ q U xx + Δ ˙ q G ( U ) x =0. (3.11)

Perform the inter product with Δ ˙ q U on both sides of (3.11) to get

Δ ˙ q U, Δ ˙ q U t + A Δ ˙ q U x , Δ ˙ q U + B Δ ˙ q U x , Δ ˙ q U x + Δ ˙ q G ( U ) x , Δ ˙ q U = B Δ ˙ q U x , Δ ˙ q U x , (3.12)

where

Δ ˙ q G ( U ) x , Δ ˙ q U [ Δ ˙ q ,U ] U x , Δ ˙ q U + U Δ ˙ q U, Δ ˙ q U x .

By integrating (3.12) with respect to x over , we deduce that

d dt Δ ˙ q U L 2 2 Δ ˙ q U L 2 Δ ˙ q U x L 2 + [ Δ ˙ q ,U ] U x L 2 Δ ˙ q U L 2 + U L Δ ˙ q U L 2 Δ ˙ q U x L 2 . (3.13)

Let ϵ>0 be a small number, dividing (3.13) by ( Δ ˙ q U L 2 2 +ϵ ) 1 2 gives

d dt ( Δ ˙ q U L 2 2 +ϵ ) 1 2 Δ ˙ q U x L 2 + [ Δ ˙ q ,U ] U x L 2 + U L Δ ˙ q U x L 2 c q ( t ) 2 q 2 ( U x L U B ˙ 2,1 1 2 + U x L U x B ˙ 2,1 1 2 )+ U L Δ ˙ q U x L 2 , (3.14)

where we used the stationary cases of estimates of commutator in [18] (Lemma 2.100, p. 112) and the sequence { c q ( t ) } satisfying c q ( t ) l 1 1 for all t[ 0, T 0 ] .

Letting ϵ0 and taking a time integration, we are lead to

U( t ) B ˙ 2,1 1 2 U 0 B ˙ 2,1 1 2 +C 0 T 0 ( U,U ) L ( U,U ) B ˙ 2,1 1 2 dτ . (3.15)

Taking the L 2 inner production (3.10) with U , it is easy to obtain

U( t ) L 2 U 0 L 2 +C 0 T 0 ( U,U ) L U,U L 2 dτ . (3.16)

Adding (3.15) to (3.16), with the aid of Proposition 2.1(1) and Lemma 2.1, we have

U( t ) B 2,1 1 2 U 0 B 2,1 1 2 +C 0 T 0 ( U,U ) L U,U B 2,1 1 2 dτ U 0 B 2,1 1 2 +C 0 T 0 ( U,U ) L U B 2,1 3 2 dτ . (3.17)

The Gronwall’s inequality implies

sup t[ 0, T 0 ] U( t ) B 2,1 3 2 U 0 B 2,1 3 2 exp( C 0 T 0 ( U,U ) L dτ ). (3.18)

So, the subsequent inequalities will hold true

0 T 0 ( U,U ) L dτ C 0 T 0 U B 2,1 3 2 dτ C T 0 sup t[ 0, T 0 ] U B 2,1 3 2 . (3.19)

This completes the proof of Proposition 3.1(ii).

Moreover, to prove that the classical solution in Proposition 3.1 is globally defined, we need to construct a priori estimates based on the dissipation mechanism generated by the Timoshenko-Fourier system. To this end, we define energy generalization in terms of N( T ) and the corresponding dissipation generalization in terms of D( T ) :

N( T ):= U L ˜ T ( B 2,1 3/2 ) ,

D( T ):= y L ˜ T 2 ( B 2,1 3/2 ) + ( v, z x ) L ˜ T 2 ( B 2,1 1/2 ) + u x L ˜ T 2 ( B 2,1 1/2 ) + θ x L ˜ T 2 ( B 2,1 3/2 )

for any time T>0 .

In the case of a=1 :

Lemma 3.1 (The dissipation for ( y, θ x ) ) If U C ˜ T ( B 2,1 3/2 ) C ˜ T 1 ( B 2,1 1/2 ) is a solution of (1.3) for any T>0 , then

E( T )+ y L ˜ T 2 ( B 2,1 3/2 ) + θ x L ˜ T 2 ( B 2,1 3/2 ) U 0 B 2,1 3/2 + N( T ) D( T ).

Proof. We start with the next equations

{ v t u x +y=0, u t v x =0, z t y x =0, y t [ σ( z ) ] x v+γy+β θ x =0, θ t k θ xx +β y x =0. (3.20)

The five equations of the above system (3.20) are each multiplied in order by v,u,σ( z )σ( 0 ),y and θ . Adding the resulting equations yields

1 2 d dt ( v 2 + u 2 +F( z )+ y 2 + θ 2 ) ( vu+( σ( z )σ( 0 ) )y+β θ y +k θ x θ ) x +γ y 2 +k θ x 2 =0, (3.21)

where

F( z )=2 0 z ( σ( η )σ( 0 ) )dη .

Noting that F( z ) is equivalent to z 2 , due to σ ( η )>0 and the following smallness assumption, then integrating about x yields the basic energy equation

1 2 d dt N 0 ( U )+γ y L 2 2 +k θ x L 2 2 =0, (3.22)

where the energy generalized function N 0 ( U ) is defined as

N 0 ( U )= ( v,u,y,θ ) L 2 2 + F( z )dx U L 2 2 .

By integrating over t[ 0,T ] and then taking the square-root of the resultant inequality, we get

U L T ( L 2 ) + 2γ y L T 2 ( L 2 ) + 2k θ x L T 2 ( L 2 ) U 0 L 2 (3.23)

for any of T>0 . Next, the dissipation rate of y, θ x is obtained by frequency localization estimation in homogeneous Chemin-Lerner space. Applying the operator Δ ˙ q ( q ) to (3.20), we have

{ Δ ˙ q v t Δ ˙ q u x + Δ ˙ q y=0, Δ ˙ q u t Δ ˙ q v x =0, Δ ˙ q z t Δ ˙ q y x =0, Δ ˙ q y t σ ( z ) Δ ˙ q z x Δ ˙ q v+γ Δ ˙ q y+β Δ ˙ q θ x =[ Δ ˙ q , σ ( z ) ] z x , Δ ˙ q θ t k Δ ˙ q θ xx +β Δ ˙ q y x =0, (3.24)

where the commutator is defined as [ f,g ]:=fggf . Multiplying (3.24) with Δ ˙ q v, Δ ˙ q u, σ ( z ) Δ ˙ q z x , Δ ˙ q y and Δ ˙ q θ , respectively, and then adding the resulting equations. We obtain

t l + x l +γ ( Δ ˙ q y ) 2 +k ( Δ ˙ q θ x ) 2 = l , (3.25)

where z t = y x , and

l := 1 2 { ( Δ ˙ q v ) 2 + ( Δ ˙ q u ) 2 + σ ( z ) ( Δ ˙ q z ) 2 + ( Δ ˙ q y ) 2 + ( Δ ˙ q θ ) 2 }, l := Δ ˙ q v Δ ˙ q u+ Δ ˙ q y Δ ˙ q θ+ σ ( z )( Δ ˙ q z )( Δ ˙ q y )+k Δ ˙ q θ x Δ ˙ q θ, l := 1 2 σ ( z ) t ( Δ ˙ q z ) 2 σ ( z ) x ( Δ ˙ q z )( Δ ˙ q y )+ Δ ˙ q y[ Δ ˙ q , σ ( z ) ] z x .

Further, integrating over x and with the aid of the Cauchy-Schwarz inequality yield

1 2 d dt N 0 [ Δ ˙ q U ]+γ Δ ˙ q y L 2 2 +k Δ ˙ q θ x L 2 2 σ ( z ) t L Δ ˙ q z L 2 2 + σ ( z ) x L Δ ˙ q z L 2 Δ ˙ q y L 2 + [ Δ ˙ q , σ ( z ) ] z x L 2 Δ ˙ q y L 2 , (3.26)

where

N 0 [ Δ ˙ q U ]= ( Δ ˙ q v, Δ ˙ q u, Δ ˙ q y, Δ ˙ q θ ) L 2 2 + σ ( z ) | Δ ˙ q z | 2 dx Δ ˙ q U L 2 2 .

From the above (1.3) and the following a priori assumption (3.69), we have

σ ( z ) t L Δ ˙ q z L 2 2 z t L Δ ˙ q z L 2 2 y x L Δ ˙ q z L 2 2 . (3.27)

Similarly,

σ ( z ) x L Δ ˙ q z L 2 Δ ˙ q y L 2 z x L Δ ˙ q z L 2 Δ ˙ q y L 2 . (3.28)

Combining (3.27)-(3.28) and integrating over t[ 0,T ] , with the aid of Young’s inequality we are led to

N 0 [ Δ ˙ q U ] + 2γ Δ ˙ q y L T 2 ( L 2 ) + 2k Δ ˙ q θ x L T 2 ( L 2 ) N 0 [ Δ ˙ q U 0 ] + ( y x , z x ) L T ( L ) ( Δ ˙ q y L T 2 ( L 2 ) + Δ ˙ q z L T 2 ( L 2 ) ) + [ Δ ˙ q , σ ( z ) ] z x L T 2 ( L 2 ) Δ ˙ q y L T 2 ( L 2 ) . (3.29)

By the commutator estimate in Proposition 2.3, it holds

[ Δ ˙ q , σ ( z ) ] z x L T 2 ( L 2 ) c q 2 3q 2 z L ˜ T ( B ˙ 2,1 3/2 ) z x L ˜ T 2 ( B ˙ 2,1 1/2 ) , (3.30)

where { c q } denotes a sequence satisfying c q l 1 1 . Therefore, we get

2 3q 2 Δ ˙ q U L T ( L 2 ) + 2γ 2 3q 2 Δ ˙ q y L T 2 ( L 2 ) + 2k 2 3q 2 Δ ˙ q θ x L T 2 ( L 2 ) Δ ˙ q U 0 L 2 + c q ( y x , z x ) L T ( B ˙ 2,1 1/2 ) ( y L ˜ T 2 ( B ˙ 2,1 3/2 ) + z x L ˜ T 2 ( B ˙ 2,1 1/2 ) ) + c q z L ˜ T ( B ˙ 2,1 3/2 ) ( y L ˜ T 2 ( B ˙ 2,1 3/2 ) z x L ˜ T 2 ( B ˙ 2,1 1/2 ) ) (3.31)

Here, we would like to point out that each { c q } may have a different form in the (3.31) equation or the inequality that emerges after, but the bound of c q l 1 1 is well satisfied. Thus, summing over q , we have

U L ˜ T ( B ˙ 2,1 3/2 ) + 2γ y L ˜ T 2 ( B ˙ 2,1 3/2 ) + 2k θ x L ˜ T 2 ( B ˙ 2,1 3/2 ) U 0 ( B ˙ 2,1 3/2 ) + ( y,z ) L ˜ T ( B ˙ 2,1 3/2 ) ( y L ˜ T 2 ( B ˙ 2,1 3/2 ) + z x L ˜ T 2 ( B ˙ 2,1 1/2 ) ). (3.32)

Finally, combining (3.23) and (3.32), we conclude from Proposition 2.2 that

N( T )+ y L ˜ T 2 ( B 2,1 3/2 ) + θ x L ˜ T 2 ( B 2,1 3/2 ) U 0 B 2,1 3/2 + N( T ) D( T ). (3.33)

Thus, the proof of Lemma 3.1 is completed.

Lemma 3.2. (The dissipation for v ) If U C ˜ T ( B 2,1 3/2 ) C ˜ T 1 ( B 2,1 1/2 ) is a solution of (1.3) for any T>0 , then

v L ˜ T 2 ( B 2,1 1/2 ) N( T )+ U 0 B 2,1 3/2 + y L ˜ T 2 ( B 2,1 3/2 ) + θ x L ˜ T 2 ( B 2,1 3/2 ) + N( T ) D( T ). (3.34)

Proof. The system (1.3) can be rewritten in the following form:

{ v t u x +y=0, u t v x =0, z t y x =0, y t z x v+γy+β θ x =g ( z ) x , θ t k θ xx +β y x =0, (3.35)

where the smooth function g( z ) is defined as

g( z )=σ( z )σ( 0 )z=O( z 2 )

and satisfies g( 0 )=0 and g ( 0 )=0 . Multiplying the first equation in (3.35) by y , the second by z , the third by u , and the fourth by v , respectively, and then adding the resulting equations together give

d dt N 1 ( U )+ v L 2 2 y L 2 2 +γ y L 2 v L 2 +β v L 2 θ x L 2 + g ( z ) x L 2 v L 2 , (3.36)

where

N 1 ( U )= ( vy+uz )dx .

Using the Young’s inequality, we can conclude that

d dt N 1 ( U )+ 1 2 v L 2 2 y L 2 2 + θ x L 2 2 + z L z x L 2 v L 2 . (3.37)

Integration of (3.37) in t[ 0,T ] leads to

v L t 2 ( L 2 ) 2 ( | N 1 ( U ) |+| N 1 ( U 0 ) | )+ y L t 2 ( L 2 ) 2 + θ x L t 2 ( L 2 ) 2 + z L t ( L ) z x L t 2 ( L 2 ) v L t 2 ( L 2 ) N ( T ) 2 + U 0 B 2,1 3/2 2 + y L T 2 ( L 2 ) 2 + θ x L T 2 ( L 2 ) 2 +N( T )D ( T ) 2 , (3.38)

for any T>0 , here we utilize the embedding property in Lemma 2.2. Then, by Young’s inequality again, we arrive at

v L T 2 ( L 2 ) N( T )+ U 0 B 2,1 3/2 + y L T 2 ( L 2 ) + θ x L T 2 ( L 2 ) + N( T ) D( T ). (3.39)

Next, Applying the operator Δ ˙ q ( q ) to (3.35) yields

{ Δ ˙ q v t Δ ˙ q u x + Δ ˙ q y=0, Δ ˙ q u t Δ ˙ q v x =0, Δ ˙ q z t Δ ˙ q y x =0, Δ ˙ q y t Δ ˙ q z x Δ ˙ q v+γ Δ ˙ q y+β Δ ˙ q θ x = Δ ˙ q g ( z ) x , Δ ˙ q θ t k Δ ˙ q θ xx +β Δ ˙ q y x =0. (3.40)

Multiplying the first equation in (3.40) by Δ ˙ q y , the second by Δ ˙ q z , the third by Δ ˙ q u , and the fourth by Δ ˙ q v , respectively, and then adding the resulting equations. we get

t 1 l + x 1 l + ( Δ ˙ q v ) 2 ( Δ ˙ q y ) 2 = 1 l , (3.41)

where

1 l := Δ ˙ q v Δ ˙ q y+ Δ ˙ q u Δ ˙ q z, 1 l := Δ ˙ q v Δ ˙ q z+ Δ ˙ q u Δ ˙ q y, 1 l :=γ Δ ˙ q v Δ ˙ q y+β Δ ˙ q v Δ ˙ q θ x Δ ˙ q v Δ ˙ q g ( z ) x .

Using the Hölder inequality and Young’s inequality gives

d dt N 1 [ Δ ˙ q U ]+ 1 2 Δ ˙ q v L 2 2 Δ ˙ q y L 2 2 + Δ ˙ q θ x L 2 2 + Δ ˙ q g ( z ) x L 2 Δ ˙ q v L 2 , (3.42)

where

N 1 [ Δ ˙ q U ]= ( Δ ˙ q v Δ ˙ q y+ Δ ˙ q u Δ ˙ q z )dx .

By integrating over t[ 0,T ] , we have

Δ ˙ q v L t 2 ( L 2 ) 2 Δ ˙ q U L T ( L 2 ) 2 + Δ ˙ q U 0 L 2 2 + Δ ˙ q y L T 2 ( L 2 ) 2 + Δ ˙ q θ x L T 2 ( L 2 ) 2 + Δ ˙ q v L T ( L 2 ) Δ ˙ q g ( z ) x L T 1 ( L 2 ) . (3.43)

In addition, Young’s inequality allows us to obtain

2 q 2 Δ ˙ q v L T 2 ( L 2 ) c q U L ˜ T ( B ˙ 2,1 1/2 ) + c q U 0 B ˙ 2,1 1/2 + c q y L ˜ T 2 ( B ˙ 2,1 1/2 ) + c q θ x L ˜ T 2 ( B ˙ 2,1 1/2 ) + c q v L ˜ T ( B ˙ 2,1 1/2 ) g ( z ) x L ˜ T 1 ( B ˙ 2,1 1/2 ) 1 2 , (3.44)

where the norm of g ( z ) x at the right end of (3.44) can be estimated by applying Lemma 2.2 and Proposition 2.4

g ( z ) x L ˜ T 1 ( B ˙ 2,1 1/2 ) 0 T g( z ) B ˙ 2,1 3/2 dt 0 T z B ˙ 2,1 3/2 2 dt z x L ˜ T 2 ( B ˙ 2,1 1/2 ) 2 . (3.45)

Thus, combining (3.44)-(3.45) and summing q , we get

v L ˜ T 2 ( B ˙ 2,1 1/2 ) U L ˜ T ( B ˙ 2,1 1/2 ) + U 0 B ˙ 2,1 1/2 + y L ˜ T 2 ( B ˙ 2,1 1/2 ) + v L ˜ T ( B ˙ 2,1 1/2 ) z x L ˜ T 2 ( B ˙ 2,1 1/2 ) + θ x L ˜ T 2 ( B ˙ 2,1 1/2 ) . (3.46)

Finally, note that (3.39) and (3.46), by proposition 2.1, we obtain

v L ˜ T 2 ( B 2,1 1/2 ) N( T )+ U 0 B 2,1 3/2 + y L ˜ T 2 ( B 2,1 3/2 ) + θ x L ˜ T 2 ( B 2,1 3/2 ) + N( T ) D( T ). (3.47)

This proves Lemma 3.2.

Lemma 3.3 (The dissipation for z x ) If U C ˜ T ( B 2,1 3/2 ) C ˜ T 1 ( B 2,1 1/2 ) is a solution of (1.3) for any T>0 , then

z x L ˜ T 2 ( B 2,1 1/2 ) N( T )+ U 0 B 2,1 3/2 + y L ˜ T 2 ( B 2,1 3/2 ) + v L ˜ T 2 ( B 2,1 1/2 ) + θ x L ˜ T 2 ( B 2,1 3/2 ) + N( T ) D( T ). (3.48)

Proof. Multiplying the third equation of (3.35) by y x and the fourth equation by z x , respectively, and then adding up the resulting equations, we have

d dt ( z x y )+ x ( y z t )+ ( z x ) 2 ( y x ) 2 = 2 l , (3.49)

where

2 l := z x ( γy+β θ x v ) z x g ( z ) x .

Integrating (3.49) over , we get

z x L 2 2 y x L 2 2 + d dt z x yd x+ | 2 l |dx ,

| 2 l |dx ϵ z x L 2 2 + C ϵ ( v,y, θ x ) L 2 2 +C z L z x L 2 2 .

That is

d dt N 2 ( U )+ z x L 2 2 y x L 2 2 +( v L 2 2 + y L 2 2 + θ x L 2 2 )+ z L z x L 2 2 , (3.50)

where

N 2 ( U ):= z x ydx .

Similar to the process for (3.39), we obtain

z x L T 2 ( L 2 ) N( T )+ U 0 B 2,1 3/2 + y L ˜ T 2 ( B 2,1 3/2 ) + v L ˜ T 2 ( B 2,1 1/2 ) + θ x L ˜ T 2 ( B 2,1 1/2 ) + N( T ) D( T ). (3.51)

Furthermore, by (3.40), it holds

{ Δ ˙ q z t Δ ˙ q y x =0, Δ ˙ q y t Δ ˙ q z x Δ ˙ q v+γ Δ ˙ q y+β Δ ˙ q θ x = Δ ˙ q g ( z ) x . (3.52)

The equations (3.52) are multiplied by Δ ˙ q y x and Δ ˙ q z x , respectively, and then utilizing the energy estimates on each block, we get

2 q 2 Δ ˙ q z x L T 2 ( L 2 ) c q ( U L ˜ T ( B 2,1 3/2 ) + U 0 B 2,1 3/2 )+ c q y x L ˜ T 2 ( B ˙ 2,1 1/2 ) + c q ϵ z x L ˜ T 2 ( B ˙ 2,1 1/2 ) + c q C ϵ ( v L ˜ T 2 ( B ˙ 2,1 1/2 ) + y L ˜ T 2 ( B ˙ 2,1 1/2 ) + θ x L ˜ T 2 ( B ˙ 2,1 1/2 ) ) + c q z x L ˜ T ( B ˙ 2,1 1/2 ) g ( z ) x L ˜ T 1 ( B ˙ 2,1 1/2 ) 1 2 . (3.53)

In addition, similar to the processes of (3.45)-(3.46), we arrive at

z x L T 2 ( B ˙ 2,1 1/2 ) U L ˜ T ( B 2,1 3/2 ) + U 0 B 2,1 3/2 + y L ˜ T 2 ( B ˙ 2,1 3/2 ) + v L ˜ T 2 ( B ˙ 2,1 1/2 ) + y L ˜ T 2 ( B ˙ 2,1 1/2 ) + θ x L ˜ T 2 ( B ˙ 2,1 3/2 ) + c q z L ˜ T ( B ˙ 2,1 3/2 ) z x L ˜ T 2 ( B ˙ 2,1 1/2 ) , (3.54)

where choosing 0<ϵ 1 2 .

Eventually, combining (3.51) and (3.54) gives Lemma 3.3.

Lemma 3.4. (The dissipation for u x ) If U C ˜ T ( B 2,1 3/2 ) C ˜ T 1 ( B 2,1 1/2 ) is a solution of (1.3) for any T>0 , then

u x L ˜ T 2 ( B 2,1 1/2 ) N( T )+ U 0 B 2,1 3/2 + v L ˜ T 2 ( B 2,1 1/2 ) + y L ˜ T 2 ( B 2,1 3/2 ) . (3.55)

Proof. Applying the nonhomogeneous operator Δ q ( q1 ) to the first and second equations of (3.35) yields

{ Δ q v t Δ q u x + Δ q y=0, Δ q u t Δ q v x =0. (3.56)

Multiplying (3.56) by Δ q u x and Δ q v x , we get

d dt N 2 [ Δ q U ]+ Δ q u x L 2 2 Δ q v x L 2 2 + Δ q u x L 2 Δ q y L 2 , (3.57)

where

N 2 [ Δ q U ]:= Δ q v Δ q u x dx .

Further, integration of (3.57) with respect to t[ 0,T ] leads to

Δ q u x L t 2 ( L 2 ) 2 | N 2 [ Δ q U ] |+ N 2 [ Δ q U 0 ]+ Δ q v x L t 2 ( L 2 ) 2 + Δ q u x L t 2 ( L 2 ) Δ q y L t 2 ( L 2 ) . (3.58)

With the aid of Young’s inequality and the embedding property in Lemma 2.2, we obtain

2 q 2 Δ q u x L T 2 ( L 2 ) c q N( T )+ c q U 0 B 2,1 3/2 + c q v L ˜ T 2 ( B 2,1 1/2 ) + c q u x L ˜ T 2 ( B 2,1 1/2 ) y L ˜ T 2 ( B 2,1 3/2 ) . (3.59)

Thus, we finish the proof of Lemma 3.4.

The case of a1 : The dissipation with respect to y, θ x , z x , u x is the same as in the case of equal waves ( a=1 ). Here we only need to show the dissipation with respect to v in the case of non-equal waves ( a1 ). It is needed to use proposition 2.1 for the topological relation between u x L ˜ T 2 ( B ˙ 2,1 1/2 ) and u x L ˜ T 2 ( B 2,1 1/2 ) to obtain the dissipative estimate of v .

Lemma 3.5. (The dissipation for v) If U C ˜ T ( B 2,1 3/2 ) C ˜ T 1 ( B 2,1 1/2 ) is a solution of (1.3) for any T>0 , then

v L ˜ T 2 ( B 2,1 1/2 ) N( T )+ U 0 B 2,1 3/2 +ε u x L ˜ T 2 ( B 2,1 1/2 ) +( 1+ C ε ) y L ˜ T 2 ( B 2,1 3/2 ) + C ε 1 θ x L ˜ T 2 ( B 2,1 3/2 ) +N( T )D( T ), (3.60)

for any ε, ε 1 >0 , C( ε ),C( ε 1 ) is a positive constant here depending on ε, ε 1 >0 .

Proof. Rewrite the system (1.3) in the following form:

{ v t u x +y=0, u t v x =0, z t a y x =0, y t a z x v+γy+β θ x =g ( z ) x , θ x k θ xx +β y x =0, (3.61)

where the smooth function g( z ) is defined by

g( z )=σ( z/a )σ( 0 ) σ ( 0 )z/a =O( z 2 )

satisfying g( 0 )=0 and g ( 0 )=0 .

First, applying the non-homogeneous frequency localization operator Δ q ( q1 ) to (3.61) leads to

{ Δ q v t Δ q u x + Δ q y=0, Δ q u t Δ q v x =0, Δ q z t a Δ q y x =0, Δ q y t a Δ q z x Δ q v+γ Δ q y+β Δ q θ x = Δ q g ( z ) x , Δ q θ x k Δ q θ xx + Δ q β y x =0. (3.62)

Subsequently, multiplying the first equation in (3.62) by Δ q y , the second equation by a Δ q z , the third equation by a Δ q u , and the fourth equation by Δ q v , and then adding up the resulting equations, we have

Δ t 3 l + x 3 l + ( Δ q v ) 2 ( Δ q y ) 2 = 3 l , (3.63)

where

3 l := Δ q v Δ q y+a Δ q u Δ q z, 3 l :=a Δ q v Δ q z+ a 2 Δ q u Δ q y, 3 l :=γ Δ q v Δ q y+β Δ q v Δ q θ x +( a 2 1 ) Δ q u x Δ q y Δ q v Δ q g ( z ) x .

Integrating (3.63) over x , using the Cauchy-Schwarz inequality, we get

d dt N 3 [ Δ q U ]+ 1 2 Δ q v L 2 2 Δ q y L 2 2 + Δ q v L 2 Δ q θ x L 2 +| a 2 1 | Δ q u x L 2 Δ q y L 2 + Δ q v L 2 Δ q g ( z ) x L 2 , (3.64)

where

N 3 [ Δ q U ]:= Δ q v Δ q y +a Δ q u Δ q z.

By integrating over t[ 0,T ] , we arrive at

Δ q v L t 2 ( L 2 ) 2 Δ q U L T ( L 2 ) 2 + Δ q U 0 L 2 2 + Δ q y L T 2 ( L 2 ) 2 + Δ q v L T 2 ( L 2 ) Δ q θ x L T 2 ( L 2 ) + Δ q y L T 2 ( L 2 ) Δ q u x L T 2 ( L 2 ) + Δ q g ( z ) x L T 2 ( L 2 ) , (3.65)

where noting the case of a1 . In addition, utilizing the Young’s inequality gives

2 q 2 Δ q v L T 2 ( L 2 ) c q U L ˜ T ( B 2,1 1/2 ) + c q U 0 B 2,1 1/2 +ε c q u x L ˜ T 2 ( B 2,1 1/2 ) + c q ( 1+ C ε ) y L ˜ T 2 ( B 2,1 3/2 ) + C ε 1 c q θ x L ˜ T 2 ( B 2,1 1/2 ) + c q g ( z ) x L ˜ T 2 ( B 2,1 1/2 ) , (3.66)

where

2 q Δ q v L T 2 ( L 2 ) Δ q θ x L T 2 ( L 2 ) ε 1 v L T 2 ( B 2,1 1/2 ) 2 + C ε 1 θ x L T 2 ( B 2,1 1/2 ) 2 ,

for ε, ε 1 >0 , where C ε , C ε 1 are position constants dependent on ε, ε 1 . Knowing the fact g ( 0 )=0 , it follows from Proposition 2.5 and Proposition 2.6 that

g ( z ) x L ˜ T 2 ( B 2,1 1/2 ) = g ( z ) z x L ˜ T 2 ( B 2,1 1/2 ) g ( z ) g ( 0 ) L ˜ T ( B 2,1 1/2 ) z x L ˜ T 2 ( B 2,1 1/2 ) z L ˜ T ( B 2,1 1/2 ) z x L ˜ T 2 ( B 2,1 1/2 ) . (3.67)

Combining the above estimations (3.66)-(3.67), and summing over q1 , we obtain

v L ˜ T 2 ( B 2,1 1/2 ) U L ˜ T ( B 2,1 1/2 ) + U 0 B 2,1 1/2 +ε u x L ˜ T 2 ( B 2,1 1/2 ) +( 1+ C ε ) y L ˜ T 2 ( B 2,1 3/2 ) + C ε 1 θ x L ˜ T 2 ( B 2,1 3/2 ) + z L ˜ T ( B 2,1 1/2 ) z x L ˜ T 2 ( B 2,1 1/2 ) . (3.68)

This gives (3.60).

From the above Lemma 3.1-Lemma 3.5, the following Proposition 3.2 holds.

Proposition 3.2. Suppose that for any T>0 , U C ˜ T ( B 2,1 3/2 ) C ˜ T 1 ( B 2,1 1/2 ) is a solution of (1.3), and there exists δ 1 >0 such that when

N( T ) δ 1 , (3.69)

then the following estimate holds

N( T )+D( T ) U 0 B 2,1 3/2 +( N( T )+ N( T ) )D( T ). (3.70)

Therefore, the following inequality holds

N( T )+D( T ) U 0 B 2,1 3/2 . (3.71)

By using the standard boot-strap argument, the proof of Theorem 1.1 is similar to the process in [20], and we omit the details here for the sake of brevity.

Acknowledgements

The first author (H. M. Cao) is supported by the National Natural Science Foundation of China (No. 12001269) and the Fundamental Research Funds for the Central Universities of China.

Conflicts of Interest

The authors declare no conflicts of interest regarding the publication of this paper.

References

[1] Timoshenko, S.P. (1921) LXVI. On the Correction for Shear of the Differential Equation for Transverse Vibrations of Prismatic Bars. The London, Edinburgh, and Dublin Philosophical Magazine and Journal of Science, 41, 744-746.[CrossRef]
[2] Timoshenko, S.P. (1922) X. On the Transverse Vibrations of Bars of Uniform Cross-section. The London, Edinburgh, and Dublin Philosophical Magazine and Journal of Science, 43, 125-131.[CrossRef]
[3] Xu, J. and Kawashima, S. (2014) Global Classical Solutions for Partially Dissipative Hyperbolic System of Balance Laws. Archive for Rational Mechanics and Analysis, 211, 513-553.[CrossRef]
[4] Kato, T. (1975) The Cauchy Problem for Quasi-Linear Symmetric Hyperbolic Systems. Archive for Rational Mechanics and Analysis, 58, 181-205.[CrossRef]
[5] Majda, A. (1984) Compressible Fluid Flow and Conservation Laws in Several Space Variables. Springer-Verlag.
[6] Abidi, H., Gui, G. and Zhang, P. (2012) On the Wellposedness of Three-Dimensional Inhomogeneous Navier-Stokes Equations in the Critical Spaces. Archive for Rational Mechanics and Analysis, 204, 189-230.[CrossRef]
[7] Danchin, R. (2000) Global Existence in Critical Spaces for Compressible Navier-Stokes Equations. Inventiones Mathematicae, 141, 579-614.[CrossRef]
[8] Haspot, B. (2011) Existence of Global Strong Solutions in Critical Spaces for Barotropic Viscous Fluids. Archive for Rational Mechanics and Analysis, 202, 427-460.[CrossRef]
[9] Paicu, M. and Zhang, P. (2012) Global Solutions to the 3-D Incompressible Inhomogeneous Navier–Stokes System. Journal of Functional Analysis, 262, 3556-3584.[CrossRef]
[10] Danchin, R. (2010) On the Well-Posedness of the Incompressible Density-Dependent Euler Equations in the Lp Framework. Journal of Differential Equations, 248, 2130-2170.
[11] Xu, J. and Yong, W.A. (2009) Relaxation-Time Limits of Non-Isentropic Hydrodynamic Models for Semiconductors. Journal of Differential Equations, 247, 1777-1795.[CrossRef]
[12] Xu, J. and Wang, Z.J. (2013) Relaxation Limit in Besov Spaces for Compressible Euler Equations. Journal de Mathématiques Pures et Appliquées, 99, 43-61.[CrossRef]
[13] Zhou, Y. (2004) Local Well-Posedness for the Incompressible Euler Equations in the Critical Besov Spaces. Annales de linstitut Fourier, 54, 773-786.[CrossRef]
[14] Xu, J., Mori, N. and Kawashima, S. (2015) Global Existence and Minimal Decay Regularity for the Timoshenko System: The Case of Non-Equal Wave Speeds. Journal of Differential Equations, 259, 5533-5553.[CrossRef]
[15] Xu, J. and Kawashima, S. (2015) The Optimal Decay Estimates on the Framework of Besov Spaces for Generally Dissipative Systems. Archive for Rational Mechanics and Analysis, 218, 275-315.[CrossRef]
[16] Ide, K. and Kawashima, S. (2008) Decay Property of Regularity-Loss Type and Nonlinear Effects for Dissipative Timoshenko System. Mathematical Models and Methods in Applied Sciences, 18, 1001-1025.
[17] Fernández Sare, H.D. and Racke, R. (2009) On the Stability of Damped Timoshenko Systems: Cattaneo versus Fourier Law. Archive for Rational Mechanics and Analysis, 194, 221-251.[CrossRef]
[18] Bahouri, H., Chemin, J.Y. and Danchin, R. (2011) Fourier Analysis and Nonlinear Partial Differential Equations. Springer-Verlag.
[19] Xu, J. (2009) Well-Posedness and Stability of Classical Solutions to the Multidimensional Full Hydrodynamic Model for Semiconductors. Communications on Pure & Applied Analysis, 8, 1073-1092.[CrossRef]
[20] Mori, N., Xu, J. and Kawashima, S. (2015) Global Existence and Optimal Decay Rates for the Timoshenko System: The Case of Equal Wave Speeds. Journal of Differential Equations, 258, 1494-1518.[CrossRef]

Copyright © 2026 by authors and Scientific Research Publishing Inc.

Creative Commons License

This work and the related PDF file are licensed under a Creative Commons Attribution 4.0 International License.