TWO IMPROVED ALGORITHMS AND IMPLEMENTATION FOR A SINGULARLY PERTURBED PROBLEM ON MOVING MESHES

This paper applies a difference scheme to a singularly perturbed problem. The authors provide two algorithms on moving mesh methods by using Richardson extrapolation which can improve the accuracy of numerical solution. In traditional algorithms of moving meshes, the initial mesh is a uniform mesh. The authors change it to Bakhvalov-Shishkin mesh, and prove that it improves efficiency by numerical experiments. Finally, the results of the two algorithms are analyzed.

A Study on CFL Conditions for the DG Solution of Conservation Laws on Adaptive Moving Meshes

The selection of time step plays a crucial role in improving stability and efficiency in the Discontinuous Galerkin(DG)solution of hyperbolic conservation laws on adaptive moving meshes that typically employs explicit stepping.A commonly used selection of time step is a direct extension based on Courant-Friedrichs-Levy(CFL)conditions established for fixed and uniform meshes.In this work,we provide a mathematical justification for those time step selection strategies used in practical adaptive DG computations.A stability analysis is presented for a moving mesh DG method for linear scalar conservation laws.Based on the analysis,a new selection strategy of the time step is proposed,which takes into consideration the coupling of theα-function(that is related to the eigenvalues of the Jacobian matrix of the flux and the mesh movement velocity)and the heights of the mesh elements.The analysis also suggests several stable combinations of the choices of theα-function in the numerical scheme and in the time step selection.Numerical results obtained with a moving mesh DG method for Burgers’and Euler equations are presented.For comparison purpose,numerical results obtained with an error-based time step-size selection strategy are also given。

Adaptive Moving Mesh Central-Upwind Schemes for Hyperbolic System of PDEs:Applications to Compressible Euler Equations and Granular Hydrodynamics

We introduce adaptive moving mesh central-upwind schemes for one-and two-dimensional hyperbolic systems of conservation and balance laws.The proposed methods consist of three steps.First,the solution is evolved by solving the studied system by the second-order semi-discrete central-upwind scheme on either the one-dimensional nonuniform grid or the two-dimensional structured quadrilateral mesh.When the evolution step is complete,the grid points are redistributed according to the moving mesh differential equation.Finally,the evolved solution is projected onto the new mesh in a conservative manner.The resulting adaptive moving mesh methods are applied to the one-and two-dimensional Euler equations of gas dynamics and granular hydrodynamics systems.Our numerical results demonstrate that in both cases,the adaptive moving mesh central-upwind schemes outperform their uniform mesh counterparts.

An Improved Moving Mesh Algorithm

We consider an iterative algorithm of mesh optimization for finite element solution, and give an improved moving mesh strategy that reduces rapidly the complexity and cost of solving variational problems. A numerical result is presented for a 2-dimensional problem by the improved algorithm.

An Analysis of the Dynamic Behavior of Damaged Reinforced Concrete Bridges under Moving Vehicle Loads by Using the Moving Mesh Technique

This work proposes a numerical investigation on the effects of damage on the structural response of Reinforced Concrete(RC)bridge structures commonly adopted in highway and railway networks.An effective three-dimensional FE-based numerical model is developed to analyze the bridge’s structural response under several damage scenarios,including the effects of moving vehicle loads.In particular,the longitudinal and transversal beams are modeled through solid finite elements,while horizontal slabs are made of shell elements.Damage phenomena are also incorporated in the numerical model according to a smeared approach consistent with Continuum Damage Mechanics(CDM).In such a context,the proposed method utilizes an advanced and efficient computational strategy for reproducing Vehicle-Bridge Interaction(VBI)effects based on a moving mesh technique consistent with the Arbitrary Lagrangian-Eulerian(ALE)formulation.The proposed model adopts a moving mesh interface for tracing the positions of the contact points between the vehicle’s wheels and the bridge slabs.Such modeling strategy avoids using extremely refined discretization for structural members,thus drastically reducing computational efforts.Vibrational analyses in terms of damage scenarios are presented to verify how the presence of damage affects the natural frequencies of the structural system.In addition,a comprehensive investigation regarding the response of the bridge under moving vehicles is developed,also providing results in terms of Dynamic Amplification Factor(DAFs)for typical design bridge variables.

