Next Article in Journal
Optimal Incentive Contract for Sales Team with Loss Aversion Preference
Next Article in Special Issue
Hysteretically Symmetrical Evolution of Elastomers-Based Vibration Isolators within α-Fractional Nonlinear Computational Dynamics
Previous Article in Journal
Transformation Properties under the Operations of the Molecular Symmetry Groups G36 and G36(EM) of Ethane H3CCH3
Previous Article in Special Issue
Two-Dimensional Finite Element in General Plane Motion Used in the Analysis of Multi-Body Systems
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

On the Partition of Energies for the Backward in Time Problem of Thermoelastic Materials with a Dipolar Structure

1
Department of Mathematics and Computer Science, Transilvania University of Brasov, 500093 Brasov, Romania
2
Department of Mechanical Engineering, Transilvania University of Brasov, 500093 Brasov, Romania
3
Center for Modeling & Computer Simulation, Research Institute, King Fahd University of Petroleum & Minerals, Dhahran 31261, Saudi Arabia
4
Department of Mathematics & Statistics, Faculty of Basic and Applied Sciences (FBAS), International Islamic University (IIUI), Islamabad 44000, Pakistan
5
College of Mathematics and Systems Science, Shandong University of Science and Technology, Qingdao 266590, Shandong, China
6
Shanghai Institute of Applied Mathematics and Mechanics, Shanghai University, Shanghai 200072, China
*
Author to whom correspondence should be addressed.
Symmetry 2019, 11(7), 863; https://doi.org/10.3390/sym11070863
Submission received: 12 June 2019 / Revised: 22 June 2019 / Accepted: 26 June 2019 / Published: 2 July 2019
(This article belongs to the Special Issue Symmetry in Applied Continuous Mechanics)

Abstract

:
We first formulate the mixed backward in time problem in the context of thermoelasticity for dipolar materials. To prove the consistency of this mixed problem, our first main result is regarding the uniqueness of the solution for this problem. This is obtained based on some auxiliary results, namely, four integral identities. The second main result is regarding the temporal behavior of our thermoelastic body with a dipolar structure. This behavior is studied by means of some relations on a partition of various parts of the energy associated to the solution of the problem.

1. Introduction

In our study, we approach a thermoelastic body having a dipolar structure. This kind of structure falls within a more general theory, namely, the theory of bodies with microstructure. The first studies in this context were published by Eringen (see, for instance, references [1,2]). One may deduce the importance of the dipolar structure due to the large number of published studies dedicated to this topic, of which we can mention [3,4,5,6,7]. As such, our present work can be considered a continuation in this respect.
A continuation of the theories of microstructure, is a theory that takes into account the voids in the materials. It is considered that the initiators of this theory were Nunziato and Cowin, in their known paper [8]. After that, the number of studies within this topic has grown impressively. We want to enumerate some of these [9,10,11,12,13,14,15,16]: The first result for the backward in time problem belongs to Serrin, who approached this problem in the context of Navier–Stokes equations (see [17]). In the paper [17], we find some uniqueness in the results with regards to the forward in time problem. After that, the number of studies dedicated to the backward in time problem has increased considerably. Of particular importance are the works [18,19,20,21,22,23,24,25,26,27]. We have to point out that the results obtained by Ciarletta in [23], and Ciarletta and Chirita in [24] were improved by Quintanilla in [25]. In addition, Quintanilla approached the question of location in time for solutions to the backward in time problem, in the context of thermoelasticity of Green and Naghdi [26,27] and the theory of porous thermoelastic bodies. The elastic porous bodies were also approached by Iovane and Passarella in [28]. The forward in time problem, in the context of theory for thermomicrostretch elastic solids, was approached by Passarella and Tibullo in [29]. It is worth noting that the idea of considering non-standard problems, in the context of the general theory of bodies having a dipolar structure, was inspired by Quintanilla and Straughan’s work [30]. In [31,32] it is proved that the nonhomogenous temperature field has a profound influence on the nanobeam mechanics. Additionally, [33] is a recent contribution on stress-driven nonlocal modeling of thermoelastic nanostructures.
Here is the plane of our study. First of all, we summarize the main equations, the initial conditions, and the boundary data of the mixed problem. Then, we prove some estimates for the gradient of classical and dipolar displacements, and for the gradient of the function of voids. In the last part of our study we prove the main result, namely, the continuous dependence of solutions—with regards to the coefficients that couple the equations describing the dipolar deformation—with the equations that describe the behavior of voids. The description of the continuous dependence was possible due to the definition of an adequate measure.

2. Basic Equations and Conditions

