A new least-squares finite element model for combined Navier–Stokes and Darcy flows in geometrically complicated domains with solid and porous boundaries

In this paper we present a novel method for linking Navier–Stokes and Darcy equations along a porous inner boundary in a flow regime which is governed by both types of these equations. The method is based on a least-squares finite element technique and uses isoparametric C¹ continuous Hermite elements for domain discretization. We show that our technique is superior to previously developed models for the combined Navier–Stokes/Darcy flows. The previous works use weighted residual finite element procedures in conjunction with C⁰ elements which are inherently incapable of linking Navier–Stokes and Darcy equations. The paper includes the application of our model to a geometrically complicated axisymmetric slurry filtration system.     

An XFEM model for cracked porous media: effects of fluid flow and heat transfer

In this work, a numerical model is developed to investigate the influence of fluid flow and heat transfer on the thermo-mechanical response of a cracked porous media. The fluid flow, governed by the Darcy’s law, is discretized with the nonconforming finite element method. Time splitting is used with the energy conservation equation to solve the fluid and the solid phases separately. A combination of Discontinuous Galerkin (DG) and multi-point flux approximation methods is used to solve the advection-diffusion heat transfer equation in the fluid phase. While the conductive heat transfers equation in the solid phase is solved using the eXtended finite element method (XFEM) to better handle the temperature discontinuities and singularities caused by the cracks. Further, the resulted temperature is used as body force to solve the thermo-mechanical problem using the XFEM. In the post processing stage, the thermal stress intensity factor is computed using the interaction integral technique at each time step and used to validate the obtained results. A good agreement was found when the results were compared with the existing ones in the literature.     

Non-Darcian convection flow in a circular duct partially filled with porous medium

The flow and heat transfer in a circular duct bounded by a porous bed is considered. The entire flow region is divided into two zones. The clean fluid region is described by the Navier–Stokes equations, while the Brinkman-extended Darcy model is used in the flow through a porous bed. In either zone the momentum and temperature equations are coupled. In order to obtain a better insight into this complex problem, the Galerkin finite element analysis with quadratic polynomial approximations is applied. The behavior of the velocity and temperature is analyzed. The shear stress and the rate of heat transfer are also obtained for various governing parameters.     

An augmented Cartesian grid method for Stokes–Darcy fluid–structure interactions

A new finite difference method based on Cartesian meshes is proposed for solving the fluid–structure interaction between a fluid flow modeled by the Stokes equations and a porous media modeled by the Darcy's law. The idea is to introduce several augmented variables along the interface between the fluid flow and the porous media so that the problem can be decoupled as several Poisson equations. The augmented variables should be chosen so that the Beavers–Joseph–Saffman and other interface conditions are satisfied. In the discretization, the augmented variables have co‐dimension one compared with that of the primitive variables and are solved through the Schur complement system. A non‐trivial analytic solution with a circular interface is constructed to check second‐order convergency of the proposed method. Numerical examples with various interfaces and parameters are also presented. Some simulations show interesting behaviors of the fluid–structure interaction between the fluid flow and the porous media. The computational framework can be applied to other multi‐phase and multi‐physics problems. Copyright © 2015 John Wiley & Sons, Ltd.     

Discontinuous finite volume methods for the stationary Stokes–Darcy problem

In this paper, we analyze discontinuous finite volume methods for the stationary Stokes–Darcy problem that models coupled fluid flow and porous media flow. The discontinuous finite volume methods are combinations of finite volume method and discontinuous Galerkin method with three interior penalty types (incomplete symmetric, nonsymmetric, and symmetric), briefly, using discontinuous functions as trial functions in the finite volume method. Optimal error estimates in broken H¹ norm are obtained for the three discontinuous finite volume methods. Optimal error estimates in the standard L² norm are derived for the symmetric interior penalty discontinuous finite volume method. Numerical experiments are presented to confirm the theoretical results with non‐matching meshes across the common interface of Stokes region and Darcy region. Copyright © 2015 John Wiley & Sons, Ltd.     

基于VOF/PLIC界面追踪方法的RTM工艺充模过程数值模拟

针对基于Darcy定律的树脂传递模塑(RTM)工艺的充模过程数值模拟的局限性,将纤维预制体内的充填流动作为两相流(树脂相和空气相)处理,在动量方程中考虑了惯性项和粘性项,采用有限体积方法(FVM)离散控制方程,并与VOF/PLIC界面追踪方法相结合,发展了求解树脂在纤维预制体内非稳态流动问题的数值模拟方法.在此基础上开发了RTM工艺的充模过程数值模拟程序,其算例的数值模拟结果与解析解或实验结果吻合良好,验证了此数值模拟方法的有效性和可靠性.     

