A Cartesian grid method for two‐phase gel dynamics on an irregular domain

We develop a numerical method for simulating models of two‐phase gel dynamics in an irregular domain using a regular Cartesian grid. The models consist of transport equations for the volume fractions of the two phases, polymer network and solvent; coupled momentum equations for the two phases; and a volume‐averaged incompressibility constraint. Multigrid with Vanka‐type box relaxation scheme is used as a preconditioner for the Krylov subspace solver (GMRES) to solve the momentum and incompressibility equations. Ghost points are used to enforce no‐slip boundary conditions for the velocity field of each phase, and no‐flux boundary conditions for the volume fractions. The behavior of the new method, including its rate of convergence, is explored through numerical experiments for a problem in which strong phase separation develops from an initially (almost) homogeneous phase distribution. We also use the method to explore situations, motivated by biology, which show that imposed boundary velocities can cause substantial redistribution of network and solvent. Copyright © 2010 John Wiley & Sons, Ltd.     

A parallel computational method for simulating two‐phase gel dynamics on a staggered grid

We develop a parallel computational algorithm for simulating models of gel dynamics where the gel is described by two phases, a networked polymer and a fluid solvent. The models consist of transport equations for the two phases, two coupled momentum equations, and a volume‐averaged incompressibility constraint. Multigrid with Vanka‐type box‐relaxation scheme is used as preconditioner for the Krylov subspace solver (GMRES) to solve the momentum and incompressibility equations. Through numerical experiments of a model problem, the efficiency, robustness and scalability of the algorithm are illustrated. Copyright © 2008 John Wiley & Sons, Ltd.     

A Model for Flow and Deformation in Unsaturated Swelling Porous Media

A thermomechanical theory for multiphase transport in unsaturated swelling porous media is developed on the basis of Hybrid Mixture Theory (saturated systems can also be modeled as a special case of this general theory). The aim is to comprehensively and non-empirically describe the effect of viscoelastic deformation on fluid transport (and vice versa) for swelling porous materials. Three phases are considered in the system: the swelling solid matrix s, liquid l, and air a. The Coleman–Noll procedure is used to obtain the restrictions on the form of the constitutive equations. The form of Darcy’s law for the fluid phase, which takes into account both Fickian and non-Fickian transport, is slightly different from the forms obtained by other researchers though all the terms have been included. When the fluid phases interact with the swelling solid porous matrix, deformation occurs. Viscoelastic large deformation of the solid matrix is investigated. A simple form of differential-integral equation is obtained for the fluid transport under isothermal conditions, which can be coupled with the deformation of the solid matrix to solve for transport in an unsaturated system. The modeling theory thus developed, which involves two-way coupling of the viscoelastic solid deformation and fluid transport, can be applied to study the processing of biopolymers, for example, soaking of foodstuffs and stress-crack predictions. Moreover, extension and modification of this modeling theory can be applied to study a vast variety of problems, such as drying of gels, consolidation of clays, drug delivery, and absorption of liquids in diapers.     

不可压饱和粘弹性多孔介质的变分原理及其无网格模拟

基于饱和多孔介质理论,在固相和液相微观不可压,固相骨架小变形且满足线性粘弹性积分型本构关系的假定下,建立了流体饱和粘弹性多孔介质动力响应的若干Gurtin型变分原理,包括Hu-Washizu变分原理.利用所建立的变分原理,导出了流体饱和粘弹性多孔介质动力响应无网格数值模拟的离散控制方程,此方程是一个关于时间的对称微分方程组,便于分析计算.作为数值例子,研究了流体饱和粘弹性多孔柱体的一维动力响应,数值结果揭示了流体饱和粘弹性多孔柱体中波的传播特性以及固相粘性的影响.     

Thermomechanics of Swelling Biopolymeric Systems

A two-scale theory for the swelling biopolymeric media is developed. At the microscale, the solid polymeric matrix interacts with the solvent through surface contact. The relaxation processes within the polymeric matrix are incorporated by modeling the solid phase as viscoelastic and the solvent phase as viscous at the mesoscale. We obtain novel equations for the total stress tensor, chemical potential of the solid phase, heat flux and the generalized Darcy's law all at the mesoscale. The constitutive relations are more general than those previously developed for the swelling colloids. The generalized Darcy's law could be used for modeling non-Fickian fluid transport over a wide range of liquid contents. The form of the generalized Fick's law is similar to that obtained in earlier works involving colloids. Using two-variable expansions, thermal gradients are coupled with the strain rate tensor for the solid phase and the deformation rate tensor for the liquid phase. This makes the experimental determination of the material coefficients easier and less ambiguous.     

