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.

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.

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.

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.

The Rate of Convergence of Augmented Lagrangian Method for Minimax Optimization Problems with Equality Constraints

The augmented Lagrangian function and the corresponding augmented Lagrangian method are constructed for solving a class of minimax optimization problems with equality constraints.We prove that,under the linear independence constraint qualification and the second-order sufficiency optimality condition for the lower level problem and the second-order sufficiency optimality condition for the minimax problem,for a given multiplier vectorμ,the rate of convergence of the augmented Lagrangian method is linear with respect to||μu-μ^(*)||and the ratio constant is proportional to 1/c when the ratio|μ-μ^(*)||/c is small enough,where c is the penalty parameter that exceeds a threshold c_(*)>O andμ^(*)is the multiplier corresponding to a local minimizer.Moreover,we prove that the sequence of multiplier vectors generated by the augmented Lagrangian method has at least Q-linear convergence if the sequence of penalty parameters(ck)is bounded and the convergence rate is superlinear if(ck)is increasing to infinity.Finally,we use a direct way to establish the rate of convergence of the augmented Lagrangian method for the minimax problem with a quadratic objective function and linear equality constraints.

Variation of the Kuroshio intrusion pathways northeast of Taiwan using the Lagrangian method

The seasonal variations of the Kuroshio intrusion pathways northeast of Taiwan were investigated using observational data from satellite-tracked sea surface drifters and a numerical particle-tracking experiment based on a high-resolution numerical ocean model. The results of sea surface drifter data observed from 1989 to 2013 indicate that the Kuroshio surface intrusion follows two distinct pathways: one is a northwestward intrusion along the northern coast of Taiwan Island, and the other is a direct intrusion near the turn of the shelf break. The former occurs primarily in the winter, while the latter exists year round. A particle-tracking experiment in the high-resolution numerical model reproduces the two observed intrusion paths by the sea surface drifters. The three-dimensional structure of the Kuroshio intrusion is revealed by the model results. The pathways, features and possible dynamic mechanisms of the subsurface intrusion are also discussed.

A Fast Augmented Lagrangian Method for Euler’s Elastica Models

In this paper,a fast algorithm for Euler’s elastica functional is proposed,in which the Euler’s elastica functional is reformulated as a constrained minimization problem.Combining the augmented Lagrangian method and operator splitting techniques,the resulting saddle-point problem is solved by a serial of subproblems.To tackle the nonlinear constraints arising in the model,a novel fixed-point-based approach is proposed so that all the subproblems either is a linear problem or has a closed-form solution.We show the good performance of our approach in terms of speed and reliability using numerous numerical examples on synthetic,real-world and medical images for image denoising,image inpainting and image zooming problems.

Transports of air particulate matters in the atmospheric boundary layer-numerical studies using Eulerian and Lagrangian methods

Transports of air particulate matters(PM) from face sources in the atmospheric boundary layer(ABL) are investigated by the Eulerian single fluid model and the Lagrangian trajectory method,respectively.Large eddy simulation is used to simulate the fluid phase for high accuracy in both two approaches.The mean and fluctuating PM concentrations,as well as instantaneous PM distributions at different downstream and height positions,are presented.Higher mean and fluctuating particle concentrations are predicted by the Eulerian approach than the Lagrangian one.For the Lagrangian method,PM distributions cluster near the ground-wall because of the preferential dispersion of inertial particles by turbulence structures in the ABL,while it cannot be obtained by the Eulerian single fluid method,because the two-phase velocity differences are neglected in the Eulerian method.

Fast Linearized Augmented Lagrangian Method for Euler’s Elastica Model

Recently,many variational models involving high order derivatives have been widely used in image processing,because they can reduce staircase effects during noise elimination.However,it is very challenging to construct efficient algo-rithms to obtain the minimizers of original high order functionals.In this paper,we propose a new linearized augmented Lagrangian method for Euler’s elastica image denoising model.We detail the procedures of finding the saddle-points of the aug-mented Lagrangian functional.Instead of solving associated linear systems by FFTor linear iterative methods(e.g.,the Gauss-Seidel method),we adopt a linearized strat-egy to get an iteration sequence so as to reduce computational cost.In addition,we give some simple complexity analysis for the proposed method.Experimental results with comparison to the previous method are supplied to demonstrate the efficiency of the proposed method,and indicate that such a linearized augmented Lagrangian method is more suitable to deal with large-sized images.

EFFECTS OF INTEGRATIONS AND ADAPTIVITY FOR THE EULERIAN-LAGRANGIAN METHOD

