A second order self-consistent IMEX method for radiation hydrodynamics

We present a second order self-consistent implicit/explicit (methods that use the combination of implicit and explicit discretizations are often referred to as IMEX (implicit/explicit) methods ,  and ) time integration technique for solving radiation hydrodynamics problems. The operators of the radiation hydrodynamics are splitted as such that the hydrodynamics equations are solved explicitly making use of the capability of well-understood explicit schemes. On the other hand, the radiation diffusion part is solved implicitly. The idea of the self-consistent IMEX method is to hybridize the implicit and explicit time discretizations in a nonlinearly consistent way to achieve second order time convergent calculations. In our self-consistent IMEX method, we solve the hydrodynamics equations inside the implicit block as part of the nonlinear function evaluation making use of the Jacobian-free Newton Krylov (JFNK) method ,  and . This is done to avoid order reductions in time convergence due to the operator splitting. We present results from several test calculations in order to validate the numerical order of our scheme. For each test, we have established second order time convergence.     

一种新型光滑粒子动力学固壁边界施加模型

由于Lagrange粒子法的本质, 固壁边界条件的施加一直是光滑粒子动力学方法的难点之一. 本文从固壁边界的物理原理出发, 应用多层虚粒子表征固壁边界, 提出了一种新型固壁边界施加模型. 将虚粒子看作流体的扩展, 计算中虚粒子密度保持不变, 压力、速度等参数通过对流体粒子的插值获得, 虚粒子有条件的参与控制方程的计算, 对流体的密度/压力产生影响, 通过压力梯度隐式地表征壁面与流体之间的作用强度并对流体粒子施加沿壁面法线方向的斥力作用, 防止流体粒子对壁面的穿透. 数值算例测试结果表明, 与现有固壁边界施加方法相比, 本文方法更加符合流体与固壁边界作用的物理原理, 可以简单、有效地施加固壁边界条件, 方便地应用于具有复杂几何边界的问题, 获得稳定的流场形态、规则的粒子秩序及良好的速度、压力等参量的分布.     

适用于等熵流动的交错拉氏Godunov方法

为消除传统单元中心型Godunov方法在求解稀疏波问题时的非物理过热现象,发展一种适用于等熵流动的交错拉氏Godunov方法.主要的特征是采用速度与热力学变量交错分布的形式,避免在单元内进行速度平均,从而消除由于动量平均过程导致的动能耗散.与传统的von Neumann型交错网格方法相比,网格的边界通量由节点处的多维黎曼求解器提供,克服了多维人工粘性选取带来的困难.为减少多维黎曼求解器在求解稀疏波问题时的非物理熵增,给出稀疏波出现的合理判据,从而保证了热力学关系式的满足.数值实验表明：该方法能很好地消除稀疏波的过热现象,同时在求解激波问题时又能保持与传统单元中心型拉氏方法相同的激波捕捉能力.     

Fully implicit 1D radiation hydrodynamics: Validation and verification

A fully implicit finite difference scheme has been developed to solve the hydrodynamic equations coupled with radiation transport. Solution of the time-dependent radiation transport equation is obtained using the discrete ordinates method and the energy flow into the Lagrangian meshes as a result of radiation interaction is fully accounted for. A tridiagonal matrix system is solved at each time step to determine the hydrodynamic variables implicitly. The results obtained from this fully implicit radiation hydrodynamics code in the planar geometry agrees well with the scaling law for radiation driven strong shock propagation in aluminium. For the point explosion problem the self similar solutions are compared with results for pure hydrodynamic case in spherical geometry. Results obtained when radiation interaction is also accounted agree with those of point explosion with heat conduction for lower input energies. Having, thus, benchmarked the code, self convergence of the method w.r.t. time step is studied in detail for both the planar and spherical problems. Spatial as well as temporal convergence rates are ?1 as expected from the difference forms of mass, momentum and energy conservation equations. This shows that the asymptotic convergence rate of the code is realized properly.     

A direct Eulerian GRP scheme for relativistic hydrodynamics: One-dimensional case

