NCCT for Micropolar Solid and Fluid Media Based on Internal Rotations and Rotation Rates with Rotational Inertial Physics: Model Problem Studies

Abstract

This paper presents model problem studies for micropolar thermoviscoelastic solids without memory and micropolar thermoviscous fluid using micropolar non-classical continuum theories (NCCT) based on internal rotations and rotation rates in which rotational inertial physics is considered in the derivation of the conservation and balance laws (CBL). The dissipation mechanism is due to strain rates as well as rotation rates. Model problems are designed to demonstrate and illustrate various significant aspects of the micropolar NCCT with rotational inertial physics considered in this paper. In case of micropolar solids, the translational and rotational waves are shown to coexist. In the absence of microconstituents (classical continuum theory, CCT) the internal rotations are a free field, hence have no influence on CCT. Absence of gradients of displacements and strains in micropolar thermoviscous fluid medium prohibits existence of translational waves as well as rotational waves even though the appearance of the mathematical model is analogous to the solids, but in terms of strain rates. It is shown that in case of micropolar thermoviscous fluids the BAM behaves more like time dependent diffusion equation i.e., like heat conduction equation in Lagrangian description. The influence of rotational inertial physics is demonstrated using BLM as well as BAM in the model problem studies.

Share and Cite:

Surana, K.S. and Kendall, J.K. (2023) NCCT for Micropolar Solid and Fluid Media Based on Internal Rotations and Rotation Rates with Rotational Inertial Physics: Model Problem Studies. Applied Mathematics, 14, 612-651. doi: 10.4236/am.2023.149037.

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

ρ 0 D 2 u k D t 2 ρ 0 F k b σ j k x j = 0 (1)

I Θ 0 ρ 0 D 2 ( i Θ k ) D t 2 ϵ i j k σ i j m j k x j = 0 (2)

ϵ j k l ( I Θ 0 ρ 0 v j ( i ω k ) m j k ) = 0 (3)

ρ 0 D e D t + q i x i tr ( [ σ s ] [ ε ˙ ] ) tr ( [ m s ] [ J ˙ s i Θ ] ) tr ( [ m a ] [ J ˙ a i Θ ] ) = 0 (4)

Additivedecompositionof σ : σ = σ s + σ a ; σ s = σ s d + σ s e (5)

σ s d = σ 0 | Ω _ I + 2 μ ε + λ ( tr ( ε ) ) I + 2 η ε ˙ + κ tr ( ε ˙ ) I (6)

m s = 2 μ m ( J s i Θ ) + 2 η m ( J ˙ s i Θ ) (7)

q = κ g (8)

In the mathematical model we have a total of 25 equations BLM(3), BAM(3), BMM(3), FLT(1), constitutive theories for σ s ( 6 ) , m s ( 6 ) , q ( 3 ) in a total of 25 dependent variables u ( 3 ) , σ s d ( 6 ) , σ a ( 3 ) , m s ( 6 ) , m a ( 3 ) , q ( 3 ) , θ ( 1 ) hence the model has closure. The last term in the entropy inequality

q : g θ σ s d : ε ˙ m s : J ˙ s i Θ m a : J ˙ a i Θ 0 (9)

must be set to zero to ensure that (9) is always satisfied. Then, in addition to (1)-(8). We must also satisfy

m a : J ˙ a i Θ = 0 (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.

D ρ ¯ D t + ρ ¯ div ( v ¯ ) = 0 ( CM ) (CM) (11)

ρ ¯ v ¯ t + ρ ¯ ( v ¯ : ¯ ) v ¯ ρ ¯ F ¯ b ¯ : σ ¯ ( 0 ) = 0 (BLM) (12)

I ¯ Θ ρ ¯ D D t ( i ω ¯ k ) ϵ i j k σ ¯ i j ( 0 ) m ¯ m k , m ( 0 ) = 0 (BAM) (13)

ϵ j k l ( I ¯ Θ ρ ¯ v ¯ j ( i ω ¯ k ) m ¯ j k ( 0 ) ) = 0 (BMM) (14)

ρ ¯ D e ¯ D t + ¯ : q ¯ σ ¯ s ( 0 ) : D ¯ m ¯ s ( 0 ) : J ¯ s i r Θ m ¯ a ( 0 ) : J ¯ a i r Θ = 0 (FLT) (15)

ρ ¯ ( D Φ ¯ D t + η ¯ D θ ¯ D t ) + q ¯ : g ¯ θ ¯ σ ¯ s ( 0 ) : D ¯ m ¯ s ( 0 ) : J ¯ s i r Θ m ¯ a ( 0 ) : J ¯ a i r Θ 0 (SLT) (16)

σ ¯ s d ( 0 ) = s ¯ | Ω _ + 2 η D ¯ + κ ( tr ( D ¯ ) ) I (17)

Additive decomposition of σ ¯ ( 0 ) : σ ¯ ( 0 ) = σ ¯ s ( 0 ) + σ ¯ a ( 0 ) ; σ ¯ s ( 0 ) = σ ¯ s d ( 0 ) + σ ¯ s e ( 0 ) (18)

m ¯ s ( 0 ) = 2 η m ( J ¯ s i r Θ ) (19)

q ¯ = κ g ¯ (20)

where

ω ¯ i = Θ ¯ i r , D ( ω ¯ i ) D t = Θ ¯ ˙ i r (21)

In this mathematical model we have CM(1), BLM(3), BAM(3), BMM(3), FLT(1) and the constitutive theories for σ ¯ s d ( 0 ) ( 6 ) , m ¯ s ( 0 ) ( 6 ) , q ¯ ( 3 ) , a total of 26 equations in 26 dependent variables ρ ¯ ( 1 ) , v ¯ ( 3 ) , σ ¯ s d ( 0 ) ( 6 ) , σ ¯ a ( 0 ) ( 3 ) , m ¯ s ( 0 ) ( 6 ) , m ¯ a ( 0 ) ( 3 ) , q ¯ ( 3 ) , θ ¯ ( 1 ) . The last term in the entropy inequality must be set to zero to ensure that (16) is not violated. Thus

m ¯ s ( 0 ) : J ¯ a i r Θ = 0 (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 σ 11 = σ s 11 = σ s d 11 and using x for x1)

ρ 0 2 u 1 t 2 ( σ s d 11 ) x = 0 (23)

I Θ 0 ρ 0 2 ( i Θ 1 ) t 2 s m 11 x = 0 (24)

σ s d 11 = ( 2 μ + λ ) u 1 x + ( 2 η + κ ) t ( u 1 x ) (25)

m s 11 = ( α 1 ) ( i Θ 1 ) x + ( α 2 ) t ( ( i Θ 1 ) x ) α 1 = μ m , α 2 = η m (26)

Substituting (25) in (23) and (26) in (24) we obtain

ρ 0 2 u 1 t 2 ( 2 μ + λ ) 2 u 1 x 2 ( 2 η + κ ) t ( 2 u 1 x 2 ) = 0 (27)

I Θ 0 ρ 0 2 ( i Θ 1 ) t 2 ( α 1 ) 2 ( i Θ 1 ) x 2 ( α 2 ) t ( 2 ( i Θ 1 ) x 2 ) = 0 (28)

for 1D case and for incompressible physics 2 μ + λ = E , modulus of elasticity, as Poison’s ratio ν = 0 and κ = 0 , hence (28) reduces to (but (28) remains same).

ρ 0 2 u 1 t 2 E 2 u 1 x 2 2 η t ( 2 u 1 x 2 ) = 0 (29)

I Θ 0 ρ 0 2 ( i Θ 1 ) t 2 ( α 1 ) 2 ( i Θ 1 ) x 2 ( α 2 ) t ( 2 ( i Θ 1 ) x 2 ) = 0 (30)

Remarks

1) Equation (29) describes translational or σ s d 11 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 u 1 , ρ 0 , E, and 2 η in (29) by i Θ 1 , I Θ 0 ρ 0 , α 1 and α 2 . Thus, (30) describes rotational waves i.e., moment ( m s 11 ) wave. The third term in (30) is due to micropolar dissipation physics which would result in amplitude decay and base elongation of m s 11 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)

