二维高速碰撞问题欧拉数值模拟的混合网格计算

提出了适用于二维平面或轴对称多介质流体力学两步欧拉数值方法中输运计算的混合网格界面处理.在一个混合网格中,将界面近似看作直线.整个方法分为3步:①用混合网格周围的8个网格的介质面积份额确定界面的法线方向;②用混合网格的介质面积份额或体积份额确定界面的直线方程;③用此直线方程求出通过网格边界的流.给出了用此方法所做的测试、数值计算及与其它算法的比较.     

一种基于MOF界面重构的二维中心型MMALE方法

发展了一种基于MOF（Moment of Fluid）界面重构的二维中心型MMALE（Multi-Material Arbitrary Lagrangian-Eulerian）方法.其中,流体力学方程组采用中心型拉氏方法进行离散求解.混合网格的热力学封闭采用Tipton压力松弛模型.混合网格内的界面重构采用MOF方法,并对MOF方法作了简化和改进.重映步采用一种基于多边形剪裁算法的精确积分守恒重映方法.计算了若干数值例子,包括二维漩涡发展问题、Sedov问题、激波与氦气泡相互作用问题、水中强激波与空气泡相互作用问题、二维RT不稳定性问题等.数值算例表明,该方法具有二阶精度,能够计算界面两侧密度比和压力比很大的问题,并且其健壮性优于交错型MMALE方法,适合计算多介质复杂流体动力学问题.     

多介质流体界面不稳定性程序

研制了二维多介质流体程序，主要包括单介质内高精度流体力学计算，多介质混合网格内各种介质输运过程和压力驰豫平衡过程计算、实际状态方程的黎曼解计算。流体计算分别采用高分辨两步PPM（Parabolic Piecewise Method）算法、TVD（Total Variation Diminishing）算法和FCT（Flux Corrected—Transport）算法，流体界面追踪采用VOF（Volume-of-Fluid）。数值求解可压缩多流体方程组和可压缩VOF方程。二维界面追踪分别采用一阶精度Youngs方法和二阶精度Elivira方法，三维界面追踪采用一阶精度Youngs方法，     

点火内爆靶辐射驱动不对称性全过程数值模拟

为满足分层掺杂点火内爆靶辐射驱动不对称性全过程物理分析的需求,在激光聚变二维总体程序LARED集成上发展了辐射输运建模下的多介质ALE方法-RTALE(Radiative Transfer Arbitrary-Lagrangian-Eulerian)。为提高多介质ALE方法的健壮性,发展了驰豫网格重构算法,该重构算法生成的新网格能自适应流场的变化。数值模拟了激波与气柱相互作用的RM不稳定性实验,模拟的气泡变形程度与试验结果基本一致,其中驰豫网格重构算法中的驰豫因子能够很好地反映流场密度梯度。基于辐射多群输运建模的LARED集成程序能够完整模拟辐射驱动不对称性条件下掺杂点火靶二维内爆过程,克服了传统ALE方法计算不下去和算不好的困难,界面变形程度也符合物理分析。     

基于MOF界面重构的多物质ALE方法

提出一种基于MOF(Moment-of-Fluid)界面重构的多物质ALE(Arbitrary Lagrangian-Eulerian)方法.流体力学方程组采用相容有限元方法进行空间离散.提出一种新的二维子网格力学模型,用来计算混合网格中的物理量经过一个拉氏步后发生的变化,混合网格内的界面重构采用MOF方法.提出一种精确积分守恒重映方法.给出数值算例,如空气和水的Riemann问题,Dukowicz问题,水中强激波与空气泡相互作用问题等.结果表明,方法具有较高的精度,能够处理物质界面和网格的大变形问题.     

三维多介质欧拉方法混合网格的有限体积计算格式

在采用Youngs界面重构技术的基础上，对三维欧拉方法混合网格的计算格式进行了研究。运用Youngs技术确定界面后，混合网格内每一部分物质一般不再是正规的六面体结构，可能是非规则的四面体、五面体、六面体或七面体。本文采用对非规则网格适应性很强的有限体积法对每一部分分别进行计算，给出了混合网格内每种物质的压力、速度、能量等的计算公式，比较有效地解决了混合网格的计算问题。     

