Fundamental validation of the finite volume particle method for 3D sloshing dynamics

The use of the finite volume particle (FVP) method is validated for three‐dimensional sloshing dynamics with a free surface by comparing with results from experiments. In the first part, two typical sloshing experiments for a single liquid phase are simulated, and slosh characteristics that include the free surface behavior and hydrodynamic pressure are reported. Moreover, the influences of the circular wall geometry and spatial resolution in the simulation are studied in a sensitivity analysis. In the second part, two sloshing problems with solid bodies are simulated to preliminarily verify the applicability of the FVP method to three‐dimensional solid bodies' motion in liquid flow. Good agreement between simulations and corresponding experiments indicates that the present FVP method well reproduces three‐dimensional sloshing behavior. Copyright © 2010 John Wiley & Sons, Ltd.     

Enhancement of the accuracy of the finite volume particle method for the simulation of incompressible flows

A finite volume particle (FVP) method for simulation of incompressible flows that provides enhanced accuracy is proposed. In this enhanced FVP method, a dummy neighbor particle is introduced for each particle in the calculation and used for the discretization of the gradient model and Laplacian model. The error‐compensating term produced by introducing the dummy neighbor particle enables higher order terms to be calculated. The proposed gradient model and Laplacian model are applied in both pressure and pressure gradient calculations. This enhanced FVP scheme provides more accurate simulations of incompressible flows. Several 2‐dimensional numerical simulations are given to confirm its enhanced performance.     

A hybrid Lagrangian–Eulerian particle‐level set method for numerical simulations of two‐fluid turbulent flows

A coupled Lagrangian interface‐tracking and Eulerian level set (LS) method is developed and implemented for numerical simulations of two‐fluid flows. In this method, the interface is identified based on the locations of notional particles and the geometrical information concerning the interface and fluid properties, such as density and viscosity, are obtained from the LS function. The LS function maintains a signed distance function without an auxiliary equation via the particle‐based Lagrangian re‐initialization technique. To assess the new hybrid method, numerical simulations of several ‘standard interface‐moving’ problems and two‐fluid laminar and turbulent flows are conducted. The numerical results are evaluated by monitoring the mass conservation, the turbulence energy spectral density function and the consistency between Eulerian and Lagrangian components. The results of our analysis indicate that the hybrid particle‐level set method can handle interfaces with complex shape change, and can accurately predict the interface values without any significant (unphysical) mass loss or gain, even in a turbulent flow. The results obtained for isotropic turbulence by the new particle‐level set method are validated by comparison with those obtained by the ‘zero Mach number’, variable‐density method. For the cases with small thermal/mass diffusivity, both methods are found to generate similar results. Analysis of the vorticity and energy equations indicates that the destabilization effect of turbulence and the stability effect of surface tension on the interface motion are strongly dependent on the density and viscosity ratios of the fluids. Copyright © 2007 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.     

Mathematical modelling of two‐phase non‐Newtonian flow in a helical pipe

Governing equations for a two‐phase 3D helical pipe flow of a non‐Newtonian fluid with large particles are derived in an orthogonal helical coordinate system. The Lagrangian approach is utilized to model solid particle trajectories. The interaction between solid particles and the fluid that carries them is accounted for by a source term in the momentum equation for the fluid. The force‐coupling method (FCM), developed by M.R. Maxey and his group, is adopted; in this method the momentum source term is no longer a Dirac delta function but is spread on a numerical mesh by using a finite‐sized envelop with a spherical Gaussian distribution. The influence of inter‐particle and particle–wall collisions is also taken into account. Copyright © 2005 John Wiley & Sons, Ltd.     

A new algorithm for surface tension model in moving particle methods

A new algorithm for the surface tension model was developed for moving particle methods. The algorithm is based on the link‐list search algorithm and the continuum surface tension (CST) model. The developed algorithm with the CST model was implemented to a kind of moving particle approach, the finite volume particle (FVP) method. The FVP method with the new algorithm was tested by oscillatory behaviour of a two‐dimensional droplet. The oscillatory period agrees well with analytical one, and the transient shape of the droplet is also in good agreement with that obtained by other numerical methods. The droplet impact on a liquid surface was also studied using the new algorithm. The deposition and splashing phenomena were clearly reproduced. Simulated spread radius of the splashing phenomena was consistent with a power law. Copyright © 2007 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.     

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.     

