Kinematics of Bulkley–Herschel fluid flow with a free surface during the filling of a channel

The non-Newtonian fluid flow with a free surface occurring during the filling of a plane channel in the gravity field is modeled. The mathematical formulation of the problem using the rheological Bulkley–Herschel model is presented. A numerical finite-difference algorithm for solving this problem is developed. A parametric investigation of the main characteristics of the process as functions of the control parameters is performed. The effect of the rheological parameters of the fluid on the distribution of the quasisolid motion zones is demonstrated.     

Pressure-driven and electroosmotic non-Newtonian flows through microporous media via lattice Boltzmann method

The lattice Boltzmann method is developed to simulate the pressure-driven flow and electroosmotic flow of non-Newtonian fluids in porous media based on the representative elementary volume scale. The flow through porous media was simulated by including the porosity into the equilibrium distribution function and adding a non-Newtonian force term to the evolution equation. The non-Newtonian behavior is considered based on the Herschel–Bulkley model. The velocity results for pressure-driven non-Newtonian flow agree well with the analytical solutions. For the electroosmotic flow, the influences of porosity, solid particle diameter, power law exponent, yield stress and electric parameters are investigated. The results demonstrate that the present lattice Boltzmann model is capable of modeling non-Newtonian flow through porous media.     

Effect of the bottom material capture and the non-Newtonian rheology on the dynamics of turbulent downslope flows

The paper is devoted to the mathematical modeling of naturally occurring downslope flows, such as snow avalanches, mudflows, and rapid landslides. The medium in motion is modeled as a non- Newtonian fluid, the non-Newtonian fluids of different types corresponding to different-in-nature flows. It is taken into account that the downslope flows capture the slope material and entrain it into the motion. The flow is assumed to be turbulent and the Lushchik–Pavel’ev–Yakubenko three-equation turbulence model is used. It is so generalized that it allows for flow unsteadiness, complicated rheological properties, the presence of a free boundary, and the mass transfer at the lower flow boundary. The effect of the bottom material capture and the nonlinear rheological properties of the medium in motion on the flow dynamics is numerically investigated.     

Non-Newtonian fluids with a yield stress

A modified Herschel–Bulkley model [E. Mitsoulis, S.S. Abdali, Flow simulation of Herschel–Bulkley fluids through extrusion dies, Can. J. Chem. Eng. 71 (1993) 147–160] predicts an infinite apparent viscosity at vanishing shear rate. Furthermore, the dimensions of one parameter depend on another parameter. In this contribution, we propose a generalized model based on earlier work by De Kee and Turcotte [D. De Kee, G. Turcotte, Viscosity of biomaterials. Chem. Eng. Commun. 6 (1980) 273–282] and on the work of Papanastasiou [T.C. Papanastasiou, Flows of materials with yield, J. Rheol. 31 (1987) 385–404] to solve the problems associated with the modified Herschel–Bulkley model. Compared to the responses of the Papanastasiou model and the modified Herschel–Bulkley model, the proposed generalized model provides the expected improvements and is capable of predicting successfully the rheological behavior (viscosity and yield stress) of Carbopol 980 dispersions.     

Application of Rheological Models in Prediction of Turbulent Slurry Flow

The paper deals with fully developed steady turbulent flow of slurry in a circular straight and smooth pipe. The Kaolin slurry consists of very fine solid particles, so the solid particles concentration, and density, and viscosity are assumed to be constant across the pipe. The mathematical model is based on the time averaged momentum equation. The problem of closure was solved by the Launder and Sharma k-ε turbulence model (Launder and Sharma, Lett Heat Mass Transf 1:131–138, 1974) but with a different turbulence damping function. The turbulence damping function, used in the mathematical model in the present paper, is that proposed by Bartosik (1997). The mathematical model uses the apparent viscosity concept and the apparent viscosity was calculated using two- and three-parameter rheological models, namely Bingham and Herschel–Bulkley. The main aim of the paper is to compare measurements and predictions of the frictional head loss and velocity distribution, taking into account two- and three-parameter rheological models, namely Bingham and Herschel–Bulkley, if the Kaolin slurry possesses low, moderate, and high yield stress. Predictions compared with measurements show an observable advantage of the Herschel–Bulkley rheological model over the Bingham model particularly if the bulk velocity decreases.     