The paper proposes a direct Eulerian generalized Riemann problem (GRP) scheme for one-dimensional relativistic hydrodynamics. It is an extension of the Eulerian GRP scheme for compressible non-relativistic hydrodynamics proposed in [M. Ben-Artzi, J.Q. Li, G. Warnecke, A direct Eulerian GRP scheme for compressible fluid flows, J. Comput. Phys. 218 (2006) 19–43]. Two main ingredients, the Riemann invariant and the Rankine–Hugoniot jump condition, are directly used to resolve the local GRP in the Eulerian formulation, and thus the crucial and delicate Lagrangian treatment in the original GRP scheme [3] can be avoided. Several numerical examples are given to demonstrate the accuracy and effectiveness of the proposed GRP scheme.     

一种改进的光滑粒子流体动力学前处理方法

为了便于对任意边界形状的计算域快速地布置均匀粒子,提出了一种改进的光滑粒子流体动力学前处理方法.该方法是在2012年Colagrossi等提出的算法基础上进行改进后得到的.Colagrossi等提出的算法能够计算一些简单外形分布比较均匀的粒子.然而当光滑长度与初始粒子间距的比值较大时该方法在计算过程中会出现较强的数值震荡问题,收敛速度慢;而且在计算过程中可能会遭遇流体粒子穿透固体壁面的问题.本文通过引入未知因素修正的平滑粒子动力学模型来提高计算稳定性,并通过对边界附近的流体粒子施加边界力来避免流体粒子穿透固体壁面.算例验证结果表明,利用改进后的光滑粒子流体动力学前处理方法能够快速地对各种边界形状的计算域分布均匀粒子,并且避免了流体粒子穿透固体壁面的问题.     

光滑粒子流体动力学方法固壁处理的一种新型排斥力模型

与传统网格法相比, 光滑粒子流体动力学方法不能直接施加壁面边界条件, 这就限制了该方法在工程中的应用.为此, 本文基于Galerkin加权余量法并结合传统排斥力方法, 推导出一种新的排斥力公式来施加壁面边界条件.该方法不含未知参数, 能在不减小边界粒子尺寸的情形下有效地防止流体粒子穿透壁面, 同时可避免邻近边界的流体粒子的速度及压力振荡. 分别通过静止液柱算例、液柱坍塌算例、容器中液体静止算例及溃坝算 例来验证本文方法的有效性, 并与传统边界处理方法进行对比, 结果表明: 本文方法克服了传统方法存在的缺陷, 是一种有效的固壁边界处理方法. 关键词： 光滑粒子流体动力学法 固壁边界 排斥力 加权余量法     

Non-analog Monte Carlo estimators for radiation momentum deposition

The standard method for calculating radiation momentum deposition in Monte Carlo simulations is the analog estimator, which tallies the change in a particle's momentum at each interaction with the matter. Unfortunately, the analog estimator can suffer from large amounts of statistical error. In this paper, we present three new non-analog techniques for estimating momentum deposition. Specifically, we use absorption, collision, and track-length estimators to evaluate a simple integral expression for momentum deposition that does not contain terms that can cause large amounts of statistical error in the analog scheme. We compare our new non-analog estimators to the analog estimator with a set of test problems that encompass a wide range of material properties and both isotropic and anisotropic scattering. In nearly all cases, the new non-analog estimators outperform the analog estimator. The track-length estimator consistently yields the highest performance gains, improving upon the analog-estimator figure of merit by factors of up to two orders of magnitude.     

高速CT阻滞产生X射线辐射的数值模拟研究

以Ar等离子体为例,应用一维辐射流体力学模型,对高速度、高密度的紧凑等离子体环(CT)撞击静态靶减速产生软X射线辐射的物理过程进行了数值模拟研究。利用分裂格式方法计算得到了描述碰撞物理过程的各力学量随时间的演化关系,以及固定CT动能时使停滞CT成为最有效辐射源的初始速度和质量的最优搭配。对计算过程中相关物理量的跟踪分析表明,在1%的误差范围内计算格式满足能量守恒方程。     

A cell-centered Lagrangian hydrodynamics scheme on general unstructured meshes in arbitrary dimension

We describe a cell-centered Godunov scheme for Lagrangian gas dynamics on general unstructured meshes in arbitrary dimension. The construction of the scheme is based upon the definition of some geometric vectors which are defined on a moving mesh. The finite volume solver is node based and compatible with the mesh displacement. We also discuss boundary conditions. Numerical results on basic 3D tests problems show the efficiency of this approach. We also consider a quasi-incompressible test problem for which our nodal solver gives very good results if compared with other Godunov solvers. We briefly discuss the compatibility with ALE and/or AMR techniques at the end of this work. We detail the coefficients of the isoparametric element in the appendix.     

