Flow of a quasi-Newtonian fluid through a planar contraction

Flows through abrupt contractions are dominated by the rapid extension experienced in passing through the contraction. Thus, it is useful to employ a fluid model which considers the extensional viscosity explicitly in its constitutive equation. In this paper, the quasi-Newtonian fluid model, which admits shear thinning and extension thickening of the viscosity depending on the local type of flow as proposed by Schunk and Scriven [P. Schunk, L. Scriven, J, Rheol 34 (1990) 1085], is applied to the numerical simulation of the flow of a dilute polyacrylamide solution through a planar 4 : 1 contraction. In this theory the extra stress tensor does not only depend on the rate of strain tensor but also on the relative rate of rotation of the fluid. The material function – the viscosity function – is allowed to depend on the invariants of these two kinematic tensors yielding a local distinction between extensional, shear or rotation dominated flow. The governing equations are discretized using a finite volume method. Different model parameters are varied and the simulation results are compared with the generalized Newtonian fluid and experimental data.     

Pore-scale modeling of viscoelastic flow in porous media using a Bautista–Manero fluid

This article examines the extensional flow and viscosity and the converging–diverging geometry as the basis of the peculiar viscoelastic behavior in porous media. The modified Bautista–Manero model, which successfully describes elasticity, thixotropic time dependency and shear-thinning, was used for modeling the flow of viscoelastic materials which also show thixotropic attributes. An algorithm, originally proposed by Philippe Tardy, that employs this model to simulate steady-state time-dependent flow was implemented in a non-Newtonian flow simulation code using pore-scale modeling. The simulation results using two topologically-complex networks confirmed the importance of the extensional flow and converging–diverging geometry on the behavior of non-Newtonian fluids in porous media. The analysis also identified a number of correct trends (qualitative and quantitative) and revealed the effect of various fluid and flow parameters on the flow process. The impact of some numerical parameters was also assessed and verified.     

Numerical study of the rheology of rigid fillers suspended in long-chain branched polymer under planar extensional flow

We report a detailed numerical study of the rheology of two-dimensional rigid fillers suspended in branched polymer melt under planar extensional flow. The polymer melt is modelled using the pom-pom constitutive equation. The numerical method uses a finite element solution of the flow in a unit cell within the self-replicating lattice for planar extensional flow identified by Kraynik and Reinelt [Int. J. Multiphase flow 18 (1992) 1045]. This is implemented using a quotient space representation that maps all points in space onto points within the unit cell. We show that the Kraynik and Reinelt cell allows simulations of suspensions under planar extensional flow to be conducted to large strains in a truly periodic cell. The influence of both the pom-pom parameters and the particle area fraction on the rheology of the suspension are investigated. We find a reduction in the degree of extension-rate thickening with the increase of both particles concentration and Weissenberg numbers in agreement with published experimental and numerical findings on other viscoelastic models. This reduction is found to be due to flow disturbance created by the particles which disrupts the alignment of backbone tube segments with the extensional axis.     

Generation of sub-micron-sized droplets,by continuous extensional flow

A method is devised to generate sub-micron-sized droplets without surfactants by using a continuous extensional flow. By using the extensional flow judiciously the physical limits on the viscosity ratio of the dispersed and dispersing media could be overcome. The strain hardening behavior of the Boger fluid under the extensional flow gave a synergic effect to generate the smaller droplets. It was shown that the drop size could be decreased to a sub-micron scale and the droplet size distribution was much narrower compared to the case of the conventional shear (stirring) method. This method is also applicable when the two phases are compatible.     

Extensional behavior influence on viscoelastic turbulent channel flow

