A three-dimensional fictitious domain method for incompressible fluid flow problems

A new Galerkin finite element method for the solution of the Navier–Stokes equations in enclosures containing internal parts which may be moving is presented. Dubbed the virtual finite element method, it is based upon optimization techniques and belongs to the class of fictitious domain methods. Only one volumetric mesh representing the enclosure without its internal parts needs to be generated. These are rather discretized using control points on which kinematic constraints are enforced and introduced into the mathematical formulation by means of Lagrange multipliers. Consequently, the meshing of the computational domain is much easier than with classical finite element approaches. First, the methodology will be presented in detail. It will then be validated in the case of the two-dimensional Couette cylinder problem for which an analytical solution is available. Finally, the three-dimensional fluid flow inside a mechanically agitated vessel will be investigated. The accuracy of the numerical results will be assessed through a comparison with experimental data and results obtained with a standard finite element method. © 1997 John Wiley & Sons, Ltd.     

On the distributed Lagrange multiplier/fictitious domain method for rigid‐particle‐laden flows: a proposition for an alternative formulation of the Lagrange multipliers

The distributed Lagrange multiplier/fictitious domain method proposed for the direct numerical simulation of particle‐laden flows is considered in this work. First, it is demonstrated that improved accuracy is obtained with a coupled numerical scheme, whereby the pressure and the Lagrange multiplier fields enforcing incompressibility and rigid body motion, respectively, are calculated and applied together. However, the convergence characteristics of the iterative solution of the coupled scheme are poor because symmetric but indefinite and poorly conditioned matrices are produced. An analysis is then presented, which suggests that the cause for the matrix pathologies lies in the interaction of the respective matrix operators enforcing incompressibility and rigid body motion. On the basis of this analysis, an alternative formulation is developed for the Lagrange multipliers, being now composed of a set of forces distributed only on the particle boundary together with a set of couples distributed within the particle core. The new formulation is tested with several types of flows with stationary or moving particles under creeping or finite Reynolds number conditions and it is demonstrated that it produces correct results and better conditioned matrices, thus enabling faster and more reliable convergence of the conjugate gradient method. The analysis and tests, therefore, support the expectation that the proposed formulation is promising and worthy of further study and improvement. Copyright © 2011 John Wiley & Sons, Ltd.     

A two‐grid fictitious domain method for direct simulation of flows involving non‐interacting particles of a very small size

The full resolution of flows involving particles whose scale is hundreds or thousands of times smaller than the size of the flow domain is a challenging problem. A naive approach would require a tremendous number of degrees of freedom in order to bridge the gap between the two spatial scales involved. The approach used in the present study employs two grids whose grid size fits the two different scales involved, one of them (the micro‐scale grid) being embedded into the other (the macro‐scale grid). Then resolving first the larger scale on the macro‐scale grid, we transfer the so obtained data to the boundary of the micro‐scale grid and solve the smaller size problem. Since the particle is moving throughout the macro‐scale domain, the micro‐scale grid is fixed at the centroid of the moving particle and therefore moves with it. In this study we combine such an approach with a fictitious domain formulation of the problem resulting in a very efficient algorithm that is also easy to implement in an existing CFD code. We validate the method against existing experimental data for a sedimenting sphere, as well as analytical results for motion of an inertia‐less ellipsoid in a shear flow. Finally, we apply the method to the flow of a high aspect ratio ellipsoid in a model of a human lung airway bifurcation. Copyright © 2009 John Wiley & Sons, Ltd.     

Adaptive non‐conformal mesh refinement and extended finite element method for viscous flow inside complex moving geometries

Two different techniques to analyze non‐Newtonian viscous flow in complex geometries with internal moving parts and narrow gaps are compared. The first technique is a non‐conforming mesh refinement approach based on the fictitious domain method (FDM), and the second one is the extended finite element method (XFEM). The refinement technique uses one fixed reference mesh, and to impose continuity across non‐conforming regions, constraints using Lagrangian multipliers are used. The size of elements locally in the high shear rate regions is reduced to increase accuracy. FDM is shown to have limitations; therefore, XFEM is applied to decouple the fluid from the internal moving rigid bodies. In XFEM, the discontinuous field variables are captured by using virtual degrees of freedom that serve as enrichment and by applying special integration over the intersected elements. The accuracy of the two methods is demonstrated by direct comparison with results of a boundary‐fitted mesh applied to a two‐dimensional cross section of a twin‐screw extruder. Compared with non‐conforming FDM, XFEM shows a considerable improvement in accuracy around the rigid body, especially in the narrow gap regions. Copyright © 2011 John Wiley & Sons, Ltd.     