An Unconventional Divergence Preserving Finite-Volume Discretization of Lagrangian Ideal MHD

We construct an unconventional divergence preserving discretization of updated Lagrangian ideal magnetohydrodynamics(MHD)over simplicial grids.The cell-centered finite-volume(FV)method employed to discretize the conservation laws of volume,momentum,and total energy is rigorously the same as the one developed to simulate hyperelasticity equations.By construction this moving mesh method ensures the compatibility between the mesh displacement and the approximation of the volume flux by means of the nodal velocity and the attached unit corner normal vector which is nothing but the partial derivative of the cell volume with respect to the node coordinate under consideration.This is precisely the definition of the compatibility with the Geometrical Conservation Law which is the cornerstone of any proper multi-dimensional moving mesh FV discretization.The momentum and the total energy fluxes are approximated utilizing the partition of cell faces into sub-faces and the concept of sub-face force which is the traction force attached to each sub-face impinging at a node.We observe that the time evolution of the magnetic field might be simply expressed in terms of the deformation gradient which characterizes the Lagrange-to-Euler mapping.In this framework,the divergence of the magnetic field is conserved with respect to time thanks to the Piola formula.Therefore,we solve the fully compatible updated Lagrangian discretization of the deformation gradient tensor for updating in a simple manner the cell-centered value of the magnetic field.Finally,the sub-face traction force is expressed in terms of the nodal velocity to ensure a semi-discrete entropy inequality within each cell.The conservation of momentum and total energy is recovered prescribing the balance of all the sub-face forces attached to the sub-faces impinging at a given node.This balance corresponds to a vectorial system satisfied by the nodal velocity.It always admits a unique solution which provides the nodal velocity.The robustness and the accuracy of this unconventional FV scheme have been demonstrated by employing various representative test cases.Finally,it is worth emphasizing that once you have an updated Lagrangian code for solving hyperelasticity you also get an almost free updated Lagrangian code for solving ideal MHD ensuring exactly the compatibility with the involution constraint for the magnetic field at the discrete level.

A Unique Modelling Strategy to Dynamically Simulate the Performance of a Lobe Pump for Industrial Applications

The performance of a newly designed tri-lobe industrial lobe pump of high capacity is simulated by using commercial CFD solver Ansys Fluent. A combination of user-defined-functions and meshing strategies is employed to capture the rotation of the lobes. The numerical model is validated by comparing the simulated results with the literature values. The processes of suction, displacement, compression and exhaust are accurately captured in the transient simulation. The fluid pressure value remains in the range of inlet pressure value till the processes of suction and displacement are over. The instantaneous process of compression is accurately captured in the simulation. The movement of a particular working chamber is traced along the gradual degree of lobe’s rotation. At five different degrees of lobe’s rotation, pressure contour plots are reported which clearly shows the pressure values inside the working chamber. Each pressure value inside the working chamber conforms to the particular process in which the working chamber is operating. Finally, the power requirement at the shaft of rotation is estimated from the simulated values. The estimated value of power requirement is 3.61 BHP FHP whereas the same calculated theoretically is 3 BHP FHP. The discrepancy is attributed to the assumption of symmetry of blower along the thickness.

Second-Order Invariant Domain Preserving ALE Approximation of Euler Equations

An invariant domain preserving arbitrary Lagrangian-Eulerian method for solving non-linear hyperbolic systems is developed.The numerical scheme is explicit in time and the approximation in space is done with continuous finite elements.The method is made invar-iant domain preserving for the Euler equations using convex limiting and is tested on vari-ous benchmarks.

Single-Step Arbitrary Lagrangian-Eulerian Discontinuous Galerkin Method for 1-D Euler Equations

We propose an explicit,single-step discontinuous Galerkin method on moving grids using the arbitrary Lagrangian-Eulerian approach for one-dimensional Euler equations.The grid is moved with the local fluid velocity modified by some smoothing,which is found to con-siderably reduce the numerical dissipation introduced by Riemann solvers.The scheme preserves constant states for any mesh motion and we also study its positivity preservation property.Local grid refinement and coarsening are performed to maintain the mesh qual-ity and avoid the appearance of very small or large cells.Second,higher order methods are developed and several test cases are provided to demonstrate the accuracy of the proposed scheme.

Simulation of Hydrodynamic Performance of Drag and Double Reverse Propeller Podded Propulsors

