A finite element method for free surface flows of incompressible fluids in three dimensions. Part II. Dynamic wetting lines

To date, few researchers have solved three‐dimensional free surface problems with dynamic wetting lines. This paper extends the free surface finite element method (FEM) described in a companion paper [Cairncross RA, Schunk PR, Baer TA, Sackinger PA, Rao RR. A finite element method for free surface flows of incompressible fluid in three dimensions. Part I. Boundary fitted mesh motion. International Journal for Numerical Methods in Fluids 2000; 33 : 375–403] to handle dynamic wetting. A generalization of the technique used in two‐dimensional modeling to circumvent double‐valued velocities at the wetting line, the so‐called kinematic paradox, is presented for a wetting line in three dimensions. This approach requires the fluid velocity normal to the contact line to be zero, the fluid velocity tangent to the contact line to be equal to the tangential component of web velocity, and the fluid velocity into the web to be zero. In addition, slip is allowed in a narrow strip along the substrate surface near the dynamic contact line. For realistic wetting line motion, a contact angle that varies with wetting speed is required because contact lines in three dimensions typically advance or recede at different rates depending upon location and/or have both advancing and receding portions. The theory is applied to capillary rise of static fluid in a corner, the initial motion of a Newtonian droplet down an inclined plane, and extrusion of a Newtonian fluid from a nozzle onto a moving substrate. The extrusion results are compared with experimental visualization. Copyright © 2000 John Wiley & Sons, Ltd.     

An adaptive Lagrangian boundary element approach for three‐dimensional transient free‐surface Stokes flow as applied to extrusion,thermoforming, and rheometry

An adaptive (Lagrangian) boundary element approach is proposed for the general three‐dimensional simulation of confined free‐surface Stokes flow. The method is stable as it includes remeshing capabilities of the deforming free surface and thus can handle large deformations. A simple algorithm is developed for mesh refinement of the deforming free‐surface mesh. Smooth transition between large and small elements is achieved without significant degradation of the aspect ratio of the elements in the mesh. Several flow problems are presented to illustrate the utility of the approach, particularly as encountered in polymer processing and rheology. These problems illustrate the transient nature of the flow during the processes of extrusion and thermoforming, the elongation of a fluid sample in an extensional rheometer, and the coating of a sphere. Surface tension effects are also explored. Copyright © 2001 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.     

Numerical simulation of liquid crystalline polymer flow into two‐dimensional thin cavity moulds

In this paper, flows of liquid crystalline polymers into two‐dimensional thin cavity moulds are simulated. The flows are modelled by Ericksen–Leslie equations of motion in the high viscosity limit. An elliptic pressure equation is derived under Hele–Shaw approximations, and the non‐isothermal natures of the flow are modelled. The equations are solved using the finite‐difference technique. A new boundary‐mapping technique is developed in this study to solve the difficulty in the finite‐difference treatment of arbitrarily shaped boundaries, which possess no natural coordinate system. This new method avoids the difficult mesh control in the body‐fitted mapping process and makes the mapping process easy to implement. It can also solve the problems caused by the uneven distribution of grid nodes in the traditional body‐fitted mapping technique. Copyright © 2008 John Wiley & Sons, Ltd.     

Three‐dimensional boundary element analysis of drop deformation in confined flow for Newtonian and viscoelastic systems

An adaptive (Lagrangian) boundary element approach is proposed for the general three‐dimensional drop deformation in confined flow. The adaptive method is stable as it includes remeshing capabilities of the deforming interface between drop and suspending fluid, and thus can handle large deformations. Both drop and surrounding fluid are viscous incompressible and can be Newtonian or viscoelastic. A boundary‐only formulation is implemented for fluids obeying the linear Jeffrey's constitutive equation. Similarly to the formulation for two‐dimensional Newtonian fluids (Khayat RE, Luciani A, Utracki LA. Boundary element analysis of planar drop deformation in confined flow. Part I. Newtonian fluids. Engineering Analysis of Boundary Elements 1997; 19 : 279), the method requires the solution of two simultaneous integral equations on the interface between the two fluids and the confining solid boundary. Although the problem is formulated for any confining geometry, the method is illustrated for a deforming drop as it is driven by the ambient flow inside a cylindrical tube. The accuracy of the method is assessed by comparison with the analytical solution for two‐phase radial spherical flow, leading to good agreement. The influence of mesh refinement is examined for a drop in simple shear flow. Copyright © 2000 John Wiley & Sons, Ltd.     

