scholarly journals Universality of shear-banding instability and crystallization in sheared granular fluid

2008 ◽  
Vol 615 ◽  
pp. 293-321 ◽  
Author(s):  
MEHEBOOB ALAM ◽  
PRIYANKA SHUKLA ◽  
STEFAN LUDING
Keyword(s):  
Constitutive Model ◽  
Constitutive Relations ◽  
Shear Banding ◽  
Transport Coefficients ◽  
Navier Stokes ◽  
Dynamic Friction ◽  
Homogeneous State ◽  
Band Formation ◽  
Uniform Shear ◽  
Uniform Shear Flow

The linear stability analysis of an uniform shear flow of granular materials is revisited using several cases of a Navier–Stokes-level constitutive model in which we incorporate the global equation of states for pressure and thermal conductivity (which are accurate up to the maximum packing density νm) and the shear viscosity is allowed to diverge at a density νμ (<νm), with all other transport coefficients diverging at νm. It is shown that the emergence of shear-banding instabilities (for perturbations having no variation along the streamwise direction), that lead to shear-band formation along the gradient direction, depends crucially on the choice of the constitutive model. In the framework of a dense constitutive model that incorporates only collisional transport mechanism, it is shown that an accurate global equation of state for pressure or a viscosity divergence at a lower density or a stronger viscosity divergence (with other transport coefficients being given by respective Enskog values that diverge at νm) can induce shear-banding instabilities, even though the original dense Enskog model is stable to such shear-banding instabilities. For any constitutive model, the onset of this shear-banding instability is tied to a universal criterion in terms of constitutive relations for viscosity and pressure, and the sheared granular flow evolves toward a state of lower ‘dynamic’ friction, leading to the shear-induced band formation, as it cannot sustain increasing dynamic friction with increasing density to stay in the homogeneous state. A similar criterion of a lower viscosity or a lower viscous-dissipation is responsible for the shear-banding state in many complex fluids.

Non-Newtonian stress, collisional dissipation and heat flux in the shear flow of inelastic disks: a reduction via Grad’s moment method

2014 ◽  
Vol 757 ◽  
pp. 251-296 ◽  
Author(s):  
Saikat Saha ◽  
Meheboob Alam
Keyword(s):  
Heat Flux ◽  
Shear Flow ◽  
Moment Tensor ◽  
Transport Coefficients ◽  
Navier Stokes ◽  
Restitution Coefficient ◽  
Moment Equations ◽  
Second Moment ◽  
Uniform Shear ◽  
Uniform Shear Flow

AbstractThe non-Newtonian stress tensor, collisional dissipation rate and heat flux in the plane shear flow of smooth inelastic disks are analysed from the Grad-level moment equations using the anisotropic Gaussian as a reference. For steady uniform shear flow, the balance equation for the second moment of velocity fluctuations is solved semi-analytically, yielding closed-form expressions for the shear viscosity $\def \xmlpi #1{}\def \mathsfbi #1{\boldsymbol {\mathsf {#1}}}\let \le =\leqslant \let \leq =\leqslant \let \ge =\geqslant \let \geq =\geqslant \def \Pr {\mathit {Pr}}\def \Fr {\mathit {Fr}}\def \Rey {\mathit {Re}}\mu $, pressure $p$, first normal stress difference ${\mathcal{N}}_1$ and dissipation rate ${\mathcal{D}}$ as functions of (i) density or area fraction $\nu $, (ii) restitution coefficient $e$, (iii) dimensionless shear rate $R$, (iv) temperature anisotropy $\eta $ (the difference between the principal eigenvalues of the second-moment tensor) and (v) angle $\phi $ between the principal directions of the shear tensor and the second-moment tensor. The last two parameters are zero at the Navier–Stokes order, recovering the known exact transport coefficients from the present analysis in the limit $\eta ,\phi \to 0$, and are therefore measures of the non-Newtonian rheology of the medium. An exact analytical solution for leading-order moment equations is given, which helped to determine the scaling relations of $R$, $\eta $ and $\phi $ with inelasticity. We show that the terms at super-Burnett order must be retained for a quantitative prediction of transport coefficients, especially at moderate to large densities for small values of the restitution coefficient ($e \ll 1$). Particle simulation data for a sheared inelastic hard-disk system are compared with theoretical results, with good agreement for $p$, $\mu $ and ${\mathcal{N}}_1$ over a range of densities spanning from the dilute to close to the freezing point. In contrast, the predictions from a constitutive model at Navier–Stokes order are found to deviate significantly from both the simulation and the moment theory even at moderate values of the restitution coefficient ($e\sim 0.9$). Lastly, a generalized Fourier law for the granular heat flux, which vanishes identically in the uniform shear state, is derived for a dilute granular gas by analysing the non-uniform shear flow via an expansion around the anisotropic Gaussian state. We show that the gradient of the deviatoric part of the kinetic stress drives a heat current and the thermal conductivity is characterized by an anisotropic second-rank tensor, for which explicit analytical expressions are given.

