High-order Lagrangian cell-centered conservative scheme on unstructured meshes

A high-order Lagrangian cell-centered conservative gas dynamics scheme is presented on unstructured meshes. A high-order piecewise pressure of the cell is intro- duced. With the high-order piecewise pressure of the cell, the high-order spatial discretiza- tion fluxes are constructed. The time discretization of the spatial fluxes is performed by means of the Taylor expansions of the spatial discretization fluxes. The vertex velocities are evaluated in a consistent manner due to an original solver located at the nodes by means of momentum conservation. Many numerical tests are presented to demonstrate the robustness and the accuracy of the scheme.

Lagrangian cell-centered conservative scheme

This paper presents a Lagrangian cell-centered conservative gas dynamics scheme. The piecewise constant pressures of cells arising from the current time sub-cell densities and the current time isentropic speed of sound are introduced. Multipling the initial cell density by the initial sub-cell volumes obtains the sub-cell Lagrangian masses, and dividing the masses by the current time sub-cell volumes gets the current time sub- cell densities. By the current time piecewise constant pressures of cells, a scheme that conserves the momentum and total energy is constructed. The vertex velocities and the numerical fluxes through the cell interfaces are computed in a consistent manner due to an original solver located at the nodes. The numerical tests are presented, which are representative for compressible flows and demonstrate the robustness and accuracy of the Lagrangian cell-centered conservative scheme.

An Updated Lagrangian Particle Hydrodynamics (ULPH)-NOSBPD Coupling Approach forModeling Fluid-Structure Interaction Problem

A fluid-structure interaction approach is proposed in this paper based onNon-Ordinary State-Based Peridynamics(NOSB-PD)and Updated Lagrangian Particle Hydrodynamics(ULPH)to simulate the fluid-structure interaction problem with large geometric deformation and material failure and solve the fluid-structure interaction problem of Newtonian fluid.In the coupled framework,the NOSB-PD theory describes the deformation and fracture of the solid material structure.ULPH is applied to describe the flow of Newtonian fluids due to its advantages in computational accuracy.The framework utilizes the advantages of NOSB-PD theory for solving discontinuous problems and ULPH theory for solving fluid problems,with good computational stability and robustness.A fluidstructure coupling algorithm using pressure as the transmission medium is established to deal with the fluidstructure interface.The dynamic model of solid structure and the PD-ULPH fluid-structure interaction model involving large deformation are verified by numerical simulations.The results agree with the analytical solution,the available experimental data,and other numerical results.Thus,the accuracy and effectiveness of the proposed method in solving the fluid-structure interaction problem are demonstrated.The fluid-structure interactionmodel based on ULPH and NOSB-PD established in this paper provides a new idea for the numerical solution of fluidstructure interaction and a promising approach for engineering design and experimental prediction.

Updated Lagrangian Particle Hydrodynamics (ULPH)Modeling of Natural Convection Problems

Natural convection is a heat transfer mechanism driven by temperature or density differences,leading to fluid motion without external influence.It occurs in various natural and engineering phenomena,influencing heat transfer,climate,and fluid mixing in industrial processes.This work aims to use the Updated Lagrangian Particle Hydrodynamics(ULPH)theory to address natural convection problems.The Navier-Stokes equation is discretized using second-order nonlocal differential operators,allowing a direct solution of the Laplace operator for temperature in the energy equation.Various numerical simulations,including cases such as natural convection in square cavities and two concentric cylinders,were conducted to validate the reliability of the model.The results demonstrate that the proposed model exhibits excellent accuracy and performance,providing a promising and effective numerical approach for natural convection problems.

A Modified Lagrange Method for Solving Convex Quadratic Optimization Problems

In this paper, a modified version of the Classical Lagrange Multiplier method is developed for convex quadratic optimization problems. The method, which is evolved from the first order derivative test for optimality of the Lagrangian function with respect to the primary variables of the problem, decomposes the solution process into two independent ones, in which the primary variables are solved for independently, and then the secondary variables, which are the Lagrange multipliers, are solved for, afterward. This is an innovation that leads to solving independently two simpler systems of equations involving the primary variables only, on one hand, and the secondary ones on the other. Solutions obtained for small sized problems (as preliminary test of the method) demonstrate that the new method is generally effective in producing the required solutions.

New Eulerian-Lagrangian Method for Salinity Calculation