A fully coupled Navier–Stokes solver for calculation of turbulent incompressible free surface flow past a ship hull

This paper deals with the calculation of free surface flow of viscous incompressible fluid around the hull of a boat moving with rectilinear motion. An original method used to avoid a large part of the theoretical problems connected with free surface boundary conditions in three‐dimensional Navier–Stokes–Reynolds equations is proposed here. The linearised system of convective equations for velocities, pressure and free surface elevation unknowns is discretised by finite differences and two methods to solve the fully coupled resulting matrix are presented here. The non‐linear convergence of fully coupled algorithm is compared with the velocity–pressure weakly coupled algorithm SIMPLER. Turbulence is taken into account through Reynolds decomposition and k–ε or k–ω model to close the equations. These two models are implemented without wall function and numerical calculations are performed up to the viscous sub‐layer. Numerical results and comparisons with experiments are presented on the Series 60 CB=0.60 ship model for a Reynolds number Rn=4.5×10⁶ and a Froude number Fn=0.316. Copyright © 1999 John Wiley & Sons, Ltd.     

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.     

Three‐dimensional transient free‐surface flow of viscous fluids inside cavities of arbitrary shape

The three‐dimensional transient free‐surface flow inside cavities of arbitrary shape is examined in this study. An adaptive (Lagrangian) boundary‐element approach is proposed for the general three‐dimensional simulation of confined free‐surface flow of viscous incompressible fluids. The method is stable as it includes remeshing capabilities of the deforming free‐surface, and thus can handle large deformations. A simple algorithm is developed for mesh refinement of the deforming free‐surface mesh. Smooth transition between large and small elements is achieved without significant degradation of the aspect ratio of the elements in the mesh. The method is used to determine the flow field and free‐surface evolution inside cubic, rectangular and cylindrical containers. These problems illustrate the transient nature of the flow during the mixing process. Surface tension effects are also explored. Copyright © 2003 John Wiley & Sons, Ltd.     

Incompressible SPH simulation of water entry of a free‐falling object

In this paper an incompressible smoothed particle hydrodynamics (Incom‐SPH) model is used to simulate the interactions between the free surface flow and a moving object. Incom‐SPH method is a two‐step semi‐implicit hydrodynamic formulation of the SPH algorithm and is capable of accurately treating the free surface deformations and impact forces during the solid–fluid interactions. For a free‐falling object, its motion is tracked by an additional Lagrangian algorithm based on Newton's law to couple with the Incom‐SPH program. The developed model is employed to investigate the water entry of a free‐falling wedge. The accuracy of the computations is validated by the good agreement in predicting the relevant hydrokinematic and hydrodynamic parameters. Finally, a numerical test is performed to study the influence of spatial resolution on the water entry features. The Incom‐SPH modeling coupled with the solid–fluid interaction algorithm could provide a promising computational tool to predict the slamming problems in coastal and offshore engineering. Copyright © 2008 John Wiley & Sons, Ltd.     

A new fully non‐hydrostatic 3D free surface flow model for water wave motions

