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 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.

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.

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.

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。

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.

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.

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.

CONVERGENCE RATES OF MOVING MESH RANNACHER METHODS FOR PDES OF ASIAN OPTIONS PRICING

This paper studies the convergence rates of a moving mesh implicit finite difference method with interpolation for partial differential equations （PDEs） with moving boundary arising in Asian option pricing. The moving mesh scheme is based on Rnnacher timestepping approach whose idea is running the implicit Euler schemes in the initial few steps and continuing with Crank-Nicolson schemes. With graded meshes for time direction and moving meshes for space direction, the fully discretized scheme is constructed using quadratic interpolation between two consecutive time level for the PDEs with moving boundary. The second-order convergence rates in both time and space are proved and numerical examples are carried out to confirm the theoretical results.

AnAdaptive Moving Mesh Method for Two-Dimensional Incompressible Viscous Flows

In this paper, we present an adaptive moving mesh technique for solvingthe incompressible viscous flows using the vorticity stream-function formulation. Themoving mesh strategy is based on the approach proposed by Li et al. [J. Comput. Phys.,170 (2001), pp. 562–588] to separate the mesh-moving and evolving PDE at each timestep. The Navier-Stokes equations are solved in the vorticity stream-function form bya finite-volume method in space, and the mesh-moving part is realized by solving theEuler-Lagrange equations to minimize a certain variation in conjunction with a moresophisticated monitor function. A conservative interpolation is used to redistributethe numerical solutions on the new meshes. This paper discusses the implementationof the periodic boundary conditions, where the physical domain is allowed to deformwith time while the computational domain remains fixed and regular throughout. Numericalresults demonstrate the accuracy and effectiveness of the proposed algorithm.

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.

A Numerical Study of Blowup in the Harmonic Map Heat Flow Using the MMPDE Moving Mesh Method

The numerical solution of the harmonic heat map flow problems with blowup in finite or infinite time is considered using an adaptive moving mesh method.A properly chosen monitor function is derived so that the moving mesh method can be used to simulate blowup and produce accurate blowup profiles which agree with formal asymptotic analysis.Moreover,the moving mesh method has finite time blowup when the underlying continuous problem does.In situations where the continuous problem has infinite time blowup,the moving mesh method exhibits finite time blowup with a blowup time tending to infinity as the number of mesh points increases.The inadequacy of a uniform mesh solution is clearly demonstrated.

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.

AN EFFICIENT MOVING MESH METHOD FOR A MODEL OF TURBULENT FLOW IN CIRCULAR TUBES

This paper presents an efficient moving problem with an optimal control constrained mesh method to solve a nonlinear singular condition. The physical problem is governed by a new model of turbulent flow in circular tubes proposed by Luo et al. using Prandtl＇s mixing-length theory. Our algorithm is formed by an outer iterative algorithm for handling the optimal control condition and an inner adaptive mesh redistribution algorithm for solving the singular governing equations. We discretize the nonlinear problem by using a upwinding approach, and the resulting nonlinear equations are solved by using the Newton- Raphson method. The mesh is generated and the grid points are moved by using the arc-length equidistribution principle. The numerical results demonstrate that proposed algorithm is effective in capturing the boundary layers associated with the turbulent model.

ANALYSIS OF MOVING MESH METHODS BASED ON GEOMETRICAL VARIABLES

Examines the moving mesh methods for solving one-dimensional time dependent partial differential equations. Introduction of the differential-algebraic formulations based on geometrical variables; Investigation of the well-posedness of the numerical approach; Discussion of some detailed numerical procedures.

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 Three Dimensional Gas-Kinetic Scheme with Moving Mesh for Low-Speed Viscous Flow Computations

The paper introduces the gas-kinetic scheme for three-dimensional(3D)flow simulation.First,under a unified coordinate transformation,the 3D gaskinetic BGK equation is transformed into a computational space with arbitrary mesh moving velocity.Second,based on the Chapman-Enskog expansion of the kinetic equation,a local solution of gas distribution function is constructed and used in a finite volume scheme.As a result,a Navier-Stokes flow solver is developed for the low speed flow computation with dynamical mesh movement.Several test cases are used to validate the 3D gas-kinetic method.The first example is a 3D cavity flow with up-moving boundary at Reynolds number 3200,where the periodic solutions are compared with the experimental measurements.Then,the flow evolution inside a rotating 3D cavity is simulated with the moving mesh method,where the solution differences between 2D and 3D simulation are explicitly presented.Finally,the scheme is applied to the falling plate study,where the unsteady plate tumbling motion inside water tank has been studied and compared with the experimental measurements.

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.

Balanced Monitoring of Flow Phenomena in Moving Mesh Methods

Adaptive moving mesh research usually focuses either on analytical deriva-tions for prescribed solutions or on pragmatic solvers with challenging physical appli-cations. In the latter case, the monitor functions that steer mesh adaptation are oftendefined in an ad-hoc way. In this paper we generalize our previously used moni-tor function to a balanced sum of any number of monitor components. This avoidsthe trial-and-error parameter fine-tuning that is often used in monitor functions. Thekey reason for the new balancing method is that the ratio between the maximum andaverage value of a monitor component should ideally be equal for all components.Vorticity as a monitor component is a good motivating example for this. Entropy alsoturns out to be a very informative monitor component. We incorporate the monitorfunction in an adaptive moving mesh higher-order finite volume solver with HLLCfluxes, which is suitable for nonlinear hyperbolic systems of conservation laws. Whenapplied to compressible gas flow it produces very sharp results for shocks and otherdiscontinuities. Moreover, it captures small instabilities (Richtmyer-Meshkov, Kelvin-Helmholtz). Thus showing the rich nature of the example problems and the effective-ness of the new monitor balancing.

A Moving-Mesh Finite Element Method and its Application to the Numerical Solution of Phase-Change Problems

A distributed Lagrangian moving-mesh finite element method is applied to problems involving changes of phase.The algorithm uses a distributed conservation principle to determine nodal mesh velocities,which are then used to move the nodes.The nodal values are obtained from an ALE(Arbitrary Lagrangian-Eulerian)equation,which represents a generalization of the original algorithm presented in Applied Numerical Mathematics,54:450–469(2005).Having described the details of the generalized algorithm it is validated on two test cases from the original paper and is then applied to one-phase and,for the first time,two-phase Stefan problems in one and two space dimensions,paying particular attention to the implementation of the interface boundary conditions.Results are presented to demonstrate the accuracy and the effectiveness of the method,including comparisons against analytical solutions where available.