ρ ^ 0 2 u ^ 1 t ^ 2 E ^ 2 u ^ 1 x ^ 2 2 η ^ t ^ ( 2 u ^ 1 x ^ 2 ) = 0 (31)

I Θ 0 ρ ^ 0 2 ( i Θ ^ 1 ) t ^ 2 ( α ^ 1 ) 2 ( i Θ ^ 1 ) x ^ 2 ( α ^ 2 ) t ^ ( 2 ( i Θ ^ 1 ) x ^ 2 ) = 0 (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 L 0 , η 0 , ( ρ 0 ) r e f , t 0 , v 0 = L 0 / t 0 = E 0 ρ 0 and E 0 = ( ρ 0 ) r e f v 0 2 as reference quantities in which v 0 is the speed of translational wave (speed of sound in CCM).

ρ 0 2 u 1 t 2 E E 0 ( ρ 0 ) r e f v 0 2 2 u 1 x 2 2 η R e t ( 2 u 1 x 2 ) = 0 (33)

where 1 / R e = η 0 L 0 ( ρ 0 ) r e f v 0 , Reynolds number. If we choose ( ρ 0 ) r e f = ρ ^ 0 , E 0 = E ^ = ( ρ 0 ) r e f v 0 2 and C d = 2 η R e as damping coefficient, then (33) can be reduced to

2 u 1 t 2 2 u 1 x 2 c d t ( 2 u 1 x 2 ) = 0 (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 ( I Θ 0 ) r e f , ( α 1 ) 0 , ( α 2 ) 0 , ω 0 = ( α 1 ) 0 / ( I Θ 0 ) r e f ( ρ 0 ) r e f and t 0 = L 0 / ω 0 , ( α 1 ) 0 = ( I Θ 0 ) r e f ( ρ 0 ) r e f ω 0 2 , then (32) can be written as

I Θ 0 ρ 0 2 ( i Θ 1 ) t 2 ( α 1 ) ( α 1 ) 0 ( I Θ 0 ) r e f ( ρ 0 ) r e f ω 0 2 2 ( i Θ 1 ) x 2 ( α 2 ) 1 R e * t ( 2 ( i Θ 1 ) x 2 ) = 0 (35)

If ( I Θ 0 ) r e f = I ^ Θ 0 , ( ρ 0 ) r e f = ρ ^ 0 , then (35) can be reduced to

2 ( i Θ 1 ) t 2 2 ( i Θ 1 ) x 2 ( α 2 ) R e * t ( 2 ( i Θ 1 ) x 2 ) = 0 (36)

in which R e * = ( I Θ 0 ) r e f ( ρ 0 ) r e f ω 0 L 0 ( α 2 ) 0 is the Reynolds number associated with dissipation due to micropolar physics. If we define c d * = α 2 R e * as the dissipation coefficient due to micropolar physics, then (36) can be written as

2 ( i Θ 1 ) t 2 2 ( i Θ 1 ) x 2 c d * t ( 2 ( i Θ 1 ) x 2 ) = 0 (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

v 0 = E 0 ( ρ 0 ) r e f t 0 = L 0 v 0 } (38)

using these reference values we can nondimensionalize (32) as follows.

I Θ 0 ρ 0 2 ( i Θ 1 ) t 2 ( α 1 ) ( α 1 ) 0 ( I Θ 0 ) r e f ( ρ 0 ) r e f v 0 2 2 ( i Θ 1 ) x 2 ( α 2 ) ( α 2 ) 0 ( I Θ 0 ) r e f ( ρ 0 ) r e f L 0 v 0 t ( 2 ( i Θ 1 ) x 2 ) = 0 (39)

If we choose ( I Θ 0 ) r e f = I ^ Θ 0 , ( ρ 0 ) r e f = ρ ^ 0 and define

α 1 * = ( α 1 ) 0 / 2 ( I Θ 0 ) r e f ( ρ 0 ) r e f v 0 2 , c d * = α 2 R e * , 1 / R e * = ( α 2 ) 0 ( I Θ 0 ) r e f ( ρ 0 ) r e f L 0 v 0

then (39) can be reduced to

2 ( i Θ 1 ) t 2 α 1 * 2 ( i Θ 1 ) x 2 c d * t ( 2 ( i Θ 1 ) x 2 ) = 0 (40)

In this case the speed of dimensionless rotational wave is α 1 * 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

v 1 t ( σ s d 11 ) x = 0 σ s d 11 u 1 x c d v 1 x = 0 v 1 u 1 t = 0 } (41)

Rotational (or moment) waves (Model RW1)

We recast (37) as

( i ω 1 ) t ( s m 11 ) x = 0 s m 11 ( i Θ 1 ) x c d * ( i ω 1 ) x = 0 i ω 1 ( i Θ 1 ) t = 0 } (42)

