High Order Semi-implicit Multistep Methods for Time-Dependent Partial Differential Equations

We consider the construction of semi-implicit linear multistep methods that can be applied to time-dependent PDEs where the separation of scales in additive form,typically used in implicit-explicit(IMEX)methods,is not possible.As shown in Boscarino et al.(J.Sci.Comput.68:975-1001,2016)for Runge-Kutta methods,these semi-implicit techniques give a great flexibility,and allow,in many cases,the construction of simple linearly implicit schemes with no need of iterative solvers.In this work,we develop a general setting for the construction of high order semi-implicit linear multistep methods and analyze their stability properties for a prototype lineal'advection-diffusion equation and in the setting of strong stability preserving(SSP)methods.Our findings are demonstrated on several examples,including nonlinear reaction-diffusion and convection-diffusion problems.

Moving-Particle Semi-Implicit Method for Vortex Patterns and Roll Damping of 2D Ship Sections

A meshless method, Moving-Particle Semi-hnplicit Method （MPS） is presented in this paper to simulate the rolling of different 2D ship sections. Sections S. S. 0.5, S.S. 5.0 and S. S. 7.0 of series 60 with CB = 0.6 are chosen for the simulation. It shows that the result of MPS is very close to results of experiments or mesh-numerical simulations. In the simulation of MPS, vortices are found periodically in bilges of ship sections. In section S. S. 5.0 and section S. S. 7.0, which are close to the middle ship, two little vortices are found at different bilges of the section, in section S. S. 0.5, which is close to the bow, only one big vortex is found at the bottom of the section, these vortices patterns are consistent with the theory of Ikeda. The distribution of shear stress and pressure on the rolling hull of ship section is calculated. When vortices are in bilges of the section, the sign clmnge of pressure can be found, but in section S. S. 0.5, there is no sign change of pressure because only one vortex in the bottom of the section. With shear stress distribution, it can be found the shear stress in bilges is bigger than that at other part of the ship section. As the free surface is considered, the shear stress of both sides near the free surface is close to zero and even sign changed.

Numerical simulation of viscous liquid sloshing by moving-particle semi-implicit method

A meshless numerical simulation method, the moving-particle semi-implicit method （MPS） is presented in this paper to study the sloshing phenomenon in ocean and naval engineering. As a meshless method, MPS uses particles to replace the mesh in traditional methods, the governing equations are discretized by virtue of the relationship of particles, and the Poisson equation of pressure is solved by incomplete Cholesky conjugate gradient method （ICCG）, the free surface is tracked by the change of numerical density. A numerical experiment of viscous liquid sloshing tank was presented and compared with the result got by the difference method with the VOF, and an additional modification step was added to make the simulation more stable. The results show that the MPS method is suitable for the simulation of viscous liquid sloshing, with the advantage in arranging the particles easily, especially on some complex curved surface.

Large deformation simulations of geomaterials using moving particle semi-implicit method

Numerical simulation tools are required to describe large deformations of geomaterials for evaluating the risk of geo-disasters. This study focused on moving particle semi-implicit(MPS) method, which is a Lagrangian gridless particle method, and investigated its performance and stability to simulate large deformation of geomaterials. A calculation method was developed using geomaterials modeled as Bingham fluids to improve the original MPS method and enhance its stability. Two numerical tests showed that results from the improved MPS method was in good agreement with the theoretical value.Furthermore, numerical simulations were calibrated by laboratory experiments. It showed that the simulation results matched well with the experimentally observed free-surface configurations for flowing sand. In addition, the model could generally predict the time-history of the impact force. The MPS method could be a useful tool to evaluate large deformation of geomaterials.

The Semi-implicit Euler Method for Stochastic Pantograph Equations with Jumps

In this paper,we present the semi-implicit Euler（SIE）numerical solution for stochastic pantograph equations with jumps and prove that the SIE approximation solution converges to the exact solution in the mean-square sense under the Local Lipschitz condition.