The unsteady performance of drag and double reverse propeller podded propulsors in open water was numerically simulated using a computational fluid dynamics （CFD） method. A moving mesh method was used to more realistically simulate propulsor working conditions, and the thrust, torque, and lateral force coefficients of both propulsors were compared and analyzed. Forces acting on different parts of the propulsors along with the flow field distribution of steady and unsteady results at different advance coefficients were compared. Moreover, the change of the lateral force and the difference between the abovementioned two methods were mainly analyzed. It was shown that the thrust and torque results of both methods were similar, with the lateral force results having the highest deviation

Novel energy dissipative method on the adaptive spatial discretization for the Allen–Cahn equation

We propose a novel energy dissipative method for the Allen–Cahn equation on nonuniform grids.For spatial discretization,the classical central difference method is utilized,while the average vector field method is applied for time discretization.Compared with the average vector field method on the uniform mesh,the proposed method can involve fewer grid points and achieve better numerical performance over long time simulation.This is due to the moving mesh method,which can concentrate the grid points more densely where the solution changes drastically.Numerical experiments are provided to illustrate the advantages of the proposed concrete adaptive energy dissipative scheme under large time and space steps over a long time.

Analytical and Numerical Investigation for the DMBBM Equation

The nonlinear dispersive modified Benjamin-Bona-Mahony(DMBBM)equation is solved numerically using adaptive moving mesh PDEs(MMPDEs)method.Indeed,the exact solution of the DMBBM equation is obtained by using the extended Jacobian elliptic function expansion method.The current methods give a wider applicability for handling nonlinear wave equations in engineering and mathematical physics.The adaptive moving mesh method is compared with exact solution by numerical examples,where the explicit solutions are known.The numerical results verify the accuracy of the proposed method.

A Transient Multidimensional CFD Approach to the Analysis of a Control Valve Using Non-Newtonian Fluids

In this paper the flow through a control directional valve is studied by means of a CFD （computational fluid-dynamics） analysis under transient operating conditions. The mesh motion is resolved on a time basis as a function of the external actuation system In the analysis, an open source fluid-dynamics code is used and both cavitation and turbulence are accounted for in the modeling. Moreover, the numerical model of the working fluid is modified in order to account also for the non-Newtonian fluids. The effects of the shear rate on the shear stress are accounted for both by using experimental measurements and correlations available in literature, such as the Herschel-Bulkley model. The analysis determines the performance of the control directional valve under different operating conditions when using either Newtonian or non-Newtonian fluids. In particular, the discharge coefficient, the recirculating regions, the flow acceleration angle and the pressure and velocity fields are investigated.

Transient Flow in Rapidly Filling Air-Entrapped Pipelines with Moving Boundaries

A mathematical model is presented for transient flow in a rapidly filling pipeline with an entrapped air pocket. The influence of transient shear stress between the pipe wall and the flowing fluid is taken into account. A coordinate transformation technique is employed to generate adaptive moving meshes for the multiphase flow system as images of the time-independent computational meshes in auxiliary domains. The method of characteristics is used to reduce the coupled nonlinear hyperbolic partial differential equations governing the motion of the filling fluid, entrapped air, and blocking fluid to ordinary differential equations. Numerical solution of resulting equations shows that the transient shear stresses have only a small damping effect on the pressure fluctuations. The peak pressure in the entrapped air pocket decreases significantly with increasing initial entrapped air volume, but decreases slightly with increasing initial entrapped air pressure.

The Error-Minimization-Based Strategy for Moving Mesh Methods

The typical elements in a numerical simulation of fluid flow using moving meshes are a time integration scheme,a rezone method in which a new mesh is defined,and a remapping(conservative interpolation)in which a solution is transferred to the new mesh.The objective of the rezone method is to move the computational mesh to improve the robustness,accuracy and eventually efficiency of the simulation.In this paper,we consider the onedimensional viscous Burgers’equation and describe a new rezone strategy which minimizes the L2 norm of error and maintains mesh smoothness.The efficiency of the proposed method is demonstrated with numerical examples.

Delaunay graph-based moving mesh method with damping functions