In this work we use in the simulation of a viscoelastic turbulent channel flow a modification of the finitely extensible of non-linear elastic dumbbells with the Peterlin approximation (FENE-P) constitutive model for dilute polymer solutions, applicable to high extensional deformations. The new feature introduced by this modification is that the free energy of the polymer (since it is assumed to be entirely entropically driven) remains always bounded (FENE-PB). The characteristics of the model under steady shear flow, pure elongational flow and transient extensional behavior are presented. It is found that the FENE-PB model is more shear thinning than FENE-P. Most importantly, it also shows a higher extensional viscosity than the FENE-P model. Although the steady-state Trouton ratio asymptotically reaches at high extensional rates the same limit as the FENE-P model, the transition from the Newtonian value is sharper and faster. We use the FENE-PB model in direct numerical simulations (DNS) of viscoelastic turbulent channel flow using spectral approximations. The results for various statistics of the flow and the polymer conformation, when compared against those obtained with the original FENE-P model and the same rheological parameters, show an enhanced polymer-induced drag reduction effect and enhanced deformation of the polymer molecules. This indicates that it is not only the asymptotic but also details from the extensional rheological behavior that matter in quantitatively specifying turbulent viscoelastic flow behavior.     

Viscous flow through microfabricated hyperbolic contractions

We study the flow of a Newtonian fluid through microfabricated hyperbolic contractions followed by a sudden expansion, with the aim of investigating the potential of this geometry to serve as an extensional microrheometer. A set of planar converging geometries, with total Hencky strains ranging from 1.0 to 3.7, were fabricated in order to produce a homogeneous extensional flow field within the contraction. The velocity field in various planes of the hyperbolic contraction was quantified by means of microparticle image velocimetry (μPIV) and the pressure drop across the converging geometry was also measured and found to vary approximately linearly with the flow rate. Additionally, an extensive range of numerical calculations were carried out using a finite-volume method to help assess the performance of this geometry as a microfluidic elongational rheometer. The measured velocity fields in the contraction and associated pressure drops compare very well (to within 10%) with the numerical predictions. For the typical dimensions used in the microfluidic devices, the steady viscous flow through the contraction is shown to be three-dimensional and it is demonstrated that regions with nearly constant strain rate can only be achieved using geometries with large total Hencky strains under Hele–Shaw (potential-like) flow conditions.     

Squeeze flow of Bingham fluids under slip with friction boundary condition

This paper develops a theoretical analysis of a Bingham fluid in slipping squeeze flow. The flow field decomposition consists in combining a central extensional flow zone in the plane of symmetry and shear flow zones near the plates. It is also considered that the slipping zone is located around a central sticking zone as previously shown from experiments. It is assumed that the shear stress at the plates is constant in the slipping zone and equals a fixed friction yield value. The squeeze force required to compress a Bingham fluid under the slipping behaviour as well as the radial evolution of the transition point between both sticking and slipping zones are finally determined.     

Computational analysis of techniques to determine extensional viscosity from entrance flows

Recent computational analysis of entrance flows (Mitsoulis et al. 1998) suggests that the entrance pressure drop is insensitive to large changes in steady extensional viscosity-a result that directly contradicts a large body of experimental work in this area. A re-examination of entrance flows using numerical simulations is presented in this work which shows that entrance pressure drops do depend on the steady extensional viscosity, provided the extension rate in the entrance flow is large enough. Numerical simulations are presented using both the strain thinning and thickening versions of the Phan-Thien–Tanner (PTT) constitutive model. Several techniques for extracting extensional viscosity from entrance pressure are applied to the results of these simulations. The resulting predictions of extensional viscosity are compared to the steady extensional viscosity curves predicted by the PTT constitutive model used to generate the simulated pressure drop curves. The analytical techniques examined here are shown to provide reasonably accurate estimates of the steady extensional viscosity. This work also clearly demonstrates the advantage of using variable power-law coefficients for the rheological properties, used as inputs to the analyses, to capture the extensional behavior at deformation rates below the power law region more accurately. Received: 23 July 1999/Accepted: 24 November 1999     

Numerical investigation of a perturbed swirling annular two-phase jet