Rotational (or moment) waves (Model RW2)

( i ω 1 ) t ( s m 11 ) x s m 11 α 1 * ( i Θ 1 ) x c d * ( i ω 1 ) x = 0 i ω 1 ( i Θ 1 ) t = 0 } (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 Δ t = 0.1 is used in all computations. The solution is computed for the first space-time strip ( 0 t Δ t ) 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 k = 2 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 Ω ¯ x t , discretization Ω ¯ x t T = e Ω ¯ x t ( e ) of the space time domain into space-time strips, discretization ( Ω ¯ x t ( 1 ) ) T = e Ω ¯ x t ( e ) of the first space-time strip Ω ¯ x t ( 1 ) and the details of stress σ s d 11 boundary conditions (for TW1) applied to the last space-time element face at x = 1 . In the numerical studies, we choose negative σ 0 (compressive) for the applied pulse at x = 1 . The details of applying BC at x = L 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, c d = 0 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 Δ t = 0.1 , the applied pulse at x = 1 reaches the impermeable boundary at x = 0 . Due to zero damping, we expect the applied pulse shape to be preserved during evolution. Figure 2(a) shows a plot of σ s d 11 versus x for various values of time. Time steps 2, 6, 10 show the incident compressive pulse entering the spatial domain at x = L 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 x = 0 as expected. Reflection of the compressive incident pulse at the impermeable boundary (at x = 0 ) 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 x = L = 1 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: σ s d 11 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 x = 0 (time steps 23 and 28 in Figure 2(a)).

In the next study for TW1, we consider dimensionless damping coefficient c d = 0.002 . Figure 2(b) shows evolution of σ s d 11 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 k = 2 (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 ( Ω ¯ x t ( i ) ) T , discretization of space-time strip Ω ¯ x t ( i ) 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 Ω ¯ x t ( i ) , 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., m s 11 pulse of time duration 2 Δ t (similar to σ s d 11 shown in Figure 1(e)). We choose p-levels of 9 in space and time, Δ t = 0.1 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 m s 11 = 1 . 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 c d * = 0 (undamped case). Incident moment pulse propagates towards the impermeable boundary without amplitude decay and base elongation. Upon reflection at x = 0 , 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 x = 0 .

Figure 3(b) shows the propagation, reflection and propagation upon reflection of the same compressive m s 11 pulse as in Figure 3(a) but in the presence of dissipation. We choose dimensionless dissipation coefficient c d * = 0.002 . From Figure 3(b) we observe almost the same behavior of the moment pulse as that of σ s d 11 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 α 1 * and the dissipation is controlled by the dimensionless dissipation coefficient c d * . We apply a moment ( m s 11 ) pulse of duration 2 Δ t at the boundary x = L = 1 .

In the first study, we choose Δ t = 0.1 , α 1 * = 2.25 and c d * = 0 (undamped case). The wave speed is 2.25 = 1.5 , thus the pulse would reach x = 0

Figure 3. Evolution of Cauchy Moment wave: s m 11 versus position x (Model RW1).

boundary in ( 1 / 1.5 ) / 0.1 = 6.67 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 α 1 * = 2.25 i.e., at the wave speed of 1.5 when c d * = 0.002 . Diminished amplitudes and elongated base of the pulse is clearly observed during evolution. This is similar to RW1 ( c d * = 0.002 ). In the next study, we apply the same moment pulse as in the previous study but choose α 1 * = 0.49 and c d * = 0 , hence rotational wave speed of 0.49 = 0.7 , thus the moment pulse will reach the boundary at x = 0 in ( 1 / 0.7 ) / 0.1 = 14.285 time steps. Figure 5(a) and Figure 5(b) show propagation and reflection of the moment wave for c d * = 0 and c d * = 0.002 . The undamped pulse ( c d * = 0 ) maintains its base and amplitude during the entire evolution and when c d * = 0.002 progressive amplitude decay and base elongation of the pulse is observed for both values of α 1 * .

Figure 4. Evolution of Cauchy Moment wave: s m 11 versus position x (Model RW2: α * = 2.25 ).

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 Δ t . In this case, we see perfect reflection of the m s 11 pulse at the 11th time step and the peak value of m s 11 doubles, same as in the case of the translational stress wave. However, when rotational wave speed is not an integer multiple of Δ t , precise arrival of m s 11 pulse at x = 0 cannot be simulated as elapsed time is always an integer multiple of Δ t . For this reason, we are not able to observe 2 ( m s 11 ) 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: s m 11 versus position x (Model RW2: α * = 0.49 ).

Remarks

1) In the model problem studies presented here the translational and rotational waves are decoupled. In R 2 and R 3 this may not be the case. Since the stress tensor is a function of [ J s d ] and the moment tensor is a function of [ J a d ] , both stress and moment tensor waves are dependent on gradients of displacement. Thus, in R 2 and R 3 we expect interaction between the two waves. That is, as shown here, rotational waves depend upon α ^ 1 and I ¯ Θ which in turn influences [ J d ] , thus [ J s d ] , hence influencing translational waves. The model problem studies in R 2 and R 3 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 R 2 and R 3 (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 x ¯ for x ¯ 1 ).

I ¯ Θ ρ ¯ ( i ω ¯ 1 ) t ( s m ¯ 11 ) x ¯ = 0 (44)

s m ¯ 11 = α ¯ 1 ( 1 2 ( i ω ¯ 1 ) x ¯ ) (45)

Substituting (45) in (44)

I ¯ Θ ρ ¯ ( i ω ¯ 1 ) t α ¯ 1 2 ( 2 ( i ω ¯ 1 ) x ¯ 2 ) = 0 (46)

Remarks

1) We clearly see that this Equation (46) will not permit a moment ( s m ¯ 11 )wave as it is not a wave equation due to the absence of 2 ( i ω ¯ 1 ) t 2 in place of ( i ω ¯ 1 ) t .

