Next Article in Journal / Special Issue
Almost Periodic Solution for Forced Perturbed Non-Instantaneous Impulsive Model
Previous Article in Journal
Coefficient Inequalities for Multivalent Janowski Type q-Starlike Functions Involving Certain Conic Domains
Previous Article in Special Issue
On (k,ψ)-Hilfer Fractional Differential Equations and Inclusions with Mixed (k,ψ)-Derivative and Integral Boundary Conditions
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Molecular-Based Q-Tensor Hydrodynamic Theory of Smectic Liquid Crystals

School of Mathematics and Statistics, Guizhou University, Guiyang 550025, China
*
Author to whom correspondence should be addressed.
Axioms 2022, 11(10), 495; https://doi.org/10.3390/axioms11100495
Submission received: 26 August 2022 / Revised: 17 September 2022 / Accepted: 19 September 2022 / Published: 23 September 2022
(This article belongs to the Special Issue Impulsive, Delay and Fractional Order Systems)

Abstract

:
The Doi–Onsager molecular theory is capable of providing a rather accurate description of the local behavior of molecules; however, its computation is extremely time-consuming, since some higher-dimensional variables are typically involved. Therefore, establishing a computable reduced model that can capture essential physical properties is an important issue. In this work, we derived a reduced Q-tensor hydrodynamic theory that described smectic phases with density variations from the Doi–Onsager molecular theory using the Bingham closure approximation. The coefficients in the tensor model were derived from those in the molecular model. The energy dissipation law was inherited from the tensor model. Some special cases for the model were also discussed.

1. Introduction

Liquid crystals, which are composed of anisotropic non-spherical molecules, are featured in the local orientational order. There are several types of liquid crystal phases that can be classified by their orientational and positional orders, such as the nematic phase, the cholesteric phase and the smectic phase. Smectic liquid crystals are characterized by orientational and partial translation orders of the anisotropic molecules. The simplest of the smectics is the smectic-A phase, which exhibits a one-dimensional lamellar order in the direction of the preferred molecular orientational axis; it can be regarded as layers of two-dimensional fluids that are stacked upon each other. Details on smectic liquid crystals can be found in [1].
There are three main theories that are used to model smectic liquid crystals: the molecular theory, the tensor theory and the vector theory. The first is a microscopic kinetic theory derived from viewpoints of statistical mechanics, and the latter two are macroscopic theories based on continuum mechanics. The primary purpose of this work was to derive a new Q-tensor dynamical theory for smectic liquid crystals from the molecular theory of the Bingham closure. Before introducing the above three theories, some notations and conventions are provided.
Notations and conventions. Throughout this paper, we adopt the Einstein summation convention on repeated indices. We denote by Q the space of symmetric traceless tensors,
Q = def { Q R 3 × 3 : Q i j = Q j i , Q i i = 0 } ,
which is endowed with the inner product Q 1 : Q 2 = Q 1 i j Q 2 i j . The set Q is a five-dimensional subspace of R 3 × 3 . The matrix norm on Q is denoted by | Q | = def t r ( Q 2 ) = Q i j Q i j . For any two tensors Q 1 , Q 2 Q , we define ( Q 1 · Q 2 ) i j = Q 1 i k Q 2 k j . We use I to denote the 3 × 3 identity tensor.
We use a small bold-faced letter m to represent a vector. We denote by m 1 m 2 m k the tensor product of vectors m 1 , m 2 , , m k , and we usually omit the symbol ⊗ for simplicity. Again, the notation · σ is defined by ( · σ ) i = j σ i j , and S 2 denotes the unit sphere of R 3 . We also employ f , i to denote i f , such as Q i j , j = j Q i j .

1.1. Molecular Theory

Based on the Virial expansion, Onsager [2] first put forward the molecular theory to study the isotropic-nematic phase transition. We use x Ω R 3 to denote the material point and f ( x , m ) to represent the number density for the number of molecules whose orientation is parallel to m S 2 at point x . Onsager’s energy functional is given by (also see [3,4])
A ε [ f ] = Ω S 2 f ( x , m ) ( ln f ( x , m ) 1 ) + 1 2 U ε ( x , m ) f ( x , m ) d m d x , U ε ( x , m ) = Ω S 2 G ε ( x , m ; x , m ) f ( x , m ) d x d m .
In (1), the first term is the entropy, and the second term represents the interaction energy between each pair of two molecules. G ε ( x , m ; x , m ) is the interaction kernel between two molecules in the configurations ( x , m ) and ( x , m ) , and it takes one of the following forms:
(i) The hard-core interaction potential:
G ε ( x , m ; x , m ) = 1 , molecule ( x , m ) is disjoint with molecule ( x , m ) , 0 , joint with each other .
(ii) The long-range Maier–Saupe interaction potential:
G ε ( x , m ; x , m ) = 1 ε 3 / 2 g x x ε α | m × m | 2 ,
where the small parameter ε represents the typical interaction distance, and g ( x ) C ( R 3 ) is the radial function with R 3 g ( x ) d x = 1 . Again, the parameter α measures the intensity of the potential. It should be noted that (3) can be regarded as a smooth approximation of (2).
The kinetic theory of liquid crystal polymers, starting from the pioneering work of Doi [5], is based on Onsager’s molecular theory. This is the so-called Doi–Onsager theory, which had been greatly successful in predicting some dynamical properties of liquid crystal polymers in a solvent, such as tumbling, wagging, flow-aligning, log-rolling [6,7,8], etc; however, it is only valid for spatially homogeneous flows. The molecular theory for nonhomogeneous flows [3,9,10,11] was subsequently developed by introducing a nonlocal intermolecular potential.
The chemical potential μ ε is defined as
μ ε = δ A ε [ f ] δ f = ln f ( x , m , t ) + U ε ( x , m , t ) .
The non-dimensional Doi–Onsager equation takes the following form [3]:
f t + v · f = ε D e · γ m m + γ ( I m m ) · ( f μ ε ) + 1 D e R · ( f R μ ε ) R · ( m × κ · m f ) ,
v t + v · v = p + γ R e Δ v + 1 γ 2 R e · ( D : m m m m f ) + 1 γ D e R e ( · τ e + F e ) ,
· v = 0 ,
where D e and R e are called the Deborah number and the Reynolds number, respectively, R = m × m is the rotational gradient operator and γ ( 0 , 1 ) is a constant. Here p is the pressure to maintain the incompressibility (6), κ = ( v ) T is the gradient tenor of the velocity field v and D = 1 2 ( κ + κ T ) , while τ e and F e represent the elastic stress and the body force, which are respectively defined as
τ e = m m × R μ ε f , F e = μ ε f .
In (4), the constants γ = L 0 D e V 0 L 2 D and γ = L 0 D e V 0 L 2 D denote the translational diffusion coefficients parallel to and normal to the orientation of the liquid crystal molecule, respectively.
The system (4)–(6) has the following energy dissipation law ([3,12]):
d d t Ω 1 2 | v | 2 d x + 1 γ D e R e A ε [ f ] = Ω ( γ R e D : D + 1 γ 2 R e ( m m : D ) 2 f + 1 γ D e 2 R e R μ ε · R μ ε f + ε ( 1 γ ) D e 2 R e μ ε · ( γ m m + γ ( I m m ) ) · μ ε f ) d x .

