NCCT for Micropolar Solid and Fluid Media Based on Internal Rotations and Rotation Rates with Rotational Inertial Physics: Model Problem Studies ()
1. Introduction, Literature Review and Scope of Current Work
In recent papers, Surana et al. [1] [2] have shown that micropolar NCCT for solid and fluid media based on internal rotations (due to deformation gradient tensor) and internal rotation rates (due to velocity gradient tensor) are thermodynamically and mathematically consistent. All other micropolar non-classical continuum theories based on Cosserat rotations and Cosserat rotation rates, those considering internal and Cosserat rotations and rotation rates, couple stress NCCT utilizing polynomial or potential functions approach for deriving constitutive theories are either in violation of thermodynamic consistency and/or mathematical consistency. In reference [3] Surana et al. extended the micropolar NCCT of reference [1] based on internal rotations for solid medium to account for rotational inertial physics. This leads to modification of balance of angular momenta (BAM) and balance of moment of moments (BMM) balance laws of reference [1] . Additionally, satisfying entropy inequality requires incorporating a constraint equation in the mathematical model consisting of CBL and the constitutive theories. Authors showed (without model problem studies) that this mathematical model permits co-existence of translational and rotational waves. In reference [4] Surana et al. extended the micropolar NCCT of reference [2] for fluids to include rotational inertial physics. This also resulted in modification of the CBL presented in reference [2] . Authors in references [3] [4] showed that the requirement of satisfying entropy inequality necessitated that a constraint equation be considered as part of the mathematical model based on CBL and the constitutive theories. Authors discussed (without model problem studies) similarity of BAM in rotation rates in the case of micropolar fluids [4] with the BAM for micropolar solids in rotations in reference [3] , but lack of rotational waves in micropolar fluids due to absence of stiffness (due to lack of strain physics). Influence of rotational inertial physics on the evolution of the IVPs in general resulting due to micropolar NCCT was not discussed. Unfortunately, there are not many published works on micropolar NCCT that consider rotational inertial physics. References [5] [6] consider flow of micropolar fluids through porous medium and dynamics of hydro-magnetic flow of micropolar fluids.
The objective of this paper is to present model problem studies for micropolar solids and micropolar fluids using the mathematical model consisting of the CBL and the constitutive theories for micropolar non-classical solids and fluids with rotational inertial physics presented in references [3] [4] . Model problem studies are intentionally kept simple (one dimensional as far as possible) so that the contributions of rotational inertial physics can be clearly illustrated and established.
2. Mathematical Models
In the following, we consider complete 3D mathematical models consisting of CBL and the constitutive theories for micropolar non-classical solid and fluids based on internal rotations and internal rotation rates with rotational inertial physics presented by Surana et al. [3] [4] . In both cases, we consider thermoviscous physics in which the dissipation mechanism is due to rate of strain (solids) and strain rate (fluids) as well as due to rate of symmetric part of rotation gradient tensor (solids) as well as due to symmetric part of rotational gradient rates (fluids).
2.1. Micropolar Thermoviscous Elastic Solid (Without Memory)
For small strain, small deformation physics the CBL and the constitutive theories presented in reference [3] modified for dissipation due to rate of strain (CCM) and rate of symmetric part of rotation gradient tensor (micropolar NCCT) are given by
(1)
(2)
(3)
(4)
(5)
(6)
(7)
(8)
In the mathematical model we have a total of 25 equations BLM(3), BAM(3), BMM(3), FLT(1), constitutive theories for
,
,
in a total of 25 dependent variables
,
,
,
,
,
,
hence the model has closure. The last term in the entropy inequality
(9)
must be set to zero to ensure that (9) is always satisfied. Then, in addition to (1)-(8). We must also satisfy
(10)
Equation (10) serves as a constraint on the mathematical Model (1)-(8).
2.2. Micropolar Thermoviscous Fluid Medium
We consider incompressible micropolar thermoviscous fluid medium. The CBL and constitutive theories of reference [4] are considered. These naturally include dissipation mechanism due to classical as well as non-classical micropolar physics. Following reference [4] we have the following for CBL and constitutive theories.
(CM) (11)
(BLM) (12)
(BAM) (13)
(BMM) (14)
(FLT) (15)
(SLT) (16)
(17)
(18)
(19)
(20)
where
(21)
In this mathematical model we have CM(1), BLM(3), BAM(3), BMM(3), FLT(1) and the constitutive theories for
,
,
, a total of 26 equations in 26 dependent variables
,
,
,
,
,
,
,
. The last term in the entropy inequality must be set to zero to ensure that (16) is not violated. Thus
(22)
must serve as constraint, hence constitutes additional condition that must be satisfied to guarantee that the entropy inequality is not violated.
3. Model Problems
In this section, we consider model problems for micropolar solid and fluid media. Model problems are intentionally kept simple so that significant features of micropolar aspects in the presence of rotational inertial physics can be demonstrated.
3.1. Translational and Rotational Waves in Micropolar Solids
The objective of this study is to demonstrate existence of pure rotational waves in micropolar solids with rotational inertial physics. Just like in classical continuum mechanics (CCM) we can show existence, propagation, reflection and interaction of translational (stress) waves, in case of micropolar solids with rotational inertial physics we can additionally show existence, propagation, reflection and interaction of rotational (moment) waves. That is, in micropolar solids with rotational inertial physics, translational and rotational waves coexist. In the following study we consider purely one dimensional case of BLM and BAM. In the studies presented here we show that BLM permits translational waves where as BAM with rotational inertial physics permits pure rotational waves.
Mathematical models
One dimensional form of BLM and balance of angular momenta with rotational inertial physics (Equations (1) and (2)) for small strain small deformation and the associated constitutive theories for thermoviscoelastic (without memory) micropolar solid with dissipation mechanism due to rate of strain (CCM) and due to rate of symmetric part of the rotation gradient tensor are given by (in the absence of body forces, body moments, initial stress and equilibrium Cauchy stress, thus
and using x for x1)
(23)
(24)
(25)
(26)
Substituting (25) in (23) and (26) in (24) we obtain
(27)
(28)
for 1D case and for incompressible physics
, modulus of elasticity, as Poison’s ratio
and
, hence (28) reduces to (but (28) remains same).
(29)
(30)
Remarks
1) Equation (29) describes translational or
stress wave in a viscous elastic medium. In the absence of the last term in (29), (29) represents translational waves in inviscid elastic medium. Dissipation (last term in (29)) results in amplitude decay and base elongation.
2) Equation (30) is analogous to equation (29) and can be obtained b replacing
,
, E, and
in (29) by
,
,
and
. Thus, (30) describes rotational waves i.e., moment (
) wave. The third term in (30) is due to micropolar dissipation physics which would result in amplitude decay and base elongation of
rotational wave.
3) We note the lack of coupling between (29) and (30) due to zero antisymmetric part of the deviatoric Cauchy stress. This is not the case in 2D and 3D applications.
4) Presence of (29) and (30) confirm coexistence of translational and rotational waves in micropolar elastic solids.
Dimensionless form of the mathematical model
We first nondimensionalize (29) and (30) before presenting their solutions. We write (29) and (30) with hat (
) on all quantities indicating that all quantities have their usual dimensions (units)
(31)
(32)
We note that in general the speed of translational wave is different than the speed of rotational wave. However, since in this case (31) and (32) are decoupled, it is possible to nondimensionalize (31) and (32) such that in the dimensionless domain the speed of prorogation of both translation and rotational waves is unity. Obviously in 2D and 3D cases it is not possible to do so. We present details in the following.
BLM: Model TW1
Choosing
,
,
,
,
and
as reference quantities in which
is the speed of translational wave (speed of sound in CCM).
(33)
where
, Reynolds number. If we choose
,
and
as damping coefficient, then (33) can be reduced to
(34)
In (34), the dimensionless speed of sound (based on reference quantities) in an inviscid medium is unity. We use (34) in the numerical studies.
BAM: Using speed of rotational waves as reference velocity: Model RW1
In this case, when nondimensionalizing BAM we choose speed of rotational wave as reference velocity. We also choose some of the same reference quantities as in BLM except we choose
,
,
,
and
,
, then (32) can be written as
(35)
If
,
, then (35) can be reduced to
(36)
in which
is the Reynolds number associated with dissipation due to micropolar physics. If we define
as the dissipation coefficient due to micropolar physics, then (36) can be written as
(37)
Equation (34) and (37) are the simplified dimensionless forms of BLM and BAM describing translational and rotational waves in the presence of dissipation mechanisms. In the numerical studies we use equations (34) and (37).
BAM: Using translation wave speed as reference velocity: Model RW2
In this case we choose
(38)
using these reference values we can nondimensionalize (32) as follows.
(39)
If we choose
,
and define
then (39) can be reduced to
(40)
In this case the speed of dimensionless rotational wave is
which may not be unity. We use (40) in the numerical studies.
Remarks
1) We consider (34), (37) and (40) for presenting numerical studies.
2) Numerical studies are presented using space-time coupled finite element method based on space-time residual functional for a space-time strip with time marching [5] in which space-time local approximation for a space-time element are p-version hierarchical with higher order global differentiability in space and time.
3) We recast (34), (37) and (40) as a system of first order PDEs for convenience of defining BCs and ICs.
First order system of PDEs
Translational (or stress) waves (Model TW1)
We recast (34) as follows
(41)
Rotational (or moment) waves (Model RW1)
We recast (37) as
(42)
Rotational (or moment) waves (Model RW2)
(43)
Numerical studies
Evolutions are computed for models TW1, RW1 and RW2 using space-time coupled finite element method for a space-time strip with time marching [7] . The space-time local approximation for a space-time element is p-version hierarchical with higher order global differentiability in space and time. Discretizations, p-levels and the minimally confirming spaces are chosen such that for each space-time strip the space-time residual functional for the discretization is O(10−8) or lower ensuring that PDEs are satisfied accurately by the computed solution. In the present studies we choose a uniform mesh of 30, nine node p-version hierarchical finite elements with p-levels of 9 in space and time for all elements of the space-time strip. We choose local approximations in space and time of class C11.
A time increment of
is used in all computations. The solution is computed for the first space-time strip (
) and then time marched to obtain evolution for desired value of
(final value of time). Since all evolutions considered here are smooth, choice of solutions of class C11 in space and time suffices. For this choice of the order of the approximation space, the space-time integrals for the discretization of the space-time strip are Riemann. The 30 element uniform discretization with p-levels of 9 for space and time and
for both space and time, yield space-time residual functional values for each space-time strip of the order O(10−8) or lower during the entire evolution, confirming good accuracy of the computed solution. In the following, we present computed evolutions for: translation wave using Model TW1 and rotations wave using Models TR1 and TR2.
Figure 1 shows schematic, space-time domain
, discretization
of the space time domain into space-time strips, discretization
of the first space-time strip
and the details of stress
boundary conditions (for TW1) applied to the last space-time element face at
. In the numerical studies, we choose negative
(compressive) for the applied pulse at
. The details of applying BC at
for mathematical models RW1 and RW2 are similar.
Translational wave: Model TW1
We consider mathematical model TW1. Evolution is computed for 30 time steps for the undamped case,
using space-time coupled finite element formulation based on space-time residual functional. Since in this case, the speed of wave propagation is one, in ten time steps with
, the applied pulse at
reaches the impermeable boundary at
. Due to zero damping, we expect the applied pulse shape to be preserved during evolution. Figure 2(a) shows a plot of
versus x for various values of time. Time steps 2, 6, 10 show the incident compressive pulse entering the spatial domain at
and propagating without amplitude decay or base elongation confirming the absence of numerical dispersion. At the 10th time step the compressive pulse is precisely at
as expected. Reflection of the compressive incident pulse at the impermeable boundary (at
) results in reflected compressive pulse at the 11th time step with double the peak value. Upon further evolution the reflected compressive pulse recovers its peak and base at the 12th time step (not shown). At the 14th time step, we see reflected pulse with the same amplitude and base as the original incident pulse propagating toward the right end. At the 20th time step the reflected compressive pulse reaches the free boundary at
and
Figure 1. Spatial domain, space-time domain, Discretization into space-time strips, First space-time strip, Boundary conditions.
Figure 2. Evolution of deviatoric stress wave:
versus position x (Model TW1).
reflects as a tensile pulse of the same amplitude and base as the original incident pulse propagating towards the boundary at
(time steps 23 and 28 in Figure 2(a)).
In the next study for TW1, we consider dimensionless damping coefficient
. Figure 2(b) shows evolution of
versus x for time steps 2, 6, 10; 11, 14, 18; 23 and 28 (same as in Figure 2(a) for the undamped case). In the incident pulse, amplitude decays and the base elongates during propagation (time steps 2, 6, 10). During reflection (time step 11), the magnitude of the pulse increases but recovers upon further evolution (time step 14). The reflected pulse reaches the free boundary and reflection from the free boundary results in a tensile pulse that propagates to the left. We observe the amplitude decay is more pronounced in the propagating incident pulse. The reflected pulses also shows some amplitude decay, but the base elongation in the reflected pulse is more pronounced. All evolutions are smooth i.e., free of oscillations. Integrated sum of squares of the space-time residual functional are of the order of O(10−8) or lower for each space-time strip ensuring that the computed solution is accurate. Choice of
(order of the approximation space in space and time) is minimally conforming for the first order system of PDEs, hence all space-time integrals for
, discretization of space-time strip
are Riemann. The space-time residual functional of the order O(10−8) ensures that the PDEs are satisfied in point wise sense over each space-time strip discretization
, hence good time accuracy of the evolution is ensured.
Rotational wave: Model RW1
We consider mathematical model RW1 defined by Equation (42). In this case, the equations are nondimensionalized using the speed of rotational wave, hence the dimensionless speed of the rotational wave is one. In the first study, we apply moment pulse i.e.,
pulse of time duration
(similar to
shown in Figure 1(e)). We choose p-levels of 9 in space and time,
and order k of the approximation space in space and time is chosen to be 2, which is minimally confirming for a first order system of PDEs. We choose peak values of
. Evolution is computed using space-time strip with time marching, Integral form in space-time finite element process is based on space-time residual functional. Figure 3(a) shows the evolution for
(undamped case). Incident moment pulse propagates towards the impermeable boundary without amplitude decay and base elongation. Upon reflection at
, the peak amplitude doubles during reflection but the pulse recovers to the original shape in the 12th time steps and continues to propagate toward the free boundary. Reflection from the free boundary results in tensile moment pulse of the same shape as original incident compressive moment pulse that continues to propagate toward the impermeable boundary at
.
Figure 3(b) shows the propagation, reflection and propagation upon reflection of the same compressive
pulse as in Figure 3(a) but in the presence of dissipation. We choose dimensionless dissipation coefficient
. From Figure 3(b) we observe almost the same behavior of the moment pulse as that of
pulse in Figure 2(b). Presence of dissipation resulting in continued amplitude decay and base elongation during the propagation. Amplitude decay is more pronounced for the incident pulse and the base elongation is more significant in the reflected pulse.
Rotational Wave: Model RW2
In this study, we consider mathematical model RW3, Equation (43). In this mathematical model rotational wave speed is
and the dissipation is controlled by the dimensionless dissipation coefficient
. We apply a moment (
) pulse of duration
at the boundary
.
In the first study, we choose
,
and
(undamped case). The wave speed is
, thus the pulse would reach
Figure 3. Evolution of Cauchy Moment wave:
versus position x (Model RW1).
boundary in
time steps. Figure 4(a) show propagation and reflection of the moment pulse. Figure 4(b) show propagation and reflection of the moment pulse for
i.e., at the wave speed of 1.5 when
. Diminished amplitudes and elongated base of the pulse is clearly observed during evolution. This is similar to RW1 (
). In the next study, we apply the same moment pulse as in the previous study but choose
and
, hence rotational wave speed of
, thus the moment pulse will reach the boundary at
in
time steps. Figure 5(a) and Figure 5(b) show propagation and reflection of the moment wave for
and
. The undamped pulse (
) maintains its base and amplitude during the entire evolution and when
progressive amplitude decay and base elongation of the pulse is observed for both values of
.
Figure 4. Evolution of Cauchy Moment wave:
versus position x (Model RW2:
).
We remark that when the rotational wave speed is one, we are able to precisely locate the pulse at the impermeable boundary as the time to reach the impermeable boundary is integer multiple of
. In this case, we see perfect reflection of the
pulse at the 11th time step and the peak value of
doubles, same as in the case of the translational stress wave. However, when rotational wave speed is not an integer multiple of
, precise arrival of
pulse at
cannot be simulated as elapsed time is always an integer multiple of
. For this reason, we are not able to observe 2 (
) in the peak values during the reflection process of the undamped case (Figure 4(a) and Figure 5(a)) when the rotational speed is not one.
Figure 5. Evolution of Cauchy Moment wave:
versus position x (Model RW2:
).
Remarks
1) In the model problem studies presented here the translational and rotational waves are decoupled. In
and
this may not be the case. Since the stress tensor is a function of
and the moment tensor is a function of
, both stress and moment tensor waves are dependent on gradients of displacement. Thus, in
and
we expect interaction between the two waves. That is, as shown here, rotational waves depend upon
and
which in turn influences
, thus
, hence influencing translational waves. The model problem studies in
and
are needed to illustrate this physics.
2) The medium with and without microconstituents is always considered isotropic and homogeneous. Thus, in this micropolar theory there is absence of wave dispersion in the micropolar medium as it is considered to be isotropic and homogeneous. This micropolar theory can only simulate wave propagation physics of translational and rotational waves and their interaction (when it exists) in
and
(future work).
3.2. Micropolar Viscous Fluids: Rotational Inertial Physics
Rotational wave?
The objective of the model problem studies presented in this section is to demonstrate the influence of rotational inertial physics in micropolar fluids in which NCCT is based on internal rotation rate physics. It is well known that in classical or micropolar viscous fluids, translational deviatoric stress waves can not exist due to the absence of elasticity. In this study, we only consider balance of angular momenta in 1D and associated constitutive theories to investigate the influence of rotational inertial physics.
Mathematical model: rotational wave?
The one dimensional form of BAM and the constitutive theory for the Cauchy moment tensor are given by (using
for
).
(44)
(45)
Substituting (45) in (44)
(46)
Remarks
1) We clearly see that this Equation (46) will not permit a moment (
)wave as it is not a wave equation due to the absence of
in place of
.
2) However it is interesting to study its solution for varying
values.
Dimensionless form of the mathematical model
Recasting (46) with hat (
) on all quantities indicating that they have their usual dimensions (or units)
(47)
and using
(48)
We can write (47) as
(49)
If we choose
,
and
then (49) can be reduced to
(50)
in which
, dimensionless dissipation coefficient. PDE (50) can be cast as a system of two first order PDEs by using moment
as auxiliary variable and we have
(51)
PDEs in (51) are helpful in defining BCs and ICs, hence are used in the computations.
Numerical Studies
Rotational wave?
We consider mathematical model (50) or (51). It is evident (50) is not a wave equation as a wave equation must have second order spatial and second order time derivatives of the dependent variable. Thus, (50) can not describe a wave, which is not a surprise, because classical as well as micropolar fluids do not possess elasticity, hence do not have stiffness, thus they cannot support deviatoric Cauchy stress wave or Cauchy moment wave.
Equation (50) is in fact a time dependent diffusion equation in which
is the diffusion coefficient. Thus low values of
correspond to high values of
and vice versa. Details of the schematic, space-time strips and discretization of a space-time strip remain the same as shown in Figure 1. We apply a negative
moment pulse of duration
on the boundary at
such that
;
and
for
(similar to BC shown in Figure 1(e)). We choose a discretization of 30 nine node p-version hierarchical space-time elements with p-levels of nine in space and time. Evolution is computed for 30 time steps using
with local approximation of class C11 in space and time, We choose two values of
. Evolutions for different values of time for the two choices of
are shown in Figure 6(a) and Figure 6(b). For both values of
, the applied pulse progressively diffuses as time elapses. We observe significantly high diffusion of the applied moment pulse for
compared to
as expected. We clearly observe lack of existence and propagation of rotational or moment waves due to lack of elasticity, demonstrating that rotational waves can not exist in micropolar fluids with rotational inertial physics.
Figure 6. Cauchy Moment
versus position
.
3.2.1. Developing Pressure Driven Flow between Parallel Plates and Developing Couette Flow
In this study, we consider developing pressure driven flow between parallel plates and developing Couette flow in which the mathematical models are based on: CCM, NCCM without rotational inertial physics and NCCM with rotational inertial physics. Only the micropolar nonclassical continuum theory based on internal rotation rates is considered [4] . Figure 7 shows a schematic of the dimensionless configuration of parallel plates.
Away from the ends A and B, the physics of the flow is purely one dimensional (in the
direction) velocity
is a dependent variable, velocity
in the
direction is zero regardless of the type of mathematical model hence, Cauchy shear stress is the only non-zero stress (details are given in the following). We present details of three mathematical models (as mentioned above): based
Figure 7. Spatial domain for pressure driven and Couette flow.
on CCM, NCCM without rotational inertial physics and NCCM with rotational inertial physics. We nondimensionalize these mathematical models using the following reference quantities and dimensionless quantities.
(52)
Mathematical model based on CCM: Model A
Balance of linear momenta in the
direction and the constitutive theory for deviatoric Cauchy stress
are the only two equations needed in this case and are given in the following (with usual dimensions for all quantities) for incompressible classical fluid.
(53)
We can nondimensionalize (53) using (52) and if we choose
,
, then the dimensionless form of (53) can be written as
(54)
In case of pressure driven flow,
is given. For non-pressure driven Couette flow,
. In (54), Re is Reynolds number and is given by (based on reference quantities)
. Equation (54) is a first order system of two PDEs in two dependent variables
and
. In this case, flow characteristics depend upon dimensionless parameter Re.
Micropolar NCCT without rotational inertial physics: Model B
In this case, the mathematical model consists of BLM in
direction, balance of angular momenta about the
direction, constitutive theory for symmetric part of deviatoric Cauchy shear stress
, constitutive theory for symmetric Cauchy moment
and definition of angular rotation rate
i.e., rotation rate about
. We can write the following (with their usual dimensions) using the CBL and the constitutive theories given for
by Equations (11)-(22), in the absence of body forces and body moments.
(55)
We can nondimensionalize (55) using (52), and if we choose
,
,
, then (55) can be written as
(56)
in which
, dimensionless material coefficient for micropolar nonclassical physics. Equations (56) is a system of five first order PDEs in five dependent variables:
,
,
,
and
.
Dimensionless parameters Re and
control classical physics and micropolar nonclassical physics.
Micropolar NCCT with rotational inertial physics: Model C
The mathematical model consists of BLM in
direction, BAM about
direction, BMM balance law, constitutive theory for symmetric part of the deviatoric Cauchy shear stress
, constitutive theory for symmetric Cauchy moment
, definition of angular rotation rate
and the constraint equations resulting from the entropy inequality.
Using the mathematical model in
(Equations (11)-(22)) and noting that in this case
,
, we can write the following for CBL and constitutive theories (in the absence of body forces and body moments) and the constraint equation
(57)
(BMM when
) (58)
(constraint equation) (59)
When
(which is the case here), BMM (58) reduces to
(60)
Equation (60) suggests that
or the nonsymmetric moment tensor
is zero, implying no micropolar physics. This of course is erroneous, hence in this case (60) cannot be used as part of the mathematical model. The constraint Equation (59) implies that either
or
for all
. This also is erroneous, thus (58) and (59) cannot be considered as a part of the mathematical model, thus Equation (57) constitutes the mathematical model. We can nondimensionalize (57) using (52), and if we choose
,
,
, then we can obtain the following from (57).
(61)
These are a system of five first order PDEs in five dependent variables:
,
,
,
and
. Re,
and
control classical physics, nonclassical physics without rotational inertial effects and nonclassical physics with rotational inertia effects, respectively. In (61),
and
are defined as
and
.
Numerical studies
We consider two model problems in this section: developing flow between parallel plates and developing Couette flow. In both model problems, we compute solutions for Model A (CCM); Model B, NCCM with micropolar physics but absence of rotational inertial physics; and Model C, micropolar NCCM with rotational inertial physics. In Model A (CCM), the only dimensionless parameter is Re which controls the flow physics. In case of Model B, Re for CCM and
for micropolar physics control the flow physics. In case of Model C, Re,
and
all three control the physics.
is associated with rotational inertial physics. For both model problem studies we present evolutions for different combinations of
and
for a fixed Re to illustrate their influence on flow physics. Computed solutions for Model B and C are always compared with Model A.
We recall that many works published by Surana et al. [8] [9] show that in micropolar NCCT presence of micro constituents offer resistance to flow i.e., increasing value of
results in diminishing flow rate. Material coefficient
, due to rotational inertial physics, related to rotational inertial physics due to microconstituents is also expected provide further resistance to flow over and beyond
. In other words, we expect increasing
as well as increasing
to provide increasing resistance to flow. Since the constitutive theories are not calibrated, we do not know actual values of
and
for various fluids with varying microconstituents. In the model problem studies presented here, we choose
and
such that we can demonstrate their relative influence on flow physics.
Developing flow between parallel plates
Figure 8 shows schematic, space-time strips and a uniform discretization of first space-time strip using 10 p-version hierarchical space-time finite elements. We consider solution of class C00 in space and time with p-level of nine in space and time,
with
. Flow is pressure driven with
.
Case I
In the first study we choose
Model A (CCM),
Model B (NCCM),
;
.
Figures 9(a)-(d) show plots of the evolution of velocity
for
,
,
and
for Models A and B. In case of Model A (CCM) the flow has no resistance due to microconstituents. In model B, with progressively increasing
the flow resistance increases resulting in diminishing flow rate. This holds true for each value of time t. At
, we almost have stationary
Figure 8. Developing flow between parallel plates, Discretization in time, Spatial discretization.
state of the developing flow. Significant reduction in flow rate with increasing
is also obvious from the stationary state in Figure 9(d). In CCM, rotation rates are a free field, hence they do not influence CCT. Figure 10(a) shows evolution of
versus
for Model A (CCM). At
, we almost have stationary sate at which
. Since
is a free field,
,
Figure 9. Developing flow between parallel plates: velocity
versus position
: Models A and B (Case I).
shown in Figure 10(b). Figure 10(c) and Figure 10(d) show evolution of
versus
and
versus
for Model B for
. Comparison of these evolutions with those in Figure 10(a) and Figure 10(b) clearly shows the influence of microconstituents on
and
as in this case
is not a free field, hence
is no longer zero.
Case II
In this study we consider the influence of both
and
for fixed Reynolds number. We consider the following
Model A (CCM),
(same as in Case I)
Model C (NCCM with rotational inertial physics),
(same as in Case I)
Discretization, p-levels and other details remain the same as in case I. Figures 11(a)-(d) and Figures 12(a)-(d) show plots of velocity
versus
for
,
; and
,
. As we mentioned
Figure 10. Developing flow between parallel plates: Model B:
versus
and
versus
.
earlier both
and
result in resistance to the flow but we expect that the degree of resistance is not the same for
and
with similar change in their value. When comparing Figures 11(a)-(d) with Figures 9(a)-(d) we note that for same values of
as in Figures 9(a)-(d) but with
, we observe further reduction in the flow rate when
. For
(Figures 12(a)-(d)) flow rate is further reduced compared to
.
Case III
In this study we consider fixed
, choose a fixed value of
and vary
to study evolution of velocity
. We choose
and
and 0.5. Remaining details are same as in case I. Figures 13(a)-(d) show plots of
versus
for
,
,
and
. At
we almost have stationary state of the flow. As expected, larger values of
offer more resistance to flow and requires more time to reach stationary state and vice versa. We note that stationary state is same for all
values. Figure 14(a) and Figure 14(b) show evolution of
versus
and
versus
(case II) for
,
and
. Comparing this with Figure 10(c) and Figure 10(d) (
,
), we clearly note reduced values of
during the evolution for
in Figure 14(a),
Figure 11. Developing flow between parallel plates: velocity
versus position
: Models A and C (Case II),
.
Cauchy moment
adjusts accordingly. This of course is due to increased resistance to flow due to
.
Developing Couette flow
The configuration of parallel plates and other details remain the same as in Figure 8. We choose
in all studies. A 10 element uniform mesh with p-levels of nine in space and time and a solution of class C00 with
is considered. At the top plate (
) a velocity of 1.0 is applied over
in a continuous and differentiable manner i.e.,
and
.
Case A
Model A (CCM):
Model B (NCCM):
;
Figures 15(a)-(c) show evolution of
versus
at
,
and
. For each value of time the velocity
for both
values is lower than from Model A (CCM) as expected. Furthermore, higher values of
yield lower values of velocity
compared to lower values of
during the entire
Figure 12. Developing flow between parallel plates: velocity
versus position
: Models A and C (Case II),
.
evolution, demonstrating increasing resistance to flow with increasing value of
during the entire evolution.
Case B
In case B we present two studies. In the first study, we choose
Model A (CCM):
Model C (NCCM):
;
;
Computations are performed using the same mesh and other details as in case A. Figures 16(a)-(d) show plots of evolution of
versus
at
and
for
and
for both values of
in each case. Once again, combination of higher values of
and
result in more reduction in the velocity due to higher resistance to flow. Velocities for
are lower than those for
for both values of
. In the second study, we choose a fixed
and vary
. Figure 17(a) and Figure 17(b) shows evolution of
versus
at
and
. Progressively increasing values of
results in progressively reduced values of velocity
. Figures 18(a)-(c) show evolution of
versus
for: Model A, CCM (
); Model B, NCCM (
,
) and Model C, NCCM (
,
and
Figure 13. Developing flow between parallel plates: velocity
versus position
: Model C (Case III).
). Corresponding
versus
graphs are shown in Figures 19(a)-(c). In Figure 18(a) and Figure 18(b) we clearly observe the change in
free field due to micropolar physics. In Figure 18(c), we observe further reduction and change in
due to
. Evolution of
versus
follows a trend opposite to
i.e.,
introduces nonzero
and introduction of
further changes evolution of existing nonzero
due to only
.
3.2.2. Remarks
1) In this work,
and
are two parameters related to the micropolar physics. We have seen in the model problem studies that both offer increasing resistance to flow with their increasing values.
2) Model problem studies show that resistance to flow due to
(
and
combined) is more pronounced compared to
.
3) Parameters
and
are properties of the micropolar medium controlled by the microconstituents, but the manner in which they exert their influence on flow physics is different. The parameter
is a measure of the
Figure 14. Developing flow between parallel plates: Model B (Case II):
and
versus
:
,
.
collective resistance offered to the flow by each microconstituents. Whereas
appears as
, suggesting that it is a volumetric or mass effect through angular acceleration. Thus
is the collective influence of a group of microconstituents in a unit volume and perhaps provides a better explanation of why
is more influential than
.
4. Summary and Conclusions
The micropolar NCCT for solids derived using internal rotations (due to deformation gradient tensor) with rotational inertial physics [3] and the micropolar
Figure 15. Developing Couette flow: velocity
versus position
: Models A and B (Case A).
Figure 16. Developing Couette flow: velocity
versus position
: Models A and C (Case B).
NCCT for fluids derived using internal rotation rates (due to velocity gradient tensor) with rotational inertial physics [4] are considered in the work presented in this paper, to present various model problem studies that illustrate the influence of rotational inertial physics in case of micropolar solids as well as micropolar fluids. In this section, we present a summary of the work and draw some conclusions from it.
1) In reference [3] , by examining the BLM and BAM for micropolar solids with NCCT based on internal rotations, authors established co-existence of translational and rotational waves in micropolar solids with small deformation, small strain physics.
2) In reference [4] , authors considered BLM and BAM of the micropolar NCCT for fluids to show that in this case neither translational nor rotational waves can exist due to absence of elasticity or stiffness due to absence of strain physics in both CCM as well as NCCM.
3) Conclusions (1) and (2) in references [3] [4] were purely based on mathematical model. Solutions of IVPs resulting from the conservation and balance laws and the constitutive theories were not presented.
Figure 17. Developing Couette flow: velocity
versus position
: (Case B, second study).
4) In this paper, we consider the mathematical models of reference [3] [4] with rotational inertial physics to construct model problems (IVPs) and present their solutions using space-time coupled finite element method. The model problems are intentionally kept simple so that significant aspects of the micropolar physics due to microconstituents and in particular influence of rotational inertial physics can be demonstrated clearly.
Figure 18. Developing Couette flow:
versus
: Models A, B and C.
Figure 19. Developing Couette flow:
versus
: Models A, B and C. Developing Couette flow:
versus
: Models A, B and C.
5) The solutions of the model problems (IVPs) are obtained by using space-time coupled finite element method based on space-time residual functional for a space-time strip with time marching. The space-time local approximation are p-version hierarchical in space and time with higher order global differentiability in space and time. In this approach, with the choice of minimally confirming spaces, when the space-time residual functional for a space-time strip is O(10−8) or lower, PDEs in the IVPs are satisfied accurately, hence the computed evolutions are almost time accurate.
6) In the present work, we consider micropolar NCCT of references [3] [4] but with dissipation mechanisms due to: strain rate (CCM) and due to rate of symmetric part of internal rotational gradient tensor (Micropolar nonclassical dissipation). While the strain rate dissipation (CCT) is viscous, the nonviscous dissipation mechanisms in micropolar NCCT is due to microconstituents and the fluid medium. These aspects were not considered in the micropolar NCCT in references [3] [4] , hence is additional new physics in the mathematical models.
7) From the one dimensional numerical studies presented for translational and rotational waves in micropolar solid medium using Models TW1, RW1 and RW2 we observe:
a) Existence and propagation of translational wave (CCM) in the absence and in the presence of strain rate dissipation. The translational wave propagates, reflects, the reflected wave propagates and reflects from the free boundary without amplitude decay and base elongation when the medium is inviscid. The same phenomenon exists in the presence of damping but with continued amplitude decay and base elongation during evolution. Amplitude decay is most pronounced for the incident wave where as base elongation is more prominent in the reflected waves.
b) In the case of Model RW1, in which the rotational wave speed is one, the evolution of rotational wave physics is exactly same as in Model TW1.
c) Model RW2 permits the choice of wave speed. Studies for wave speed faster than one and slower than one with and without dissipation are presented. In 2D and 3D applications choice of
using
or
permits translational wave speed of one or rotational wave speed of one, but both wave speeds cannot be one. In this study, we show that wave speed different than one only influences when the propagating wave reaches the boundaries. The studies conclusively demonstrate that rotational inertial physics in micropolar solids is essential for the existence of rotational waves in micropolar solids.
d) In the studies presented here, rotational and translational waves are decoupled but coexist. In
and
this is not the case (Section 3.1.1). Since translational waves depend upon gradients of displacements in
and the rotational waves depend upon the gradients of displacements in
, both
and
being due to
, we expect coupling in the two in
and
.
e) Both micropolar and non-micropolar media are assumed isotropic and homogeneous, thus this micropolar theory can not simulate wave dispersion in micropolar media due to microconstituents. This theory can only simulate the influence of micropolar physics on wave propagation physics without dispersion.
8) In the case of micropolar fluid, neither translational nor rotational waves exist. The 1D form of BAM in this case is a time dependent diffusion equation. Two numerical studies are presented for two values of diffusion coefficient. From the studies, we clearly see higher values of the diffusion coefficient result in faster diffusion of the applied pulse.
9) Numerical studies presented for pressure driven developing flow between parallel plates show:
a) Micropolar physics with increasing
(
) results in progressively increasing resistance to flow and progressively decreasing flow rate. This holds during the entire evolution for each value of time.
b) We have shown that
, which controls rotational inertial physics, also offers resistance to flow. Increasing values of
results in decreasing flow rate.
c) Thus, in micropolar fluids with both
and
, the resulting flow rate is decreased. Since the constitutive model is not calibrated their relative influence on the flow rate is difficult to ascertain. Based on the values of
and
used,
influences flow rate more than
. We have discussed the physics related to
and
in Section 3.2.2. Both
and
offer resistance to flow but the physics of the resistance mechanisms is different in the two cases.
d) The significance of these studies is that we can conclusively see that rotational inertial physics in micropolar nonclassical continuum theories for fluids offers further resistance to the motion of the fluids over and beyond rotation rates or rotation rate gradients, thus reducing velocities and flow rates.
10) Numerical studies presented for Couette flow confirms the findings reported in item (9).
11) Presence of rotational inertial physics in fluids only offers added resistance to fluid motion. It can not possibly result in rotational waves as micropolar fluids have no elasticity associated with micropolar physics. In the case of micropolar fluids, BAM is not a wave equation in rotation rate.
12) Evolution of
versus
and
versus
(Figure 10(a) and Figure 10(b)) for Model A (CCM) and their comparison with similar evolution when
(
) i.e., model B shown in Figure 10(c) and Figure 10(d) and those for
and
in Figure 14(a) and Figure 14(b) clearly demonstrate how the free field
and zero
in Figure 10(a) and Figure 10(b) are affected by the presence of microconstituents without and with rotational inertial physics. These studies show that increasing
(
) reduces
field due to resistance offered by the microconstituents. When both
and
are nonzero, resistance offered by the microconstituents increases resulting in further reduction in
field. The consequences of reducing
field are increase in
as shown in the graphs.
Acknowledgements
First author is grateful for his endowed professorship and the department of mechanical engineering of the University of Kansas for providing financial support to the second author. The computational facilities provided by the Computational Mechanics Laboratory of the mechanical engineering department are also acknowledged.
Nomenclatures
Diffusion coefficient in Model RW2
,
Material coefficient for constitutive theory for Cauchy moment tensor
Nonclassical physics material coefficient
Dimensionless nonclassical physics material coefficient
Rotational inertial physics nonclassical material coefficient
Dimensionless dissipation coefficient fluids
Damping coefficient Model TW1
Damping coefficient Model RW1, RW2
Modulus of elasticity
Reference modulus of elasticity
Specific internal energy in Lagrangian description
Specific internal energy in Eulerian description
,
,
Permutation tensor
Strain tensor in Lagrangian description
Strain tensor in Eulerian description
Strain rate tensor in Lagrangian description
Strain rate tensor in Eulerian description
Viscosity in Lagrangian description
Viscosity in Eulerian description
Heat vector in Lagrangian description
Heat vector in Eulerian description
Internal rotational inertia in Lagrangian description
Internal rotational inertia in Eulerian description
Deformation gradient tensor in Lagrangian description
Symmetric part of deformation gradient tensor in Lagrangian description
Skew Symmetric part of deformation gradient tensor in Lagrangian description
Displacement gradient tensor in Lagrangian description
Symmetric part of displacement gradient tensor in Lagrangian description
Skew symmetric part of displacement gradient tensor in Lagrangian description
,
Length, reference length
,
Velocity gradient tensor in Lagrangian description
Velocity gradient tensor in Eulerian description
Symmetric Cauchy moment tensor in Lagrangian description component 11
Symmetric Cauchy moment tensor in Eulerian description component 11
Symmetric Cauchy moment tensor in Eulerian description component 23
Skew symmetric Cauchy moment tensor in Eulerian description component 23
,
Internal rotational velocity or rotation rates
Internal rotation rate in Lagrangian description about the x1 axis
Reference rotational velocity
Internal rotation rate in Eulerian description about the x1 axis
Internal rotation rate in Eulerian description about the x3 axis
Space-time domain
Discretization of space-time domain
ith space-time strip
Discretization of ith space-time strip
Space-time element e
Reference Density
Density in Lagrangian description
Density in Eulerian description
Pressure in Eulerian description
Cauchy stress tensor in Lagrangian description
Symmetric Cauchy stress tensor in Lagrangian description
Skew symmetric Cauchy stress tensor in Lagrangian description
Symmetric equilibrium Cauchy stress tensor in Lagrangian description
Symmetric deviatoric Cauchy stress tensor in Lagrangian description
Symmetric deviatoric contravariant Cauchy stress in Eulerian description component 12
Antisymmetric deviatoric contravariant Cauchy stress in Eulerian description component 12
Time
Reference time
Reference stress
,
,
Internal or classical rotations in Lagrangian description
,
,
Internal or classical rotations in Eulerian description
Internal or classical rotation rates in Lagrangian description
Internal or classical rotation rates in Eulerian description
Internal rotation about the x1 axis
Internal rotation about the x3 axis
,
,
Velocity components in x1, x2 and x3 directions
Reference velocity
,
,
Velocity vector in Lagrangian description
,
,
Velocity vector in Eulerian descriptions
,
,
Cartesian Coordinates
,
,
Cartesian Coordinates
Cartesian Coordinates
Cartesian Coordinates
Abbreviations
BAM: Balance of angular momentum
BLM: Balance of linear momentum
BMM: Balance of moment of moments
BVP: Boundary value problem
CBL: Conservation and balance laws
CCM: Classical continuum mechanics
CCT: Classical continuum theory
CM: Conservation of mass
FLT: First law of thermodynamics
IVP: Initial value problem
NCCM: Non-classical continuum mechanics
NCCT: Non-classical continuum theory
ODE: Ordinary differential equation
PDE: Partial differential equation
SLT: Second law of thermodynamics