2) However it is interesting to study its solution for varying I ¯ Θ values.

Dimensionless form of the mathematical model

Recasting (46) with hat ( ) on all quantities indicating that they have their usual dimensions (or units)

I ¯ ^ Θ ρ ¯ ^ ( i ω ¯ ^ 1 ) t ^ α ¯ ^ 1 2 ( 2 ( i ω ¯ ^ 1 ) x ¯ ^ 2 ) = 0 (47)

and using

I ¯ Θ = I ¯ ^ Θ / I Θ 0 , ρ ¯ = ρ ¯ ^ / ρ 0 , x ¯ = x ¯ ^ / L 0 ω ¯ i 1 = ω ¯ ^ i 1 / ω 0 , t 0 = L 0 / ω 0 , and ω 0 = ( α 1 ) 0 / I Θ 0 ρ 0 , α ¯ 1 = α ¯ ^ 1 / ( α 1 ) 0 } (48)

We can write (47) as

I ¯ Θ ρ ¯ ( i ω ¯ 1 ) t ( α ¯ 1 2 ) ( α 1 ) 0 L 0 I Θ 0 ρ 0 ( 2 ( i ω ¯ 1 ) x ¯ 2 ) = 0 (49)

If we choose I Θ 0 = I ¯ ^ Θ , ρ 0 = ρ ¯ ^ and ( α 1 ) 0 = α ¯ ^ 1 then (49) can be reduced to

( i ω ¯ 1 ) t c ¯ 2 ( 2 ( i ω ¯ 1 ) x ¯ 2 ) = 0 (50)

in which c ¯ 2 = ( α ¯ ^ 1 ) / ( 2 L 0 I ¯ ^ Θ ρ ¯ ^ ) , dimensionless dissipation coefficient. PDE (50) can be cast as a system of two first order PDEs by using moment s m ¯ 11 as auxiliary variable and we have

( i ω ¯ 1 ) t ( s m ¯ 11 ) x ¯ = 0 s m ¯ 11 = c ¯ 2 ( i ω ¯ 1 ) x ¯ } (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 c ¯ 2 1 I ¯ Θ is the diffusion coefficient. Thus low values of I ¯ Θ correspond to high values of c ¯ 2 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 s m ¯ 11 moment pulse of duration 2 Δ t on the boundary at x ¯ = 1 such that s m ¯ 11 [ 0, 1 ] t [ 0, Δ t ] ; s m ¯ 11 [ 1,0 ] t [ Δ t ,2 Δ t ] and s m ¯ 11 = 0 for t 2 Δ t (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 Δ t = 0.1 with local approximation of class C11 in space and time, We choose two values of c ¯ 2 = 0.0001 , 0.001 . Evolutions for different values of time for the two choices of c ¯ 2 are shown in Figure 6(a) and Figure 6(b). For both values of c ¯ 2 , the applied pulse progressively diffuses as time elapses. We observe significantly high diffusion of the applied moment pulse for c ¯ 2 = 0.001 compared to c ¯ 2 = 0.0001 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 s m ¯ 11 versus position x ¯ .

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 y ¯ direction) velocity u ¯ is a dependent variable, velocity v ¯ in the y ¯ 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.

x ¯ = x ¯ ^ / L 0 , y ¯ = y ¯ ^ / L 0 , η ¯ = η ¯ ^ / η 0 , α ¯ = α ¯ ^ / α 0 m ¯ 23 ( 0 ) = m ¯ ^ 23 ( 0 ) / m 0 , σ ¯ s d 21 ( 0 ) = σ ¯ ^ s d 21 ( 0 ) / τ 0 , σ ¯ a 21 ( 0 ) = σ ¯ ^ a 21 ( 0 ) / τ 0 m 0 = τ 0 L 0 , τ 0 = ρ 0 v 0 2 = p 0 ( characteristickineticenergy ) u ¯ = u ¯ ^ / v 0 , v ¯ = v ¯ ^ / v 0 , t 0 = L 0 / v 0 ω ¯ i 3 = ω ¯ ^ i 3 / ω 0 , ω 0 = t 0 , I ¯ Θ = I ¯ ^ Θ / I Θ 0 } (52)

Mathematical model based on CCM: Model A

Balance of linear momenta in the x ¯ direction and the constitutive theory for deviatoric Cauchy stress σ ¯ d 12 ( 0 ) = σ ¯ d 21 ( 0 ) 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.