In our paper, we approach a thermoelastic body having a dipolar structure. We will use an anisotropic body, which is situated in a regular domain D, included in the physical space E 3 , that is, the three-dimensional Euclidean space. Consider that the boundary of the domain is a piecewise smooth surface D . The closure of D is usual denoted by D ¯ , D ¯ = D D . An orthonormal system of references is introduced, and then tensors and vectors have components with Latin subscripts over 1, 2, 3. Typical conventions for summation over repeated indices and for derivation operations are implied. So, a subscript preceded by a comma is for a partial derivative with regards to corresponding spatial coordinate; while a superposed dot is for a derivative with regards to time variable. All the functions we use are assumed to be sufficiently regular as necessary. Additionally, if there is no possibility of confusion, then the dependence of function with regards to its spatial or time variables will be omitted. The evolution of the body with a dipolar structure will be described with the help of the following specific variables:
u i ( x , t ) , ϕ i j ( x , t ) , θ ( x , t ) , ( x , t ) D × [ 0 , t 0 ) .
Here, we denoted by u i the components of the displacement vector field, by ϕ i j the components of the dipolar displacement tensor field, and by θ the absolute temperature.
Using the above variables u i ( x , t ) , and ϕ i j ( x , t ) we will introduce the components of the tensors of strain, namely, ε i j , κ i j , and χ i j k , as follows:
2 ε i j = u j , i + u i , j , κ i j = u j , i ϕ i j , χ i j k = ϕ i j , k .
All our considerations are made within a linear theory, therefore it is natural to consider that the Helmholtz’s free energy is a quadratic form with regards to its independent constitutive variables. The Helmholtz’s free energy in the reference configuration will be denoted by W. So, in accordance with the principle of conservation of energy, we develop in series the function W and we keep the terms only until the second order. Because the reference state was assumed to be free of loadings, we deduce that the Helmholtz’s free energy per mass can be considered of the following form (see [30]):
W = 1 2 A i j m n ε i j ε m n + D i j m n ε i j κ m n + F i j m n r ε i j χ m n r + 1 2 B i j m n κ i j κ m n + G i j m n r κ i j χ m n r + 1 2 C i j k m n r χ i j k χ m n r a i j ε i j θ b i j κ i j θ c i j k χ i j k θ 1 2 c θ 2 .
We will use this form of free energy used in the entropy production inequality and deduce the motion equations. In addition, from the same inequality, the constitutive equations are obtained. These equations express the tensors of stress with the help of the tensors of deformation. We will denote the components of the stress measures by τ i j , η i j , and μ i j k . In this way, the constitutive equations establish a connection between the tensors τ i j , η i j , μ i j k and the tensors ε i j , κ i j , χ i j k .
We will use a procedure similar to that used by Green and Rivlin in [6], so that considering the Helmholtz’ free energy (3) we deduce the next constitutive equations
τ i j = W ε i j = A i j m n ε m n + D m n i j κ m n + F m n r i j χ m n r a i j θ , η i j = W κ i j = D i j m n ε m n + B i j m n κ m n + G i j m n r χ m n r b i j θ , μ i j k = W χ i j k = F i j k m n ε m n + G m n i j k κ m n + C i j k m n r χ m n r c i j k θ , η = W θ = a i j ε i j + b i j κ i j + c i j k χ i j k + c θ ,
which are satisfied in D × [ 0 , t 0 ) . Here, we denoted by η the entropy per unit mass.
For the vector of heat flux, having the components q i we have a classical constitutive relation, namely,
q i = K i j θ , j ,
where K i j is the thermal conductivity symmetric tensor.
Also, we can deduce the main equations that govern the thermoelasticity of bodies with a dipolar structure, namely (see [5,6]):
-
the motion equations:
τ i j + η i j , j + ρ f i = ρ u ¨ i , μ i j k , i + η j k + ρ g j k = I k r ϕ ¨ j r ;
-
the equation of energy:
ρ T 0 η ˙ = q i , i + ρ r .
The signification of the notations that we introduced in preceding equations is as follows: ρ , the density of mass, which is a constant; I i j , the symmetric tensor of microinertia; k, the intrinsic inertia; ε i j , κ i j , χ i j k , the strain tensors; τ i j , η i j , μ i j k , the stress tensors; f i , the body forces; g j k , the dipolar charges; A i j m n , B i j m n , . . . , a i j , the functions that describe the properties of the material in terms of elasticity. Suppose the following symmetry relations take place:
A i j m n = A j i m n = A m n i j , B i j m n = B m n i j , a i j = a j i , C i j k m n r = C m n r i j k , F i j k m n = F i j k n m , D i j m n = D i j n m .
Assuming that there are no supply terms and taking into account the constitutive Equations (4) and (5) and the kinematic Equation (2), Equations (6) and (7) become
ρ u ¨ i = C i j m n + G i j m n u n , m + G m n i j + B i j m n u n , m ϕ m n + F m n r i j + D i j m n r ϕ n r , m a i j + b i j θ , j ,
I k r ϕ ¨ j r = F i j k m n u n , m + D m n i j k u n , m ϕ m n + A i j k m n r ϕ n r , m c i j k θ , i + G j k m n u m , n + B j k m n u n , m ϕ m n + D j k m n r ϕ n r , m b j k θ ,
K i j θ , j , i = T 0 a i j u ˙ i , j + b i j u ˙ j , i ϕ ˙ i j + c i j k ϕ ˙ i j , k + c θ ˙ .
From now, we will assume that the Equations (2), (4), and (9) will be satisfied on the interval ( , 0 ] .
The outward unit normal to the surface D has the components n i . With the help of this normal we can define the surface traction’s of components t i , the surface couple of components μ j k , and the flux of heat, q. All of this makes sense in every point of regularity of the boundary D and has the following expressions
t i = τ i j + η i j n j , μ j k = μ i j k n i , q = q i n i .
In close relation to these surface tractions, we consider the following homogeneous boundary conditions:
u i ( x , t ) = 0 , ( x , t ) D u × ( , 0 ] , t i = 0 , ( x , t ) D u c × ( , 0 ] , ϕ i j ( x , t ) = 0 , ( x , t ) D ϕ × ( , 0 ] , m j k = 0 , ( x , t ) D ϕ c × ( , 0 ] , θ ( x , t ) = 0 , ( x , t ) D θ × ( , 0 ] , q = 0 , ( x , t ) D θ c × ( , 0 ] ,
where the surfaces D u , D ϕ , D θ , and its complements D u c , D ϕ c , D θ c are subsurfaces of the border D , which are subject to the following restrictions:
D ¯ u D u c = D ¯ ϕ D ϕ c = D ¯ θ D θ c = D , D u D u c = D ϕ D ϕ c = D θ D θ c = .
We still have to add the final restrictions. So, on the closed domain D ¯ we have:
u i ( x , 0 ) = u i 0 ( x ) , u ˙ i ( x , 0 ) = u i 1 ( x ) , θ ( x , 0 ) = θ 0 ( x ) , ϕ i j ( x , 0 ) = ϕ i j 0 ( x ) , ϕ ˙ i j ( x , 0 ) = ϕ i j 1 ( x ) ,
where u i 0 ( x ) , u i 1 ( x ) , ϕ i j 0 ( x ) , ϕ i j 1 ( x ) , and θ 0 ( x ) are continuous prescribed functions in all points where they are defined. Additionally, these functions are assumed be compatible with conditions (11) on the appropriate subsets of D .
Let us consider the internal energy density Ψ (see [30]), which has the following expression:
Ψ = 1 2 A i j m n ε i j ε m n + D i j m n ε i j κ m n + F i j m n r ε i j χ m n r + 1 2 B i j m n κ i j κ m n + G i j m n r κ i j χ m n r + 1 2 C i j k m n r χ i j k χ m n r .
P is denoted the so-called boundary-final value problem, which consists of Equation (9), the boundary restrictions (11), and the final data (12).
To obtain the results we have proposed, we will have to impose some conditions on the functions we are dealing with.
So, if J m ( x ) is the minimum eigenvalue of the inertia tensor I i j ( x ) , then we need to assume that J m and ρ are continuous functions and the constitutive coefficients are of class C 1 ( D ) . We also assume that:
(a)
ρ ( x ) a 1 , J m ( x ) a 2 , c ( x ) c 0 , where a 1 , a 2 , c 0 are real positive constants;
(b)
the tensor K i j is positive definite;
(c)
the internal energy density Ψ is a positive definite quadratic form.
Based on hypothesis (b), we deduce that there exist two positive numbers, K m and K M , so that
K m θ , i θ , j K i j θ , i θ , j K M θ , i θ , j ,
and, as a consequence of the hypothesis (c), we can find the positive constants M 1 and M 2 so that the next inequality is satisfied:
M 1 2 ε i j ε i j + κ i j κ i j + χ i j k χ i j k Ψ M 2 2 ε i j ε i j + κ i j κ i j + χ i j k χ i j k .
These hypotheses are not considered as very restrictive, as they are commonly imposed in mechanics of continuous media.
It is not difficult to equate our boundary-final value problem P with a boundary-initial problem, denoted by P , by a convenient change of variables. In this regard, we set h ( t ) = h ( t ) , for t = t . But, to simplify writing, we will give up the sign “prime” so that the P problem will be defined by the following conditions and equations:
-
the motion Equations (9a) and (9b), satisfied in D × [ 0 , ) ;
-
the energy equation:
K i j θ , j , i = T 0 a i j u ˙ i , j + b i j u ˙ j , i ϕ ˙ i j + c i j k ϕ ˙ i j , k + c θ ˙ , in D × [ 0 , ) ;
-
the geometric Equation (2), satisfied in D × [ 0 , ) ;
-
the constitutive Equation (4), satisfied in D × [ 0 , ) ;
-
the initial conditions (11), satisfied in D ¯ ;
-
the boundary conditions:
u i ( x , t ) = 0 , ( x , t ) D u × [ 0 , ) , t i = 0 , ( x , t ) D u c × [ 0 , ) , ϕ i j ( x , t ) = 0 , ( x , t ) D ϕ × [ 0 , ) , m j k = 0 , ( x , t ) D ϕ c × [ 0 , ) , θ ( x , t ) = 0 , ( x , t ) D θ × [ 0 , ) , q = 0 , ( x , t ) D θ c × [ 0 , ) .

