Finite element computations of viscoelastic two-phase flows using local projection stabilization

A three-field local projection stabilized (LPS) finite element method is developed for computations of a three-dimensional axisymmetric buoyancy driven liquid drop rising in a liquid column where one of the liquid is viscoelastic. The two-phase flow is described by the time-dependent incompressible Navier-Stokes equations, whereas the viscoelasticity is modeled by the Giesekus constitutive equation in a time-dependent domain. The arbitrary Lagrangian-Eulerian (ALE) formulation with finite elements is used to solve the governing equations in the time-dependent domain. Interface-resolved moving meshes in ALE allows to incorporate the interfacial tension force and jumps in the material parameters accurately. A one-level LPS based on an enriched approximation space and a discontinuous projection space is used to stabilize the numerical scheme. A comprehensive numerical investigation is performed for a Newtonian drop rising in a viscoelastic fluid column and a viscoelastic drop rising in a Newtonian fluid column. The influence of the viscosity ratio, Newtonian solvent ratio, Giesekus mobility factor, and the Eötvös number on the drop dynamics are analyzed. The numerical study shows that beyond a critical Capillary number, a Newtonian drop rising in a viscoelastic fluid column experiences an extended trailing edge with a cusp-like shape and also exhibits a negative wake phenomena. However, a viscoelastic drop rising in a Newtonian fluid column develops an indentation around the rear stagnation point with a dimpled shape.     

Some exact solutions of the ferrohydrodynamics equations

The theoretical study of nonisothermal flows of magnetizable liquids presents serious matheical difficulties, which are intensified as compared to to the study of normal liquids by the necessity of simultaneous solution of both the hydrodynamics and Maxwell's equations, with corresponding boundary conditions for the magnetic field. Thus, in most cases problems of this type are solved by neglecting the effect of the liquid's nonisothermal state on the field distribution within the liquid, and also, as a rule, with use of an approximate solution for Maxwell's equations and fulfillment of the boundary conditions for the field [1–3]. The present study will present easily realizable practical formulations of the problem which permit exact one-dimensional solutions of the complete system of the equations of thermomechanic s of electrically nonconductive incompressible Newtonian magnetizable liquids with constant transfer coefficients. A common feature of the formulations is the presence of a longitudinal temperature gradient along the boundaries along which liquid motion is accomplished. Plane-parallel convective flows of this type in a conventional liquid and their stability were considered in [4–6].Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 3, pp. 126–133, May–June, 1979.     

Selective withdrawal of polymer solutions: Experiments

Selective withdrawal refers to the process of drawing one or both components of stratified fluids through a tube placed near their interface. This paper reports an experimental study of selective withdrawal of viscous and viscoelastic liquids under air. The key mechanism of interest is how the viscoelasticity in the bulk liquid affects the evolution of the free surface. This is investigated by comparing the interfacial behavior between a Newtonian silicone oil and two dilute polymer solutions. While the surface undergoes smooth and gradual deformation for Newtonian liquids, for the polymer solutions there is a critical transition where the surface forms a cusp from which an air jet emanates toward the suction tube. This transition shows a hysteresis when the flow rate or location of the tube is varied. In the subcritical state, the surface of polymer solutions deform much more than its Newtonian counterpart but the deformation is more localized. The interfacial behavior of the polymer solutions can be attributed to the large polymer stress that develops under the surface because of predominantly extensional deformation.     

Jet with variable fluid properties: Free jet and dissipative jet

An analysis of the laminar jet of an incompressible Newtonian fluid emerging from a narrow slot or a circular hole, where the physical properties like viscosity and thermal conductivity depends upon the temperature, is given. Both the cases: the case of In the absence of viscous heat dissipation and the case of In the presence of viscous heat dissipation are considered. The governing partial differential equations of the flow problem are transformed into the ordinary differential equations by group theoretic technique. The Runge–Kutta method is applied to obtained numerical solution of the transformed ordinary differential equations.     

Numerical simulation of shear-thinning flow problems in mixing vessels

In many applications of chemical and process engineering numerous important flow situations appear, in which the shear-thinning properties of the fluid dominate the normal stress effects and time-dependent elastic effects. A mixed finite element method for steady incompressible flow is presented taking as a basis the model of a generalized Newtonian fluid. An augmented Lagrangian functional is constructed corresponding to the equations of motion and the continuity constraint together with appropriate boundary conditions. The treatment of the resulting nonlinear system of equations by the Newton-Raphson scheme is made for general models of the viscosity function. The linear subproblems are solved by Uzawa's algorithm. The method is applied to the numerical simulation of various mixing problems in cylindrical unbaffled vessels. The computations were performed for a real polymer liquid (CMC in water), where the viscosity data were found experimentally and were fitted to a suitable mathematical model. Several numerical results are discussed and as far as possible compared with experimental data.This work was supported by Deutsche Forschungsgemeinschaft.     

Nonisothermal model of glass fiber drawing stability