A new fully non‐hydrostatic model is presented by simulating three‐dimensional free surface flow on a vertical boundary‐fitted coordinate system. A projection method, known as pressure correction technique, is employed to solve the incompressible Euler equations. A new grid arrangement is proposed under a horizontal Cartesian grid framework and vertical boundary‐fitted coordinate system. The resulting model is relatively simple. Moreover, the discretized Poisson equation for pressure correction is symmetric and positive definite, and thus it can be solved effectively by the preconditioned conjugate gradient method. Several test cases of surface wave motion are used to demonstrate the capabilities and numerical stability of the model. Comparisons between numerical results and analytical or experimental data are presented. It is shown that the proposed model could accurately and effectively resolve the motion of short waves with only two layers, where wave shoaling, nonlinearity, dispersion, refraction, and diffraction phenomena occur. Copyright © 2010 John Wiley & Sons, Ltd.     

Algorithms for interface treatment and load computation in embedded boundary methods for fluid and fluid–structure interaction problems

Embedded boundary methods for CFD (computational fluid dynamics) simplify a number of issues. These range from meshing the fluid domain, to designing and implementing Eulerian‐based algorithms for fluid–structure applications featuring large structural motions and/or deformations. Unfortunately, embedded boundary methods also complicate other issues such as the treatment of the wall boundary conditions in general, and fluid–structure transmission conditions in particular. This paper focuses on this aspect of the problem in the context of compressible flows, the finite volume method for the fluid, and the finite element method for the structure. First, it presents a numerical method for treating simultaneously the fluid pressure and velocity conditions on static and dynamic embedded interfaces. This method is based on the exact solution of local, one‐dimensional, fluid–structure Riemann problems. Next, it describes two consistent and conservative approaches for computing the flow‐induced loads on rigid and flexible embedded structures. The first approach reconstructs the interfaces within the CFD solver. The second one represents them as zero level sets, and works instead with surrogate fluid/structure interfaces. For example, the surrogate interfaces obtained simply by joining contiguous segments of the boundary surfaces of the fluid control volumes that are the closest to the zero level sets are explored in this work. All numerical algorithms presented in this paper are applicable with any embedding CFD mesh, whether it is structured or unstructured. Their performance is illustrated by their application to the solution of three‐dimensional fluid–structure interaction problems associated with the fields of aeronautics and underwater implosion. Copyright © 2011 John Wiley & Sons, Ltd.     

RKDG有限元方法计算流体与刚体耦合

用Level set方法配合Runge-Kutta discontinuous Galerkin (RKDG)有限元方法求解流体与刚体耦合问题。用RKDG有限元方法求解欧拉方程,通过求解Level set方程对界面进行追踪,并用推广的Ghost fluid方法对流刚界面进行处理。数值实验表明,该方法具有较高的分辨率。由于该方法不需要对移动网格进行处理,因此可以处理任意形状的拓扑问题,并且很容易推广到三维。     

Solution of incompressible flows with or without a free surface using the finite volume method on unstructured triangular meshes

An incompressible Navier–Stokes solver based on a cell‐centre finite volume formulation for unstructured triangular meshes is developed and tested. The solution methodology makes use of pseudocompressibility, whereby the convective terms are computed using a Godunov‐type second‐order upwind finite volume formulation. The evolution of the solution in time is obtained by subiterating the equations in pseudotime for each physical time step, with the pseudotime step set equal to infinity. For flows with a free surface the computational mesh is fitted to the free surface boundary at each time step, with the free surface elevation satisfying a kinematic boundary condition. A ‘leakage coefficient’, ε, is introduced for the calculation of flows with a free surface in order to control the leakage of flow through the free surface. This allows the assumption of stationarity of mesh points to be made during the course of pseudotime iteration. The solver is tested by comparing the output with a wide range of documented published results, both for flows with and without a free surface. The presented results show that the solver is robust. © 1999 John Wiley & Sons, Ltd.     

An immersed boundary method for unstructured meshes in depth averaged shallow water models