Three‐dimensional remeshed smoothed particle hydrodynamics for the simulation of isotropic turbulence

We present a remeshed particle‐mesh method for the simulation of three‐dimensional compressible turbulent flow. The method is related to the meshfree smoothed particle hydrodynamics method, but the present method introduces a mesh for efficient calculation of the pressure gradient, and laminar and turbulent diffusion. In addition, the mesh is used to remesh (reorganise uniformly) the particles to ensure a regular particle distribution and convergence of the method. The accuracy of the presented methodology is tested for a number of benchmark problems involving two‐ and three‐dimensional Taylor‐Green flow, thin double shear layer, and three‐dimensional isotropic turbulence. Two models were implemented, direct numerical simulations, and Smagorinsky model. Taking advantage of the Lagrangian advection, and the finite difference efficiency, the method is capable of providing quality simulations while maintaining its robustness and versatility.     

A three‐dimensional vortex particle‐in‐cell method for vortex motions in the vicinity of a wall

A new vortex particle‐in‐cell method for the simulation of three‐dimensional unsteady incompressible viscous flow is presented. The projection of the vortex strengths onto the mesh is based on volume interpolation. The convection of vorticity is treated as a Lagrangian move operation but one where the velocity of each particle is interpolated from an Eulerian mesh solution of velocity–Poisson equations. The change in vorticity due to diffusion is also computed on the Eulerian mesh and projected back to the particles. Where diffusive fluxes cause vorticity to enter a cell not already containing any particles new particles are created. The surface vorticity and the cancellation of tangential velocity at the plate are related by the Neumann conditions. The basic framework for implementation of the procedure is also introduced where the solution update comprises a sequence of two fractional steps. The method is applied to a problem where an unsteady boundary layer develops under the impact of a vortex ring and comparison is made with the experimental and numerical literature. Copyright © 2001 John Wiley & Sons, Ltd.     

A direct numerical simulation of axisymmetric liquid–gas two‐phase laminar flows in a pipe

An accurate finite‐volume‐based numerical method for the simulation of an isothermal two‐phase flow, consisting of a liquid slug translating in a non‐reacting gas in a circular pipe is presented. This method is built on a sharp interface concept and developed on an Eulerian Cartesian fixed‐grid with a cut‐cell scheme and marker points to track the moving interface. The unsteady, axisymmetric Navier–Stokes equations in both liquid and gas phases are solved separately. The mass continuity and momentum flux conditions are explicitly matched at the true surface phase boundary to determine the interface shape and movement. A quadratic curve fitting algorithm with marker points is used to yield smooth and accurate information of the interface curvatures. It is uniquely demonstrated for the first time with the current method that conservation of mass is strictly enforced for continuous infusion of flow into the domain of computation. The method has been used to compute the velocity and pressure fields and the deformation of the liquid core. It is also shown that the current method is capable of producing accurate results for a wide range of Reynolds number, Re, Weber number, We, and large property jumps at the interface. Copyright © 2010 John Wiley & Sons, Ltd.     

A high‐precision unstructured adaptive mesh technique for gas–liquid two‐phase flows

