Effect of tangential interface motion on the viscous instability in fluid flow past flexible surfaces

The stability of linear shear flow of a Newtonian fluid past a flexible membrane is analysed in the limit of low Reynolds number as well as in the intermediate Reynolds number regime for two different membrane models. The objective of this paper is to demonstrate the importance of tangential motion in the membrane on the stability characteristics of the shear flow. The first model assumes the wall to be a “spring-backed” plate membrane, and the displacement of the wall is phenomenologically related in a linear manner to the change in the fluid stresses at the wall. In the second model, the membrane is assumed to be a two-dimensional compressible viscoelastic sheet of infinitesimal thickness, in which the constitutive relation for the shear stress contains an elastic part that depends on the local displacement field and a viscous component that depends on the local velocity in the membrane. The stability characteristics of the laminar flow in the limit of low are crucially dependent on the tangential motion in the membrane wall. In both cases, the flow is stable in the low Reynolds number limit in the absence of tangential motion in the membrane. However, the presence of tangential motion in the membrane destabilises the shear flow even in the absence of fluid inertia. In this case, the non-dimensional velocity (Λ_t) required for unstable fluctuations is proportional to the wavenumber k ( Λ _t∼k) in the plate membrane type of wall while it scales as k² in the viscoelastic membrane type of wall ( Λ _t∼k ²) in the limit k→ 0. The results of the low Reynolds number analysis are extended numerically to the intermediate Reynolds number regime for the case of a viscoelastic membrane. The numerical results show that for a given set of wall parameters, the flow is unstable only in a finite range of Reynolds number, and it is stable in the limit of large Reynolds number. Received 8 November 2000 and Received in final form 20 March 2001     

Asymptotic analysis of wall modes in a flexible tube

The stability of wall modes in a flexible tube of radius R surrounded by a viscoelastic material in the region R < r < H R in the high Reynolds number limit is studied using asymptotic techniques. The fluid is a Newtonian fluid, while the wall material is modeled as an incompressible visco-elastic solid. In the limit of high Reynolds number, the vorticity of the wall modes is confined to a region of thickness in the fluid near the wall of the tube, where the small parameter , and the Reynolds number is , and are the fluid density and viscosity, and V is the maximum fluid velocity. The regime is considered in the asymptotic analysis, where G is the shear modulus of the wall material. In this limit, the ratio of the normal stress and normal displacement in the wall, , is only a function of H and scaled wave number . There are multiple solutions for the growth rate which depend on the parameter .In the limit , which is equivalent to using a zero normal stress boundary condition for the fluid, all the roots have negative real parts, indicating that the wall modes are stable. In the limit , which corresponds to the flow in a rigid tube, the stable roots of previous studies on the flow in a rigid tube are recovered. In addition, there is one root in the limit which does not reduce to any of the rigid tube solutions determined previously. The decay rate of this solution decreases proportional to in the limit , and the frequency increases proportional to . Received: 5 November 1997 / Revised: 10 March 1998 / Accepted: 29 April 1998     

微多孔介质非牛顿流体格子Boltzmann模拟

根据修正Blake-Kozeny方程对Herschel-Bulkley模型的非牛顿流体在微多孔介质中流动的定量描述,通过增加外力项建立了非牛顿流体在宏观多孔介质中的格子Boltzmann模型,并将该模型应用到平板多孔介质中非牛顿流体的电渗驱动流模拟。对多孔介质渗透性参数和非牛顿流体流变性参数的影响进行模拟和分析,获得了各种参数对多孔介质中非牛顿流体流动速度的影响规律。     

Elastic waves in non-Newtonian (Maxwell) fluid-saturated porous media

Abstract

This paper studies the elastic waves in non-Newtonian (Maxwell) fluid-saturated porous media with the nonzero boundary slip velocity for pore size distribution. The coefficient bF _m (ω) that measures the deviation from Poiseuille flow friction in such media is presented. Based on this coefficient, we investigate the properties of elastic waves by calculating their phase velocities and attenuation coefficients as functions of frequency and the behaviour of the dynamic permeability. The study shows that the pore size distribution removes oscillations in all physical quantities in the non-Newtonian regime. Consideration of the nonzero boundary slip effect in non-Newtonian (Maxwell) fluid-saturated porous media results in (a) an overall increase of the dynamic permeability, (b) an increase of phase velocities of fast Biot waves and shear waves except in the low frequency domain and an overall increase of phase velocity of slow Biot waves and (c) an overall increase of the attenuation of three Biot waves in the intermediate frequency domain except in the deeply non-Newtonian regime. The study also shows that the attenuation coefficient of slow Biot waves is small in the deeply non-Newtonian regime at higher frequency, which encourages us to detect slow Biot waves in oil-saturated porous rock.     

