Research

# Global existence and uniform decay for the one-dimensional model of thermodiffusion with second sound

Ming Zhang1* and Yuming Qin2

Author Affiliations

1 College of Information Science and Technology, Donghua University, Shanghai, 201620, P.R. China

2 Department of Applied Mathematics, Donghua University, Shanghai, 201620, P.R. China

For all author emails, please log on.

Boundary Value Problems 2013, 2013:222  doi:10.1186/1687-2770-2013-222

 Received: 19 June 2013 Accepted: 4 September 2013 Published: 7 November 2013

This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

### Abstract

In this paper, we investigate an initial boundary value problem for the one-dimensional linear model of thermodiffusion with second sound in a bounded region. Using the semigroup approach, boundary control and the multiplier method, we obtain the existence of global solutions and the uniform decay estimates for the energy.

MSC: 35B40, 35M13, 35Q79.

##### Keywords:
thermodiffusion; second sound; global existence; exponential decay

### 1 Introduction

In this paper, we investigate the global existence and uniform decay rate of the energy for solutions for the one-dimensional model of thermodiffusion with second sound:

ρ u t t ( λ + 2 μ ) u x x + γ 1 θ 1 x + γ 2 θ 2 x = 0 in  ( 0 , 1 ) × ( 0 , + ) , (1.1)

c θ 1 t + k q 1 x + γ 1 u t x + d θ 2 t = 0 in  ( 0 , 1 ) × ( 0 , + ) , (1.2)

n θ 2 t + D q 2 x + γ 2 u t x + d θ 1 t = 0 in  ( 0 , 1 ) × ( 0 , + ) , (1.3)

τ 1 q 1 t + q 1 + k θ 1 x = 0 in  ( 0 , 1 ) × ( 0 , + ) , (1.4)

τ 2 q 2 t + q 2 + D θ 2 x = 0 in  ( 0 , 1 ) × ( 0 , + ) , (1.5)

together with the initial conditions