Adaptive mesh techniques are used widely in the numerical simulations of fluid flows, and the simulation results with high accuracies are obtained by appropriate mesh adaptations. However, gas–liquid two‐phase flows are still difficult to be simulated on adaptive meshes, especially on unstructured adaptive meshes, because the physical phenomena near gas–liquid interfaces are highly complicated and in general, not modeled appropriately on adaptive meshes. In this paper, a high‐precision unstructured adaptive mesh technique for gas–liquid two‐phase flows is developed and verified/validated. In the unstructured adaptive mesh technique, the PLIC algorithm is employed to simulate interfacial dynamic behaviors and, therefore, the reconstruction method for the interfaces in refined cells is developed, which satisfies the gas and liquid volume conservations and geometrical conservations of interfaces. In addition, the physics‐based consideration is performed on the momentum calculations near interfaces, and the calculation method with gas and liquid momentum conservations is developed. For verification, the slotted‐disk revolution problem is solved. As a result, the unstructured adaptive mesh technique succeeds in reproducing the slotted‐disk shape accurately and well maintaining the shape after one full‐revolution. The dam‐break problem is also simulated and the momentum conservative calculation method succeeds in providing physically appropriate results, which show good agreements with experimental data. Therefore, it is confirmed that the developed unstructured adaptive mesh technique is very efficient to simulate gas–liquid two‐phase flows accurately. Copyright © 2010 John Wiley & Sons, Ltd.     

Simulating gas–liquid multiphase flow using meshless Lagrangian method

A multiphase flow model has been established based on a moving particle semi‐implicit method. A surface tension model is introduced to the particle method to improve the numerical accuracy and stability. Several computational techniques are employed to simplify the numerical procedure and further improve the accuracy. A particle fraction multiphase flow model is developed and verified by a two‐phase Poiseuille flow. The multiphase surface tension model is discussed in detail, and an ethanol drop case is introduced to verify the surface tension model. A simple dam break is simulated to demonstrate the improvements with various modifications in particle method along with a new boundary condition. Finally, we simulate several bubble rising cases to show the capacity of this new model in simulating gas–liquid multiphase flow with large density ratio difference between phases. The comparisons among numerical results of mesh‐based model, experimental data, and the present model, indicate that the new multiphase particle method is acceptable in gas–liquid multiphase fluids simulation. Copyright © 2014 John Wiley & Sons, Ltd.     

An immersed‐boundary finite‐volume method for direct simulation of flows with suspended paramagnetic particles

In the present study, we have proposed an immersed‐boundary finite‐volume method for the direct numerical simulation of flows with inertialess paramagnetic particles suspended in a nonmagnetic fluid under an external magnetic field without the need for any model such as the dipole–dipole interaction. In the proposed method, the magnetic field (or force) is described by the numerical solution of the Maxwell equation without current, where the smoothed representation technique is employed to tackle the discontinuity of magnetic permeability across the particle–fluid interface. The flow field, on the other hand, is described by the solution of the continuity and momentum equations, where the discrete‐forcing‐based immersed‐boundary method is employed to satisfy the no‐slip condition at the interface. To validate the method, we performed numerical simulations on the two‐dimensional motion of two and three paramagnetic particles in a nonmagnetic fluid subjected to an external uniform magnetic field and then compared the results with the existing finite‐element and semi‐analytical solutions. Comparison shows that the proposed method is robust in the direct simulation of such magnetic particulate flows. This method can be extended to more general flows without difficulty: three‐dimensional particulate flows, flows with a great number of particles, or flows under an arbitrary external magnetic field. Copyright © 2010 John Wiley & Sons, Ltd.     

Fitted finite element discretization of two‐phase Stokes flow

We propose a novel fitted finite element method for two‐phase Stokes flow problems that uses piecewise linear finite elements to approximate the moving interface. The method can be shown to be unconditionally stable. Moreover, spherical stationary solutions are captured exactly by the numerical approximation. In addition, the meshes describing the discrete interface in general do not deteriorate in time, which means that in numerical simulations, a smoothing or a remeshing of the interface mesh is not necessary. We present several numerical experiments for our numerical method, which demonstrate the accuracy and robustness of the proposed algorithm. Copyright © 2016 John Wiley & Sons, Ltd.     

Simplex space‐time meshes in two‐phase flow simulations

In this paper, we present the numerical solution of 2‐phase flow problems of engineering significance with a space‐time finite element method that allows for local temporal refinement. Arbitrary temporal refinement is applied to preselected regions of the mesh and is governed by a quantity that is part of the solution process, namely, the interface position in 2‐phase flow. Because of local effects such as surface tension, jumps in material properties, etc, the interface can in general be considered a region that requires high flexibility and high resolution, both in space and in time. The new method, which leads to tetrahedral (for 2D problems) and pentatope (for 3D problems) meshes, offers an efficient yet accurate approach to the underlying 2‐phase flow problems.     

