Clausius–Duhem inequality

Laws

Conservations
Energy Mass Momentum
Inequalities
Clausius–Duhem (entropy)

Solid mechanics

Contact mechanics (frictional)

Fluid mechanics

Fluids
Statics · Dynamics Archimedes' principle · Bernoulli's principle Navier–Stokes equations Poiseuille equation · Pascal's law Viscosity (Newtonian · non-Newtonian) Buoyancy · Mixing · Pressure
Liquids
Surface tension Capillary action
Gases
Atmosphere Boyle's law Charles's law Gay-Lussac's law Combined gas law
Plasma

Rheology

Smart fluids
Viscoelasticity Rheometry Rheometer
Magnetorheological Electrorheological Ferrofluids

Scientists

The Clausius–Duhem inequality^[1]^[2] is a way of expressing the second law of thermodynamics that is used in continuum mechanics. This inequality is particularly useful in determining whether the constitutive relation of a material is thermodynamically allowable.^[3]

This inequality is a statement concerning the irreversibility of natural processes, especially when energy dissipation is involved. It was named after the German physicist Rudolf Clausius and French physicist Pierre Duhem.

Clausius–Duhem inequality in terms of the specific entropy

The Clausius–Duhem inequality can be expressed in integral form as

{\cfrac {d}{dt}}\left(\int _{\Omega }\rho ~\eta ~{\text{dV}}\right)\geq \int _{\partial \Omega }\rho ~\eta ~(u_{n}-\mathbf {v} \cdot \mathbf {n} )~{\text{dA}}-\int _{\partial \Omega }{\cfrac {\mathbf {q} \cdot \mathbf {n} }{T}}~{\text{dA}}+\int _{\Omega }{\cfrac {\rho ~s}{T}}~{\text{dV}}.

In this equation $t\,$ is the time, $\Omega \,$ represents a body and the integration is over the volume of the body, $\partial \Omega \,$ represents the surface of the body, $\rho \,$ is the mass density of the body, $\eta \,$ is the specific entropy (entropy per unit mass), $u_{n}\,$ is the normal velocity of $\partial \Omega \,$ , $\mathbf {v}$ is the velocity of particles inside $\Omega \,$ , $\mathbf {n}$ is the unit normal to the surface, $\mathbf {q}$ is the heat flux vector, $s\,$ is an energy source per unit mass, and $T\,$ is the absolute temperature. All the variables are functions of a material point at $\mathbf {x}$ at time $t\,$ .

In differential form the Clausius–Duhem inequality can be written as

\rho ~{\dot {\eta }}\geq -{\boldsymbol {\nabla }}\cdot \left({\cfrac {\mathbf {q} }{T}}\right)+{\cfrac {\rho ~s}{T}}

where ${\dot {\eta }}$ is the time derivative of $\eta \,$ and ${\boldsymbol {\nabla }}\cdot (\mathbf {a} )$ is the divergence of the vector $\mathbf {a}$ .

Proof

Assume that

\Omega

is an arbitrary fixed control volume. Then

$u_{n}=0$ and the derivative can be taken inside the integral to give

\int _{\Omega }{\frac {\partial }{\partial t}}(\rho ~\eta )~{\text{dV}}\geq -\int _{\partial \Omega }\rho ~\eta ~(\mathbf {v} \cdot \mathbf {n} )~{\text{dA}}-\int _{\partial \Omega }{\cfrac {\mathbf {q} \cdot \mathbf {n} }{T}}~{\text{dA}}+\int _{\Omega }{\cfrac {\rho ~s}{T}}~{\text{dV}}.

Using the divergence theorem, we get

\int _{\Omega }{\frac {\partial }{\partial t}}(\rho ~\eta )~{\text{dV}}\geq -\int _{\Omega }{\boldsymbol {\nabla }}\cdot (\rho ~\eta ~\mathbf {v} )~{\text{dV}}-\int _{\Omega }{\boldsymbol {\nabla }}\cdot \left({\cfrac {\mathbf {q} }{T}}\right)~{\text{dV}}+\int _{\Omega }{\cfrac {\rho ~s}{T}}~{\text{dV}}.

Since $\Omega$ is arbitrary, we must have

{\frac {\partial }{\partial t}}(\rho ~\eta )\geq -{\boldsymbol {\nabla }}\cdot (\rho ~\eta ~\mathbf {v} )-{\boldsymbol {\nabla }}\cdot \left({\cfrac {\mathbf {q} }{T}}\right)+{\cfrac {\rho ~s}{T}}.

