Considering droplet phenomena at low Mach numbers,large differences in the magnitude of the occurring characteristic waves are presented.As acoustic phenomena often play a minor role in such applications,classical exp...Considering droplet phenomena at low Mach numbers,large differences in the magnitude of the occurring characteristic waves are presented.As acoustic phenomena often play a minor role in such applications,classical explicit schemes which resolve these waves suffer from a very restrictive timestep restriction.In this work,a novel scheme based on a specific level set ghost fluid method and an implicit-explicit(IMEX)flux splitting is proposed to overcome this timestep restriction.A fully implicit narrow band around the sharp phase interface is combined with a splitting of the convective and acoustic phenomena away from the interface.In this part of the domain,the IMEX Runge-Kutta time discretization and the high order discontinuous Galerkin spectral element method are applied to achieve high accuracies in the bulk phases.It is shown that for low Mach numbers a significant gain in computational time can be achieved compared to a fully explicit method.Applica-tions to typical droplet dynamic phenomena validate the proposed method and illustrate its capabilities.展开更多
The ghost fluid method and the level set function is used to simulate the Richtmyer Meshkov instability of multi material flow interface in two dimensions. To keep the density profile from smearing out, the isobaric f...The ghost fluid method and the level set function is used to simulate the Richtmyer Meshkov instability of multi material flow interface in two dimensions. To keep the density profile from smearing out, the isobaric fix technique is used in numerical algorithms.展开更多
In this work,the modified ghost fluid method is developed to deal with 2D compressible fluid interacting with elastic solid in an Euler-Lagrange coupled system.In applying the modified Ghost Fluid Method to treat the ...In this work,the modified ghost fluid method is developed to deal with 2D compressible fluid interacting with elastic solid in an Euler-Lagrange coupled system.In applying the modified Ghost Fluid Method to treat the fluid-elastic solid coupling,the Navier equations for elastic solid are cast into a system similar to the Euler equations but in Lagrangian coordinates.Furthermore,to take into account the influence of material deformation and nonlinear wave interaction at the interface,an Euler-Lagrange Riemann problem is constructed and solved approximately along the normal direction of the interface to predict the interfacial status and then define the ghost fluid and ghost solid states.Numerical tests are presented to verify the resultant method.展开更多
In this work,the modified Ghost Fluid Method is further developed to ap-ply to compressible fluid coupled to deformable structure,where the pressure in the structure or flow can vary from an initial extremely high mag...In this work,the modified Ghost Fluid Method is further developed to ap-ply to compressible fluid coupled to deformable structure,where the pressure in the structure or flow can vary from an initial extremely high magnitude(such that the solid medium can be under plastic compression)to a subsequently very low quantity(so that cavitation can occur in the fluid).New techniques are also developed in the definition of the ghost fluid status when the structure is under plastic deformation or when the flow is under cavitation next to the structure.Numerical results show that the improved MGFM for treatment of the fluid-deformable structure coupling works efficiently for all pressure ranges and is capable of simulating cavitation evolution and cavitation re-loading in conjunction with the employment of the isentropic one-fluid cavitation model.展开更多
Recent years the modify ghost fluid method （MGFM） and the real ghost fluid method （RGFM） based on Riemann problem have been developed for multimedium compressible flows. According to authors, these methods have on...Recent years the modify ghost fluid method （MGFM） and the real ghost fluid method （RGFM） based on Riemann problem have been developed for multimedium compressible flows. According to authors, these methods have only been used with the level set technique to track the interface. In this paper, we combine the MCFM and the RGFM respectively with front tracking method, for which the fluid interfaces are explicitly tracked by connected points. The method is tested with some one-dimensional problems, and its applicability is also studied. Furthermore, in order to capture the interface more accurately, especially for strong shock impacting on interface, a shock monitor is proposed to determine the initial states of the Riemann problem. The present method is applied to various one- dimensional problems involving strong shock-interface interaction. An extension of the present method to two dimension is also introduced and preliminary results are given.展开更多
The original ghost fluid method (GFM) developed in [13] and the modifiedGFM (MGFM) in [26] have provided a simple and yet flexible way to treat twomediumflow problems. The original GFM and MGFM make the material inter...The original ghost fluid method (GFM) developed in [13] and the modifiedGFM (MGFM) in [26] have provided a simple and yet flexible way to treat twomediumflow problems. The original GFM and MGFM make the material interface"invisible" during computations and the calculations are carried out as for a singlemedium such that its extension to multi-dimensions becomes fairly straightforward.The Runge-Kutta discontinuous Galerkin (RKDG) method for solving hyperbolic conservationlaws is a high order accurate finite element method employing the usefulfeatures from high resolution finite volume schemes, such as the exact or approximateRiemann solvers, TVD Runge-Kutta time discretizations, and limiters. In this paper,we investigate using RKDG finite element methods for two-medium flow simulationsin one and two dimensions in which the moving material interfaces is treated via nonconservativemethods based on the original GFM and MGFM. Numerical results forboth gas-gas and gas-water flows are provided to show the characteristic behaviors ofthese combinations.展开更多
The modified ghost fluid method(MGFM)provides a robust and efficient interface treatment for various multi-medium flow simulations and some particular fluid-structure interaction(FSI)simulations.However,this methodolo...The modified ghost fluid method(MGFM)provides a robust and efficient interface treatment for various multi-medium flow simulations and some particular fluid-structure interaction(FSI)simulations.However,this methodology for one specific class of FSI problems,where the structure is plate,remains to be developed.This work is devoted to extending the MGFM to treat compressible fluid coupled with a thin elastic plate.In order to take into account the influence of simultaneous interaction at the interface,a fluid-plate coupling system is constructed at each time step and solved approximately to predict the interfacial states.Then,ghost fluid states and plate load can be defined by utilizing the obtained interfacial states.A type of acceleration strategy in the coupling process is presented to pursue higher efficiency.Several one-dimensional examples are used to highlight the utility of this method over looselycoupled method and validate the acceleration techniques.Especially,this method is applied to compute the underwater explosions(UNDEX)near thin elastic plates.Evolution of strong shock impacting on the thin elastic plate and dynamic response of the plate are investigated.Numerical results disclose that this methodology for treatment of the fluid-plate coupling indeed works conveniently and accurately for different structural flexibilities and is capable of efficiently simulating the processes of UNDEX with the employment of the acceleration strategy.展开更多
A conservative modification to the ghost fluid method(GFM)is developed for compressible multiphase flows.The motivation is to eliminate or reduce the conservation error of the GFM without affecting its performance.We ...A conservative modification to the ghost fluid method(GFM)is developed for compressible multiphase flows.The motivation is to eliminate or reduce the conservation error of the GFM without affecting its performance.We track the conservative variables near the material interface and use this information to modify the numerical solution for an interfacing cell when the interface has passed the cell.The modification procedure can be used on the GFM with any base schemes.In this paper we use the fifth order finite difference WENO scheme for the spatial discretization and the third order TVD Runge-Kutta method for the time discretization.The level set method is used to capture the interface.Numerical experiments show that the method is at least mass and momentum conservative and is in general comparable in numerical resolution with the original GFM.展开更多
The modified ghost fluid method(MGFM)has been shown to be robust and efficient when being applied to multi-medium compressible flows.In this paper,we rigorously analyze the optimal error estimation of the MGFM when it...The modified ghost fluid method(MGFM)has been shown to be robust and efficient when being applied to multi-medium compressible flows.In this paper,we rigorously analyze the optimal error estimation of the MGFM when it is applied to the multi-fluid Riemann problem.By analyzing the properties of the MGFM and the approximate Riemann problem solver(ARPS),we show that the interfacial status provided by the MGFM can achieve“third-order accuracy”in the sense of comparing to the exact solution of the Riemann problem,regardless of the solution type.In addition,our analysis further reveals that the ARPS based on a doubled shock structure in the MGFM is suitable for almost any conditions for predicting the interfacial status,and that the“natural”approach of“third-order accuracy”is practically less useful.Various examples are presented to validate the conclusions made.展开更多
Two interface capturing methods are studied for multi fluid flows, governed by the stiffened gas equation of state. The mixture type interface capturing algorithm uses a simple volume fraction model Euler equations wr...Two interface capturing methods are studied for multi fluid flows, governed by the stiffened gas equation of state. The mixture type interface capturing algorithm uses a simple volume fraction model Euler equations written in a quasi conservative form, which is solved by a standard high resolution piecewise parabolic method (PPM) with multi fluid Riemann solver. The level set interface capturing method uses a narrow band ghost fluid method (GFM) with no numerical smearing. Several examples are presented and compared for one and two dimensions, which show the feasibility of the two methods applied to various multi fluid problems.展开更多
In the present paper, a three-dimensional （3D） Eulerian technique for the 3D numerical simulation of high-velocity impact problems is proposed. In the Eulerian framework, a complete 3D conservation element and solut...In the present paper, a three-dimensional （3D） Eulerian technique for the 3D numerical simulation of high-velocity impact problems is proposed. In the Eulerian framework, a complete 3D conservation element and solution element scheme for conservative hyperbolic governing equations with source terms is given. A modified ghost fluid method is proposed for the treatment of the boundary conditions. Numerical simulations of the Taylor bar problem and the ricochet phenomenon of a sphere impacting a plate target at an angle of 60~ are carried out. The numerical results are in good agreement with the corresponding experimental observations. It is proved that our computational technique is feasible for analyzing 3D high-velocity impact problems.展开更多
Enhancement of two fluid mixing was numerically studied by tracking the multi fluid interfaces. Level set equations were used to capture the interfaces, and flow field was obtained by upwind TVD scheme to solve 2D Eul...Enhancement of two fluid mixing was numerically studied by tracking the multi fluid interfaces. Level set equations were used to capture the interfaces, and flow field was obtained by upwind TVD scheme to solve 2D Eulerian equations. The boundary conditions at interface of two fluids are determined by Ghost fluid method (GFM). The distributions of fluid parameters, such as pressure and density, were got at different time steps. The results show that the method presented in this paper can track the density discontinuity perfectly. Superior to previous results, the density discontinuity remains sharper. Also, the mixing of fluids can be greatly enhanced by setting disturbances along the initial fluid interfaces.展开更多
In the paper, the numerical simulation of interface problems for multiple material fluids is studied. The level set function is designed to capture the location of the material interface. For multi-dimensional and mul...In the paper, the numerical simulation of interface problems for multiple material fluids is studied. The level set function is designed to capture the location of the material interface. For multi-dimensional and multi-material fluids, the modified ghost fluid method needs a Riemann solution to renew the variable states near the interface. Here we present a new convenient and effective algorithm for solving the Riemann problem in the normal direction. The extrapolated variables are populated by Taylor series expansions in the direction. The anti-diffusive high order WENO difference scheme with the limiter is adopted for the numerical simulation. Finally we implement a series of numerical experiments of multi-material flows. The obtained results are satisfying, compared to those by other methods.展开更多
基金support provided by the Deutsche Forschun-gsgemeinschaft(DFG,German Research Foundation)through the project GRK 2160/1“Droplet Interaction Technologies”and through the project no.457811052
文摘Considering droplet phenomena at low Mach numbers,large differences in the magnitude of the occurring characteristic waves are presented.As acoustic phenomena often play a minor role in such applications,classical explicit schemes which resolve these waves suffer from a very restrictive timestep restriction.In this work,a novel scheme based on a specific level set ghost fluid method and an implicit-explicit(IMEX)flux splitting is proposed to overcome this timestep restriction.A fully implicit narrow band around the sharp phase interface is combined with a splitting of the convective and acoustic phenomena away from the interface.In this part of the domain,the IMEX Runge-Kutta time discretization and the high order discontinuous Galerkin spectral element method are applied to achieve high accuracies in the bulk phases.It is shown that for low Mach numbers a significant gain in computational time can be achieved compared to a fully explicit method.Applica-tions to typical droplet dynamic phenomena validate the proposed method and illustrate its capabilities.
文摘The ghost fluid method and the level set function is used to simulate the Richtmyer Meshkov instability of multi material flow interface in two dimensions. To keep the density profile from smearing out, the isobaric fix technique is used in numerical algorithms.
基金This research was partially supported by the National Natural Science Foundation of China(NSFC)(10871018,10931004)the National Key Laboratory of Explosion Science and Technology,Beijing Institute of Technology(KFJJ08-7).
文摘In this work,the modified ghost fluid method is developed to deal with 2D compressible fluid interacting with elastic solid in an Euler-Lagrange coupled system.In applying the modified Ghost Fluid Method to treat the fluid-elastic solid coupling,the Navier equations for elastic solid are cast into a system similar to the Euler equations but in Lagrangian coordinates.Furthermore,to take into account the influence of material deformation and nonlinear wave interaction at the interface,an Euler-Lagrange Riemann problem is constructed and solved approximately along the normal direction of the interface to predict the interfacial status and then define the ghost fluid and ghost solid states.Numerical tests are presented to verify the resultant method.
文摘In this work,the modified Ghost Fluid Method is further developed to ap-ply to compressible fluid coupled to deformable structure,where the pressure in the structure or flow can vary from an initial extremely high magnitude(such that the solid medium can be under plastic compression)to a subsequently very low quantity(so that cavitation can occur in the fluid).New techniques are also developed in the definition of the ghost fluid status when the structure is under plastic deformation or when the flow is under cavitation next to the structure.Numerical results show that the improved MGFM for treatment of the fluid-deformable structure coupling works efficiently for all pressure ranges and is capable of simulating cavitation evolution and cavitation re-loading in conjunction with the employment of the isentropic one-fluid cavitation model.
基金supported by National Science Foundation of China (10576015)
文摘Recent years the modify ghost fluid method （MGFM） and the real ghost fluid method （RGFM） based on Riemann problem have been developed for multimedium compressible flows. According to authors, these methods have only been used with the level set technique to track the interface. In this paper, we combine the MCFM and the RGFM respectively with front tracking method, for which the fluid interfaces are explicitly tracked by connected points. The method is tested with some one-dimensional problems, and its applicability is also studied. Furthermore, in order to capture the interface more accurately, especially for strong shock impacting on interface, a shock monitor is proposed to determine the initial states of the Riemann problem. The present method is applied to various one- dimensional problems involving strong shock-interface interaction. An extension of the present method to two dimension is also introduced and preliminary results are given.
基金NSFC grant 10671091Nanjing University Talent Development Foundation and SRF for ROCS,SEM.Additional support was provided by NUS Research Project R-265-000-118-112 while he was in residence at the Department of Mechanical Engineering,National University of Singapore,Singapore 119260.
文摘The original ghost fluid method (GFM) developed in [13] and the modifiedGFM (MGFM) in [26] have provided a simple and yet flexible way to treat twomediumflow problems. The original GFM and MGFM make the material interface"invisible" during computations and the calculations are carried out as for a singlemedium such that its extension to multi-dimensions becomes fairly straightforward.The Runge-Kutta discontinuous Galerkin (RKDG) method for solving hyperbolic conservationlaws is a high order accurate finite element method employing the usefulfeatures from high resolution finite volume schemes, such as the exact or approximateRiemann solvers, TVD Runge-Kutta time discretizations, and limiters. In this paper,we investigate using RKDG finite element methods for two-medium flow simulationsin one and two dimensions in which the moving material interfaces is treated via nonconservativemethods based on the original GFM and MGFM. Numerical results forboth gas-gas and gas-water flows are provided to show the characteristic behaviors ofthese combinations.
基金the National Natural Science Foundation of China(Nos.11201442 and 10931004)。
文摘The modified ghost fluid method(MGFM)provides a robust and efficient interface treatment for various multi-medium flow simulations and some particular fluid-structure interaction(FSI)simulations.However,this methodology for one specific class of FSI problems,where the structure is plate,remains to be developed.This work is devoted to extending the MGFM to treat compressible fluid coupled with a thin elastic plate.In order to take into account the influence of simultaneous interaction at the interface,a fluid-plate coupling system is constructed at each time step and solved approximately to predict the interfacial states.Then,ghost fluid states and plate load can be defined by utilizing the obtained interfacial states.A type of acceleration strategy in the coupling process is presented to pursue higher efficiency.Several one-dimensional examples are used to highlight the utility of this method over looselycoupled method and validate the acceleration techniques.Especially,this method is applied to compute the underwater explosions(UNDEX)near thin elastic plates.Evolution of strong shock impacting on the thin elastic plate and dynamic response of the plate are investigated.Numerical results disclose that this methodology for treatment of the fluid-plate coupling indeed works conveniently and accurately for different structural flexibilities and is capable of efficiently simulating the processes of UNDEX with the employment of the acceleration strategy.
基金supported by NSFC 10531080,10972230,and 973 project 2005CB321703 and 2010CB731505supported by ARO grant W911NF-08-1-0520 and NSF grant DMS-0809086.
文摘A conservative modification to the ghost fluid method(GFM)is developed for compressible multiphase flows.The motivation is to eliminate or reduce the conservation error of the GFM without affecting its performance.We track the conservative variables near the material interface and use this information to modify the numerical solution for an interfacing cell when the interface has passed the cell.The modification procedure can be used on the GFM with any base schemes.In this paper we use the fifth order finite difference WENO scheme for the spatial discretization and the third order TVD Runge-Kutta method for the time discretization.The level set method is used to capture the interface.Numerical experiments show that the method is at least mass and momentum conservative and is in general comparable in numerical resolution with the original GFM.
基金supported under the National Natural Science Foundation of China(No.10871018)the funding of National Key Lab of Explosion Science and Technology(No.KFJJ08-7).
文摘The modified ghost fluid method(MGFM)has been shown to be robust and efficient when being applied to multi-medium compressible flows.In this paper,we rigorously analyze the optimal error estimation of the MGFM when it is applied to the multi-fluid Riemann problem.By analyzing the properties of the MGFM and the approximate Riemann problem solver(ARPS),we show that the interfacial status provided by the MGFM can achieve“third-order accuracy”in the sense of comparing to the exact solution of the Riemann problem,regardless of the solution type.In addition,our analysis further reveals that the ARPS based on a doubled shock structure in the MGFM is suitable for almost any conditions for predicting the interfacial status,and that the“natural”approach of“third-order accuracy”is practically less useful.Various examples are presented to validate the conclusions made.
文摘Two interface capturing methods are studied for multi fluid flows, governed by the stiffened gas equation of state. The mixture type interface capturing algorithm uses a simple volume fraction model Euler equations written in a quasi conservative form, which is solved by a standard high resolution piecewise parabolic method (PPM) with multi fluid Riemann solver. The level set interface capturing method uses a narrow band ghost fluid method (GFM) with no numerical smearing. Several examples are presented and compared for one and two dimensions, which show the feasibility of the two methods applied to various multi fluid problems.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.10732010,10972010,and 11332002)
文摘In the present paper, a three-dimensional （3D） Eulerian technique for the 3D numerical simulation of high-velocity impact problems is proposed. In the Eulerian framework, a complete 3D conservation element and solution element scheme for conservative hyperbolic governing equations with source terms is given. A modified ghost fluid method is proposed for the treatment of the boundary conditions. Numerical simulations of the Taylor bar problem and the ricochet phenomenon of a sphere impacting a plate target at an angle of 60~ are carried out. The numerical results are in good agreement with the corresponding experimental observations. It is proved that our computational technique is feasible for analyzing 3D high-velocity impact problems.
文摘Enhancement of two fluid mixing was numerically studied by tracking the multi fluid interfaces. Level set equations were used to capture the interfaces, and flow field was obtained by upwind TVD scheme to solve 2D Eulerian equations. The boundary conditions at interface of two fluids are determined by Ghost fluid method (GFM). The distributions of fluid parameters, such as pressure and density, were got at different time steps. The results show that the method presented in this paper can track the density discontinuity perfectly. Superior to previous results, the density discontinuity remains sharper. Also, the mixing of fluids can be greatly enhanced by setting disturbances along the initial fluid interfaces.
文摘In the paper, the numerical simulation of interface problems for multiple material fluids is studied. The level set function is designed to capture the location of the material interface. For multi-dimensional and multi-material fluids, the modified ghost fluid method needs a Riemann solution to renew the variable states near the interface. Here we present a new convenient and effective algorithm for solving the Riemann problem in the normal direction. The extrapolated variables are populated by Taylor series expansions in the direction. The anti-diffusive high order WENO difference scheme with the limiter is adopted for the numerical simulation. Finally we implement a series of numerical experiments of multi-material flows. The obtained results are satisfying, compared to those by other methods.