Modeling of the Installation and Stability of Gel Barriers in Main Fractures

A quasi-one-dimensional nonstationary model for the shutoff of water flow into main fractures using Herschel–Bulkley viscoplastic fluid is presented. The mathematical model, developed in a one-dimensional isothermal approximation with the hydrodynamic parameters (pressure and velocity) averaged over the cross section of the fracture, can be used to determine the optimal technological parameters of the shutoff process and the size of gel barriers to provide their stability on exposure to intense filtration flow after the shutoff of the flow into the fracture. The ranges of flow rates which provide shutoff of a fracture of fixed width for the selected gel are determined. The stability of the installed gel barriers is estimated.     

Three‐dimensional lubrication flow of a Herschel–Bulkley fluid

In this paper three‐dimensional lubrication flow of grease is analysed numerically. The lubrication flow configuration is formed by two ellipsoid rollers. The load is assumed to be light enough for the lubrication mode to be purely hydrodynamic. The fluid behaviour is modelled using the Herschel–Bulkley model, and a two‐dimensional modified Reynolds equation is derived. The numerical solutions are obtained by using a hybrid spectral/iterative technique and the Galerkin projection scheme. The effects of the material and geometrical parameters on pressure distribution are emphasized in the study. The investigation is conducted for a situation where the two ellipsoids are fully immersed in a grease lubricant. The effect of the geometry on the pressure distribution is determined by varying the ratio of the semi‐axes and the minimum gap of the two rollers, respectively. The effect of the material parameters is examined by varying the power‐law index and yield stress. It is found that the pressure distribution is strongly influenced by the shape of the rollers, the size of the minimum gap of the rollers and the rheological parameters. Copyright © 2005 John Wiley & Sons, Ltd.     

The extrusion of a model yield stress fluid imaged by MRI velocimetry

Extrusion tests were performed by forcing a well-characterized model yield stress fluid from a cylindrical cartridge through various cylindrical extrusion dies using a variety of different piston velocities. In this study the Bingham number within the die ranged from 0.1 to 10. MRI techniques allowed for the non-invasive determination of the local velocity within the extruded material in the range [0.015; 20 mm s^?1]. The velocity profile within a very long die was determined by MRI and agreed very well with the analytical results for the flow of a Herschel–Bulkley fluid within a conduit using parameters determined from independent rheometrical tests, validating both the rheological approach and the accuracy of the MRI techniques. Although the velocity was determined by MRI in the upper and lower zones separately, the intersection of these zones showed great agreement, providing an entire view of the extrusion process. In the range of Bingham number studied, the velocity field for a given contraction ratio appeared similar when scaled by the piston velocity, with a dimpled acceleration zone above the die and lateral dead zones varying negligibly with the piston velocity. For a further analysis the experimental results were compared with the results of numerical simulations. Finite element simulations using an elastic solids model were performed to provide this comparison. It was found that this model did well in representing the characteristics of extrusion flow seen in the experiments; an aspect that was not present in the biviscous simulations. The MRI results show that for the range of values studied, both the piston velocity and the contraction ratio have little effect on the characteristics of the flow, including the size and location of the apparent dead zones. It was found that with an appropriate scaling the central, longitudinal velocity follows a master curve. A decreasing contraction ratio, on the other hand, appears to increase the size of the weak velocity region, in contrast with the simulation results.     

Rheological characterization and constitutive modeling of bread dough

