Parallel finite element simulation of mooring forces on floating objects

The coupling between the equations governing the free‐surface flows, the six degrees of freedom non‐linear rigid body dynamics, the linear elasticity equations for mesh‐moving and the cables has resulted in a fluid‐structure interaction technology capable of simulating mooring forces on floating objects. The finite element solution strategy is based on a combination approach derived from fixed‐mesh and moving‐mesh techniques. Here, the free‐surface flow simulations are based on the Navier–Stokes equations written for two incompressible fluids where the impact of one fluid on the other one is extremely small. An interface function with two distinct values is used to locate the position of the free‐surface. The stabilized finite element formulations are written and integrated in an arbitrary Lagrangian–Eulerian domain. This allows us to handle the motion of the time dependent geometries. Forces and momentums exerted on the floating object by both water and hawsers are calculated and used to update the position of the floating object in time. In the mesh moving scheme, we assume that the computational domain is made of elastic materials. The linear elasticity equations are solved to obtain the displacements for each computational node. The non‐linear rigid body dynamics equations are coupled with the governing equations of fluid flow and are solved simultaneously to update the position of the floating object. The numerical examples includes a 3D simulation of water waves impacting on a moored floating box and a model boat and simulation of floating object under water constrained with a cable. Copyright © 2003 John Wiley & Sons, Ltd.     

On the usage of NURBS as interface representation in free‐surface flows

When simulating free‐surface flows using the finite element method, there are many cases where the governing equations require information which must be derived from the available discretized geometry. Examples are curvature or normal vectors. The accurate computation of this information directly from the finite element mesh often requires a high degree of refinement—which is not necessarily required to obtain an accurate flow solution. As a remedy and an option to be able to use coarser meshes, the representation of the free surface using non‐uniform rational B‐splines (NURBS) curves or surfaces is investigated in this work. The advantages of a NURBS parameterization in comparison with the standard approach are discussed. In addition, it is explored how the pressure jump resulting from surface tension effects can be handled using doubled interface nodes. Numerical examples include the computation of surface tension in a two‐phase flow as well as the computation of normal vectors as a basis for mesh deformation methods. For these examples, the improvement of the numerical solution compared with the standard approaches on identical meshes is shown. Copyright © 2011 John Wiley & Sons, Ltd.     

Computation of turbulent free‐surface flows around modern ships

This paper presents the calculated results for three classes of typical modern ships in modelling of ship‐generated waves. Simulations of turbulent free‐surface flows around ships are performed in a numerical water tank, based on the FINFLO‐RANS SHIP solver developed at Helsinki University of Technology. The Reynolds‐averaged Navier–Stokes (RANS) equations with the artificial compressibility and the non‐linear free‐surface boundary conditions are discretized by means of a cell‐centred finite‐volume scheme. The convergence performance is improved with the multigrid method. A free surface is tracked using a moving mesh technology, in which the non‐linear free‐surface boundary conditions are given on the actual location of the free surface. Test cases recommended are a container ship, a US Navy combatant and a tanker. The calculated results are compared with the experimental data available in the literature in terms of the wave profiles, wave pattern, and turbulent flow fields for two turbulence models, Chien's low Reynolds number k–εmodel and Baldwin–Lomax's model. Furthermore, the convergence performance, the grid refinement study and the effect of turbulence models on the waves have been investigated. Additionally, comparison of two types of the dynamic free‐surface boundary conditions is made. Copyright © 2003 John Wiley& Sons, Ltd.     

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 adaptive numerical scheme for solving incompressible 2‐phase and free‐surface flows

In this paper, we present a numerical scheme for solving 2‐phase or free‐surface flows. Here, the interface/free surface is modeled using the level‐set formulation, and the underlying mesh is adapted at each iteration of the flow solver. This adaptation allows us to obtain a precise approximation for the interface/free‐surface location. In addition, it enables us to solve the time‐discretized fluid equation only in the fluid domain in the case of free‐surface problems. Fluids here are considered incompressible. Therefore, their motion is described by the incompressible Navier‐Stokes equation, which is temporally discretized using the method of characteristics and is solved at each time iteration by a first‐order Lagrange‐Galerkin method. The level‐set function representing the interface/free surface satisfies an advection equation that is also solved using the method of characteristics. The algorithm is completed by some intermediate steps like the construction of a convenient initial level‐set function (redistancing) as well as the construction of a convenient flow for the level‐set advection equation. Numerical results are presented for both bifluid and free‐surface problems.     

Local mesh adaptation technique for front tracking problems