{ u ( x , 0 ) = u 0 ( x ) , u t ( x , 0 ) = u 1 ( x ) , θ 1 ( x , 0 ) = θ 1 0 ( x ) , θ 2 ( x , 0 ) = θ 2 0 ( x ) , q 1 ( x , 0 ) = q 1 0 ( x ) , q 2 ( x , 0 ) = q 2 0 ( x ) , (1.6)

and the boundary conditions

u ( x , t ) | x = 0 , 1 = θ 1 ( x , t ) | x = 0 , 1 = θ 2 ( x , t ) | x = 0 , 1 = 0 , (1.7)

where u, θ 1 , and q 1 are the displacement, temperature, and heat flux, θ 2 , and q 2 are the chemical potentials and the associated flux. The boundary conditions (1.7) model a rigidly clamped medium with temperature and chemical potentials held constant on the boundary.

Here, we denote by λ, μ the material constants, ρ the density, γ 1 , γ 2 the coefficients of thermal and diffusion dilatation, k, D the coefficients of thermal conductivity, n, c, d the coefficients of thermodiffusion, and τ 1 , τ 2 the (in general very small) relaxation time. All the coefficients above are positive constants and satisfy the condition

n c d 2 > 0 . (1.8)

The classical thermodiffusion equations were first given by Nowacki [1,2] in 1971. The equations describe the process of thermodiffusion in a solid body (see, e.g., [1-5]):

{ ρ u t t ( λ + 2 μ ) u x x + γ 1 θ 1 x + γ 2 θ 2 x = 0 , c θ 1 t k θ 1 x x + γ 1 u t x + d θ 2 t = 0 , n θ 2 t D θ 2 x x + γ 2 u t x + d θ 1 t = 0 . (1.9)

There are many results about the classical thermodiffusion equations. By the method of integral transformations and integral equations, Nowacki [2], Podstrigach [6] and Fichera [7] investigated the initial boundary value problem for the linear homogeneous system. Gawinecki [8] proved the existence, uniqueness and regularity of solutions to an initial boundary value problem for the linear system of thermodiffusion in a solid body. Szymaniec [5] proved the L p - L q time decay estimates along the conjugate line for the solutions of the linear thermodiffusion system. Using the results from [5], Szymaniec [9] obtained the global existence and uniqueness of small data solutions to the Cauchy problem of nonlinear thermodiffusion equations in a solid body. Using the semigroup approach and the multiplier method, Qin et al.[4] obtained the global existence and exponential stability of solutions for homogeneous, nonhomogeneous and semilinear thermodiffusion equations subject to various boundary conditions. Liu and Reissig [3] studied the Cauchy problem for one-dimensional models of thermodiffusion and explained qualitative properties of solutions and showed which part of the model has a dominant influence on wellposedness, propagation of singularities, L p - L q decay estimates on the conjugate line and the diffusion phenomenon.

If we neglect the diffusion in (1.9), then we obtain the classical thermo-elasticity equations. Today models of type I (classical model of thermo-elasticity), of type II (thermal wave), of type III (visco-elastic damping) or second sound present some classification of models of thermo-elasticity (see, e.g., [3,10,11]). By considerations of the total energy equation and comparisons with the models of classical thermo-elasticity and thermodiffusion, we shall propose the linear one-dimensional model of thermodiffusion with second sound as mentioned above. Due to our knowledge, there exist no results for thermodiffusion models with second sound.

Our paper is organized as follows. In Section 2, we present some notations and the main result. Section 3 is devoted to the proof of the main result.

### 2 Notations and main result

Let Ω = ( 0 , 1 ) and

L 2 ( Ω ) = { y L 2 ( Ω ) | 0 1 y ( x ) d x = 0 } , H 1 ( Ω ) = { y H 1 ( Ω ) | 0 1 y ( x ) d x = 0 } . (2.1)

The associated first-order and second-order energy is defined by

E 1 ( t ) E 1 ( t ; u x , u t , θ 1 , θ 2 , q 1 , q 2 ) : = 1 2 0 1 ( ( λ + 2 μ ) u x 2 + ρ u t 2 + c θ 1 2 + n θ 2 2 + τ 1 q 1 2 + τ 2 q 2 2 + 2 d θ 1 θ 2 ) d x , (2.2)

E 2 ( t ) : = E 1 ( t ; u x t , u t t , θ 1 t , θ 2 t , q 1 t , q 2 t ) . (2.3)

The energy E ( t ) is defined by

E ( t ) : = E 1 ( t ) + E 2 ( t ) . (2.4)

Our main result reads as follows.

Theorem 2.1Assume that u 0 H 2 ( Ω ) H 0 1 ( Ω ) , u 1 H 0 1 ( Ω ) , θ 1 0 , θ 2 0 H 2 ( Ω ) H 0 1 ( Ω ) , q 1 0 , q 2 0 H 2 ( Ω ) H 1 ( Ω ) . Then problem (1.1)-(1.7) has a unique global solution such that

u ( t ) C 2 ( [ 0 , + ) , L 2 ( Ω ) ) C 1 ( [ 0 , + ) , H 0 1 ( Ω ) ) C ( [ 0 , + ) , H 2 ( Ω ) H 0 1 ( Ω ) ) , θ 1 ( t ) , θ 2 ( t ) C 1 ( [ 0 , + ) , L 2 ( Ω ) ) C ( [ 0 , + ) , H 2 ( Ω ) H 0 1 ( Ω ) ) , q 1 ( t ) , q 2 ( t ) C 1 ( [ 0 , + ) , H 1 ( Ω ) ) C ( [ 0 , + ) , H 2 ( Ω ) H 1 ( Ω ) ) .

Moreover, the associated energy E ( t ) defined by (2.4) decays exponentially, i.e., there exist positive constants c 0 and C 0 such that

E ( t ) C 0 e c 0 t E ( 0 ) , t > 0 . (2.5)

Remark 2.1 If the initial value ( u 0 , u 1 , θ 1 0 , θ 2 0 , q 1 0 , q 2 0 ) D ( A n ) , n N , ( D ( A n ) will be defined later), then the solution ( u , u t , θ 1 , θ 2 , q 1 , q 2 ) C ( [ 0 , + ) , D ( A n ) ) , and problem (1.1)-(1.7) yields higher regularity in t.

### 3 Proof of the main result

We shall divide the proof into two steps: in Step 1, we shall use the semigroup approach to prove the existence of global solutions and the Remark 2.1; Step 2 is devoted to proving the uniform decay of the energy by the boundary control and the multiplier method.

Step 1. Existence of global solutions.

The proof is based on the semigroup approach (see [4,12]) that can be used to reduce problem (1.1)-(1.7) to an abstract initial value problem for a first-order evolution equation. In order to choose proper space for (1.1)-(1.7), we shall consider the static system associated with them (see [4]). Considering the energy and the property of operator A, we can choose the following state space and the domain of operator A for problem (1.1)-(1.7):

H = H 0 1 ( Ω ) × ( L 2 ( Ω ) ) 3 × ( H 1 ( Ω ) ) 2 , (3.1)

D ( A ) = ( H 2 ( Ω ) H 0 1 ( Ω ) ) × H 0 1 ( Ω ) × ( H 2 ( Ω ) H 0 1 ( Ω ) ) 2 × ( H 2 ( Ω ) H 1 ( Ω ) ) 2 . (3.2)

Using the same method as in [4,12], we can prove that the operator A generates a C 0 -semigroup of contractions on the Hilbert space H. Define

D ( A k ) = { y | y D ( A k 1 ) , A u D ( A k 1 ) , k N } .

Then by Theorem 2.3.1 of [13] about the existence and regularities of solutions, we can complete the proof.

Step 2. Uniform decay of the energy.

In this section, we shall assume the existence of solutions in the Sobolev spaces that we need for our computations. The proof of uniform decay is difficult. It is necessary to construct a suitable Lyapunov function and to combine various techniques from energy method, multiplier approaches and boundary control (see [10,11]). We mainly refer to Racke [11] for the approaches of thermo-elastic models with second sound.

Multiplying (1.1) by u t , (1.2) by θ 1 , (1.3) by θ 2 , (1.4) by q 1 , and (1.5) by q 2 in L 2 , respectively, and summing up the results, yields

d d t E 1 ( t ) = 0 1 ( q 1 2 + q 2 2 ) d x . (3.3)

Similarly, we can get

d d t E 2 ( t ) = 0 1 ( q 1 t 2 + q 2 t 2 ) d x . (3.4)

Multiplying (1.1) by 1 λ + 2 μ u x x in L 2 , we get

2 3 0 1 u x x 2 d x + ρ λ + 2 μ d d t ( 0 1 u t x u x d x ) ρ 0 1 u t x 2 d x λ + 2 μ + 3 0 1 ( γ 1 2 θ 1 x 2 + γ 2 2 θ 2 x 2 ) d x 2 ( λ + 2 μ ) 2 . (3.5)

Multiplying (1.2) by 3 ρ 2 ( λ + 2 μ ) γ 1 u t x , (1.3) by 3 ρ 2 ( λ + 2 μ ) γ 2 u t x in L 2 and summing them up, yields

3 ρ λ + 2 μ 0 1 u t x 2 d x d d t G 1 ( t ) + 1 6 0 1 u x x 2 d x + h 1 0 1 θ 1 x 2 d x + h 2 0 1 θ 2 x 2 d x + h 3 0 1 θ 1 x θ 2 x d x + 81 ρ 2 k 4 ( λ + 2 μ ) 2 γ 1 2 0 1 q 1 t 2 d x + 81 ρ 2 D 4 ( λ + 2 μ ) 2 γ 2 2 0 1 q 2 t 2 d x 3 ρ k 2 ( λ + 2 μ ) γ 1 [ q 1 u t x ] x = 0 x = 1 3 ρ D 2 ( λ + 2 μ ) γ 2 [ q 2 u t x ] x = 0 x = 1 , (3.6)

where

G 1 ( t ) : = 3 ρ k 2 ( λ + 2 μ ) 2 γ 1 0 1 ( ρ k γ 1 q 1 u t t + ρ D γ 2 q 2 u t t τ 1 q 1 q 1 t q 1 2 q 2 2 ρ k + ρ D k D q 1 q 2 τ 2 q 2 q 2 t k τ 2 D q 1 q 2 t D τ 1 k q 2 q 1 t ) d x 3 ρ 2 ( λ + 2 μ ) γ 1 γ 2 × 0 1 [ c γ 2 + d γ 1 k ( q 1 u t + τ 1 q 1 t u t ) + d γ 2 + n γ 1 D ( q 2 u t + τ 2 q 2 t u t ) ] d x , h 3 : = 3 λ + 2 μ ( c γ 2 2 γ 1 + n γ 1 2 γ 2 + d ) , h 1 : = 3 2 ( c λ + 2 μ + d γ 1 γ 2 ) + 81 4 ( c 2 γ 1 2 + d 2 γ 2 2 ) , h 2 : = 3 2 ( n λ + 2 μ + d γ 2 γ 1 ) + 81 4 ( d 2 γ 1 2 + n 2 γ 2 2 ) .

Combining (3.5) with (3.6), we get

2 ρ λ + 2 μ 0 1 u t x 2 d x + 1 2 0 1 u x x 2 d x + d d t ( G 1 ( t ) + ρ λ + 2 μ 0 1 u t x u x d x ) 81 ρ 2 k 4 ( λ + 2 μ ) 2 γ 1 2 0 1 q 1 t 2 d x + ( h 1 + h 3 2 + 3 γ 1 2 2 ( λ + 2 μ ) 2 ) 0 1 θ 1 x 2 d x + 81 ρ 2 D 4 ( λ + 2 μ ) 2 γ 2 2 0 1 q 2 t 2 d x + ( h 2 + h 3 2 + 3 γ 2 2 2 ( λ + 2 μ ) 2 ) 0 1 θ 2 x 2 d x 3 ρ k 2 ( λ + 2 μ ) γ 1 [ q 1 u t x ] x = 0 x = 1 3 ρ D 2 ( λ + 2 μ ) γ 2 [ q 2 u t x ] x = 0 x = 1 . (3.7)

Now, we conclude from (1.1), (1.4), (1.5) and Poincaré inequality

0 1 ( u t 2 + u t t 2 + θ 1 2 + θ 2 2 ) d x 2 ( λ + 2 μ ) 2 ρ 2 0 1 u x x 2 d x + 2 D ( 2 γ 2 2 ρ 2 + 1 π 2 ) 0 1 ( q 2 2 + τ 2 2 q 2 t 2 ) d x + 1 π 2 0 1 u t x 2 d x + 2 k ( 2 γ 1 2 ρ 2 + 1 π 2 ) 0 1 ( q 1 2 + τ 1 2 q 1 t 2 ) d x . (3.8)

Multiplying (1.1) by u in L 2 , we obtain

λ + 2 μ 2 0 1 u x 2 d x 3 2 π 2 ( λ + 2 μ ) 0 1 ( ρ 2 u t t 2 + γ 1 2 θ 1 x 2 + γ 2 2 θ 2 x 2 ) d x . (3.9)

From (1.2) and (1.3), we get

θ 1 t = 1 n c d 2 ( n k q 1 x d D q 2 x + ( n γ 1 d γ 2 ) u t x ) , (3.10)

θ 2 t = 1 n c d 2 ( c D q 2 x d k q 1 x + ( c γ 2 d γ 1 ) u t x ) . (3.11)

Multiplying (3.10) by θ 1 t , (3.11) by θ 2 t in L 2 , and summing up the results, we get

0 1 ( θ 1 t 2 + θ 2 t 2 ) d x d d t G 2 ( t ) 2 k ( n 2 + d 2 ) ( n c d 2 ) 2 0 1 q 1 t 2 d x + ( n γ 1 d γ 2 ) 2 + ( c γ 2 d γ 1 ) 2 ( n c d 2 ) 2 0 1 u t x 2 d x + 0 1 ( θ 1 x 2 + θ 2 x 2 ) d x + 2 D ( c 2 + d 2 ) ( n c d 2 ) 2 0 1 q 2 t 2 d x , (3.12)

where G 2 ( t ) : = 2 n c d 2 0 1 ( n k q 1 θ 1 x d D q 2 θ 1 x + c D q 2 θ 2 x d k q 1 θ 2 x ) d x .

The boundary terms are estimated as follows.

| [ 3 ρ k 2 ( λ + 2 μ ) γ 1 q 1 u t x ] x = 0 x = 1 | 9 ρ 2 k ( | q 1 ( 1 ) | 2 + | q 1 ( 0 ) | 2 ) 8 ( λ + 2 μ ) 2 γ 1 2 ε ˆ + ε ˆ ( | u t x ( 1 ) | 2 + | u t x ( 0 ) | 2 ) 2 , (3.13)

| [ 3 ρ D 2 ( λ + 2 μ ) γ 2 q 2 u t x ] x = 0 x = 1 | 9 ρ 2 D ( | q 2 ( 1 ) | 2 + | q 2 ( 0 ) | 2 ) 8 ( λ + 2 μ ) 2 γ 2 2 ε ˆ + ε ˆ ( | u t x ( 1 ) | 2 + | u t x ( 0 ) | 2 ) 2 , (3.14)

for some ε ˆ > 0 ,

| q 1 ( 1 ) | 2 + | q 1 ( 0 ) | 2 2 ( 1 + 1 ε ˆ 2 ) 0 1 q 1 2 d x + 4 ε ˆ 2 k 0 1 c 2 θ 1 t 2 + d 2 θ 2 t 2 + γ 1 2 u t x 2 d x , (3.15)

| q 2 ( 1 ) | 2 + | q 2 ( 0 ) | 2 2 ( 1 + 1 ε ˆ 2 ) 0 1 q 2 2 d x + 4 ε ˆ 2 D 0 1 d 2 θ 1 t 2 + n 2 θ 2 t 2 + γ 2 2 u t x 2 d x . (3.16)

Combining (3.13)-(3.16), we get

| [ 3 ρ k 2 ( λ + 2 μ ) γ 1 q 1 u t x ] x = 0 x = 1 | + | [ 3 ρ D 2 ( λ + 2 μ ) γ 2 q 2 u t x ] x = 0 x = 1 | 9 ρ 2 ( 1 + ε ˆ 2 ) 4 ( λ + 2 μ ) 2 ε ˆ 3 0 1 ( k γ 1 2 q 1 2 + D γ 2 2 q 2 2 ) d x + ε ˆ ( | u t x ( 1 ) | 2 + | u t x ( 0 ) | 2 ) + 9 ρ 2 ε ˆ 2 ( λ + 2 μ ) 2 0 1 [ ( c 2 γ 1 2 + d 2 γ 2 2 ) θ 1 t 2 + ( d 2 γ 1 2 + n 2 γ 2 2 ) θ 2 t 2 + 2 u t x 2 ] d x . (3.17)

Differentiating (1.1) with respect to t and multiplying by φ ( x ) u t x , where

φ ( x ) : = 1 2 x , (3.18)

we obtain

d d t ( 0 1 ρ u t t φ u t x d x ) ρ 0 1 u t t 2 d x + ( λ + 2 μ ) 2 ( u t x 2 ( 1 ) + u t x 2 ( 0 ) ) ( λ + 2 μ ) 0 1 u t x 2 d x + γ 1 0 1 θ 1 t x φ u t x d x + γ 2 0 1 θ 2 t x φ u t x d x = 0 . (3.19)

Multiplications of (1.2) by φ θ 1 t x and (1.3) by φ θ 2 t x yield

c 0 1 θ 1 t 2 d x k d d t ( 0 1 q 1 x φ θ 1 x d x ) + k 0 1 q 1 t x φ θ 1 x d x = γ 1 0 1 u t x φ θ 1 t x d x + d 0 1 θ 2 t φ θ 1 t x d x , n 0 1 θ 2 t 2 d x D d d t ( 0 1 q 2 x φ θ 2 x d x ) + D 0 1 q 2 t x φ θ 2 x d x = γ 2 0 1 u t x φ θ 2 t x d x + d 0 1 θ 1 t φ θ 2 t x d x ,

which implies, using (1.4) and (1.5),

c 0 1 θ 1 t 2 d x k d d t ( 0 1 q 1 x φ θ 1 x d x ) k τ 1 0 1 q 1 x φ θ 1 x d x + k 2 τ 1 ( θ 1 x 2 ( 1 ) + θ 1 x 2 ( 0 ) ) = k τ 1 0 1 θ 1 x 2 d x + γ 1 0 1 u t x φ θ 1 t x d x + d 0 1 θ 2 t φ θ 1 t x d x , (3.20)

n 0 1 θ 2 t 2 d x D d d t ( 0 1 q 2 x φ θ 2 x d x ) D τ 1 0 1 q 2 x φ θ 2 x d x + D 2 τ 2 ( θ 2 x 2 ( 1 ) + θ 2 x 2 ( 0 ) ) = D τ 2 0 1 θ 2 x 2 d x + γ 2 0 1 u t x φ θ 2 t x d x + d 0 1 θ 1 t φ θ 2 t x d x = 0 . (3.21)

Combining (3.19)-(3.21), we conclude

d d t ( 0 1 ρ u t t φ u t x k q 1 x φ θ 1 x D q 2 x φ θ 2 x d x ) + λ + 2 μ 2 ( u t x 2 ( 1 ) + u t x 2 ( 0 ) ) ρ 0 1 u t t 2 d x + ( λ + 2 μ ) 0 1 u t x 2 d x + c 0 1 θ 1 t 2 d x + n 0 1 θ 2 t 2 d x + 0 1 q 1 x 2 d x + k τ 1 0 1 θ 1 x 2 d x + D τ 2 0 1 θ 2 x 2 d x + k 4 τ 1 2 0 1 θ 1 x 2 d x + D 4 τ 2 2 0 1 θ 2 x 2 d x + 0 1 q 2 x 2 d x . (3.22)

Using (1.2) and (1.3), we get

d d t ( 2 ε ˆ λ + 2 μ 0 1 ρ u t t φ u t x k q 1 x φ θ 1 x D q 2 x φ θ 2 x d x ) + ε ˆ ( u t x 2 ( 1 ) + u t x 2 ( 0 ) ) ( 2 ε ˆ + 4 ε ˆ γ 1 2 ( λ + 2 μ ) k + 4 ε ˆ γ 2 2 ( λ + 2 μ ) D ) 0 1 u t x 2 d x + ( 2 ε ˆ c λ + 2 μ + 4 ε ˆ c 2 ( λ + 2 μ ) k + 4 ε ˆ d 2 ( λ + 2 μ ) D ) 0 1 θ 1 t 2 d x + 2 ε ˆ ρ λ + 2 μ 0 1 u t t 2 d x + ( 2 ε ˆ n λ + 2 μ + 4 ε ˆ d 2 ( λ + 2 μ ) k + 4 ε ˆ n 2 ( λ + 2 μ ) D ) 0 1 θ 2 t 2 d x + ( 2 ε ˆ k ( λ + 2 μ ) τ 1 + ε ˆ k 2 ( λ + 2 μ ) τ 1 2 ) 0 1 θ 1 x 2 d x + ( 2 ε ˆ D ( λ + 2 μ ) τ 2 + ε ˆ D 2 ( λ + 2 μ ) τ 2 2 ) 0 1 θ 2 x 2 d x . (3.23)

With (3.17) and (3.23), we can estimate

| [ 3 ρ k 2 ( λ + 2 μ ) γ 1 q 1 u t x ] x = 0 x = 1 | + | [ 3 ρ D 2 ( λ + 2 μ ) γ 2 q 2 u t x ] x = 0 x = 1 | 9 ρ 2 ( 1 + ε ˆ 2 ) 4 ( λ + 2 μ ) 2 ε ˆ 3 0 1 ( k γ 1 2 q 1 2 + D γ 2 2 q 2 2 ) d x + 2 ε ˆ ρ λ + 2 μ 0 1 u t t 2 d x d d t G 3 ( t ) + ε ˆ ( 9 ρ 2 2 ( λ + 2 μ ) 2 ( c 2 γ 1 2 + d 2 γ 2 2 ) + 2 c λ + 2 μ + 4 c 2 ( λ + 2 μ ) k + 4 d 2 ( λ + 2 μ ) D ) 0 1 θ 1 t 2 d x + ε ˆ ( 9 ρ 2 2 ( λ + 2 μ ) 2 ( d 2 γ 1 2 + n 2 γ 2 2 ) + 2 n λ + 2 μ + 4 d 2 ( λ + 2 μ ) k + 4 n 2 ( λ + 2 μ ) D ) 0 1 θ 2 t 2 d x + ε ˆ ( 2 k ( λ + 2 μ ) τ 1 + k 2 ( λ + 2 μ ) τ 1 2 ) 0 1 θ 1 x 2 d x + ε ˆ ( 2 D ( λ + 2 μ ) τ 2 + D 2 ( λ + 2 μ ) τ 2 2 ) 0 1 θ 2 x 2 d x + ε ˆ ( 9 ρ 2 ( λ + 2 μ ) 2 + 2 + 4 γ 1 2 ( λ + 2 μ ) k + 4 γ 2 2 ( λ + 2 μ ) D ) 0 1 u t x 2 d x , (3.24)

where G 3 : = 2 ε ˆ λ + 2 μ 0 1 ρ u t t φ u t x k q 1 x φ θ 1 x D q 2 x φ θ 2 x d x .

Define ξ : = ρ ( n c d 2 ) 2 4 ( λ + 2 μ ) [ ( n γ 1 d γ 2 ) 2 + ( c γ 2 d γ 1 ) 2 ] . Multiplying both sides of (3.12) by ξ and combining the result with (3.7) and (3.24), we obtain for sufficiently small ε ˆ the estimate

ρ 0 1 u t x 2 d x λ + 2 μ + 1 4 0 1 u x x 2 d x + ξ 2 0 1 ( θ 1 t 2 + θ 2 t 2 ) d x + d G ( t ) d t C 1 0 1 ( q 1 2 + q 2 2 + q 1 t 2 + q 2 t 2 ) d x , (3.25)

where G ( t ) = G 1 ( t ) + ρ λ + 2 μ 0 1 u t x u x d x ξ G 2 ( t ) + G 3 ( t ) , and C 1 = C 1 ( ε ˆ , λ , μ , γ 1 , γ 2 , k , D , n , c , d , τ 1 , τ 2 ) . Now, we can define the desired Lyapunov functional F ( t ) . For ε > 0 , to be determined later on, let

F ( t ) : = 1 ε E ( t ) + G ( t ) . (3.26)

Then we conclude from (3.3), (3.4), and (3.15)

d d t F ( t ) 1 ε 0 1 ( q 1 2 + q 2 2 + q 1 t 2 + q 2 t 2 ) d x ρ λ + 2 μ 0 1 u t x 2 d x 1 4 0 1 u x x 2 d x ξ 2 0 1 ( θ 1 t 2 + θ 2 t 2 ) d x + C 1 0 1 ( q 1 2 + q 2 2 + q 1 t 2 + q 2 t 2 ) d x . (3.27)

By using (3.8), we arrive at

C 2 ε 0 1 ( q 1 2 + q 2 2 + q 1 t 2 + q 2 t 2 ) d x ε π 2 0 1 u t x 2 d x 2 ε ( λ + 2 μ ) 2 ρ 2 0 1 u x x 2 d x ε 0 1 ( u t 2 + u t t 2 + θ 1 2 + θ 2 2 ) d x , (3.28)

while (3.9) yields

C 3 ε 2 0 1 u x x 2 d x C 4 ε 2 0 1 ( q 1 2 + q 2 2 + q 1 t 2 + q 2 t 2 ) d x ε 2 λ + 2 μ 2 0 1 u x 2 d x . (3.29)

Combining (3.27)-(3.29), we conclude

d d t F ( t ) ( 1 ε C 1 C 2 ε C 4 ε 2 ) 0 1 ( q 1 2 + q 2 2 + q 1 t 2 + q 2 t 2 ) d x ε 2 λ + 2 μ 2 0 1 u x 2 d x ( 2 ρ λ + 2 μ ε π 2 ) 0 1 u t x 2 d x [ 1 4 2 ε ( λ + 2 μ ) 2 ρ 2 C 3 ε 2 ] 0 1 u x x 2 d x ξ 2 0 1 ( θ 1 t 2 + θ 2 t 2 ) d x ε 0 1 ( u t 2 + u t t 2 + θ 1 2 + θ 2 2 ) d x . (3.30)

We choose 0 < ε 1 such that all terms on the right-hand side of (3.30) become negative,

i.e. , ε ε 1 : = min { ρ π 2 λ + 2 μ , ρ 2 32 ( λ + 2 μ ) 2 , 1 4 1 C 3 , 1 2 ( C 1 + C 2 + C 4 ) } . (3.31)

Choosing ε as in (3.31), we obtain from (3.30)

d d t F ( t ) d 1 0 1 ( u x 2 + u t 2 + u t x 2 + u t t 2 + θ 1 2 + θ 2 2 + θ 1 t 2 + θ 2 t 2 + q 1 2 + q 2 2 + q 1 t 2 + q 2 t 2 ) d x

with

d 1 : = min { ε 2 , 1 8 , ρ λ + 2 μ , ( λ + 2 μ ) ε 2 4 , ξ 4 , 1 2 ε } , (3.32)

which implies

d d t F ( t ) d 2 E ( t ) , d 2 : = d 1 2 min { 1 ρ , 1 λ + 2 μ , 1 c , 1 n , 1 τ 1 , 1 τ 2 } . (3.33)

There exist positive constants C 5 , C 6 and ε 2 such that for any ε ε 2 and t 0 , it holds

C 5 E ( t ) F ( t ) C 6 E ( t ) , (3.34)

where C 5 , C 6 are determined later on. In fact,

| G ( t ) | C 5 E ( t ) (3.35)

with C 5 : = max { 2 i 1 τ 1 , 2 i 2 τ 2 , 2 i 3 τ 1 , 2 i 4 τ 2 , 2 i 5 ρ , 2 i 6 λ + 2 μ , 2 i 7 ρ , 2 i 8 c , 2 i 9 n , ρ ( λ + 2 μ ) 2 } ,

i 1 : = 3 ρ 4 ( λ + 2 μ ) 2 [ ρ k γ 1 + ( τ 1 + 2 ) + k + D k D + k τ 2 D ] + 3 ρ ( c γ 2 + d γ 1 ) 4 ( λ + 2 μ ) γ 1 γ 2 k + ( n + d ) k n c d 2 + [ n k + d D n c d 2 + ε ˆ k λ + 2 μ ] 2 k , i 2 : = 3 ρ 4 ( λ + 2 μ ) 2 [ ρ D γ 2 + ( τ 2 + 2 ) + k + D k D + D τ 1 k ] + 3 ρ ( n γ 1 + d γ 2 ) 4 ( λ + 2 μ ) γ 1 γ 2 D + ( c + d ) D n c d 2 + [ c D + d k n c d 2 + ε ˆ D λ + 2 μ ] 2 D , i 3 : = 3 ρ τ 1 4 ( λ + 2 μ ) 2 ( 1 + D k ) + 3 ρ τ 1 ( c γ 2 + d γ 1 ) 4 ( λ + 2 μ ) γ 1 γ 2 k + [ n k + d D n c d 2 + ε ˆ k λ + 2 μ ] 2 τ 1 2 k , i 4 : = 3 ρ τ 2 4 ( λ + 2 μ ) 2 ( 1 + k D ) + 3 ρ τ 2 ( d γ 2 + n γ 1 ) 4 ( λ + 2 μ ) γ 1 γ 2 D + [ c D + d k n c d 2 + ε ˆ D λ + 2 μ ] 2 τ 2 2 D , i 5 : = 3 ρ ( 1 + τ 1 ) ( d γ 1 + c γ 2 ) 4 ( λ + 2 μ ) γ 1 γ 2 k + 3 ρ ( 1 + τ 2 ) ( d γ 2 + n γ 1 ) 4 ( λ + 2 μ ) γ 1 γ 2 D , i 6 : = 1 λ + 2 μ [ ρ ( 1 + 2 ε ˆ ) 2 + 2 ε ˆ γ 1 2 k + 2 ε ˆ γ 2 2 D ] , i 7 = ρ λ + 2 μ [ 3 k 4 ( λ + 2 μ ) γ 1 + 3 D 4 ( λ + 2 μ ) γ 2 + ε ˆ ] , i 8 : = 2 ε ˆ c 2 ( λ + 2 μ ) k + 2 ε ˆ d 2 ( λ + 2 μ ) D , i 9 = 2 ε ˆ d 2 ( λ + 2 μ ) k + 2 ε ˆ n 2 ( λ + 2 μ ) D .

At this point, we choose ε ε 2 : = 1 2 C 5 , C 6 : = 1 ε + C 5 . Finally, we choose ε : = min { ε 1 , ε 2 } . Thus, we have the validity of (3.34). Combining (3.33) with (3.34), we get

d d t F ( t ) d 0 F ( t ) , d 0 : = d 2 C 6 . (3.36)

Hence, it follows from (3.36), F ( t ) e d 0 t F ( 0 ) . Applying (3.34) again, we can conclude (2.5) with C 0 : = C 6 C 5 . The proof is complete.

### Competing interests

The authors declare that they have no competing interests.

### Authors’ contributions

The paper is a joint work of all authors who contributed equally to the final version of the paper. All authors read and approved the final manuscript.

### Acknowledgements

This paper was in part supported by the NNSF of China with contract numbers 11031003, 11271066 and a grant from Shanghai Education Commission 13ZZ048.

### References

1. Nowacki, W: Certain problem of thermodiffusion. Arch. Mech.. 23, 731–754 (1971)

2. Nowacki, W: Dynamical problem of thermodiffusion in solids II. Bull. Pol. Acad. Sci., Tech. Sci.. 22, 205–211 (1974)

3. Liu, Y, Reissig, M: Models of thermodiffusion in 1D. Math. Methods Appl. Sci. (in press) (2013)

(in press)

Publisher Full Text

4. Qin, Y, Zhang, M, Feng, B, Li, H: Global existence and asymptotic behavior of solutions for thermodiffusion equations. J. Math. Anal. Appl.. 408, 140–153 (2013). Publisher Full Text

5. Szymaniec, A: L p - L q time decay estimates for the solution of the linear partial differential equations of thermodiffusion. Appl. Math.. 37, 143–170 (2010)

6. Podstrigach, Y: Differential equations of the problem of thermodiffusion in an isotropic deformable solid. Dopov Akad. Nauk Ukr. RSR. 2, 169–172 (1961)

7. Fichera, G: Uniqueness, existence and estimate of the solution in the dynamical problem of thermodiffusion in an elastic solid. Arch. Ration. Mech. Anal.. 26, 903–920 (1974)

8. Gawinecki, JA, Sierpinski, K: Existence, uniqueness and regularity of the solution of the first boundary initial value problem for the equation of the thermodiffusion in a solid body. Bull. Pol. Acad. Sci., Tech. Sci.. 30, 541–547 (1982)

9. Szymaniec, A: Global solution to the Cauchy problem of nonlinear thermodiffusion in a solid body. Appl. Math.. 37, 437–458 (2010)

10. Jiang, S, Racke, R: Evolution Equations in Thermoelasticity, Chapman & Hall/CRC, Boca Raton (2000)

11. Racke, R: Thermoelasticity with second sound-exponential stability in linear and non-linear 1-d. Math. Methods Appl. Sci.. 25, 409–441 (2002). Publisher Full Text

12. Liu, Z, Zheng, S: Semigroups Associated with Dissipative Systems, Chapman & Hall/CRC, Boca Raton (1999)

13. Zheng, S: Nonlinear Evolution Equations, CRC Press, Boca Raton (2004)