Magnetohydrodynamic flow of a viscoelastic fluid

The non-linear differential equation for the magnetohydrodynamic Poiseuille flow of Phan-Thein-Tanner (PTT) conducting fluid is derived. Using the homotopy analysis method (HAM), the series solution is developed and its convergence is discussed. Also, the results are presented graphically and the effects of non-dimensional parameters on the flow field are analyzed. The results obtained reveal many interesting behaviors that warrant further study on the equations related to non-Newtonian fluid phenomena, especially the shear-thinning phenomena. Shear thinning reduces the wall shear stress.     

Plasticity and dynamical heterogeneity in driven glassy materials

Many amorphous glassy materials exhibit complex spatio-temporal mechanical response and rheology, characterized by an intermittent stress strain response and a fluctuating velocity profile. Under quasistatic and athermal deformation protocols this heterogeneous plastic flow was shown to be composed of plastic events of various sizes, ranging from local quadrupolar plastic rearrangements to system spanning shear bands. In this paper, through numerical study of a 2D Lennard-Jones amorphous solid, we generalize the study of the heterogeneous dynamics of glassy materials to the finite shear rate ( [(g)\dot] \dot{{\gamma}} 1 \neq 0 and temperature case (T 1 \neq 0 . In practice, we choose an effectively athermal limit (T ∼ 0 and focus on the influence of shear rate on the rheology of the glass. In line with previous works we find that the model Lennard-Jones glass follows the rheological behavior of a yield stress fluid with a Herschel-Bulkley response of the form, s \sigma = s_Y \sigma_{{Y}}^{} + c ₁ [(g)\dot]^b \dot{{\gamma}}^{{\beta}}_{} . The global mechanical response obtained through the use of Molecular Dynamics is shown to converge in the limit [(g)\dot] \dot{{\gamma}} ? \rightarrow 0 to the quasistatic limit obtained with an energy minimization protocol. The detailed analysis of the plastic deformation at different shear rates shows that the glass follows different flow regimes. At sufficiently low shear rates the mechanical response reaches a shear-rate-independent regime that exhibits all the characteristics of the quasistatic response (finite-size effects, cascades of plastic rearrangements, yield stress, ...). At intermediate shear rates the rheological properties are determined by the externally applied shear rate and the response deviates from the quasistatic limit. Finally at higher shear the system reaches a shear-rate-independent homogeneous regime. The existence of these three regimes is also confirmed by the detailed analysis of the atomic motion. The computation of the four-point correlation function shows that the transition from the shear-rate-dominated to the quasistatic regime is accompanied by the growth of a dynamical cooperativity length scale x \xi that is shown to diverge with shear rate as x \xi μ \propto [(g)\dot]^-n \dot{{\gamma}}^{{-\nu}}_{} , with n \nu ∼ 0.2 -0.3. This scaling is compared with the prediction of a simple model that assumes the diffusive propagation of plastic events.     

Analysis of ferroelectric switching in sodium nitrite:poly(vinyl alcohol) nanocomposite films

The Fourier transform infrared (FTIR) spectra and switching current response in sodium nitrite:poly(vinyl alcohol) nanocomposite films have been studied as a function of composition of NaNO₂. The switching current data fitted well to infinite-grain model (IGM) in the region t<t _s and to finite-grain model (FGM) in the region t≥t _s . The microscopic parameters like the dimensionality, the domain wall velocity, and the nucleation rate have been evaluated which provide more physical insight of the switching phenomena in the composite films. The polarization current and nucleation rate are optimum in 50 wt.% composite film and have been discussed in terms of grain size and strain variations with the composition. The effect of applied field and pulse width variation on the switching behavior of 50 wt.% composition has also been studied. The exponential field dependence of the domain wall velocity and the nucleation rate indicate that nucleation mechanism is responsible for switching phenomena in the composite films. The writing pulse width affects significantly on the switching behavior of the composite films.     

Thermal boundary layer for non-Newtonian fluid

Summary Analytical and numerical solutions for the momentum and thermal boundary layer equations of a non-Newtonian power law fluid are presented. The flow is assumed to be under the influence of an external magnetic fieldB (x) applied perpendicular to the surface and an electric fieldE(x) perpendicular toB(x) and the direction of the longitudinal velocity in the boundary layer. For the power law fluid it is assumed that the shear stress is proportional to then-th power of the velocity gradient andn is called the flow index. The variations of the velocity fieldf′, the temperature field θ, the shear stress on the surfaceτ _W , the displacement thicknessδ ₁ and the momentum thicknessδ ₂ with the magnetic-field parameter γ, the flow indexn, the heat transfer indexS and the Prandtl number Pr are studied. It is found that, if the outer flow velocityU(x) (potential flow) is proportional to the arc lengthx raised to a powerm, then the similarity solution for the thermal boundary layer equation is possible only whenm=1/3, which represents flow past a wedge of included angle π/2. It is established that the temperature of the wedge increases with the increase of γ, Pr,S and the decrease ofn. In general the magnetic field can be used as a heater for the surface of the wedge.     

Analysis of peristaltic flow of Jeffrey six constant nano fluid in a vertical non-uniform tube

Peristaltic flow of non-Newtonian nano fluid through a non-uniform surface has been investigated in this paper. The fluid motion along the wall of the surface is caused by the sinusoidal wave traveling with constant speed. The governing equations are converted into cylindrical coordinate system and assuming low Reynolds number and long wave length partial differential equations are simplified. Analytically solutions of the problem are obtained by utilizing the homotopy perturbation method (HPM). In order to insight the impact of embedded parameters on temperature, concentration and velocity some graphs are plotted for different peristaltic waves. At the end, some observations were made from the graphical presentation that velocity, pressure rise and nano particle concentration are increasing function of thermophoresis parameter N_t while temperature and frictional forces show opposite trend.     

Simulation of the Second Grade Fluid Model for Blood Flow through a Tapered Artery with a Stenosis

We analyze the blood flow through a tapered artery, assuming the blood to be a second order fluid model. The resulting nonlinear implicit system of partial differential equations is solved by the perturbation method. The expressions for shear stress, velocity, flow rate, wall shear stress and longitudinal impedance are obtained. The physical behavior of different parameters is also discussed, as are trapping phenomena.     

Study on blood flow containing nanoparticles through porous arteries in presence of magnetic field using analytical methods

In this paper, flow analysis for a third grade non-Newtonian blood in porous arteries in presence of magnetic field is simulated analytically and numerically. Blood is considered as the third grade non-Newtonian fluid containing nanoparticles. Collocation Method (CM) and Optimal Homotopy Asymptotic Method (OHAM) are used to solve the Partial Differential Equation (PDE) governing equation which a good agreement between them was observed in the results. The influences of the some physical parameters such as Brownian motion parameter, pressure gradient and thermophoresis parameter, etc. on temperature, velocity and nanoparticles concentration profiles are considered. For instance, increasing the thermophoresis parameter (N_t) caused an increase in temperature values in whole domain and an increase in nanoparticles concentration near the inner wall.     

NMR velocimetry studies of the steady-shear rheology of a concentrated hard-sphere colloidal system

NMR velocimetry has been used to observe the steady-shear rheological behaviour of a concentrated suspension of hard-sphere like 370 nm diameter PMMA core-shell latex particles at the volume fraction Φ = 0.46, the liquid core of the spheres rendering possible NMR observation of the particles themselves. Rheological measurements in a cone-and-plate geometry indicate that when aged (i.e. left at rest for two weeks), the material exhibits yield stress behaviour at very low shear rates. For shear rates greater than 1 s ^{- 1} a transition to liquid-like behaviour was observed, leading to a rejuvenated fluid state which exhibits shear-thinning behaviour over a wide range of shear rates. A similar yield stress behaviour was reflected in NMR velocimetry measurements in a Couette geometry, where the solid-to liquid transition could be clearly observed. Under steady-state flow, the fluid state inside the radius at which yield stress was observed, exhibited shear-thinning behaviour with a power law exponent n slowly approaching unity with increasing shear rate. This behaviour has some similarities with a model of Derec et al. in which aging and rejuvenation effects compete. Substantial wall slip was observed both at the inner and at the outer wall, an effect which disappeared as the shear rate was increased. No radial particle migration from the high-shear region at the inner wall was observed.     

Strengthened opposition control for skin-friction reduction in wall-bounded turbulent flows

An opposition control scheme with strengthened control input is proposed and tested in turbulent channel flows at friction Reynolds number Re_τ = 180 by direct numerical simulations. When the detection plane is located at less than 20 wall units, the drag reduction rate can be greatly enhanced by increasing the control amplitude parameter. The maximum drag reduction rate achieved in the present study is around 33%, which is much higher than the best value of 25% reported in literature. The strengthened control can be more efficient to attain a given drag reduction rate. Based on the total shear stress at the virtual wall established between the real wall and the detection plane by the control, a new friction velocity is proposed and the corresponding coordinate transform is made. Scaled by the proposed friction velocity, the wall-normal velocity fluctuation and the Reynolds shear stress of the controlled flows are collapsed well with those of the uncontrolled flow in the new coordinate. Based on the similarity, a relation between drag reduction rate and the effectiveness of the virtual wall is deduced, which disclosed that the elevation and residual Reynolds shear stress at the virtual wall are the key parameters to determine the drag reduction rate. The conclusion are also validated at Re_τ = 395 and 590. The decrease of the drag reduction rate with the increase of the Reynolds number is attributed to the enhanced residual Reynolds shear stress at the virtual wall.     

Large eddy simulations and parameterisation of roughness element orientation and flow direction effects in rough wall boundary layers

We conduct a series of large eddy simulations (LES) of turbulent boundary layers over arrays of cuboidal roughness elements at arbitrary orientation angles (non-frontal orientations with the incident flow). Flow response to changing roughness orientation is systematically studied at two ground coverage densities, λ_p = 0.06 and 0.11. As expected, the effective roughness heights z_o measured from LES are higher for λ_p = 0.11 than for λ_p = 0.06, although appreciable changes both in z_o and wall shear stress (friction velocity) are observed at both ground coverage densities as the roughness orientation angle changes. This suggests the necessity of accounting for detailed rough wall topology (including more information than just λ_p, λ_f) when relating rough wall morphology to its aerodynamic properties. To this end, a recently developed analytical rough wall parameterisation is used to predict the aerodynamic properties of the simulated rough surfaces. In this rough wall model, wake interactions among roughness elements are explicitly modelled using the concept of sheltering height and exponential attenuation coefficient. As a result, the parameterisation is responsive to detailed ground roughness arrangements and flow conditions, including roughness height variations, element orientation, incident flow direction, transverse displacements, etc. Model-predicted effective roughness heights, wall stress, mean velocity at the height of the roughness, and in some cases displacement height, are compared against the LES measurements from this study as well as numerical/experiment measurements from other authors. The predictions from the model are found to agree well with the measurements both in trends and in absolute values, thus extending the applicability of the analytical rough wall model to more general surfaces than those previously tested.     

Scaling study of the combustion performance of gasben gas rocket injectors

To obtain the key subelements that may influence the scaling of gas-gas injector combustor performance, the combustion performance subelements in a liquid propellant rocket engine combustor are initially analysed based on the results of a previous study on the scaling of a gas-gas combustion flowfield. Analysis indicates that inner wall friction loss and heat-flux loss are two key issues in gaining the scaling criterion of the combustion performance. The similarity conditions of the inner wall friction loss and heat-flux loss in a gas-gas combustion chamber are obtained by theoretical analyses. Then the theoretical scaling criterion was obtained for the combustion performance, but it proved to be impractical. The criterion conditions, the wall friction and the heat flux are further analysed in detail to obtain the specific engineering scaling criterion of the combustion performance. The results indicate that when the inner flowfields in the combustors are similar, the combustor wall shear stress will have similar distributions qualitatively and will be directly proportional to p_c^0.8 d_t^-0.2 quantitatively. In addition, the combustion peformance will remain unchanged. Furthermore, multi-element injector chambers with different geometric sizes and at different pressures are numerically simulated and the wall shear stress and combustion efficiencies are solved and compared with each other. A multi-element injector chamber is designed and hot-fire tested at several chamber pressures and the combustion performances are measured in a total of nine hot-fire tests. The numerical and experimental results verified the similarities among combustor wall shear stress and combustion performances at different chamber pressures and geometries, with the criterion applied.     

Tracer dispersion in power law fluids flow through porous media: Evidence of a cross-over from a logarithmic to a power law behaviour

An analytical model is presented to describe the dispersion of tracers in a power-law fluid flowing through a statistically homogeneous and isotropic porous medium. The model is an extension of Saffman's approach to the case of non-Newtonian fluids. It is shown that an effective viscosity depending on the pressure gradient and on the characteristics of the fluid, must be introduced to satisfy Darcy's law. An analytical expression of the longitudinal dispersivity is given as a function of the Peclet number Pe and of the power-law index n that characterizes the dependence of the viscosity on the shear rate . As the flow velocity increases the dispersivity obeys an asymptotic power law: . This asymptotic regime is achieved at moderate Peclet numbers with strongly non-Newtonian fluids and on the contrary at very large values when n goes to 1 ( for n=0.9). This reflects the cross-over from a scaling behaviour for towards a logarithmic behaviour predicted for Newtonian fluids (n=1). Received: 22 July 1997 / Revised and Accepted: 2 July 1998     

Heat Transfer in an Upper Convected Maxwell Fluid with Fluid Particle Suspension

An analysis is carried out to study the magnetohydrodynamic (MHD) flow and heat transfer characteristics of an electrically conducting dusty non-Newtonian fluid, namely, the upper convected Maxwell (UCM) fluid over a stretching sheet. The stretching velocity and the temperature at the surface are assumed to vary linearly with the distance from the origin. Using a similarity transformation, the governing nonlinear partial differential equations of the model problem are transformed into coupled non-linear ordinary differential equations and the equations are solved numerically by a second order finite difference implicit method known as the Keller-box method. Comparisons with the available results in the literature are presented as a special case. The effects of the physical parameters on the fluid velocity, the velocity of the dust particle, the density of the dust particle, the fluid temperature, the dust-phase temperature, the skin friction, and the wall-temperature gradient are presented through tables and graphs. It is observed that, Maxwell fluid reduces the wall-shear stress. Also, the fluid particle interaction reduces the fluid temperature in the boundary layer. Furthermore, the results obtained for the flow and heat transfer characteristics reveal many interesting behaviors that warrant further study on the non-Newtonian fluid flow phenomena, especially the dusty UCM fluid flow phenomena.     

HEAT AND MASS TRANSFER AND WALL SHEAR STRESS IN VERTICAL GAS-LIQUID FLOW

Results of an experimental investigation of heat and mass transfer and wall shear stress at gas-liquid flow in a vertical tube are presented. Local wall shear stress and mass transfer coefficients were measured by an electrochemical method. Experiments were performed in the range of Reynolds number variation with respect to liquid Rc_i, = 8.5 × 10³-5.4 × 10⁴, gas Re_g = 3 × 10³-1.4 × 10⁵, pressure 0.1-1 MPa. The relationship between heat and mass transfer and wall shear at gas-liquid flows is shown to exist. The results of measuring heat and mass transfer coefficients are generalized by formulas applied to calculate heat and mass transfer in single-phase turbulent flow.     

Investigation of scaling laws in a turbulent boundary layer flow with adverse pressure gradient using PIV

We present an experimental investigation and data analysis of a turbulent boundary layer flow at a significant adverse pressure gradient at Reynolds number up to Re_θ = 10, 000. We combine large-scale particle image velocimetry (PIV) with microscopic PIV for measuring the near wall region including the viscous sublayer. We investigate scaling laws for the mean velocity and for the total shear stress in the inner part of the boundary layer. In the inner part the mean velocity can be fitted by a log-law. In the outer part of the inner layer the log-law ceases to be valid. Instead, a modified log-law provides a good fit, which is given in terms of the pressure gradient parameter and a parameter for the mean inertial effects. Finally we describe and assess a simple quantitative model for the total shear stress distribution which is local in wall-normal direction without streamwise history effects.     

Peculiarities of deformation of compliant coatings with increased strength

Influence of elasticity module of coating material on the parameters of hard compliant coatings deformation has been analysed. Calculation using two-dimensional model has shown that maximum coating deformation is achieved at the ratio of flow rate U to the parameter C _t ⁰ = (E/3ρ)^0.5 approximately equal to 2.5, however, velocity of wall surface motion has first local maximum at U/C _t ⁰≈1. The range of coating parameters’ values at which compromise between its hardness and intensity of interaction with turbulent flow is provided has been determined. For rubbery materials with Poisson coefficient of about 0.5, correlations of the flow velocity and parameter C _t ⁰ shall be in the range 1÷1.5. It is shown that at such parameters, the mean square value of the coating surface deflection/inflection is less than the viscous sublayer thickness, its correlation with the wavelength is very small and equals (1÷5)·10⁻⁴. Such form of deformed surface fundamentally differs from the parameters of the wave wall in Kendall’s experiments which results are used for calculation of inverse influence of wall deformation on the flow. It was assumed that solid compliant coatings do not cause instability of interaction with the ambient flow.