Coupled Solvent and Heat Transport of a Mixture of Swelling Porous Particles and Fluids: Single Time-Scale Problem

A three-spatial scale, single time-scale model for both moisture and heat transport is developed for an unsaturated swelling porous media from first principles within a mixture theoretic framework. On the smallest (micro) scale, the system consists of macromolecules (clay particles, polymers, etc.) and a solvating liquid (vicinal fluid), each of which are viewed as individual phases or nonoverlapping continua occupying distinct regions of space and satisfying the classical field equations. These equations are homogenized forming overlaying continua on the intermediate (meso) scale via hybrid mixture theory (HMT). On the mesoscale the homogenized swelling particles consisting of the homogenized vicinal fluid and colloid are then mixed with two bulk phase fluids: the bulk solvent and its vapor. At this scale, there exists three nonoverlapping continua occupying distinct regions of space. On the largest (macro) scale the saturated homogenized particles, bulk liquid and vapor solvent, are again homogenized forming four overlaying continua: doubly homogenized vicinal fluid, doubly homogenized macromolecules, and singly homogenized bulk liquid and vapor phases. Two constitutive theories are developed, one at the mesoscale and the other at the macroscale. Both are developed via the Coleman and Noll method of exploiting the entropy inequality coupled with linearization about equilibrium. The macroscale constitutive theory does not rely upon the mesoscale theory as is common in other upscaling methods. The energy equation on either the mesoscale or macroscale generalizes de Vries classical theory of heat and moisture transport. The momentum balance allows for flow of fluid via volume fraction gradients, pressure gradients, external force fields, and temperature gradients.     

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.     

A thermo-mechanically coupled theory for fluid permeation in elastomeric materials: Application to thermally responsive gels

An elastomeric gel is a cross-linked polymer network swollen with a solvent, and certain gels can undergo large reversible volume changes as they are cycled about a critical temperature. We have developed a continuum-level theory to describe the coupled mechanical deformation, fluid permeation, and heat transfer of such thermally responsive gels. In discussing special constitutive equations we limit our attention to isotropic materials, and consider a model based on a Flory–Huggins model for the free energy change due to mixing of the fluid with the polymer network, coupled with a non-Gaussian statistical–mechanical model for the change in configurational entropy—a model which accounts for the limited extensibility of polymer chains. We have numerically implemented our theory in a finite element program. We show that our theory is capable of simulating swelling, squeezing of fluid by applied mechanical forces, and thermally responsive swelling/de-swelling of such materials.     

Finite volume multigrid method of the planar contraction flow of a viscoelastic fluid

This paper reports on a numerical algorithm for the steady flow of viscoelastic fluid. The conservative and constitutive equations are solved using the finite volume method (FVM) with a hybrid scheme for the velocities and first‐order upwind approximation for the viscoelastic stress. A non‐uniform staggered grid system is used. The iterative SIMPLE algorithm is employed to relax the coupled momentum and continuity equations. The non‐linear algebraic equations over the flow domain are solved iteratively by the symmetrical coupled Gauss–Seidel (SCGS) method. In both, the full approximation storage (FAS) multigrid algorithm is used. An Oldroyd‐B fluid model was selected for the calculation. Results are reported for planar 4:1 abrupt contraction at various Weissenberg numbers. The solutions are found to be stable and smooth. The solutions show that at high Weissenberg number the domain must be long enough. The convergence of the method has been verified with grid refinement. All the calculations have been performed on a PC equipped with a Pentium III processor at 550 MHz. Copyright © 2001 John Wiley & Sons, Ltd.     

Surfactant spreading on a thin weakly viscoelastic film

A mathematical model is presented for surfactant-driven thin weakly viscoelastic film flows on a flat, impermeable plane. The Oldroyd-B constitutive relation is used to model the viscoelastic fluid. Lubrication theory and a perturbation expansion in powers of the Weissenberg number (We) are employed, which give rise to non-linear coupled evolution equations governing the transport of insoluble surfactant and thin liquid film thickness. Spreading on a Newtonian film is recovered to leading order and corrections to viscoelasticity are obtained at order We. These equations are solved numerically over a wide range of viscosity ratio (ratio of solvent viscosity to the sum of solvent and polymeric viscosities), pre-existing surfactant level and Peclet number (Pe). The effect of viscoelasticity on surfactant transport and fluid flow is investigated and the mechanisms underlying this effect are explored. Shear stress, streamwise normal stress and the temporal rate of change of extra shear stress generated from gradients in surfactant concentration dominate thin viscoelastic film flows whereas only shear stresses play a role in Newtonian thin film flows. Our results also reveal that, for weak viscoelasticity, the influence of viscosity ratio on the evolution of surfactant concentration and film thickness can be significant and varies considerably, depending on the concentration of pre-existing surfactant and surfactant surface diffusivity.     