A numerical model is developed for the simulation of moving interfaces in viscous incompressible flows. The model is based on the finite element method with a pseudo-concentration technique to track the front. Since a Eulerian approach is chosen, the interface is advected by the flow through a fixed mesh. Therefore, material discontinuity across the interface cannot be described accurately. To remedy this problem, the model has been supplemented with a local mesh adaptation technique. This latter consists in updating the mesh at each time step to the interface position, such that element boundaries lie along the front. It has been implemented for unstructured triangular finite element meshes. The outcome of this technique is that it allows an accurate treatment of material discontinuity across the interface and, if necessary, a modelling of interface phenomena such as surface tension by using specific boundary elements. For illustration, two examples are computed and presented in this paper: the broken dam problem and the Rayleigh–Taylor instability. Good agreement has been obtained in the comparison of the numerical results with theory or available experimental data. © 1998 John Wiley & Sons, Ltd.     

Cahn–Hilliard modeling of particles suspended in two‐phase flows

In this paper, we present a model for the dynamics of particles suspended in two‐phase flows by coupling the Cahn–Hilliard theory with the extended finite element method (XFEM). In the Cahn–Hilliard model the interface is considered to have a small but finite thickness, which circumvents explicit tracking of the interface. For the direct numerical simulation of particle‐suspended flows, we incorporate an XFEM, in which the particle domain is decoupled from the fluid domain. To cope with the movement of the particles, a temporary ALE scheme is used for the mapping of field variables at the previous time levels onto the computational mesh at the current time level. By combining the Cahn–Hilliard model with the XFEM, the particle motion at an interface can be simulated on a fixed Eulerian mesh without any need of re‐meshing. The model is general, but to demonstrate and validate the technique, here the dynamics of a single particle at a fluid–fluid interface is studied. First, we apply a small disturbance on a particle resting at an interface between two fluids, and investigate the particle movement towards its equilibrium position. In particular, we are interested in the effect of interfacial thickness, surface tension, particle size and viscosity ratio of two fluids on the particle movement towards its equilibrium position. Finally, we show the movement of a particle passing through multiple layers of fluids. Copyright © 2011 John Wiley & Sons, Ltd.     

Coupling methods between finite element–based Boussinesq‐type wave and particle‐based free‐surface flow models

We develop one‐way coupling methods between a Boussinesq‐type wave model based on the discontinuous Galerkin finite element method and a free‐surface flow model based on a mesh‐free particle method to strike a balance between accuracy and computational cost. In our proposed model, computation of the wave model in the global domain is conducted first, and the nonconstant velocity profiles in the vertical direction are reproduced by using its results. Computation of the free‐surface flow is performed in a local domain included within the global domain with interface boundaries that move along the reproduced velocity field in a Lagrangian fashion. To represent the moving interfaces, we used a polygon wall boundary model for mesh‐free particle methods. Verification and validation tests of our proposed model are performed, and results obtained by the model are compared with theoretical values and experimental results to show its accuracy and applicability.     

Numerical simulation of transient planar flow of a viscoelastic material with two moving free surfaces

In this study, we examine the numerical simulation of transient viscoelastic flows with two moving free surfaces. A modified Galerkin finite element method is implemented to the two-dimensional non-steady motion of the fluid of the Oldroyd-B type. The fluid is initially placed between two parallel plates and bounded by two straight free boundaries. In this Lagrangian finite element method, the spatial mesh deforms in time along with the moving free boundaries. The unknown shape of the free surfaces is determined with the flow field u, v, τ, p by the deformable finite element method, combined with a predictor-corrector scheme in an uncoupled fashion. The moving free surfaces and fluid motion of both Newtonian and non-Newtonian flows are investigated. The results include the influence of surface tension, fluid inertia and elasticity.     

NUMERICAL SIMULATION OF INCOMPRESSIBLE FLOWS WITH MOVING INTERFACES

A numerical model has been developed for the 2D simulation of free surface flows or, more generally speaking, moving interface ones. The bulk fluids on both sides of the interface are taken into account in simulating the incompressible laminar flow state. In the case of heat transfer the whole system, i.e. walls as well as possible obstacles, is considered. This model is based on finite element analysis with an Eulerian approach and an unstructured fixed mesh. A special technique to localize the interface allows its temporal evolution through this mesh. Several numerical examples are presented to demonstrate the capabilities of the model. © 1997 by John Wiley & Sons, Ltd.     

Anisotropic adaptive stabilized finite element solver for RANS models

