An arbitrary Lagrangian–Eulerian finite element approach to non‐steady state turbulent fluid flow with application to mould filling in casting

This paper presents a two‐dimensional Lagrangian–Eulerian finite element approach of non‐steady state turbulent fluid flows with free surfaces. The proposed model is based on a velocity–pressure finite element Navier–Stokes solver, including an augmented Lagrangian technique and an iterative resolution of Uzawa type. Turbulent effects are taken into account with the k–ε two‐equation statistical model. Mesh updating is carried out through an arbitrary Lagrangian–Eulerian (ALE) method in order to describe properly the free surface evolution. Three comparisons between experimental and numerical results illustrate the efficiency of the method. The first one is turbulent flow in an academic geometry, the second one is a mould filling in effective casting conditions and the third one is a precise confrontation to a water model. Copyright © 2000 John Wiley & Sons, Ltd.     

Three-dimensional arbitrary Lagrangian–Eulerian numerical prediction method for non-linear free surface oscillation

A numerical prediction method has been proposed to predict non-linear free surface oscillation in an arbitrarily-shaped three-dimensional container. The liquid motions are described with Navier–Stokes equations rather than Laplace equations which are derived by assuming the velocity potential. The profile of a liquid surface is precisely represented with the three-dimensional curvilinear co-ordinates which are regenerated in each computational step on the basis of the arbitrary Lagrangian–Eulerian (ALE) formulation. In the transformed space, the governing equations are discretized on a Lagrangian scheme with sufficient numerical accuracy and the boundary conditions near the liquid surface are implemented in a complete manner. In order to confirm the applicability of the present computational technique, numerical simulations are conducted for the free oscillations of viscid and inviscid liquids and for highly non-linear oscillation. In addition, non-linear sloshing motions caused by horizontal and vertical excitations and a transition from non-linear sloshing to swirling are numerically predicted in three-dimensional cylindrical containers. Conclusively, it is shown that these sloshing motions associated with high non-linearity are reasonably predicted with the present numerical technique. © 1998 John Wiley & Sons, Ltd.     

A Lagrangian–Eulerian model of particle dispersion in a turbulent plane mixing layer

A Lagrangian–Eulerian model for the dispersion of solid particles in a two‐dimensional, incompressible, turbulent flow is reported and validated. Prediction of the continuous phase is done by solving an Eulerian model using a control‐volume finite element method (CVFEM). A Lagrangian model is also applied, using a Runge–Kutta method to obtain the particle trajectories. The effect of fluid turbulence upon particle dispersion is taken into consideration through a simple stochastic approach. Validation tests are performed by comparing predictions for both phases in a particle‐laden, plane mixing layer airflow with corresponding measurements formerly reported by other authors. Even though some limitations are detected in the calculation of particle dispersion, on the whole the validation results are rather successful. Copyright © 2002 John Wiley & Sons, Ltd.     

Developing a unified FVE‐ALE approach to solve unsteady fluid flow with moving boundaries

In this study, an arbitrary Lagrangian–Eulerian (ALE) approach is incorporated with a mixed finite‐volume–element (FVE) method to establish a novel moving boundary method for simulating unsteady incompressible flow on non‐stationary meshes. The method collects the advantages of both finite‐volume and finite‐element (FE) methods as well as the ALE approach in a unified algorithm. In this regard, the convection terms are treated at the cell faces using a physical‐influence upwinding scheme, while the diffusion terms are treated using bilinear FE shape functions. On the other hand, the performance of ALE approach is improved by using the Laplace method to improve the hybrid grids, involving triangular and quadrilateral elements, either partially or entirely. The use of hybrid FE grids facilitates this achievement. To show the robustness of the unified algorithm, we examine both the first‐ and the second‐order temporal stencils. The accuracy and performance of the extended method are evaluated via simulating the unsteady flow fields around a fixed cylinder, a transversely oscillating cylinder, and in a channel with an indented wall. The numerical results presented demonstrate significant accuracy benefits for the new hybrid method on coarse meshes and where large time steps are taken. Of importance, the current method yields the second‐order temporal accuracy when the second‐order stencil is used to discretize the unsteady terms. Copyright © 2009 John Wiley & Sons, Ltd.     

Simulation of liquid sloshing in curved‐wall containers with arbitrary Lagrangian–Eulerian method