A distributed Lagrange multiplier/fictitious domain method for flows around moving rigid bodies: Application to particulate flow

This article discusses the application of a Lagrange multiplier‐based fictitious domain method to the numerical simulation of incompressible viscous flow modeled by the Navier–Stokes equations around moving rigid bodies; the rigid body motions are due to hydrodynamical forces and gravity. The solution method combines finite element approximations, time discretization by operator splitting and conjugate gradient algorithms for the solution of the linearly constrained quadratic minimization problems coming from the splitting method. The study concludes with the presentation of numerical results concerning four test problems, namely the simulation of an incompressible viscous flow around a NACA0012 airfoil with a fixed center but free to rotate, then the sedimentation of 200 and 1008 cylinders in a two‐dimensional channel, and finally the sedimentation of two spherical balls in a rectangular cylinder. Copyright © 1999 John Wiley & Sons, Ltd.     

A lattice Boltzmann fictitious domain method for modeling red blood cell deformation and multiple‐cell hydrodynamic interactions in flow

To model red blood cell (RBC) deformation and multiple‐cell interactions in flow, the recently developed technique derived from the lattice Boltzmann method and the distributed Lagrange multiplier/fictitious domain method is extended to employ the mesoscopic network model for simulations of RBCs in flow. The flow is simulated by the lattice Boltzmann method with an external force, while the network model is used for modeling RBC deformation. The fluid–RBC interactions are enforced by the Lagrange multiplier. To validate parameters of the RBC network model, stretching tests on both coarse and fine meshes are performed and compared with the corresponding experimental data. Furthermore, RBC deformation in pipe and shear flows is simulated, revealing the capacity of the current method for modeling RBC deformation in various flows. Moreover, hydrodynamic interactions between two RBCs are studied in pipe flow. Numerical results illustrate that the leading cell always has a larger flow velocity and deformation, while the following cells move slower and deform less.Copyright © 2013 John Wiley & Sons, Ltd.     

Solving high Reynolds‐number viscous flows by the general BEM and domain decomposition method

In this paper, the domain decomposition method (DDM) and the general boundary element method (GBEM) are applied to solve the laminar viscous flow in a driven square cavity, governed by the exact Navier–Stokes equations. The convergent numerical results at high Reynolds number Re = 7500 are obtained. We find that the DDM can considerably improve the efficiency of the GBEM, and that the combination of the domain decomposition techniques and the parallel computation can further greatly improve the efficiency of the GBEM. This verifies the great potential of the GBEM for strongly non‐linear problems in science and engineering. Copyright © 2004 John Wiley & Sons, Ltd.     

Semi-analytical finite element method for fictitious crack model in fracture mechanics of concrete

Based on the Hamiltonian governing equations of plane elasticity for sectorial domain, the variable separation and eigenfunction expansion techniques were employed to develop a novel analytical finite element for the fictitious crack model in fracture mechanics of concrete. The new analytical element can be implemented into FEM program systems to solve fictitious crack propagation problems for concrete cracked plates with arbitrary shapes and loads. Numerical results indicate that the method is more efficient and accurate than ordinary finite element method.     

A second‐order time accurate finite element method for quasi‐incompressible viscous flows

A finite element method for quasi‐incompressible viscous flows is presented. An equation for pressure is derived from a second‐order time accurate Taylor–Galerkin procedure that combines the mass and the momentum conservation laws. At each time step, once the pressure has been determined, the velocity field is computed solving discretized equations obtained from another second‐order time accurate scheme and a least‐squares minimization of spatial momentum residuals. The terms that stabilize the finite element method (controlling wiggles and circumventing the Babuska–Brezzi condition) arise naturally from the process, rather than being introduced a priori in the variational formulation. A comparison between the present second‐order accurate method and our previous first‐order accurate formulation is shown. The method is also demonstrated in the computation of the leaky‐lid driven cavity flow and in the simulation of a crossflow past a circular cylinder. In both cases, good agreement with previously published experimental and computational results has been obtained. Copyright © 2010 John Wiley & Sons, Ltd.     