A difference scheme in curvilinear coordinates is put forward for calculation of salinity in estuaries and coastal waters, which is based on Eulerian-Lagrangian method. It combines first-order and second-order Lagrangian interpolation to reduce numerical dispersion and oscillation. And the length of the curvilinear grid is also considered in the interpolation. Then the scheme is used in estuary, coast and ocean model, and several numerical experiments for the Yangtze Estuary and the Hangzhou Bay are conducted to test it. These experiments show that it is suitable for simulations of salinity in estuaries and coastal waters with the models using curvilinear coordinates.

The discontinuous Petrov-Galerkin method for one-dimensional compressible Euler equations in the Lagrangian coordinate

In this paper, a Petrov-Galerkin scheme named the Runge-Kutta control volume （RKCV） discontinuous finite ele- ment method is constructed to solve the one-dimensional compressible Euler equations in the Lagrangian coordinate. Its advantages include preservation of the local conservation and a high resolution. Compared with the Runge-Kutta discon- tinuous Galerkin （RKDG） method, the RKCV method is easier to implement. Moreover, the advantages of the RKCV and the Lagrangian methods are combined in the new method. Several numerical examples are given to illustrate the accuracy and the reliability of the algorithm.

An RKDG finite element method for the one-dimensional inviscid compressible gas dynamics equations in a Lagrangian coordinate

In this paper,Runge-Kutta Discontinuous Galerkin（RKDG） finite element method is presented to solve the onedimensional inviscid compressible gas dynamic equations in a Lagrangian coordinate.The equations are discretized by the DG method in space and the temporal discretization is accomplished by the total variation diminishing Runge-Kutta method.A limiter based on the characteristic field decomposition is applied to maintain stability and non-oscillatory property of the RKDG method.For multi-medium fluid simulation,the two cells adjacent to the interface are treated differently from other cells.At first,a linear Riemann solver is applied to calculate the numerical ？ux at the interface.Numerical examples show that there is some oscillation in the vicinity of the interface.Then a nonlinear Riemann solver based on the characteristic formulation of the equation and the discontinuity relations is adopted to calculate the numerical ？ux at the interface,which suppresses the oscillation successfully.Several single-medium and multi-medium fluid examples are given to demonstrate the reliability and efficiency of the algorithm.

An efficient formulation based on the Lagrangian method for contact–impact analysis of flexible multi-body system

In this paper,an efficien formulation based on the Lagrangian method is presented to investigate the contact–impact problems of f exible multi-body systems.Generally,the penalty method and the Hertz contact law are the most commonly used methods in engineering applications.However,these methods are highly dependent on various non-physical parameters,which have great effects on the simulation results.Moreover,a tremendous number of degrees of freedom in the contact–impact problems will influenc thenumericalefficien ysignificantl.Withtheconsideration of these two problems,a formulation combining the component mode synthesis method and the Lagrangian method is presented to investigate the contact–impact problems in fl xible multi-body system numerically.Meanwhile,the finit element meshing laws of the contact bodies will be studied preliminarily.A numerical example with experimental verificatio will certify the reliability of the presented formulationincontact–impactanalysis.Furthermore,aseries of numerical investigations explain how great the influenc of the finit element meshing has on the simulation results.Finally the limitations of the element size in different regions are summarized to satisfy both the accuracy and efficien y.

Hybrid N-order Lagrangian Interpolation Eulerian-Lagrangian Method for Salinity Calculation

The Eulerian？Lagrangian method（ELM） has been used by many ocean models as the solution of the advection equation,but the numerical error caused by interpolation imposes restriction on its accuracy.In the present study,hybrid N-order Lagrangian interpolation ELM（Li ELM） is put forward in which the N-order Lagrangian interpolation is used at first,then the lower order Lagrangian interpolation is applied in the points where the interpolation results are abnormally higher or lower.The calculation results of a step-shaped salinity advection model are analyzed,which show that higher order（N=3？8） Li ELM can reduce the mean numerical error of salinity calculation,but the numerical oscillation error is still significant.Even number order Li ELM makes larger numerical oscillation error than its adjacent odd number order Li ELM.Hybrid N-order Li ELM can remove numerical oscillation,and it significantly reduces the mean numerical error when N is even and the current is in fixed direction,while it makes less effect on mean numerical error when N is odd or the current direction changes periodically.Hybrid odd number order Li ELM makes less mean numerical error than its adjacent even number order Li ELM when the current is in the fixed direction,while the mean numerical error decreases as N increases when the current direction changes periodically,so odd number of N may be better for application.Among various types of Hybrid N-order Li ELM,the scheme reducing N-order directly to 1st-order may be the optimal for synthetic selection of accuracy and computational efficiency.