Topology optimization for unsteady flow

A computational methodology for optimizing the conceptual layout of unsteady flow problems at low Reynolds numbers is presented. The geometry of the design is described by the spatial distribution of a fictitious material with continuously varying porosity. The flow is predicted by a stabilized finite element formulation of the incompressible Navier–Stokes equations. A Brinkman penalization is used to enforce zero‐velocities in solid material. The resulting parameter optimization problem is solved by a non‐linear programming method. The paper studies the feasibility of the material interpolation approach for optimizing the topology of unsteady flow problems. The derivation of the governing equations and the adjoint sensitivity analysis are presented. A design‐dependent stabilization scheme is introduced to mitigate numerical instabilities in porous material. The emergence of non‐physical artifacts in the optimized material distribution is observed and linked to an insufficient resolution of the flow field and an improper representation of the pressure field within solid material by the Brinkman penalization. Two numerical examples demonstrate that the designs optimized for unsteady flow differ significantly from their steady‐state counterparts. Copyright © 2011 John Wiley & Sons, Ltd.     

Estimation of the dynamic permeability of an assembly of permeable spherical porous particles using the cell model

This paper is concerned with the theoretical estimation of the dynamic permeability of an assembly of permeable porous particles subject to an oscillatory Stokes flow. A cell model is used to approximate the hydrodynamic interactions of particles. The flow field inside the representative porous sphere is governed by Darcy’s law and that within the fluid region by the unsteady Stokes equations. Faxén’s laws for drag and torque exerted on the representative particle are derived, and the results are compared in special cases with the existing literature. The dynamic permeability of the bed of particles is obtained and the variation in dynamic permeability with various parameters such as porosity, Darcy number and, frequency is analyzed.     

A least-squares front-tracking finite element method analysis of phase change with natural convection

This paper focuses on the numerical modelling of phase-change processes with natural convection. In particular, two-dimensional solidification and melting problems are studied for pure metals using an energy preserving deforming finite element model. The transient Navier–Stokes equations for incompressible fluid flow are solved simultaneously with the transient heat flow equations and the Stefan condition. A least-squares variational finite element method formulation is implemented for both the heat flow and fluid flow equations. The Boussinesq approximation is used to generate the bulk fluid motion in the melt. The mesh motion and mesh generation schemes are performed dynamically using a transfinite mapping. The consistent penalty method is used for modelling incompressibility. The effect of natural convection on the solid/liquid interface motion, the solidification rate and the temperature gradients is found to be important. The proposed method does not possess some of the false diffusion problems associated with the standard Galerkin formulations and it is shown to produce accurate numerical solutions for convection dominated phase-change problems.     

Simulation of 3D flow in porous media by boundary element method

In the present paper problem of natural convection in a cubic porous cavity is studied numerically, using an algorithm based on a combination of single domain and subdomain boundary element method (BEM). The modified Navier–Stokes equations (Brinkman-extended Darcy formulation with inertial term included) were adopted to model fluid flow in porous media, coupled with the energy equation using the Boussinesq approximation. The governing equations are transformed by the velocity–vorticity variables formulation which separates the computation scheme into kinematic and kinetic parts. The kinematics equation, vorticity transport equation and energy equation are solved by the subdomain BEM, while the boundary vorticity values, needed as a boundary conditions for the vorticity transport equation, are calculated by single domain BEM solution of the kinematics equation. Computations are performed for steady state cases, for a range of Darcy numbers from 10^?6 to 10^?1, and porous Rayleigh numbers ranging from 50 to 1000. The heat flux through the cavity and the flow fields are analyzed for different cases of governing parameters and compared to the results in some published studies.     

A numerical contact algorithm in saturated porous media with the extended finite element method

In this paper, a coupled hydro-mechanical formulation is developed for deformable porous media subjected to crack interfaces in the framework of extended finite element method. Governing equations of the porous medium consist of the momentum balance of the bulk together with the momentum balance and continuity equations of the fluid phase, known as formulation. The discontinuity in fractured porous medium is modeled for both opening and closing modes that results in the fluid flow within the fracture, and/or contact behavior at the crack edges. The fluid flow through the fracture is assumed to be viscous and is modeled by employing the Darcy law in which the permeability of fracture is obtained using the cubic law. The contact condition in fractured porous medium is handled by taking the advantage from two different algorithms of LATIN method and penalty algorithm. The effect of contact on fluid phase is employed by considering no leak-off from/into the porous medium. The nonlinearity of coupled equations produced due to opening and closing modes is carried out using an iterative algorithm in the Newton–Raphson procedure. Finally, several numerical examples are solved to illustrate the performance of proposed X-FEM method for hydro-mechanical behavior of fractured porous media with opening and closing modes.     