Interaction of a deformable solid with two-phase flows: An Eulerian-based numerical model for fluid-structure interaction using the level contour reconstruction method

We describe the formulation of a method for fluid-structure interaction involving the coupling of moving and/or flexible solid structures with multiphase flows in the framework of the Level Contour Reconstruction Method. We present an Eulerian-based numerical procedure for tracking the motion and interaction of a liquid-gas interface with a fluid-solid interface in the Lagrangian frame together with the evaluation of the fluid transport equations coupled to those for the solid transport, namely the left Cauchy-Green strain tensor field, in the Eulerian frame. To prevent excessive dissipation due to the convective nature of the solid transport equation, a simple incompressibility constraint for the strain field is enforced. A single grid structure is used for both the fluid and solid phases which allows for a simple and natural coupling of the fluid and solid dynamics. Several benchmark tests are performed to show the accuracy of the numerical method and which demonstrate accurate results compared to several of those in the existing literature. In particular we show that surface tension effects including contact line dynamics on the deforming solid phase can be properly simulated. The three-phase interaction of a droplet impacting on a flexible cantilever is investigated in detail. The simulations follow the detailed motion of the droplet impact (and subsequent deformation, breakup, and fall trajectory) along with the motion of the deformable solid cantilever due to its own weight as well as due to the force of the droplet impact.     

Nonlinear waves in incompressible viscoelastic Maxwell medium

In this paper we study two-dimensional flows of incompressible viscoelastic Maxwell media with Jaumann corotational derivative in the rheological constitutive law. In the general case, due to the incompressibility condition, the equations of motion have both real and complex characteristics. Group properties of this system are studied. On this basis, two submodels of the Maxwell model are selected, which can be reduced to hyperbolic ones. More precisely, we consider plane shear flow between two parallel planes and Couette type flow caused by the inertial cylinder rotation. As a result, we obtain the closed systems of three equations of mixed type, which describe nonlinear transverse waves in an incompressible Maxwell fluid. It is demonstrated that discontinuities can develop in elastic media even from smooth initial data. Stability of shocks in the Maxwell fluid with and without retardation time is discussed.     

A finite element analysis of laminar unsteady flow of viscoelastic fluids through channels with non-uniform cross-sections

The effects of non-Newtonian behaviour of a fluid and unsteadiness on flow in a channel with non-uniform cross-section have been investigated. The rheological behaviour of the fluid is assumed to be described by the constitutive equation of a viscoelastic fluid obeying the Oldroyd-B model. The finite element method is used to analyse the flow. The novel features of the present method are the adoption of the velocity correction technique for the momentum equations and of the two-step explicit scheme for the extra stress equations. This approach makes the computational scheme simple in algorithmic structure, which therefore implies that the present technique is capable of handling large-scale problems. The scheme is completed by the introduction of balancing tensor diffusivity (wherever necessary) in the momentum equations. It is important to mention that the proper boundary condition for pressure (at the outlet) has been developed to solve the pressure Poisson equation, and then the results for velocity, pressure and extra stress fields have been computed for different values of the Weissenberg number, viscosity due to elasticity, etc. Finally, it is pertinent to point out that the present numerical scheme, along with the proper boundary condition for pressure developed here, demonstrates its versatility and suitability for analysing the unsteady flow of viscoelastic fluid through a channel with non-uniform cross-section.     

A numerical continuous model for the hydrodynamics of fluid particle systems

In order to understand the hydrodynamic interactions that can appear in a fluid particle motion, an original method based on the equations governing the motion of two immiscible fluids has been developed. These momentum equations are solved for both the fluid and solid phases. The solid phase is assumed to be a fluid phase with physical properties, such as its behaviour can be assimilated to that of pseudo‐rigid particles. The only unknowns are the velocity and the pressure defined in both phases. The unsteady two‐dimensional momentum equations are solved by using a staggered finite volume formulation and a projection method. The transport of each particle is solved by using a second‐order explicit scheme. The physical model and the numerical method are presented, and the method is validated through experimental measurements and numerical results concerning the flow around a circular cylinder. Good agreement is observed in most cases. The method is then applied to study the trajectory of one settling particle initially off‐centred between two parallel walls and the corresponding wake effects. Different particle trajectories related to particulate Reynolds numbers are presented and commented. A two‐body interaction problem is investigated too. This method allows the simulation of the transport of particles in a dilute suspension in reasonable time. One of the important features of this method is the computational cost that scales linearly with the number of particles. Copyright © 1999 John Wiley & Sons, Ltd.     