There are many challenges in the numerical simulation of liquid sloshing in horizontal cylinders and spherical containers using the finite element method of arbitrary Lagrangian–Eulerian (ALE) formulation: tracking the motion of the free surface with the contact points, defining the mesh velocity on the curved wall boundary and updating the computational mesh. In order to keep the contact points slipping along the curved side wall, the shape vector in each time advancement is defined to modify the kinematical boundary conditions on the free surface. A special function is introduced to automatically smooth the nodal velocities on the curved wall boundary based on the liquid nodal velocities. The elliptic partial differential equation with Dirichlet boundary conditions can directly rezone the inner nodal velocities in more than a single freedom. The incremental fractional step method is introduced to solve the finite element liquid equations. The numerical results that stemmed from the algorithm show good agreement with experimental phenomena, which demonstrates that the ALE method provides an efficient computing scheme in moving curved wall boundaries. This method can be extended to 3D cases by improving the technique to compute the shape vector. Copyright © 2007 John Wiley & Sons, Ltd.     

High‐order implicit Runge–Kutta time integrators for fluid‐structure interactions

This paper presents an approach to develop high‐order, temporally accurate, finite element approximations of fluid‐structure interaction (FSI) problems. The proposed numerical method uses an implicit monolithic formulation in which the same implicit Runge–Kutta (IRK) temporal integrator is used for the incompressible flow, the structural equations undergoing large displacements, and the coupling terms at the fluid‐solid interface. In this context of stiff interaction problems, the fully implicit one‐step approach presented is an original alternative to traditional multistep or explicit one‐step finite element approaches. The numerical scheme takes advantage of an arbitrary Lagrangian–Eulerian formulation of the equations designed to satisfy the geometric conservation law and to guarantee that the high‐order temporal accuracy of the IRK time integrators observed on fixed meshes is preserved on arbitrary Lagrangian–Eulerian deforming meshes. A thorough review of the literature reveals that in most previous works, high‐order time accuracy (higher than second order) is seldom achieved for FSI problems. We present thorough time‐step refinement studies for a rigid oscillating‐airfoil on deforming meshes to confirm the time accuracy on the extracted aerodynamics reactions of IRK time integrators up to fifth order. Efficiency of the proposed approach is then tested on a stiff FSI problem of flow‐induced vibrations of a flexible strip. The time‐step refinement studies indicate the following: stability of the proposed approach is always observed even with large time step and spurious oscillations on the structure are avoided without added damping. While higher order IRK schemes require more memory than classical schemes (implicit Euler), they are faster for a given level of temporal accuracy in two dimensions. Copyright © 2015 John Wiley & Sons, Ltd.     

High‐order time integrators for front‐tracking finite‐element analysis of viscous free‐surface flows

This paper proposes implicit Runge–Kutta (IRK) time integrators to improve the accuracy of a front‐tracking finite‐element method for viscous free‐surface flow predictions. In the front‐tracking approach, the modeling equations must be solved on a moving domain, which is usually performed using an arbitrary Lagrangian–Eulerian (ALE) frame of reference. One of the main difficulties associated with the ALE formulation is related to the accuracy of the time integration procedure. Indeed, most formulations reported in the literature are limited to second‐order accurate time integrators at best. In this paper, we present a finite‐element ALE formulation in which a consistent evaluation of the mesh velocity and its divergence guarantees satisfaction of the discrete geometrical conservation law. More importantly, it also ensures that the high‐order fixed mesh temporal accuracy of time integrators is preserved on deforming grids. It is combined with the use of a family of L‐stable IRK time integrators for the incompressible Navier–Stokes equations to yield high‐order time‐accurate free‐surface simulations. This is demonstrated in the paper using the method of manufactured solution in space and time as recommended in Verification and Validation. In particular, we report up to fifth‐order accuracy in time. The proposed free‐surface front‐tracking approach is then validated against cases of practical interest such as sloshing in a tank, solitary waves propagation, and coupled interaction between a wave and a submerged cylinder. Copyright © 2015 John Wiley & Sons, Ltd.     

A projection scheme for incompressible multiphase flow using adaptive Eulerian grid: 3D validation

A three‐dimensional finite element method for incompressible multiphase flows with capillary interfaces is developed based on a (formally) second‐order projection scheme. The discretization is on a fixed (Eulerian) reference grid with an edge‐based local h‐refinement in the neighbourhood of the interfaces. The fluid phases are identified and advected using the level‐set function. The reference grid is then temporarily reconnected around the interface to maintain optimal interpolations accounting for the singularities of the primary variables. Using a time splitting procedure, the convection substep is integrated with an explicit scheme. The remaining generalized Stokes problem is solved by means of a pressure‐stabilized projection. This method is simple and efficient, as demonstrated by a wide range of difficult free‐surface validation problems, considered in the paper. Copyright © 2005 John Wiley & Sons, Ltd.     

Handling contacts in an Eulerian frame: a finite element approach for fluid structures with contacts

ABSTRACT

Eulerian variational formulations for deformable solids, with or without fluids around them, end up, after implicit time discretisation, as large non-linear systems for the velocities in the moving domains. Handling moving domains and moving boundaries requires careful meshing procedures; on the other hand, the detection of contact is particularly simple with a distance function. Then at every time step, a variational inequality can be used to update the velocities. This article gives new implementation details and two new complex simulations: a very soft bouncing ball in an axisymmetric flow and a disk hit by a club.     