Draw resonance is caused by a constant speed winder that leads to non-constant axial forces (Schultz, 1984). The well studied isothermal Newtonian fiber drawing predicts very modest critical draw ratios (around 20, much less than the typical production draw ratios for glass fibers of 10³ – 10⁵). The nonisothermal fiber drawing model presented here shows that cooling along the spin line strongly stabilizes the process. However, we show that the conclusion of Shah and Pearson (1972a,b) that non-isothermal Newtonian fiber spinning is unconditionally stable is based on non-converged numerical results. The choice of viscosity-temperature correlation function has a strong influence in determining the stability of the process. While viscoelasticity generally has an adverse effect on the stability, low viscoelasticity in the presence of extensional thinning helps to slightly improve the maximum critical draw ratio.Dedicated to the memory of Professor Tasos C. Papanastasiou     

Numerical investigation on a transition prediction model

A new transition prediction model is introduced,which couples the intermittency effect into the turbulence transport equations and takes the characteristics of fluid transition into consideration to mimic the exact process of transition.Test cases include a two-dimensional incompressible plate and a two-dimensional NACA0012 airfoil.Performance of this transition model for incompressible flows is studied,with numerical results consistent to experimental data.The requirement of grid resolution for this transition model is also studied.     

An ILU smoother for the incompressible Navier-Stokes equations in general co-ordinates

The solution of the incompressible Navier-Stokes equations in general two- and three-dimensional domains using a multigrid method is considered. Because a great variety of boundary-fitted grids may occur, robustness is at a premium. Therefore a new ILU smoother called CILU (collective ILU) is described, based on r-transformations. In CILU the matrix that is factorized is block-structured, with blocks corresponding to the set of physical variables. A multigrid algorithm using CILU as smoother is investigated. The performance of the algorithm in two and three dimensions is assessed by numerical experments. The results show that CILU is a good smoother for the incompressible Navier-Stokes equations discretized on general non-orthogonal curvilinear grids.     

On Solutions of Some Non-Linear Differential Equations Arising in Newtonian and Non-Newtonian Fluids

Some steady as well as unsteady solutions of the equations of motion for an incompressible Newtonian and non-Newtonian (second-grade) fluids are obtained by applying different methods including the Lie symmetry group method. The flows considered are axially symmetric with the swirling motion, and the governing equations for second-grade fluid flow have been modeled. Expressions for streamlines, velocity and vorticity components are constructed explicitly in each case. Exact analytical solutions in second-grade fluid are obtained and compared with the corresponding viscous solutions.     

Well-posedness of the problem of fiber suspension flows

The well-posedness of the equations governing the flow of fiber suspensions is studied. The fluid is assumed to be Newtonian and incompressible, and the presence of fibers is accounted for through the use of second- and fourth-order orientation tensors, which model the effects of the orientation of fibers in an averaged sense. The fourth-order orientation tensor is expressed in terms of the second-order tensor through various closure relations. It is shown that the linear closure relation leads to anomalous behavior, in that the rest state of the fluid is unstable, in the sense of Liapounov, for certain ranges of the fiber particle number. No such anomalies arise in the case of quadratic and hybrid closure relations. For the quadratic closure relation, it is shown that a unique solution exists locally in time for small data.     

Formation of radially expanding liquid sheet by impinging two round jets

A thin circular liquid sheet can be formed by impinging two identical round jets against each other. The liquid sheet expands to a certain critical radial distance and breaks. The unsteady process of the formation and breakup of the liquid sheet in the ambient gas is simulated numerically. Both liquid and gas are treated as incompressible Newtonian fluids. The flow considered is axisymmetric. The liquid-gas interface is modeled with a level set function. A finite difference scheme is used to solve the governing Navier-Stokes equations with physical boundary conditions. The numerical results show how a thin circular sheet can be formed and break at its circular edge in slow motion. The sheet continues to thin as it expands radially. Hence, the Weber number decreases radially. The Weber number is defined as ρu 2 h/σ, where ρ and σ are, respectively, the liquid density and the surface tension, and u and h are, respectively, the average velocity and the half sheet thickness at a local radial location in the liquid sheet. The numerical results show that the sheet indeed terminates at a radial location, where the Weber number reaches one as observed in experiments. The spatio-temporal linear theory predicts that the breakup is initiated by the sinuous mode at the critical Weber number We c =1, below which the absolute instability occurs. The other independent mode called the varicose mode grows more slowly than the sinuous mode according to the linear theory. However, our numerical results show that the varicose mode actually overtakes the sinuous mode during the nonlinear evolution, and is responsible for the final breakup. The linear theory predicts the nature of disturbance waves correctly only at the onset of the instability, but cannot predict the exact consequence of the instability.     

Front‐tracking by the level‐set and the volume penalization methods in a two‐phase microfluidic network