The representation of geometries as buildings, flood barriers or dikes in free surface flow models implies tedious and time‐consuming operations in order to define accurately the shape of these objects when using a body fitted numerical mesh. The immersed boundary method is an alternative way to define solid bodies inside the computational domain without the need of fitting the mesh boundaries to the shape of the object. In the direct forcing immersed boundary method, a solid body is represented by a grid of Lagrangian markers, which define its shape and which are independent from the fluid Eulerian mesh. This paper presents a new implementation of the immersed boundary method in an unstructured finite volume solver for the 2D shallow water equations. Moving least‐squares is used to transmit information between the grid of Lagrangian markers and the fluid Eulerian mesh. The performance of the proposed implementation is analysed in three test cases involving different flow conditions: the flow around a spur dike, a dam break flow with an isolated obstacle and the flow around an array of obstacles. A very good agreement between the classic body fitted approach and the immersed boundary method was found. The differences between the results obtained with both methods are less relevant than the errors because of the intrinsic shallow water assumptions. Copyright © 2015 John Wiley & Sons, Ltd.     

Implicit Calculations of an Aeroelasticity Problem

The purpose of this work is to show that a linearized implicit scheme for the flow resolution can be an efficient and accurate method for solving fluid-structure interaction. The fluid is modeled by the Euler equations in two dimensions and the structure by a one (free piston) or a two (NACA0012 airfoil) degrees of freedom system. The schemes are developed using a finite volume/finite element formulation and, stating the moving boundary problem in the space-time domain, the Riemann solver is generalized in a suitable manner. Assuming a modal decomposition for the structure's response, an analytical solution to the equation of motion is obtained.

The effects of the linearized implicit scheme on the aeroelastic response are demonstrated on the free piston and the NACA 0012 airfoil problems. In the latter case, we focus on the capability of the linearized implicit scheme to accurately predict the stability limit of the coupled response (wing flutter analysis). Although the above analysis is performed using a rigid transformation, a robust moving mesh strategy is presented for more general 2-D and 3-D deformations.     

A coupled boundary element–finite difference model of surface wave motion over a wall turbulent flow

An effective numerical technique is presented to model turbulent motion of a standing surface wave in a tank. The equations of motion for turbulent boundary layers at the solid surfaces are coupled with the potential flow in the bulk of the fluid, and a mixed BEM–finite difference technique is used to model the wave motion and the corresponding boundary layer flow. A mixing‐length theory is used for turbulence modelling. The model results are in good agreement with previous physical and numerical experiments. Although the technique is presented for a standing surface wave, it can be easily applied to other free surface problems. Copyright © 2005 John Wiley & Sons, Ltd.     

SPH for 3D floating bodies using variable mass particle distribution

A floating body with substantial heave motion is a challenging fluid–structure interaction problem for numerical simulation. In this paper we develop SPH in three dimensions to include variable particle mass distribution using an arbitrary Lagrange–Eulerian formulation with an embedded Riemann solver. A wedge or cone in initially still water is forced to move with a displacement equal to the surface elevation of a focused wave group. A two‐dimensional wedge case is used to evaluate two forms of repulsive‐force boundary condition on the body; the force depending on the normal distance from the object surface produced closer agreement with the experiment. For a three‐dimensional heaving cone the comparison between SPH and experiment shows excellent agreement for the force and free surface for motion with low peak spectral frequencies while for a higher peak frequency the agreement is reasonable in terms of phase and magnitude, but a small discrepancy appears at the troughs in the motion. Capturing the entire three‐dimensional flow field using an initially uniform particle distribution with sufficiently fine resolution requires an extremely large number of particles and consequently large computing resource. To mitigate this issue, we employ a variable mass distribution with fine resolution around the body. Using a refined mass distribution in a preselected area avoids the need for a dynamic particle refinement scheme and leads to a computational speedup of more than 600% or much improved results for a given number of particles. SPH with variable mass distribution is then applied to a single heaving‐float wave energy converter, the ‘Manchester Bobber’, in extreme waves and compared with experiments in a wave tank. The SPH simulations are presented for two cases: a single degree‐of‐freedom system with motion restricted to the vertical direction and with general motion allowing six degrees‐of‐freedom. The motion predicted for the float with general motion is in much closer agreement with experimental data than the vertically constrained system. Using variable particle mass distribution is shown to produce close agreement with a computation time 20% of that required with a uniformly fine resolution. Copyright © 2012 John Wiley & Sons, Ltd.     