1.2. Vector Theory

In this theory, the configuration at the point x Ω is described by the director ( x ) S 2 , a unit vector. The corresponding Oseen–Frank free energy takes the form [13,14],
F OF = Ω ( K 1 2 ( · n ) 2 + K 2 2 ( n · ( × n ) ) 2 + K 3 2 | n × ( × n ) | 2 + K 2 + K 4 2 tr ( n ) 2 ( · n ) 2 ) d x ,
where K 1 , , K 4 are the Frank elastic constants. The free energy (8) is only suitable for characterizing uniaxial nematic phases.
In order to characterize the smectic phase, besides the director n ( x ) S 2 , de Gennes [15] introduced a complex-valued order parameter,
Ψ ( x ) = ψ ( x ) e i q φ ( x ) .
The molecular mass density is defined by
δ ( x ) = ρ 0 ( x ) + 1 2 Ψ ( x ) + Ψ * ( x ) = ρ 0 ( x ) + ρ ( x ) cos q φ ( x ) ,
where ρ 0 is a locally uniform mass density, ρ ( x ) is the mass density of the smectic layer and φ parameterizes the layers so that φ is the direction of a layer’s normal. In addition, q is the wave number and 2 π / q is the layer thickness. If Ψ 0 then the smectics can become the nematics.
In (8), if we increase the terms with density variations, then it becomes the famous Chen–Lubensky model [16]:
F CL = F Ψ + F OF ,
where F Ψ contains terms up to the second order derivatives of the density Ψ ( x ) ,
F Ψ = 1 2 Ω ( a Ψ 2 + D 1 [ ( n · ) 2 Ψ ] 2 C 1 ( n · Ψ ) 2 + C 1 2 4 D 1 Ψ 2 + C 2 ( δ i j n i n j ) i Ψ j Ψ + D 2 [ ( δ i j n i n j ) i j Ψ ] 2 ) d x ,
where a = a 0 [ ( T T N A ) / T N A ] , and T N A is the transition temperature from the nematic state to the smectic-A state. In (10), the coefficients C 1 and D 1 determine the horizontal period, C 2 and D 2 determine the perpendicular period.
The hydrodynamic theory of nematic phases was established by Ericksen [17] and Leslie [18] in the 1960s. Based on the Ericksen–Leslie theory, the works [1,19,20] established the linear dynamical continuum theory of smectic-A liquid crystals. E [21] first set up a general nonlinear continuum theory for smectic-A liquid crystals that was applicable to situations with large deformations and non-trivial flows. Motivated by the approach of E, Stewart [22] also presented a nonlinear continuum theory for the dynamics of smectic-A liquid crystals, which included the usual Oseen constraint and the permeation.

1.3. Tensor Theory

