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.展开更多
In this paper,several arbitrary Lagrangian-Eulerian discontinuous Galerkin(ALE-DG)methods are presented for Korteweg-de Vries(KdV)type equations on moving meshes.Based on the L^(2) conservation law of KdV equations,we...In this paper,several arbitrary Lagrangian-Eulerian discontinuous Galerkin(ALE-DG)methods are presented for Korteweg-de Vries(KdV)type equations on moving meshes.Based on the L^(2) conservation law of KdV equations,we adopt the conservative and dissipative numerical fuxes for the nonlinear convection and linear dispersive terms,respectively.Thus,one conservative and three dissipative ALE-DG schemes are proposed for the equations.The invariant preserving property for the conservative scheme and the corresponding dissipative properties for the other three dissipative schemes are all presented and proved in this paper.In addition,the L^(2)-norm error estimates are also proved for two schemes,whose numerical fuxes for the linear dispersive term are both dissipative type.More precisely,when choosing the approximation space with the piecewise kth degree polynomials,the error estimate provides the kth order of convergence rate in L^(2)-norm for the scheme with the conservative numerical fuxes applied for the nonlinear convection term.Furthermore,the(k+1∕2)th order of accuracy can be proved for the ALE-DG scheme with dissipative numerical fuxes applied for the convection term.Moreover,a Hamiltonian conservative ALE-DG scheme is also presented based on the conservation of the Hamiltonian for KdV equations.Numerical examples are shown to demonstrate the accuracy and capability of the moving mesh ALE-DG methods and compare with stationary DG methods.展开更多
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 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.展开更多
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.展开更多
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 paper the numerical method for solution of an aeroelastic model describing the interactions of air flow with vocal folds is described.The flow is modelled by the incompressible Navier-Stokes equations spatiall...In this paper the numerical method for solution of an aeroelastic model describing the interactions of air flow with vocal folds is described.The flow is modelled by the incompressible Navier-Stokes equations spatially discretizedwith the aid of the stabilized finite element method.The motion of the computational domain is treated with the aid of the Arbitrary Lagrangian Eulerian method.The structure dynamics is replaced by a mechanically equivalent system with the two degrees of freedom governed by a system of ordinary differential equations and discretized in time with the aid of an implicit multistep method and strongly coupled with the flow model.The influence of inlet/outlet boundary conditions is studied and the numerical analysis is performed and compared to the related results from literature.展开更多
基金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.
基金This work was supported by the National Numerical Windtunnel Project NNW2019ZT4-B08Science Challenge Project TZZT2019-A2.3the National Natural Science Foundation of China Grant no.11871449.
文摘In this paper,several arbitrary Lagrangian-Eulerian discontinuous Galerkin(ALE-DG)methods are presented for Korteweg-de Vries(KdV)type equations on moving meshes.Based on the L^(2) conservation law of KdV equations,we adopt the conservative and dissipative numerical fuxes for the nonlinear convection and linear dispersive terms,respectively.Thus,one conservative and three dissipative ALE-DG schemes are proposed for the equations.The invariant preserving property for the conservative scheme and the corresponding dissipative properties for the other three dissipative schemes are all presented and proved in this paper.In addition,the L^(2)-norm error estimates are also proved for two schemes,whose numerical fuxes for the linear dispersive term are both dissipative type.More precisely,when choosing the approximation space with the piecewise kth degree polynomials,the error estimate provides the kth order of convergence rate in L^(2)-norm for the scheme with the conservative numerical fuxes applied for the nonlinear convection term.Furthermore,the(k+1∕2)th order of accuracy can be proved for the ALE-DG scheme with dissipative numerical fuxes applied for the convection term.Moreover,a Hamiltonian conservative ALE-DG scheme is also presented based on the conservation of the Hamiltonian for KdV equations.Numerical examples are shown to demonstrate the accuracy and capability of the moving mesh ALE-DG methods and compare with stationary DG methods.
文摘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.
基金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.
文摘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.
基金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.
基金This research was supported by the Project OC 09019“Modelling of voice production based on biomechanics”within the program COST of the Ministry of Education of the Czech Republic,under grants No.201/08/0012 and No.P101/11/0207 of the Grant Agency of the Czech Republic and the Research Plan MSM6840770003 of the Ministry of Education of the Czech Republic.
文摘In this paper the numerical method for solution of an aeroelastic model describing the interactions of air flow with vocal folds is described.The flow is modelled by the incompressible Navier-Stokes equations spatially discretizedwith the aid of the stabilized finite element method.The motion of the computational domain is treated with the aid of the Arbitrary Lagrangian Eulerian method.The structure dynamics is replaced by a mechanically equivalent system with the two degrees of freedom governed by a system of ordinary differential equations and discretized in time with the aid of an implicit multistep method and strongly coupled with the flow model.The influence of inlet/outlet boundary conditions is studied and the numerical analysis is performed and compared to the related results from literature.