GPU accelerated simulations of three-dimensional flow of power-law fluids in a driven cube

Newtonian fluid flow in two- and three-dimensional cavities with a moving wall has been studied extensively in a number of previous works. However, relatively a fewer number of studies have considered the motion of non-Newtonian fluids such as shear thinning and shear thickening power law fluids. In this paper, we have simulated the three-dimensional, non-Newtonian flow of a power law fluid in a cubic cavity driven by shear from the top wall. We have used an in-house developed fractional step code, implemented on a Graphics Processor Unit. Three Reynolds numbers have been studied with power law index set to 0.5, 1.0 and 1.5. The flow patterns, viscosity distributions and velocity profiles are presented for Reynolds numbers of 100, 400 and 1000. All three Reynolds numbers are found to yield steady state flows. Tabulated values of velocity are given for the nine cases studied, including the Newtonian cases.     

Anomalous predictions of pressure drop and heat transfer in ducts of arbitrary cross-section with modified power law fluids

 The apparent viscosities of purely viscous non-Newtonian fluids are shear rate dependent. At low shear rates, many of such fluids exhibit Newtonian behaviour while at higher shear rates non-Newtonian, power law characteristics exist. Between these two ranges, the fluid's viscous properties are neither Newtonian or power law. Utilizing an apparent viscosity constitutive equation called the “Modified Power Law” which accounts for the above behavior, solutions have been obtained for forced convection flows. A shear rate similarity parameter is identified which specifies both the shear rate range for a given fluid and set of operating conditions and the appropriate solution for that range. The results of numerical solutions for the friction factor–Reynolds number product and for the Nusselt number as a function of a dimensionless shear rate parameter have been presented for forced fully developed laminer duct flows of different cross-sections with modified power law fluids. Experimental data is also presented showing the suitability of the “Modified Power Law” constitutive equation to represent the apparent viscosity of various polymer solutions. Received on 21 August 2000     

Displacement of a Newtonian fluid by a non-Newtonian fluid in a porous medium

This paper presents an analytical Buckley-Leverett-type solution for one-dimensibnal immiscible displacement of a Newtonian fluid by a non-Newtonian fluid in porous media. The non-Newtonian fluid viscosity is assumed to be a function of the flow potential gradient and the non-Newtonian phase saturation. To apply this method to field problems a practical procedure has been developed which is based on the analytical solution and is similar to the graphic technique of Welge. Our solution can be regarded as an extension of the Buckley-Leverett method to Non-Newtonian fluids. The analytical result reveals how the saturation profile and the displacement efficiency are controlled not only by the relative permeabilities, as in the Buckley-Leverett solution, but also by the inherent complexities of the non-Newtonian fluid. Two examples of the application of the solution are given. One application is the verification of a numerical model, which has been developed for simulation of flow of immiscible non-Newtonian and Newtonian fluids in porous media. Excellent agreement between the numerical and analytical results has been obtained using a power-law non-Newtonian fluid. Another application is to examine the effects of non-Newtonian behavior on immiscible displacement of a Newtonian fluid by a power-law non-Newtonian fluid.     

A unified compatibility method for exact solutions of non-linear flow models of Newtonian and non-Newtonian fluids

Criteria are established for higher order ordinary differential equations to be compatible with lower order ordinary differential equations. Necessary and sufficient compatibility conditions are derived which can be used to construct exact solutions of higher order ordinary differential equations subject to lower order equations. We provide the connection to generalized groups through conditional symmetries. Using this approach of compatibility and generalized groups, new exact solutions of non-linear flow problems arising in the study of Newtonian and non-Newtonian fluids are derived. The ansatz approach for obtaining exact solutions for non-linear flow models of Newtonian and non-Newtonian fluids is unified with the application of the compatibility and generalized group criteria.     

Parallel domain decomposition method for finite element approximation of 3D steady state non‐Newtonian fluids

We introduce a stabilized finite element method for the 3D non‐Newtonian Navier–Stokes equations and a parallel domain decomposition method for solving the sparse system of nonlinear equations arising from the discretization. Non‐Newtonian flow problems are, generally speaking, more challenging than Newtonian flows because the nonlinearities are not only in the convection term but also in the viscosity term, which depends on the shear rate. Many good iterative methods and preconditioning techniques that work well for the Newtonian flows do not work well for the non‐Newtonian flows. We employ a Galerkin/least squares finite element method, with stabilization parameters adjusted to count the non‐Newtonian effect, to discretize the equations, and the resulting highly nonlinear system of equations is solved by a Newton–Krylov–Schwarz algorithm. In this study, we apply the proposed method to some inelastic power‐law fluid flows through the eccentric annuli with inner cylinder rotation and investigate the robustness of the method with respect to some physical parameters, including the power‐law index and the Reynolds number ratios. We then report the superlinear speedup achieved by the domain decomposition algorithm on a computer with up to 512 processors. Copyright © 2015 John Wiley & Sons, Ltd.     