A Novel Staggered Semi-implicit Space-Time Discontinuous Galerkin Method for the Incompressible Navier-Stokes Equations

A new high-order accurate staggered semi-implicit space-time discontinuous Galerkin(DG)method is presented for the simulation of viscous incompressible flows on unstructured triangular grids in two space dimensions.The staggered DG scheme defines the discrete pressure on the primal triangular mesh,while the discrete velocity is defined on a staggered edge-based dual quadrilateral mesh.In this paper,a new pair of equal-order-interpolation velocity-pressure finite elements is proposed.On the primary triangular mesh(the pressure elements),the basis functions are piecewise polynomials of degree N and are allowed to jump on the boundaries of each triangle.On the dual mesh instead(the velocity elements),the basis functions consist in the union of piecewise polynomials of degree N on the two subtriangles that compose each quadrilateral and are allowed to jump only on the dual element boundaries,while they are continuous inside.In other words,the basis functions on the dual mesh arc built by continuous finite elements on the subtriangles.This choice allows the construction of an efficient,quadrature-free and memory saving algorithm.In our coupled space-time pressure correction formulation for the incompressible Navier-Stokes equations,the arbitrary high order of accuracy in time is achieved through tire use of time-dependent test and basis functions,in combination with simple and efficient Picard iterations.Several numerical tests on classical benchmarks confirm that the proposed method outperforms existing staggered semi-implicit space-time DG schemes,not only from a computer memory point of view,but also concerning the computational time.

Numerical Simulation of the Solitary Wave Interacting with an Elastic Structure Using MPS-FEM Coupled Method

Fluid-Structure Interaction(FSI) caused by fluid impacting onto a flexible structure commonly occurs in naval architecture and ocean engineering. Research on the problem of wave-structure interaction is important to ensure the safety of offshore structures. This paper presents the Moving Particle Semi-implicit and Finite Element Coupled Method(MPS-FEM) to simulate FSI problems. The Moving Particle Semi-implicit(MPS) method is used to calculate the fluid domain, while the Finite Element Method(FEM) is used to address the structure domain. The scheme for the coupling of MPS and FEM is introduced first. Then, numerical validation and convergent study are performed to verify the accuracy of the solver for solitary wave generation and FSI problems. The interaction between the solitary wave and an elastic structure is investigated by using the MPS-FEM coupled method.

Numerical Simulation of Sloshing with Large Deforming Free Surface by MPS-LES Method

Moving particle semi-implicit （MPS） method is a fully Lagrangian particle method which can easily solve problems with violent free surface. Although it has demonstrated its advantage in ocean engineering applications, it still has some defects to be improved. In this paper, MPS method is extended to the large eddy simulation （LES） by coupling with a sub-particle-scale （SPS） turbulence model. The SPS turbulence model turns into the Reynolds stress terms in the filtered momentum equation, and the Smagorinsky model is introduced to describe the Reynolds stress terms. Although MPS method has the advantage in the simulation of the free surface flow, a lot of non-free surface particles are treated as free surface particles in the original MPS model. In this paper, we use a new free surface tracing method and the key point is ＂neighbor particle＂. In this new method, the zone around each particle is divided into eight parts, and the particle will be treated as a free surface particle as long as there are no ＂neighbor particles＂ in any two parts of the zone. As the number density parameter judging method has a high efficiency for the free surface particles tracing, we combine it with the neighbor detected method. First, we select out the particles which may be mistreated with high probabilities by using the number density parameter judging method. And then we deal with these particles with the neighbor detected method. By doing this, the new mixed free surface tracing method can reduce the mistreatment problem efficiently. The serious pressure fluctuation is an obvious defect in MPS method, and therefore an area-time average technique is used in this paper to remove the pressure fluctuation with a quite good result. With these improvements, the modified MPS-LES method is applied to simulate liquid sloshing problems with large deforming free surface. Results show that the modified MPS-LES method can simulate the large deforming free surface easily. It can not only capture the large impact pressure accurately on rolling tank wall but also can generate all physical phenomena successfully. The good agreement between numerical and experimental results proves that the modified MPS-LES method is a good CFD methodology in free surface flow simulations.