Simulation of sheet metal extrusion processes with Arbitrary Lagrangian-Eulerian method

An Arbitrary Lagrangian-Eulerian(ALE) method was employed to simulate the sheet metal extrusion process,aiming at avoiding mesh distortion and improving the computational accuracy.The method was implemented based on MSC/MARC by using a fractional step method,i.e.a Lagrangian step followed by an Euler step.The Lagrangian step was a pure updated Lagrangian calculation and the Euler step was performed using mesh smoothing and remapping scheme.Due to the extreme distortion of deformation domain,it was almost impossible to complete the whole simulation with only one mesh topology.Therefore,global remeshing combined with the ALE method was used in the simulation work.Based on the numerical model of the process,some deformation features of the sheet metal extrusion process,such as distribution of localized equivalent plastic strain,and shrinkage cavity,were revealed.Furthermore,the differences between conventional extrusion and sheet metal extrusion process were also analyzed.

Improved Numerical Computing Method for the 3D Tidally Induced Lagrangian Residual Current and Its Application in a Model Bay with a Longitudinal Topography

An improved method for computing the three-dimensional(3 D)first-order Lagrangian residual velocity(uL)is estab-lished.The method computes tidal body force using the harmonic constants of the zeroth-order tidal current.Compared with using the tidal-averaging method to compute the tidal body force,the proposed method filters out the clutter other than the single-frequency tidal input from the open boundary and obtains uL that is more consistent with the analytic solution.Based on the new method,uL is calculated for a wide bay with a longitudinal topography.The strength and pattern of uL are mostly determined by the parts of the tidal body force related to the vertical mixing of the Stokes’drift and the Coriolis effect,with a minor contribution from the advection effect.The geometrical shape of the bay can influence uL through the topographic gradient.The magnitude of uL increases with the increases in tidal energy input and vertical eddy viscosity and decreases in terms of the bottom friction coefficient.

An Augmented Lagrangian based Semismooth Newton Method for a Class of Bilinear Programming Problems

This paper proposes a semismooth Newton method for a class of bilinear programming problems(BLPs)based on the augmented Lagrangian,in which the BLPs are reformulated as a system of nonlinear equations with original variables and Lagrange multipliers.Without strict complementarity,the convergence of the method is studied by means of theories of semismooth analysis under the linear independence constraint qualification and strong second order sufficient condition.At last,numerical results are reported to show the performance of the proposed method.

Augmented Lagrangian Methods for Numerical Solutions to Higher Order Differential Equations

A large number of problems in engineering can be formulated as the optimization of certain functionals. In this paper, we present an algorithm that uses the augmented Lagrangian methods for finding numerical solutions to engineering problems. These engineering problems are described by differential equations with boundary values and are formulated as optimization of some functionals. The algorithm achieves its simplicity and versatility by choosing linear equality relations recursively for the augmented Lagrangian associated with an optimization problem. We demonstrate the formulation of an optimization functional for a 4th order nonlinear differential equation with boundary values. We also derive the associated augmented Lagrangian for this 4th order differential equation. Numerical test results are included that match up with well-established experimental outcomes. These numerical results indicate that the new algorithm is fully capable of producing accurate and stable solutions to differential equations.

An Optimization Model for the Strip-packing Problem and Its Augmented Lagrangian Method

This paper formulates a two-dimensional strip packing problem as a non- linear programming （NLP） problem and establishes the first-order optimality conditions for the NLP problem. A numerical algorithm for solving this NLP problem is given to find exact solutions to strip-packing problems involving up to 10 items. Approximate solutions can be found for big-sized problems by decomposing the set of items into small-sized blocks of which each block adopts the proposed numerical algorithm. Numerical results show that the approximate solutions to big-sized problems obtained by this method are superior to those by NFDH, FFDH and BFDH approaches.

Deviation of the Lagrangian particle tracing method in the evaluation of the Southern Hemisphere annual subduction rate

The classical Lagrangian particle tracing method is widely used in the evaluation of the ocean annual subduction rate.However,our analysis indicates that in addition to neglecting the effect of mixing,there are two possible deviations in the method:one is an overestimation due to not considering that the amount of subducted water at the source location may be inadequate during the late winter of the first year when the mixed layer becomes shallow;the other one is an underestimation due to the neglect of the effective subduction caused by strong vertical pumping.Quantitative analysis shows that these two deviations mainly exist in the low-latitude subduction areas of the South Pacific and South Atlantic.The two deviations have very similar distribution areas and can partially off set each other.However,the overall deviation is still large,and the maximum relative deviation ratio can reach 50%;therefore,it cannot be ignored.