Accurate numerical estimation of interphase momentum transfer in Lagrangian–Eulerian simulations of dispersed two-phase flows

The Lagrangian–Eulerian (LE) approach is used in many computational methods to simulate two-way coupled dispersed two-phase flows. These include averaged equation solvers, as well as direct numerical simulations (DNS) and large-eddy simulations (LES) that approximate the dispersed-phase particles (or droplets or bubbles) as point sources. Accurate calculation of the interphase momentum transfer term in LE simulations is crucial for predicting qualitatively correct physical behavior, as well as for quantitative comparison with experiments. Numerical error in the interphase momentum transfer calculation arises from both forward interpolation/approximation of fluid velocity at grid nodes to particle locations, and from backward estimation of the interphase momentum transfer term at particle locations to grid nodes. A novel test that admits an analytical form for the interphase momentum transfer term is devised to test the accuracy of the following numerical schemes: (1) fourth-order Lagrange Polynomial Interpolation (LPI-4), (3) Piecewise Cubic Approximation (PCA), (3) second-order Lagrange Polynomial Interpolation (LPI-2) which is basically linear interpolation, and (4) a Two-Stage Estimation algorithm (TSE). A number of tests are performed to systematically characterize the effects of varying the particle velocity variance, the distribution of particle positions, and fluid velocity field spectrum on estimation of the mean interphase momentum transfer term. Numerical error resulting from backward estimation is decomposed into statistical and deterministic (bias and discretization) components, and their convergence with number of particles and grid resolution is characterized. It is found that when the interphase momentum transfer is computed using values for these numerical parameters typically encountered in the literature, it can incur errors as high as 80% for the LPI-4 scheme, whereas TSE incurs a maximum error of 20%. The tests reveal that using multiple independent simulations and higher number of particles per cell are required for accurate estimation using current algorithms. The study motivates further testing of LE numerical methods, and the development of better algorithms for computing interphase transfer terms.     

Flow simulation on moving boundary‐fitted grids and application to fluid–structure interaction problems

We present a method for the parallel numerical simulation of transient three‐dimensional fluid–structure interaction problems. Here, we consider the interaction of incompressible flow in the fluid domain and linear elastic deformation in the solid domain. The coupled problem is tackled by an approach based on the classical alternating Schwarz method with non‐overlapping subdomains, the subproblems are solved alternatingly and the coupling conditions are realized via the exchange of boundary conditions. The elasticity problem is solved by a standard linear finite element method. A main issue is that the flow solver has to be able to handle time‐dependent domains. To this end, we present a technique to solve the incompressible Navier–Stokes equation in three‐dimensional domains with moving boundaries. This numerical method is a generalization of a finite volume discretization using curvilinear coordinates to time‐dependent coordinate transformations. It corresponds to a discretization of the arbitrary Lagrangian–Eulerian formulation of the Navier–Stokes equations. Here the grid velocity is treated in such a way that the so‐called Geometric Conservation Law is implicitly satisfied. Altogether, our approach results in a scheme which is an extension of the well‐known MAC‐method to a staggered mesh in moving boundary‐fitted coordinates which uses grid‐dependent velocity components as the primary variables. To validate our method, we present some numerical results which show that second‐order convergence in space is obtained on moving grids. Finally, we give the results of a fully coupled fluid–structure interaction problem. It turns out that already a simple explicit coupling with one iteration of the Schwarz method, i.e. one solution of the fluid problem and one solution of the elasticity problem per time step, yields a convergent, simple, yet efficient overall method for fluid–structure interaction problems. Copyright © 2005 John Wiley & Sons, Ltd.     

A meshless finite point method for three‐dimensional analysis of compressible flow problems involving moving boundaries and adaptivity

A finite point method for solving compressible flow problems involving moving boundaries and adaptivity is presented. The numerical methodology is based on an upwind‐biased discretization of the Euler equations, written in arbitrary Lagrangian–Eulerian form and integrated in time by means of a dual‐time steeping technique. In order to exploit the meshless potential of the method, a domain deformation approach based on the spring network analogy is implemented, and h‐adaptivity is also employed in the computations. Typical movable boundary problems in transonic flow regime are solved to assess the performance of the proposed technique. In addition, an application to a fluid–structure interaction problem involving static aeroelasticity illustrates the capability of the method to deal with practical engineering analyses. The computational cost and multi‐core performance of the proposed technique is also discussed through the examples provided. Copyright © 2013 John Wiley & Sons, Ltd.     

Monotonicity in high‐order curvilinear finite element arbitrary Lagrangian–Eulerian remap