Numerical Simulation of Sloshing Using the MPS-FSI Method with Large Eddy Simulation

A numerical model has been developed to study sloshing of turbulent flow in a tank with elastic baffles. The Moving-Particle Semi-implicit method(MPS) is a kind of meshless Lagrangian calculation method. The large eddy simulation(LES) approach is employed to model the turbulence by using the Smagorinsky Sub-Particle Scale(SPS)closure model. This paper uses MPS-FSI method with LES to simulate the interaction between free surface flow and a thin elastic baffle in sloshing. Then, the numerical model is validated, and the numerical solution has good agreement with experimental data for sloshing in a tank with elastic baffles. Furthermore, under external excitations,the MPS is applied to viscous laminar flow and turbulent flow, with both the deformation of elastic baffles and the wave height of the free surface are compared with each other. Besides, the impact pressure with/without baffles and wave height of free surface are investigated and discussed in detail. Finally, preliminary simulations are carried out in the damage problem of elastic baffles, taking the advantage of the MPS-FSI method in computations of the fluid–structure interaction with large deformation.

Computational-Analysis of the Non-Isothermal Dynamics of the Gravity-Driven Flow of Viscoelastic-Fluid-Based Nanofluids Down an Inclined Plane

The paper explores the gravity-driven flow of the thin film of a viscoelastic-fluid-based nanofluids(VFBN)along an inclined plane under non-isothermal conditions and subjected to convective cooling at the free-surface.The Newton’s law of cooling is used to model the convective heat-exchange with the ambient at the free-surface.The Giesekus viscoelastic constitutive model,with appropriate modifications to account for non-isothermal effects,is employed to describe the polymeric effects.The unsteady and coupled non-linear partial differential equations(PDEs)describing the model problem are obtained and solved via efficient semi-implicit numerical schemes based on finite difference methods(FDM)implemented in Matlab.The response of the VFBN velocity,temperature,thermal-conductivity and polymeric-stresses to variations in the volume-fraction of embedded nanoparticles is investigated.It is shown that these quantities all increase as the nanoparticle volume-fraction becomes higher.

UNCONDITIONAL CONVERGENCE AND ERROR ESTIMATES OF A FULLY DISCRETE FINITE ELEMENT METHOD FOR THE MICROPOLAR NAVIER-STOKES EQUATIONS

In this paper,we consider the initial-boundary value problem(IBVP)for the micropolar Naviers-Stokes equations(MNSE)and analyze a first order fully discrete mixed finite element scheme.We first establish some regularity results for the solution of MNSE,which seem to be not available in the literature.Next,we study a semi-implicit time-discrete scheme for the MNSE and prove L2-H1 error estimates for the time discrete solution.Furthermore,certain regularity results for the time discrete solution are establishes rigorously.Based on these regularity results,we prove the unconditional L2-H1 error estimates for the finite element solution of MNSE.Finally,some numerical examples are carried out to demonstrate both accuracy and efficiency of the fully discrete finite element scheme.

Numerical simulation of liquid sloshing in a spherical tank by MPS method

This paper investigates the sloshing phenomena in a spherical liquid tank using the moving particle semi-implicit(MPS)method,a crucial study in fluid dynamics.Distinct from previous research focused on rectangular or LNG tanks,this work explores the unique motion patterns inherent to spherical geometries.The accuracy of our in-house MPS solver MLParticle-SJTU is validated against experimental data and finite volume method(FVM).And the MPS method reveals a closer alignment with experimental outcomes,which suggests that MPS method is particularly effective for modeling complex,non-linear fluid behaviors.Then the fluid’s response to excitation at its natural frequency is simulated,showcasing vigorous sloshing and rotational motion.Detailed analyses of the fluid motion are conducted by drawing streamline diagrams,velocity vector diagrams,and vorticity maps.The fluid’s motion response is explored using both time-domain and frequency-domain curves of the fluid centroid,as well as the sloshing force.