Expanding out

{\frac {\partial \rho }{\partial t}}~\eta +\rho ~{\frac {\partial \eta }{\partial t}}\geq -{\boldsymbol {\nabla }}(\rho _{\eta })\cdot \mathbf {v} -\rho ~\eta ~({\boldsymbol {\nabla }}\cdot \mathbf {v} )-{\boldsymbol {\nabla }}\cdot \left({\cfrac {\mathbf {q} }{T}}\right)+{\cfrac {\rho ~s}{T}}

or,

{\frac {\partial \rho }{\partial t}}~\eta +\rho ~{\frac {\partial \eta }{\partial t}}\geq -\eta ~{\boldsymbol {\nabla }}\rho \cdot \mathbf {v} -\rho ~{\boldsymbol {\nabla }}\eta \cdot \mathbf {v} -\rho ~\eta ~({\boldsymbol {\nabla }}\cdot \mathbf {v} )-{\boldsymbol {\nabla }}\cdot \left({\cfrac {\mathbf {q} }{T}}\right)+{\cfrac {\rho ~s}{T}}

or,

\left({\frac {\partial \rho }{\partial t}}+{\boldsymbol {\nabla }}\rho \cdot \mathbf {v} +\rho ~{\boldsymbol {\nabla }}\cdot \mathbf {v} \right)~\eta +\rho ~\left({\frac {\partial \eta }{\partial t}}+{\boldsymbol {\nabla }}\eta \cdot \mathbf {v} \right)\geq -{\boldsymbol {\nabla }}\cdot \left({\cfrac {\mathbf {q} }{T}}\right)+{\cfrac {\rho ~s}{T}}.

Now, the material time derivatives of $\rho$ and $\eta$ are given by

{\dot {\rho }}={\frac {\partial \rho }{\partial t}}+{\boldsymbol {\nabla }}\rho \cdot \mathbf {v} ~;~~{\dot {\eta }}={\frac {\partial \eta }{\partial t}}+{\boldsymbol {\nabla }}\eta \cdot \mathbf {v} .

Therefore,

\left({\dot {\rho }}+\rho ~{\boldsymbol {\nabla }}\cdot \mathbf {v} \right)~\eta +\rho ~{\dot {\eta }}\geq -{\boldsymbol {\nabla }}\cdot \left({\cfrac {\mathbf {q} }{T}}\right)+{\cfrac {\rho ~s}{T}}.

From the conservation of mass ${\dot {\rho }}+\rho ~{\boldsymbol {\nabla }}\cdot \mathbf {v} =0$ . Hence,

{\rho ~{\dot {\eta }}\geq -{\boldsymbol {\nabla }}\cdot \left({\cfrac {\mathbf {q} }{T}}\right)+{\cfrac {\rho ~s}{T}}.}

Clausius–Duhem inequality in terms of specific internal energy

The inequality can be expressed in terms of the internal energy as

\rho ~({\dot {e}}-T~{\dot {\eta }})-{\boldsymbol {\sigma }}:{\boldsymbol {\nabla }}\mathbf {v} \leq -{\cfrac {\mathbf {q} \cdot {\boldsymbol {\nabla }}T}{T}}

where ${\dot {e}}$ is the time derivative of the specific internal energy $e\,$ (the internal energy per unit mass), ${\boldsymbol {\sigma }}$ is the Cauchy stress, and ${\boldsymbol {\nabla }}\mathbf {v}$ is the gradient of the velocity. This inequality incorporates the balance of energy and the balance of linear and angular momentum into the expression for the Clausius–Duhem inequality.

Proof

Using the identity

${\boldsymbol {\nabla }}\cdot (\varphi ~\mathbf {v} )=\varphi ~{\boldsymbol {\nabla }}\cdot \mathbf {v} +\mathbf {v} \cdot {\boldsymbol {\nabla }}\varphi$ in the Clausius–Duhem inequality, we get

\rho ~{\dot {\eta }}\geq -{\boldsymbol {\nabla }}\cdot \left({\cfrac {\mathbf {q} }{T}}\right)+{\cfrac {\rho ~s}{T}}\qquad {\text{or}}\qquad \rho ~{\dot {\eta }}\geq -{\cfrac {1}{T}}~{\boldsymbol {\nabla }}\cdot \mathbf {q} -\mathbf {q} \cdot {\boldsymbol {\nabla }}\left({\cfrac {1}{T}}\right)+{\cfrac {\rho ~s}{T}}.

