Axisymmetric stagnation point flow and mass transfer for power-law fluids II. Moderate schmidt number correction

The theory for axisymmetric stagnation point flow of power-law fluids has been extended to include the correction terms for convective diffusion at moderate Schmidt numbers. The dimensionless mass transfer rate is expressed as an asymptotic series that is valid for Re^{(1 ? n)/3(1 + n)}Sc^{?$\text{1}\text{3}$} < 1. The result can be used to predict accurate diffusion coefficients for dilute species in fluids with specified power-law characteristics.     

OVERLAPPING CONTROL VOLUME APPROACH FOR CONVECTION–DIFFUSION PROBLEMS

This paper introduces a finite volume method to solve 2D steady state convection–diffusion problems on structured non-orthogonal grids. Overlapping control volumes (OCV) are used to discretize the physical domain and the governing equations are solved without transformation. An isoparametric formulation is used to compute diffusion and for upwinding. Four test problems are solved using this and other schemes. The modelling of diffusion in OCV seems very effective even on distorted meshes. The convection modelling in OCV is found to be second-order-accurate, like QUICK, on regular meshes. Although its accuracy is slightly inferior to the latter on rectangular grids, its faster convergence gives it a better overall performance. On non-orthogonal grids, OCV gives better accuracy for a large and practical range of Peclet numbers than does QUICK applied to the transformed equations using the conventional five-point diffusion modelling. The results obtained also demonstrate that the scheme reduces false diffusion to a considerable extent in comparison with the power-law scheme.     

Steady flow of a power-law non-Newtonian fluid across an unconfined square cylinder

A two-dimensional flow of a non-Newtonian power-law fluid directed normally to a horizontal cylinder with a square cross section is considered in the present paper. The problem is investigated numerically with a finite volume method by using the commercial code Ansys Fluent with a very large computational domain so that the flow could be considered unbounded. The investigation covers the power-law index from 0.1 to 2.0 and the Reynolds number range from 0.001 to 45.000. It is found that the drag coefficient for low Reynolds numbers and low power-law index (n ≤ 0.5) obeys the relationship C_D = A/Re. An equation for the quantity A as a function of the power-law index is derived. The drag coefficient becomes almost independent of the power-law index at high Reynolds numbers and the wake length changes nonlinearly with the Reynolds number and power-law index.     

Heat and mass transfer to a reacting particle in a gas flow for arbitrary temperature dependence of the transport coefficients

A study is made of simultaneous heat and mass transfer to a reacting particle of any shape in a translational (and shear) flow of a viscous heat conducting compressible gas for which the thermal conductivity and diffusion coefficient, and also the specific heat depend on the temperature. The first two terms of the asymptotic expansion with respect to the small Reynolds number are obtained for the mean Sherwood and Nusselt numbers. The case of a power-law temperature dependence of the gas viscosity is considered in detail.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 111–119, January–February, 1984.     

Non-Newtonian Behavior of Ballistic Gelatin at High Shear Rates

A coordinated modeling and experimental effort to investigate the shear stress-shear strain rate response of ballistic gelatin is presented. A power-law constitutive model that captures non-Newtonian shear-thickening behavior, the evolution of viscosity, and the momentum diffusion at high shear rates is adopted. A simple asymptotic relationship between the maximum wall shear stress and the maximum striking wall velocity is derived in the high diffusion rate regime for a shear flow between two parallel plates. Experimental investigation is conducted on double lap-shear test fixture with gelatin specimens of different thicknesses subjected to high strain rate input on the inner surface, generated by a polymer split Hopkinson pressure bar. This test fixture allows measurement of transmitted shear stress as well as visualization of momentum diffusion through gelatin when imaged by a high speed camera. Gelatin specimens of various thicknesses were used for extracting the power-law model parameters. It is found that ballistic gelatin behaves as a shear-thickening fluid at high shear rates with a power-law exponent of 2.22.     

Variable viscosity and slip velocity effects on the flow and heat transfer of a power-law fluid over a non-linearly stretching surface with heat flux and thermal radiation