A stabilized volume‐averaging finite element method for flow in porous media and binary alloy solidification processes

A stabilized equal‐order velocity–pressure finite element algorithm is presented for the analysis of flow in porous media and in the solidification of binary alloys. The adopted governing macroscopic conservation equations of momentum, energy and species transport are derived from their microscopic counterparts using the volume‐averaging method. The analysis is performed in a single domain with a fixed numerical grid. The fluid flow scheme developed includes SUPG (streamline‐upwind/Petrov–Galerkin), PSPG (pressure stabilizing/Petrov–Galerkin) and DSPG (Darcy stabilizing/Petrov–Galerkin) stabilization terms in a variable porosity medium. For the energy and species equations a classical SUPG‐based finite element method is employed. The developed algorithms were tested extensively with bilinear elements and were shown to perform stably and with nearly quadratic convergence in high Rayleigh number flows in varying porosity media. Examples are shown in natural and double diffusive convection in porous media and in the directional solidification of a binary‐alloy. Copyright © 2004 John Wiley & Sons, Ltd.     

Immersed stress method for fluid–structure interaction using anisotropic mesh adaptation

This paper presents advancements toward a monolithic solution procedure and anisotropic mesh adaptation for the numerical solution of fluid–structure interaction with complex geometry. First, a new stabilized three‐field stress, velocity, and pressure finite element formulation is presented for modeling the interaction between the fluid (laminar or turbulent) and the rigid body. The presence of the structure will be taken into account by means of an extra stress in the Navier–Stokes equations. The system is solved using a finite element variational multiscale method. We combine this method with anisotropic mesh adaptation to ensure an accurate capturing of the discontinuities at the fluid–solid interface. We assess the behavior and accuracy of the proposed formulation in the simulation of 2D and 3D time‐dependent numerical examples such as the flow past a circular cylinder and turbulent flows behind an immersed helicopter in a forward flight. Copyright © 2013 John Wiley & Sons, Ltd.     

An iterative BEM/FVM protocol for steady-state multi-dimensional conjugate heat transfer in compressible flows

A BEM-based temperature forward/flux back (TFFB) coupling algorithm is developed to solve the conjugate heat transfer (CHT), which arises naturally in analysis of systems exposed to a convective environment. Here, heat conduction within a structure is coupled to heat transfer to the external fluid, which is convecting heat into or out of the solid structure. There are two basic approaches for solving coupled fluid-structural systems. The first is a direct coupling where the solution of the different fields is solved simultaneously in one large set of equations. The second approach is a loose coupling strategy where each set of field equations is solved to provide boundary conditions for the other. The equations are solved in turn until an iterated convergence criterion is met at the fluid–solid interface. The loose coupling strategy is particularly attractive when coupling auxiliary field equations to computational fluid dynamics codes. We adopt the latter method in which the BEM is used to solve heat conduction inside a structure which is exposed to a convective field which in turn is resolved by solving the Navier–Stokes equations by finite volume methods. Interface of flux and temperature is enforced at the solid/fluid interface.     

Lagrangian finite element analysis of Newtonian fluid flows

A fully Lagrangian finite element method for the analysis of Newtonian flows is developed. The approach furnishes, in effect, a Lagrangian implementation of the compressible Navier–Stokes equations. As the flow proceeds, the mesh is maintained undistorted through continuous and adaptive remeshing of the fluid mass. The principal advantage of the present approach lies in the treatment of boundary conditions at material surfaces such as free boundaries, fluid/fluid or fluid/solid interfaces. In contrast to Eulerian approaches, boundary conditions are enforced at material surfaces ab initio and therefore require no special attention. Consistent tangents are obtained for Lagrangian implicit analysis of a Newtonian fluid flow which may exhibit compressibility effects. The accuracy of the approach is assessed by comparison of the solution for a sloshing problem with existing numerical results and its versatility demonstrated through a simulation of wave breaking. The finite element mesh is maintained undistorted throughout the computation by recourse to frequent and adaptive remeshing © 1998 John Wiley & Sons, Ltd.     