Quadric SFDI for Laplacian Discretisation in Lagrangian Meshless Methods

In the Lagrangian meshless(particle)methods,such as the smoothed particle hydrodynamics(SPH),moving particle semi-implicit(MPS)method and meshless local Petrov-Galerkin method based on Rankine source solution(MLPG_R),the Laplacian discretisation is often required in order to solve the governing equations and/or estimate physical quantities(such as the viscous stresses).In some meshless applications,the Laplacians are also needed as stabilisation operators to enhance the pressure calculation.The particles in the Lagrangian methods move following the material velocity,yielding a disordered(random)particle distribution even though they may be distributed uniformly in the initial state.Different schemes have been developed for a direct estimation of second derivatives using finite difference,kernel integrations and weighted/moving least square method.Some of the schemes suffer from a poor convergent rate.Some have a better convergent rate but require inversions of high order matrices,yielding high computational costs.This paper presents a quadric semi-analytical finite-difference interpolation(QSFDI)scheme,which can achieve the same degree of the convergent rate as the best schemes available to date but requires the inversion of significant lower-order matrices,i.e.3×3 for 3D cases,compared with 6×6 or 10×10 in the schemes with the best convergent rate.Systematic patch tests have been carried out for either estimating the Laplacian of given functions or solving Poisson’s equations.The convergence,accuracy and robustness of the present schemes are compared with the existing schemes.It will show that the present scheme requires considerably less computational time to achieve the same accuracy as the best schemes available in literatures,particularly for estimating the Laplacian of given functions.

Lagrangian Relaxation Method for Multiobjective Optimization Methods: Solution Approaches

This paper introduces the Lagrangian relaxation method to solve multiobjective optimization problems. It is often required to use the appropriate technique to determine the Lagrangian multipliers in the relaxation method that leads to finding the optimal solution to the problem. Our analysis aims to find a suitable technique to generate Lagrangian multipliers, and later these multipliers are used in the relaxation method to solve Multiobjective optimization problems. We propose a search-based technique to generate Lagrange multipliers. In our paper, we choose a suitable and well-known scalarization method that transforms the original multiobjective into a scalar objective optimization problem. Later, we solve this scalar objective problem using Lagrangian relaxation techniques. We use Brute force techniques to sort optimum solutions. Finally, we analyze the results, and efficient methods are recommended.

A hybrid subcell-remapping algorithm for staggered multi-material arbitrary Lagrangian-Eulerian methods

A new flux-based hybrid subcell-remapping algorithm for staggered multimaterial arbitrary Lagrangian-Eulerian (MMALE) methods is presented. This new method is an effective generalization of the original subcell-remapping method to the multi-material regime (LOUBERE, R. and SHASHKOV,M. A subcell remapping method on staggered polygonal grids for arbitrary-Lagrangian-Eulerian methods. Journal of Computational Physics, 209, 105–138 (2005)). A complete remapping procedure of all fluid quantities is described detailedly in this paper. In the pure material regions, remapping of mass and internal energy is performed by using the original subcell-remapping method. In the regions near the material interfaces, remapping of mass and internal energy is performed with the intersection-based fluxes where intersections are performed between the swept regions and pure material polygons in the Lagrangian mesh, and an approximate approach is then introduced for constructing the subcell mass fluxes. In remapping of the subcell momentum, the mass fluxes are used to construct the momentum fluxes by multiplying a reconstructed velocity in the swept region. The nodal velocity is then conservatively recovered. Some numerical examples simulated in the full MMALE regime and several purely cyclic remapping examples are presented to prove the properties of the remapping method.

The Approximated Semi-Lagrangian WENO Methods Based on Flux Vector Splitting for Hyperbolic Conservation Laws

The paper is devised to combine the approximated semi-Lagrange weighted essentially non-oscillatory scheme and flux vector splitting. The approximated finite volume semi-Lagrange that is weighted essentially non-oscillatory scheme with Roe flux had been proposed. The methods using Roe speed to construct the flux probably generates entropy-violating solutions. More seriously, the methods maybe perform numerical instability in two-dimensional cases. A robust and simply remedy is to use a global flux splitting to substitute Roe flux. The combination is tested by several numerical examples. In addition, the comparisons of computing time and resolution between the classical weighted essentially non-oscillatory scheme (WENOJS-LF) and the semi-Lagrange weighted essentially non-oscillatory scheme (WENOEL-LF) which is presented (both combining with the flux vector splitting).