Analysis of gas flow over a cylinder at all Knudsen numbers based on non-newtonian constitutive models

2020 ◽  
Vol 34 (14n16) ◽  
pp. 2040076
Author(s):  
Zhen-Yu Yuan ◽  
Zhong-Zheng Jiang ◽  
Wen-Wen Zhao ◽  
Wei-Fang Chen
Keyword(s):  
Constitutive Model ◽  
Gas Flow ◽  
Constitutive Relations ◽  
Constitutive Models ◽  
Direct Simulation Monte Carlo ◽  
Navier Stokes ◽  
Dsmc Simulation ◽  
Nonequilibrium Flows ◽  
Simulation Results ◽  
Better Than

This paper is focused on the gas properties over a cylinder from continuum to rarefied regimes based on the non-Newtonian constitutive model. This new constitutive model is first derived from Eu’s nonequilibrium ensemble method, which is intended for accurate description of nonequilibrium flows. Some assumptions and simplifications are made during the establishing progress of the new constitutive model by both Eu and Myong. To verify its accuracy, temperature contours and skin frictions around the cylinder are simulated by this new model. The inflow Mach number is equal to 10 and the Knudsen number ranges from 0.002 to 0.05. All simulation results are compared with Navier–Stokes (NS) and the direct simulation Monte Carlo (DSMC) methods in detail. The comparisons of friction around the surface show that the non-Newtonian constitutive models are better than the linear constitutive relations of NS equations for the prediction of nonequilibrium flow and much more close to DSMC simulation results.

Asymptotic analysis of a linearized trailing edge flow

1972 ◽  
Vol 6 (3) ◽  
pp. 327-347 ◽  
Author(s):  
K. Capell
Keyword(s):  
Boundary Layer ◽  
Stokes Equations ◽  
Navier Stokes ◽  
Wake Region ◽  
Trailing Edge ◽  
Navier Stokes Equations ◽  
Uniform Shear ◽  
Uniform Shear Flow ◽  
Appl Math ◽  
Edge Flow

An Oseén type linearization of the Navier-Stokes equations is made with respect to a uniform shear flow at the trailing edge of a flat plate. Asymptotic expansions are obtained to describe a symmetrical merging flow for distances from the trailing edge that are, in a certain sense, large. Expansions for three regions are found:(i) a wake region,(ii) an inviscid region, and(iii) an upstream lower order boundary layer.The results are compared with those of Hakkinen and O'Neil (Douglas Aircraft Co. Report, 1967) and Stewartson (Proc. Roy. Soc. Ser. A 306 (1968)). They are further related to the results of Stewartson (Mathematika 16 (1969)) and Messiter (SIAM J. Appl. Math. 18 (1970)).

Normal stress differences, their origin and constitutive relations for a sheared granular fluid

