Eulerian method is a main numerical simulation method in elastoplastic hydrodynamics, which is suitable for the problems with multi-component and large deformation. As the feature of the problems to be simulated, such...Eulerian method is a main numerical simulation method in elastoplastic hydrodynamics, which is suitable for the problems with multi-component and large deformation. As the feature of the problems to be simulated, such as detonation and penetration, the dynamic parallel method (DPM) is designed to adjust the computational domain dynamically to get better load balance. Dynamic parallel method can be separated into two parts: one is division of initial computational domain and location of the data, the other is expansion of the computational domain and adjustment of the data location. DPM program can greatly shorten computational time and be preferable in simulating actual problems. The speedup of the DPM program is linear in parallel test. DPM can be popularized to parallel program of other multi-component high dimension Eulerian methods naturally.展开更多
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.展开更多
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.展开更多
We present a simple and efficient strategy for the acceleration of explicit Eulerian methods for multidimensional hyperbolic systems of conservation laws.The strategy is based on the Galilean invariance of dynamic equ...We present a simple and efficient strategy for the acceleration of explicit Eulerian methods for multidimensional hyperbolic systems of conservation laws.The strategy is based on the Galilean invariance of dynamic equations and optimization of the reference frame,in which the equations are numerically solved.The optimal reference frame moves(locally in time)with the average characteristic speed of the system,and,in this sense,the resulting method is quasi-Lagrangian.This leads to the acceleration of the numerical computations thanks to the optimal CFL condition and automatic adjustment of the computational domain to the evolving part of the solution.We show that our quasi-Lagrangian acceleration procedure may also reduce the numerical dissipation of the underlying Eulerian method.This leads to a significantly enhanced resolution,especially in the supersonic case.We demonstrate a great potential of the proposed method on a number of numerical examples.展开更多
In this paper, we incorporate fuzzy mathematics approach into the Eulerian method to simulate three dimensional multi-material interfaces. In particular, we propose a fuzzy interface treatment for describing interface...In this paper, we incorporate fuzzy mathematics approach into the Eulerian method to simulate three dimensional multi-material interfaces. In particular, we propose a fuzzy interface treatment for describing interfaces, designing transport plans, and computing transport quantities. Using a set of three-dimensional inviscid isothermal elastic-plastic hydrodynamic equations, we simulate shaped charge jet in different filled dynamite structures. Strain and stress have been under consideration in simulations. Numerical results demonstrate that the fuzzy interface treatment is correct and efficient for three-dimensional multi-material problems.展开更多
n the area of naval architecture and ocean engineering,the research about the underwater xplosion problem is of great significance.To achieve prolonged simulation of near-free surface underwater explosion,the underwat...n the area of naval architecture and ocean engineering,the research about the underwater xplosion problem is of great significance.To achieve prolonged simulation of near-free surface underwater explosion,the underwater explosion transient numerical model is established in this paper based on compressible Eulerian finite element method(EFEM).Compared with Geers Hunter formula,EFEM is availably validated by simulating the free-field underwater xplosion case.Then,the bubble pulsation and flow field dynamic characteristics of the cases with different underwater explosive depth are compared in this work.Lastly,the height of the water hump and the pressure of flow flied are analyzed quantitatively through the simulation results.展开更多
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 st...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.展开更多
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 Lagrang...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 passenger side airbags(PAB)are usually larger than the driver airbags.Therefore,the inflator of PAB is more powerful with high mass rate.In this paper,an Arbitrary Lagrangian-Eulerian(ALE)method based computationa...The passenger side airbags(PAB)are usually larger than the driver airbags.Therefore,the inflator of PAB is more powerful with high mass rate.In this paper,an Arbitrary Lagrangian-Eulerian(ALE)method based computational method is developed to simulate the deployment of a PAB.The tank test is used to test the property of the inflator.Through comparison of numerical and experimental results,the ALE method is validated.Based on a failed airbag test,a smaller sub-airbag is placed inside PAB to disperse the gas flow to directions which are less damaging.By applying dynamic relaxation,the initial mesh corresponding to the experimental terms is obtained.The results indicate that the interior pressure and impact force coincide with the test data,and the method in this paper is capable of capturing airbag deploying process of the PAB module accurately.展开更多
Combining the strengths of Lagrangian and Eulerian descriptions,the coupled Lagrangian–Eulerian methods play an increasingly important role in various subjects.This work reviews their development and application in o...Combining the strengths of Lagrangian and Eulerian descriptions,the coupled Lagrangian–Eulerian methods play an increasingly important role in various subjects.This work reviews their development and application in ocean engineering.Initially,we briefly outline the advantages and disadvantages of the Lagrangian and Eulerian descriptions and the main characteristics of the coupled Lagrangian–Eulerian approach.Then,following the developmental trajectory of these methods,the fundamental formulations and the frameworks of various approaches,including the arbitrary Lagrangian–Eulerian finite element method,the particle-in-cell method,the material point method,and the recently developed Lagrangian–Eulerian stabilized collocation method,are detailedly reviewed.In addition,the article reviews the research progress of these methods with applications in ocean hydrodynamics,focusing on free surface flows,numerical wave generation,wave overturning and breaking,interactions between waves and coastal structures,fluid–rigid body interactions,fluid–elastic body interactions,multiphase flow problems and visualization of ocean flows,etc.Furthermore,the latest research advancements in the numerical stability,accuracy,efficiency,and consistency of the coupled Lagrangian–Eulerian particle methods are reviewed;these advancements enable efficient and highly accurate simulation of complicated multiphysics problems in ocean and coastal engineering.By building on these works,the current challenges and future directions of the hybrid Lagrangian–Eulerian particle methods are summarized.展开更多
The interaction between geogrid and soil is crucial for the stability of geogrid-reinforced earth structure. In finite element (FE) analysis, geogrids are usually assumed as beam or truss elements, and the interaction...The interaction between geogrid and soil is crucial for the stability of geogrid-reinforced earth structure. In finite element (FE) analysis, geogrids are usually assumed as beam or truss elements, and the interaction between geogrid and soil is considered as Coulomb friction resistance, which cannot reflect the true stress and displacement developed in the reinforcement. And the traditional Lagrangian elements used to simulate soil always become highly distorted and lose accuracy in high-stress blocks. An improved geogrid model that can produce shear resistance and passive resistance and a soil model using the Eulerian technique, in combination with the coupled Eulerian-Lagrangian (CEL) method, are used to analyze the interaction between geogrid and soil of reinforced foundation test in ABAQUS. The stress in the backfill, resistance of geogrid, and settlement of foundation were computed and the results of analysis agree well with the experimental results. This simulation method is of referential value for FE analysis of reinforced earth structure.展开更多
基金Sponsored by State Key Laboratory of Computational Physics Fundation(9140C690101070C69)
文摘Eulerian method is a main numerical simulation method in elastoplastic hydrodynamics, which is suitable for the problems with multi-component and large deformation. As the feature of the problems to be simulated, such as detonation and penetration, the dynamic parallel method (DPM) is designed to adjust the computational domain dynamically to get better load balance. Dynamic parallel method can be separated into two parts: one is division of initial computational domain and location of the data, the other is expansion of the computational domain and adjustment of the data location. DPM program can greatly shorten computational time and be preferable in simulating actual problems. The speedup of the DPM program is linear in parallel test. DPM can be popularized to parallel program of other multi-component high dimension Eulerian methods naturally.
基金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.
基金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 work of A.Kurganov was supported in part by NSF grants DMS-0310585 and DMS-0610430.
文摘We present a simple and efficient strategy for the acceleration of explicit Eulerian methods for multidimensional hyperbolic systems of conservation laws.The strategy is based on the Galilean invariance of dynamic equations and optimization of the reference frame,in which the equations are numerically solved.The optimal reference frame moves(locally in time)with the average characteristic speed of the system,and,in this sense,the resulting method is quasi-Lagrangian.This leads to the acceleration of the numerical computations thanks to the optimal CFL condition and automatic adjustment of the computational domain to the evolving part of the solution.We show that our quasi-Lagrangian acceleration procedure may also reduce the numerical dissipation of the underlying Eulerian method.This leads to a significantly enhanced resolution,especially in the supersonic case.We demonstrate a great potential of the proposed method on a number of numerical examples.
基金the National Natural Science Foundation of China (Grant No.10272023).
文摘In this paper, we incorporate fuzzy mathematics approach into the Eulerian method to simulate three dimensional multi-material interfaces. In particular, we propose a fuzzy interface treatment for describing interfaces, designing transport plans, and computing transport quantities. Using a set of three-dimensional inviscid isothermal elastic-plastic hydrodynamic equations, we simulate shaped charge jet in different filled dynamite structures. Strain and stress have been under consideration in simulations. Numerical results demonstrate that the fuzzy interface treatment is correct and efficient for three-dimensional multi-material problems.
基金The authors would like to acknowledge the support of the National Natural Science Foundation of China(Grant 11672081)the Industrial Technology Development Program(Grants JCKY2018604C010 and JCKY2017604C002).Finally,Thanks for the help of Zu-Hui Li during writing the paper.
文摘n the area of naval architecture and ocean engineering,the research about the underwater xplosion problem is of great significance.To achieve prolonged simulation of near-free surface underwater explosion,the underwater explosion transient numerical model is established in this paper based on compressible Eulerian finite element method(EFEM).Compared with Geers Hunter formula,EFEM is availably validated by simulating the free-field underwater xplosion case.Then,the bubble pulsation and flow field dynamic characteristics of the cases with different underwater explosive depth are compared in this work.Lastly,the height of the water hump and the pressure of flow flied are analyzed quantitatively through the simulation results.
基金financially supported by the National Natural Science Foundation of China(Grant Nos.40906044 and 41076048)the Fundamental Research Funds for the Central Universities Project(Grant No.2011B05714)
文摘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.
基金This project was supported by the Major State Basic Research Program under Contract Grant No. G1999043803the University Fund for Mainstay Teachers of State Ministry of Education and the Opening Fund of Open Laboratory of Marine Dynamic Process and Sa
文摘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.
基金supported by the National Basic Research Program of China("973"Project)(Grant Nos.2013CB036101 and 2010CB832704)the National Natural Science Foundation of China(Grant Nos.51221961,51279030 and 51309040)
文摘The passenger side airbags(PAB)are usually larger than the driver airbags.Therefore,the inflator of PAB is more powerful with high mass rate.In this paper,an Arbitrary Lagrangian-Eulerian(ALE)method based computational method is developed to simulate the deployment of a PAB.The tank test is used to test the property of the inflator.Through comparison of numerical and experimental results,the ALE method is validated.Based on a failed airbag test,a smaller sub-airbag is placed inside PAB to disperse the gas flow to directions which are less damaging.By applying dynamic relaxation,the initial mesh corresponding to the experimental terms is obtained.The results indicate that the interior pressure and impact force coincide with the test data,and the method in this paper is capable of capturing airbag deploying process of the PAB module accurately.
基金the support received from the Laoshan Laboratory(No.LSKJ202202000)the National Natural Science Foundation of China(Grant Nos.12032002,U22A20256,and 12302253)the Natural Science Foundation of Beijing(No.L212023)for partially funding this work.
文摘Combining the strengths of Lagrangian and Eulerian descriptions,the coupled Lagrangian–Eulerian methods play an increasingly important role in various subjects.This work reviews their development and application in ocean engineering.Initially,we briefly outline the advantages and disadvantages of the Lagrangian and Eulerian descriptions and the main characteristics of the coupled Lagrangian–Eulerian approach.Then,following the developmental trajectory of these methods,the fundamental formulations and the frameworks of various approaches,including the arbitrary Lagrangian–Eulerian finite element method,the particle-in-cell method,the material point method,and the recently developed Lagrangian–Eulerian stabilized collocation method,are detailedly reviewed.In addition,the article reviews the research progress of these methods with applications in ocean hydrodynamics,focusing on free surface flows,numerical wave generation,wave overturning and breaking,interactions between waves and coastal structures,fluid–rigid body interactions,fluid–elastic body interactions,multiphase flow problems and visualization of ocean flows,etc.Furthermore,the latest research advancements in the numerical stability,accuracy,efficiency,and consistency of the coupled Lagrangian–Eulerian particle methods are reviewed;these advancements enable efficient and highly accurate simulation of complicated multiphysics problems in ocean and coastal engineering.By building on these works,the current challenges and future directions of the hybrid Lagrangian–Eulerian particle methods are summarized.
基金Supported by National Natural Science Foundation of China (No. 50678032)
文摘The interaction between geogrid and soil is crucial for the stability of geogrid-reinforced earth structure. In finite element (FE) analysis, geogrids are usually assumed as beam or truss elements, and the interaction between geogrid and soil is considered as Coulomb friction resistance, which cannot reflect the true stress and displacement developed in the reinforcement. And the traditional Lagrangian elements used to simulate soil always become highly distorted and lose accuracy in high-stress blocks. An improved geogrid model that can produce shear resistance and passive resistance and a soil model using the Eulerian technique, in combination with the coupled Eulerian-Lagrangian (CEL) method, are used to analyze the interaction between geogrid and soil of reinforced foundation test in ABAQUS. The stress in the backfill, resistance of geogrid, and settlement of foundation were computed and the results of analysis agree well with the experimental results. This simulation method is of referential value for FE analysis of reinforced earth structure.