Interface‐fitted moving mesh method for axisymmetric two‐phase flow in microchannels

A boundary‐fitted moving mesh scheme is presented for the simulation of two‐phase flow in two‐dimensional and axisymmetric geometries. The incompressible Navier‐Stokes equations are solved using the finite element method, and the mini element is used to satisfy the inf‐sup condition. The interface between the phases is represented explicitly by an interface adapted mesh, thus allowing a sharp transition of the fluid properties. Surface tension is modelled as a volume force and is discretized in a consistent manner, thus allowing to obtain exact equilibrium (up to rounding errors) with the pressure gradient. This is demonstrated for a spherical droplet moving in a constant flow field. The curvature of the interface, required for the surface tension term, is efficiently computed with simple but very accurate geometric formulas. An adaptive moving mesh technique, where smoothing mesh velocities and remeshing are used to preserve the mesh quality, is developed and presented. Mesh refinement strategies, allowing tailoring of the refinement of the computational mesh, are also discussed. Accuracy and robustness of the present method are demonstrated on several validation test cases. The method is developed with the prospect of being applied to microfluidic flows and the simulation of microchannel evaporators used for electronics cooling. Therefore, the simulation results for the flow of a bubble in a microchannel are presented and compared to experimental data.     

Theoretical analysis for achieving high‐order spatial accuracy in Lagrangian/Eulerian source terms

In a fully coupled Lagrangian/Eulerian two‐phase calculation, the source terms from computational particles must be agglomerated to nearby gas‐phase nodes. Existing methods are capable of accomplishing this particle‐to‐gas coupling with second‐order accuracy. However, higher‐order methods would be useful for applications such as two‐phase direct numerical simulation and large eddy simulation. A theoretical basis is provided for producing high spatial accuracy in particle‐to‐gas source terms with low computational cost. The present work derives fourth‐ and sixth‐order accurate methods, and the procedure for even higher accuracy is discussed. The theory is also expanded to include two‐ and three‐dimensional calculations. One‐ and two‐dimensional tests are used to demonstrate the convergence of this method and to highlight problems with statistical noise. Finally, the potential for application in computational fluid dynamics codes is discussed. It is concluded that high‐order kernels have practical benefits only under limited ranges of statistical and spatial resolution. Additionally, convergence demonstrations with full CFD codes will be extremely difficult due to the worsening of statistical errors with increasing mesh resolution. Copyright © 2006 John Wiley & Sons, Ltd.     

Finite‐volume multi‐stage schemes for shallow‐water flow simulations

A finite‐volume multi‐stage (FMUSTA) scheme is proposed for simulating the free‐surface shallow‐water flows with the hydraulic shocks. On the basis of the multi‐stage (MUSTA) method, the original Riemann problem is transformed to an independent MUSTA mesh. The local Lax–Friedrichs scheme is then adopted for solving the solution of the Riemann problem at the cell interface on the MUSTA mesh. The resulting first‐order monotonic FMUSTA scheme, which does not require the use of the eigenstructure and the special treatment of entropy fixes, has the generality as well as simplicity. In order to achieve the high‐resolution property, the monotonic upstream schemes for conservation laws (MUSCL) method are used. For modeling shallow‐water flows with source terms, the surface gradient method (SGM) is adopted. The proposed schemes are verified using the simulations of six shallow‐water problems, including the 1D idealized dam breaking, the steady transcritical flow over a hump, the 2D oblique hydraulic jump, the circular dam breaking and two dam‐break experiments. The simulated results by the proposed schemes are in satisfactory agreement with the exact solutions and experimental data. It is demonstrated that the proposed FMUSTA schemes have superior overall numerical accuracy among the schemes tested such as the commonly adopted Roe and HLL schemes. Copyright © 2007 John Wiley & Sons, Ltd.