2016 ◽  
Vol 795 ◽  
pp. 549-580 ◽  
Author(s):  
Saikat Saha ◽  
Meheboob Alam
Keyword(s):  
Normal Stress ◽  
Constitutive Relations ◽  
Moment Tensor ◽  
Normal Stress Difference ◽  
Stress Difference ◽  
Uniform Shear ◽  
Uniform Shear Flow ◽  
Sign Change ◽  
Normal Stress Differences ◽  
Dilute Limit

The rheology of the steady uniform shear flow of smooth inelastic spheres is analysed by choosing the anisotropic/triaxial Gaussian as the single-particle distribution function. An exact solution of the balance equation for the second-moment tensor of velocity fluctuations, truncated at the ‘Burnett order’ (second order in the shear rate), is derived, leading to analytical expressions for the first and second ($\unicode[STIX]{x1D615}_{1}$ and $\unicode[STIX]{x1D615}_{2}$) normal stress differences and other transport coefficients as functions of density (i.e. the volume fraction of particles), restitution coefficient and other control parameters. Moreover, the perturbation solution at fourth order in the shear rate is obtained which helped to assess the range of validity of Burnett-order constitutive relations. Theoretical expressions for both $\unicode[STIX]{x1D615}_{1}$ and $\unicode[STIX]{x1D615}_{2}$ and those for pressure and shear viscosity agree well with particle simulation data for the uniform shear flow of inelastic hard spheres for a large range of volume fractions spanning from the dilute regime to close to the freezing-point density (${\it\nu}\sim 0.5$). While the first normal stress difference $\unicode[STIX]{x1D615}_{1}$ is found to be positive in the dilute limit and decreases monotonically to zero in the dense limit, the second normal stress difference $\unicode[STIX]{x1D615}_{2}$ is negative and positive in the dilute and dense limits, respectively, and undergoes a sign change at a finite density due to the sign change of its kinetic component. It is shown that the origin of $\unicode[STIX]{x1D615}_{1}$ is tied to the non-coaxiality (${\it\phi}\neq 0$) between the eigendirections of the second-moment tensor $\unicode[STIX]{x1D648}$ and those of the shear tensor $\unicode[STIX]{x1D63F}$. In contrast, the origin of $\unicode[STIX]{x1D615}_{2}$ in the dilute limit is tied to the ‘excess’ temperature ($T_{z}^{ex}=T-T_{z}$, where $T_{z}$ and $T$ are the $z$-component and the average of the granular temperature, respectively) along the mean vorticity ($z$) direction, whereas its origin in the dense limit is tied to the imposed shear field.

scholarly journals Alternative derivation of the higher-order constitutive model for six-parameter elastic shells

2021 ◽  
Vol 72 (2) ◽  
Author(s):  
Mircea Bîrsan
Keyword(s):  
Energy Density ◽  
Constitutive Model ◽  
Strain Energy Density ◽  
Strain Energy ◽  
Shell Theory ◽  
Constitutive Relations ◽  
Three Dimensional ◽  
Classical Shell Theory ◽  
Three Dimensional Elasticity ◽  
General Method

AbstractIn this paper, we present a general method to derive the explicit constitutive relations for isotropic elastic 6-parameter shells made from a Cosserat material. The dimensional reduction procedure extends the methods of the classical shell theory to the case of Cosserat shells. Starting from the three-dimensional Cosserat parent model, we perform the integration over the thickness and obtain a consistent shell model of order $$ O(h^5) $$ O ( h 5 ) with respect to the shell thickness h. We derive the explicit form of the strain energy density for 6-parameter (Cosserat) shells, in which the constitutive coefficients are expressed in terms of the three-dimensional elasticity constants and depend on the initial curvature of the shell. The obtained form of the shell strain energy density is compared with other previous variants from the literature, and the advantages of our constitutive model are discussed.

Formation of Vortex Street in Laminar Boundary Layer