Numerical modelling of void inclusions in porous media

Numerical modelling of porous flow in a low‐permeability matrix with high‐permeability inclusions is a challenging task because the large ratio of permeabilities ill‐conditions the finite element system of equations. We propose a coupled model where Darcy flow is used for the porous matrix and potential flow is used for the inclusions. We discuss appropriate interface conditions in detail and show that the head drop in the inclusions can be prescribed in a very simple way. Algorithmic aspects are treated in full detail. Numerical examples show that this coupled approach precludes ill‐conditioning and is more efficient than heterogeneous Darcy flow. Copyright © 2003 John Wiley & Sons, Ltd.     

Topology optimization for stationary fluid–structure interaction problems using a new monolithic formulation

This paper outlines a new procedure for topology optimization in the steady‐state fluid–structure interaction (FSI) problem. A review of current topology optimization methods highlights the difficulties in alternating between the two distinct sets of governing equations for fluid and structure dynamics (hereafter, the fluid and structural equations, respectively) and in imposing coupling boundary conditions between the separated fluid and solid domains. To overcome these difficulties, we propose an alternative monolithic procedure employing a unified domain rather than separated domains, which is not computationally efficient. In the proposed analysis procedure, the spatial differential operator of the fluid and structural equations for a deformed configuration is transformed into that for an undeformed configuration with the help of the deformation gradient tensor. For the coupling boundary conditions, the divergence of the pressure and the Darcy damping force are inserted into the solid and fluid equations, respectively. The proposed method is validated in several benchmark analysis problems. Topology optimization in the FSI problem is then made possible by interpolating Young's modulus, the fluid pressure of the modified solid equation, and the inverse permeability from the damping force with respect to the design variables. Copyright © 2009 John Wiley & Sons, Ltd.     

A rheological interface model and its space–time finite element formulation for fluid–structure interaction

This contribution discusses extended physical interface models for fluid–structure interaction problems and investigates their phenomenological effects on the behavior of coupled systems by numerical simulation. Besides the various types of friction at the fluid–structure interface the most interesting phenomena are related to effects due to additional interface stiffness and damping. The paper introduces extended models at the fluid–structure interface on the basis of rheological devices (Hooke, Newton, Kelvin, Maxwell, Zener). The interface is decomposed into a Lagrangian layer for the solid‐like part and an Eulerian layer for the fluid‐like part. The mechanical model for fluid–structure interaction is based on the equations of rigid body dynamics for the structural part and the incompressible Navier–Stokes equations for viscous flow. The resulting weighted residual form uses the interface velocity and interface tractions in both layers in addition to the field variables for fluid and structure. The weak formulation of the whole coupled system is discretized using space–time finite elements with a discontinuous Galerkin method for time‐integration leading to a monolithic algebraic system. The deforming fluid domain is taken into account by deformable space–time finite elements and a pseudo‐structure approach for mesh motion. The sensitivity of coupled systems to modification of the interface model and its parameters is investigated by numerical simulation of flow induced vibrations of a spring supported fluid‐immersed cylinder. It is shown that the presented rheological interface model allows to influence flow‐induced vibrations. Copyright © 2010 John Wiley & Sons, Ltd.     

Mixed finite element method for saturated poroelastoplastic media at large strains

A mixed finite element for hydro‐dynamic analysis in saturated porous media in the frame of the Biot theory is proposed. Displacements, effective stresses, strains for the solid phase and pressure, pressure gradients, and Darcy velocities for the fluid phase are interpolated as independent variables. The weak form of the governing equations of coupled hydro‐dynamic problems in saturated porous media within the element are given on the basis of the Hu–Washizu three‐field variational principle. In light of the stabilized one point quadrature super‐convergent element developed in solid continuum, the interpolation approximation modes for the primary unknowns and their spatial derivatives of the solid and the fluid phases within the element are assumed independently. The proposed mixed finite element formulation is derived. The non‐linear version of the element formulation is further derived with particular consideration of pressure‐dependent non‐associated plasticity. The return mapping algorithm for the integration of the rate constitutive equation, the consistent elastoplastic tangent modulus matrix and the element tangent stiffness matrix are developed. For geometrical non‐linearity, the co‐rotational formulation approach is used. Numerical results demonstrate the capability and the performance of the proposed element in modelling progressive failure characterized by strain localization due to strain softening in poroelastoplastic media subjected to dynamic loading at large strain. Copyright © 2003 John Wiley & Sons, Ltd.