ρ ¯ ^ u ¯ ^ t ^ + p ¯ ^ x ¯ ^ ( σ ¯ ^ s d 21 ( 0 ) ) y ¯ ^ = 0 σ ¯ ^ s d 21 ( 0 ) = η ¯ ^ u ¯ ^ y ¯ ^ } (53)

We can nondimensionalize (53) using (52) and if we choose ρ 0 = ρ ¯ ^ , η 0 = η ¯ ^ , then the dimensionless form of (53) can be written as

u ¯ t + p ¯ x ¯ ( σ ¯ s d 21 ( 0 ) ) y ¯ = 0 σ ¯ s d 21 ( 0 ) = 1 R e u ¯ y ¯ } (54)

In case of pressure driven flow, p ¯ x ¯ is given. For non-pressure driven Couette flow, p ¯ x ¯ = 0 . In (54), Re is Reynolds number and is given by (based on reference quantities) R e = L 0 ρ 0 v 0 η 0 . Equation (54) is a first order system of two PDEs in two dependent variables u ¯ and σ ¯ s d 21 ( 0 ) . 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 x ¯ direction, balance of angular momenta about the z ¯ direction, constitutive theory for symmetric part of deviatoric Cauchy shear stress σ ¯ s d 12 ( 0 ) = σ ¯ s d 21 ( 0 ) , constitutive theory for symmetric Cauchy moment s m ¯ 23 and definition of angular rotation rate i ω ¯ 3 i.e., rotation rate about z ¯ . We can write the following (with their usual dimensions) using the CBL and the constitutive theories given for 3 by Equations (11)-(22), in the absence of body forces and body moments.

ρ ¯ ^ u ¯ ^ t ^ + p ¯ ^ x ¯ ( σ ¯ ^ s d 21 ( 0 ) ) y ¯ ^ ( σ ¯ ^ a d 21 ( 0 ) ) y ¯ ^ = 0 ( m ¯ ^ s 23 ( 0 ) ) y ¯ ^ + 2 ( σ ¯ ^ a d 21 ( 0 ) ) = 0 σ ¯ ^ s d 12 ( 0 ) = σ ¯ ^ s d 21 ( 0 ) = η ¯ ^ u ¯ ^ y ¯ ^ m ¯ ^ s 23 ( 0 ) = α ¯ ^ 2 ( i ω ¯ ^ 3 ) y ¯ ω ¯ ^ i 3 = 1 2 u ¯ ^ y ¯ ^ } (55)

We can nondimensionalize (55) using (52), and if we choose I Θ 0 = I ¯ ^ Θ , ρ 0 = ρ ¯ ^ , η 0 = η ¯ ^ , then (55) can be written as

ρ ¯ u ¯ t + p ¯ x ¯ ( σ ¯ s d 21 ( 0 ) ) y ¯ ( σ ¯ a d 21 ( 0 ) ) y ¯ = 0 ( m ¯ s 23 ( 0 ) ) y ¯ + 2 ( σ ¯ a d 21 ( 0 ) ) = 0 σ ¯ s d 12 ( 0 ) = σ ¯ s d 21 ( 0 ) = 1 R e u ¯ y ¯ m ¯ s 23 ( 0 ) = α ¯ 2 ( i ω ¯ 3 ) y ¯ ω ¯ i 3 = 1 2 u ¯ y ¯ } (56)

in which α ¯ = α ¯ ^ / τ 0 L 0 2 t 0 , dimensionless material coefficient for micropolar nonclassical physics. Equations (56) is a system of five first order PDEs in five dependent variables: u ¯ , σ ¯ s d 21 ( 0 ) , σ ¯ a d 21 ( 0 ) , s m ¯ 23 ( 0 ) and i ω ¯ 3 .

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 x ¯ direction, BAM about z ¯ direction, BMM balance law, constitutive theory for symmetric part of the deviatoric Cauchy shear stress σ ¯ s d 12 ( 0 ) = σ ¯ s d 21 ( 0 ) , constitutive theory for symmetric Cauchy moment s m ¯ 23 ( 0 ) , definition of angular rotation rate i ω ¯ 3 and the constraint equations resulting from the entropy inequality.

Using the mathematical model in 3 (Equations (11)-(22)) and noting that in this case v ¯ = 0 , w ¯ = 0 , we can write the following for CBL and constitutive theories (in the absence of body forces and body moments) and the constraint equation

ρ ¯ ^ u ¯ ^ t ^ + p ¯ ^ x ¯ ( σ ¯ ^ s d 21 ( 0 ) ) y ¯ ^ ( σ ¯ ^ a d 21 ( 0 ) ) y ¯ ^ = 0 ρ ¯ ^ ( I ¯ ^ Θ ) ( ( ω ¯ ^ i 3 ) t ^ ) ( m ¯ ^ s 23 ( 0 ) ) y ¯ ^ ( m ¯ ^ a 23 ( 0 ) ) y ¯ ^ + 2 ( σ ¯ ^ a d 21 ( 0 ) ) = 0 σ ¯ ^ s d 12 ( 0 ) = σ ¯ ^ s d 21 ( 0 ) = η ¯ ^ u ¯ ^ y ¯ ^ m ¯ ^ s 23 ( 0 ) = α ¯ ^ 2 ( i ω ¯ ^ 3 ) y ¯ ω ¯ ^ i 3 = 1 2 u ¯ ^ y ¯ ^ } (57)

I ¯ ^ Θ v ¯ ^ ( ω ¯ ^ i 3 ) m ¯ ^ s 23 ( 0 ) m ¯ ^ a 23 ( 0 ) = 0 (BMM when v ¯ 0 ) (58)

s m ¯ ^ 23 ( 0 ) ( i ω ¯ ^ 3 ) = 0 (constraint equation) (59)

