1. Introduction
As is well known, the Second Law can be expressed in terms of a deterministic inequality (e.g., [
1])
where
is the irreversible part of the entropy increment
. The later quantity can at any given time be written as the sum
in which
is the reversible part of
S, while
θ is the absolute temperature.
While random fluctuations are negligible on macroscales, the Second Law gets spontaneously violated on very small (molecular) scales, as expressed by the so-called
fluctuation theorem (
FT) which gives the relative probability of observing processes that have positive (
A) and negative (
) total dissipation in non-equilibrium systems [
2,
3]:
Here
is the total dissipation for a trajectory
of
N particles originating at
and evolving for a time
t:
The integral in (3) involves an instantaneous dissipation function:
It follows that the Second Law is correct as either an ensemble, or a temporal, or a volume average. Clearly, the FT represents a stochastic generalization of the Second Law. Depending on the particular assumptions/conditions, there exist similar versions of the FT [
2]. Note that experiments validating the FT have been conducted [
4,
5]. The fluctuations are Gaussian.
There are three types of classical physics phenomena where spontaneous violations occur–viscous flow [
6], heat conduction [
7], and electrical resistance [
3]—and in this paper we review the recently introduced extensions of continuum mechanics incorporating the first two of these. See also Figures 1 and 2 in [
8] for molecular fluid dynamics simulations showing such violations.
2. Irreversible Entropy Evolves as a Submartingale
In view of the random fluctuations,
is a stochastic process with a specific type of memory effect to be examined as follows. First, every stochastic process is defined with reference to a probability space
, where Ω is the sample space (of elementary events),
is the
field, and
the probability measure, the argument
being employed to indicate an elementary event as well as the random character of
. We now switch from a continuous (
t) to a discrete (
n) time parametrization
the reason for this switch being that the analytical aspects of discrete-time stochastic processes are simpler than those of continuous-time processes; the integral in (4) is replaced by a summation, while the derivative in (5) is understood in a finite-difference sense.
Our growing knowledge of the process
at the successive times (i.e., its history) is represented by a so-called
filtration on Ω: a sequence
of sub-sigma fields of
such that for all time instants
,
. In view of (5), we observe that this inequality is satisfied
which indicates that
is a
submartingale. On the technical side dictated by the probability theory, (7) has to be accompanied by two more conditions: (i)
is a filtration and
is adapted to
; (ii) for each
n,
is integrable.
If the ≥ sign in (7) were replaced by an equality sign, we would have a so-called
martingale. In fact, this observation acquires more light in view of the so-called
Doob decomposition [
9] saying that any submartingale is the sum of a martingale (
M) and an increasing process (
G): Let
be a submartingale relative to the filtration
. Then there exists a martingale
and a process
such that
- (i)
M is a martingale relative to ;
- (ii)
G is an increasing process: almost everywhere;
- (iii)
is -measurable ;
- (iv)
.
In [
10] we have employed an analogous (Doob–Meyer decomposition) theorem in continuous time, also giving a unique decomposition of a submartingale into a martingale and a “drift” process. The discrete time case should be sufficient for most continuum physics applications, while allowing a simpler analytical treatment.
Since we are interested in continuum mechanics/physics, we write Equations (1) and (2) in terms of the rates of entropy densities (
,
,
)
where
(
) is the heat flux. As we are focusing on continuum physics, we use equivalent notations of subscript-type (
) and symbolic-type (
) for tensors. While an overdot indicates a material time derivative corresponding to a total differential of a given quantity, a superscript * is used to indicate a time rate of a quantity which is not a state function [
1,
11].
While the inequality in (8)3 is assumed to hold instantaneously (i.e., for ) in conventional continuum mechanics/physics, in fact it holds only for statistical averages (such as (7) above), or for sufficiently long time averages, or for sufficiently large systems. We are interested in extending the conventional continuum mechanics/physics so as to account for violations of the Second Law according to FT.
3. Stochastic Thermomechanics
3.1. Basic Framework
Since the fluctuation theorem involves a dissipation function, the natural framework to generalize the continuum mechanics is one employing the thermodynamics with internal variables (TIV) [
1], where the key role is played by the internal energy
u and the entropy production rate
. Thus,
is the functional of strain
(assuming for simplicity that we work with small elastic strains) and entropy
s, satisfying the first law of thermodynamics (the energy balance)
where
ρ is the mass density,
is the Cauchy stress,
is the deformation rate, and
is the heat flux. Note that the first term on the right is the specific power of deformation. By a well-known partial Legendre transform, the free energy is
Following the standard derivations of TIV [
11], on account of (8)
, (8)
, (9), and (10), we obtain the relation
With
ψ playing the role of a potential for the entropy (
s) and the quasi-conservative Cauchy stress (
), we have
The stress is the dissipative part of the total stress ; it appears in the discussion of irreversible phenomena immediately below.
If we were to assume (8)
, we would obtain the standard Clausius–Duhem inequality
where we employ a standard representation of
as a scalar product of two conjugate vectors
If we admit violations of the Second Law according to the FT discussed in the first section of this paper, the dissipation function
is neither strictly non-negative for all times, nor deterministic. Rather,
is a scale-dependent stochastic functional
, which satisfies the Clausius–Duhem inequality only in the ensemble average sense
The explicit dependence of
ϕ on
ω (
) indicates the stochastic character of
ϕ. The basic model, sufficiently general for our purposes and valid for linear viscous (Stokes-type) and thermal (Fourier-type) fluids, involves an additive decomposition of
ϕ into two parts: one mechanical (
) and another thermal (
):
Now, the randomness itself is introduced by taking each of these parts as a sum of a deterministic mean (which must be positive-valued) plus a fluctuating part which may spontaneously take negative values, see (19) below.
3.2. Axioms of Continuum Thermomechanics
In view of the preceding discussion, the fundamental rôle in TIV [
1] is played by the free energy (
ψ) and the dissipation function (
ϕ), the latter quantity, upon (temporal, volume, or statistical) averaging, giving the positive entropy production rate. Clearly, TIV lends itself to an extension replacing the Second Law by the FT. Notably, an amazingly wide range of continuum constitutive behaviors can be obtained from the pair
[
1,
11,
12].
By contrast, in the rational continuum mechanics (RCM) one works with a quartet of the stress state
, heat flux
, free energy density
ψ, and entropy
s, all of which are taken as functionals of the history at a continuum point. As outlined in [
13], the axioms (also called
principles) of RCM need to be modified so as to admit a negative entropy production:
1. The Axiom of Causality: “
The future state of the system depends solely on the probabilities of events in the past”. That is, quoting Evans & Searles [
3]: “
the probability of subsequent events can be predicted from the probabilities of finding initial phases and a knowledge of preceding changes in the applied field and environment of the system.”
2. The FT is derived from the Axiom of Causality. While the entropy
is a submartingale in time (
), it is a random field in spatial domain
:
3. The Second Law in its conventional deterministic form is obtained–upon spatial, statistical, or time averaging–as a special case of the fluctuation theorem, and the entire field of continuum mechanics is recovered as a special case.
4. The Axiom of Determinism: “at any instant t, the value of at depends on the whole history of the thermokinetic process up to the time t”. The choice of the thermokinetic process depends on the particular physics involved, and may take the form of a classical or non-classical (e.g., micropolar) theory or a stochastic field.
5. The Axiom of Local Action is to be replaced by the scale dependence of the adopted continuum approximation for, clearly, that axiom makes no clear reference to the microstructure of the medium.
6. The Axiom of Equipresence (all the constitutive quantities depend a priori on the same variables) is to be abandoned since the violation of Second Law may occur in one physical process present in constitutive relations, not all.
3.3. Random Fields
In view of the preceding discussion, the deterministic continuum theory must be replaced by a stochastic one: working with TIV, the free energy density
ψ and the dissipation function
ϕ are real-valued random fields over the material (
) and time (
) domains:
For example, considering the thermal dissipation, we haveo
Given the spontaneous, randomly occurring violations of the Second Law, the dissipation function reads
In the case of linear Fourier-type conductivity, (20) becomes more explicit with
Here
involves the deterministic thermal conduction
(which is positive definite) and
, with
M being the martingale modeling the random fluctuations according to the Doob decomposition (iv) in
Section 2. Clearly, the randomness residing in
allows the total thermal conductivity
to sometimes become negative-definite, since
is not required to be positive-definite, thus signifying the violations of the Second Law. More specifically,
(where
is a real vector space) is a second-order, rank 2 tensor random field (e.g., [
14,
15])
In view of the Gaussian character of nanoscale fluctuations (recall
Section 1),
is a Gaussian tensor random field, making negative values possible.
By analogy, the same type of approach may be used to introduce fluctuations in mechanical dissipation having spontaneous violations of the Second Law, recall (15).
4. Thermoviscous Fluids
We start from the free energy and the dissipation function
and assume the macroscopic incompressibility
The free energy function
ψ (taken per unit volume), assuming no elastic response but the presence of a relaxation time
, is
Here and are the free energy and entropy in the reference state, μ is the shear elastic modulus, is the specific heat at constant strain, and ϑ is the temperature difference from the reference temperature . The last term on the right hand side is taken by analogy to the thermoelasticity with two relaxation times in the next section, so as to retain the Fourier-type heat conduction, but to obey the hyperbolic (finite speed) heat propagation. In this section, an overdot denotes a material derivative for absolute tensors (like the temperature gradient and deformation rate) and an Oldroyd derivative for tensor densities (like the heat flux and stress tensor).
The free energy being a potential for quasi-conservative stresses
and the entropy
s, we find
The first relation above reflects the presence of a hydrostatic pressure and the fact that ψ does not depend on the strain since the fluid under consideration does not have any elasticity. The relation (26) is immediately identified as the constitutive equation for entropy.
Henceforth, in view of the fluid’s incompressibility in (24)
and to simplify the notation,
is taken as the deviatoric part of the quasi-conservative stress tensor. As always in TIV [
1], there holds the relation
where
is the dissipative stress, which now also represents the deviatoric part only.
This Ansatz leads to the Clausius-Duhem inequality in the form
where
is the dissipative stress, which is now equal to the total stress
. As discussed earlier, the inequality (28) may spontaneously be violated according to the FT.
On account of the equality in (28), we take the specific (per unit volume) dissipation
ϕ as a functional of temperature gradient
and deformation rate
:
Therefore, the inequality in (28) may be stated in terms of the scalar product:
.
Next, for the entropy production rate we adopt the dissipation functional in the space of velocity
:
where
k is the Fourier conductivity and
H the fluid viscosity, both parameters modeling an isotropic medium. The Equation (
30) is a special case of (29), with both processes being effectively compound [
1]. Using the thermodynamic orthogonality, (30) yields
Collecting the three parts of the constitutive law–mechanical (Stokesian), thermal (Fourier law), and entropy- and introducing the FT-type fluctuations, we have
Here, on account of
Section 3, we have introduced the martingale-type fluctuation
in thermal response (32)
and an analogous one (
) in the viscous response (32)
. Relation (32)
shows that, while the Fourier-type law holds, there is a relaxation effect involved in the entropy. If the scale of the elementary volume
of continuum mechanics approximation is very small, that element will see violations of the Second Law according to [
7]. In this model the heat is conducted with finite speeds—i.e., not infinite speeds as would be the case with
. In other words, instead of having a parabolic (diffusion) equation for temperature, we have (by application of the energy balance (24) and the entropy-temperature relation (32)
)
Here we have also used the approximation of small temperature fluctuations. In effect,
ϑ (and also
θ) is governed by the telegraph (damped hyperbolic) Equation (
33), whose limiting case (for
) is the conventional (parabolic) heat conduction equation.
6. Closure
The situations where violations of the Second Law are relevant involve very small length scales (extremely thin wavefronts, nano-channels, nano-rods, ...) and short time scales (although times up to 3 s have been observed for cholesteric liquids [
3]). Thus, a natural setting in which to expect such violations is that of poromechanics involving nanoscale channels, where the Second Law may apply only on average in the fluid within the channels (viscous flow and heat conduction) and the solid skeleton (heat conduction). The third possible phenomenon (electrical resistance) is one more possibility, that still needs to be worked out. As a reference, we start from the classical poromechanics obeying the Second Law [
19], the Clausius-Duhem inequality is written in terms of the irreversible entropy production
(
) taking the form
where three possible contributions to dissipation are identified: (1) fluid dissipation:
; (2) thermal dissipation:
; (3) skeleton dissipation:
. Conventionally, each of these contributions to dissipation is assumed to satisfy its own second law inequality. It now follows that, in the case of poromechanics describing phenomena on very small space and time scales, the spontaneous violations of the second law (43) can occur in either one, two or three processes.