激光X光转换研究中解三温方程的两种差分格式

在激光Ｘ光转换规律的研究中必须认真讨论辐射的计算，本文提出了解电子，离子，光子三温相脱离的能量方程的两种差分格式，即隐式格式的整全迭代求解和分裂格式解法，文中介绍了能量交换项与其它项分开计算的分步方法，并考察了对温度变化的影响，数值计算表明结果是正确有效的。     

电磁内爆X射线数值计算

 叙述一维辐射磁流体力学方程组（三温）的数值模拟,用于电磁内爆X射线计算。动量方程、三温方程组采用分裂法求解,以克服强耦合非线性可能引起的数值不稳定性。采用隐式差分格式避免小时间步长限制。隐式格式采用Newton-Raphson法迭代求解。使用了分层网格、加权平均,迭代初值公式,自动调整时间步长,奇点处理等特殊处理方法。最后给出Pegasus l装置一个实验模型的计算结果,获得15 J的X射线能量,峰值功率为16 GW。     

二维三温能量方程的九点差分格式及其迭代解法

给出了激光驱动内爆数值模拟中二维三温能量方程的九点差分格式及其迭代解法,并给出了九点差分格式与五点差分格式对比计算结果,对一维和二维物理模型进行了数值模拟,得到了令人满意的数值结果。     

三温方程中的电子-辐射能量交换项的计算

讨论了激光靶耦合数值模拟中电子-辐射能量交换项的处理。为了克服它给三温方程求解带来的困难,在非平衡区(T_e>T_r)和接近平衡区(T_e~T_r)采用了不同的能量交换项形式。此外,由于接近平衡时的能量交换项用辐射平均自由程表示,避免了此时计算束缚电子占据概率不准确的问题。     

辐射扩散计算方法若干研究进展

辐射流体力学研究辐射的传输对流体运动的影响,并在此条件下研究流体的运动规律.实际应用问题中辐射流体力学所描述的是非常复杂的物理过程,数值模拟是主要的研究手段之一.模拟通常采用流体计算和辐射计算分裂求解的方法.讨论求解辐射扩散方程时迫切需要解决的一些计算方法问题,包括大变形网格上扩散计算格式与非线性迭代方法,并简要介绍部分研究进展.     

数值模拟辐射驱动2维内爆压缩过程

采用2维三温非平衡辐射流体力学程序LARED-H数值模拟辐射驱动内爆过程。针对2维三温能量方程九点差分格式离散后的线性方程组,采用了高效的Krysolv子空间迭代解法,改进了代数解法器。将1维间接驱动内爆总体程序CFJ与LARED-H程序的计算结果进行比对,验证了LARED-H程序数值模拟1维内爆问题的正确性。并数值模拟了不同腔长辐射温度源驱动下的2维靶球运动,数值结果显示:随着腔长的增加,高压缩内爆燃料区分别被压缩成香肠形、球形和铁饼形,数值模拟结果与神光Ⅱ的实验结果定性上相同。     

�����ά���·���������ѧ�������������Ǩ�Ʒ�

在三维三温辐射流体力学方程组求解过程中,为解决复杂的等离子体流场和难度大的三维网格计算,提出了用物理量迁移法求解方程组,并对三温能量方程采用了一种改进的显示算法。由此获得了较好的数值结果。     

辐射扩散方程的自适应坐标变换方法

借鉴自适应坐标变换流体力学方法和其他一些有益工作，研究辐射扩散方程的自适应坐标变换方法，并考虑辐射扩散和流体力学耦合的辐射流体力学问题。对辐射扩散方程和一般二维对流扩散方程，从积分和微分形式出发，导出自适应坐标变换下的等价形式，对数值离散格式进行初步讨论。对二维三温辐射扩散方程，采用时空有限体积方法，进行数值离散；利用局部函数近似离散扩散算子。利用合理的进程分裂，简化耦合问题的计算。     

不同辐射建模下球腔的整体数值模拟

不同辐射建模对于腔内辐射场描述的精确程度不同,需要分析不同建模对腔内辐射温度的影响。开展了三温建模与辐射多群输运建模下LARED集成程序数值模拟两孔球型黑腔模型,实现了球腔的完整数值模拟。数值模拟结果表明,三温与辐射多群输运模拟的等离子体状态接近,辐射温度存在差异。物理分析显示辐射温度差异的主要原因是使用的辐射不透明度,修改辐射不透明度参数后的三温计算结果与输运计算符合更好,从而可以用三温建模更快更准确地估计出所需的激光能量和功率。     

金属氧化物半导体场效应管长期辐射效应的数值模拟

将粒子输运的蒙特卡罗方法与器件数值模拟的有限体积法相耦合来模拟典型金属氧化物半导体场效应管(MOSFET)的长期辐射效应。二氧化硅中的陷阱电荷及硅中的自由电子和空穴均使用漂移扩散模型来描述,入射粒子的能量沉积可作为源项耦合至漂移扩散模型方程,并根据有限体积法得到控制方程的离散格式,方程的数值解即为MOSFET的长期辐射响应结果。使用该方法模拟了MOSFET受射线粒子辐照后的阈值电压漂移与关态漏电流现象。结果表明,耦合方法适用于典型半导体器件长期辐射效应模拟,其阈值电压漂移及漏电流计算结果与文献符合较好。     

A high order kinetic flux-vector splitting method for the reduced five-equation model of compressible two-fluid flows