Two‐phase immiscible fluids in a two‐dimensional micro‐channels network are considered. The incompressible Stokes equations are used to describe the Newtonian fluid flow, while the Oldroyd‐B rheological model is used to capture the viscoelastic behavior. In order to perform numerical simulations in a complex geometry like a micro‐channels network, the volume penalization method is implemented. To follow the interface between the two fluids, the level‐set method is used, and the dynamics of the contact line is modeled by Cox law. Numerical results show the ability of the method to simulate two‐phase flows and to follow properly the contact line between the two immiscible fluids. Finally, simulations with realistic parameters are performed to show the difference when a Newtonian fluid is pushed by a viscoelastic fluid instead of a Newtonian one. Copyright © 2015 John Wiley & Sons, Ltd.     

Flow of a viscous liquid in a bent tube of finite length that undergoes deformation

A study is made of the closed flow of a viscous incompressible Newtonian liquid in a deformed tube with a section used as a peristaltic pump connected to a given hydraulic load. The time-average flow rate of the liquid (blood) in the external section is determined, together with the pressure difference across the load and the condition for the absence of hemolysis.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 2, pp. 177–179, March–April, 1981.     

Fundamental Solution of the System of Equations of Steady Oscillations in the Theory of Fluid-Saturated Porous Media

In the present paper the linear theory of the liquid-saturated porous medium consisting of a microscopically incompressible solid skeleton containing microscopically incompressible liquid is considered. The fundamental solution of the system of linear coupled partial differential equations of the steady oscillations of the porous solids is constructed in terms of elementary functions and some basic properties are established.     

Oscillations of a gas bubble in a non-Newtonian liquid under the action of an acoustic field

The evolution of the radius of a spherical cavitation bubble in an incompressible non-Newtonian liquid under the action of an external acoustic field is investigated. Non-Newtonian liquids having relaxation properties and also pseudoplastic and dilatant liquids with powerlaw equation of state are studied. The equations for the oscillation of the gas bubble are derived, the stability of its radial oscillation and its spherical form are investigated, and formulas are given for the characteristic frequency of oscillations of the cavitation hollow in a relaxing liquid. The equations are integrated numerically. It is shown that in a relaxing non-Newtonian liquid the viscosity may lead to the instability of the radial oscillations and the spherical form of the bubble. The results obtained here are compared with the behavior of a gas bubble in a Newtonian liquid.     

Exact solutions of equations of rotationally symmetric motion of an ideal incompressible liquid

A partially invariant solution of the Euler equations is considered, where the vertical component of velocity is a function of the vertical coordinate and time, whereas the remaining components of velocity and pressure are independent of the polar angle in a cylindrical coordinate system. Using the classification of equations obtained by analysis of an overdetermined system, we consider two hyperbolic systems: the first one describes the motion of a cylindrical layer of an ideal incompressible liquid under a punch, and the second system allows obtaining solutions in a halfcylinder with singularities at the axis of symmetry. A class of new exact solutions is obtained, which describe vortex motion of an ideal incompressible liquid, including the motion with singularities (sources of vortices) located along the axis of symmetry.     

Unsteady effects in weakly perturbed pulsatile channel flows

Pulsatile incompressible laminar flow in a plane channel with slight asymmetric deformation of the walls is considered on the basis of an asymptotic approach to the solution of the Navier-Stokes equations at large Reynolds numbers. It is shown that for extended deformations the change in the direction of the undisturbed flow in the wall regions is accompanied by a sharp increase in the amplitude of the perturbations; this is an essentially unsteady process. For flows with a small positive friction stress two classes of eigensolutions are obtained in the quasisteady approximation.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 2, pp. 96–104, March–April, 1987.The author is grateful to V. V. Sychev and A. I. Ruban for a number of useful consultations in the course of the work.     

A Local Model for the Energy Transfer in a Saturated Static Mixture

In the present work the transient energy transfer in a nonsaturated porous medium is studied, using a mixture theory viewpoint. The porous matrix is assumed homogeneous, rigid and isotropic, while the fluid is a Newtonian incompressible one and both are assumed static. Since the homogeneous matrix is not saturated, gradients of concentration are present. The porous medium and the fluid (a liquid) will be regarded as continuous constituents of a mixture that will have also a third constituent, an inert gas, assumed with zero mass density and thermal conductivity. The problem is described by a set of two partial differential equations which represent the energy balances for the fluid and the solid constituents. Isovalues for these two constituents are plotted, considering representative time instants and selected values for the energy equations coefficients and for the saturation.     

Isothermal extrusion of non-dilute fiber suspensions

The extrusion of a rod-like fiber suspension is a Newtonian solvent, as a first step to the fast and inexpensive production of composite materials, is investigated. The analysis is carried out by means of an integral constitutive equation for a non-dilute suspension, streamlined finite element for liquid with memory, and Newton iteration of nonlinear integro-differential equations. The predictions show substantial differences between dilute and nondilute fiber suspension regarding the processing conditions (pressure drop, velocity distribution, die-swell) and the resulting fiber orientation. Nondilute fiber suspensions exhibit substantial shear-thinning and negligible elasticity as evidenced by the small die-swell, and fiber concentration viscosity-thickening as evidenced by the large pressure drop. The fiber orientation is computed by solving the orientation distribution function along selected streamlines of the complex velocity field. It is shown that the fiber orientation far downstream can be made independent of the random fiber orientation at the inlet.     

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.