The turbulent flow of non-Newtonian fluids in the absence of anomalous wall effects

For Newtonian fluids, the engineering predictions for pressure drop in turbulent pipe flow are well established. However, in the case of non-Newtonian liquids, only a few design techniques have been proposed and these do not share a common basis with the approach for Newtonian systems. This present work attempts to provide a common basis for both Newtonian and non-Newtonian systems in situations where anomalous wall effects are absent. Previously published experimental data suggest that if the Reynolds number is calculated on the basis of the apparent viscosity at the wall then the standard Newtonian correlations can be used for the prediction of pressure drop. The use of the wall viscosity in defining the Reynolds number also serves as a test for anomalous behaviour. Any departure of the experimental data from the Newtonian turbulent friction factor correlation indicates anomalous behaviour.     

Separating forces in blade coating of viscous and viscoelastic liquids

Coating of viscous and viscoelastic liquids is examined both theoretically and experimentally. A single simple geometry, a blade over a rotating roll, is considered. A perturbation solution to the Navier-Stokes equations yields a lubrication theory with first order corrections for curvature and inertia. A numerical solutions by the Finite Element Method (FEM) is compared to the analytical solutions. For Newtonian fluids, agreement between these mathematical models, and data on blade loading, is quite good.The effect of a non-Newtonian viscosity is explored by adopting a purely viscous power law model. The zeroth-order (lubrication) equations are solved by the method of Steidler and Horowitz, and predictions for coating thickness and blade loading agree quite well with those obtained from a FEM solution of the full equations of motion for a power law fluid. Data on blade loading, obtained using a strongly elastic polymer solution, are compared to these mathematical models, and discrepancies are noted.     

Pulsatile flows of Leslie–Ericksen liquid crystals

Capillary pulsatile flows of calamitic (rod-like) and discotic nematic liquid crystals are analyzed using the Leslie–Ericksen equations for low-molar mass liquid crystals, using computational, analytical, and scaling methods. The dependence of flow-enhancement and power requirement on frequency, amplitude, pressure drop wave-form, molecular geometry is characterized. The unique roles of orientation-dependent local viscosity and backflow (orientation-driven flow) on flow-enhancement and power requirement are elucidated. The local viscosity effect is shown to be a significant factor in flow-enhancement at all pressure drops, but only affects power requirement at higher pressure drops. Backflow has weak effects on flow-enhancement and large effects on power requirements at low average pressure drops. Amplitude, frequency, and molecular geometry effects are clearly manifested through viscosity and backflow. A detailed comparison with predictions for power law fluids shows a clear correspondence between these non-Newtonian fluids and nematic liquid crystals. The unique distinguishing feature of pulsatile flows of liquid crystals is found to be backflow, such that power increases with increasing frequency, a featured that does not exist in other non-Newtonian fluids due to lack of a strong flow driven by restructuring/re-orientation processes. Future use of these new results may include measurements of viscoelastic parameters that control backflow.     

The Effect of Double Dispersion on Natural Convection Heat and Mass Transfer in a Non-Newtonian Fluid Saturated Non-Darcy Porous Medium

The effect of power law index parameter of the non-Newtonian fluid on free convection heat and mass transfer from a vertical wall is analyzed by considering double dispersion in a non-Darcy porous medium with constant wall temperature and concentration conditions. The Ostwald–de Waele power law model is used to characterize the non-Newtonian fluid behavior. In this case a similarity solution is possible. The variation of heat and mass transfer coefficients with the governing parameters such as power law index, thermal and solutal dispersion parameters, inertia parameter, buoyancy ratio, and the Lewis number is discussed for a wide range of values of these parameters.     

基于舌/上颚微间隙下流体流动行为研究

为探讨口腔环境下流体的流动行为，采用数值方法与流变试验深入研究舌/上颚微间隙下流体流量的影响因素. 建立舌/上颚微间隙的简化模型及Reynolds方程，通过数值方法获取微间隙下流量变化；在DHR-2流变仪上研究非牛顿流体的黏度与剪切率的变化，探讨牛顿流体和非牛顿流体的流量影响. 结果表明:牛顿流体流量平方的倒数同载荷和黏度比值和时间均呈线性函数关系；所制备的非牛顿流体近似为幂律流体，其黏度随脂肪含量的增加而增大，而非牛顿流体流量率先高于后低于等效牛顿流体，其研究结果将为特定人群功能产品的研发提供技术支持.      