Towards development of enhanced fully-Lagrangian mesh-free computational methods for fluid-structure interaction

Simulation of incompressible fluid flow-elastic structure interactions is targeted by using fully-Lagrangian mesh-free computational methods. A projection-based fluid model（moving particle semi-implicit（MPS）） is coupled with either a Newtonian or a Hamiltonian Lagrangian structure model（MPS or HMPS） in a mathematically-physically consistent manner. The fluid model is founded on the solution of Navier-Stokes and continuity equations. The structure models are configured either in the framework of Newtonian mechanics on the basis of conservation of linear and angular momenta, or Hamiltonian mechanics on the basis of variational principle for incompressible elastodynamics. A set of enhanced schemes are incorporated for projection-based fluid model（Enhanced MPS）, thus, the developed coupled solvers for fluid structure interaction（FSI） are referred to as Enhanced MPS-MPS and Enhanced MPS-HMPS. Besides, two smoothed particle hydrodynamics（SPH）-based FSI solvers, being developed by the authors, are considered and their potential applicability and comparable performance are briefly discussed in comparison with MPS-based FSI solvers. The SPH-based FSI solvers are established through coupling of projection-based incompressible SPH（ISPH） fluid model and SPH-based Newtonian/Hamiltonian structure models, leading to Enhanced ISPH-SPH and Enhanced ISPH-HSPH. A comparative study is carried out on the performances of the FSI solvers through a set of benchmark tests, including hydrostatic water column on an elastic plate,high speed impact of an elastic aluminum beam, hydroelastic slamming of a marine panel and dam break with elastic gate.

Dynamic analysis of dam-reservoir-foundation interaction using finite difference technique

Time domain dynamic analysis of inclined dam-reservoir-foundation interaction was conducted using finite difference method （FDM）. The Timoshenko beam theory and the Euler-Bemoulli beam theory were implemented to draw out governing equation of beam. The interactions between the dam and the soil were modeled by using a translational spring and a rotational spring. A Sommerfeld＇s radiation condition at the infinity boundary of the fluid domain was adopted. The effects of the reservoir bottom absorption and surface waves on the dam-reservoir-foundation interaction due to the earthquake were studied. To avoid the instability of solution, a semi-implicit scheme was used for the discretization of the governing equation of dam and an explicit scheme was used for the discretization of the governing equation of fluid. The results show that as the slope of upstream dam increases, the hydrodynamic pressure on the dam is reduced. Moreover, when the Timoshenko beam theory is used, the system response increases.

Semi-Implicit Spectral Deferred Correction Method Based on the Invariant Energy Quadratization Approach for Phase Field Problems

This paper presents a high order time discretization method by combining the semi-implicit spectral deferred correction method with energy stable linear schemes to simulate a series of phase field problems.We start with the linear scheme,which is based on the invariant energy quadratization approach and is proved to be linear unconditionally energy stable.The scheme also takes advantage of avoiding nonlinear iteration and the restriction of time step to guarantee the nonlinear system uniquely solvable.Moreover,the scheme leads to linear algebraic system to solve at each iteration,and we employ the multigrid solver to solve it efficiently.Numerical re-sults are given to illustrate that the combination of local discontinuous Galerkin(LDG)spatial discretization and the high order temporal scheme is a practical,accurate and efficient simulation tool when solving phase field problems.Namely,we can obtain high order accuracy in both time and space by solving some simple linear algebraic equations.

A semi-implicit three-step method based on SUPG finite element formulation for flow in lid driven cavities with different geometries

A numerical algorithm using a bilinear or linear finite element and semi-implicit three-step method is presented for the analysis of incompressible viscous fluid problems. The streamline upwind/Petrov-Galerkin (SUPG) stabilization scheme is used for the formulation of the Navier-Stokes equations. For the spatial discretization, the convection term is treated explicitly, while the viscous term is treated implicitly, and for the temporal discretization, a three-step method is employed. The present method is applied to simulate the lid driven cavity problems with different geometries at low and high Reynolds numbers. The results compared with other numerical experiments are found to be feasible and satisfactory.