When v ¯ = 0 (which is the case here), BMM (58) reduces to

m ¯ ^ s 23 + m ¯ ^ a 23 = 0 (60)

Equation (60) suggests that m ¯ ^ a 23 ( 0 ) = m ¯ ^ s 23 ( 0 ) or the nonsymmetric moment tensor m ¯ ^ 23 ( 0 ) 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 s m ¯ ^ 23 ( 0 ) = 0 or i ω ¯ ^ 3 = 0 for all x ¯ , t Ω ¯ x t . 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 I Θ 0 = I ¯ ^ Θ , ρ 0 = ρ ¯ ^ , η 0 = η ¯ ^ , then we can obtain the following from (57).

ρ ¯ u ¯ t + p ¯ x ¯ ( σ ¯ s d 21 ( 0 ) ) y ¯ ( σ ¯ a d 21 ( 0 ) ) y ¯ = 0 β ¯ ( ω ¯ i 3 ) t ( m ¯ s 23 ( 0 ) ) y ¯ + 2 ( σ ¯ a d 21 ( 0 ) ) = 0 σ ¯ s d 12 ( 0 ) = σ ¯ s d 21 ( 0 ) = 1 R e u ¯ y ¯ m ¯ s 23 ( 0 ) = α ¯ 2 ( i ω ¯ 3 ) y ¯ ω ¯ i 3 = 1 2 u ¯ y ¯ } (61)

These are a system of five first order PDEs in five dependent variables: u ¯ , σ ¯ s d 21 ( 0 ) , σ ¯ a d 21 ( 0 ) , s m ¯ 23 ( 0 ) and i ω ¯ 3 . 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 α ¯ = α ¯ ^ / τ 0 L 0 2 t 0 and β ¯ = I Θ 0 / L 0 2 .

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, α ¯ 1 and β ¯ 1 all three control the physics. β ¯ 1 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, R e = 100 with Δ t = 1.0 . Flow is pressure driven with p ¯ x ¯ = 0.1 .

Case I

In the first study we choose

Model A (CCM), R e = 100

Model B (NCCM), R e = 100 ; α ¯ = 0.0001 , 0.001 .

Figures 9(a)-(d) show plots of the evolution of velocity u ¯ for t = 5 Δ t , 10 Δ t , 20 Δ t and 40 Δ t 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 t = 40 Δ t , 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 i ω ¯ 3 versus y ¯ for Model A (CCM). At t = 40 Δ t , we almost have stationary sate at which ( i ω ¯ 3 ) y ¯ = c o n s t a n t . Since i ω ¯ 3 is a free field, s m ¯ 23 = 0 ,

Figure 9. Developing flow between parallel plates: velocity u ¯ versus position y ¯ : Models A and B (Case I).

shown in Figure 10(b). Figure 10(c) and Figure 10(d) show evolution of i ω ¯ 3 versus y ¯ and s m ¯ 23 versus y ¯ for Model B for α ¯ = 0.001 . Comparison of these evolutions with those in Figure 10(a) and Figure 10(b) clearly shows the influence of microconstituents on i ω ¯ 3 and s m ¯ 23 as in this case i ω ¯ 3 is not a free field, hence s m ¯ 23 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), R e = 100 (same as in Case I)

Model C (NCCM with rotational inertial physics),

R e = 100 (same as in Case I)

α ¯ = 0.0001 α ¯ = 0.001 } β ¯ = 0.05 , 0.9

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 u ¯ versus y ¯ for α ¯ = 0.0001 , 0.001 , β ¯ = 0.05 ; and α ¯ = 0.0001 , 0.001 , β ¯ = 0.9 . As we mentioned

Figure 10. Developing flow between parallel plates: Model B: i ω ¯ 3 versus y ¯ and s m ¯ 23 versus y ¯ .

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 β ¯ = 0 , we observe further reduction in the flow rate when β ¯ = 0.05 . For β ¯ = 0.9 (Figures 12(a)-(d)) flow rate is further reduced compared to β ¯ = 0.05 .

Case III

In this study we consider fixed R e = 100 , choose a fixed value of α = 0.01 and vary β ¯ to study evolution of velocity u ¯ . We choose Δ t = 1.0 and β ¯ = 0.05 , 0.2 and 0.5. Remaining details are same as in case I. Figures 13(a)-(d) show plots of u ¯ versus y ¯ for t = 5 Δ t , 10 Δ t , 20 Δ t and 40 Δ t . At t = 40 Δ t 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 i ω ¯ 3 versus y ¯ and s m ¯ 23 versus y ¯ (case II) for R e = 100 , α ¯ = 0.001 and β ¯ = 0.9 . Comparing this with Figure 10(c) and Figure 10(d) ( α ¯ = 0.001 , β ¯ = 0.0 ), we clearly note reduced values of i ω ¯ 3 during the evolution for β ¯ = 0.9 in Figure 14(a),

Figure 11. Developing flow between parallel plates: velocity u ¯ versus position y ¯ : Models A and C (Case II), β ¯ = 0.05 .

Cauchy moment s m ¯ 23 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 R e = 5 in all studies. A 10 element uniform mesh with p-levels of nine in space and time and a solution of class C00 with Δ t = 0.01 is considered. At the top plate ( y ¯ = 1.0 ) a velocity of 1.0 is applied over Δ t in a continuous and differentiable manner i.e., u ¯ [ 0,1 ] t [ 0, Δ t ] and u = 1.0 t Δ t .

Case A

Model A (CCM): R e = 5

Model B (NCCM): R e = 5 ; α ¯ = 0.0001 , 0.001

