一维球几何菱形差分中子输运方程的扩散综合加速解法

本文导出了一维球几何定态中子输运方程菱形格式的扩散综合加速方程,并给出了差分公式。所给出的加速方法可以加速菱形格式的输运方程的迭代求解。并给出了部分模型的数值计算结果。     

自适应有限元方法和后验误差估计

自适应有限元方法是科学研究和工程设计领域中非常有效的一种求解偏微分方程的数值计算方法。这种方法是为了以尽可能小的代价取得尽可能好的计算效果。后验误差估计是实现自适应有限元计算的关键性手段。文章综合介绍了自适应有限元方法和后验误差估计在求解椭圆型方程、抛物型方程和双曲型方程方面所取得的比较新的成就。     

综合核方法求解辐射输运问题的求积组选取

简要介绍求解辐射输运方程的综合核方法,分析计算误差和收敛性,提出新的求积组和误差修正方法,提高综合核方法的计算精度.通过对基准问题的计算比对表明,采用提出的求积组并通过误差修正,综合核方法在低阶时的结果具有较高的计算精度.     

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

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

湍流扩散燃烧的数值研究-PDF方法和火焰面模型的性能比较

采用标量概率密度函数(PDF)方法、稳态和非稳态火焰面模型三种方法对一个值班湍流CH_4/O_2/N_2射流扩散火焰(Sandia Flame D)进行数值计算,以比较不同燃烧模型的性能。PDF方法通过计算反应标量的PDF输运方程来得到标量分布,而火焰面模型只求解单标量混合物分数的PDF方程,组分和温度分布通过火焰面方程的求解或者火焰面数据库的插值得到。计算结果和实验数据对比表明PDF方法计算结果最好但计算量相当大,稳态火焰面模型则反之。综合而言,非稳态火焰面模型的预测结果相对稳态模型有了非常大的改进,而计算量仍然容易接受,非常适合工程应用。     

用改进PT方程计算低温流体汽化潜热

本文用改进ＰＴ方程对２２种常见低温流体的汽化潜热进行了计算．计算结果与原始ＰＴ方程、原始ＰＲ方程和改进ＰＲ方程计算结果的比较表明，用改进型ＰＴ方程所得的结果与实验值的偏差明显小于原始ＰＴ方程、ＰＲ方程和改进ＰＲ方程，特别是该方程对强极性物质、低温流体的汽化潜热计算，优越性更为显著．     

用改进的PT方程计算低温流体混合物的汽液平衡

本用改进的ＰＴ状态方程对２１种常见的低温流体混合物的汽液平衡进行了计算，并与原始ＰＴ方程，ＰＲ方程琢ＰＲＯＭＩＸ软件的计算结果进行了比较。结果表明，用改进的ＰＴ方程计算所得结果与实验值的偏差优于原始ＰＴ方程，ＰＲ方程及ＰＲＯＭＩＸ软件的计算结果。     

基于响应面模型的高超声速低温喷管优化

 建立了高超声速低温(HYLTE)喷管的响应面模型,样本点设计矩阵采用正交试验设计法构造,样本数据通过3维粘性NS方程计算喷管光腔耦合段的多组分混合流场获得。基于HYLTE喷管响应面模型,对喷管进行了多目标优化,定义了综合性能参数,优化后喷管混合性能显著改善,综合性能明显提高。计算表明,基于响应面模型的优化策略,能够满足喷管优化对计算精度和计算量的要求。运用多学科设计优化软件框架iSIGHT作为建模和优化的辅助工具,提高了优化设计的实现效率。     

测量激光聚变靶丸内燃料气体总量的对比状态法

推导了用对比参数表示的平衡蒸汽压力方程和平衡蒸汽密度方程，这些方程适用于D2在临界点以下的行为，计算表明，用对比状态计算露点温度下的蒸汽密度比露点法准确，用BWR方程计算充气压力比用理想气体方程合理。     

高空核爆炸火球的二维辐射流体力学计算

本文介绍了高空核爆炸火球的二维辐射流体力学计算的主要计算方法和部分计算结果,用柱对称最大熵变Eddington因子P-1近似方程描述了高空核爆炸的辐射输运过程,讨论了P-1近似的边界条件,对能量方程和P-1近似方程的耦合,给出了一个新的差分格式,提高了耦合计算的精度,显著地减少了计算量,同时有效地克服了辐射输运计算中出现的非物理振荡现象。     

