辐射与物质强耦合问题的二维输运数值模拟

分析确定性输运数值方法模拟辐射与物质强耦合问题的难点,给出在任意四边形网格上求解该类问题的简单隅角平衡算法,并应用灰体输运综合加速算法提高输运方程的源迭代收敛速度,最后给出数值例子.     

强爆炸火球数值模拟中的算子分裂方法

采用算子分裂方法将辐射流体力学方程组分裂为对流项和刚性源项,设计一种高效求解刚性源项方程组的数值方法.数值实验表明:该方法对时间步长不敏感,计算精度能够满足工程计算要求.在不对方程组做任何近似的情况下,数值给出了较长时间内火球冲击波阵面压力峰值及阵面位置的变化,结果与实际强爆炸中的经验公式吻合较好.     

二维柱对称辐射输运基准模型及程序考核

对二维辐射扩散和输运模型问题,分别构造具有含外源项的精确解,解具有一定的物理特征且形式简单,并且辐射输运方程的解可以直接退化到辐射扩散的情形.给出LARED-R程序求解扩散模型数值解及分析.     

激光与靶非平衡耦合的数值模拟

对激光与高Z靶耦合主要物理过程,其中包括激光吸收,电子热传导,等离子体运动,非平衡原子物理,X光转换和辐射输运等现象的数值模拟作了全面的概括性描述。提出了总体微分方程组,建立了差分计算格式和求解方法,给出了程序数值计算的结果。     

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

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

粒子输运离散纵标方程基于界面修正的并行计算方法

为了改造粒子输运方程求解的隐式格式,研究设计适应大型并行计算机的并行计算方法,介绍一类求解粒子输运方程离散纵标方程组的基于界面修正的源迭代并行计算方法.应用空间区域分解,在子区域内界面处首先采用迎风显式差分格式进行预估,构造子区域的入射边界条件,然后,在各个子区域内部进行源迭代求解隐式离散纵标方程组.在源迭代过程中,在内界面入射边界处采用隐式格式进行界面修正.数值算例表明该并行计算方法在精度、并行度、简单性诸方面均具有良好的性质.     

吸气式激光推进中激光能量沉积过程的数值模拟

 结合辐射输运方程,在流体力学方程组的能量方程中加入包括空气吸收的激光能量以及高温气体向周围辐射损失的能量源项,转化为辐射流体力学方程组,建立了用于模拟吸气式激光推进中能量沉积过程的物理力学模型和计算方法。该辐射流体力学计算程序可以很好地模拟激光能量沉积过程中空气对激光能量的吸收、等离子体对激光的屏蔽作用以及激光维持的爆轰波的传播规律,计算得到激光能量的沉积效率约为57%,激光维持的爆轰波的传播速度与同等条件下的理论和实验结果吻合得较好。     

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

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

���·��̵ķ��Ѹ�ʽ�����ڽ�յ������廷��ײ�о��е�Ӧ��

从Lagrange观点出发,采用分裂格式法求解一维辐射流体力学方程组中的三温方程,用于紧凑等离子体环与靶碰撞产生的X射线辐射的数值模拟研究。提出了求解电子、离子、辐射场三温相脱离的能量方程的隐式差分格式,介绍了能量交换项与压力做功及热传导项分开计算的分裂格式方法,数值模拟得到了其对温度变化的影响。     

超强激光与等离子体作用中的能量输运

应用多光子非线性Compton散射模型和数值计算方法,研究了Compton散射对超强激光与等离子体作用中能量输运的影响,提出了将Compton散射光和入射超强光作为电子能量输运的新机制,给出了电子热传导新模型和能量输运数值计算结果。结果表明:散射使等离子体中Weibel不稳定性和自生磁场增强效应导致耦合光传输方向的电子密度显著减小,更多激光能量以热流形式分布在横向方向。散射使电子吸收能量的时间缩短和自生磁场线性阶段最大增长率增大效应导致等离子体表面处沿耦合激光横向方向的热流几乎被完全限制,电子在激光传输方向的能量显著增加。     

Resonant electronic energy transfer from excitons confined in silicon nanocrystals to oxygen molecules

We demonstrate efficient resonant energy transfer from excitons confined in silicon nanocrystals to molecular oxygen (MO). Quenching of photoluminescence (PL) of silicon nanocrystals by MO physisorbed on their surface is found to be most efficient when the energy of excitons coincides with triplet-singlet splitting energy of oxygen molecules. The dependence of PL quenching efficiency on nanocrystal surface termination is consistent with short-range resonant electron exchange mechanism of energy transfer. A highly developed surface of silicon nanocrystal assemblies and a long radiative lifetime of excitons are favorable for achieving a high efficiency of this process.     

Radiative transfer in finite cylindrical media using transformed integral equations

A modified direct integration method is presented to solve three-dimensional radiative transfer in emitting, absorbing and linear-anisotropic scattering finite cylindrical media. This scheme effectively avoids an integral singularity in the coupled Fredholm type integral equations of radiative transfer. The scheme leads to faster and more accurate results, which are needed in combined mode and non-gray problems. The calculated incident radiation and heat fluxes agree well with published results by discrete ordinates method. Using the transformed integral equations, the effects of boundary emission and reflection can also be easily handled.     

辐射传输的格子Boltzmann模拟

辐射动力学理论是描述辐射传输是重要手段,基于此,本文建立了辐射能和辐射动量守恒方程,并基于Chapman-Enskog多尺度展开方法实现了从辐射传输Boltzmann方程到宏观方程的推导,进而建立了适于一维辐射传输的2分量格子Boltzmann模型。数值结果与精确解吻合较好,表明本文提出的LBM方法具有很好的准确性和稳...     