Figures 15(a)-(c) show evolution of u ¯ versus y ¯ at t = 2 Δ t , 5 Δ t and 10 Δ t . For each value of time the velocity u ¯ for both α ¯ values is lower than from Model A (CCM) as expected. Furthermore, higher values of α ¯ yield lower values of velocity u ¯ compared to lower values of α ¯ during the entire

Figure 12. Developing flow between parallel plates: velocity u ¯ versus position y ¯ : Models A and C (Case II), β ¯ = 0.9 .

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): R e = 5

Model C (NCCM): R e = 5 ; α ¯ = 0.01 , 0.1 ; β ¯ = 0.05 , 0.2

Computations are performed using the same mesh and other details as in case A. Figures 16(a)-(d) show plots of evolution of u ¯ versus y ¯ at t = 5 Δ t and t = 20 Δ t for β ¯ = 0.05 and β ¯ = 0.2 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 β ¯ = 0.2 are lower than those for β ¯ = 0.05 for both values of α ¯ . In the second study, we choose a fixed α ¯ = 0.1 and vary β ¯ = 0.1 , 0.2 , 0.3 . Figure 17(a) and Figure 17(b) shows evolution of u ¯ versus y ¯ at t = Δ t and t = 5 Δ t . Progressively increasing values of β ¯ results in progressively reduced values of velocity u ¯ . Figures 18(a)-(c) show evolution of i ω ¯ 3 versus y ¯ for: Model A, CCM ( R e = 5 ); Model B, NCCM ( R e = 5 , α ¯ = 0.001 ) and Model C, NCCM ( R e = 5 , α ¯ = 0.001 and

Figure 13. Developing flow between parallel plates: velocity u ¯ versus position y ¯ : Model C (Case III).

β ¯ = 0.05 ). Corresponding s m ¯ 23 versus y ¯ graphs are shown in Figures 19(a)-(c). In Figure 18(a) and Figure 18(b) we clearly observe the change in i ω ¯ 3 free field due to micropolar physics. In Figure 18(c), we observe further reduction and change in i ω ¯ 3 due to β ¯ = 0.05 . Evolution of s m ¯ 23 versus y ¯ follows a trend opposite to i ω ¯ 3 i.e., α ¯ 0 introduces nonzero s m ¯ 23 and introduction of β ¯ further changes evolution of existing nonzero s m ¯ 23 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 β ¯ ( I ¯ Θ and ρ ¯ combined) is more pronounced compared to α ¯ .

3) Parameters α ¯ and I ¯ Θ 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): i ω ¯ 3 and s m ¯ 23 versus y ¯ : α ¯ = 0.001 , β ¯ = 0.9 .

collective resistance offered to the flow by each microconstituents. Whereas I ¯ Θ appears as I ¯ Θ ρ ¯ , suggesting that it is a volumetric or mass effect through angular acceleration. Thus I ¯ Θ ρ ¯ 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 u ¯ versus position y ¯ : Models A and B (Case A).

