Euler-Lagrange耦合计算的GEL方法

研究了一种Euler-Lagrange耦合数值方法——GEL（Ghost-fluid Euler-Lagrange）方法，编写了GEL二维计算程序。     

基于激光靶耦合问题的二阶精度重映算法

在现有二维激光靶耦合程序LARED-H中，用流体力学方程组描述等离子体的运动，流体计算以拉氏计算为主，经常由于网格的大变形导致数值计算提前终止。ALE方法是目前解决大变形问题中比较流行的方法，在二维空腔计算中，希望能在变形较大的地方采用欧拉区来解决大变形问题。高精度的重映算法是ALE算法中最重要的部分之一，本工作的目的就是研究适用于激光靶耦合问题的高精度重映算法。     

激光靶耦合问题的二阶精度重映算法研究

 任意拉格朗日欧拉方法被认为是解决2维流体力学大变形数值模拟的有效途径之一。以激光黑腔靶耦合问题为背景，结合质量流的物理意义，提出了一种二阶守恒物理量重映的算法，并利用1维激光靶耦合程序对方法及其可行性进行了考察，结合物理问题调整网格，用较少的网格得到了与较多网格拉氏计算同样的结果。为算法在2维LARED-H程序中的应用奠定了基础。     

MOC/SN耦合三维中子输运程序KYCORE开发与初步验证

KYCORE程序是中国核动力研究设计院开发的径向MOC（特征线方法）与轴向SN耦合三维中子输运程序。KYCORE将二维MOC与一维SN通过角通量实现高精度耦合,并通过粗网有限差分实现快速收敛,是目前可工程化应用于三维中子输运计算中精度最高的方法之一。介绍了2D/1D计算与加速理论,并通过与蒙特卡罗程序的计算对比,数值验证了KYCORE三维中子计算的准确性与高效性。     

阻抗梯度飞片加载下的超高速发射二维数值模拟方法

阻抗梯度飞片准等熵加载和超高速发射的二维数值模拟,在计算方法上集中体现了多介质、多界面、大变形、高密度比等特点.采用多介质流体高精度PPM计算方法,以VOF为基础研制MFPPM2二维计算程序,数值模拟Sandia实验室的超高速发射实验模型,获得了与Sandia实验室数值计算一致的结果.     

二维拉氏自适应流体动力学（LAD2D）软件研制