The flow and heat transfer of a non-Newtonian power-law fluid over a non-linearly stretching surface has been studied numerically under conditions of constant heat flux and thermal radiation and evaluated for the effect of wall slip. The governing partial differential equations are transformed into a set of coupled non-linear ordinary differential equations which are using appropriate boundary conditions for various physical parameters. The remaining set of ordinary differential equations is solved numerically by fourth-order Runge–Kutta method using the shooting technique. The effects of the viscosity, the slip velocity, the radiation parameter, power-law index, and the Prandtl number on the flow and temperature profiles are presented. Moreover, the local skin friction and Nusselt numbers are presented. Comparison of numerical results is made with the earlier published results under limiting cases.     

The Dean instability in power-law and Bingham fluids in a curved rectangular duct

The laminar flow of power-law and yield-stress fluids in 180° curved channels of rectangular cross section was studied experimentally and numerically in order to understand the effect of rheological fluid behavior on the Dean instability that appears beyond a critical condition in the flow. This leads to the apparition of Dean vortices that differ from the two corner vortices created by the channel wall curvature.Flow visualizations showed that the Dean vortices develop first in the near-wall zone on the concave (outer) wall, where the shear rate is higher and the viscosity weaker; then they penetrate into the centre of the channel cross section where power-law fluids have high viscosity and Bingham fluids are unyielded in laminar flow. Based on the complete formation on the concave wall of the new pairs of counter-rotating vortices (Dean vortices), the critical value of the Dean number decreases as the power-law index increases for the power-law fluids, and the Bingham number decreases for the Bingham fluids. For power-law fluids, a diagram of critical Dean numbers, based on the number of Dean vortices formed, was established for different axial positions. For the same flow conditions, the critical Dean number obtained using the axial velocity gradient criterion was smaller then that obtained with the visualization technique.     

Simulation of two generalised Newtonian fluids flow in micropore with dead end

This research deals with the numerical simulation of Carreau and power-law fluids flow in an open capillary of a reservoir. The capillary is connected to a dead end. The finite volume method (FVM) on a structured and co-located grid has been used. The numerical method has been validated through the comparison of numerical results against the analytical solutions of power-law fluid flow in a planar channel. The effects of fluids, the operating conditions and the aspect ratio of dead end at the low Reynolds (Re) numbers on the oil sweeping from the dead end are investigated. The simulation results show that by increasing the power-law exponent in the case of power-law fluids, the swept depth in the dead end increases. However, according to the results, the effect of Re number on the flow pattern and the oil sweeping from the dead end is insignificant at the investigated conditions. In the case of Carreau model, at the conditions investigated, the swept area increases as the power-law exponent increases, but the Reynolds number has still minor effects on the flow pattern. Also, as the aspect ratio of dead end increases, the sweep efficiency increases.     

A new numerical scheme for convective dominated heat transfer in a porous medium with strong temperature gradients

A new numerical scheme, theimplicit correction scheme, has been developed for heat transfer in a porous medium with strong temperature gradients. The scheme includes diffusion, convection and transverse heat transfer processes. By using correction coefficients which are based on transverse heat transfer, the effects of convection enthalpy flow and diffusion are modified. Under suitable limiting conditions, the implicit correction scheme can be reduced to the central-difference, upwind, or power-law scheme. The correction scheme is shown to be especially useful in calculations of the thermal effectiveness of the regenerator in Stirling cycle refrigeration.     

Comparative study of forced convection of a power-law fluid in a channel with a built-in square cylinder

A detailed comparison between the lattice Boltzmann method and the finite element method is presented for an incompressible steady laminar flow and heat transfer of a power-law fluid past a square cylinder between two parallel plates. Computations are performed for three different blockage ratios (ratios of the square side length to the channel width) and different values of the power-law index n covering both pseudo-plastic fluids (n < 1) and dilatant fluids (n > 1). The methodology is validated against the exact solution. The local and averaged Nusselt numbers are also presented. The results show that the relatively simple lattice Boltzmann method is a good alternative to the finite element method for analyzing non-Newtonian fluids.     