Implementation of the moving particle semi-implicit method on GPU

The Moving Particle Semi-implicit (MPS) method performs well in simulating violent free surface flow and hence becomes popular in the area of fluid flow simulation. However, the implementations of searching neighbouring particles and solving the large sparse matrix equations (Poisson-type equation) are very time-consuming. In order to utilize the tremendous power of parallel computation of Graphics Processing Units (GPU), this study has developed a GPU-based MPS model employing the Compute Unified Device Architecture (CUDA) on NVIDIA GTX 280. The efficient neighbourhood particle searching is done through an indirect method and the Poisson-type pressure equation is solved by the Bi-Conjugate Gradient (BiCG) method. Four different optimization levels for the present general parallel GPU-based MPS model are demonstrated. In addition, the elaborate optimization of GPU code is also discussed. A benchmark problem of dam-breaking flow is simulated using both codes of the present GPU-based MPS and the original CPU-based MPS. The comparisons between them show that the GPU-based MPS model outperforms 26 times the traditional CPU model.

Control of Vortex Shedding and Drag Reduction through Dual Splitter Plates Attached to a Square Cylinder

In this paper we have made a numerical study on the control of vortex shedding and drag reduction of a cylinder by attaching thin splitter plates. The wake structure of the cylinder of square cross-section with attached splitter plates is analyzed for a range of Reynolds number, based on the incident stream and height of the cylinder, in the laminar range. The Navier-Stokes equations governing the flow are solved by the control volume method over a staggered grid arrangement. We have used the semi-implicit method for pressure-linked equation （SIMPLE） algorithm for computation. Our results show that the presence of a splitter plate upstream of the cylinder reduces the drag, but it has a small impact on the vortex shedding frequency when the plate length is beyond 1.5 time the height of the cylinder. The presence of a downstream splitter plate dampens the vortex shedding frequency. The entrainment of fluid into the inner side of the separated shear layers is obstructed by the downstream splitter plate. Our results suggest that by attaching in-line splitter plates both upstream and downstream of the cylinder, the vortex shedding can be suppressed, as well as a reduction in drag be obtained. We made a parametric study to determine the optimal length of these splitter plates so as to achieve low drag and low vortex shedding frequency.

Numerical simulations of faraday waves in cylindrical and hexagonal tanks based on MPS method

When a partially loaded liquid container vibrates along the vertical direction,the liquid inside will oscillate regularly,which is called Faraday wave.In some cases,the wave form of the Faraday wave is stable and smooth,and sometimes there is violent wave breaking and liquid splashing.In this paper,the Faraday waves inside the cylindrical tank and the hexagonal tanks are simulated by the in-house solver MLParticle-SJTU base on the moving particle semi-implicit(MPS)method.The surface tension model is used to better model the free surfaces with large deformations.Phenomena such as wave breaking and liquid splashing are well captured and simulated.The results show that the waveforms are significantly different at different excitation frequencies.And the tank shape also has an obvious effect on the waveform.

Research on Numerical Wave Tank Based on the Improved Moving-Particle Semi-Implicit Method with Large Eddy Simulation

Moving-particle semi-implicit(MPS) method is a new mesh-free numerical method based on Lagrangian particle. In this paper, MPS method is applied to the study on numerical wave tank. For the purpose of simulating numerical wave, we combine the MPS method with large eddy simulation(LES) which can simulate the turbulence in the flow. The intense pressure fluctuation is a significant shortcoming in MPS method. So, we improve the original MPS method by using a new pressure Poisson equation to ease the pressure fluctuation. Divergencefree condition representing fluid incompressible is used to calculate pressure smoothly. Then, area-time average technique is used to deal with the calculation. With these improvements, the modified MPS-LES method is applied to the simulation of numerical wave. As a contrast, we also use the original MPS-LES method to simulate the wave in a numerical wave tank. The result shows that the new method is better than the original MPS-LES method.