A fictitious domain approach for the Stokes problem based on the extended finite element method

In the present work, we propose to extend to the Stokes problem a fictitious domain approach inspired by extended finite element method and studied for the Poisson problem in a paper of Renard and Haslinger of 2009. The method allows computations in domains whose boundaries do not match. A mixed FEM is used for the fluid flow. The interface between the fluid and the structure is localized by a level‐set function. Dirichlet boundary conditions are taken into account using Lagrange multiplier. A stabilization term is introduced to improve the approximation of the normal trace of the Cauchy stress tensor at the interface and avoid the inf‐sup condition between the spaces for the velocity and the Lagrange multiplier. Convergence analysis is given, and several numerical tests are performed to illustrate the capabilities of the method. Copyright © 2013 John Wiley & Sons, Ltd.     

A domain decomposition method for modelling Stokes flow in porous materials

An algorithm is presented for solving the Stokes equation in large disordered two‐dimensional porous domains. In this work, it is applied to random packings of discs, but the geometry can be essentially arbitrary. The approach includes the subdivision of the domain and a subsequent application of boundary integral equations to the subdomains. This gives a block diagonal matrix with sparse off‐block components that arise from shared variables on internal subdomain boundaries. The global problem is solved using a biconjugate gradient routine with preconditioning. Results show that the effectiveness of the preconditioner is strongly affected by the subdomain structure, from which a methodology is proposed for the domain decomposition step. A minimum is observed in the solution time versus subdomain size, which is governed by the time required for preconditioning, the time for vector multiplications in the biconjugate gradient routine, the iterative convergence rate and issues related to memory allocation. The method is demonstrated on various domains including a random 1000‐particle domain. The solution can be used for efficient recovery of point velocities, which is discussed in the context of stochastic modelling of solute transport. Copyright © 2002 John Wiley & Sons, Ltd.     

An assessment of a parallel,finite element method for three‐dimensional,moving‐boundary flows driven by capillarity for simulation of viscous sintering

A parallel, finite element method is presented for the computation of three‐dimensional, free‐surface flows where surface tension effects are significant. The method employs an unstructured tetrahedral mesh, a front‐tracking arbitrary Lagrangian–Eulerian formulation, and fully implicit time integration. Interior mesh motion is accomplished via pseudo‐solid mesh deformation. Surface tension effects are incorporated directly into the momentum equation boundary conditions using surface identities that circumvent the need to compute second derivatives of the surface shape, resulting in a robust representation of capillary phenomena. Sample results are shown for the viscous sintering of glassy ceramic particles. The most serious performance issue is error arising from mesh distortion when boundary motion is significant. This effect can be severe enough to stop the calculations; some simple strategies for improving performance are tested. Copyright © 2001 John Wiley & Sons, Ltd.     

A fictitious domain/mortar element method for fluid–structure interaction

A new method for the computational analysis of fluid–structure interaction of a Newtonian fluid with slender bodies is developed. It combines ideas of the fictitious domain and the mortar element method by imposing continuity of the velocity field along an interface by means of Lagrange multipliers. The key advantage of the method is that it circumvents the need for complicated mesh movement strategies common in arbitrary Lagrangian–Eulerian (ALE) methods, usually used for this purpose. Copyright © 2001 John Wiley & Sons, Ltd.     

Extended finite element method for viscous flow inside complex three‐dimensional geometries with moving internal boundaries

A three‐dimensional extended finite element method is presented to simulate Stokes flow in complex geometries with internal moving parts. Instead of re‐meshing the flow domain, the kinematics of the internal objects are imposed on the conservation equations using a constraint, implemented with a Lagrangian multiplier. To capture discontinuities of field variables, such as pressure and velocity, on the intersected elements, XFEM is used. To validate our method, it is first applied to a relatively simple problem, that is, the flow around a cylinder in a channel. The results are verified by comparing with a boundary‐fitted solution. After validation of the model and its implementation, the three‐dimensional flow in a twin‐screw extruder is simulated and the results are compared with experimental data from literature. XFEM shows very good accuracy for complex geometries with internal moving parts and narrow gap regions where the shear rate is orders of magnitude higher than in other regions. Copyright © 2011 John Wiley & Sons, Ltd.     