The fluid–structure interaction and aerodynamic shape optimization usually involve the moving or deforming boundaries, thus the dynamic mesh techniques are the key techniques to cope with such deformation. A novel dynamic mesh method was developed based on the Delaunay graph in this paper. According to the Delaunay graph, the mesh points were divided into groups. In each group, a factor ranging from 0 to 1 was calculated based on the area/volume ratio. By introducing a proper function for this factor, this method can control the mesh quality with high efficiency. Several test cases were compared with other dynamic mesh methods regarding mesh quality and CPU time, such as radial basis function method and Delaunay graph mapping method.

AnAdaptive Moving Mesh Method for Two-Dimensional Relativistic Hydrodynamics

This paper extends the adaptive moving mesh method developed by Tang and Tang[36]to two-dimensional(2D)relativistic hydrodynamic(RHD)equations.The algorithm consists of two“independent”parts:the time evolution of the RHD equations and the(static)mesh iteration redistribution.In the first part,the RHD equations are discretized by using a high resolution finite volume scheme on the fixed but nonuniform meshes without the full characteristic decomposition of the governing equations.The second part is an iterative procedure.In each iteration,the mesh points are first redistributed,and then the cell averages of the conservative variables are remapped onto the new mesh in a conservative way.Several numerical examples are given to demonstrate the accuracy and effectiveness of the proposed method.

ON A MOVING MESH METHOD FOR SOLVING PARTIAL INTEGRO-DIFFERENTIAL EQUATIONS

This paper develops and analyzes a moving mesh finite difference method for solving partial integro-differential equations. First, the time-dependent mapping of the coordinate transformation is approximated by a a piecewise linear function in time. Then, piecewise quadratic polynomial in space and an efficient method to discretize the memory term of the equation is designed using the moving mesh approach. In each time slice, a simple piecewise constant approximation of the integrand is used, and thus a quadrature is constructed for the memory term. The central finite difference scheme for space and the backward Euler scheme for time are used. The paper proves that the accumulation of the quadrature error is uniformly bounded and that the convergence of the method is second order in space and first order in time. Numerical experiments are carried out to confirm the theoretical predictions.

A Well-Balanced Positivity-Preserving Quasi-Lagrange Moving Mesh DG Method for the Shallow Water Equations

A high-order, well-balanced, positivity-preserving quasi-Lagrange movingmesh DG method is presented for the shallow water equations with non-flat bottomtopography. The well-balance property is crucial to the ability of a scheme to simulate perturbation waves over the lake-at-rest steady state such as waves on a lake ortsunami waves in the deep ocean. The method combines a quasi-Lagrange movingmesh DG method, a hydrostatic reconstruction technique, and a change of unknownvariables. The strategies in the use of slope limiting, positivity-preservation limiting,and change of variables to ensure the well-balance and positivity-preserving properties are discussed. Compared to rezoning-type methods, the current method treatsmesh movement continuously in time and has the advantages that it does not need tointerpolate flow variables from the old mesh to the new one and places no constraintfor the choice of a update scheme for the bottom topography on the new mesh. A selection of one- and two-dimensional examples are presented to demonstrate the wellbalance property, positivity preservation, and high-order accuracy of the method andits ability to adapt the mesh according to features in the flow and bottom topography.

Lagrangian-Eulerian One-Step WENO Finite Volume Schemes on Unstructured Triangular Meshes

In this article we present a new class of high order accurate Arbitrary-Eulerian-Lagrangian(ALE)one-step WENO finite volume schemes for solving nonlinear hyperbolic systems of conservation laws on moving two dimensional unstructured triangular meshes.A WENO reconstruction algorithm is used to achieve high order accuracy in space and a high order one-step time discretization is achieved by using the local space-time Galerkin predictor proposed in[25].For that purpose,a new element-local weak formulation of the governing PDE is adopted on moving space-time elements.The space-time basis and test functions are obtained considering Lagrange interpolation polynomials passing through a predefined set of nodes.Moreover,a polynomial mapping defined by the same local space-time basis functions as the weak solution of the PDE is used to map the moving physical space-time element onto a space-time reference element.To maintain algorithmic simplicity,the final ALE one-step finite volume scheme uses moving triangular meshes with straight edges.This is possible in the ALE framework,which allows a local mesh velocity that is different from the local fluid velocity.We present numerical convergence rates for the schemes presented in this paper up to sixth order of accuracy in space and time and show some classical numerical test problems for the two-dimensional Euler equations of compressible gas dynamics.