二维三温能量方程的部分Newton迭代法

符尚武等对二维三温能量方程提出了一种9点差分格式，它适用于任意三维网格。但他们用非线性块Gauss—Seidel方法求解所得的非线性代数方程组，收敛得非常慢且经常不得不因为迭代某些次数后仍不收敛而缩小时间步长。     

二维三温能量方程的差分格式及平面靶数值模拟

通过统一二维直角坐标与柱坐标下三温能量方程的差分格式,实现了LARED-H(laser radiationelectron dynamic holehurm)程序既可以模拟正入射激光平面靶过程,也可模拟斜入射平面靶问题,从而扩展了LARED-H程序的功能.使用发展后的LARED-H程序进行了直角坐标系下线聚焦打靶耦合过程数值模拟,给出了不同入射角度下靶等离子体的膨胀过程.     

二维三温能量方程组的高效代数解法

针对二维三温能量方程九点格式离散后形成的非线性方程组,研制了高效求解的代数解法器.主要思想是在部分Newton-Krylov(PNK)方法和Jacobi矩阵自由的Newton-Krylov(JFNK)方法的框架下,结合非精确Newton类方法和预条件Krylov子空间方法进行高效求解.数值结果显示,PNK方法比非线性块Gauss-Seidel方法快6倍以上,在PNK框架下比较了3种预条件子和4种Krylov子空间方法,得出不同组合的最佳方案.还比较了JFNK方法和PNK方法.     

二维三温耦合差分方程组的解法

针对二维三温耦合热传导差分议程组的特点,给出并比较了块G－S算法与ICCG算法,数值结果表明,ICCG算法的收敛速度比块G－S算法的收敛速度可以快数千倍,并且,在二维激光黑腔靶耦合的数值计算中,如果能够实现向量化计算,使用ICCG方法将会大大提高数值模拟的计算效率。     

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

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

二维三温能量方程组离散求解的两个新预处理技术

二维三温能量方程离散后得到的稀疏线性代数方程组中,系数矩阵各行的对角占优性相差十分悬殊,矩阵元素相差也十分大.针对前一问题,提出了改善对角占优性的一个新比例化方法.针对后一问题,利用每次舍弃前计算多个行的技术提出了多行ILUT预条件方法.最后,将对角占优性改善技术、多行ILUT与对角元比例化技术、RCM排序联合使用于实际的能量方程离散求解中,取得了较好的加速效果.     

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

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

基于变分原理的二维热传导方程差分格式

研究二维热传导方程的差分数值模拟.用变分原理在不规则结构网格上建立热流通量形式的差分格式.将热流通量作为未知函数求泛函极值,并与温度函数联立求解.克服通常九点格式用插值方法计算网格边界上的热传导系数和网格结点上的温度所引入的误差.     

二维三温流体力学计算中时间步长的自动控制

研究了二维三温流体力学计算中时间步长的控制问题,提出了控制时间步长的多种约束条件:既考虑了显式流体力学离散方程的CFL(Courant-Friedrichs-Lewy)条件,又考虑了隐式三温能量方程温度的相对变化,还考虑了Lagrange网格密度(或体积)的相对变化以及三温能量方程的迭代次数等.在计算过程中这些约束条件可以随时自动地改变时间步长,以最经济和合理的时间步长完成计算.最后给出了数值实验结果.     

解二维中子输运问题的三角网格配置法

本文研究了二维轴对称中子输运(包括各向同性散射和各向异性散射)问题的一种有限元数值方法。给出了一系列的数值计算结果,并与国内外的SN、间断有限元方法的数值结果及实验结果进行比较,取得了令人满意的结果。     

三维流线上风方法

1引言在流场的数值模拟中，为了消除高雷诺数流动中数值的不稳定和非物理振荡，普遍采用能反映流动物理本质的上风方法。在有限元法中，上风法主要有三种表达形式：（1）采用有上风特性的权函数，如SUPG法山；你对流项和扩散项分别离散，对流项直接采用上游点作为积分点，如MonotoneStreamelineUpwind法［‘］及SkewedPositly，CoeficientUpwind法［’］；由用特征法处理非定常的对流扩散方程，使其具有上风特性，采用传统的Galerkin法离散方程，如TaylorGalerkin法＊。ttice及SchniPke门提出的一种流线上风方法能方便地应用于已有的…     

解流体力学方程组的一种隐式完全守恒差分格式

对Lagrange非守恒流体力学方程组给出了一种隐式完全守恒差分格式,既保证了质量、动量和总能量守恒的差分近似,又能满足内能与动能的平衡特性,提高了数值解的精度。并用该格式对两个可压缩理想流体模型进行了数值计算,并与其它差分格式作了比较。     

色散方程的一类具任意稳定性的显格式

本文对色散方程u_r=au_xxx构造了一类中间层包含四个结点,带两个参数m和θ的三层显式差分格式。当m和θ满足一定的关系时,其稳定性条件为|γ|≤(m+1)/(4(m-1))(|m|>1),从而当取m充分接近1时,可得到任意大的稳定性条件,并且保持截断误差阶不变。数值例子验证了理论分析的结果。     

常微分方程边值问题的高阶三对角OCI差分法

本文给出了二阶线性常微分方程两点边值问题(ODETPBVP)的高阶差分格式构造的基本思想,推导出六阶三对角OCI差分格式,并对端点有奇异性的方程进行了极限值处理,消去了奇异性,对边界层问题采用了非均匀网格上的六阶三对角OCI差分格式。通过大量的数值比较实验表明,这种高阶三对角OCI差分格式能很好地求解奇异性问题,固有不稳定性问题,奇异摄动问题,对生不稳定性问题和振荡性问题。     

求解扩散方程的一类交替分组四点方法

给出了一种交替分组四点格式,该方法绝对稳定,具有并行本性,便于在并行计算机上直接使用.并进行了数值试验.     

三维内肋管内插入螺旋扭带的强化传热实验

本文分别以水和乙二醇为工质，在Ｒｅ数范围为：６００～４００００，Ｐｒ数范围为：５．５～１１０之间，对四根分别插入三种不同扭率螺旋扭带的三维内助管内的换热和流阻特性进行了实验研究。结果表明：三维内肋管内加装扭带的强化传热技术适用于低Ｒｅ数下高Ｐｒ数工质的管内对流换热强化。根据实验值得到了流阻和换热关联式。     

三维对流扩散方程非等距网格上的四阶紧致格式及其多重网格方法

提出了数值求解三维变系数对流扩散方程非等距网格上的四阶精度19点紧致差分格式,为了提高求解效率,采用多重网格方法求解高精度格式所形成的大型代数方程组。数值实验结果表明本文方法对于不同的网格雷诺数问题,在精确性、稳定性和减少计算工作量方面均明显优于7点中心差分格式。     

放射性递次衰变方程的E因子解法

从放射性衰这的一般物理模型出发，利用Ｅ因子函数及其递推关系，给出了求解放射性递次衰变方程的一种新方法，这种解法数学上容易推导，物理上直观易懂，便于计算机程序上实现。