A two-dimensional hybrid numerical model, FEM-LES-VOF, for free surface flows is proposed in this study, which is a combination of three-step Taylor-Galerkin finite element method, large eddy simulation with the Smago...A two-dimensional hybrid numerical model, FEM-LES-VOF, for free surface flows is proposed in this study, which is a combination of three-step Taylor-Galerkin finite element method, large eddy simulation with the Smagorinsky sub-grid model and Computational Lagrangian-Eulerian Advection Remap Volume of Fluid （CLEAR-VOF） method. The present FEM-LES-VOF model allows the fluid flows involving violent free surface and turbulence subject to complex boundary configuration to be simulated in a straightforward manner with unstructured grids in the context of finite element method. Numerical simulation of a benchmark problem of dam breaking is conducted to verify the present model. Comparisons with experimental data show that the proposed model works well and is capable of producing reliable predictions for free surface flows. Using the FEM-LES-VOF model, the free surface flow over a semi-circular obstruction is investigated. The simulation results are compared with available experimental and numerical results. Good performance of the FEM-LES-VOF model is demonstrated again. Moreover, the numerical studies show that the turbulence plays an important role in the evolution of free surface when the reflected wave propagates upstream during the fluid flow passing the submerged obstacle.展开更多
A novel VOF-type volume-tracking method for two-dimensional free-surface flows based on the unstructured triangular mesh is presented. Owing to the inherent merit of the unstructured triangular mesh in fitting curved ...A novel VOF-type volume-tracking method for two-dimensional free-surface flows based on the unstructured triangular mesh is presented. Owing to the inherent merit of the unstructured triangular mesh in fitting curved boundaries, this method can handle the free-surface problems with complex geometries accurately and directly, without introducing any complicated boundary treatment or artificial diffusion. The method solves the volume transport equation geometrically through the Modified Lagrangian-Eulerian Re-map （MLER） method, which is applied to advective fluid volumes. Moreover, the PLIC method is adopted to give a second-order reconstructed interface approximation. To validate this method, two advection tests were performed for the establishment of the accuracy and convergence rate of the solutions. Numerical results for these complex tests provide convincing evidence for the excellent solution quality and fidelity of the method.展开更多
An axially variable-length solid element with eight nodes is proposed by integrating the arbitrary Lagrangian-Eulerian (ALE) formulation and the absolute nodal coordinate formulation (ANCF). In addition to the nodal p...An axially variable-length solid element with eight nodes is proposed by integrating the arbitrary Lagrangian-Eulerian (ALE) formulation and the absolute nodal coordinate formulation (ANCF). In addition to the nodal positions and slopes of eight nodes, two material coordinates in the axial direction are used as the generalized coordinates. As a consequence, the nodes in the ALE-ANCF are not associated with any specific material points and the axial length of the solid element can be varied over time. These two material coordinates give rise to a variable mass matrix and an additional inertial force vector. Computationally efficient formulae of the additional inertial forces and elastic forces, as well as their Jacobians, are also derived. The dynamic equation of a flexible multibody system (FMBS) with variable-length bodies is presented. The maximum and minimum lengths of the boundary elements of an FMBS have to be appropriately defined to ensure accuracy and non-singularity when solving the dynamic equation. Three numerical examples of static and dynamic problems are given to validate the variable-length solid elements of ALE-ANCF and show their capability.展开更多
A computational procedure is developed to solve the problems of coupled motion of a structure and a viscous incompressible fluid. In order to incorporate the effect of the moving surface of the structure as well as th...A computational procedure is developed to solve the problems of coupled motion of a structure and a viscous incompressible fluid. In order to incorporate the effect of the moving surface of the structure as well as the free surface motion, the arbitrary Lagrangian-Eulerian formulation is employed as the basis of the finite element spatial discretization. For numerical integration in time, the fraction,step method is used. This method is useful because one can use the same linear interpolation function for both velocity and pressure. The method is applied to the nonlinear interaction of a structure and a tuned liquid damper. All computations are performed with a personal computer.展开更多
An effective computational method is developed for dynamic analysis offluid-structure interaction problems involving large-amplitude sloshing of the fluid andlarge-displacement motion of the structure. The structure i...An effective computational method is developed for dynamic analysis offluid-structure interaction problems involving large-amplitude sloshing of the fluid andlarge-displacement motion of the structure. The structure is modeled as a rigid container supportedby a system consisting of springs and dashpots. The motion of the fluid is decomposed into twoparts: the large-displacement motion with the container and the large-amplitude sloshing relative tothe container. The former is conveniently dealt with by defining a container-fixed noninertiallocal frame, while the latter is easily handled by adopting an ALE kinematical description. Thisleads to an easy and accurate treatment of both the fluid-structure interface and the fluid freesurface without producing excessive distortion of the computational mesh. The coupling between thefluid and the structure is accomplished through the coupling matrices that can be easilyestablished. Two numerical examples, including a TLD-structure system and a simplified liquid-loadedvehicle system, are presented to demonstrate the effectiveness and reliability of the proposedmethod. The present work can also be applied to simulate fluid-structure problems incorporatingmultibody systems and several fluid domains.展开更多
In this paper,we present a new vertex-centered arbitrary LagrangianEulerian(ALE)finite volume scheme for two-dimensional compressible flow.In our scheme,the momentum equation is discretized on the vertex control volum...In this paper,we present a new vertex-centered arbitrary LagrangianEulerian(ALE)finite volume scheme for two-dimensional compressible flow.In our scheme,the momentum equation is discretized on the vertex control volume,while the mass equation and the energy equation are discretized on the sub-cells which are included in the vertex control volume.We attain the average of the fluid velocity on the vertex control volume directly by solving the conservation equations.Then we can obtain the fluid velocity at vertex with the reconstructed polynomial of the velocity.This fluid velocity is chosen as the mesh velocity,which makes the mesh move in a Lagrangian manner.Two WENO(Weighted Essentially Non-Oscillatory)reconstructions for the density(the total energy)and the velocity are used to make our scheme achieve the anticipated accuracy.Compared with the general vertexcentered schemes,our scheme with the new approach for the space discretization can simulate some multi-material flows which do not involve large deformations.In addition,our scheme has good robustness,and some numerical examples are presented to demonstrate the anticipated accuracy and the good properties of our scheme.展开更多
The hemodynamic mechanism of rolling manipulation (RM) of traditional Chinese medical massage (TCMM) is investigated. An axisymmetrical nonlinear model and an arbitrary Lagrangian-Eulerian finite element method (ALE-F...The hemodynamic mechanism of rolling manipulation (RM) of traditional Chinese medical massage (TCMM) is investigated. An axisymmetrical nonlinear model and an arbitrary Lagrangian-Eulerian finite element method (ALE-FEM) with rezoning algorithm were introduced to study the viscous flow through an axisymmetrical rigid tube with axially moving stenosis to simulate the rolling manipulation. Flow rate and wall shear stress were obtained by solving complete Navier-Stokes equations numerically. The numerical results show that the stenosis moving frequency, namely the frequency of rolling manipulation, has great effect on the disturbance of flow and the wall shear stress. The stenosis coefficient, which characterizes the severity of the stenosis, another adjustable parameter in rolling manipulation, also shows the significant effect on flow rate and wall shear stress. These numerical results may provide some data that can be taken into consideration when massage is used in clinic.展开更多
Free surface flow problems involving large free motions are analysed using finite element techniques. In solving these problems an Arbitrary Lagrangian-Eulerian(ALE)kinematical description of the fluid domain is adopt...Free surface flow problems involving large free motions are analysed using finite element techniques. In solving these problems an Arbitrary Lagrangian-Eulerian(ALE)kinematical description of the fluid domain is adopted, in which the nodal points can be displaced independently of the fluid motion. A new mesh tracing method is proposed in this paper. To confirm the effectiveness of the new method, solitary wave propagation is analysed and the numerical results are compared with the analytical results. The behaviour of the viscous fluid flow with a free surface is expressed by the unsteady Navier-Stokes equation. For numerical integration in time the velocity correction fractional step method is used.展开更多
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.展开更多
This paper is concerned with the numerical simulation of the transient effect of an inertialess Boger flow past a confined circular cylinder and the comparison of predictions with particle image velocimetry (PIV) meas...This paper is concerned with the numerical simulation of the transient effect of an inertialess Boger flow past a confined circular cylinder and the comparison of predictions with particle image velocimetry (PIV) measurements given by Shiang et al.. Dynamic simulation based on the Oldroyd-B constitutive model was carried out using a Lagrangian-Eulerian algorithm. The evolution of velocity field was obtained for the flow at two Deborah (De) numbers, i.e. De = 1.2 and 3.0. At low De, the flow reached steady state rapidly, and showed a symmetric flow regime. However, at high De, the time required to reach steady flow behind the cylinder increased significantly, and the distribution of the velocity field appears to be asymmetric with respect to the stagnation line. Fairly good agreement between the numerical results and the experimental observations is reported. It can be concluded that both the experimental measurements and the present simulations indicate that the elasticity of the polymeric flow strongly affect the flow regime of viscoelastic flow around a confined cylinder.展开更多
Sheet bulk metal forming processes have been widely developed to the facilitate manufacture of complicated 3D parts. However, there is still not enough know-how available. In this paper, as one of the typical sheet bu...Sheet bulk metal forming processes have been widely developed to the facilitate manufacture of complicated 3D parts. However, there is still not enough know-how available. In this paper, as one of the typical sheet bulk metal forming processes, the sheet metal extrusion process was studied. A reasonable finite element method （FEM） model of sheet metal extrusion process taking the influence of flow-stress curve with wide range of plastic strain and ductile damage into consideration was established and simulated by an arbitrary Lagrangian-Eulerian （ALE） FEM implemented in MSC.Marc. Validated by comparing the results with experiment, some phenomenological characteristics, such as metal flow behavior, shrinkage cavity, and the influence of different combinations of diameter of punch, diameter of extrusion outlet, and diameter of pre-punched hole were analyzed and concluded, which can be used as theoretical fundamental for the design of the sheet metal extrusion process.展开更多
利用Arbitrary Lagrangian-Eulerian(ALE)有限元方法求解爆炸冲击过程中的流固耦合问题:采用ALE算法描述流体和炸药模型,采用Lagrangain方法描述舵结构模型,不同介质间的界面采用接触罚函数耦合算法。应用有限元软件ANSYS/LS-DYNA仿真...利用Arbitrary Lagrangian-Eulerian(ALE)有限元方法求解爆炸冲击过程中的流固耦合问题:采用ALE算法描述流体和炸药模型,采用Lagrangain方法描述舵结构模型,不同介质间的界面采用接触罚函数耦合算法。应用有限元软件ANSYS/LS-DYNA仿真模拟舵遭受TNT炸药爆炸冲击作用的全过程,得到舵的应力云图、位移云图,典型位置压力时间历程,加速度时间曲线等冲击响应。计算结果表明:利用A L E方法可以预估舵在水下爆炸冲击载荷作用下的损伤情况,为舵的抗冲击设计提供依据。展开更多
The Positive Crankcase Ventilation (PCV) system in a car engine is designed to lower the pressure in the crankcase, which otherwise could lead to oil leaks and seal damage. The rotation of crankshaft in the crankcase ...The Positive Crankcase Ventilation (PCV) system in a car engine is designed to lower the pressure in the crankcase, which otherwise could lead to oil leaks and seal damage. The rotation of crankshaft in the crankcase causes the churn up of oil which conducts to occurrence of oil droplets which in turn may end in the PCV exhaust air intended to be re-injected in the engine admission. The oil catch can (OCC) is a device designed to trap these oil droplets, allowing the air to escape from the crankcase with the lowest content of oil as possible and thus, reducing the generation and emission of extra pollutants during the combustion of the air-fuel mixture. The main purpose of this paper is to optimize the design of a typical OCC used in many commercial cars by varying the length of its inner tube and the relative position of the outlet from radial to tangential fitting to the can body. For this purpose, CFD parametric analysis is performed to compute a one-way coupled Lagrangian-Eulerian two-phase flow simulation of the engine oil droplets driven by the air flow stream running through the device. The study was performed using the finite volume method with second-order spatial discretization scheme on governing equations in the Solid Works-EFD CFD platform. The turbulence was modelled using the k-? model with wall functions. Numerical results have proved that maximum efficiency is obtained for the longest inner tube and the tangential position of the outlet;however, it is recommended further investigation to assess the potential erosion on the bottom of the can under such a design configuration.展开更多
In this work the feasibility of a numerical wave tank using a hybrid particle-mesh method is investigated.Based on the fluid implicit particle method（FLIP）a formulation for the hybrid method is presented for incompr...In this work the feasibility of a numerical wave tank using a hybrid particle-mesh method is investigated.Based on the fluid implicit particle method（FLIP）a formulation for the hybrid method is presented for incompressible multiphase flows involving large density jumps and wave generating boundaries.The performance of the method is assessed for a standing wave and for the generation and propagation of a solitary wave over a flat and a sloping bed.A comparison is made with results obtained with a well-established SPH package.The tests demonstrate that the method is a promising and attractive tool for simulating the nearshore propagation of waves.展开更多
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 continuou...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.展开更多
The arbitrary Lagrangian-Eulerian(ALE)method is widely used in the field of compressible multi-material and multi-phase flow problems.In order to implement the indirect ALE approach for the simulation of compressible ...The arbitrary Lagrangian-Eulerian(ALE)method is widely used in the field of compressible multi-material and multi-phase flow problems.In order to implement the indirect ALE approach for the simulation of compressible flow in the context of high order discontinuous Galerkin(DG)discretizations,we present a high order positivity-preserving DG remapping method based on a moving mesh solver in this paper.This remapping method is based on the ALE-DG method developed by Klingenberg et al.[17,18]to solve the trivial equation∂u/∂t=0 on a moving mesh,which is the old mesh before remapping at t=0 and is the new mesh after remapping at t=T.An appropriate selection of the final pseudo-time T can always satisfy the relatively mild smoothness requirement(Lipschitz continuity)on the mesh movement velocity,which guarantees the high order accuracy of the remapping procedure.We use a multi-resolution weighted essentially non-oscillatory(WENO)limiter which can keep the essentially non-oscillatory property near strong discontinuities while maintaining high order accuracy in smooth regions.We further employ an effective linear scaling limiter to preserve the positivity of the relevant physical variables without sacrificing conservation and the original high order accuracy.Numerical experiments are provided to illustrate the high order accuracy,essentially non-oscillatory performance and positivity-preserving of our remapping algorithm.In addition,the performance of the ALE simulation based on the DG framework with our remapping algorithm is examined in one-and two-dimensional Euler equations.展开更多
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.展开更多
Interaction between fuel and air in a combustion chamber is one of the main drivers of the mixing process.Experimentally,flow visualizations are limited by high droplet density in the spray.Numerically,the ability of ...Interaction between fuel and air in a combustion chamber is one of the main drivers of the mixing process.Experimentally,flow visualizations are limited by high droplet density in the spray.Numerically,the ability of large eddy simulations(LES)to resolve large scales of flow offers good perspectives on capturing flow structures issued from the interaction between the Lagrangian(fuel droplets)and Eulerian(ambient gas)phases.This study examined these interactions first during a single injection using 3D and 2D criteria for both phases.As for the 3D criteria,the spray shape was analyzed in parallel to the Q-criteria applied to the Eulerian phase,making it possible to relate the spray deformations to some specific Eulerian structures.Secondly,2D criteria were the fuel mass-fraction field and Eulerian streamlines,both taken in the mid-plane of the spray.This last analysis allows for identifying certain mechanisms involved in the Eulerian phase’s structure generation and relates it to high fuel-concentration areas in the fuel mass-fraction visualizations.展开更多
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.展开更多
基金the National Natural Science Foundation of China (Grant No. 50409015)the Program forChangjiang Scholars and Innovative Research Team inUniversity (Grant No. IRT0420) the 40th ChinaPostdoctoral Science Foundation
文摘A two-dimensional hybrid numerical model, FEM-LES-VOF, for free surface flows is proposed in this study, which is a combination of three-step Taylor-Galerkin finite element method, large eddy simulation with the Smagorinsky sub-grid model and Computational Lagrangian-Eulerian Advection Remap Volume of Fluid （CLEAR-VOF） method. The present FEM-LES-VOF model allows the fluid flows involving violent free surface and turbulence subject to complex boundary configuration to be simulated in a straightforward manner with unstructured grids in the context of finite element method. Numerical simulation of a benchmark problem of dam breaking is conducted to verify the present model. Comparisons with experimental data show that the proposed model works well and is capable of producing reliable predictions for free surface flows. Using the FEM-LES-VOF model, the free surface flow over a semi-circular obstruction is investigated. The simulation results are compared with available experimental and numerical results. Good performance of the FEM-LES-VOF model is demonstrated again. Moreover, the numerical studies show that the turbulence plays an important role in the evolution of free surface when the reflected wave propagates upstream during the fluid flow passing the submerged obstacle.
文摘A novel VOF-type volume-tracking method for two-dimensional free-surface flows based on the unstructured triangular mesh is presented. Owing to the inherent merit of the unstructured triangular mesh in fitting curved boundaries, this method can handle the free-surface problems with complex geometries accurately and directly, without introducing any complicated boundary treatment or artificial diffusion. The method solves the volume transport equation geometrically through the Modified Lagrangian-Eulerian Re-map （MLER） method, which is applied to advective fluid volumes. Moreover, the PLIC method is adopted to give a second-order reconstructed interface approximation. To validate this method, two advection tests were performed for the establishment of the accuracy and convergence rate of the solutions. Numerical results for these complex tests provide convincing evidence for the excellent solution quality and fidelity of the method.
基金the National Natural Science Foundation of China (Grants 11521062, 11722216)the 111 China Project (Grant B16003)+1 种基金Postgraduate Research and Practice Innovation Program of Jiangsu Province (Grant KYCX17_0226)China Scholarship Council.
文摘An axially variable-length solid element with eight nodes is proposed by integrating the arbitrary Lagrangian-Eulerian (ALE) formulation and the absolute nodal coordinate formulation (ANCF). In addition to the nodal positions and slopes of eight nodes, two material coordinates in the axial direction are used as the generalized coordinates. As a consequence, the nodes in the ALE-ANCF are not associated with any specific material points and the axial length of the solid element can be varied over time. These two material coordinates give rise to a variable mass matrix and an additional inertial force vector. Computationally efficient formulae of the additional inertial forces and elastic forces, as well as their Jacobians, are also derived. The dynamic equation of a flexible multibody system (FMBS) with variable-length bodies is presented. The maximum and minimum lengths of the boundary elements of an FMBS have to be appropriately defined to ensure accuracy and non-singularity when solving the dynamic equation. Three numerical examples of static and dynamic problems are given to validate the variable-length solid elements of ALE-ANCF and show their capability.
文摘A computational procedure is developed to solve the problems of coupled motion of a structure and a viscous incompressible fluid. In order to incorporate the effect of the moving surface of the structure as well as the free surface motion, the arbitrary Lagrangian-Eulerian formulation is employed as the basis of the finite element spatial discretization. For numerical integration in time, the fraction,step method is used. This method is useful because one can use the same linear interpolation function for both velocity and pressure. The method is applied to the nonlinear interaction of a structure and a tuned liquid damper. All computations are performed with a personal computer.
基金This project is supported by National 863 Hi-Tech Project Foundation (No. 2002AA411030).
文摘An effective computational method is developed for dynamic analysis offluid-structure interaction problems involving large-amplitude sloshing of the fluid andlarge-displacement motion of the structure. The structure is modeled as a rigid container supportedby a system consisting of springs and dashpots. The motion of the fluid is decomposed into twoparts: the large-displacement motion with the container and the large-amplitude sloshing relative tothe container. The former is conveniently dealt with by defining a container-fixed noninertiallocal frame, while the latter is easily handled by adopting an ALE kinematical description. Thisleads to an easy and accurate treatment of both the fluid-structure interface and the fluid freesurface without producing excessive distortion of the computational mesh. The coupling between thefluid and the structure is accomplished through the coupling matrices that can be easilyestablished. Two numerical examples, including a TLD-structure system and a simplified liquid-loadedvehicle system, are presented to demonstrate the effectiveness and reliability of the proposedmethod. The present work can also be applied to simulate fluid-structure problems incorporatingmultibody systems and several fluid domains.
基金supported by Natural Science Foundation of Guangdong province of China(Grant No.2018A030310038)National Natural Science Foundation of China(Grant Nos.11571002,11772067,11702028 and 12071046)。
文摘In this paper,we present a new vertex-centered arbitrary LagrangianEulerian(ALE)finite volume scheme for two-dimensional compressible flow.In our scheme,the momentum equation is discretized on the vertex control volume,while the mass equation and the energy equation are discretized on the sub-cells which are included in the vertex control volume.We attain the average of the fluid velocity on the vertex control volume directly by solving the conservation equations.Then we can obtain the fluid velocity at vertex with the reconstructed polynomial of the velocity.This fluid velocity is chosen as the mesh velocity,which makes the mesh move in a Lagrangian manner.Two WENO(Weighted Essentially Non-Oscillatory)reconstructions for the density(the total energy)and the velocity are used to make our scheme achieve the anticipated accuracy.Compared with the general vertexcentered schemes,our scheme with the new approach for the space discretization can simulate some multi-material flows which do not involve large deformations.In addition,our scheme has good robustness,and some numerical examples are presented to demonstrate the anticipated accuracy and the good properties of our scheme.
基金Project supported by the National Natural Science Foundation of China (No. 30070951)
文摘The hemodynamic mechanism of rolling manipulation (RM) of traditional Chinese medical massage (TCMM) is investigated. An axisymmetrical nonlinear model and an arbitrary Lagrangian-Eulerian finite element method (ALE-FEM) with rezoning algorithm were introduced to study the viscous flow through an axisymmetrical rigid tube with axially moving stenosis to simulate the rolling manipulation. Flow rate and wall shear stress were obtained by solving complete Navier-Stokes equations numerically. The numerical results show that the stenosis moving frequency, namely the frequency of rolling manipulation, has great effect on the disturbance of flow and the wall shear stress. The stenosis coefficient, which characterizes the severity of the stenosis, another adjustable parameter in rolling manipulation, also shows the significant effect on flow rate and wall shear stress. These numerical results may provide some data that can be taken into consideration when massage is used in clinic.
文摘Free surface flow problems involving large free motions are analysed using finite element techniques. In solving these problems an Arbitrary Lagrangian-Eulerian(ALE)kinematical description of the fluid domain is adopted, in which the nodal points can be displaced independently of the fluid motion. A new mesh tracing method is proposed in this paper. To confirm the effectiveness of the new method, solitary wave propagation is analysed and the numerical results are compared with the analytical results. The behaviour of the viscous fluid flow with a free surface is expressed by the unsteady Navier-Stokes equation. For numerical integration in time the velocity correction fractional step method is used.
基金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.
基金This work is supported by the National Natural Science Foundation of China (No. 29634030) and subsidized by the Special Funds for Major State Basic Research Projects (G1999064800).
文摘This paper is concerned with the numerical simulation of the transient effect of an inertialess Boger flow past a confined circular cylinder and the comparison of predictions with particle image velocimetry (PIV) measurements given by Shiang et al.. Dynamic simulation based on the Oldroyd-B constitutive model was carried out using a Lagrangian-Eulerian algorithm. The evolution of velocity field was obtained for the flow at two Deborah (De) numbers, i.e. De = 1.2 and 3.0. At low De, the flow reached steady state rapidly, and showed a symmetric flow regime. However, at high De, the time required to reach steady flow behind the cylinder increased significantly, and the distribution of the velocity field appears to be asymmetric with respect to the stagnation line. Fairly good agreement between the numerical results and the experimental observations is reported. It can be concluded that both the experimental measurements and the present simulations indicate that the elasticity of the polymeric flow strongly affect the flow regime of viscoelastic flow around a confined cylinder.
基金supported by National Science & Technology Major Project of China (No. 2009ZX04014-073)National Natural Science Foundation of China (No. 50975175)
文摘Sheet bulk metal forming processes have been widely developed to the facilitate manufacture of complicated 3D parts. However, there is still not enough know-how available. In this paper, as one of the typical sheet bulk metal forming processes, the sheet metal extrusion process was studied. A reasonable finite element method （FEM） model of sheet metal extrusion process taking the influence of flow-stress curve with wide range of plastic strain and ductile damage into consideration was established and simulated by an arbitrary Lagrangian-Eulerian （ALE） FEM implemented in MSC.Marc. Validated by comparing the results with experiment, some phenomenological characteristics, such as metal flow behavior, shrinkage cavity, and the influence of different combinations of diameter of punch, diameter of extrusion outlet, and diameter of pre-punched hole were analyzed and concluded, which can be used as theoretical fundamental for the design of the sheet metal extrusion process.
文摘利用Arbitrary Lagrangian-Eulerian(ALE)有限元方法求解爆炸冲击过程中的流固耦合问题:采用ALE算法描述流体和炸药模型,采用Lagrangain方法描述舵结构模型,不同介质间的界面采用接触罚函数耦合算法。应用有限元软件ANSYS/LS-DYNA仿真模拟舵遭受TNT炸药爆炸冲击作用的全过程,得到舵的应力云图、位移云图,典型位置压力时间历程,加速度时间曲线等冲击响应。计算结果表明:利用A L E方法可以预估舵在水下爆炸冲击载荷作用下的损伤情况,为舵的抗冲击设计提供依据。
文摘The Positive Crankcase Ventilation (PCV) system in a car engine is designed to lower the pressure in the crankcase, which otherwise could lead to oil leaks and seal damage. The rotation of crankshaft in the crankcase causes the churn up of oil which conducts to occurrence of oil droplets which in turn may end in the PCV exhaust air intended to be re-injected in the engine admission. The oil catch can (OCC) is a device designed to trap these oil droplets, allowing the air to escape from the crankcase with the lowest content of oil as possible and thus, reducing the generation and emission of extra pollutants during the combustion of the air-fuel mixture. The main purpose of this paper is to optimize the design of a typical OCC used in many commercial cars by varying the length of its inner tube and the relative position of the outlet from radial to tangential fitting to the can body. For this purpose, CFD parametric analysis is performed to compute a one-way coupled Lagrangian-Eulerian two-phase flow simulation of the engine oil droplets driven by the air flow stream running through the device. The study was performed using the finite volume method with second-order spatial discretization scheme on governing equations in the Solid Works-EFD CFD platform. The turbulence was modelled using the k-? model with wall functions. Numerical results have proved that maximum efficiency is obtained for the longest inner tube and the tangential position of the outlet;however, it is recommended further investigation to assess the potential erosion on the bottom of the can under such a design configuration.
基金The NWO(Netherlands Organisation for Scientific Research)is acknowledged for their support through the JMBC-EM Graduate Programme research grant
文摘In this work the feasibility of a numerical wave tank using a hybrid particle-mesh method is investigated.Based on the fluid implicit particle method（FLIP）a formulation for the hybrid method is presented for incompressible multiphase flows involving large density jumps and wave generating boundaries.The performance of the method is assessed for a standing wave and for the generation and propagation of a solitary wave over a flat and a sloping bed.A comparison is made with results obtained with a well-established SPH package.The tests demonstrate that the method is a promising and attractive tool for simulating the nearshore propagation of waves.
基金supported in part by a“Computational R&D in Support of Stockpile Stewardship”Grant from Lawrence Livermore National Laboratorythe National Science Foundation Grants DMS-1619892+2 种基金the Air Force Office of Scientifc Research,USAF,under Grant/contract number FA9955012-0358the Army Research Office under Grant/contract number W911NF-15-1-0517the Spanish MCINN under Project PGC2018-097565-B-I00
文摘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.
基金supported in part by NSFC grant 12031001National Key R&D Program of China No.2023YFA1009003supported in part by NSF grant DMS-2010107.
文摘The arbitrary Lagrangian-Eulerian(ALE)method is widely used in the field of compressible multi-material and multi-phase flow problems.In order to implement the indirect ALE approach for the simulation of compressible flow in the context of high order discontinuous Galerkin(DG)discretizations,we present a high order positivity-preserving DG remapping method based on a moving mesh solver in this paper.This remapping method is based on the ALE-DG method developed by Klingenberg et al.[17,18]to solve the trivial equation∂u/∂t=0 on a moving mesh,which is the old mesh before remapping at t=0 and is the new mesh after remapping at t=T.An appropriate selection of the final pseudo-time T can always satisfy the relatively mild smoothness requirement(Lipschitz continuity)on the mesh movement velocity,which guarantees the high order accuracy of the remapping procedure.We use a multi-resolution weighted essentially non-oscillatory(WENO)limiter which can keep the essentially non-oscillatory property near strong discontinuities while maintaining high order accuracy in smooth regions.We further employ an effective linear scaling limiter to preserve the positivity of the relevant physical variables without sacrificing conservation and the original high order accuracy.Numerical experiments are provided to illustrate the high order accuracy,essentially non-oscillatory performance and positivity-preserving of our remapping algorithm.In addition,the performance of the ALE simulation based on the DG framework with our remapping algorithm is examined in one-and two-dimensional Euler equations.
基金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.
文摘Interaction between fuel and air in a combustion chamber is one of the main drivers of the mixing process.Experimentally,flow visualizations are limited by high droplet density in the spray.Numerically,the ability of large eddy simulations(LES)to resolve large scales of flow offers good perspectives on capturing flow structures issued from the interaction between the Lagrangian(fuel droplets)and Eulerian(ambient gas)phases.This study examined these interactions first during a single injection using 3D and 2D criteria for both phases.As for the 3D criteria,the spray shape was analyzed in parallel to the Q-criteria applied to the Eulerian phase,making it possible to relate the spray deformations to some specific Eulerian structures.Secondly,2D criteria were the fuel mass-fraction field and Eulerian streamlines,both taken in the mid-plane of the spray.This last analysis allows for identifying certain mechanisms involved in the Eulerian phase’s structure generation and relates it to high fuel-concentration areas in the fuel mass-fraction visualizations.
基金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.