Now, using index notation with respect to a Cartesian coordinate system $\mathbf {e} _{j}$ ,

{\boldsymbol {\nabla }}\left({\cfrac {1}{T}}\right)={\frac {\partial }{\partial x_{j}}}\left(T^{-1}\right)~\mathbf {e} _{j}=-\left(T^{-2}\right)~{\frac {\partial T}{\partial x_{j}}}~\mathbf {e} _{j}=-{\cfrac {1}{T^{2}}}~{\boldsymbol {\nabla }}T.

Hence,

\rho ~{\dot {\eta }}\geq -{\cfrac {1}{T}}~{\boldsymbol {\nabla }}\cdot \mathbf {q} +{\cfrac {1}{T^{2}}}~\mathbf {q} \cdot {\boldsymbol {\nabla }}T+{\cfrac {\rho ~s}{T}}\qquad {\text{or}}\qquad \rho ~{\dot {\eta }}\geq -{\cfrac {1}{T}}\left({\boldsymbol {\nabla }}\cdot \mathbf {q} -\rho ~s\right)+{\cfrac {1}{T^{2}}}~\mathbf {q} \cdot {\boldsymbol {\nabla }}T.

From the balance of energy

\rho ~{\dot {e}}-{\boldsymbol {\sigma }}:{\boldsymbol {\nabla }}\mathbf {v} +{\boldsymbol {\nabla }}\cdot \mathbf {q} -\rho ~s=0\qquad \implies \qquad \rho ~{\dot {e}}-{\boldsymbol {\sigma }}:{\boldsymbol {\nabla }}\mathbf {v} =-({\boldsymbol {\nabla }}\cdot \mathbf {q} -\rho ~s).

Therefore,

\rho ~{\dot {\eta }}\geq {\cfrac {1}{T}}\left(\rho ~{\dot {e}}-{\boldsymbol {\sigma }}:{\boldsymbol {\nabla }}\mathbf {v} \right)+{\cfrac {1}{T^{2}}}~\mathbf {q} \cdot {\boldsymbol {\nabla }}T\qquad \implies \qquad \rho ~{\dot {\eta }}~T\geq \rho ~{\dot {e}}-{\boldsymbol {\sigma }}:{\boldsymbol {\nabla }}\mathbf {v} +{\cfrac {\mathbf {q} \cdot {\boldsymbol {\nabla }}T}{T}}.

Rearranging,

{\rho ~({\dot {e}}-T~{\dot {\eta }})-{\boldsymbol {\sigma }}:{\boldsymbol {\nabla }}\mathbf {v} \leq -{\cfrac {\mathbf {q} \cdot {\boldsymbol {\nabla }}T}{T}}\qquad \qquad \square }

Dissipation

The quantity

{\mathcal {D}}:=\rho ~(T~{\dot {\eta }}-{\dot {e}})+{\boldsymbol {\sigma }}:{\boldsymbol {\nabla }}\mathbf {v} -{\cfrac {\mathbf {q} \cdot {\boldsymbol {\nabla }}T}{T}}\geq 0

is called the dissipation which is defined as the rate of internal entropy production per unit volume times the absolute temperature. Hence the Clausius–Duhem inequality is also called the dissipation inequality. In a real material, the dissipation is always greater than zero.

References

↑ Truesdell, Clifford (1952), "The Mechanical foundations of elasticity and fluid dynamics", Journal of Rational Mechanics and Analysis, 1: 125–300 .
↑ Truesdell, Clifford & Toupin, Richard (1960), "The Classical Field Theories of Mechanics", Handbuch der Physik, III, Berlin: Springer .
↑ Frémond, M. (2006), "The Clausius–Duhem Inequality, an Interesting and Productive Inequality", Nonsmooth Mechanics and Analysis, Advances in mechanics and mathematics, 12, New York: Springer, pp. 107–118, doi:10.1007/0-387-29195-4_10, ISBN 0-387-29196-2 .

External links

Memories of Clifford Truesdell by Bernard D. Coleman, Journal of Elasticity, 2003.
Thoughts on Thermomechanics by Walter Noll, 2008.

This article is issued from Wikipedia - version of the 6/14/2016. The text is available under the Creative Commons Attribution/Share Alike but additional terms may apply for the media files.