The remap phase in arbitrary Lagrangian–Eulerian (ALE) hydrodynamics involves the transfer of field quantities defined on a post‐Lagrangian mesh to some new mesh, usually generated by a mesh optimization algorithm. This problem is often posed in terms of transporting (or advecting) some state variable from the old mesh to the new mesh over a fictitious time interval. It is imperative that this remap process be monotonic, that is, not generate any new extrema in the field variables. It is well known that the only linear methods that are guaranteed to be monotonic for such problems are first‐order accurate; however, much work has been performed in developing non‐linear methods, which blend both high and low (first) order solutions to achieve monotonicity and preserve high‐order accuracy when the field is sufficiently smooth. In this paper, we present a set of methods for enforcing monotonicity targeting high‐order discontinuous Galerkin methods for advection equations in the context of high‐order curvilinear ALE hydrodynamics. Published 2014. This article is a U.S. Government work and is in the public domain in the USA.     

Prolonged simulation of near-free surface underwater explosion based on Eulerian finite element method

In the area of naval architecture and ocean engineering, the research about the underwater explosion problem is of great significance. To achieve prolonged simulation of near-free surface underwater explosion, the underwater explosion transient numerical model is established in this paper based on compressible Eulerian finite element method(EFEM). Compared with Geers–Hunter formula, EFEM is availably validated by simulating the free-field underwater explosion case. Then, the bubble pulsation and flow field dynamic characteristics of the cases with different underwater explosive depth are compared in this work. Lastly, the height of the water hump and the pressure of flow flied are analyzed quantitatively through the simulation results.     

Two‐dimensional non‐Newtonian injection molding with a new control volume FEM/volume of fluid method

A new method based on volume of fluid for interface tracking in the simulation of injection molding is presented. The proposed method is comprised of two main stages: accumulation and distribution of the volume fraction. In the first stage the equation for the volume fraction with a noninterfacial flux condition is solved. In the second stage the accumulated volume of fluid that arises as a consequence of the application of the first one is dispersed. This procedure guarantees that the fluid fills the available space without dispersion of the interface. The mathematical model is based on two‐phase transport equations that are numerically integrated through the control volume finite element method. The numerical results for the interface position are successfully verified with analytical results and numerical data available in the literature for one‐dimensional and two‐dimensional domains. The transient position of the advance fronts showed an effective and consistent simulation of an injection molding process. The nondispersive volume of fluid method here proposed is implemented for the simulation of nonisothermal injection molding in two‐dimensional cavities. The obtained results are represented as transient interface positions, isotherms and pressure distributions during the injection molding of low density polyethylene. Copyright © 2012 John Wiley & Sons, Ltd.     

Quantitative comparison of Taylor flow simulations based on sharp‐interface and diffuse‐interface models

A pressure‐driven flow of elongated bullet‐shaped bubbles in a narrow channel is known as Taylor flow or bubble‐train flow. This process is of relevance in various applications of chemical engineering. In this paper, we describe a typical simplified experimental setting, with surface tension, density and viscosity as prescribed input parameters. We compare a sharp‐interface model based on a moving grid aligned with the bubble boundary (ALE coordinates) and a diffuse‐interface model where the bubble shape is implicitly given by a phase‐field function. Four independent implementations based on the two modeling approaches are introduced and described briefly. Besides the simulation of the bubble shapes, we compare some resulting quantities such as pressure difference and film widths within the implementations and to existing analytical and experimental results. The simulations were conducted in 2D and 3D (rotationally symmetric). Good accordance of the results indicate the applicability and the usability of all approaches. Differences between the models and their implementations are visible but in no contradiction to theoretical results. Copyright © 2013 John Wiley & Sons, Ltd.     

An accurate and efficient semi‐implicit method for section‐averaged free‐surface flow modelling

An accurate, efficient and robust numerical method for the solution of the section‐averaged De St. Venant equations of open channel flow is presented and discussed. The method consists in a semi‐implicit, finite‐volume discretization of the continuity equation capable to deal with arbitrary cross‐section geometry and in a semi‐implicit, finite‐difference discretization of the momentum equation. By using a proper semi‐Lagrangian discretization of the momentum equation, a highly efficient scheme that is particularly suitable for subcritical regimes is derived. Accurate solutions are obtained in all regimes, except in presence of strong unsteady shocks as in dam‐break cases. By using a suitable upwind, Eulerian discretization of the same equation, instead, a scheme capable of describing accurately also unsteady shocks can be obtained, although this scheme requires to comply with a more restrictive stability condition. The formulation of the two approaches allows a unified implementation and an easy switch between the two. The code is verified in a wide range of idealized test cases, highlighting its accuracy and efficiency characteristics, especially for long time range simulations of subcritical river flow. Finally, a model validation on field data is presented, concerning simulations of a flooding event of the Adige river. Copyright © 2009 John Wiley & Sons, Ltd.