A numerical model for multigroup radiation hydrodynamics

We present in this paper a multigroup model for radiation hydrodynamics to account for variations of the gas opacity as a function of frequency. The entropy closure model (M₁) is applied to multigroup radiation transfer in a radiation hydrodynamics code. In difference from the previous grey model, we are able to reproduce the crucial effects of frequency-variable gas opacities, a situation omnipresent in physics and astrophysics. We also account for the energy exchange between neighbouring groups which is important in flows with strong velocity divergence. These terms were computed using a finite volume method in the frequency domain. The radiative transfer aspect of the method was first tested separately for global consistency (reversion to grey model) and against a well-established kinetic model through Marshak wave tests with frequency-dependent opacities. Very good agreement between the multigroup M₁ and kinetic models was observed in all tests. The successful coupling of the multigroup radiative transfer to the hydrodynamics was then confirmed through a second series of tests. Finally, the model was linked to a database of opacities for a Xe gas in order to simulate realistic multigroup radiative shocks in Xe. The differences with the previous grey models are discussed.     

Simulation of surface tension in 2D and 3D with smoothed particle hydrodynamics method

The methods for simulating surface tension with smoothed particle hydrodynamics (SPH) method in two dimensions and three dimensions are developed. In 2D surface tension model, the SPH particle on the boundary in 2D is detected dynamically according to the algorithm developed by Dilts [G.A. Dilts, Moving least-squares particle hydrodynamics II: conservation and boundaries, International Journal for Numerical Methods in Engineering 48 (2000) 1503–1524]. The boundary curve in 2D is reconstructed locally with Lagrangian interpolation polynomial. In 3D surface tension model, the SPH particle on the boundary in 3D is detected dynamically according to the algorithm developed by Haque and Dilts [A. Haque, G.A. Dilts, Three-dimensional boundary detection for particle methods, Journal of Computational Physics 226 (2007) 1710–1730]. The boundary surface in 3D is reconstructed locally with moving least squares (MLS) method. By transforming the coordinate system, it is guaranteed that the interface function is one-valued in the local coordinate system. The normal vector and curvature of the boundary surface are calculated according to the reconstructed boundary surface and then surface tension force can be calculated. Surface tension force acts only on the boundary particle. Density correction is applied to the boundary particle in order to remove the boundary inconsistency. The surface tension models in 2D and 3D have been applied to benchmark tests for surface tension. The ability of the current method applying to the simulation of surface tension in 2D and 3D is proved.     

光滑粒子动力学方法中粒子分布与数值稳定性分析

光滑粒子动力学(SPH)作为一种拉格朗日型无网格粒子方法,已经成功地应用于包括含多相流动界面以及移动边界的可压缩和不可压缩流体运动的研究中.通过对Poiseuille流动的深入研究,探索了SPH方法中粒子分布对计算精度的影响,揭示了一种因为粒子不规则分布而导致的数值不稳定现象.研究显示,这种数值不稳定性起源于SPH方法粒子近似过程中的不连续性.使用了一种新的粒子近似格式以确保SPH方法中粒子近似的连续性.计算结果表明,这种新的粒子近似格式对于规则和不规则的粒子分布都能得到稳定精度的结果.     

渐近保持PIC算法在等离子体粒子模拟中的应用

针对准中性Vlasov-Poisson系统中多尺度物理参数共存的问题开展研究,采用渐近保持格式推导出与Vlasov方程相耦合的Poisson方程。数值模拟了一维等离子体在真空中的膨胀,比较了渐近保持PIC算法与传统PIC算法模拟结果的异同。结果表明：两种算法在满足多尺度中小尺度量限制的前提下,其结果保持高度一致;且在克服小尺度量束缚时,渐近保持PIC算法依然可以正确地模拟等离子体的物理行为,大大节省了计算时间。     

充模过程中熔接痕的改进光滑粒子动力学方法模拟与预测