Aerodynamic characteristics of various geometries are predicted using a finite element formulation coupled with several numerical techniques to ensure stability and accuracy of the method. First, an edge‐based error estimator and anisotropic mesh adaptation are used to detect automatically all flow features under the constraint of a fixed number of elements, thus controlling the computational cost. A variational multiscale‐stabilized finite element method is used to solve the incompressible Navier‐Stokes equations. Finally, the Spalart‐Allmaras turbulence model is solved using the streamline upwind Petrov‐Galerkin method. This paper is meant to show that the combination of anisotropic unsteady mesh adaptation with stabilized finite element methods provides an adequate framework for solving turbulent flows at high Reynolds numbers. The proposed method was validated on several test cases by confrontation with literature of both numerical and experimental results, in terms of accuracy on the prediction of the drag and lift coefficients as well as their evolution in time for unsteady cases.     

CFD simulation of two‐ and three‐dimensional free‐surface flow

The paper describes the implementation of moving‐mesh and free‐surface capabilities within a 3‐d finite‐volume Reynolds‐averaged‐Navier–Stokes solver, using surface‐conforming multi‐block structured meshes. The free‐surface kinematic condition can be applied in two ways: enforcing zero net mass flux or solving the kinematic equation by a finite‐difference method. The free surface is best defined by intermediate control points rather than the mesh vertices. Application of the dynamic boundary condition to the piezometric pressure at these points provides a hydrostatic restoring force which helps to eliminate any unnatural free‐surface undulations. The implementation of time‐marching methods on moving grids are described in some detail and it is shown that a second‐order scheme must be applied in both scalar‐transport and free‐surface equations if flows driven by free‐surface height variations are to be computed without significant wave attenuation using a modest number of time steps. Computations of five flows of theoretical and practical interest—forced motion in a pump, linear waves in a tank, quasi‐1d flow over a ramp, solitary wave interaction with a submerged obstacle and 3‐d flow about a surface‐penetrating cylinder—are described to illustrate the capabilities of our code and methods. Copyright © 2003 John Wiley & Sons, Ltd.     

Prediction of extrudate swell in polymer melt extrusion using an Arbitrary Lagrangian Eulerian (ALE) based finite element method

Accurate prediction of extrudate (die) swell in polymer melt extrusion is important as this helps in appropriate die design for profile extrusion applications. Extrudate swell prediction has shown significant difficulties due to two key reasons. The first is the appropriate representation of the constitutive behavior of the polymer melt. The second is regarding the simulation of the free surface, which requires special techniques in the traditionally used Eulerian framework. In this paper we propose a method for simulation of extrudate swell using an Arbitrary Lagrangian Eulerian (ALE) technique based finite element formulation. The ALE technique provides advantages of both Lagrangian and Eulerian frameworks by allowing the computational mesh to move in an arbitrary manner, independent of the material motion. In the present method, a fractional-step ALE technique is employed in which the Lagrangian phase of material motion and convection arising out of mesh motion are decoupled. In the first step, the relevant flow and constitutive equations are solved in Lagrangian framework. The simpler representation of polymer constitutive equations in a Lagrangian framework avoids the difficulties associated with convective terms thereby resulting in a robust numerical formulation besides allowing for natural evolution of the free surface with the flow. In the second step, mesh is moved in ALE mode and the associated convection of the variables due to relative motion of the mesh is performed using a Godunov type scheme. While the mesh is fixed in space in the die region, the nodal points of the mesh on the extrudate free surface are allowed to move normal to flow direction with special rules to facilitate the simulation of swell. A differential exponential Phan Thien Tanner (PTT) model is used to represent the constitutive behavior of the melt. Using this method we simulate extrudate swell in planar and axisymmetric extrusion with abrupt contraction ahead of the die exit. This geometry allows the extrudate to have significant memory for shorter die lengths and acts as a good test for swell predictions. We demonstrate that our predictions of extrudate swell match well with reported experimental and numerical simulations.     

Numerical simulation of 3D free surface flows,with multiple incompressible immiscible phases. Applications to impulse waves

A numerical method for the solution to the density‐dependent incompressible Navier–Stokes equations modeling the flow of N immiscible incompressible liquid phases with a free surface is proposed. It allows to model the flow of an arbitrary number of liquid phases together with an additional vacuum phase separated with a free surface. It is based on a volume‐of‐fluid approach involving N indicator functions (one per phase, identified by its density) that guarantees mass conservation within each phase. An additional indicator function for the whole liquid domain allows to treat boundary conditions at the interface between the liquid domain and a vacuum. The system of partial differential equations is solved by implicit operator splitting at each time step: first, transport equations are solved by a forward characteristics method on a fine Cartesian grid to predict the new location of each liquid phase; second, a generalized Stokes problem with a density‐dependent viscosity is solved with a FEM on a coarser mesh of the liquid domain. A novel algorithm ensuring the maximum principle and limiting the numerical diffusion for the transport of the N phases is validated on benchmark flows. Then, we focus on a novel application and compare the numerical and physical simulations of impulse waves, that is, waves generated at the free surface of a water basin initially at rest after the impact of a denser phase. A particularly useful application in hydraulic engineering is to predict the effects of a landslide‐generated impulse wave in a reservoir. Copyright © 2014 John Wiley & Sons, Ltd.     