三维弹塑性流体力学自适应欧拉方法研究

 研究了三维多介质弹塑性流体力学欧拉数值方法的网格自适应技术,编制了可计算多介质弹塑性流体力学问题的三维自适应欧拉程序,解决了计算弹塑性问题时应用网格自适应方法所遇到的一系列关键问题,如网格细分和合并规则,父网格或子网格物理量的填充,幻影网格的选择、填充及存储方案,插值函数的选择及如何选择时间步长等问题。给出了侵彻、爆轰等弹塑性问题的数值算例,验证了方法的有效性。     

二维欧拉多介质流体力学数值模拟-TVD＋VOF方法

用TVD格式结合VOF界面处理方法编制了二维多介质高分辨欧拉程序，以解决冲击波和多介质界面处理。程序包括单介质网格高精度流体力学计算、多介质网格内界面重构、各种介质输运和压力驰豫平衡过程。其中单介质网格的计算采用Harten二阶TVD格式结合MacCormark方法计算含有源项的非齐次守恒定律方程组，通过4节点限制函数保证格式单调。多介质网格采用Youngs方法构造界面，采用x，y方向分裂格式计算体积份额输运，再根据体积份额输运计算质量、动量和能量的输运，最后利用等熵条件计算各种介质的压力驰豫平衡过程。     

多介质流界面高精度自适应欧拉算法

采用欧拉网格自适应算法捕捉多介质流界面,获得了高精度界面特征,对不同物质引入不同位标函数跟踪界面运动,将位标函数方程与流体动力学方程非耦合求解,在笛卡尔坐标系中运用二阶精度有限体积算法,在保持流场守恒条件下,采用多层网格级对笛卡尔网格嵌套细化,实现了多介质流物质界面的高精度自适应跟踪.方法逻辑简单,大大节省了CPU时间,且能够对局部参数急剧变化的流场(如激波)进行自适应跟踪.     

混合Krylov子空间算法及其应用

给出了一种适合二维三温辐射流体力学能量方程的大型稀疏线性代数方程组的混合迭代算法.计算结果显示,该算法解二维三温辐射流体力学能量方程的大型稀疏线性代数方程组比原有算法快4倍左右;原有算法不收敛时,该算法收敛;各物理量也符合得很好.     

Interchange-turbulence-based radial transport model for SOLPS-ITER: A COMPASS case study

Mean-field plasma edge transport codes such as SOLPS-ITER heavily rely on ad-hoc radial diffusion coefficients to approximately model anomalous transport. Such coefficients are experimentally determined and vary between different machines, and also depend on the operational regime and plasma location within the same device. Therefore, to match experimental data the modeller is required to manually tune several free parameters in expensive simulations, and the code's predictive capabilities are significantly downgraded. As a solution, a new model has been developed for SOLPS-ITER, solving an additional transport equation for the turbulent kinetic energy k, derived by consistently time-averaging the Braginskii equations, and including a diffusive closure for the anomalous particle flux. This closure model relates the anomalous diffusion coefficient to the local k value. The resulting equation structure and its closure are inspired by TOKAM2D isothermal interchange turbulence simulation results. Within this model, fewer and hopefully more universal free parameters are retained, thus improving the code's predictive capabilities. The new model has been tested on a COMPASS case for which upstream plasma profiles were available. Experimental data and a reference solution, obtained by matching the profiles through manual tuning of radial diffusivities, have been used to estimate the parameters of our new transport model. A ballooned particle diffusivity profile is retrieved by the new radial transport model, thanks to the proposed interchange drive. The obtained upstream profiles qualitatively agree with the experiment and prove the new model is a promising first attempt to be further refined.     

Measurement of thermal parameters of a heat insulating material using infrared thermography

The article presents results of research developing methods for determining thermal parameters of a thermal insulating material. This method applies periodic heating as an excitation and an infrared camera is used to measure the temperature distribution on the surface of the tested material. The usefulness of known analytical solution of the inverse problem was examined in simulation study, using a three-dimensional model of the heat diffusion phenomenon in the sample of the material under test. To solve the coefficient inverse problem an approach using an artificial neural network is proposed. The measurements were performed on an experimental setup equipped with a ThermaCAM PM 595 infrared camera and a frame grabber. The experiment allowed verification of the chosen 3-D model of the heat diffusion phenomenon and proved suitability of the proposed test method.     

Thermal diffusivity measurement of insulating material using infrared thermography