The state of the nematic liquid crystals is described by a macroscopic tensor value order parameter Q ( x ) , which is a symmetric and traceless 3 × 3 matrix. Physically, it can be interpreted as the second-order traceless moment of the orientational distribution function f, that is,
Q ( x , m ) = S 2 ( m m 1 3 I ) f ( x , m ) d m .
Under this interpretation, the so-called physical constraint is that the eigenvalues of Q should satisfy
λ i ( Q ) ( 1 3 , 2 3 ) , for 1 i 3 .
The classic Landau–de Gennes energy functional takes the following general form [1],
F L G ( Q , Q ) = Ω { a 2 t r ( Q 2 ) b 3 t r ( Q 3 ) + c 4 ( t r ( Q 2 ) 2 bulk energy + 1 2 L 1 | Q | 2 + L 2 Q i j , j Q i k , k + L 3 Q i j , k Q i k , j + L 4 Q i j Q k l , i Q k l , j elastic energy } d x ,
where a , b , c are material-dependent and temperature-dependent non-negative constants and L i ( i = 1 , 2 , 3 , 4 ) are material-dependent elastic constants. We refer to [1,23] for more details.
There are two types of Q-tensor hydrodynamics for nematic liquid crystal. The first type is based on variational methods, such as the Beris–Edwards model [24] and the Qian–Sheng model [25]. The second type is derived from the molecular kinetic theory of closure approximations, which is called the molecular-based Q-tensor dynamical theory [4]. In such models, the evolution of Q is derived from the evolution of the probability density function f ( x , m ) by the relation (11); however, the equations of lower-order moments of f ( x , m ) obtained from kinetic theory are coupled with higher-order moments, which causes a problem. So, we need to express the higher-order moment as a function of the lower-order moment using the closure approximation.
Since the above three theories are derived from different physical considerations, exploring the connection between the different theories is not only of significance in mathematics but it is also directly related to many physical properties. In this respect, Wang, Zhang and Zhang [12] presented a rigorous derivation from the Doi–Onsager molecular theory to the Ericksen–Leslie theory. Based on the same idea, Li, Wang and Zhang [26] rigorously derived the Ericksen–Leslie theory starting from the molecular-based Q-tensor theory for nematic phases. By the Hilbert expansion, Li and Wang [27] rigorously justified the connection between the inertial Qian–Sheng theory and the full Ericksen–Leslie theory. Based on the moment method, Degond, Frouvelle and Liu [28] presented the convergence of the Doi–Onsager theory to the Ericksen–Leslie theory. These works were mainly concerned with the hydrodynamic theory of nematic phases. For a detailed survey on modeling, analysis and computation of liquid crystals, the work of [29] can be refered to.
To our knowledge, Q-tensor hydrodynamics had not be applied to smectic phases prior to this work. In this paper, by using the Bingham closure, we aimed to derive new Q-tensor hydrodynamics that could describe the essential physical properties of smectic liquid crystal flows. Meanwhile, the new reduced Q-tensor model would be computable for future numerical works.

2. Derivation of Q -Tensor Hydrodynamics for Smectic Liquid Crystals

Before deriving the new Q-tensor hydrodynamics for smectic phases, we introduce the Bingham closure and the previous hydrodynamics for nematic phases established in [4].
For a given density function f ( x , m ) satisfying
S 2 f ( x , m ) d m = c ( x ) , S 2 ( m m 1 3 I ) f ( x , m ) d m = c ( x ) Q ( x ) ,
the Bingham closure uses the quasi-equilibrium distribution [4],
f Q = c ( x ) Z Q exp ( B Q ( x ) : m m ) , Z Q = S 2 exp ( B Q ( x ) : m m ) d m ,
to approximate f. Notice that if the density c ( x ) = c is constant, then the state is the nematic phase; if c ( x ) is varied, then the state is the smectic phase. Again, B Q depends on Q and can be determined by the following relation (see [4,26,30])
S 2 ( m m 1 3 I ) f Q d m = Q .
It can be seen from [26,30] that B Q can be uniquely determined for Q and satisfies the physical constraint. Then, the fourth-order moment and the sixth-order moment of f can be viewed as functions of Q.
Now, we introduce the molecular-based Q-tensor dynamical model for nematic liquid crystals presented in [4]. For a given free energy functional F ( Q , Q ) , we define
μ Q = δ F ( Q , Q ) δ Q .
Here the energy functional F ( Q , Q ) , derived from Onsager’s molecular theory, is given by ([4])
F ( Q , Q ) = F b ( Q ) + F e ( Q , Q ) ,
where the bulk energy F b and the elastic distortion energy F e are respectively expressed as
F b ( Q ) = L 0 Ω ln Z Q + Q : B Q α 2 | Q | 2 d x , F e ( Q , Q ) = ε 2 Ω { L 1 | Q | 2 + L 2 i ( Q i k ) j ( Q j k ) + i ( Q j k ) j ( Q i k ) + L 3 | Q ( 4 ) | 2 + L 4 i ( Q i k l m ( 4 ) ) j ( Q j k l m ( 4 ) ) + i ( Q j k l m ( 4 ) ) j ( Q i k l m ( 4 ) ) + L 5 i ( Q i j k l ( 4 ) ) j ( Q k l ) } d x ,
where the coefficients L i ( i = 0 , , 5 ) are explicitly expressed as functions of molecular parameters, and Q ( 4 ) = Q ( 4 ) ( Q ) is the fourth order symmetric traceless moment of f Q . The parameter ε appears in F e due to the fact that the ratios between the coefficients in F e and the ones in F b are at the order of the square of the molecule length.
We introduce two linear operators
M Q ( A ) = 1 3 A + Q · A A : M Q ( 4 ) , N Q ( A ) α β = i { [ γ ( M α β k l ( 4 ) δ i j 1 3 δ α β Q k l δ i j ) + ( γ γ ) ( M α β k l i j ( 6 ) 1 3 δ α β M k l i j ( 4 ) ) ] A k l } ,
where the fourth-order moment M Q ( 4 ) and the sixth-order moment M Q ( 6 ) are respectively defined as follows:
M Q ( 4 ) = S 2 m m m m f Q d m , M Q ( 6 ) = S 2 m m m m m m f Q d m .
Based on Doi–Onsager’s dynamical theory, the molecular-based Q-tensor model for the nematic phase [4] is given by
Q t + v · Q = ε D e N Q ( μ Q ) 2 D e M Q ( μ Q ) + M Q T ( μ Q ) + M Q ( v ) + M Q T ( v ) , v t + v · v = p + γ R e Δ v + 1 γ 2 R e · ( D : M Q ( 4 ) )
+ 1 γ D e R e 2 · M Q ( μ Q ) + μ Q : Q ,
· v = 0 ,
where the term N Q ( μ Q ) represents the translational diffusion. The system (15)–(17) obeys the following basic energy dissipative law (see [4]):
d d t 1 2 Ω | v | 2 d x + 1 γ R e D e F ( Q , Q ) = Ω ( γ R e | v | 2 + 1 γ 2 R e D : M Q ( 4 ) : D ε ( 1 γ ) R e D e 2 μ Q : N ( μ Q ) + 4 ( 1 γ ) R e D e 2 μ Q : M Q ( μ Q ) ) d x .

2.1. From the Doi–Onsager Model to the New Q-Tensor Dynamical Model

The energy functional of the molecular-based Q-tensor static model in a modified version from the works [4,31] is
F [ c , Q ] k B T c 0 = F bulk + F elasitc , 1 + F elasitc , 2 ,
which is measured by the product of the Boltzmann constant k B and the absolute temperature T, and c = c ( x ) is the density of molecules. In (18), F bulk contains the entropy and the quadratic terms of the order parameters, F elasitc , 1 and F elasitc , 2 contain the first derivative terms and the second derivative terms of the order parameters. They are respectively given by
F bulk = Ω c ln c + B Q : Q ln Z + π c 0 L 3 η 2 A 1 c 2 + A 2 | c Q | 2 d x , F elasitc , 1 = π c 0 L 5 η 2 Ω [ G 1 | c | 2 + G 2 | ( c Q ) | 2 + G 3 i ( c Q i j ) j c + G 4 i ( c Q i k j ( c Q j k ) ] d x , F elasitc , 2 = π c 0 L 7 η 2 Ω [ H 1 | 2 c | 2 + H 2 | 2 ( c Q ) | 2 + H 3 i j ( c Q i j ) k k ( c ) + H 4 i k ( c Q i l ) j k ( c Q j l ) + H 5 i j ( c Q i j ) k l ( Q k l ) ] d x ,
where the constant coefficients A i , G i , H i 1 only depend on the dimensionless molecular parameter η = D / L , where D and L are the molecular radius and molecular length, respectively. The coefficient c 0 is the average density of molecules in the region Ω , which is defined by
c 0 = 1 | Ω | Ω c ( x ) d x .
Let x * = x / L 0 , where L 0 is the typical size of the flow region, then the nondimensional form of (18) can be expressed by
F [ c , Q ] = ε 3 2 Ω { c ( ln c + B Q : Q ln Z Q ) + α 2 ( A 1 c 2 + A 2 | c Q | 2 ) ε 2 α G 1 | c | 2 + G 2 | ( c Q ) | 2 + G 3 i ( c Q i j ) j c + G 4 i ( c Q i k ) j ( c Q j k ) + ε 2 2 α [ H 1 | 2 c | 2 + H 2 | 2 ( c Q ) | 2 + H 3 i j ( c Q i j ) k k ( c ) + H 4 i k ( c Q i l ) j k ( c Q j l ) + H 5 i j ( c Q i j ) k l ( c Q k l ) ] } d x .
A direct calculation of the first variational derivative of (19) yields
δ F δ ( c Q ) = ε 3 2 { B Q A 2 α c Q ε 2 α [ 2 G 2 Δ ( c Q ) + G 3 ( i j c 1 3 δ i j Δ c ) + G 4 ( Q j k , j i + Q j i , j k ) ] + ε 2 2 α [ 2 H 2 Δ 2 ( c Q ) + H 3 ( i j Δ c 1 3 δ i j Δ 2 c ) + H 4 ( j i Δ ( c Q j l ) + j l Δ ( c Q j i ) ) + H 5 Q k l , k l i j + Q k l , k l j i ] } ,
δ F δ c = ε 3 2 { ln c + 1 ln Z Q + A 1 α c ε 2 α 2 G 1 Δ c + G 3 ( c Q i j ) , i j + ε 2 2 α 2 H 1 Δ 2 c + H 3 i j Δ ( c Q i j ) } .
In light of the above two variational derivatives, (20) and (21), we can define
μ ^ = μ c Q : m m + μ c , μ c Q = δ F δ ( c Q ) , μ c = δ F δ c .
For convenience, we also define the following linear operators:
S Q ( b ) = i ( γ γ ) c M i j ( 2 ) + c γ δ i j j b ,
B Q ( A ) = i ( γ γ ) c M i j k l ( 4 ) + γ δ i j c M k l ( 2 ) j A k l ,
K Q ( b ) α β = i ( γ γ ) c M α β i j ( 4 ) 1 3 c δ α β M i j ( 2 ) + γ c δ i j Q α β j b ,
N Q ( A ) α β = i { [ ( γ γ ) c M α β k l i j ( 6 ) 1 3 δ α β c M k l i j ( 4 )
+ γ c M α β k l ( 4 ) δ i j 1 3 δ α β c Q k l δ i j ] j A k l } .
Lemma 1.
For any vector field u defined on S 2 , it holds that
S 2 ( m m 1 3 I ) R · ( u f ) d m = m × u m + m m × u f .
The proof of the lemma, which can be found in [3], is easy.
Based on the Doi–Onsager dynamical model (4)–(6), we are now in a position to derive the new Q-tensor dynamical model for smectic liquid crystals. We first derive the time evolution equation for c Q . Multiplying Equation (4) by m m 1 3 I and integrating over a unit ball S 2 , and using the Bingham closure, we have
t c Q + v · c Q = ε D e S 2 ( m m 1 3 I ) · γ m m + γ ( I m m ) · ( f Q μ ^ ) d m I 1 + 1 D e S 2 ( m m 1 3 I ) R · ( f Q R μ ^ ) d m I 2 S 2 ( m m 1 3 I ) R · ( m × κ · m f Q ) d m I 3 .
By a straightforward calculation, we infer that
I 1 = ε D e S 2 ( m m 1 3 I ) · γ m m + γ ( I m m ) · ( f Q μ ^ ) d m = ε D e · S 2 ( γ γ ) ( m m 1 3 I ) m m + γ ( m m 1 3 I ) I · ( μ c ) f Q d m + ε D e · S 2 ( γ γ ) ( m m 1 3 I ) m m + γ ( m m 1 3 I ) I · ( μ c Q ) f Q d m = ε D e N Q μ c Q + K Q μ c .
According to Lemma 1, taking u = μ c Q : R ( m m ) , we have
I 2 = m × μ c Q : R ( m m ) m + m m × μ c Q : R ( m m ) f Q = 2 D e 2 3 c μ c Q + μ c Q · c Q + c Q · μ c Q 2 μ c Q : c M Q ( 4 ) = 2 D e c M Q ( μ c Q ) + M Q T ( μ c Q ) .
Similarly, by Lemma 1 and taking u = m × κ · m , we obtain
I 3 = m × ( m × κ · m ) m + m m × ( m × κ · m ) f Q = 2 ( m × κ · m ) m m ( κ · m ) m m ( κ · m ) f Q = 2 D : m m m m f Q D · m m f Q + Ω · m m f Q m m f Q · ( D + Ω ) = c 2 D : M Q ( 4 ) D · ( Q + 1 3 I ) + Ω · ( Q + 1 3 I ) ( Q + 1 3 I ) · ( D + Ω ) = c M Q ( v ) + M Q T ( v ) .
Therefore, substituting (28)–(30) into (27) we obtain the time evolution equation of c Q
t c Q + v · c Q = ε D e N Q μ c Q + K Q μ c 2 D e c M Q ( μ c Q ) + M Q T ( μ c Q ) + c M Q ( v ) + M Q T ( v ) .
Replacing f with f Q in the elastic stress τ e of (5), it follows that
m m × R μ ^ f Q = 2 c M Q ( μ c Q ) .
In a similar way, the force F e can be also rewritten as
μ ^ f Q = c μ c c Q : μ c Q .
Thus, from Equation (5) we know that the evolution equation of the fluid velocity v is given by
v t + v · v = p + γ R e Δ v + 1 γ 2 R e · D : c M Q ( 4 ) s + 1 γ D e R e 2 · c M Q μ c Q c μ c + c Q : μ c Q .
Finally, from Equation (4) we derive the evolution equation of the molecular density c ( x ) . Integrating over a unit ball S 2 for two sides of (4) gives
c t + v · c = ε D e S 2 · γ m m + γ ( I m m ) · ( f Q μ ^ ) d m J 1 + 1 D e S 2 R · ( f Q R μ ^ ) d m J 2 S 2 R · ( m × κ · m f Q ) d m J 3 .
A straightforward calculation leads to
J 1 = ε D e i ( γ γ ) c M i j k l ( 4 ) + γ δ i j c M k l ( 2 ) j ( μ c Q ) k l + ( γ γ ) c M i j ( 2 ) + γ c δ i j j μ c = ε D e B Q ( μ c Q ) + S Q ( μ c ) .
Note that an important property of the rotational operator R is the formula of integration by parts [5],
S 2 A ( m ) R B ( m ) d m = S 2 [ R A ( m ) ] B ( m ) d m .
Then, it is easy to see from (35) that
J 2 = 0 , J 3 = 0 .
Thus, plugging (34) and (36) into (33), the evolution equation of c ( x ) is given by
c t + v · c = ε D e B Q ( μ c Q ) + S Q ( μ c ) .
To sum up, we obtain the new molecular-based Q-tensor dynamical system of smectic liquid crystals:
t c Q + v · c Q = ε D e N Q μ c Q + K Q μ c 2 D e c M Q ( μ c Q ) + M Q T ( μ c Q ) + c M Q ( v ) + M Q T ( v ) ,
c t + v · c = ε D e B Q ( μ c Q ) + S Q ( μ c ) ,
v t + v · v = p + γ R e Δ v + 1 γ 2 R e · D : c M Q ( 4 ) + 1 γ D e R e 2 · c M Q μ c Q c μ c + c Q : μ c Q ,
· v = 0 ,
where the operator M Q is expressed by (14), the other operators S Q , B Q , K Q , N Q are expressed by (22)–(26), respectively.

2.2. Energy Dissipation Law

One of the most important properties of the Bingham closure approximation is that it maintains the energy dissipation of the kinetic model. We derived the energy dissipation relation of the new molecular-based Q-tensor system (37)–(40). Combining the first variational derivative (20)–(21) and the dynamical system (37)–(40), it follows that
d d t F [ c , Q ] = Ω δ F δ ( c Q ) : ( c Q ) t + δ F δ c c t d x = Ω { ε D e μ c Q : N Q μ c Q + K Q μ c 2 D e μ c Q : c M Q ( μ c Q ) + M Q T ( μ c Q ) + μ c Q : c M Q ( v ) + M Q T ( v ) v · c Q + ε D e μ c B Q ( μ c Q ) + S Q ( μ c ) D e ε v · c } d x ,
d d t Ω 1 2 | v | 2 d x = Ω { γ R e | D | 2 1 γ 2 R e D : c M Q ( 4 ) : D + 1 γ D e R e 2 c M Q μ c Q : v c v · μ c + c Q : ( v · ) μ c Q } d x .
It is easy to see that M Q ( A ) : B = M Q ( B ) : A , if A or B is a symmetric tensor. Then, we deduce from (41) and (42) that
d d t 1 D e R e F [ c , Q ] + Ω 1 2 ( 1 γ ) | v | 2 d x = Ω { ε R e D e 2 μ c Q : N Q μ c Q + K Q μ c 4 R e D e 2 μ c Q : c M Q μ c Q + 1 R e D e μ c Q : 2 c M Q ( v ) v · c Q γ R e ( 1 γ ) | D | 2 1 2 R e D : c M Q ( 4 ) : D + 1 R e D e 2 c M Q μ c Q : v c v · μ c + c Q : ( v · ) μ c Q + ε R e D e 2 μ c B Q ( μ c Q ) + S Q ( μ c ) D e ε v · c } d x = Ω { ε R e D e 2 μ c Q : N Q μ c Q 4 R e D e 2 μ c Q : c M Q μ c Q γ R e ( 1 γ ) | D | 2 1 2 R e D : c M Q ( 4 ) : D + ε R e D e 2 μ c B Q ( μ c Q ) + μ c S Q μ c + μ c Q : K Q μ c } d x .
Note that j μ ^ = j μ c Q : m m + j μ c , we have
0 Ω μ ^ · ( γ γ ) m m + γ I · μ ^ f Q d x = Ω { S 2 i μ c ( γ γ ) m i m j + γ δ i j j μ c f Q d m + S 2 i μ c ( γ γ ) m i m j + γ δ i j ( m k m l 1 3 δ k l ) j ( μ c Q ) k l f Q d m + S 2 i ( μ c Q ) α β ( m α m β 1 3 δ α β ) ( γ γ ) m i m j + γ δ i j j μ c f Q d m + S 2 i ( μ c Q ) α β ( m α m β 1 3 δ α β ) ( γ γ ) m i m j + γ δ i j ( m k m l 1 3 δ k l ) j ( μ c Q ) k l f Q d m } d x = Ω { μ c i ( γ γ ) c M i j ( 2 ) + γ δ i j c j μ c + μ c i ( γ γ ) c M i j k l ( 4 ) + γ δ i j c M k l ( 2 ) j ( μ c Q ) k l + ( μ c Q ) α β i ( γ γ ) c M α β i j ( 4 ) 1 3 c δ α β M i j ( 2 ) + γ c δ i j Q α β j μ c + ( μ c Q ) α β i { [ ( γ γ ) ( c M α β k l i j ( 6 ) 1 3 δ α β c M k l i j ( 4 ) ) + γ ( c M α β k l ( 4 ) δ i j 1 3 δ α β c Q k l δ i j ) ] j ( μ c Q ) k l } } d x = Ω μ c S Q ( μ c Q ) + μ c Q B Q ( μ c Q ) + μ c Q : K Q ( μ c ) + μ c Q : N Q ( μ c Q ) d x .
Therefore, from (43) and (44), the system (37)–(40) fulfills the following energy dissipation law:
d d t Ω 1 2 | v | 2 d x + 1 γ D e R e F ( c , Q ) = Ω { γ R e D : D + 1 γ 2 R e D : c M Q ( 4 ) : D + 4 ( 1 γ ) D e 2 R e μ c Q : c M Q μ c Q + ε ( 1 γ ) D e 2 R e μ ^ · ( γ γ ) m m + γ I · μ ^ f Q } d x .
Similarly to system (15)–(17), the parameters in the system (37)–(40) have clear physical significance but not are phenomenological. The coefficients A i , G i and H i are also explicitly calculated in terms of physical molecular parameters. Another important feature of the new system (37)–(40) is that the translational and rotational diffusions are still maintained.

3. Discussion on the New Q -Tensor Model

In this section, we give some special cases for the molecular-based Q-tensor dynamical system (37)–(40).
Case I. If the order parameter tensor Q and the molecular density c ( x ) are all equal to zero, i.e., Q ( x , t ) = 0 , c ( x , t ) = 0 , then the liquid crystal system will become the isotropic phase. At the moment, the system (37)–(40) is degenerated into
v t + v · v = p + γ R e Δ v ,
· v = 0 .
The system (46) and (47) is just the standard incompressible Navier–Stokes equation for ordinary isotropic liquids. Here γ R e represents the viscous coefficient of isotropic flows.
Case II. If the order parameters Q and c satisfy
Q ( x , t ) = 0 , c ( x , t ) 0 ,
then it follows that
F [ c ] = ε 3 2 c ln c 4 π + 1 2 α A 1 c 2 ε G 1 | c | 2 + ε 2 H 1 | 2 c | 2 d x , μ c = δ F δ c = ε 3 2 ln c 4 π + 1 + α A 1 c + ε G 1 Δ c + ε 2 H 1 Δ 2 c , μ c Q = 0 .
Consequently, the system (37)–(40) can be reduced to
v t + v · v = p + γ R e Δ v + 1 γ 2 R e · ε 2 D e K 0 ( μ c ) + 1 3 c D 1 γ D e R e c μ c ,
c t + v · c = ε D e S 0 ( μ c ) ,
· v = 0 ,
where K 0 ( b ) = K Q ( b ) | Q = 0 , S 0 ( b ) = S Q ( b ) | Q = 0 . We further simplify the above system (48). Note that B Q = 0 , Z Q = 4 π if Q = 0 . Let m = ( sin θ cos φ , sin θ sin φ , cos θ ) T , then we have
M i j ( 2 ) | Q = 0 = 1 4 π S 2 m i m j d m = 1 4 π 0 2 π d φ 0 π d θ sin θ m i m j = 1 3 δ i j ,
M α β i j ( 4 ) | Q = 0 = 1 4 π S 2 m α m β m i m j d m = 1 4 π S 2 1 7 ( m α m β δ i j ) s y m 1 35 ( δ α β δ i j ) s y m d m = 1 4 π 2 7 · 4 π 3 ( δ α β δ i j ) s y m 4 π 35 ( δ α β δ i j ) s y m = 1 15 ( δ α β δ i j ) s y m .
Thus, using the definition of S Q ( b ) and K Q ( b ) , it can be seen from (51) and (52) that
S 0 ( b ) = i ( γ γ ) c M i j ( 2 ) | Q = 0 + c γ δ i j j b = ( 1 3 γ + 2 3 γ ) i ( c j b ) δ i j = def γ ˜ · ( c b ) , K 0 ( b ) α β = i ( γ γ ) c M α β i j ( 4 ) | Q = 0 1 3 c δ α β M i j ( 2 ) | Q = 0 j b = i ( γ γ ) 1 15 ( δ α β δ i j ) s y m 1 9 δ α β δ i j c j b = def γ ^ 1 ( δ α i δ β j + δ α j δ β i ) + γ ^ 2 δ α β δ i j i ( c j b ) = γ ^ 1 α ( c β b ) + β ( c α b ) + γ ^ 2 δ α β · ( c b ) ,
where γ ˜ = 1 3 γ + 2 3 γ , γ ^ 1 = 1 15 ( γ γ ) , γ ^ 2 = 2 45 ( γ γ ) . Let
τ α β ε = γ ^ 1 α ( c β μ c ) + β ( c α μ c ) + γ ^ 2 δ α β · ( c μ c ) + 2 3 D e ε c D α β ,
then the system (48)–(50) becomes the following more simplified incompressible system for the fluid velocity v ( x , t ) and the molecular density c ( x , t ) :
v t + v · v = p + γ R e Δ v + ( 1 γ ) ε 4 D e R e · τ ε 1 γ D e R e c μ c , c t + v · c = γ ˜ ε D e · ( c μ c ) , · v = 0 .
Case III. Suppose that the molecular density c ( x , t ) is constant (particularly c ( x , t ) = 1 ) and Q ( x , t ) 0 , then the system will become the hydrodynamics of the nematic phases. Consequently, we have
μ ^ = μ Q : m m , μ c = 0 , μ Q = δ F δ Q .
The non-dimensional energy functional (19) can be reduced to
F [ Q ] = Ω { B Q : Q ln Z Q + α 2 A 2 | Q | 2 ε 2 α G 2 | Q | 2 + G 4 i Q i k j Q j k + ε 2 2 α H 2 | 2 Q | 2 + H 4 i k Q i l j k Q j l + H 5 i j Q i j k l Q k l } d x .
It should be noted that there some second derivative terms of the order parameter tensor Q appear in (53). For this reason we needed to truncate at the fourth order derivatives of the Taylor expansion for the molecular order parameter f ( x , m ) with regards to x , when deriving the molecular-based Q-tensor static model for smectic liquid crystals using the Bingham closure. If we only considered describing the nematic liquid crystals, then we did not need such high-order Taylor expansions. Therefore, after removing the second order derivative terms of Q in (53), the corresponding energy functional is
F [ Q ] = Ω B Q : Q ln Z Q + α 2 A 2 | Q | 2 ε 2 α G 2 | Q | 2 + G 4 i Q i k j Q j k d x .
The functional (54) is a special case of the free energy functional (13) with the coefficients L 0 = 1 , L 3 = L 4 = L 5 = 0 . Hence, under the condition of incompressible flows, the system (37)–(40) is degenerated into the molecular-based Q-tensor dynamical model (15)–(17) that describes the nematic liquid crystal flows.

4. Conclusions

In this paper, based on the Doi–Onsager molecular theory, a novel Q-tensor hydrodynamic model for smectic phases was derived using the Bingham closure. The model demonstrated some new features: its coefficients were all expressed as those in the molecular model, and the energy dissipation was maintained. The hydrodynamics could be further reduced to the version of nematic liquid crystals if the density was constant. The new model not only described many interesting physical phenomena of smectic liquid crystals, but it also included complex structures of equations. More importantly, the new model will be a computable model when considering the quasi-entropy proposed by recent work [32], since quasi-entropy is an elementary function that maintains the same asymptotic behavior and underlying physics of the Bingham closure. Therefore, we aim to investigate this model from numerical computations and well-posedness analysis in future works. Furthermore, exploring the connection between the new Q-tensor hydrodynamics and the Ericksen–Leslie model with density variations will be a rather challenging problem.

Author Contributions

The contributions of all authors are equal. All authors have read and agreed to the published version of the manuscript.

Funding

Sirui Li is supported by the Growth Foundation for Youth Science and Technology Talent of the Educational Commission of Guizhou Province of China under grant No. [2021]087.

Data Availability Statement

Not applicable.

Acknowledgments

We would like to thank the editors and the referees for their valuable comments and suggestions for improving our paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. De Gennes, P.G. The Physics of Liquid Crystals; Clarendon Press: Oxford, UK, 1974. [Google Scholar]
  2. Onsager, L. The effects of shape on the interaction of colloidal particles. Ann. N. Y. Acad. Sci. 1949, 51, 627–659. [Google Scholar] [CrossRef]
  3. Zhang, P. A molecular kinetic theory of inhomogeneous liquid crystal flow and the small Deborah number limit. Methods Appl. Anal. 2006, 13, 181–198. [Google Scholar]
  4. Han, J.; Luo, Y.; Wang, W.; Zhang, P.; Zhang, Z. From microscopic theory to macroscopic theory: Systematic study on modeling for liquid crystals. Arch. Ration. Mech. Anal. 2015, 215, 741–809. [Google Scholar] [CrossRef]
  5. Doi, M.; Edwards, S.F. The Theory of Polymer Dynamics; Oxford University Press: Oxford, UK, 1986. [Google Scholar]
  6. Forest, M.G.; Wang, Q.; Zhou, R. The weak shear kinetic phase diagram for nematic polymers. Rheol. Acta 2004, 43, 17–37. [Google Scholar] [CrossRef]
  7. Forest, M.G.; Wang, Q.; Zhou, R. The flow-phase diagram of Doi-Hess theory for sheared nematic polymers II: Finite shear rates. Rheol. Acta 1961, 41, 80–93. [Google Scholar] [CrossRef]
  8. Larson, R.G.; Öttinger, H. Effect of molecular elasticity on out-of-plane orientitations in shearing flows of liquid-crystalline polymers. Macromolecules 1991, 24, 6270–6282. [Google Scholar] [CrossRef]
  9. Wang, Q.; Weinan, E.; Liu, C.; Zhang, P. Kinetic theories for flows of nonhomogeneous rodlike liquid crystalline polymers with a nonlocal intermolecular potential. Phys. Rev. E 2002, 65, 051504. [Google Scholar] [CrossRef]
  10. Marrucci, G.; Greco, F. The elastic constants of Maier-Saupe rodlike molecule nematics. Mol. Cryst. Liq. Cryst. 1991, 206, 17–30. [Google Scholar] [CrossRef]
  11. Wang, Q. A hydrodynamic theory for solutions of nonhomogeneous nematic liquid crystalline polymers of different configurations. J. Chem. Phys. 2002, 116, 9120–9136. [Google Scholar] [CrossRef]
  12. Wang, W.; Zhang, P.; Zhang, Z. The small Deborah number limit of the Doi-Onsager equation to the Ericksen-Leslie equation. Comm. Pure Appl. Math. 2015, 68, 1326–1398. [Google Scholar] [CrossRef]
  13. Frank, F.C. On the theory of liquid crystals. Trans. Soc. Rheol. 1961, 5, 22–34. [Google Scholar]
  14. Oseen, C.W. The theory of liquid crystals. Trans. Faraday Soc. 1933, 29, 883–899. [Google Scholar] [CrossRef]
  15. De Gennes, P. Short range order effects in the isotropic phase of nematics and cholesterics. Mol. Cryst. Liq. Cryst. 1971, 12, 193–214. [Google Scholar] [CrossRef]
  16. Chen, J.; Lubensky, T.C. Landau-Ginzburg mean-field theory for the nematic to smectic-C and nematic to smectic-A phase transitions. Phys. Rev. A 1976, 14, 1202–1207. [Google Scholar] [CrossRef]
  17. Ericksen, J. Conservation laws for liquid crystals. Trans. Soc. Rheol. 1961, 5, 22–34. [Google Scholar] [CrossRef]
  18. Leslie, F.M. Some constitutive equations for liquid crystals. Arch. Ration. Mech. Anal. 1968, 28, 265–283. [Google Scholar] [CrossRef]
  19. Ahmadi, G. A continuum theory of smectic A liquid crystals. J. Rheol. 1982, 26, 535–556. [Google Scholar] [CrossRef]
  20. Martin, P.C.; Parodi, O.; Pershan, P.S. Unified hydrodynamic theory for crystals, liquid crystals and normal fluids. Phys. Rev. A 1972, 6, 2401–2420. [Google Scholar] [CrossRef]
  21. Weinan, E. Nonlinear continuum theory of smectic-A liquid crystals. Arch. Rational. Mech. Anal. 1997, 137, 343–360. [Google Scholar] [CrossRef]
  22. Stewart, I.W. Dynamic theory for smectic A liquid crystals. Continuum Mech. Thermodyn. 2007, 18, 343–360. [Google Scholar] [CrossRef]
  23. Mottram, N.J.; Newton, C. Introduction to Q-Tensor Theory; Research Report; Department of Mathematics, University of Strathclyde: Glasgow, UK, 2004; Volume 10. [Google Scholar]
  24. Beris, A.N.; Edwards, B.J. Thermodynamics of Flowing Systems with Internal Microstructure; Oxford Engrg, Sci. Ser. 36; Oxford University Press: Oxford, UK; New York, NY, USA, 1994. [Google Scholar]
  25. Qian, T.; Sheng, P. Generalized hydrodynamic equations for nematic liquid crystals. Phys. Rev. E 1998, 58, 7475–7485. [Google Scholar] [CrossRef]
  26. Li, S.; Wang, W.; Zhang, P. Local posedness and small Deborah limit of a molecule-based Q-tensor system. Discrete. Contin Dyn. Syst.-B 2015, 20, 2611–2655. [Google Scholar] [CrossRef]
  27. Li, S.; Wang, W. Rigorous justification of the uniaxial limit from the Qian-Sheng inertial Q-tensor theory to the Ericksen-Lesile theory. SIAM J. Math. Anal. 2020, 52, 4421–4468. [Google Scholar] [CrossRef]
  28. Degond, P.; Frouvelle, A.; Liu, J.G. From kinetic to fluid models of liquid crystals by the moment method. Kinet. Relat. Model. 2022, 15, 417–465. [Google Scholar] [CrossRef]
  29. Wang, W.; Zhang, L.; Zhang, P. Modeling and computation of liquid crystals. Acta Numer. 2021, 30, 765–851. [Google Scholar] [CrossRef]
  30. Ball, J.M.; Majumdar, A. Nematic liquid crystals: From Maier-Saupe to a continuum theory. Mol. Cryst. Liq. Cryst. 2010, 525, 1–11. [Google Scholar] [CrossRef]
  31. Mei, S.; Zhang, P. On a molecular based Q-tensor model for liquid crystals with density variations. Multiscale Model. Simul. 2015, 13, 977–1000. [Google Scholar] [CrossRef]
  32. Xu, J. Quasi-entropy by log-determinant covariance matrix and application to liquid crystals. Phys. Nonlinear Phenom. 2022, 435, 133308. [Google Scholar] [CrossRef]
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Feng, X.; Li, S. A Molecular-Based Q-Tensor Hydrodynamic Theory of Smectic Liquid Crystals. Axioms 2022, 11, 495. https://doi.org/10.3390/axioms11100495

AMA Style

Feng X, Li S. A Molecular-Based Q-Tensor Hydrodynamic Theory of Smectic Liquid Crystals. Axioms. 2022; 11(10):495. https://doi.org/10.3390/axioms11100495

Chicago/Turabian Style

Feng, Xinxin, and Sirui Li. 2022. "A Molecular-Based Q-Tensor Hydrodynamic Theory of Smectic Liquid Crystals" Axioms 11, no. 10: 495. https://doi.org/10.3390/axioms11100495

APA Style

Feng, X., & Li, S. (2022). A Molecular-Based Q-Tensor Hydrodynamic Theory of Smectic Liquid Crystals. Axioms, 11(10), 495. https://doi.org/10.3390/axioms11100495

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