A new modification of the immersed‐boundary method for simulating flows with complex moving boundaries

In this paper, a new immersed‐boundary method for simulating flows over complex immersed, moving boundaries is presented. The flow is computed on a fixed Cartesian mesh and the solid boundaries are allowed to move freely through the mesh. The present method is based on a finite‐difference approach on a staggered mesh together with a fractional‐step method. It must be noted that the immersed boundary is generally not coincident with the position of the solution variables on the grid, therefore, an appropriate strategy is needed to construct a relationship between the curved boundary and the grid points nearby. Furthermore, a momentum forcing is added on the body boundaries and also inside the body to satisfy the no‐slip boundary condition. The immersed boundary is represented by a series of interfacial markers, and the markers are also used as Lagrangian forcing points. A linear interpolation is then used to scale the Lagrangian forcing from the interfacial markers to the corresponding grid points nearby. This treatment of the immersed‐boundary is used to simulate several problems, which have been validated with previous experimental results in the open literature, verifying the accuracy of the present method. Copyright © 2006 John Wiley & Sons, Ltd.     

The solution of viscous incompressible jet flows using non‐staggered boundary fitted co‐ordinate methods

A new approach for the solution of the steady incompressible Navier–Stokes equations in a domain bounded in part by a free surface is presented. The procedure is based on the finite difference technique, with the non‐staggered grid fractional step method used to solve the flow equations written in terms of primitive variables. The physical domain is transformed to a rectangle by means of a numerical mapping technique. In order to design an effective free solution scheme, we distinguish between flows dominated by surface tension and those dominated by inertia and viscosity. When the surface tension effect is insignificant we used the kinematic condition to update the surface; whereas, in the opposite case, we used the normal stress condition to obtain the free surface boundary. Results obtained with the improved boundary conditions for a plane Newtonian jet are found to compare well with the available two‐dimensional numerical solutions for Reynolds numbers, up to Re=100, and Capillary numbers in the range of 0≤Ca<1000. Copyright © 2001 John Wiley & Sons, Ltd.     

A Q2Q1 finite element/level‐set method for simulating two‐phase flows with surface tension

A Q2Q1 (quadratic velocity/linear pressure) finite element/level‐set method was proposed for simulating incompressible two‐phase flows with surface tension. The Navier–Stokes equations were solved using the Q2Q1 integrated FEM, and the level‐set variable was linearly interpolated using a ‘pseudo’ Q2Q1 finite element when calculating the density and viscosity of a fluid to avoid an unbounded density/viscosity. The advection of the level‐set function was calculated through the Taylor–Galerkin method, and the direct approach method is employed for reinitialization. The proposed method was tested by solving several benchmark problems including rising bubbles exhibiting a large density difference and the surface tension effect. The numerical results of the rising bubbles were compared with the existing results to validate the benchmark quantities such as the centroid, circularity, and rising velocity. Furthermore, we focused our attention mainly on mass conservation and time‐step. We observed that the present method represented a convergence rate between 1.0 and 1.5 orders in terms of mass conservation and provided more stable solutions even when using a larger time‐step than the critical time‐step that was imposed because of the explicit treatment of surface tension. Copyright © 2011 John Wiley & Sons, Ltd.     

Computations of a laminar backward‐facing step flow at Re=800 with a spectral domain decomposition method

The two‐dimensional laminar incompressible flow over a backward‐facing step is computed using a spectral domain decomposition approach. A minimum number of subdomains (two) is used; high resolution being achieved by increasing the order of the basis Chebyshev polynomial. Results for the case of a Reynolds number of 800 are presented and compared in detail with benchmark computations. Stable accurate steady flow solutions were obtained using substantially fewer nodes than in previously reported simulations. In addition, the problem of outflow boundary conditions was examined on a shortened domain. Because of their more global nature, spectral methods are particularly sensitive to imposed boundary conditions, which may be exploited in examining the effect of artificial (non‐physical) outflow boundary conditions. Two widely used set of conditions were tested: pseudo stress‐free conditions and zero normal gradient conditions. Contrary to previous results using the finite volume approach, the latter is found to yield a qualitatively erroneous yet stable flow‐field. Copyright © 1999 John Wiley & Sons, Ltd.