The least-square gridless method was extended to simulate the compressible multi-material flows. The algorithm was accomplished to solve the Arbitrary Lagrange-Euler( ALE) formulation. The local least-square curve fit...The least-square gridless method was extended to simulate the compressible multi-material flows. The algorithm was accomplished to solve the Arbitrary Lagrange-Euler( ALE) formulation. The local least-square curve fits was adopted to approximate the spatial derivatives of a point on the base of the points in its circular support domain,and the basis function was linear. The HLLC( Harten-Lax-van Leer-Contact) scheme was used to calculate the inviscid flux. On the material interfaces,the gridless points were endued with a dual definition corresponding to different materials. The moving velocity of the interface points was updated by solving the Riemann problem. The interface boundary condition was built by using the Ghost Fluid Method( GFM).Computations were performed for several one and two dimensional typical examples. The numerical results show that the interface and the shock wave are well captured,which proves the effectiveness of gridless method in dealing with multi-material flow problems.展开更多
Hydrocodes are necessary numerical tools in the fields of implosion and high-velocity impact,which often involve large deformations with changing-topology interfaces.It is very difficult for Lagrangian or Simplified A...Hydrocodes are necessary numerical tools in the fields of implosion and high-velocity impact,which often involve large deformations with changing-topology interfaces.It is very difficult for Lagrangian or Simplified Arbitrary Lagrangian-Eulerian(SALE)codes to tackle these kinds of large-deformation problems,so a stag-gered Multi-Material ALE(MMALE)code is developed in this paper,which is the explicit time-marching Lagrange plus remap type.We use the Moment Of Fluid(MOF)method to reconstruct the interfaces of multi-material cells and present an adaptive bisection method to search for the global minimum value of the nonlinear objective function.To keep the Lagrangian computations as long as possible,we develop a ro-bust rezoning method named as Combined Rezoning Method(CRM)to generate the convex,smooth grids for the large-deformation domain.Regarding the staggered remap phase,we use two methods to remap the variables of Lagrangian mesh to the rezoned one.One is the first-order intersection-based remapping method that doesn’t limit the distances between the rezoned and Lagrangian meshes,so it can be used in the applications of wide scope.The other one is the conservative second-order flux-based remapping method developed by Kucharika and Shashkov[22]that requires the rezoned element to locate in its adjacent old elements.Numerical results of triple point problem show that the result of first-order remapping method using ALE computations is gradually convergent to that of second-order remapping method using Eulerian computations with the decrease of rezoning,thereby telling us that MMALE computations should be performed as few as possible to reduce the errors of the interface reconstruction and the remapping.Numerical results provide a clear evidence of the robustness and the accuracy of this MMALE scheme,and that our MMALE code is powerful for the large-deformation problems.展开更多
文摘The least-square gridless method was extended to simulate the compressible multi-material flows. The algorithm was accomplished to solve the Arbitrary Lagrange-Euler( ALE) formulation. The local least-square curve fits was adopted to approximate the spatial derivatives of a point on the base of the points in its circular support domain,and the basis function was linear. The HLLC( Harten-Lax-van Leer-Contact) scheme was used to calculate the inviscid flux. On the material interfaces,the gridless points were endued with a dual definition corresponding to different materials. The moving velocity of the interface points was updated by solving the Riemann problem. The interface boundary condition was built by using the Ghost Fluid Method( GFM).Computations were performed for several one and two dimensional typical examples. The numerical results show that the interface and the shock wave are well captured,which proves the effectiveness of gridless method in dealing with multi-material flow problems.
基金This work was performed under the auspices of National Natural Science Founda-tion of China(No.11702030)NSAF(No.U1630247).
文摘Hydrocodes are necessary numerical tools in the fields of implosion and high-velocity impact,which often involve large deformations with changing-topology interfaces.It is very difficult for Lagrangian or Simplified Arbitrary Lagrangian-Eulerian(SALE)codes to tackle these kinds of large-deformation problems,so a stag-gered Multi-Material ALE(MMALE)code is developed in this paper,which is the explicit time-marching Lagrange plus remap type.We use the Moment Of Fluid(MOF)method to reconstruct the interfaces of multi-material cells and present an adaptive bisection method to search for the global minimum value of the nonlinear objective function.To keep the Lagrangian computations as long as possible,we develop a ro-bust rezoning method named as Combined Rezoning Method(CRM)to generate the convex,smooth grids for the large-deformation domain.Regarding the staggered remap phase,we use two methods to remap the variables of Lagrangian mesh to the rezoned one.One is the first-order intersection-based remapping method that doesn’t limit the distances between the rezoned and Lagrangian meshes,so it can be used in the applications of wide scope.The other one is the conservative second-order flux-based remapping method developed by Kucharika and Shashkov[22]that requires the rezoned element to locate in its adjacent old elements.Numerical results of triple point problem show that the result of first-order remapping method using ALE computations is gradually convergent to that of second-order remapping method using Eulerian computations with the decrease of rezoning,thereby telling us that MMALE computations should be performed as few as possible to reduce the errors of the interface reconstruction and the remapping.Numerical results provide a clear evidence of the robustness and the accuracy of this MMALE scheme,and that our MMALE code is powerful for the large-deformation problems.