A swirling annular gas–liquid two-phase jet flow system has been investigated by solving the compressible, time-dependent, non-dimensional Navier–Stokes equations using highly accurate numerical methods. The mathematical formulation for the flow system is based on an Eulerian approach with mixed-fluid treatment while an adjusted volume of fluid method is utilised to account for the gas compressibility. Surface tension effects are captured by a continuum surface force model. Swirling motion is applied at the inlet while a small helical perturbation is also applied to initiate the instability. Three-dimensional spatial direct numerical simulation has been performed with parallelisation of the code based on domain decomposition. The results show that the flow is characterised by a geometrical recirculation zone adjacent to the nozzle exit and by a central recirculation zone further downstream. Swirl enhances the flow instability and vorticity and promotes liquid dispersion in the cross-streamwise directions. A dynamic precessing vortex core is developed demonstrating that the growth of such a vortex in annular configurations can be initiated even at low swirl numbers, in agreement with experimental findings. Analysis of the averaged results revealed the existence of a geometrical recirculation zone and a swirl induced central recirculation zone in the flow field.     

Numerical simulations of Boger fluids through different contraction configurations for the development of a measuring system for extensional viscosity

This paper reports the flow behaviour of Newtonian and Boger fluids through various axisymmetric contraction configurations by means of numerical predictions. A principal aim has been to evaluate the geometrical design choice of the hyperbolic contraction flow. The FENE-CR model has been used to reflect the behaviour of Boger fluids, with constant shear viscosity, finite (yet large) extensional viscosity and less than quadratic first normal stress difference. Numerical calculations have been performed on six different contraction configurations to evaluate an optimized geometry for measuring extensional viscosity in uniaxial extensional flow. The influence of a sharp or rounded recess-corner on the nozzle has also been investigated. Few commercial measuring systems are currently available for measurement of the extensional rheology of medium-viscosity fluids, such as foods and other biological systems. In this context, a technique based on the hyperbolic contraction flow would be a suitable alternative. The pressure drop, the velocity field, the first normal stress difference and the strain rate across the geometry have each been evaluated for Newtonian and Boger fluids. This numerical study has shown that the hyperbolic configuration is superior to the other geometry choices in achieving a constant extension rate. In this hyperbolic configuration, no vortices are formed, the measuring range is broader and the strain rate is constant throughout the geometric domain, unlike in the alternative configurations tested. The difference between sharp and rounded recess-corner configurations proved to be negligible and a rise in excess pressure drop (epd) for increasing deformation rates has been observed.     

On one model of cyclic ejections of magma in the process of the explosive volcanic eruption

The dynamics of the structure of a cavitating magma flow behind the decompression wave front is experimentally studied by the method of hydrodynamic shock tubes. It is demonstrated that a discrete system of intensely cavitating zones with alternation of low and high densities of the gas phase can be formed at a certain regime of shock-wave loading of the examined fluid sample. Based on the results of a numerical analysis of the formation of an anomalous zone in the cavitating magma flow with anomalously high flow characteristics exceeding the values of these characteristics outside this zone at least by an order of magnitude, a model of an instantaneous transformation of the cavitating magma in the anomalous zone to a gas-droplet system, its ejections, and formation of a free surface on the interface is proposed. A numerical analysis shows that the characteristic wave structure and the anomalous saturation zone are fairly rapidly reconstructed in the vicinity of this free surface of the flow part remaining in the conduit after the ejections, and the above-mentioned jumps of flow characteristics are again formed in the anomalous zone.     

Steady non-swirling axisymmetric flows: flow classification and rheological consequences

For steady non-swirling axisymmetric flow () of an incompressible fluid two invariants of the rate of strain dyadic D are introduced, which directly enter into the expression for D. This being the case they - in conjunction with the vorticity - allow a flow classification into strong and weak flows. For a generalized Newtonian fluid an expression for the viscosity function is listed, which reduces for model fluids to correct results in shearing and, respectively, extensional flow. A possible modification of is proposed, which involves the relative vorticity as well (quasi-Newtonian fluid), since this allows to adjust itself to the local nature of the flow. As such it should prove useful for numerical calculation. Received April 23, 1998     