We present a high order kinetic flux-vector splitting (KFVS) scheme for the numerical solution of a conservative interface-capturing five-equation model of compressible two-fluid flows. This model was initially introduced by Wackers and Koren (2004) [21]. The flow equations are the bulk equations, combined with mass and energy equations for one of the two fluids. The latter equation contains a source term in order to account for the energy exchange. We numerically investigate both one- and two-dimensional flow models. The proposed numerical scheme is based on the direct splitting of macroscopic flux functions of the system of equations. In two space dimensions the scheme is derived in a usual dimensionally split manner. The second order accuracy of the scheme is achieved by using MUSCL-type initial reconstruction and Runge–Kutta time stepping method. For validation, the results of our scheme are compared with those from the high resolution central scheme of Nessyahu and Tadmor [14]. The accuracy, efficiency and simplicity of the KFVS scheme demonstrate its potential for modeling two-phase flows.     

Algebraic Splitting for Incompressible Navier–Stokes Equations

Fully discretized incompressible Navier–Stokes equations are solved by splitting the algebraic system with an approximate factorization. This splitting affects the temporal convergence order of velocity and pressure. The splitting error is proportional to the pressure variable, and a simple analysis shows that the original convergence order of the time-integration scheme can be retained by solving for incremental pressure. The combination of splitting and incremental pressure is shown to be equivalent to an error-correcting method using the full pressure. In numerical experiments employing a third-order time-integration scheme and various orders for the pressure increment, the splitting error is shown to control the convergence order, and the full order of the scheme is recaptured for both velocity and pressure. The difference between perturbing the momentum or the continuity equation is also explored.     

A gradient stable scheme for a phase field model for the moving contact line problem

In this paper, an efficient numerical scheme is designed for a phase field model for the moving contact line problem, which consists of a coupled system of the Cahn–Hilliard and Navier–Stokes equations with the generalized Navier boundary condition [1], [2], [4]. The nonlinear version of the scheme is semi-implicit in time and is based on a convex splitting of the Cahn–Hilliard free energy (including the boundary energy) together with a projection method for the Navier–Stokes equations. We show, under certain conditions, the scheme has the total energy decaying property and is unconditionally stable. The linearized scheme is easy to implement and introduces only mild CFL time constraint. Numerical tests are carried out to verify the accuracy and stability of the scheme. The behavior of the solution near the contact line is examined. It is verified that, when the interface intersects with the boundary, the consistent splitting scheme [21], [22] for the Navier Stokes equations has the better accuracy for pressure.     

流场/电磁场的数值优化计算

采用矢通量分裂方法计算了绕翼型的时域电磁散射场特性及雷达散射截面积(RCS)。采用VanLeer矢通量分裂格式求解Euler方程方法计算了绕翼型的气动特性。并用一种简单而有效的数值优化方法对流场解和电磁场解进行了综合优化计算     

轴对称可压缩流的统一坐标系

将Hui等人提出的统一坐标系推广到轴对称欧拉方程,讨论统一坐标系下的轴对称欧拉方程的解法,并通过数值算例证明轴对称统一坐标系的优越性.     

A consistent-splitting approach to computing stiff steady-state reacting flows with adaptive chemistry

Splitting techniques have been used extensively for computing reacting flows with detailed chemistry. Nevertheless, there are still some open questions with respect to efficiency and the error introduced by splitting. In this paper, the accuracy and effectiveness of split-operator methods for computing steady-state reacting flows are determined. A fully coupled scheme is described together with two splitting schemes: a standard Strang-splitting scheme and a consistent-splitting scheme, all with implicit transport computations. The effect of splitting errors on the convergence and solution accuracy is investigated analytically using a one-dimensional scalar equation. The accuracy with respect to the original discretized equations is tested for an H₂/O₂ burner flame. Finally, consistent splitting is combined with an adaptive chemistry approach to compute three partially premixed laminar methane flames using detailed chemistry (217 reactions). The calculations confirm that HCO radical concentration is an excellent surrogate for heat release rate.     

用代数多重网格法求解一维分裂格式的Euler方程

提出了代数多重网格法(AMG)的一种新算法。新算法改进了插值公式和粗网格方程,并把它应用到求解一维的分裂格式Euler方程。数值结果表明,对于具有高CFL条件数的Euler方程,代数多重网格法可以求解;对于Gaus-Seidel方法求解不能收敛的代数方程组,代数多重网格法求解可以收敛。新算法改进了代数多重网格法的收敛性和扩展了它的应用范围,数值结果表明了它的有效性和强壮性。     

Multi-Symplectic Wavelet Collocation Method for Maxwell's Equations

In this paper, we develop a multi-symplectic wavelet collocation method for three-dimensional (3-D) Maxwell's equations. For the multi-symplectic formulation of the equations, wavelet collocation method based on autocorrelation functions is applied for spatial discretization and appropriate symplectic scheme is employed for time integration. Theoretical analysis shows that the proposed method is multi-symplectic, unconditionally stable and energy-preserving under periodic boundary conditions. The numerical dispersion relation is investigated. Combined with splitting scheme, an explicit splitting symplectic wavelet collocation method is also constructed. Numerical experiments illustrate that the proposed methods are efficient, have high spatial accuracy and can preserve energy conservation laws exactly.     

A LAX-WENDROFF TYPE SCHEME WITH PREDICTOR-CORRECTOR

In this paper a Lax-Wendroff type scheme with predictor-corrector is presented for solving the hyperbolic equations of fluid flow. The predictor solutions without introducing nonphysical oscillation are obtained by the method of flux vector splitting. The numerical flux of the corrector step is an anti-diffusive term, on which contraints are imposed to decrease the oscillation. Numerical experiments are presented to illustrate the performance of our proposed scheme.     

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

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