An explicit formulation for the evolution of nonlinear surface waves interacting with a submerged body

An explicit formulation to study nonlinear waves interacting with a submerged body in an ideal fluid of infinite depth is presented. The formulation allows one to decompose the nonlinear wave–body interaction problem into body and free‐surface problems. After the decomposition, the body problem satisfies a modified body boundary condition in an unbounded fluid domain, while the free‐surface problem satisfies modified nonlinear free‐surface boundary conditions. It is then shown that the nonlinear free‐surface problem can be further reduced to a closed system of two nonlinear evolution equations expanded in infinite series for the free‐surface elevation and the velocity potential at the free surface. For numerical experiments, the body problem is solved using a distribution of singularities along the body surface and the system of evolution equations, truncated at third order in wave steepness, is then solved using a pseudo‐spectral method based on the fast Fourier transform. A circular cylinder translating steadily near the free surface is considered and it is found that our numerical solutions show excellent agreement with the fully nonlinear solution using a boundary integral method. We further validate our solutions for a submerged circular cylinder oscillating vertically or fixed under incoming nonlinear waves with other analytical and numerical results. Copyright © 2007 John Wiley & Sons, Ltd.     

A variationally consistent finite element approach to the two‐fluid internal contact‐line problem

A new method for simulating incompressible viscous fluid flow involving moving internal contact lines is presented. The steady state interface shape is determined by a variationally consistent formulation of the surface tension contribution to the equations of motion adapted to the case of internal contact lines through the application of a global force balance compatibility condition that consistently removes the pressure indeterminacy. The Crouzeix–Raviart element is chosen to capture the pressure discontinuity at the two‐fluid interface. The resulting discrete equations are solved by an iterative procedure which displays fast convergence characteristics for small capillary numbers. Numerical results for the case of the steady movement of a fluid meniscus in a two‐dimensional channel are presented. Copyright © 2000 John Wiley & Sons, Ltd.     

A hybrid immersed boundary and material point method for simulating 3D fluid–structure interaction problems

A numerical method is developed for solving the 3D, unsteady, incompressible Navier–Stokes equations in curvilinear coordinates containing immersed boundaries (IBs) of arbitrary geometrical complexity moving and deforming under forces acting on the body. Since simulations of flow in complex geometries with deformable surfaces require special treatment, the present approach combines a hybrid immersed boundary method (HIBM) for handling complex moving boundaries and a material point method (MPM) for resolving structural stresses and movement. This combined HIBM & MPM approach is presented as an effective approach for solving fluid–structure interaction (FSI) problems. In the HIBM, a curvilinear grid is defined and the variable values at grid points adjacent to a boundary are forced or interpolated to satisfy the boundary conditions. The MPM is used for solving the equations of solid structure and communicates with the fluid through appropriate interface‐boundary conditions. The governing flow equations are discretized on a non‐staggered grid layout using second‐order accurate finite‐difference formulas. The discrete equations are integrated in time via a second‐order accurate dual time stepping, artificial compressibility scheme. Unstructured, triangular meshes are employed to discretize the complex surface of the IBs. The nodes of the surface mesh constitute a set of Lagrangian control points used for tracking the motion of the flexible body. The equations of the solid body are integrated in time via the MPM. At every instant in time, the influence of the body on the flow is accounted for by applying boundary conditions at stationary curvilinear grid nodes located in the exterior but in the immediate vicinity of the body by reconstructing the solution along the local normal to the body surface. The influence of the fluid on the body is defined through pressure and shear stresses acting on the surface of the body. The HIBM & MPM approach is validated for FSI problems by solving for a falling rigid and flexible sphere in a fluid‐filled channel. The behavior of a capsule in a shear flow was also examined. Agreement with the published results is excellent. Copyright © 2007 John Wiley & Sons, Ltd.