The effect of step-stretch parameters on capillary breakup extensional rheology (CaBER) measurements

Extensional rheometry has only recently been developed into a commercially available tool with the introduction of the capillary breakup extensional rheometer (CaBER). CaBER is currently being used to measure the transient extensional viscosity evolution of mid to low-viscosity viscoelastic fluids. The elegance of capillary breakup extensional experiments lies in the simplicity of the procedure. An initial step-stretch is applied to generate a fluid filament. What follows is a self-driven uniaxial extensional flow in which surface tension is balanced by the extensional stresses resulting from the capillary thinning of the liquid bridge. In this paper, we describe the results from a series of experiments in which the step-stretch parameters of final length, and the extension rate of the stretch were varied and their effects on the measured extensional viscosity and extensional relaxation time were recorded. To focus on the parameter effects, well-characterized surfactant wormlike micelle solutions, polymer solutions, and immiscible polymer blends were used to include a range of characteristic relaxation times and morphologies. Our experimental results demonstrate a strong dependence of extensional rheology on step-stretch conditions for both wormlike micelle solutions and immiscible polymer blends. Both the extensional viscosity and extensional relaxation time of the wormlike micelle solutions were found to decrease with increasing extension rate and strain of the step-stretch. For the case of the immiscible polymer blends, fast step-stretches were found to result in droplet deformation and an overshoot in the extensional viscosity which increased with increasing strain rates. Conversely, the polymer solutions tested were found to be insensitive to step-stretch parameters. In addition, numerical simulations were performed using the appropriate constitutive models to assist in both the interpretation of the CaBER results and the optimization of the experimental protocol. From our results, it is clear that any rheological results obtained using the CaBER technique must be properly considered in the context of the stretch parameters and the effects that preconditioning has on viscoelastic fluids.     

Swirling and Impinging Effects in an Annular Nonpremixed Jet Flame

The effects of swirl and downstream wall confinement on an annular nonpremixed flame were investigated using direct numerical simulation (DNS). Fully three-dimensional parallel DNS was performed employing high-order numerical methods and high-fidelity boundary conditions to solve governing equations for variable-density flow and finite-rate Arrhenius chemistry. Three swirl numbers have been examined: 0 (without swirl), 0.4 and 0.8, while the effects of downstream wall confinement have been examined for swirl numbers 0 and 0.4. Results have been presented in terms of instantaneous and time-averaged flow quantities, which have also been analysed using energy spectra and proper orthogonal decomposition (POD). Effects of swirl on the fluid dynamic behaviour of the annular nonpremixed flame were found to be significant. The fluid dynamic behaviour of the flame is greatly affected by the interaction between the geometrical recirculation zone (GRZ) near the jet nozzle exit due to the annular configuration, the central recirculation zone (CRZ) associated with swirl, the unsteady vortical structures in the jet column due to the shear instability, and the downstream wall confinement. Depending on the degree of swirl, the GRZ near the burner mouth and the CRZ may co-exist or one zone may be overwhelmed by another. At a moderate swirl number, the co-existence leads to a flame with strong reaction attached to the burner mouth; while at a high swirl number, the CRZ dominates over the GRZ. The precessing vortex core was observed to exist in the swirling flow fields. The Nusselt number distribution of the annular impinging flames differs from that of round impinging jets. The POD analysis revealed that wall effects on the flow field are mainly associated with the higher mode numbers.     

THE SIMULATION OF SOME MODEL VISCOELASTIC EXTENSIONAL FLOWS