Jeffery-Hamel flow of non-Newtonian fluid with nonlinear viscosity and wall friction

A Jeffery-Hamel(J-H) flow model of the non-Newtonian fluid type inside a convergent wedge(inclined walls) with a wall friction is derived by a nonlinear ordinary differential equation with appropriate boundary conditions based on similarity relationships. Unlike the usual power law model, this paper develops nonlinear viscosity based only on a tangential coordinate function due to the radial geometry shape. Two kinds of solutions are developed, i.e., analytical and semi-analytical(numerical) solutions with suitable assumptions. As a result of the parametric examination, it has been found that the Newtonian normalized velocity gradually decreases with the tangential direction progress. Also, an increase in the friction coefficient leads to a decrease in the normalized Newtonian velocity profile values. However, an increase in the Reynolds number causes an increase in the normalized velocity function values. Additionally, for the small values of wedge semi-angle, the present solutions are in good agreement with the previous results in the literature.     

Computational Modelling of Non-Newtonian Effects on Flow in Channels with Moving Wall Indentations

Non-Newtonian effects in a channel with moving wall indentations are assessed numerically by a finite volume method for solving the unsteady incompressible Navier-Stokes equations and using a power-law model exhibiting shear thinning viscosity and Casson's model as the constitutive equations for the non-Newtonian fluid. The computations show that for a non-Newtonian fluid, there are differences in the velocity profiles and in the structure and size of the reversed flow regions as compared with the corresponding Newtonian fluid. The comparison of non-Newtonian and Newtonian wall shear stress reveals a slight decrease in the magnitude on the average for the non-Newtonian case, eventually resulting in the strength of the “wave train” being slightly weaker than those corresponding to a Newtonian fluid.     

Numerical simulation of the non-Newtonian blood flow through a mechanical aortic valve

This work focuses on the comparison between Newtonian and non-Newtonian blood flows through a bileaflet mechanical heart valve in the aortic root. The blood, in fact, is a concentrated suspension of cells, mainly red blood cells, in a Newtonian matrix, the plasma, and consequently its overall behavior is that of a non-Newtonian fluid owing to the action of the cells’ membrane on the fluid part. The common practice, however, assumes the blood in large vessels as a Newtonian fluid since the shear rate is generally high and the effective viscosity becomes independent of the former. In this paper, we show that this is not always the case even in the aorta, the largest artery of the systemic circulation, owing to the pulsatile and transitional nature of the flow. Unexpectedly, for most of the pulsating cycle and in a large part of the fluid volume, the shear rate is smaller than the threshold level for the blood to display a constant effective viscosity and its shear thinning character might affect the system dynamics. A direct inspection of the various flow features has shown that the valve dynamics, the transvalvular pressure drop and the large-scale features of the flow are very similar for the Newtonian and non-Newtonian fluid models. On the other hand, the mechanical damage of the red blood cells (hemolysis), induced by the altered stress values in the flow, is larger for the non-Newtonian fluid model than for the Newtonian one.     

Resistance to uniaxial extensional flow of fibre suspensions

The tensile stress due to resistance to uniaxial extensional flow of fibre suspensions in Newtonian and non-Newtonian fluids has been measured using the filament stretching technique. It has been found that addition of fibres to a Newtonian fluid increases the extensional viscosity. The steady state results agree with Bachelors theory and the stress growth behaviour is qualitatively predicted by the theory of Dinh and Armstrong. Experimental results from this work have also shown that the behaviour of a fibre suspension in viscoelastic fluid is qualitatively described by Fans equation. The added fibres increase the extensional stress growth coefficient of the viscoelastic fluid at low strain but have marginal effect on the fluid after the onset of strain-hardening.     

The stratified flow of gas and non-newtonian liquid in horizontal pipes

Predictions of pressure drop and holdup are presented for the stratified flow of gas and non-Newtonian liquid obeying the Ostwald-de Waele power law model. The model of Taitel & Dukler (1976) for gas/Newtonian liquid flow is extended to liquids possessing either shear-thinning or shear-thickening laminar flow behaviour and computed results are given for flow behaviour indices in the range $\text{0.1 ≤ n ≤ 2}$. In particular, conditions are defined for drag reduction of the liquid flow by the presence of the gas. It is concluded that drag reduction occurs over the largest ranges of liquid and gas flow rates at the lowest $\text{n}$ values, provided that liquid flow remains laminar, but that maximum drag reduction may be expected for shear-thickening liquids with $\text{n}$ values of 2 or greater. Ratios of the liquid flow rate in the presence of gas to that for liquid flow alone under a constant pressure gradient are also presented. These ratios frequently exceed unity and are greatest for highly shear-thinning liquids.Although the Taitel & Dukler approach is consistent with experiments on gas/Newtonian liquid flow, and, in addition, appears to be valid for immiscible Newtonian liquid-liquid systems, provided that the viscosity ratio of the two phases is at least five, experiments are required to confirm its applicability for gas/non-Newtonian systems.     