A free surface finite element model for low Froude number mould filling problems on fixed meshes

The simulation of low Froude number mould filling problems on fixed meshes presents significant difficulties. As the Froude number decreases, the coupling between the position of the interface and the resulting flow field increases. The usual two‐phase flow model provides poor results for such flow. In order to overcome the difficulties, a free surface model that applies boundary conditions at the interface accurately is used. Moreover, the use of wall laws on curved boundaries also fails in the case of low Froude number flows. To solve this second problem, we combine wall laws with ‘do nothing’ boundary conditions. A special feature of our approach is that ‘do nothing’ boundary conditions are only applied in the normal direction. These two key ingredients together with the Level Set method allow us to simulate three‐dimensional mould filling problems borrowed directly from the foundry. Copyright © 2010 John Wiley & Sons, Ltd.     

Numerical simulation of unsteady multidimensional free surface motions by level set method

This paper presents a numerical method that couples the incompressible Navier–Stokes equations with the level set method in a curvilinear co‐ordinate system for study of free surface flows. The finite volume method is used to discretize the governing equations on a non‐staggered grid with a four‐step fractional step method. The free surface flow problem is converted into a two‐phase flow system on a fixed grid in which the free surface is implicitly captured by the zero level set. We compare different numerical schemes for advection of the level set function in a generalized curvilinear format, including the third order quadratic upwind interpolation for convective kinematics (QUICK) scheme, and the second and third order essentially non‐oscillatory (ENO) schemes. The level set equations of evolution and reinitialization are validated with benchmark cases, e.g. a stationary circle, a rotating slotted disk and stretching of a circular fluid element. The coupled system is then applied to a travelling solitary wave, and two‐ and three‐dimensional dam breaking problems. Some interesting free surface phenomena are revealed by the computational results, such as, the large free surface vortices, air entrapment and splashing of the water surge front. The computational results are in excellent agreement with theoretical predictions and experimental data, where they are available. Copyright © 2003 John Wiley & Sons, Ltd.     

A two‐phase flow interface capturing finite element method

A new interface capturing algorithm is proposed for the finite element simulation of two‐phase flows. It relies on the solution of an advection equation for the interface between the two phases by a streamline upwind Petrov–Galerkin (SUPG) scheme combined with an adaptive mesh refinement procedure and a filtering technique. This method is illustrated in the case of a Rayleigh–Taylor two‐phase flow problem governed by the Stokes equations. Copyright © 2006 John Wiley & Sons, Ltd.     

A weakly compressible MPS method for modeling of open‐boundary free‐surface flow

A mesh‐free particle method, based on the moving particle semi‐implicit (MPS) interaction model, has been developed for the simulation of two‐dimensional open‐boundary free‐surface flows. The incompressibility model in the original MPS has been replaced with a weakly incompressible model. The effect of this replacement on the efficiency and accuracy of the model has been investigated. The new inflow–outflow boundary conditions along with the particle recycling strategy proposed in this study extend the application of the model to open‐boundary problems. The final model is able to simulate open‐boundary free surface flow in cases of large deformation and fragmentation of free surface. The models and proposed algorithms have been validated and applied to sample problems. The results confirm the model's efficiency and accuracy. Copyright © 2009 John Wiley & Sons, Ltd.     

An algorithm for rotating axisymmetric flows: model,validation and industrial applications

The study of axisymmetric flows is of interest not only from an academic point of view, due to the existence of exact solutions of Navier–Stokes equations, but also from an industrial point of view, since these kind of flows are frequently found in several applications. In the present work the development and implementation of a finite element algorithm to solve Navier–Stokes equations with axisymmetric geometry and boundary conditions is presented. Such algorithm allows the simulation of flows with tangential velocity, including free surface flows, for both laminar and turbulent conditions. Pseudo‐concentration technique is used to model the free surface (or the interface between two fluids) and the k–ε model is employed to take into account turbulent effects. The finite element model is validated by comparisons with analytical solutions of Navier–Stokes equations and experimental measurements. Two different industrial applications are presented. Copyright © 2005 John Wiley & Sons, Ltd.