The direct simulation Monte Carlo(DSMC) method is the most mature and wildly used approach for nonequilibrium gas flow simulation.The phenomenological nature of this method brings flexibility to the computation algori...The direct simulation Monte Carlo(DSMC) method is the most mature and wildly used approach for nonequilibrium gas flow simulation.The phenomenological nature of this method brings flexibility to the computation algorithms.In this study,the theoretical foundations to decouple the molecular motion and collision within a time step are discussed in detail,which can be treated as criterions for the DSMC algorithms.Based on the theoretical developments,an improved local time stepping scheme is proposed,which specifies the movement time attribute and the collision time attribute for each representative particle.A free flow about a sphere body is considered as an example,which is compared with the calculations using the published local time stepping technique.The results show that the improved local time scheme is valid and is promising in realizing flow structures with strong variations.展开更多
In this paper,a new sharp-interface approach to simulate compressible multiphase flows is proposed.The new scheme consists of a high order WENO finite volume scheme for solving the Euler equations coupled with a high ...In this paper,a new sharp-interface approach to simulate compressible multiphase flows is proposed.The new scheme consists of a high order WENO finite volume scheme for solving the Euler equations coupled with a high order pathconservative discontinuous Galerkin finite element scheme to evolve an indicator function that tracks the material interface.At the interface our method applies ghost cells to compute the numerical flux,as the ghost fluid method.However,unlike the original ghost fluid scheme of Fedkiw et al.[15],the state of the ghost fluid is derived from an approximate-state Riemann solver,similar to the approach proposed in[25],but based on a much simpler formulation.Our formulation leads only to one single scalar nonlinear algebraic equation that has to be solved at the interface,instead of the system used in[25].Away from the interface,we use the new general Osher-type flux recently proposed by Dumbser and Toro[13],which is a simple but complete Riemann solver,applicable to general hyperbolic conservation laws.The time integration is performed using a fully-discrete one-step scheme,based on the approaches recently proposed in[5,7].This allows us to evolve the system also with time-accurate local time stepping.Due to the sub-cell resolution and the subsequent more restrictive time-step constraint of the DG scheme,a local evolution for the indicator function is applied,which is matched with the finite volume scheme for the solution of the Euler equations that runs with a larger time step.The use of a locally optimal time step avoids the introduction of excessive numerical diffusion in the finite volume scheme.Two different fluids have been used,namely an ideal gas and a weakly compressible fluid modeled by the Tait equation.Several tests have been computed to assess the accuracy and the performance of the new high order scheme.A verification of our algorithm has been carefully carried out using exact solutions as well as a comparison with other numerical reference solutions.The material interface is resolved sharply and accurately without spurious oscillations in the pressure field.展开更多
不连续伽辽金时域有限元法(discontinuous Galerkin-finite element time domain,DG-FETD)便于处理多尺度电磁问题,但是由于精细结构或高介电参数媒质的存在,考虑到稳定性条件,整体时间步长的选取会受到最小剖分尺寸的影响,导致计算效...不连续伽辽金时域有限元法(discontinuous Galerkin-finite element time domain,DG-FETD)便于处理多尺度电磁问题,但是由于精细结构或高介电参数媒质的存在,考虑到稳定性条件,整体时间步长的选取会受到最小剖分尺寸的影响,导致计算效率降低。针对该问题,在传统局部时间步长(local time stepping,LTS)技术的基础上提出一种基于蛙跳格式的任意比例LTS技术,该方法在求解多尺度问题时,减少了迭代所需时间,提高了不同时间步长选取的灵活性,同时随着未知量增加,其优势更加明显,结合数值算例,验证了该方法的正确性和有效性。展开更多
文摘The direct simulation Monte Carlo(DSMC) method is the most mature and wildly used approach for nonequilibrium gas flow simulation.The phenomenological nature of this method brings flexibility to the computation algorithms.In this study,the theoretical foundations to decouple the molecular motion and collision within a time step are discussed in detail,which can be treated as criterions for the DSMC algorithms.Based on the theoretical developments,an improved local time stepping scheme is proposed,which specifies the movement time attribute and the collision time attribute for each representative particle.A free flow about a sphere body is considered as an example,which is compared with the calculations using the published local time stepping technique.The results show that the improved local time scheme is valid and is promising in realizing flow structures with strong variations.
文摘In this paper,a new sharp-interface approach to simulate compressible multiphase flows is proposed.The new scheme consists of a high order WENO finite volume scheme for solving the Euler equations coupled with a high order pathconservative discontinuous Galerkin finite element scheme to evolve an indicator function that tracks the material interface.At the interface our method applies ghost cells to compute the numerical flux,as the ghost fluid method.However,unlike the original ghost fluid scheme of Fedkiw et al.[15],the state of the ghost fluid is derived from an approximate-state Riemann solver,similar to the approach proposed in[25],but based on a much simpler formulation.Our formulation leads only to one single scalar nonlinear algebraic equation that has to be solved at the interface,instead of the system used in[25].Away from the interface,we use the new general Osher-type flux recently proposed by Dumbser and Toro[13],which is a simple but complete Riemann solver,applicable to general hyperbolic conservation laws.The time integration is performed using a fully-discrete one-step scheme,based on the approaches recently proposed in[5,7].This allows us to evolve the system also with time-accurate local time stepping.Due to the sub-cell resolution and the subsequent more restrictive time-step constraint of the DG scheme,a local evolution for the indicator function is applied,which is matched with the finite volume scheme for the solution of the Euler equations that runs with a larger time step.The use of a locally optimal time step avoids the introduction of excessive numerical diffusion in the finite volume scheme.Two different fluids have been used,namely an ideal gas and a weakly compressible fluid modeled by the Tait equation.Several tests have been computed to assess the accuracy and the performance of the new high order scheme.A verification of our algorithm has been carefully carried out using exact solutions as well as a comparison with other numerical reference solutions.The material interface is resolved sharply and accurately without spurious oscillations in the pressure field.
文摘不连续伽辽金时域有限元法(discontinuous Galerkin-finite element time domain,DG-FETD)便于处理多尺度电磁问题,但是由于精细结构或高介电参数媒质的存在,考虑到稳定性条件,整体时间步长的选取会受到最小剖分尺寸的影响,导致计算效率降低。针对该问题,在传统局部时间步长(local time stepping,LTS)技术的基础上提出一种基于蛙跳格式的任意比例LTS技术,该方法在求解多尺度问题时,减少了迭代所需时间,提高了不同时间步长选取的灵活性,同时随着未知量增加,其优势更加明显,结合数值算例,验证了该方法的正确性和有效性。