The article presents the results of research developing methods for determining the coefficient of thermal diffusivity of thermal insulating material. This method applies a periodic heating as an excitation and an infrared camera is used to measure the temperature distribution on the surface of tested material. In simulation study, the usefulness of known analytical solution of the inverse problem was examined using a three-dimensional model of the phenomenon of heat diffusion in the sample of tested material. To solve the coefficient inverse problem, an approach using artificial neural network is proposed. The measurements were performed on an experimental setup equipped with a ThermaCAM PM 595 infrared camera and frame grabber. The experiment allowed to verify the chosen 3D model of heat diffusion phenomenon and to determine suitability of the proposed test method.     

间断有限元方法求解一维非平衡辐射扩散方程

研究一维非平衡辐射扩散方程的数值方法.通过求解间断系数热传导方程的广义黎曼问题,得到一种带加权数值流量,基于该数值流量构造了一类新型的间断有限元方法.在时间离散上采用向后Euler方法,形成的非线性方程组采用Picard迭代求解.数值试验表明该方法具有捕捉大梯度的能力,而且能适应扩散系数间断的情形.     

热辐射输运问题的高效蒙特卡罗模拟方法

惯性约束聚变研究中,热辐射光子在介质中的输运以及热辐射光子与介质的相互作用是重要研究课题,蒙特卡罗方法是该类问题的重要研究手段之一.隐式蒙特卡罗方法虽然能正确地模拟热辐射在介质中的输运过程,但当模拟重介质(材料的吸收系数大)问题时,该方法花费的计算时间将变得很长,导致模拟效率很低.本文以离散扩散蒙特卡罗方法为基础,开发了"离散扩散蒙特卡罗方法辐射输运模拟程序",可以较好地解决重介质区的计算效率问题,但是离散扩散蒙卡罗方法在模拟轻介质区时精度不够高.辐射输运问题中通常既有轻介质也有重介质,为了能同时解决蒙特卡罗方法模拟的效率和精度问题,本文研究了离散扩散蒙特卡罗方法与隐式蒙特卡罗方法相结合的模拟方法,并提出了新的扩散区与输运区界面处理方法,研制了混合蒙特卡罗方法的辐射输运模拟程序.典型辐射输运问题模拟显示:在模拟重介质问题时,该程序能大幅缩短模拟时间,且能取得与隐式蒙特卡罗方法一致的结果;在模拟轻重介质均存在的问题时,与隐式蒙特卡罗方法相比,混合蒙特卡罗方法的模拟精度与其相当且计算效率同样能够得到显著提升.     

两相湍流色噪声扩维法PDF模型及数值模拟

由颗粒运动的朗之万方程出发,对流体脉动速度采用扩维方法,得到两个不同层次的PDF输运方程.通过对颗粒运动方程求解和高斯分布假设,解决PDF方程的封闭问题,获得颗粒二阶矩模型,然后将颗粒应力方程简化成代数方程,建立代数应力模型.将对流扩散方程的有限分析法运用到求解两相流模型中,对壁面两相射流进行数值模拟,对比分析数值结果与实验结果.     

引入源-样品距离影响的主动中子多重性质量反演公式北大核心CSCD

主动中子多重性计数测量方法是常用的核材料质量无损测量方法,已广泛应用于核材料衡算、核安保测量与军控核查等领域。我们通过对JMCT中子-光子输运程序的二次开发,实现了对经典点模型铀样品质量估算实验的数值模拟,并提出了改进的铀样品质量计算公式。该算法可以显著降低本实验中源-样品耦合与源中子反照等作用对铀样品质量估算精度的影响。建立了主动中子多重性计数测量探测系统模型和32个铀样品半球壳模型,模拟得到了与铀样品距离不同的DT源和AmLi源主动中子多重性计数,利用数值模拟手段检验了质量估算算法的有效性。数值模拟结果表明,改进的铀质量估算算法可以使质量估算的平均偏差率降低到10%以下。     

Grey and multigroup radiation transport models for two-dimensional stochastic media with material temperature coupling

Current theories for approximating the effects of stochastic media on radiation transport assume very limited physics such as one dimension, constant grey opacities, and no material energy balance equation. When applied to more complex physical problems, the standard theory fails to match the results from direct numerical simulations. This work presents the first direct numerical simulations of multigroup radiation transport coupled to a material temperature equation in a 2D stochastic medium that are compared to closures proposed by various authors. After extending it from grey to multigroup physics, one closure that is not commonly used successfully models the results in dilute systems where one material comprises less than 5% of the total. This closure is more accurate for related grey transport problems than it is for the multigroup problem. When the specific heats are material- and temperature-dependent, it is much more difficult to fit the direct numerical solutions with an approximate closure.