3. Main Result

We first establish some integral identities regarding a solution u = u i , ϕ i j , θ of the mixed problem P . These will be useful in obtaining the important results of our study.
Proposition 1.
If the ordered array u = u i , ϕ i j , θ satisfies the mixed problem P , then the following equality takes place:
B 1 2 ρ u ˙ i ( t ) u ˙ i ( t ) + I j k ϕ ˙ j m ( t ) ϕ ˙ k m ( t ) + Ψ ( t ) + 1 2 c θ 2 ( t ) d V = B ρ u ˙ i ( 0 ) u ˙ i ( 0 ) + I j k ϕ ˙ j m ( 0 ) ϕ ˙ k m ( 0 ) + Ψ ( 0 ) + 1 2 c θ 2 ( 0 ) d V + 0 t D 1 T 0 K i j θ , i ( τ ) θ , j ( τ ) d V d τ , t [ 0 , ) .
Proof. 
Taking into account the kinematics compatibility relations (2) and the differential conditions of equilibrium (9a) and (9b), we obtain the following equality:
1 2 t ρ u ˙ i ( t ) u ˙ i ( t ) + I j k ϕ ˙ j m ( t ) ϕ ˙ k m ( t ) = τ i j + η i j u ˙ j + μ i j k ϕ ˙ j k , i τ i j ε ˙ i j + η i j κ ˙ i j + μ i j k χ ˙ i j k .
Taking into account the constitutive Equation (4), the symmetry relations (8), and the expression of the internal energy density Ψ from (13), the last parentheses in the right-hand side of (19) becomes
τ i j ε ˙ i j + η i j κ ˙ i j + μ i j k χ ˙ i j k = t Ψ + 1 2 c θ 2 + 1 T 0 q j θ , j 1 T 0 K i j θ , i θ , i .
We substitute Equation (20) into Equation (19), then the resulting equality is integrated on cylinder [ 0 , t ] × D . If we use the theorem of divergence and consider the conditions to the limit (17), we are led to equality (18), such that the proof of the proposition is finished. □
In a similar way, one can demonstrate the identity that follows, as a complement to identity (18):
D 1 2 ρ u ˙ i ( t ) u ˙ i ( t ) + I j k ϕ ˙ j m ( t ) ϕ ˙ k m ( t ) + Ψ ( t ) 1 2 c θ 2 ( t ) d V = D ρ u ˙ i ( 0 ) u ˙ i ( 0 ) + I j k ϕ ˙ j m ( 0 ) ϕ ˙ k m ( 0 ) + Ψ ( 0 ) 1 2 c θ 2 ( 0 ) d V 0 t D u ˙ i ( τ ) a i j + b j i θ ( τ ) , j + ϕ ˙ i j ( τ ) c i j k θ ( τ ) , k b i j ϕ ˙ i j ( τ ) θ ( τ ) + 1 T 0 K i j θ , i ( τ ) θ , j ( τ ) d V d τ , t [ 0 , ) .
To simplify writing, we enter the notation
2 F ( x , y ) = A i j m n ε i j ( x ) ε m n ( y ) + D i j m n ε i j ( x ) κ m n ( y ) + ε i j ( y ) κ m n ( x ) + F i j m n r ε i j ( x ) χ m n r ( y ) + ε i j ( y ) χ m n r ( x ) + B i j m n κ i j ( x ) κ m n ( y ) + G i j m n r κ i j ( x ) χ m n r ( y ) + κ i j ( y ) χ m n r ( x ) + C i j k m n r χ i j k ( x ) χ m n r ( y ) .
By using the symmetry relations (8), from (22) we deduce
F ( x , y ) = F ( y , x ) .
By direct substitution in (22) and taking into account Equation (13), we also obtain
F ( τ , τ ) = Ψ ( τ ) .
Now, we can prove a result similar to (18), but in the case of homogeneous initial conditions.
Proposition 2.
Consider a solution u i , ϕ i j , θ of the problem backward in time P , which corresponds to null initial data. Then, the next identity takes place
B 1 2 ρ u ˙ i ( t ) u ˙ i ( t ) + I j k ϕ ˙ j m ( t ) ϕ ˙ k m ( t ) 1 2 c θ 2 ( t ) d V = B 1 2 A i j m n ε i j ( t ) ε m n ( t ) + D i j m n ε i j ( t ) κ m n ( t ) + F i j m n r ε i j ( t ) χ m n r ( t ) + 1 2 B i j m n κ i j ( t ) κ m n ( t ) + G i j m n r κ i j ( t ) χ m n r ( t ) + 1 2 C i j k m n r χ i j k ( t ) χ m n r ( t ) d V ,
for all t [ 0 , ) .
Proof. 
By direct calculations, for a fixed t ( 0 , ) , we get the identity:
t ρ u ˙ i ( τ ) u ˙ i ( 2 t τ ) + I j k ϕ ˙ j m ( τ ) ϕ ˙ k m ( 2 t τ ) c θ ( τ ) θ ( 2 t τ ) = ρ u ¨ i ( τ ) u ˙ i ( 2 t τ ) + I j k ϕ ¨ j m ( τ ) ϕ ˙ k m ( 2 t τ ) + c θ ( τ ) θ ˙ ( 2 t τ ) ρ u ˙ i ( τ ) u ¨ i ( 2 t τ ) + I j k ϕ ˙ j m ( τ ) ϕ ¨ k m ( 2 t τ ) c θ ˙ ( τ ) θ ( 2 t τ ) .
Taking into account the kinematic Equation (2), the constitutive Equation (4), the motion Equation (9), and the symmetry relations (8), we are led to the equality
t ρ u ˙ i ( τ ) u ˙ i ( 2 t τ ) + I j k ϕ ˙ j m ( τ ) ϕ ˙ k m ( 2 t τ ) c θ ( τ ) θ ( 2 t τ ) = τ i j + η i j ( τ ) u ˙ i ( 2 t τ ) τ i j + η i j ( 2 t τ ) u ˙ j ( τ ) + μ i j k ϕ ¨ j k ( τ ) ϕ ˙ j k ( 2 t τ ) μ i j k ϕ ¨ j k ( 2 t τ ) ϕ ˙ j k ( τ ) 1 T 0 θ ( τ ) q i ( 2 t τ ) + 1 T 0 θ ( 2 t τ ) q i ( τ ) , i + F ( τ , 2 t τ ) ,
the function F ( . , . ) being defined in (22).
We just need to integrate this equality into [ 0 , t ] × D , to keep in mind that the initial data are null and to use the definition (13), so we get the desired equality (25) and the proof of proposition is complete. □
The following two propositions are also useful in establishing the main outcomes of our study.
Proposition 3.
Consider that the ordered array u = u i , ϕ i j , θ satisfies the mixed problem P . Then, the following equality takes place
D ρ u i ( t ) u ˙ i ( t ) + I j k ϕ j m ( t ) ϕ ˙ k m ( t ) 1 2 T 0 K i j 0 t θ ( τ ) d τ , i 0 t θ ( τ ) d τ , j d V = D ρ u i ( 0 ) u ˙ i ( 0 ) + I j k ϕ j m ( 0 ) ϕ ˙ k m ( 0 ) d V + 0 t D ρ η ( 0 ) θ ( τ ) d V d τ + 0 t D ρ u ˙ i ( t ) u ˙ i ( t ) + I j k ϕ ˙ j m ( t ) ϕ ˙ k m ( t ) 2 Ψ ( τ ) c θ 2 ( τ ) d V d τ ,
for all t [ 0 , ) .
Proof. 
We will consider the constitutive relations (4), the geometric Equation (2), and the motion Equations (9a) and (9b), we deduce
t ρ u i ( t ) u ˙ i ( t ) + I j k ϕ j m ( t ) ϕ ˙ k m ( t ) = ρ u ˙ i ( t ) u ˙ i ( t ) + I j k ϕ ˙ j m ( t ) ϕ ˙ k m ( t ) + τ i j ( t ) + η i j ( t ) u j ( t ) + μ i j k ( t ) ϕ j k ( t ) , i τ i j ( t ) + η i j ( t ) ε i j ( t ) + η i j ( t ) κ i j ( t ) + μ i j k ( t ) χ i j k ( t ) .
Taking into account definition (13), the last parentheses from the right-hand side of (29) can be restated in the following form:
τ i j ( t ) ε i j ( t ) + η i j ( t ) κ i j ( t ) + μ i j k ( t ) χ i j k ( t ) = 1 2 c θ 2 ( t ) + Ψ ( t ) + 1 T 0 θ ( t ) 0 t q i ( τ ) d τ , i 1 T 0 K i j 0 t θ ( τ ) d τ , i 0 t θ ( τ ) d τ , i ρ η ( 0 ) θ ( t ) .
Here, the expression of the entropy η is obtained by integrating its equation of evolution (19), with regards to time variable.
We integrate the equality (30) on [ 0 , t ] × D , and take into account the data on the border (17) and the theorem of divergence. As such, we arrive at the proposed equality (29), and the demonstration of proposition is over. □
Proposition 4.
Let us consider a solution u i , ϕ i j , θ of the mixed problem P . Then, the following equality takes place
2 D ρ u i ( t ) u ˙ i ( t ) + I j k ϕ j m ( t ) ϕ ˙ k m ( t ) 1 2 T 0 K i j 0 t θ ( τ ) d τ , i 0 t θ ( τ ) d τ , j d V = D ρ u ˙ i ( 0 ) u i ( 2 t ) + I j k ϕ ˙ j m ( 0 ) ϕ k m ( 2 t ) d V + D ρ u i ( 0 ) u ˙ i ( 2 t ) + I j k ϕ j m ( 0 ) ϕ ˙ k m ( 2 t ) d V 0 t D ρ η ( 0 ) θ ( t + τ ) θ ( t τ ) d V d τ ,
for all t [ 0 , ) .
Proof. 
We will take into account the geometric relations (2) and the motion Equations (9a) and (9b), we deduce the following identity
t ρ u ˙ i ( t + τ ) u i ( t τ ) + u i ( t + τ ) u ˙ i ( t τ ) + t I j k ϕ ˙ j m ( t + τ ) ϕ k m ( t τ ) + ϕ j m ( t + τ ) ϕ ˙ k m ( t τ ) = τ i j + η i j ( t + τ ) u i ( t τ ) τ i j + η i j ( t τ ) u j ( t + τ ) + μ i j k ( t + τ ) ϕ j k ( t τ ) μ i j k ( t τ ) ϕ j k ( t + τ ) , i τ i j + η i j ( t + τ ) ε i j ( t τ ) τ i j + η i j ( t τ ) ε i j ( t + τ ) μ i j k ( t + τ ) χ i j k ( t τ ) μ i j k ( t τ ) χ j k ( t + τ ) .
On the other hand, by using the symmetry Equation (8) and the constitutive relations (4), the last two brackets receive the following form:
τ i j + η i j ( t + τ ) ε i j ( t τ ) τ i j + η i j ( t τ ) ε i j ( t + τ ) + μ i j k ( t + τ ) χ i j k ( t τ ) μ i j k ( t τ ) χ j k ( t + τ ) = 1 T 0 θ ( t + τ ) 0 t τ q i ( s ) d s θ ( t τ ) 0 t + τ q i ( s ) d s , i 1 T 0 K i j 0 t + τ θ ( s ) d s , i 0 t τ θ ( s ) d s , j 0 t + τ θ ( s ) d s , i 0 t τ θ ( s ) d s , j ρ η ( 0 ) θ ( t + τ ) θ ( t τ ) .
Here, the expression of the entropy η is obtained by integrating its equation of evolution (19), with regards to time variable.
Let us integrate the equality (33) on [ 0 , t ] × D , and take into account the border data (17) and the theorem of divergence. As such, we arrive at the proposed equality (31) and the demonstration of the proposition is over. □
If we combine the results from Equations (25) and (31), then we obtain a new useful equality
2 B ρ u ˙ i ( t ) u ˙ i ( t ) + I j k ϕ ˙ j m ( t ) ϕ ˙ k m ( t ) c θ 2 ( t ) d V = 2 D ρ u i ( 0 ) u ˙ i ( 0 ) + I j k ϕ j m ( 0 ) ϕ ˙ k m ( 0 ) d V + D ρ u ˙ i ( 0 ) u i ( 2 t ) + I j k ϕ ˙ j m ( 0 ) ϕ k m ( 2 t ) d V + D ρ u i ( 0 ) u ˙ i ( 2 t ) + I j k ϕ j m ( 0 ) ϕ ˙ k m ( 2 t ) d V 0 t D ρ η ( 0 ) 2 θ ( τ ) + θ ( t + τ ) θ ( t τ ) d V d τ ,
for all t [ 0 , ) .
Based on the previously demonstrated integral identities, we are able to address the main results of our study. First, we established a result of uniqueness for the solution of the backward in time problem. As a consequence, we approach the question of localization of the solutions of the backward in time problem.
Theorem 1.
At most, an ordered array u = u i , ϕ i j , θ can satisfy the equations and conditions of the backward problem P .
Proof. 
As usual, we will assume, by absurdum, that the problem would admit two solutions. The difference of the two solutions is also a solution, because the problem P is a linear one. Suffice it to show that this difference is null. For this we have to show that the problem P , for which the boundary and initial data are null, admits the null solution. It is clear that for the difference of two solutions, the boundary and initial conditions become homogeneous.
To simplify writing, we introduce the function M as a measure of the solution, defined by
M ( t ) = D ε 2 ρ u ˙ i ( t ) u ˙ i ( t ) + I j k ϕ ˙ j m ( t ) ϕ ˙ k m ( t ) + ( ε + 2 ) Ψ ( t ) + ε 2 c θ 2 ( t ) .
Here, ε is a small positive number.
Based on assumptions (a), (b), and (c), it can be deduced that the function M is positive.
Because the initial data are zero, the identity (18) received the simpler form:
B 1 2 ρ u ˙ i ( t ) u ˙ i ( t ) + I j k ϕ ˙ j m ( t ) ϕ ˙ k m ( t ) + Ψ ( t ) + 1 2 c θ 2 ( t ) d V = 0 t D 1 T 0 K i j θ , i ( τ ) θ , j ( τ ) d V d τ , t [ 0 , ) .
Analogously, the identity (21) becomes
B 1 2 ρ u ˙ i ( t ) u ˙ i ( t ) + I j k ϕ ˙ j m ( t ) ϕ ˙ k m ( t ) + Ψ ( t ) 1 2 c θ 2 ( t ) d V = 0 t D u ˙ i ( τ ) a i j + b j i θ ( τ ) , j + ϕ ˙ i j ( τ ) c i j k θ ( τ ) , k b i j ϕ ˙ i j ( τ ) θ ( τ ) + 1 T 0 K i j θ , i ( τ ) θ , j ( τ ) d V d τ , t [ 0 , ) .
If we use Equations (36) and (37), then the function M from (35) becomes
M ( t ) = 0 t D 2 u ˙ i ( τ ) a i j + b j i θ ( τ ) , j + ϕ ˙ i j ( τ ) c i j k θ ( τ ) , k b i j ϕ ˙ i j ( τ ) θ ( τ ) + 1 ε T 0 K i j θ , i ( τ ) θ , j ( τ ) d V d τ , t [ 0 , ) .
By direct derivation with regards to the variable t in (37), we are led to the following equality:
d M ( t ) d t = 2 D u ˙ i ( τ ) a i j + b j i θ ( τ ) , j + ϕ ˙ i j ( τ ) c i j k θ ( τ ) , k b i j ϕ ˙ i j ( τ ) θ ( τ ) + 1 ε 2 T 0 K i j θ , i ( τ ) θ , j ( τ ) d V , t [ 0 , ) .
With the help of Schwarz’ inequality and by using the arithmetic–geometric mean inequality, from (39) we can deduce the following inequality:
d M ( t ) d t C 1 D ρ u ˙ i ( t ) u ˙ i ( t ) + I j k ϕ ˙ j m ( t ) ϕ ˙ k m ( t ) + c θ 2 ( t ) d V + δ 1 + ε T 0 D K i j θ , i ( τ ) θ , j ( τ ) d V , t [ 0 , ) .
Now, we take into account that the internal energy density Ψ is a positive definite quadratic form, according to the hypothesis (c), using the definition (35) of the function M and choose δ 1 ε . Then, from (40) we obtain the inequality:
d M ( t ) d t C 1 ε D ε ρ u ˙ i ( t ) u ˙ i ( t ) + I j k ϕ ˙ j m ( t ) ϕ ˙ k m ( t ) + c θ 2 ( t ) d V C 1 ε M ,
and it is then clear that a solution to this inequality meets the next inequality:
0 M ( t ) M ( 0 ) e C 1 / ε .
We recall that for the difference of the two supposed solutions, the initial conditions are homogeneous, then we have M ( 0 ) = 0 , so from (42) we get
M ( t ) = 0 , t [ 0 , )
and this together with the assumptions leads to the conclusion that our problem has only the solution
u i ( t ) = 0 , ϕ i j ( t ) = 0 , θ ( t ) = 0 , t [ 0 , ) ,
and the proof of Theorem 1 is concluded. □
Our final result is dedicated to the partition of various energies associated with the solution of the backward in time problem P * . We recall that this problem consists of the equations of motion (9), the constitutive relations (4), the kinematic Equation (2), the initial data (12), and the boundary restrictions in their homogeneous form (11).
First, using the procedure outlined at the end of Section 2, we transform the boundary-final value problem P * into the boundary-initial value problem P . In this way, in what follows, we will be able to make the considerations only on the problem P . Let us denote by T the set of of those thermoelastodynamic processes defined in the cylinder ( , 0 ] × D , which satisfy the restriction
D 1 T 0 K i j θ , i ( τ ) θ , j ( τ ) d V C 1 ,
for all t [ 0 , ) . Here, C 1 is a given positive constant.
We will introduce the known Cesaro means necessary to evaluate the various types of energies that can be attached to a solution to the problem P . So, if u i , ϕ i j , θ is a solution of the mixed problem, then the Cesaro means are:
K ( t ) = 1 t 0 t D ρ u ˙ i ( τ ) u ˙ i ( τ ) + I j k ϕ ˙ j m ( τ ) ϕ ˙ j m ( τ ) d V d τ ; S ( t ) = 1 t 0 t D 1 2 A i j m n ε i j ( τ ) ε m n ( τ ) + D i j m n ε i j ( τ ) κ m n ( τ ) + F i j m n r ε i j ( τ ) χ m n r ( τ ) + 1 2 B i j m n κ i j ( τ ) κ m n ( τ ) + G i j m n r κ i j ( τ ) χ m n r ( τ ) + 1 2 C i j k m n r χ i j k ( τ ) χ m n r ( τ ) d V d τ ; R ( t ) = 1 t 0 t D 1 2 c θ 2 ( τ ) d V d τ ; D ( t ) = 1 t 0 t 0 τ D 1 T 0 K i j θ , i ( s ) θ , j ( s ) d V d s d τ .
In the particular case when meas( D u ) = 0 and meas( D ϕ ) = 0 , it can be determined a rigid displacement, a rigid dipolar displacement, and a null temperature, which satisfy the equations of motion (9), the constitutive relations (4), the kinematic relations (2), and verify the homogeneous boundary conditions (17). In this case, the initial data can be decomposed as follows:
u i 0 = u i + V i 0 , u ˙ i 0 = u ˙ i + V ˙ i 0 , ϕ i j 0 = ϕ i j + Ψ i j 0 , ϕ ˙ i j 0 = ϕ ˙ i j + Ψ ˙ i j 0 .
The rigid displacements u i , u ˙ i , ϕ i j , ϕ ˙ i j can be computed using the functions
W i ( ω ) = D ρ ω i D v , W j k ( ψ ) = D I m n ψ m j ψ n k d V ,
such that we have
W i ( V i 0 ) = 0 , W i ( V ˙ i 0 ) = 0 , W j k ( Ψ i j 0 ) = 0 , W j k ( Ψ ˙ i j 0 ) = 0 .
Together with known notation C 1 ( D ) , we will use the notation W n ( D ) for a Sobolev space defined on the domain D, and W n ( D ) = [ W n ( D ) ] 3 .
Other new notations:
C 1 ( D ) = u i , ϕ i j C 1 ( D ) 3 × C 1 ( D ) 9 : u i = 0 on D u , ϕ i j = 0 on D ϕ ; if meas ( D u ) = 0 and meas ( D ϕ ) = 0 , then W i ( u i ) = 0 , W j k ( ϕ i j 0 ) = 0 ; C ˜ 1 ( D ) = θ C 1 ( D ) : θ = 0 on D θ ; W 1 ( D ) the completion of C 1 ( D ) ; W ˜ 1 ( D ) the completion of C ˜ 1 ( D ) ,
the completion is by means of the original norm of the respective Sobolev space. Based on definition (13), and taking into account the hypothesis (15), we can obtain the following inequality (see [34]):
D Ψ ( u ) d V M 1 2 D ρ u i u i + I j k ϕ j m ϕ k m d V ,
for any u = u i , ϕ i j W 1 ( D ) . M 1 > 0 is a convenient chosen constant.
On the other hand, if we take into account the hypothesis (14), we can obtain the next inequality of Poincare’ type:
D K i j θ , i θ , j d V M 2 D θ 2 d V ,
for any θ W ˜ 1 ( D ) . M 2 > 0 is a convenient chosen constant.
Let us consider a solution u i , ϕ i j , θ to the problem P in the particular case when meas( D u ) = 0 and meas( D ϕ ) = 0 . We can represent this solution in the following form:
u i ( t , x ) = u i ( t , x ) + u i ( t , x ) + t u ˙ i ( t , x ) , θ ( t , x ) = ϑ ( t , x ) , ϕ i j ( t , x ) = ψ i j ( t , x ) + ϕ i j ( t , x ) + t ϕ ˙ i j ( t , x ) , ( t , x ) [ 0 , ) × D ,
in which u i , ψ i j , ϑ W 1 ( D ) × W ˜ 1 ( D ) . In addition, the ordered array u i , ψ i j , ϑ satisfies the problem P , which corresponds to the following initial data:
u i ( 0 , x ) = V i ( x ) , u ˙ i ( 0 , x ) = V ˙ i ( x ) , θ ( t , x ) = ϑ ( x ) , ϕ i j ( 0 , x ) = Ψ i j ( x ) , ϕ ˙ i j ( 0 , x ) = Ψ ˙ i j ( x ) , x D .
We now have everything prepared to address the problem of the equipartition of the various types of energies associated with the solution to the problem P .
Theorem 2.
Consider a solution u i , ϕ i j , θ of the backward in time problem P . If the initial data satisfy the following conditions:
u = u i W 1 ( D ) , u ˙ = u ˙ i W 0 ( D ) , ϕ = ϕ i j W 1 ( D ) , ϕ ˙ = ϕ ˙ i j W 0 ( D ) , θ W 0 ( D ) ,
then the following three statements are true:
i) The thermal component of energy, R, vanishes as t :
lim t R ( t ) = 0 ;
ii) if meas( D u ) = 0 and meas( D ϕ ) = 0 , then we have
lim t K ( t ) = lim t S ( t ) + 1 2 D ρ u ˙ i u ˙ i + I j k ϕ ˙ j m ϕ ˙ k m d V ,
lim t D ( t ) = 2 lim t K ( t ) 1 2 D ρ u ˙ i u ˙ i + I j k ϕ ˙ j m ϕ ˙ k m d V E ( 0 ) = 2 lim t S ( t ) + 1 2 D ρ u ˙ i u ˙ i + I j k ϕ ˙ j m ϕ ˙ k m d V E ( 0 ) ;
iii) if meas( D u ) 0 or meas( D ϕ ) 0 , then we have
lim t K ( t ) = lim t S ( t ) , lim t D ( t ) = 2 lim t K ( t ) E ( 0 ) = 2 lim t S ( t ) E ( 0 ) ,
where, to simplify writing in (52) and (53), we used the notation:
E ( t ) = D 1 2 ρ u ˙ i ( t ) u ˙ i ( t ) + I j k ϕ ˙ j m ( t ) ϕ ˙ k m ( t ) + Ψ ( t ) + 1 2 c θ 2 ( t ) d V ,
the internal energy density Ψ being defined in (13).
Proof. 
Taking into account the relations (18), (54), and (44) we immediately deduce the following equality:
K ( t ) + S ( t ) + R ( t ) = D ( t ) + E ( 0 ) , t ( 0 , ) .
Now, we consider the restrictions (43) and (48), the definitions (44), and the notation (54) in order to deduce the following estimation:
R ( t ) T 0 C 1 2 M 2 t max x D { c ( x ) } , t ( 0 , ) .
Clearly, since c ( x ) is a continuous function, we deduce that max x D { c ( x ) } is bounded, such that after we pass to the limit in (56) for t , we obtain the result (50).
If we consider the identity (34) and use the notations (44), we are led to the the following identity:
K ( t ) S ( t ) R ( t ) = 1 2 t D ρ u i ( 0 ) u ˙ i ( 0 ) + I j k ϕ j m ( 0 ) ϕ ˙ k m ( 0 ) d V + 1 4 t D ρ u ˙ i ( 0 ) u i ( 2 t ) + I j k ϕ ˙ j m ( 0 ) ϕ k m ( 2 t ) d V + 1 4 t D ρ u i ( 0 ) u ˙ i ( 2 t ) + I j k ϕ j m ( 0 ) ϕ ˙ k m ( 2 t ) d V 1 4 t 0 t D ρ η ( 0 ) 2 θ ( τ ) + θ ( t + τ ) θ ( t τ ) d V d τ ,
for all t [ 0 , ) .
On the other hand, taking into account the hypothesis (15), the inequalities (43) and (48), and the identity (18) we obtain the estimates
D ρ u ˙ i ( t ) u ˙ i ( t ) d V 2 ( C 1 + E ( 0 ) ) , D I j k ϕ ˙ j m ( t ) ϕ ˙ k m ( t ) d V 2 ( C 1 + E ( 0 ) ) , c 0 D ρ θ 2 ( t ) d V 2 ( C 1 + E ( 0 ) ) , D 2 Ψ ( t ) d V 2 ( C 1 + E ( 0 ) ) ,
where c 0 is from hypothesis (a), C 1 is from (43), and E from (54).
In (57) we use the Schwarz’s inequality and taking into account the estimates (58), we pass the limit—for t —so we get the equality
lim t K ( t ) = lim t S ( t ) + lim t 1 4 t D ρ u ˙ i ( 0 ) u i ( 2 t ) + I j k ϕ ˙ j m ( 0 ) ϕ ˙ k m ( 2 t ) d V .
In this way, it will be easy to demonstrate the relation (53), if we show that the integral from the right-hand side of the identity (59) is bounded.
For this aim, we will use the fact that meas( D u ) 0 or meas( D ϕ ) 0 and u i , ψ i j W 1 ( D ) . Furthermore, we consider the relations (47), (54), and (18) in order to get the following estimates:
M 1 D ρ u i ( t ) u i ( t ) d V D 2 Ψ ( t ) d V 2 ( C 1 + E ( 0 ) ) , M 1 D I j k ϕ j m ( t ) ϕ k m ( t ) d V D 2 Ψ ( t ) d V 2 ( C 1 + E ( 0 ) ) ,
such that it is easy, after we apply the Schwarz’s inequality, to deduce
lim t 1 4 t D ρ u ˙ i ( 0 ) u i ( 2 t ) + I j k ϕ ˙ j m ( 0 ) ϕ k m ( 2 t ) d V = 0 .
From (59) and (60) we obtain the first result from (53), and after we consider the equality (55), the second equality from (53) is proven.
Finally, we will prove the equalities (51) and (52). In this regard, we take into account that meas( D u ) = 0 and meas( D ϕ ) 0 such that with the help of the decompositions (45) and (49), and conditions (46), we are led to the equality
1 4 t D ρ u ˙ i ( 0 ) u i ( 2 t ) + I j k ϕ ˙ j m ( 0 ) ϕ k m ( 2 t ) d V = 1 4 t D ρ u ˙ i u i ( 2 t ) d V + 1 4 t D ρ u ˙ i + V ˙ i 0 u i ( 2 t ) d V + 1 2 D ρ u ˙ i u ˙ i d V + 1 4 t D I j k ϕ ˙ j m ϕ k m d V + 1 4 t D I j k ϕ ˙ j m + Ψ ˙ j m 0 ψ j k ( 2 t ) d V + 1 2 D I j k ϕ ˙ j m ϕ ˙ k m d V .
Let us observe that relations (47), (54), and (18) involve the estimates
M 1 D ρ u i ( t ) u i ( t ) d V 2 ( C 1 + E ( 0 ) ) , M 1 D I j k ψ j m ( t ) ψ k m ( t ) d V 2 ( C 1 + E ( 0 ) ) ,
such that from (61) we are led to the equality
lim t 1 4 t D ρ u ˙ i ( 0 ) u i ( 2 t ) + I j k ϕ ˙ j m ( 0 ) ϕ k m ( 2 t ) d V = 1 2 D ρ u ˙ i u ˙ i d V + 1 2 D I j k ϕ ˙ j m ϕ ˙ k m d V .
Now, we substitute (62) into Equation (59) and then obtain the relation (51). Lastly, we consider the relations (50), (51), and (55) in order to obtain (52). With this, the proof of the theorem is completed. □