二维拉氏自适应流体动力学软件LAD2D，是采用建立在拉氏自适应结构和非结构网格上的有限体积格式，可以计算平面二维和柱对称二维多物质大变形弹塑性流体动力学问题。LAD2D软件系统主要由5部分组成：主控程序、数据模块、前处理模块、主体计算模块、网格模块和后处理模块。其中主体计算采用了结构网格与非结构网格联合使用的拉氏网格体系，计算格式采用了有限体积格式。网格模块包括网格生成、自适应网格加密（AMR技术）和网格重分技术，以及网格改变后物理量守恒重映技术。LAD2D软件系统由主体程序、二维网格生成程序（GRID2D）、二维自适应网格加密程序（AMR2D）、二维自适应程序（ADAPT2D）f[I-维物理量重映程序（REMAP2D）组成。     

FCT技术在ALE方法中的初步应用

ALE方法在一定程度上克服了纯Lagrange和纯Euler方法的缺点，成为计算大变形问题比较有效的方法。在模拟流固耦合、ICF和磁流体动力学等方面有广泛的应用。对于ALE方法来说，计算网格间物理量的输运（或物理量的重映）是十分关键的。在计算网格间的输运时，用二阶精度的FCT技术，使ALE方法程序计算大变形能力和计算精度都有一定程度的提高。用添加FCT技术后的程序对Taylor杆问题和炸药驱动飞片问题进行数值模拟，Taylor杆问题的计算结果与用LS-DYNA得到的计算结果一致。炸药驱动问题的计算结果与理论分析及实验结果符合较好，比修改前的程序计算变形能力有较大的提高。     

基于旋转场曲率的二维剪切梁单元建模

对二维剪切梁单元进行研究,利用平面旋转场理论推导了精确曲率模型.采用几何精确梁理论构建了剪切梁单元弹性力矩阵.通过绝对节点坐标方法建立了系统的非线性动力学方程,提出基于旋转场曲率的二维剪切梁单元,并分别引入经典二维剪切梁单元和基于位移场曲率的二维剪切梁单元进行比较研究.首先,静力学分析证明了所提模型的正确性;其次,特征频率分析验证了模型可与理论解符合,收敛精度高,并且能准确地预测单元固有频率对应的振型;最后,在非线性动力学问题上,通过与ANSYS结果对比分析,证明了该模型可有效处理柔性大变形问题,并且与经典二维剪切梁单元相比具有缓解剪切闭锁的优势.因此,本文提出的基于旋转场曲率的二维剪切梁单元在处理几何非线性问题中具有较大的应用潜力.     

SPH方法对金属表面微射流的数值模拟

 光滑粒子流体动力学（Smoothed Particle Hydrodynamics）方法是一种无网格的拉氏计算方法,具有计算格式简单、易于计算大变形问题等优势。用二维SPH程序模拟了冲击加载下金属铝表面的沟槽状缺陷产生微射流的过程,得到了喷射物的总质量、最大速度和质量-速度曲线。还计算了相同深度、不同夹角的沟槽微射流,并与实验结果进行了比较。结果表明：计算得到的喷射物总质量和最大速度与实验结果符合较好,SPH方法对于估算金属表面的微射流提供了一种有效的数值模拟手段。     

有限元方法在二维六角形组件堆芯临界问题计算中的应用

应用有限元方法对二维六角组件维芯的少群中子扩散临界计算进行了研究，并编制了二维少群中子扩散方程有限元计算程序２ＤＦＭＥ。以快堆３００ＭＷ－ＬＭＦＢＲ临界问题为例，作了数值计算。且与国内外有限元程序和有限差分程序ＣＩＴＡＴＩＯＮ的计算结果进行了分析比较。     

一个三维台阶流全N-S方程的求解程序

本文用SIMPLE方法,设计了一个求解不可压三维N-S方程的程序,作为算例,求解了一个三维台阶流,并和二维台阶流结果作一比较,分析了三维壁面效应。     

The method of lines solution of discrete ordinates method for non-grey media

A radiation code based on method of lines (MOL) solution of discrete ordinates method (DOM) for radiative heat transfer in non-grey absorbing-emitting media was developed by incorporation of a gas spectral radiative property model, namely wide band correlated-k (WBCK) model, which is compatible with MOL solution of DOM. Predictive accuracy of the code was evaluated by applying it to 1-D parallel plate and 2-D axisymmetric cylindrical enclosure problems containing absorbing-emitting medium and benchmarking its predictions against line-by-line solutions available in the literature. Comparisons reveal that MOL solution of DOM with WBCK model produces accurate results for radiative heat fluxes and source terms and can be used with confidence in conjunction with computational fluid dynamics codes based on the same approach.     

Numerical Simulation of 2-D Radiation-Drive Ignition Implosion Process

A Lagrangian compatible radiation hydrodynamic algorithm and the nuclear dynamics computing module are developed and implemented in the LARED Integration code, which is a radiation hydrodynamic code based on the 2-D cylindrical coordinates for the numerical simulation of the indirect-drive Inertial Confined Fusion. A number of 1-D and 2-D ignition implosion numerical simulations by using the improved LARED Integration code (ILARED) are presented which show that the 1-D numerical results are consistent with those computed by the 1-D radiation hydrodynamic code RDMG, while the simulation results of the 2-D low-mode radiative asymmetry and hydrodynamic instability growth,according to the physical analysis and anticipation, are satisfactory. The capsules driven by the sources from SGII experiments are also simulated by ILARED, and the fuel shapes agree well with the experimental results. The numerical simulations demonstrate that ILARED can be used in the simulation of the 1-D and 2-D ignition capsule implosion using the multi-group diffusion model for radiation.     

用非结构网格与欧拉方程计算复杂区域的二维流动

提出用Delaunay三角化方法生成非结构网格的一种过程。所生成的网格可用于复杂多连通域内的可压流计算。采用Euler方程和格心有限体积法,研制出程序,给出了算例。     

OCDMA系统中2-D码和3-D码

光码分多址(OCDMA)技术主要应用在接入网中。设计大容量的地址码是该技术实用化的前提。本文在构造2 D素数码的基础上提出了一种用于OCDMA系统的空域/频域/时域3 D地址码。这种码基于素数序列运算。理论分析和数学计算结果表明，3 D码比2 D码的容量更大，系统带宽效率更高，误码率(BER)更低，表现出的性能更好。     

Three-dimensional modelling of inductively coupled plasma torches

A three-dimensional model has been developed for simulating the behaviour of inductively coupled plasma torches (ICPTs), using customized CFD commercial code FLUENT ^?. The helicoidal coil is taken into account in its actual 3-D shape, showing the effects of its non-axisymmetry on the plasma discharge. Steady state, continuity, momentum and energy equations are solved for argon optically thin plasmas under the assumptions of LTE and laminar flow. The electromagnetic field is obtained by solving the 3-D vector potential equation on a grid extending outside the torch region. In order to evaluate the importance of various 3-D effects on calculated plasma temperature and flow fields, comparisons of our new results with the ones obtainable from conventional 2-D models and from an improved 2-D model that includes 3-D coil effects are presented. The presence of wall temperature hot spots due to plasma discharge displacement from the torch axis is evidenced, while the use of the new 3-D code for optimization of induction coil geometry and plasma gas inlet features is foreseen. Received 5 September 2002 Published online 13 December 2002 RID="a" ID="a"e-mail: colombo@ciram.ing.unibo.it     

靶与光斑尺寸对激光冲量影响的数值模拟

采用二维雷诺平均N-S方程,数值模拟研究了大气条件下短脉冲激光与固体靶相互作用所产生等离子体的动力学过程。采用k-ε两方程模型用于湍流的数值模拟,分别利用ROE格式和二阶中心格式对对流通量和粘性通量进行离散处理;用高斯-赛德尔隐式格式对方程进行时间推进求解。数值模拟给出了激光引发靶蒸气等离子体侧向膨胀、稀疏等二维流体动力学过程的物理图像,讨论了靶与光斑尺寸对脉冲激光冲量的影响。结果表明,不同宽度固体靶受到的激光冲量有很大差异,固体靶宽度越大,受到的激光冲量也越大。     

二维磁流体力学程序MFCG Ⅱ与MFCG Ⅰ计算结果的对比分析

 根据磁场、电流密度和感应电流的解析公式研制了爆炸磁通压缩发生器充电和爆磁压缩过程的计算程序MFCG Ⅱ。在相同条件下用该程序和二维磁流体力学程序MFCG Ⅰ分别计算相同模型所得结果符合较好。因此用MFCG Ⅱ可分析、解释MFCG实验中有关的物理问题。     

Complex surface three-dimensional shape measurement method based on defocused Gray code plus phase-shifting

Binary pattern defocused projection method can overcome the nonlinear gamma of the projector, as well as filter high harmonics and high-frequency noise. However, high-accuracy three-dimensional (3-D) shape measurement of complex surface using it still remains a challenge. Therefore, a novel Gray code plus phase-shifting method based on defocusing is proposed to solve the problem. The edges of Gray code patterns become blurred owing to defocus, which makes the recovery of accurate Gray code patterns difficult. To solve this problem, the positive and inverse Gray code patterns are projected to obtain threshold values, which are used to achieve the binarization of Gray code patterns. This method is robust and suitable for different defocus levels. Compared with the traditional Gray code plus phase-shifting method, the experimental results prove the effectiveness and feasibility of the proposed method.     

Modeling of humps formation during deep-penetration laser welding

A 3-D transient modeling based on the numerical resolution of the fluid flow ant the heat transfer equations is developed for deep-penetration laser welding at high-welding speed regime, which results in the humping phenomenon. The physical mechanisms included in our model concern matter melting, vaporization inducing a recoil pressure, and resolidification. The implementation of developed procedures called User Defined Functions (UDFs) working interactively with the CFD Fluent code and a dynamic mesh method allowed us to treat the problem with specific and complex boundary conditions. The recoil pressure, fusion, resolidification, and the temperature dependence of the physical properties were thus taken into account. As a result, regular humps could be observed after resolidification on the weld seam.