1980 ◽  
Vol 47 (2) ◽  
pp. 227-233 ◽  
Author(s):  
M. Kiya ◽  
M. Arie
Keyword(s):  
Boundary Layer ◽  
Shear Flow ◽  
Laminar Boundary Layer ◽  
Solid Wall ◽  
Vortex Sheet ◽  
Vortex Street ◽  
Bluff Bodies ◽  
Uniform Shear ◽  
Uniform Shear Flow ◽  
Vortex Patterns

Main features of the formation of vortex street from free shear layers emanating from two-dimensional bluff bodies placed in uniform shear flow which is a model of a laminar boundary layer along a solid wall. This problem is concerned with the mechanism governing transition induced by small bluff bodies suspended in a laminar boundary layer. Calculations show that the background vorticity of shear flow promotes the rolling up of the vortex sheet of the same sign whereas it decelerates that of the vortex sheet of the opposite sign. The steady configuration of the conventional Karman vortex street is not possible in shear flow. Theoretical vortex patterns are experimentally examined by a flow-visualization technique.

The formation mechanism and shedding frequency of vortices from a sphere in uniform shear flow

1995 ◽  
Vol 287 ◽  
pp. 151-171 ◽  
Author(s):  
Hiroshi Sakamoto ◽  
Hiroyuki Haniu
Keyword(s):  
Reynolds Number ◽  
Shear Flow ◽  
Formation Mechanism ◽  
Vortex Shedding ◽  
Uniform Flow ◽  
Sphere Diameter ◽  
Uniform Shear ◽  
Uniform Shear Flow ◽  
Vortex Loops ◽  
Shear Parameter

Experiments to investigate the formation mechanism and frequency of vortex shedding from a sphere in uniform shear flow were conducted in a water channel using flow visualization and velocity measurement. The Reynolds number, defined in terms of the sphere diameter and approach velocity at its centre, ranged from 200 to 3000. The shear parameter K, defined as the transverse velocity gradient of the shear flow non-dimensionalized by the above two parameters, was varied from 0 to 0.25. The critical Reynolds number beyond which vortex shedding from the sphere occurred was found to be lower than that for uniform flow and decreased approximately linearly with increasing shear parameter. Also, the Strouhal number of the hairpin-shaped vortex loops became larger than that for uniform flow and increased as the shear parameter increased.The formation mechanism and the structure of vortex shedding were examined on the basis of series of photographs and subsequent image processing using computer graphics. The range of Reynolds number in the present investigation, extending up to 3000, could be classified into three regions on the basis of this study, and it was observed that the wake configuration did not differ substantially from that for uniform flow. Also, unlike the detachment point of vortex loops in uniform flow, which was irregularly located along the circumference of the sphere, the detachment point in shear flow was always on the high-velocity side.

The heat/mass transfer to a finite strip at small Péclet numbers

1978 ◽  
Vol 86 (1) ◽  
pp. 49-65 ◽  
Author(s):  
R. C. Ackerberg ◽  
R. D. Patel ◽  
S. K. Gupta
Keyword(s):  
Heat Transfer ◽  
Mass Transfer ◽  
Shear Flow ◽  
Sodium Hydroxide Solution ◽  
Flow Direction ◽  
Mathematical Problem ◽  
Limiting Current ◽  
Transfer Rates ◽  
Uniform Shear ◽  
Uniform Shear Flow

The problem of heat transfer (or mass transfer at low transfer rates) to a strip of finite length in a uniform shear flow is considered. For small values of the Péclet number (based on wall shear rate and strip length), diffusion in the flow direction cannot be neglected as in the classical Leveque solution. The mathematical problem is solved by the method of matched asymptotic expansions and expressions for the local and overall dimensionless heat-transfer rate from the strip are found. Experimental data on wall mass-transfer rates in a tube at small Péclet numbers have been obtained by the well-known limiting-current method using potassium ferrocyanide and potassium ferricyanide in sodium hydroxide solution. The Schmidt number is large, so that a uniform shear flow can be assumed near the wall. Experimental results are compared with our theoretical predictions and the work of others, and the agreement is found to be excellent.

Sign in / Sign up

Export Citation Format

Share Document