泄漏迭代法在求解中子扩散方程中的应用

本文着重介绍了用于求解中子扩散方程的泄漏迭代方法及其边界条件的解析处理,并将此方法用于反应堆设计的计算,结果表明,这种泄漏迭代方法可以满足反应堆的设计要求。     

THE TRACK SIMULATIONS OF CHARGED PARTICLES IN THE REAL MAGNETIC FIELDS

A simulation method for trajectories of charged particles in the real magnetic fields is introduced for the forward spectrometer in experiments of high energy physics.By using the numerical integrations to the Lorentz equation,simulation tracks with good accuracy have been obtained.     

玻璃各温度段热膨胀系数的一种计算方法

本文提出一种计算玻璃各温度段热膨胀系数的方法。该方法是根据同类玻璃的组成及实验数据,利用加和性原理,写成线性方程组,利用计算机求其计算系数,然后,再求已知配方玻璃的热膨胀系数或已知各温度段的热膨胀系数,求玻璃的配方组成。     

关于推广的组态混合法(Ⅰ)

本文提出了一个借助无规位相近似法的解,应用推广的组态混合法以求得更精确解的方法,导出了有关矩阵元的一组明显表达式。此外并证明了推广的组态混合法是自洽的,给出了有关的格林函数所满足的运动方程及展开式,根据后者对推广的组态混合法的解的性质进行了讨论。     

抛物方程基于自然边界归化的耦合法

将冯康和余德浩提出的自然边界归化方法^[1~4]应用于求解抛物方程初边值外区域问题,提出一种自然边界元与有限元耦合算法。先将控制方程对时间进行离散化,得到关于时间步长的离散化格式,给出圆外域上的自然积分方程,基于此研究抛物方程无界区域问题的自然边界元与有限元耦合法,最后给出相应的数值例子。     

撕裂模研究中一类奇异边值问题的数值解法

在数值求解自洽螺旋平衡方程的基础上,计算了撕裂模不稳定性对电流分布的依赖关系,考虑了外螺旋场对安全因子的影响。我们采用了动态交替隐式法来加速平衡方程的收敛,用差分方法有效地处理了具有两个奇点的奇异边值问题。     

自由界面等离子体平衡的变分原理

在自由界面情况下用求位能汛函极值的变分原理推出了等离子体平衡所需满足的平衡方程和边界条件;其泛函的Euler方程即为具有自由界面条件的磁面函数方程。也给出了便于计算的变分泛函,它相当于等值面边界问题。对于导体壁为简单几何形状的情况,用Rit法进行了数值计算。 关键词：     

气体放电击穿过程的物理和数值研究

本文对低气压（１０＾－２Ｐａ）热阴极气体放电的击穿过程给出了物理描述和相应的双流体数学型，并发展了一种选择和调整未知初始条件的有效算法，可以通过伴随试射法得到对初始条件十分敏感的非线性两点边值常微分方程组的数值解，从而给出这类气体放电中击穿过程的定量描述。     

化学反应体系中涨落的时间空间关联(Ⅰ)——涨落、扩散和波

本文着重分析了如何描述化学反应体系中密度涨落的空间关联,对照反应扩散方程,我们在密度涨落分布函数所满足的Master方程中引入了扩散项或压力项,对无限介质,解Master方程,得到了一系列波,第一支波是与反应扩散方程的解一致,描述了平均密度起伏在空间的传播,在气体中就是声波,其它各支波,则描述了局部密度涨落高阶矩的变化在空间的传播(平均密度不改变),其实就是局部机率分布函数畸变后的传播过程,我们称之为涨落波,文中还讨论了如何求解线性初值问题。 关键词：     

非平衡系统Master方程的稳定性

针对非平衡统计中出现的多元线性Master方程,利用“熵产生”和“剩余熵产生”的概念讨论了Master方程在线性平衡区和非线性远离平衡区的稳定性问题。从而得到与Prigo-gine宏观热力学理论中一致的结果。此外还提出了Master方程所决定的“概率流”的概念,给出了概率流分解的具体解析表达式。 关键词：