Viscosity measurements with a magnetoviscometer in the zero-shear and transition region of polypropylenes

Streaming of a non-Newtonian fluid around a sphere is of special importance not only for measuring viscosities with falling spheres, but also for many problems connected with polymer processing. Using the mentioned measuring principle, attention has to be paid to the following points: The sphere is moving in a fluid (melt) of finite extension which requires the application of wall and perhaps end corrections. These are possibly not the same for Newtonian and non-Newtonian fluids. To calculate the viscosity with the help of Stokes law the steady-state velocity is necessary, and it is essential, how long it takes the sphere to reach it. To compare our results with other data, an average shear rate has to be calculated, since there is no uniform shear rate around the sphere. Velocities being very low in our experiments result in very small Reynolds numbers (Re < 10^–3), which allows the application of Stokes law practically without corrections.The experiments were performed at zero shear and in the transition region above. It turned out, that it is usually not possible to extrapolate from shear-dependent viscosity data to zero-shear viscosity.Dedicated to Prof. A. Neckel on the occasion of his 60th birthday     

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.     

Normal stress effects in viscoelastic fluid lubrication

Non-Newtonian flow effects are evaluated in a slider-bearing configuration. The material model taken is that of the Coleman—Noll second-order fluid. An explicit result is given for the portion of the bearing load supported by the non-Newtonian normal stresses as well as that portion resulting from the usual lubrication theory (Newtonian effect). Particular attention is given to the non-Newtonian effect of a high-polymer additive applied to a Newtonian base stock. The non-Newtonian effect has a particular dependence on the bearing geometry as well as a dependence on the relaxation time of the addtive and the amount by which the additive increases the viscosity. The strength of the non-Newtonian effect is assessed in realistic conditions of bearing operation. We find that under certain conditions the non-Newtonian effect could provide a significant load-supporting capability. However, with slight changes in the conditions of the bearing operation, the non-Newtonian load support is negligible. These results are interpreted and qualified with respect to the limitations of the second-order theory, which does not include shear thinning effects.     

Non-Newtonian viscous oscillating free surface jets, and a new strain-rate dependent viscosity form for flows experiencing low strain rates

A model for oscillating free surface jet flow of a fluid from an elliptical orifice, together with experimental measurements, can be exploited to characterize the elongational viscosity of non-Newtonian inelastic fluids. The oscillating jet flow is predominantly elongational, with a small strain that oscillates rapidly between large and zero strain rates. We find that to reproduce the experimentally observed steady oscillating jet flow in model simulations, the assumed form of the non-Newtonian viscosity as a function of strain rate must have zero gradient, i.e., be Newtonian, at zero strain rate (a behavior exhibited, in general, by real inelastic fluids). We demonstrate that the Cross, Carreau, Prandtl-Eyring, and Powell-Eyring forms, although they have finite viscosity at zero strain rate, have either nonzero or even unbounded gradient at zero, and hence are unable to model oscillating jet behavior. We propose a new non-Newtonian viscous form which has all of the desirable features of existing forms (high and low strain rate plateaus, with adjustable location and steepness of the transition) and the additional feature of Newtonian behavior at low strain rates. Received: 7 February 2000 Accepted: 31 October 2000     

The Falkner–Skan Wedge Flows of Power-Law Fluids Embedded in a Porous Medium

The non-Darcy flow characteristics of power-law non-Newtonian fluids past a wedge embedded in a porous medium have been studied. The governing equations are converted to a system of first-order ordinary differential equations by means of a local similarity transformation and have been solved numerically, for a number of parameter combinations of wedge angle parameter m, power-law index of the non-Newtonian fluids n, first-order resistance A and second-order resistance B, using a fourth-order Runge–Kutta integration scheme with the Newton–Raphson shooting method. Velocity and shear stress at the body surface are presented for a range of the above parameters. These results are also compared with the corresponding flow problems for a Newtonian fluid. Numerical results show that for the case of the constant wedge angle and material parameter A, the local skin friction coefficient is lower for a dilatant fluid as compared with the pseudo-plastic or Newtonian fluids.