4. Conclusions

We first formulate the mixed backward in time problem in the context of thermoelasticity for dipolar materials. To prove the consistency of this mixed problem, our first main result is regarding the uniqueness of the solution for this problem. This is obtained based on some auxiliary results, namely, four integral identities. The second main result is regarding the temporal behavior of our thermoelastic body with dipolar structure. This behavior is studied by means of some relations on partition of various parts of the energy associated to the solution of the problem. After we introduce the Cesaro means for all parts of the total energy, we can evaluate the asymptotic partition of these parts. We have to say that the kinetic energy and potential energy become asymptotically equal when the variable time tends to infinity.

Author Contributions

All four authors conceived the framework and structured the whole manuscript, checked the results and completed the revision of the paper. The authors have equally contributed to the elaboration of this manuscript.

Funding

This research received no external funding.

Acknowledgments

The authors would like to express their gratitude to the anonymous reviewers for their constructive comments and suggestions, which led to an improved form of the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Eringen, A.C. Theory of thermo-microstretch elastic solids. Int. J. Eng. Sci. 1990, 28, 1291–1301. [Google Scholar] [CrossRef]
  2. Eringen, A.C. Microcontinuum Field Theories; Springer: New York, NY, USA,, 1990. [Google Scholar]
  3. Hassan, M.; Marin, M.; Ellahi, R.; Alamri, S.Z. Exploration of convective heat transfer and flow characteristics synthesis by Cu–Ag/water hybrid-nanofluids. Heat Transf. Res. 2018, 49, 1837–1848. [Google Scholar] [CrossRef]
  4. Marin, M.; Ellahi, R.; Chirila, A. On solutions of Saint-Venant’s problem for elastic dipolar bodies with voids. Carpathian J. Math. 2017, 33, 219–232. [Google Scholar]
  5. Mindlin, R.D. Micro-structure in linear elasticity. Arch. Ration. Mech. Anal. 1964, 51–78. [Google Scholar] [CrossRef]
  6. Green, A.E.; Rivlin, R.S. Multipolar continuum mechanics. Arch. Ration. Mech. Anal. 1964, 17, 113–147. [Google Scholar] [CrossRef]
  7. Fried, E.; Gurtin, M.E. Thermomechanics of the interface between a body and its environment. Contin. Mech. Therm. 2007, 19, 253–271. [Google Scholar] [CrossRef]
  8. Nunziato, J.W.; Cowin, S.C. A nonlinear theory of materials with voids. Arch. Ration. Mech. Anal. 1979, 72, 175–201. [Google Scholar] [CrossRef]
  9. Cowin, S.C.; Nunziato, J.W. Linear elastic materials with voids. J. Elast. 1983, 13, 125–147. [Google Scholar] [CrossRef]
  10. Goodman, M.A.; Cowin, S.C. A continuum theory of granular material. Arch. Ration. Mech. Anal. 1971, 44, 249–266. [Google Scholar] [CrossRef]
  11. Othman, M.I.A.; Marin, M. Effect of thermal loading due to laser pulse on thermo-elastic porous medium under G-N theory. Results Phys. 2017, 7, 3863–3872. [Google Scholar] [CrossRef]
  12. Abbas, I.A. A GN model based upon two-temperature generalized thermoelastic theory in an unbounded medium with a spherical cavity. Appl. Math. Comput. 2014, 245, 108–115. [Google Scholar] [CrossRef]
  13. Abbas, I.A. Eigenvalue approach for an unbounded medium with a spherical cavity based upon two-temperature generalized thermoelastic theory. J. Mech. Sci. Technol. 2014, 28, 4193–4198. [Google Scholar] [CrossRef]
  14. Othman, M.I.A.; Hasona, W.M.; Abd-Elaziz, E.M. Effect of Rotation on Micropolar Generalized Thermoelasticity with Two-Temperatures using a Dual-Phase-Lag Model. Can. J. Phys. 2014, 92, 149–158. [Google Scholar] [CrossRef]
  15. Marin, M. Cesaro means in thermoelasticity of dipolar bodies. Acta Mech. 1997, 122, 155–168. [Google Scholar] [CrossRef]
  16. Marin, M.; Craciun, E.M. Uniqueness results for a boundary value problem in dipolar thermoelasticity to model composite materials. Compos. Part B Eng. 2017, 126, 27–37. [Google Scholar] [CrossRef]
  17. Serrin, J. The Initial Value Problem for the Navier-Stokes Equations. In Nonlinear Problems Proc. Sympos.; Univ. of Wisconsin Press: Madison, WI, USA, 1963; pp. 69–98. [Google Scholar]
  18. Ciarletta, M.; Scalia, A. Some Results in Linear Theory of Thermomicrostretch Elastic Solids. Meccanica 2004, 39, 191–206. [Google Scholar] [CrossRef]
  19. Knops, R.J.; Payne, L.E. On the Stability of Solutions of the Navier-Stokes Equations Backward in Time. Arch. Ration. Mech. Anal. 1968, 29, 331–335. [Google Scholar] [CrossRef]
  20. Galdi, G.P.; Straughan, B. Stability of Solutions of the Navier-Stokes Equations Backward in Time. Arch. Ration. Mech. Anal. 1988, 101, 107–114. [Google Scholar] [CrossRef]
  21. Payne, L.; Straughan, B. Improperly Posed and Nonstandard Problems for Parabolic Partial Differential Equations. In Elasticity: Mathematical Methods and Applications; Ellis-Horwood Pub: Chichester, UK, 1990; pp. 273–300. [Google Scholar]
  22. Ames, K.A.; Payne, L.E. Stabilizing Solutions of the Equations of Dynamical Linear Thermoelasticity Backward in Time. Stab. Appl. Anal. Contin. Media 1991, 1, 243–260. [Google Scholar]
  23. Ciarletta, M. On the Uniqueness and Continuous Dependence of Solutions in Dynamical Thermoelasticity Backward in Time. J. Therm. Stress. 2002, 25, 969–984. [Google Scholar] [CrossRef]
  24. Ciarletta, M.; Chirita, S. Spatial Behavior in Linear Thermoelasticity Backward in Time. In Proceedings of the Fourth International Congress on Thermal Stresses, Osaka, Japan, 8–11 June 2001; pp. 485–488. [Google Scholar]
  25. Quintanilla, R. Impossibility of Localization in Linear Thermoelasticity with Voids. Mech. Res. Commun. 2007, 34, 522–527. [Google Scholar] [CrossRef]
  26. Green, A.E.; Naghdi, P.M. On Undamped Heat Waves in an Elastic Solid. J. Therm. Stress. 1992, 15, 253–264. [Google Scholar] [CrossRef]
  27. Green, A.E.; Naghdi, P.M. Thermoelasticity without Energy Dissipation. J. Elast. 1993, 31, 189–208. [Google Scholar] [CrossRef]
  28. Iovane, G.; Passarella, F. Saint-Venant’s Principle in Dynamic Porous Thermoelastic Media with Memory for Heat Flux. J. Therm. Stress. 2004, 27, 983–999. [Google Scholar] [CrossRef]
  29. Passarella, F.; Tibullo, V. Some Results in Linear Theory of Thermoelasticity Backward in Time for Microstretch Materials. J. Therm. Stress. 2010, 33, 559–576. [Google Scholar] [CrossRef]
  30. Quintanilla, R.; Straughan, B. Energy Bounds for Some Non-standard Problems in Thermoelasticity. Proc. R. Soc. A Math. Phys. Eng. Sci. 2005, 461, 1–15. [Google Scholar] [CrossRef]
  31. De Sciarra, F.M.; Salerno, M. On thermodynamic functions in thermoelasticity without energy dissipation. Eur. J. Mech. A Solids 2014, 46, 84–95. [Google Scholar] [CrossRef]
  32. Canadija, M.; Barretta, R.; de Sciarra, F.M. A gradient elasticity model of Bernoulli-Euler nanobeams in non-isothermal environments. Eur. J. Mech. A Solids 2016, 55, 243–255. [Google Scholar] [CrossRef]
  33. Barretta, R.; Canadija, M.; Luciano, R.; de Sciarra, F.M. Stress-driven modeling of nonlocal thermoelastic behavior of nanobeams. Int. J. Eng. Sci. 2018, 126, 53–67. [Google Scholar] [CrossRef] [Green Version]
  34. Hlavacek, I.; Hlavacek, M. On the Existence and Uniqueness of Solution and Some Variational Principles in Linear Theories of Elasticity with Couple-Stresses. I. Cosserat Continuum. Appl. Math. 1969, 14, 387–410. [Google Scholar]

Share and Cite

MDPI and ACS Style

Marin, M.; Vlase, S.; Ellahi, R.; Bhatti, M.M. On the Partition of Energies for the Backward in Time Problem of Thermoelastic Materials with a Dipolar Structure. Symmetry 2019, 11, 863. https://doi.org/10.3390/sym11070863

AMA Style

Marin M, Vlase S, Ellahi R, Bhatti MM. On the Partition of Energies for the Backward in Time Problem of Thermoelastic Materials with a Dipolar Structure. Symmetry. 2019; 11(7):863. https://doi.org/10.3390/sym11070863

Chicago/Turabian Style

Marin, M., S. Vlase, R. Ellahi, and M.M. Bhatti. 2019. "On the Partition of Energies for the Backward in Time Problem of Thermoelastic Materials with a Dipolar Structure" Symmetry 11, no. 7: 863. https://doi.org/10.3390/sym11070863

APA Style

Marin, M., Vlase, S., Ellahi, R., & Bhatti, M. M. (2019). On the Partition of Energies for the Backward in Time Problem of Thermoelastic Materials with a Dipolar Structure. Symmetry, 11(7), 863. https://doi.org/10.3390/sym11070863

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop