水平不变海洋声道中的WKBZ简正波方法

本文利用本征函数的ＷＫＢＺ近似和海面相移修正发展了一种新的计算海洋声场的数值方法，该方法具有形式简明、易于计算、精度满足实用要求的优点．本文给出了ＷＫＢＺ方法应用于北太平洋声道中声传播的一些数值结果，这些结果表明该方法是一种计算分层海洋声道会聚区声场的快速而精确的数值方法．     

多层线型分割法分析渐变型平板介质波导

提出了一种求解渐变折射率平板光波导的传播常数和场强分布的数值计算方法，计算了几类典型的渐变折射率波导的传播常数数值，并与其它方法进行了比较，结果表明，该方法计算的精度极高，速度较快，有实用价值。     

高空核电磁脉冲模拟波形的双指数函数拟合法

 许多标准和公开出版物中都用双指数函数描述高空核电磁脉冲典型波形。通过数值方法，研究了双指数函数一项重要的性质。根据该性质，详细讨论了高空核电磁脉冲模拟波形数值拟合中，双指数函数特征参数与脉冲峰值、前沿、后沿以及半宽等物理参数的关系，从而提出一种简单有效的脉冲参数计算方法。通过IEC标准中规定的高空核电磁脉冲参数的估计与一个实测高空核电磁脉冲模拟波形的数值拟合，验证了该方法的有效性和可靠性。该方法能够适应多种峰值、前沿和半宽的高宽核电磁脉冲模拟波形的数值拟合。     

束流输运理论的新表述（II）——线性力的检验

 用数值方法在恒定和周期线性力情况下对纯粹和混合系统检验了我们在文献[1]中所提出的输运理论框架的可靠性和优越性。结果表明，该理论的数值计算和常用的粒子跟踪模拟完全一致，但计算时间少，误差涨落小。     

二维各向异性浮雕型光栅的矢量衍射分析

利用严格模式理论对二维电磁各向异性浮雕型介质光栅的衍射特性进行了理论分析，并对其退化情形进行讨论。文中采用反射-透射系数矩阵算法和校正傅里叶展开规则以提高数值计算的稳定性、收敛性和计算效率。通过对光栅进行数值计算表明，该方法在可靠性、稳定性和收敛性方面是令人满意的。     

数值求解瞬态温度场的有限体积法

本文介绍了求解三维瞬态热传导的有限体积法,给出了控制方程和数值方法。通过一些算例对该方法进行了验证,数值试验结果表明该方法计算稳定、可靠。     

煤灰沉积的传热过程模型及其数值研究

计算流体力学（ＣＦＤ）方法的应用在锅炉设计或燃烧设备的改造过程中有着十分重要的作用．本文研究了实际燃烧过程中普遍存在的煤灰沉积现象对数值计算结果的影响，提出了描述煤灰沉积的新型传热模型，比较了新模型采用前后数值计算结果与实测数据的差异，从而验证了该模型的合理性，提高了对炉内积灰、结渣过程数值描述的精度．     

对旋轴流水轮机的非定常数值模拟研究

为了研究采用对旋转轮的新型轴流式水轮机的非定常性能,应用三维非定常数值计算方法,对该新型水轮机进行了数值模拟研究.数值计算中,采用了全三维全流道的数值计算方法,基于DES(Detached Eddy Simulation)方法获得了对旋轴流水轮机的非定常流场和能量性能.计算结果显示,该新型水轮机在一个转动周期的不同时刻,其内部流场存在显著的非定常特性.     

激光与靶相互作用体汽化模型的数值计算方法

利用克劳修斯-克拉珀龙(Clausius Clapeyron)方程描述相变过程,利用von Neumann差分方法构造辐射流体动力学方程组的计算方法,数值求解了激光辐照金属材料的热效应、汽化过程和激光维持爆震波,并同作者前期的工作作了比较.结果说明,该计算方法对激光与靶材料相互作用体汽化模型的数值计算是行之有效的,且该汽化模型为相关问题的二维、三维问题的数值计算提供了简便有效的途径.     

光脉冲传输模拟的快速数值差分递推算法及其应用研究

利用零邻域的马克劳林展开对归一化非线性薛定谔方程的频域差分形式进行分析，得出一个能够同时考虑光学媒质的色散作用以及非线性克尔作用的时域快速数值差分递推公式.选取若干算例，将运用该公式的数值计算结果与已知的解析结果和传统分步傅里叶方法所得的计算结果做了相应的比较，结果表明这种快速数值差分递推公式不但拥有相当快的计算速度，也有很高的计算精度，而且在物理上符合光脉冲在光学媒质中传输时色散和非线性作用同时施加影响的客观实际，这说明它是研究非线性光学媒质中光脉冲传输的一种科学、合理而快速有效的数值计算方法. 关键词： 快速数值差分递推公式 零邻域马克劳林展开 非线性薛定谔方程 光脉冲传输     

轴对称渐变型类周期慢波结构的色散特性

 运用场匹配法和傅里叶级数理论，提出一种原则上可数值求解任意轴对称渐变型类周期慢波结构色散特性的方法。采用该方法编制了计算渐变型波纹波导和渐变型盘荷波导色散曲线的Matlab程序，详细分析并讨论了这两类典型渐变型类周期慢波结构的色散特性。数值计算结果与多维全电磁模拟软件模拟结果的数据吻合度较高，验证了该数值算法的可靠性。另外，该方法具有较强的普适性和扩展性，也可退化到任意轴对称周期慢波结构色散特性的求解，为慢波结构的设计提供一种简单有效的途径。     

激光匙孔焊接过程两传热数值计算

对高能量密度激光匙孔焊接过程的二维准稳态传热进行了计算方法的研究。数学模型的 质量、动量和能量守恒为依据，数值计算在曲线坐标系上进行。提出了位置预测－－修正技术，以便准确地捕捉焊过程中固－液交界面的形状。并应用该计算方法对一典型的ＡＩＳ１３０４不锈钢板的匙妆过程中进行了算例验证，结果表明该计算方法进合理可行的。     

双重弹性波波场分离数值模拟及相关理论证明

提出一种等效的双重弹性波波场分离数值模拟方法,用于模拟纯纵波和纯横波分离模式的质点振动速度、位移以及散度场和旋度场,并将该方法应用于全弹性波波动方程数值模拟中.同时,详细推导双重弹性波波场分离波动方程的高阶交错网格有限差分数值计算公式及其稳定性条件、数值频散关系和完全匹配层(PML)吸收边界条件.理论分析和数值计算均表明,该方法可以实现高精度双重弹性波波场分离数值模拟,且纯纵波和纯横波得到完全分离,边界吸收效果较好.与前人工作相比,存储量和计算时间均得到有效改善,数值计算结果进一步验证了该方法的优越性.     

电磁脉冲作用下地面铺设屏蔽电缆蒙皮电流分布的时域计算方法研究

为克服传输线理论方法求解电磁脉冲电缆耦合问题存在的不足，对电磁脉冲电缆耦合问题进行了时域求解方法的研究。引入了适用于电磁脉冲电缆耦合问题数值求解的基于细线散射的时域有限差分(FDTD)法，将Gedney's PML进行了扩展，使其可以用于各向同性有耗计算域截断问题的处理。利用建立的数值方法，模拟了电磁脉冲作用下地面铺设屏蔽电缆引起的蒙皮感应电流分布规律及波形特征，并利用辐射波电磁脉冲模拟器对该问题进行了实验研究。计算与实验结果在误差范围内的一致性证明了数值结果的可靠性，以及该时域计算方法在处理电磁脉冲电缆耦合问题的可靠性。     

任意各向异性介质光栅的严格模式理论分析

将“逆规则”傅里叶级数展开法和反射透射系数阵递推算法进行推广。对任意面形及任意入射角的电磁各向异性介质光栅衍射特性进行研究，该方法物理概念清晰，数值计算稳定收敛，计算速度快，并与其它方法的计算结果进行了比较。     

多层线型分割法分析渐变型平板介质波导

提出了一种求解渐变折射率平板光波导的传播常数和场强分布的数值计算方法,计算了几类典型的渐变折射率波导的传播常数数值,并与其它方法进行了比较,结果表明,该方法计算的精度极高,速度较快,有实用价值.     

数值模拟铅蓄电池在SFUDS循环下的行为

本文应用传热传质理论，提供了一种用于模拟电动汽车用铅酸蓄电池工作的计算机数值计算方法。该数值方法 采用高级计算流体动力学技术，结合铅酸蓄电池内部的电化学过程和传质过程进行计算，可用于预测电动汽车用铅酸蓄电池在ＳＦＵＤＳ循环下的非稳态行为，如酸浓度，电极的多孔率，活性物质利用率，以及电池在放电、休息、充电的循环过 程中的充电状态等各种性能。模拟结果显示电池的负极有５０％的活性物质未充分利用。     

自然对流换热的高精度非结构网格有限体积法

用非结构网格有限体积法求解自然对流换热时,传统的对流项离散格式难以兼顾数值精度与计算效率,我们发展了一种耦合高精度格式的延迟修正方法,用于对流项的离散.高Re数下方腔驱动流数值计算验证了该方法具有较高的计算精度和较好的稳定性.Boussinesq流体的自然对流换热数值模拟,表明该方法能有效克服高Ra数时数值计算发散,可准确捕捉自然对流换热问题中不同偏心率下的等温线和流线分布特征.     

一种计算色度学三刺激值加权表的新方法

三刺激值计算是求取各种色度量的基础。首先回顾了国际照明委员会三刺激值的基本定义及其数值计算方法；其次简单描述和分析了相关工业界常用的三种ASTM加权表及Venable最优加权表；最后，基于Venable的工作，从减小光谱反射比测量误差的影响出发，提出了一种在形式上不同于Venable的最优加权系数选择准则，给出了一种新的计算最优加权表的方法。该方法与现有的计算加权表方法相比，具有不需要通过迭代求解、易于程序设计从而比较简单且有可能作为标准方法予以推广应用等特点。数值仿真表明用此方法计算的加权表较ASTM三种加权表更准确。可以期望，给出的计算最优加权表的方法将会广泛应用于相关工业领域。     

一种分段式γ扫描系统探测效率响应函数的数值计算方法

基于SGS探测系统中放射性核废物桶、准直器及HPGe探测器的空间几何分布建立数学模型。对低、中放射性核废物桶中点源137Cs的探测效率响应函数进行了数值公式推导，然后利用蒙特卡罗工具包MCNP模拟计算了相应的效率函数。将这两种方法模拟的结果进行了对比分析发现:数值计算方法的数据普遍高于MCNP模拟的数据，最大相对误差达143.26%，误差平均值为37.15%。对数值计算方法进行修正后，最大值相对误差为17.22%，平均值为4.54%。在误差范围内，可以认为该修正方法有效。