Efficient weakly-radiative wireless energy transfer: An EIT-like approach

Inspired by a quantum interference phenomenon known in the atomic physics community as electromagnetically induced transparency (EIT), we propose an efficient weakly radiative wireless energy transfer scheme between two identical classical resonant objects, strongly coupled to an intermediate classical resonant object of substantially different properties, but with the same resonance frequency. The transfer mechanism essentially makes use of the adiabatic evolution of an instantaneous (so called “dark”) eigenstate of the coupled 3-object system. Our analysis is based on temporal coupled mode theory (CMT), and is general enough to be valid for various possible sorts of coupling, including the resonant inductive coupling on which witricity-type wireless energy transfer is based. We show that in certain parameter regimes of interest, this scheme can be more efficient, and/or less radiative than other, more conventional approaches. A concrete example of wireless energy transfer between capacitively-loaded metallic loops is illustrated at the beginning, as a motivation for the more general case. We also explore the performance of the currently proposed EIT-like scheme, in terms of improving efficiency and reducing radiation, as the relevant parameters of the system are varied.     

Transient coupled heat transfer in a rectangular medium with black surfaces

The main objective of this paper is to extend to two-dimensional (2-D) medium the ray tracing-node analyzing method, which has already been successfully used to solve one-dimensional (1-D) problem of coupled heat transfer in a semitransparent medium. For simplicity, an infinitely long rectangular semitransparent medium with four black opaque surfaces is chosen as our studying object. A control volume method in the implicit scheme is adopted for discretizing the partial transient energy equation. In combination with spectral band model, the radiative heat source term is calculated using the radiative transfer coefficients (RTCs), which are deduced by the ray tracing method. The Partankar's linearization method is used to linearize the radiative source term and the opaque boundary condition, and the linearized equations are solved by the ADI method. Effects of absorption coefficient, refractive index and conductivity on transient cooling process in the 2-D gray rectangular medium are investigated under the condition that the radiation and convection processes cool one side of the rectangular medium while heat the remaining three sides.     

表面不透明三层漫反射散射性介质耦合换热

本文采用射线踪迹、节点分析法研究了三层吸收、各向同性散射性介质层内的一维辐射和导热瞬态耦合换热,复合层表面不透明漫反射,介质层交界面半透明漫反射,且半透明漫反射交界面的反射率采用Fresnel反射定律确定。采用一层和二层辐射能量传递模型跟踪辐射能量在三层介质内的传递,从而推导出辐射传递系数。运用辐射传递系数求解辐射源项,在辐射对流边界条件下、采用全隐格式求解瞬态能量方程,并从机理上研究了辐射和导热耦合换热过程。     

Transient coupled radiation-conduction in infinite semitransparent cylinders

A method is developed to analyze the transient coupled radiation-conduction in infinite semitransparent cylinders surrounded by isothermal black walls. The radiative heat source term is calculated by the radiative transfer coefficients and the transient energy equation is solved by an implicit finite difference method. The radiative transfer coefficients are deduced by use of the ray tracing method in combination with the Hottel and Sarofim zonal method. The effects of the related parameters on the transient radiative heat source and temperature distribution are analyzed. It is found that the peak of the dimensionless radial radiative heat source can be located at the interior shell of the cylinder with small optical thickness when heated by the surrounding irradiation. Treating the volume radiation as a surface radiation will result in large errors of transient temperature distribution for the cylinder with small optical thickness.     

一维瞬态辐射传输的有限元法模拟

随着超短脉冲激光的快速发展,吸收散射性介质内的瞬态辐射传输引起了人们的广泛关注.本文基于离散坐标法和最小二乘有限元法(LSFEM),提出了模拟多维吸收散射性介质内瞬态辐射传输的数值模型.该模型有效地克服了在标准Galerkin有限元法(GFEM)中发生的伪振荡现象,在时间步长较大的情况下仍然可以得到光滑无振荡的解.而且,最小二乘法产生的求解系数矩阵是对称正定的,与GFEM中的系数矩阵相比,仅需要存储一半的非零系数,可以应用许多高效的迭代求解方法进行求解.为了检验模型,本文研究了一维吸收散射性介质内瞬态辐射传输问题,其结果与蒙特卡洛法(MCM)和积分模型法(IE)的结果进行了比较,结果证实:本文的方法可以精确、高效地模拟参与性介质内的瞬态辐射传输.     

用射线踪迹法解黑表面矩形介质耦合换热

本文用射线踪迹-节点分析法研究了二维黑体表面矩形、各向同性散射半透明介质内辐射与导热瞬态耦合换热。采用全隐格式的有限差分法离散二维瞬态微分能量方程,用辐射传递系数来表示辐射源项,结合谱带模型并采用射线踪迹法求解辐射传递系数。采用Patankar线性化方法将辐射源项及不透明边界条件线性化,并采用附加源项法处理边界条件,运用ADI方法求解名以上的线性化方程组,从而解得二维矩形介质内的瞬态温度分布。     

Radiative characteristics of high-frequency electrodeless helium lamps

Based on the system of equations describing the populations of highly excited helium atomic levels and the electron density and energy balance and on the Maxwell equations, we develop a model for an electrodeless high-frequency discharge in helium. We suggest a method of self-consistent calculation in a plasma-field system. The model developed is used to calculate the radiative characteristics of high-frequency electrodeless lamps. The derived dependences of line intensities on gas pressure, lamp radius, and discharge power are in good agreement with the experimental data.