提出了一种黏弹性流体的改进光滑粒子动力学(SPH)方法以试探性地模拟和预测黏弹性FENE-P熔体充模过程中熔接痕的形态演化. 首先基于SPH方法建立了聚合物流动的宏微观耦合模型, 同时提出了黏弹性流体的改进SPH离散格式. 随后, 通过模拟一些基准算例验证了改进的SPH方法模拟聚合物宏微观耦合问题的有效性及收敛性, 以及所提出的黏弹性温度模型的有效性. 最后, 模拟了环型腔内的充模过程, 试探性地展示了充模过程中微观分子的变形过程. 同时采用顺序热流道技术模拟了多浇口C形腔内的充模过程, 并与其他数值结果做比较. 数值结果表明: 对于大制件多浇口充模过程, 顺序热流道技术能够改善甚至消除充模过程中的熔接痕.     

A study of the radiation of a Nd:YAG laser-MIG hybrid plasma

To study the coupling of a Nd:YAG laser beam and a MIG arc with the spectrum of the plasma, a new hollow probe method is used to collect radiation from specific points within the arc plasma. With the probe, the spatial distributions of the MIG arc radiation with and without the Nd:YAG laser hybrid plasma are collected and analyzed by a fiber spectrometer. Through analysis of the spectrum, radiation in the ultraviolet (Fe II dominated), visible (Fe I dominated) and infra (Ar I dominated) spectral zones are acquired for comparison of the difference between the MIG and Laser-MIG hybrid welding processes. The electronic temperature is calculated using the Boltzmann plot method and the electronic density is calculated using the Stark broadening method. The results show that the Laser-MIG hybrid welding processes cause the plasma energy to focus on the center of the welding arc and approach the welding pool. An ionizing duct abundant with Fe ions will be formed near the position of the impinging laser beam, which makes the welding arc stable. The electronic temperature of the Laser-MIG hybrid welding arc is a little higher than that of the MIG. The electronic density of Laser-MIG hybrid welding is higher than that of the MIG.     

基于光滑粒子动力学方法的液滴冲击固壁面问题数值模拟

采用改进的光滑粒子动力学(SPH)方法对液滴冲击固壁面问题进行了数值模拟. 为了提高传统SPH方法的计算精度和数值稳定性, 在传统的SPH方法的基础上对粒子方法中的密度和核梯度进行了修正, 采用了考虑黎曼解法的SPH流体控制方程, 构造了一种新型的粒子间相互作用力(IIF)模型来模拟表面张力的影响. 应用改进的SPH方法对液滴冲击固壁面问题进行了数值模拟. 计算结果表明：新型的IIF 模型能够较好地模拟表面张力的影响, 改进的SPH方法能够精细地描述液滴与固壁面相互作用过程中液滴的内部压力场演变和自由面形态变化, 液滴的铺展因子随初始韦伯数的增大而增大, 数值模拟结果与实验得到的结果基本一致. 关键词： 液滴 固壁面 光滑粒子动力学 表面张力     

液滴溅落问题的光滑粒子动力学模拟

对传统的光滑粒子动力学方法进行了改进, 改进的光滑粒子动力学方法对传统粒子方法中的核近似式和粒子近似式进行了修正, 采用Riemann 算法求解光滑粒子动力学流体控制方程, 添加了表面张力的计算程序, 考虑了表面张力对液滴溅落的影响. 应用改进的光滑粒子动力学方法对液滴静止状态下冲击液面的飞溅过程进行了数值模拟. 计算结果表明, 改进的光滑粒子动力学方法能够有效地描述液滴溅落液面的动力学特性和自由表面变化特征, 能够得到稳定精度的结果.     

A Multiphase Godunov Method for Compressible Multifluid and Multiphase Flows

We propose a new model and a solution method for two-phase compressible flows. The model involves six equations obtained from conservation principles applied to each phase, completed by a seventh equation for the evolution of the volume fraction. This equation is necessary to close the overall system. The model is valid for fluid mixtures, as well as for pure fluids. The system of partial differential equations is hyperbolic. Hyperbolicity is obtained because each phase is considered to be compressible. Two difficulties arise for the solution: one of the equations is written in non-conservative form; non-conservative terms exist in the momentum and energy equations. We propose robust and accurate discretisation of these terms. The method solves the same system at each mesh point with the same algorithm. It allows the simulation of interface problems between pure fluids as well as multiphase mixtures. Several test cases where fluids have compressible behavior are shown as well as some other test problems where one of the phases is incompressible. The method provides reliable results, is able to compute strong shock waves, and deals with complex equations of state.