Modelling of Power-law Fluid Flow Through Porous Media Using Smoothed Particle Hydrodynamics

The flow of non-Newtonian fluids through two-dimensional porous media is analyzed at the pore scale using the smoothed particle hydrodynamics (SPH) method. A fully explicit projection method is used to simulate incompressible flow. This study focuses on a shear-thinning power-law model (n < 1), though the method is sufficiently general to include other stress-shear rate relationships. The capabilities of the proposed method are demonstrated by analyzing a Poiseuille problem at low Reynolds numbers. Two test cases are also solved to evaluate validity of Darcy’s law for power-law fluids and to investigate the effect of anisotropy at the pore scale. Results show that the proposed algorithm can accurately simulate non-Newtonian fluid flows in porous media.     

Non-Darcy Free Convection of Power-Law Fluids Over a Two-Dimensional Body Embedded in a Porous Medium

A boundary layer analysis was presented to study the non-Darcy-free convection of a power-law fluid over a non-isothermal two-dimensional body embedded in a porous medium. The Ostwald-de Waele power-law model was used to characterize the non-Newtonian fluid behavior. Similarity solutions were obtained with variations in surface temperature or surface heat flux. In view of the fact that most of the non-Newtonian fluids have large Prandtl numbers, this study was directed toward such fluids. The effects of the porous medium parameters, k ₁ and k ₂, body shape parameter, m, and surface thermal variations parameter, p, as well as the power-law index, n, were examined.     

MHD Flow and heat transfer of non-Newtonian power-law fluid with diffusion and chemical reaction on a moving cylinder

The problem of flow and heat transfer of an electrically conducting non-Newtonian fluid over a continuously moving cylinder in the presence of a uniform magnetic field is analyzed for the case of power-law variation in the temperature and concentration at the cylinder surface. A diffusion equation with a chemical reaction source term is taken into account. The governing non-similar partial differential equation are solved numerically by employing shooting method. The effects of various parameters on the velocity, temperature and concentration profiles as well as the heat and mass transfer rate from the cylinder surface to the surrounding fluid are presented graphically and in tabulated form.     

A micromechanical model for polycrystalline creep with grain boundary cavitation

Summary  A micromechanical model is developed to describe effects such as combined power-law creep and diffusion, grain boundary sliding and cavitation in polycrystals. Several aspects of creep-constrained cavitation are taken into account such as diffusion in a cage of creeping matrix material and cavitating facets in a cage of creeping grains. Grain boundary sliding is modelled by distributed micro-shearcracks. It is shown that the different physical mechanisms and their interactions are functions of a well-defined material parameter λ, which can be related to the material length scale L introduced by Rice. Received 18 January 2000; accepted for publication 17 May 2000     

Thermal radiation effect on mixed convection heat and mass transfer of a non-Newtonian fluid over a vertical surface embedded in a porous medium in the presence of thermal diffusion and diffusion-thermo effects

Thermal radiation, thermal diffusion, and diffusion-thermo effects on heat and mass transfer by mixed convection of non-Newtonian power-law fluids over a vertical permeable surface embedded in a saturated porous medium are investigated. The governing equations describing the problem are non-dimensionalized and transformed into a non-similar form. The transformed equations are solved by using the local non-similarity method combined with the shooting technique. The effects of the physical parameters of the problem on the fluid temperature and concentration are illustrated graphically and analyzed. Also, the effects of the pertinent parameters on the local Nusselt number and the local Sherwood number are presented.     

Soret and Dufour Effects on Double-Diffusive Free Convection Over a Vertical Truncated Cone in Porous Media with Variable Wall Heat and Mass Fluxes

