In conventional modeling of a cable-pulley system, the cable must be finely meshed with Lagrangian elements for valid contact detections with pulleys, leading to extremely low efficiency. The sliding joint method base...In conventional modeling of a cable-pulley system, the cable must be finely meshed with Lagrangian elements for valid contact detections with pulleys, leading to extremely low efficiency. The sliding joint method based on the arbitrary-Lagrangian- Eulerian (ALE) formulation still lacks an efficient cable element, and in particular, modeling of friction between a sliding joint and the cable has not been studied. This paper presents efficient multi-body modeling of a cable-pulley system with friction. A variable- length cable element with a node movable along the cable, which is described with ALE, is developed to mesh the cable. A transitional cable element is then proposed to model the contact part of the cable by fixing its two nodes to the two corresponding locations of the pulley. Friction of the cable-pulley is derived as a simple law of tension decay and embedded in the multi-body system modeling. It is simplified as a generalized friction force acting only on the arc-length coordinate. This approach can use a rough mesh on the cable, and is free of contact detections, thus significantly saving computation time. Several examples are presented to validate the proposed method, and show its effectiveness in real engineering applications.展开更多
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-remapp...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.展开更多
In this paper,we present the negative norm estimates for the arbitrary Lagrangian-Eulerian discontinuous Galerkin(ALE-DG)method solving nonlinear hyperbolic equations with smooth solutions.The smoothness-increasing ac...In this paper,we present the negative norm estimates for the arbitrary Lagrangian-Eulerian discontinuous Galerkin(ALE-DG)method solving nonlinear hyperbolic equations with smooth solutions.The smoothness-increasing accuracy-conserving(SIAC)filter is a post-processing technique to enhance the accuracy of the discontinuous Galerkin(DG)solutions.This work is the essential step to extend the SIAC filter to the moving mesh for nonlinear problems.By the post-processing theory,the negative norm estimates are vital to get the superconvergence error estimates of the solutions after post-processing in the L2 norm.Although the SIAC filter has been extended to nonuniform mesh,the analysis of fil-tered solutions on the nonuniform mesh is complicated.We prove superconvergence error estimates in the negative norm for the ALE-DG method on moving meshes.The main dif-ficulties of the analysis are the terms in the ALE-DG scheme brought by the grid velocity field,and the time-dependent function space.The mapping from time-dependent cells to reference cells is very crucial in the proof.The numerical results also confirm the theoreti-cal proof.展开更多
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 mod...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.展开更多
A hybrid Lagrangian - Eulerian (HLE) method is developed for sea ice dynamics, which combines the high computational efficiency of finite difference method (FDM) with the high numerical accuracy of smoothed partic...A hybrid Lagrangian - Eulerian (HLE) method is developed for sea ice dynamics, which combines the high computational efficiency of finite difference method (FDM) with the high numerical accuracy of smoothed particle hydrodynamics (SPH). In this HLE model, the sea ice cover is represented by a group of Lagrangian ice particles with their own thicknesses and concentrations. These ice variables are interpolated to the Eularian gird nodes using the Gaussian interpolation function. The FDM is used to determine the ice velocities at Eulerian grid nodes, and the velocities of Lagrangian ice particles are interpolated from these grid velocities with the Gaussian function also. The thicknesses and concentrations of ice particles are determined based on their new locations. With the HLE numerical model, the ice ridging process in a rectangular basin is simulated, and the simulated results are validated with the analytical solution. This method is also applied to the simulation of sea ice dynamics in a vortex wind field. At last, this HLE model is applied to the Bohai Sea, and the simulated concentration, thickness and velocity match the satellite images and the field observed data well.展开更多
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 M...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.展开更多
针对海上设施的安全防护问题,提出了一种网-桁架式海上拦截装置,目的是在有效拦截来袭小艇撞击的同时尽量降低对拦截装置本身的损伤.为验证这一装置的有效性,基于显式动力学和欧拉-拉格朗日耦合方法对拦截装置拦截小艇过程进行数值模拟...针对海上设施的安全防护问题,提出了一种网-桁架式海上拦截装置,目的是在有效拦截来袭小艇撞击的同时尽量降低对拦截装置本身的损伤.为验证这一装置的有效性,基于显式动力学和欧拉-拉格朗日耦合方法对拦截装置拦截小艇过程进行数值模拟,设定小艇垂直撞击支撑柱、支撑柱间隙和45°撞击支撑柱、支撑柱间隙4种工况,每种工况包含10、20、30 m/s 3种速度.通过对这12组仿真的碰撞力、能量以及破损状态分析,全方位分析了该装置的防撞性能和对小艇的拦截效果.所有工况下该装置都能完成对小艇的拦截,说明这种装置具有优秀的拦截效果.展开更多
In this article we present a new family of high order accurate Arbitrary Lagrangian-Eulerian one-step WENO finite volume schemes for the solution of stiff hyperbolic balance laws.High order accuracy in space is obtain...In this article we present a new family of high order accurate Arbitrary Lagrangian-Eulerian one-step WENO finite volume schemes for the solution of stiff hyperbolic balance laws.High order accuracy in space is obtained with a standard WENO reconstruction algorithm and high order in time is obtained using the local space-time discontinuous Galerkinmethod recently proposed in[20].In the Lagrangian framework considered here,the local space-time DG predictor is based on a weak formulation of the governing PDE on a moving space-time element.For the spacetime basis and test functions we use Lagrange interpolation polynomials defined by tensor-product Gauss-Legendre quadrature points.The moving space-time elements are mapped to a reference element using an isoparametric approach,i.e.the spacetime mapping is defined by the same basis functions as the weak solution of the PDE.We show some computational examples in one space-dimension for non-stiff and for stiff balance laws,in particular for the Euler equations of compressible gas dynamics,for the resistive relativistic MHD equations,and for the relativistic radiation hydrodynamics equations.Numerical convergence results are presented for the stiff case up to sixth order of accuracy in space and time and for the non-stiff case up to eighth order of accuracy in space and time.展开更多
基金Project supported by the National Natural Science Foundation of China(No.11302114)the Major State Basic Research Development Program(No.2012CB821203)the Guangdong Provincial Key Laboratory Construction Project of China(No.2011A060901026)
文摘In conventional modeling of a cable-pulley system, the cable must be finely meshed with Lagrangian elements for valid contact detections with pulleys, leading to extremely low efficiency. The sliding joint method based on the arbitrary-Lagrangian- Eulerian (ALE) formulation still lacks an efficient cable element, and in particular, modeling of friction between a sliding joint and the cable has not been studied. This paper presents efficient multi-body modeling of a cable-pulley system with friction. A variable- length cable element with a node movable along the cable, which is described with ALE, is developed to mesh the cable. A transitional cable element is then proposed to model the contact part of the cable by fixing its two nodes to the two corresponding locations of the pulley. Friction of the cable-pulley is derived as a simple law of tension decay and embedded in the multi-body system modeling. It is simplified as a generalized friction force acting only on the arc-length coordinate. This approach can use a rough mesh on the cable, and is free of contact detections, thus significantly saving computation time. Several examples are presented to validate the proposed method, and show its effectiveness in real engineering applications.
基金Project supported by the China Postdoctoral Science Foundation(No.2017M610823)
文摘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 fellowship of China Postdoctoral Science Foundation,no:2020TQ0030.Y.Xu:Research supported by National Numerical Windtunnel Project NNW2019ZT4-B08+1 种基金Science Challenge Project TZZT2019-A2.3NSFC Grants 11722112,12071455.X.Li:Research supported by NSFC Grant 11801062.
文摘In this paper,we present the negative norm estimates for the arbitrary Lagrangian-Eulerian discontinuous Galerkin(ALE-DG)method solving nonlinear hyperbolic equations with smooth solutions.The smoothness-increasing accuracy-conserving(SIAC)filter is a post-processing technique to enhance the accuracy of the discontinuous Galerkin(DG)solutions.This work is the essential step to extend the SIAC filter to the moving mesh for nonlinear problems.By the post-processing theory,the negative norm estimates are vital to get the superconvergence error estimates of the solutions after post-processing in the L2 norm.Although the SIAC filter has been extended to nonuniform mesh,the analysis of fil-tered solutions on the nonuniform mesh is complicated.We prove superconvergence error estimates in the negative norm for the ALE-DG method on moving meshes.The main dif-ficulties of the analysis are the terms in the ALE-DG scheme brought by the grid velocity field,and the time-dependent function space.The mapping from time-dependent cells to reference cells is very crucial in the proof.The numerical results also confirm the theoreti-cal proof.
文摘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.
基金The study was supported by the National Natural Science Foundation of China under contract No.10772041the State 0ceamic Administration Key Laboratory for Ploar Science of China under contract No.KP 2007004.
文摘A hybrid Lagrangian - Eulerian (HLE) method is developed for sea ice dynamics, which combines the high computational efficiency of finite difference method (FDM) with the high numerical accuracy of smoothed particle hydrodynamics (SPH). In this HLE model, the sea ice cover is represented by a group of Lagrangian ice particles with their own thicknesses and concentrations. These ice variables are interpolated to the Eularian gird nodes using the Gaussian interpolation function. The FDM is used to determine the ice velocities at Eulerian grid nodes, and the velocities of Lagrangian ice particles are interpolated from these grid velocities with the Gaussian function also. The thicknesses and concentrations of ice particles are determined based on their new locations. With the HLE numerical model, the ice ridging process in a rectangular basin is simulated, and the simulated results are validated with the analytical solution. This method is also applied to the simulation of sea ice dynamics in a vortex wind field. At last, this HLE model is applied to the Bohai Sea, and the simulated concentration, thickness and velocity match the satellite images and the field observed data well.
基金Project(50505027) supported by the National Natural Science Foundation of ChinaProject(20070248056) supported by the Research Fund for the Doctoral Program of Higher Education of China
文摘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.
文摘针对海上设施的安全防护问题,提出了一种网-桁架式海上拦截装置,目的是在有效拦截来袭小艇撞击的同时尽量降低对拦截装置本身的损伤.为验证这一装置的有效性,基于显式动力学和欧拉-拉格朗日耦合方法对拦截装置拦截小艇过程进行数值模拟,设定小艇垂直撞击支撑柱、支撑柱间隙和45°撞击支撑柱、支撑柱间隙4种工况,每种工况包含10、20、30 m/s 3种速度.通过对这12组仿真的碰撞力、能量以及破损状态分析,全方位分析了该装置的防撞性能和对小艇的拦截效果.所有工况下该装置都能完成对小艇的拦截,说明这种装置具有优秀的拦截效果.
基金the European Research Council under the European Union’s Seventh Framework Programme(FP7/2007-2013)under the research project STiMulUs,ERC Grant agreement no.278267.
文摘In this article we present a new family of high order accurate Arbitrary Lagrangian-Eulerian one-step WENO finite volume schemes for the solution of stiff hyperbolic balance laws.High order accuracy in space is obtained with a standard WENO reconstruction algorithm and high order in time is obtained using the local space-time discontinuous Galerkinmethod recently proposed in[20].In the Lagrangian framework considered here,the local space-time DG predictor is based on a weak formulation of the governing PDE on a moving space-time element.For the spacetime basis and test functions we use Lagrange interpolation polynomials defined by tensor-product Gauss-Legendre quadrature points.The moving space-time elements are mapped to a reference element using an isoparametric approach,i.e.the spacetime mapping is defined by the same basis functions as the weak solution of the PDE.We show some computational examples in one space-dimension for non-stiff and for stiff balance laws,in particular for the Euler equations of compressible gas dynamics,for the resistive relativistic MHD equations,and for the relativistic radiation hydrodynamics equations.Numerical convergence results are presented for the stiff case up to sixth order of accuracy in space and time and for the non-stiff case up to eighth order of accuracy in space and time.