The numerical simulation of three model viscoelastic extensional flows is considered: sink flow, model draw-down and conical section draw-down. A transient finite element scheme with a pressure correction method is employed to analyse the numerical treatment of such problems for Oldroyd- Band Phan-Thien/Tanner constitutive models. Both decoupled and coupled formulations are compared for these highly convective flows and effective mechanisms are proposed for removing numerical oscillations in the temporally developing solution. In pure viscoelastic extensional flow from an initial stress-free state, the maximum stress level attained decreases with increase in material relaxation time. When this is followed by stress relaxation, as in conical section draw-down, increasing the relaxation time inhibits stress decay.     

Numerical investigation of transient contraction flows for worm-like micellar systems using Bautista–Manero models

This study is concerned with the numerical modelling of the Modified Bautista–Manero (MBM) model, for both steady-state and transient solutions in planar 4:1 contraction flow. This model was proposed to represent the structured composition and behaviour of worm-like micellar systems which have importance in industrial oil-reservoir recovery applications. A parameter sensitivity analysis for the rheology of this model is presented in both transient and steady response, covering pure shear and uniaxial extension. In addition, some features in evolutionary flow-structure are demonstrated in contraction flows due to the influence and imposition of start-up transient boundary conditions. The different effects of various model parameter choices are described through transient field response, stress and viscosity fields in the contraction flow setting. Distinction may be drawn between fluid response in the strong/moderate extension hardening regimes by matching both steady-state and transient shear and extensional viscosity peaks, contrasting between micellar (MBM) models against network-based counterparts Phan-Thien/Tanner (PTT). Simulations are performed with a hybrid finite volume/element algorithm. The momentum and continuity equations are solved by a Taylor–Galerkin/pressure-correction finite element method, whilst the constitutive equation is dealt with by a cell-vertex finite volume algorithm.     

Effect of spatial discretization schemes on numerical solutions of viscoelastic fluid flows

This study examines the effect of discretization schemes for the convection term in the constitutive equation on numerical solutions of viscoelastic fluid flows. For this purpose, a temporally evolving mixing layer, a two-dimensional vortex pair interacting with a wall, and a fully developed turbulent channel flow are selected as test cases, and eight different discretization schemes are considered. Among them, the first-order upwind difference scheme (UD) and artificial diffusion scheme (AD), which are commonly used in the literature, show most stable and smooth solutions even for highly extensional flows. However, the stress fields are smeared too much by these schemes and the corresponding flow fields are quite different from those obtained by higher-order upwind difference schemes. Among higher-order upwind difference schemes investigated in this study, a third-order compact upwind difference scheme (CUD3) with locally added AD shows stable and most accurate solutions for highly extensional flows even at relatively high Weissenberg numbers.     

Efficient Thermo-Chemistry Tabulation for Non-Premixed Combustion at High-Pressure Conditions

Propellant injection and turbulent combustion in high-pressure engines is often dominated by real-gas effects. However, previous studies suggested that the departure of the fluid properties from an ideal gas behavior has only a limited effect on the laminar flame structure. This is due to the fact that chemical reactions take place in the flame zone where the temperature is sufficiently high and molecular interactions are negligible, i.e., the ideal gas assumption is valid. On the other hand, various experimental and numerical studies of injection processes at high-pressure conditions demonstrated that real-gas effects can have a strong impact on the turbulent flow. Mixing is influenced by the rapid change of fluid properties. In this work, we exploit the gap in the fidelity of the thermodynamics model needed to describe the laminar flame structure and that needed to describe the turbulent flow field. We then propose a new real-gas flamelet model with increased numerical performance. The computational cost of the new formulation is not significantly higher than that of an ideal gas simulation. The performance of the method is analyzed and the error that is introduced by our assumptions is assessed by comparison to more complete modeling. Finally, the method is used to simulate a turbulent jet flame emanating from a coaxial injector at supercritical pressure and cryogenic oxidizer temperature. The results are compared with experimental OH^? images giving evidence of the suitability of the present method.