This work studies the Soret and Dufour effects on the double-diffusive free convection over a downward-pointing vertical truncated cone with variable wall heat and mass fluxes in fluid-saturated porous media. A coordinate transformation is used to derive the nondimensional boundary-layer governing equations, and the obtained nonsimilar equations are then solved by the cubic spline collocation method. Results for local surface temperature and the local surface concentration are presented as functions of Soret parameters, Dufour parameters, power-law exponents, buoyancy ratios, and Lewis numbers. Results show that increasing the Dufour parameter tends to increase the local surface temperature, while it tends to decrease the local surface concentration. An increase in the Soret number leads to a decrease in the local surface temperature for buoyancy assisting flows, while it leads to an increase in the local surface temperature for buoyancy opposing flows. Increasing the Soret number tends to increase the local surface concentration. Moreover, the local surface temperature and the local surface concentration of the truncated cones with higher power-law exponents are lower than those with lower exponents.     

Steady Herschel–Bulkley fluid flow in three-dimensional expansions

In this paper we study steady flow of Herschel–Bulkley fluids in a canonical three-dimensional expansion. The fluid behavior was modeled using a regularized continuous constitutive relation, and the flow was obtained numerically using a mixed-Galerkin finite element formulation with a Newton–Raphson iteration procedure coupled to an iterative solver. Results for the topology of the yielded and unyielded regions, and recirculation zones as a function of the Reynolds and Bingham numbers and the power-law exponent, are presented and discussed for a 2:1 and a 4:1 expansion ratio. The results reveal the strong interplay between the Bingham and Reynolds numbers and their influence on the formation and break up of stagnant zones in the corner of the expansion and on the size and location of core regions.     

Transient natural convection along a vertical clyinder with variable surface temperature and mass diffusion

A numerical study is carried out for thermal and concentration driven transient natural convection adjacent to a vertical cylinder. The temperature and concentration level at the cylinder surface are assumed to vary as power-law type functions, with exponents n and m respectively in the streamwise co-ordinate. The governing boundary layer equations are converted into a non-dimensional form. A Crank-Nicolson type of implicit finite-difference method is used to solve the governing non-linear set of equations. Numerical results are obtained and presented with various thermal and mass Grashof numbers and power law variations. Transient effects of velocity, temperature and concentration are analyzed. Local and average skin-friction, Nusselt number and Sherwood number are shown graphically.     

Flow of non-newtonian particulate suspension with a compressible particle phase

Continuum equations for a two-phase fluid-particle flow are developed and applied to the problem of steady, laminar flow over an infinite porous flat plate. Both phases are assumed to behave as non-Newtonian power-law fluids. The effects of particle-particle interaction and diffusion of particles are taken into account in the mathematical model. In addition, the particle phase is assumed to have a non-uniform density distribution. The resulting governing equations are nondimensionalized and solved numerically subject to appropriate boundary conditions using an iterative, implicit finite-difference method. Graphical results for the displacement thicknesses and the skin-friction coefficients for both the fluid and particle phases are presented and discussed to elucidate interesting features of the solutions.     

复杂介质中任意阶频率依赖耗散声波的分数阶导数模型

实验现象表明,声波在复杂介质中传播时,其衰减往往呈现频率的任意次幂律依赖现象.鉴于复杂介质的力学和物理性质的记忆性和长程相关性,频率幂律依赖的声波衰减现象难以用经典的声波方程描述,因为经典的阻尼波方程和近似热黏性波方程只能分别描述与频率无关和频率二次方依赖的声衰减.近年来,带有分数阶导数项的声波方程已被成功用于描述这一声衰减现象.基于课题组对声波衰减分数阶导数建模的研究,对已有的分数阶导数声波方程的研究进展及获得的成果做一个系统的综述,重点讨论这些模型的力学本构、统计力学解释等.简述了软物质中声波传播的时间分数阶导数唯象模型和本构模型,空间分数阶导数唯象模型和本构模型,并深入讨论了各种模型之间的联系与区别:介绍了分数阶导数声波模型在多孔介质中的成功应用,该部分内容涉及了均匀和非均匀多孔介质,刚性固体骨架和可变形固体骨架多孔介质等;通过空间分数阶扩散方程与Levy稳定分布之间的联系,给出了频率幂律依赖指数的变化区间为[0,2]的统计力学解释.最后,讨论了声波传播耗散行为的分数阶导数建模领域仍然存在的问题,并对今后的研究方向进行了探讨和展望.