Figure 16. Developing Couette flow: velocity u ¯ versus position y ¯ : 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 u ¯ versus position y ¯ : (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: i ω ¯ 3 versus y ¯ : Models A, B and C.

Figure 19. Developing Couette flow: i ω ¯ 3 versus y ¯ : Models A, B and C. Developing Couette flow: s m ¯ 23 versus y ¯ : 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(108) 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 t 0 using v 0 / L 0 or i ω 0 / L 0 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 R 2 and R 3 this is not the case (Section 3.1.1). Since translational waves depend upon gradients of displacements in [ s d J ] and the rotational waves depend upon the gradients of displacements in [ a d J ] , both [ s d J ] and [ a d J ] being due to [ J ] , we expect coupling in the two in R 2 and R 3 .

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 α ¯ ( β ¯ = 0 ) 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 i ω ¯ 3 versus y ¯ and s m ¯ 23 versus y ¯ (Figure 10(a) and Figure 10(b)) for Model A (CCM) and their comparison with similar evolution when α ¯ = 0.001 ( β ¯ = 0.0 ) i.e., model B shown in Figure 10(c) and Figure 10(d) and those for α ¯ = 0.001 and β ¯ = 0.9 in Figure 14(a) and Figure 14(b) clearly demonstrate how the free field i ω ¯ 3 and zero s m ¯ 23 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 α ¯ ( β ¯ = 0.0 ) reduces i ω ¯ 3 field due to resistance offered by the microconstituents. When both α ¯ and β ¯ are nonzero, resistance offered by the microconstituents increases resulting in further reduction in i ω ¯ 3 field. The consequences of reducing i ω ¯ 3 field are increase in s m ¯ 23 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

α 1 * Diffusion coefficient in Model RW2

α 1 , α 2 Material coefficient for constitutive theory for Cauchy moment tensor

α ¯ Nonclassical physics material coefficient

α 0 Dimensionless nonclassical physics material coefficient

β ¯ Rotational inertial physics nonclassical material coefficient

c ¯ 2 Dimensionless dissipation coefficient fluids

c d Damping coefficient Model TW1

c d * Damping coefficient Model RW1, RW2

E Modulus of elasticity

E 0 Reference modulus of elasticity

e Specific internal energy in Lagrangian description

e ¯ Specific internal energy in Eulerian description

ϵ , ϵ , ϵ i j k 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

g Heat vector in Lagrangian description

g ¯ Heat vector in Eulerian description

I Θ Internal rotational inertia in Lagrangian description

I ¯ Θ Internal rotational inertia in Eulerian description

J Deformation gradient tensor in Lagrangian description

J s Symmetric part of deformation gradient tensor in Lagrangian description

J a Skew Symmetric part of deformation gradient tensor in Lagrangian description

J d Displacement gradient tensor in Lagrangian description

J s d Symmetric part of displacement gradient tensor in Lagrangian description

J a d Skew symmetric part of displacement gradient tensor in Lagrangian description

L , L 0 Length, reference length

L , J ˙ Velocity gradient tensor in Lagrangian description

L ¯ Velocity gradient tensor in Eulerian description

s m 11 Symmetric Cauchy moment tensor in Lagrangian description component 11

s m ¯ 11 Symmetric Cauchy moment tensor in Eulerian description component 11

s m ¯ 23 Symmetric Cauchy moment tensor in Eulerian description component 23

a m ¯ 23 Skew symmetric Cauchy moment tensor in Eulerian description component 23

i ω , i ω k Internal rotational velocity or rotation rates

i ω 1 Internal rotation rate in Lagrangian description about the x1 axis

ω 0 Reference rotational velocity

i ω ¯ 1 Internal rotation rate in Eulerian description about the x1 axis

i ω ¯ 3 Internal rotation rate in Eulerian description about the x3 axis

Ω ¯ x t Space-time domain

Ω ¯ x t T Discretization of space-time domain Ω ¯ x t

Ω ¯ x t ( i ) ith space-time strip

( Ω ¯ x t ( i ) ) T Discretization of ith space-time strip

Ω ¯ x t e Space-time element e

ρ 0 Reference Density

ρ Density in Lagrangian description

ρ ¯ Density in Eulerian description

p ¯ Pressure in Eulerian description

σ Cauchy stress tensor in Lagrangian description

σ s Symmetric Cauchy stress tensor in Lagrangian description

σ a Skew symmetric Cauchy stress tensor in Lagrangian description

σ s e Symmetric equilibrium Cauchy stress tensor in Lagrangian description

σ s d Symmetric deviatoric Cauchy stress tensor in Lagrangian description

σ ¯ s d 21 ( 0 ) = σ ¯ s d 12 ( 0 ) Symmetric deviatoric contravariant Cauchy stress in Eulerian description component 12

σ ¯ a d 21 ( 0 ) = σ ¯ a d 12 ( 0 ) Antisymmetric deviatoric contravariant Cauchy stress in Eulerian description component 12

t Time

t 0 Reference time

τ 0 Reference stress

i Θ , i Θ j , { i Θ } Internal or classical rotations in Lagrangian description

i Θ ¯ , i Θ ¯ j , { i Θ ¯ } Internal or classical rotations in Eulerian description

i r Θ Internal or classical rotation rates in Lagrangian description

i r Θ ¯ Internal or classical rotation rates in Eulerian description

i Θ 1 Internal rotation about the x1 axis

i Θ 3 Internal rotation about the x3 axis

u ¯ , v ¯ , w ¯ Velocity components in x1, x2 and x3 directions

v 0 Reference velocity

v , v i , { v } Velocity vector in Lagrangian description

v ¯ , v ¯ i , { v ¯ } Velocity vector in Eulerian descriptions

x , x i , { x } Cartesian Coordinates

x ¯ , x ¯ i , { x ¯ } Cartesian Coordinates

y Cartesian Coordinates

y ¯ 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

Conflicts of Interest

The authors declare no conflicts of interest regarding the publication of this paper.

References

[1] Surana, K.S. and Mathi, S.S.C. (2023) Thermodynamic Consistency of Nonclassical Continuum Theories for Solid Continua Incorporating Rotations. Continuum Mechanics and Thermodynamics, 35, 17-59.
https://doi.org/10.1007/s00161-022-01163-y
[2] Surana, K.S. and Carranza, C.H. (2022) Nonclassical Continuum Theories for Fluent Media Incorporating Rotation Rates and Their Thermodynamic Consistency. ZAMM—Journal of Applied Mathematics and Mechanics, 103, e202200079.
https://doi.org/10.1002/zamm.202200079
[3] Surana, K.S. and Kendall, J.K. (2020) Existence of Rotational Waves in Non-Classical Thermoelastic Solid Continua Incorporating Internal Rotations. Continuum Mechanics and Thermodynamics, 32, 1659-1683.
https://doi.org/10.1007/s00161-020-00872-6
[4] Surana, K.S. and Kendall, J.K. (2022) Rotational Inertial Physics in Non-Classical Thermoviscous Fluent Continua Incorporating Internal Rotation Rates. Applied Mathematics, 13, 453-487.
https://doi.org/10.4236/am.2022.136030
[5] Islam, M., Nasrin, S. and Alam, M. (2020) Unsteady Electromagnetic Free Convection Micropolar Fluid Flow through a Porous Medium along a Vertical Porous Plate. Open Journal of Applied Sciences, 10, 701-718.
https://doi.org/10.4236/ojapps.2020.1011049
[6] Koriko, O., Omowaye, A., Popoola, A., Oreyeni, T., Adegbite, A., Oni, E., Omokhuale, E. and Altine, M. (2022) Insight into Dynamics of Hydromagnetic Flow of Micropolar Fluid Containing Nanoparticles and Gyrotactic Microorganisms at Weak and Strong Concentrations of Microelements: Homotopy Analysis Method. American Journal of Computational Mathematics, 12, 267-282.
https://doi.org/10.4236/ajcm.2022.122017
[7] Surana, K.S. and Reddy, J.N. (2017) The Finite Element Method for Initial Value Problems. CRC/Taylor and Francis, Boca Raton, FL.
https://doi.org/10.1201/b22512
[8] Surana, K.S., Powell, M.J. and Reddy, J.N. (2015) A More Complete Thermodynamic Framework for Fluent Continua. Journal of Thermal Engineering, 1, 460-475.
https://doi.org/10.18186/jte.00314
[9] Surana, K.S., Powell, M.J. and Reddy, J.N. (2015) A Polar Continuum Theory for Fluent Continua. International Journal of Engineering Research and Industrial Applications, 8, 107-146.

Copyright © 2024 by authors and Scientific Research Publishing Inc.

Creative Commons License

This work and the related PDF file are licensed under a Creative Commons Attribution 4.0 International License.