Bread dough (a flour–water system) has been rheologically characterized using a parallel-plate, an extensional, and a capillary rheometer at room temperature. Based on the linear and nonlinear viscoelastic and viscoplastic data, two constitutive equations have been applied, namely a viscoplastic Herschel–Bulkley model and a viscoelastoplastic K–BKZ model with a yield stress. For cases where time effects are unimportant, the viscoplastic Herschel–Bulkley model can be used. For cases where transient effects are important, it is more appropriate to use the K-BKZ model with the addition of a yield stress. Finally, the wall slip behavior of dough was studied in capillary flow, and an appropriate slip law was formulated. These models characterize the rheological behavior of bread dough and constitute the basic ingredients for flow simulation of dough processing, such as extrusion, calendering, and rolling.     

Analysis of the impact of a sharp indenter

The present work investigates the impact of a sharp indenter at low impact velocities. A one-dimensional model is developed by assuming that the variation of indentation load as a function of depth under dynamic conditions has the same parabolic form (Kick's Law) as under static conditions. The motion of the indenter as it indents and rebounds from the target is described. Predictions are made of the peak indentation depth, residual indentation depth, contact time, and rebound velocity as functions of the impact velocity, indenter mass and target properties. Finite element simulations were carried out to assess the validity of the model for elastoplastic materials. For rate-independent materials agreement with the model was good provided the impact velocity did not exceed certain critical values. For rate-dependent materials the relationship between load and depth in the impact problem is no longer parabolic and the model predictions cannot be applied to this case. The rate-dependent case can be solved by incorporating the relationship between the motion of the indenter and the dynamic flow properties of the material into the equation of motion for the indenter.     

Nonlinear Dynamics of Turbulent Thermals in Shear Flow

The nonlinear integral model of a turbulent thermal is extended to the case of the horizontal component of its motion relative to the medium (e.g., thermal floating-up in shear flow). In contrast to traditional models, the possibility of a heat source in the thermal is taken into account. For a piecewise constant vertical profile of the horizontal velocity of the medium and a constant vertical velocity shear, analytical solutions are obtained which describe different modes of dynamics of thermals. The nonlinear interaction between the horizontal and vertical components of thermal motion is studied because each of the components influences the rate of entrainment of the surrounding medium, i.e., the growth rate of the thermal size and, hence, its mobility. It is shown that the enhancement of the entrainment of the medium due to the interaction between the thermal and the cross flow can lead to a significant decrease in the mobility of the thermal.     

An augmented Lagrangian approach to simulating yield stress fluid flows around a spherical gas bubble

We simulated the flow of a yield stress fluid around a gas bubble using an augmented Lagrange approach. The piecewise linear equal‐order finite elements for both the velocity and the pressure approximations proposed and analyzed by Latché and Vola in 2004 were applied. A mesh adaptive strategy based on this element‐pair choice was also proposed to render the yield surfaces of desired resolution. The corresponding numerical scheme was formulated for general Herschel–Bulkley fluids. Numerical results on Bingham fluid flows around a slowly rising spherical gas bubble were provided to validate the proposed algorithm. Copyright © 2011 John Wiley & Sons, Ltd.     

波浪作用下直立结构物附近强湍动掺气流体运动的数值模拟

波浪破碎卷入气体易对建筑物受力产生压力振荡, 了解波浪作用下建筑物附近掺气水流的运动特性是精确计算建筑物受力的前提. 基于OpenFOAM开源程序包和修正速度入口造波方法建立三维数值波浪水槽, 模型采用S-A IDDES湍流模型进行湍流封闭, 并采用修正的VOF 方法捕捉自由液面, 数值模拟了规则波在1:10的光滑斜坡上与直立结构物的相互作用过程, 重点分析了结构物附近的水动力和掺气水流运动特性. 结果表明, 建立的数值模型能精确地捕捉波浪作用下直立结构物附近的自由液面的变化以及气泡输运过程, 较好地描述气体卷入所形成的气腔形态以及多气腔之间的融合、分裂等过程; 波浪与直立结构物相互作用产生强湍动掺气水流, 其运动过程十分复杂; 掺气流体输运过程中水气界面周围一直伴随着涡的存在, 其中, 气泡的分裂与周围正负涡量剪切作用密切相关, 且其输运轨迹主要受周围流场的影响; 研究揭示了结构物附近湍动能与掺气特性的关系, 发现波浪作用下直立结构物附近湍动能的分布与掺气水流特征参数(气泡数量、空隙率)整体呈现一定的线性关系.      

Using the Euler–Lagrange variational principle to obtain flow relations for generalized Newtonian fluids

The Euler–Lagrange variational principle is used to obtain analytical and numerical flow relations in cylindrical tubes. The method is based on minimizing the total stress in theflow duct using the fluid constitutive relation between stress and rate of strain. Newtonian and non-Newtonian fluid models, which include power law, Bingham, Herschel–Bulkley, Carreau, and Cross, are used for demonstration.     

Dynamic coupling of fluid and structural mechanics for simulating particle motion and interaction in high speed compressible gas particle laden flow

In this work, structural finite element analyses of particles moving and interacting within high speed compressible flow are directly coupled to computational fluid dynamics and heat transfer analyses to provide more detailed and improved simulations of particle laden flow under these operating conditions. For a given solid material model, stresses and displacements throughout the solid body are determined with the particle–particle contact following an element to element local spring force model and local fluid induced forces directly calculated from the finite volume flow solution. Plasticity and particle deformation common in such a flow regime can be incorporated in a more rigorous manner than typical discrete element models where structural conditions are not directly modeled. Using the developed techniques, simulations of normal collisions between two 1 mm radius particles with initial particle velocities of 50–150 m/s are conducted with different levels of pressure driven gas flow moving normal to the initial particle motion for elastic and elastic–plastic with strain hardening based solid material models. In this manner, the relationships between the collision velocity, the material behavior models, and the fluid flow and the particle motion and deformation can be investigated. The elastic–plastic material behavior results in post collision velocities 16–50% of their pre-collision values while the elastic-based particle collisions nearly regained their initial velocity upon rebound. The elastic–plastic material models produce contact forces less than half of those for elastic collisions, longer contact times, and greater particle deformation. Fluid flow forces affect the particle motion even at high collision speeds regardless of the solid material behavior model. With the elastic models, the collision force varied little with the strength of the gas flow driver. For the elastic–plastic models, the larger particle deformation and the resulting increasingly asymmetric loading lead to growing differences in the collision force magnitudes and directions as the gas flow strength increased. The coupled finite volume flow and finite element structural analyses provide a capability to capture the interdependencies between the interaction of the particles, the particle deformation, the fluid flow and the particle motion.     

Rheology of pulp suspensions using ultrasonic Doppler velocimetry

Conventional rheometry coupled with local velocity measurements (ultrasonic Doppler velocimetry) are used to study the flow behaviour of various commercial pulp fibre suspensions at fibre mass concentrations ranging from 1 to 5 wt.%. Experimental data obtained using a stress-controlled rheometer by implementing a vane in large cup geometry exhibits apparent yield stress values which are lower than those predicted before mainly due to existence of apparent slip. Pulp suspensions exhibit shear-thinning behaviour up to a high shear rate value after which Newtonian behaviour prevails. Local velocity measurements prove the existence of significant wall slippage at the vane surface. The velocimetry technique is also used to study the influence of pH and lignin content on the flow behaviour of pulp suspensions. The Herschel–Bulkley constitutive equation is used to fit the local steady-state velocity profiles and to predict the steady-state flow curves obtained by conventional rheometry. Consistency between the various sets of data is found for all suspensions studied, including apparent yield stress, apparent wall slip and complete flow curves.     

Permanent Waves in Slow Free-Surface Flow of a Herschel–Bulkley fluid

Unsteady flow of a viscoplastic fluid on an inclined plane is examined. The fluid is described by the three-parameter Herschel–Bulkley constitutive equation. The set of equations governing the flow is presented, recovering earlier results for a Bingham fluid and steady uniform motion. A permanent wave solution is then derived, and the relation between wave speed and flow depth is discussed. It is shown that more types of gravity currents are possible than in a Newtonian fluid; these include some cases of flows propagating up a slope. The speed of permanent waves is derived and the possible surface profiles are illustrated as functions of the flow behavior index.     

Multi-scale defect interactions in high-rate failure of brittle materials,Part II: Application to design of protection materials

Micromechanics based damage models, such as the model presented in Part I of this 2 part series (Tonge and Ramesh, 2015), have the potential to suggest promising directions for materials design. However, to reach their full potential these models must demonstrate that they capture the relevant physical processes. In this work, we apply the multiscale material model described in Tonge and Ramesh (2015) to ballistic impacts on the advanced ceramic boron carbide and suggest possible directions for improving the performance of boron carbide under impact conditions. We simulate both dynamic uniaxial compression and simplified ballistic loading geometries to demonstrate that the material model captures the relevant physics in these problems and to interrogate the sensitivity of the simulation results to some of the model input parameters. Under dynamic compression, we show that the simulated peak strength is sensitive to the maximum crack growth velocity and the flaw distribution, while the stress collapse portion of the test is partially influenced by the granular flow behavior of the fully damaged material. From simulations of simplified ballistic impact, we suggest that the total amount of granular flow (a possible performance metric) can be reduced by either a larger granular flow slope (more angular fragments) or a larger granular flow timescale (larger fragments). We then discuss the implications for materials design.     

Experimental investigation of turbulent entrainment in an annulus with moving sidewalls

The establishment of a turbulent mixed layer in a two-layer stratified shear flow, and the rate of entrainment into that layer were studied experimentally in a modified annulus. The modification of the conventional annulus was made by replacing the upper rotating screen with inner rotating sidewalls, extending over the upper half of the channel, so that the flow in the upper layer was nearly uniform and almost laminar, while the bottom layer was quiescent. Vertical density profile measurements were conducted using single electrode conductivity probes. The flow was visualized during the various stages of the experiment using the hydrogen bubble technique.After the start of the sidewalls rotation, the upper layer accelerates from rest, and consequently a transition process is taking place during which the initial density interface between the two layers is developed into a turbulent mixed layer. This turbulent layer is bounded by two sharp interfaces, each separating it from an outer non-turbulent zone. The generation of this five-layer structure seemed to be dominated by instabilities activated by the velocity difference between the upper and lower layer.Once a turbulent mixed layer is formed, entrainment of nonturbulent fluid into that layer is taking place causing its thickness to increase continuously. Depending on the overall Richardson number, based on the channel width, the slope of the entrainment law curve was found to have two different values, each indicating the dominance of a different source of turbulent energy production. For relatively low Richardson numbers, the slope is close to -1.8, implying that the velocity shear across each interface contributes significantly to the entrainment. On the other hand, for larger Richardson numbers the slope is about -1.25, in agreement with previous results of shear-free entrainment experiments.The measured velocity profiles indicate that as long as the mixed layer is not too thick, the radial inhomogeneities are small and the flow may be considered as nearly one-dimensional. It seems, therefore, that for the understanding of entrainment processes occurring in realistic stratified flows, the modified annulus is a more reliable tool than the conventional one.     

Criteria for the appearance of recirculating and non-stationary regimes behind a cylinder in a viscoplastic fluid

The appearance of a recirculation zone and the formation of non-stationary vortices behind a cylinder in the unconfined flow of a Herschel–Bulkley fluid have been studied by numerical simulation. The Herschel–Bulkley constitutive equation was regularised by using the Papanastasiou model. Special attention was paid to determining the numerical parameters and comparing them to existing results. The influence of the Oldroyd number and power-law index on flow morphology and, in particular, on the unyielded zones was studied over a wide spectrum (0 ≤ Od ≤ 10) and (0.3 ≤ n ≤ 1.8). It was seen that the greater the Oldroyd number, the greater the critical Reynolds numbers and Strouhal number for the two flow regimes. The influence of the power-law index is more complex.