This paper provides an analysis on the effects of exact and inexact integrations on stability, convergence, numerical diffusion, and numerical oscillations for the Eulerian- Lagrangian method （ELM）. In the finite element ELM, when more accurate integrations are used for the right-hand-side, less numerical diffusion is introduced and better approximation is obtained. When linear interpolation is used for numerical integrations, the resulting ELM is shown to be unconditionally stable and of first-order accuracy. When Gauss quadrature is used, conditional stability and second-order accuracy are established under some mild constraints for the convection-diffusion problems. Finally, numerical experiments demonstrate that more accurate integrations lead to better approximation, and spatial adaptivity can substantially reduce numerical oscillations and smearing that often occur in the ELM when inexact numerical integrations are used.

LAGRANGIAN METHOD USED IN THE RESEARCH OF ATMOSPHERIC CIRCULATION

Based on the conservation of entropy and potential vorticity in adiabatic atmospheric motion without the consideration of friction,calculation is made of the trajectory of a particle on an isentropic surface by use of the data of FGGE III-b.Results of several calculation schemes of the trajectory discussed show that the local data interpolation and Runge-Kutta time-integral scheme is the best.The calculated trajectory reflects the large-scale atmospheric motion only and the small-scale motion emerges as a deviation term of the calculated trajectory.And then the outbreak and propagation of planetary wave are studied by means of the deformation of a material line,with the result showing that the material line can be tracked in the trop- osphere only in a few days,beyond which the interaction between the small-scale waves and large-scale motion leads to its dramatical twisting and deformation.Therefore,the Lagrangian method is assumed to be an effective means of diagnostic research in the nonlinear intcraction in atmospheric circulation,in addition to the general study of the atmospheric circulation.

Augmented Lagrangian Methods for p-Harmonic Flows with the Generalized Penalization Terms and Application to Image Processing

In this paper,we propose a generalized penalization technique and a convex constraint minimization approach for the p-harmonic flow problem following the ideas in[Kang&March,IEEE T.Image Process.,16(2007),2251–2261].We use fast algorithms to solve the subproblems,such as the dual projection methods,primal-dual methods and augmented Lagrangian methods.With a special penalization term,some special algorithms are presented.Numerical experiments are given to demonstrate the performance of the proposed methods.We successfully show that our algorithms are effective and efficient due to two reasons:the solver for subproblem is fast in essence and there is no need to solve the subproblem accurately(even 2 inner iterations of the subproblem are enough).It is also observed that better PSNR values are produced using the new algorithms.

Augmented Lagrangian Methods for Convex Matrix Optimization Problems

In this paper,we provide some gentle introductions to the recent advance in augmented Lagrangian methods for solving large-scale convex matrix optimization problems(cMOP).Specifically,we reviewed two types of sufficient conditions for ensuring the quadratic growth conditions of a class of constrained convex matrix optimization problems regularized by nonsmooth spectral functions.Under a mild quadratic growth condition on the dual of cMOP,we further discussed the R-superlinear convergence of the Karush-Kuhn-Tucker(KKT)residuals of the sequence generated by the augmented Lagrangian methods(ALM)for solving convex matrix optimization problems.Implementation details of the ALM for solving core convex matrix optimization problems are also provided.

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.

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.

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.

SHOCK AND BOUNDARY STRUCTURE FORMATION BY SPECTRAL-LAGRANGIAN METHODS FOR THE INHOMOGENEOUS BOLTZMANN TRANSPORT EQUATION

The numerical approximation of the Spectral-Lagrangian scheme developed by the authors in [30] for a wide range of homogeneous non-linear Boltzmann type equations is extended to the space inhomogeneous case and several shock problems are benchmark. Recognizing that the Boltzmann equation is an important tool in the analysis of formation of shock and boundary layer structures, we present the computational algorithm in Section 3.3 and perform a numerical study case in shock tube geometries well modeled in for ID in x times 3D in v in Section 4. The classic Riemann problem is numerically analyzed for Knudsen numbers close to continuum. The shock tube problem of Aoki et al [2], where the wall temperature is suddenly increased or decreased, is also studied. We consider the problem of heat transfer between two parallel plates with diffusive boundary conditions for a range of Knudsen numbers from close to continuum to a highly rarefied state. Finally, the classical infinite shock tube problem that generates a non-moving shock wave is studied. The point worth noting in this example is that the flow in the final case turns from a supersonic flow to a subsonic flow across the shock.

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.

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.