Modified incompressible SPH method for simulating free surface problems

An incompressible smoothed particle hydrodynamics (I-SPH) formulation is presented to simulate free surface incompressible fluid problems. The governing equations are mass and momentum conservation that are solved in a Lagrangian form using a two-step fractional method. In the first step, velocity field is computed without enforcing incompressibility. In the second step, a Poisson equation of pressure is used to satisfy incompressibility condition. The source term in the Poisson equation for the pressure is approximated, based on the SPH continuity equation, by an interpolation summation involving the relative velocities between a reference particle and its neighboring particles. A new form of source term for the Poisson equation is proposed and also a modified Poisson equation of pressure is used to satisfy incompressibility condition of free surface particles. By employing these corrections, the stability and accuracy of SPH method are improved. In order to show the ability of SPH method to simulate fluid mechanical problems, this method is used to simulate four test problems such as 2-D dam-break and wave propagation.     

An incompressible SPH method for simulation of unsteady viscoelastic free-surface flows

In this paper, an incompressible smoothed particle hydrodynamics (SPH) method is presented to solve unsteady free-surface flows. Both Newtonian and viscoelastic fluids are considered. In the case of viscoelastic fluids, both the Maxwell and Oldroyd-B models are investigated. The proposed SPH method uses a Poisson pressure equation to satisfy the incompressibility constraints. The solution algorithm is an explicit predictor-corrector scheme and employs an adaptive smoothing length based on density variations. To alleviate the numerical difficulties encountered when fluid is highly stretched, an artificial stress term is incorporated into the momentum equation which reduces the risk of unrealistic fractures in the material. Two challenging test cases, the impacting drop and the jet buckling problems, are solved to demonstrate the capability of the proposed scheme in handling viscoelastic flows with complex free surfaces. The jet buckling test case was solved for a wide range of Weissenberg numbers. It was shown that in all cases the method is stable and fairly accurate and agrees well with the available data.     

FAST ITERATIVE SOLVERS FOR THE DISCRETIZED INCOMPRESSIBLE NAVIER–STOKES EQUATIONS

In this paper some iterative solution methods of the GMRES type for the discretized Navier–Stokes equations are treated. The discretization combined with a pressure correction scheme leads to two different types of systems of linear equations: the momentum system and the pressure system. These systems may be coupled to one or more transport equations. For every system we specify a particular ILU-type preconditioner and show how to vectorize these preconditions. Finally, some numerical experiments to show the efficiency of the proposed methods are presented.     

Finite Elastic Deformations in Liquid-Saturated and Empty Porous Solids

Based on the Theory of Porous Media (TPM), a formulation of a fluid-saturated porous solid is presented where both constituents, the solid and the fluid, are assumed to be materially incompressible. Therefore, the so-called point of compaction exists. This deformation state is reached when all pores are closed and any further volume compression is impossible due to the incompressibility constraint of the solid skeleton material. To describe this effect, a new finite elasticity law is developed on the basis of a hyperelastic strain energy function, thus governing the constraint of material incompressibility for the solid material. Furthermore, a power function to describe deformation dependent permeability effects is introduced.After the spatial discretization of the governing field equations within the finite element method, a differential algebraic system in time arises due to the incompressibility constraint of both constituents. For the efficient numerical treatment of the strongly coupled nonlinear solid-fluid problem, a consistent linearization of the weak forms of the governing equations with respect to the unknowns must be used.     

不可压饱和粘弹性多孔介质固结问题的Gurtin型变分原理

基于多孔介质理论,在固相骨架和孔隙流体微观不可压,固相骨架小变形且满足线性粘弹性积分型本构关系的假定下,利用卷积积分的性质,本文首先建立了以固相骨架位移、孔隙流体相对速度和孔隙流体压力为宗量的流体饱和粘弹性多孔介质固结问题的一个Gurtin型变分原理.其次,利用Lagrange乘子法解除相关的变分约束条件,建立了流体饱和粘弹性多孔介质固结问题的若干广义Gurtin型变分原理,包括第三类的Hu-Washizu型变分原理.最后,简单讨论了等价初边值问题的相应变分原理.这些Gurtin型变分原理的建立不仅丰富了饱和粘弹性多孔介质的相关理论,而且为相关数值模拟方法,如有限元法、无网格法等的建立奠定了理论基础.     

Viscoelastic MHD flow boundary layer over a stretching surface with viscous and ohmic dissipations

In this study the momentum and heat transfer characteristics in an incompressible electrically conducting viscoelastic boundary layer fluid flow over a linear stretching sheet are considered. Highly non-linear momentum and thermal boundary layer equations are reduced to set of nonlinear ordinary differential equations by appropriate transformation.