美国恢复核试验可能性分析(下)

3、科技进步为进一步提高和改善武器性能提供的可能性为了在不进行核试验的环境下，维持和确保美国核武器的安全性、可靠性和性能，美国能源部根据总统的指示和国会授权法，于1995年制定并开始实施“核武库维护与管理计划”(SSMP)，提出要研发先进的超级计算机；开发三维高分辨率、高保真、全系     

可压缩多介质流动问题的高精度数值模拟方法

在可压缩流动问题的数值模拟领域,激波的高分辨率计算已取得重要进展。但是包含物质界面的可压缩多介质流动的数值模拟还存在诸多数值挑战,主要表现为界面处数值耗散过大和非物理振荡等问题。界面处流体性质的不连续性是造成可压缩多介质流动问题物理建模与数值方法困难的主要原因。为了建立一套高效的可压缩多介质流动问题的高精度数值模拟方案,本文从数值框架的选择、非守恒方程相容离散、高精度有界格式构造、界面压缩、界面-单介质分区计算方法等多个维度出发,综述近几年我们在可压缩多介质流动问题高精度数值模拟方法方面的研究进展。通过上述多个维度的工作,我们建立了一套适用于可压缩多介质流动问题的低耗散、基本无数值振荡的高精度欧拉数值方法,并成功应用于可压缩多介质大变形流动和界面不稳定性诱导湍流混合等问题的数值模拟。相关数值方法研究成果已集成至武器物理等领域工程数值模拟软件中,为相关工程任务提供了重要技术支撑。     

一种中心型间断有限元MMALE方法

针对多介质可压缩流体动力学问题,提出了一种单元中心型二维MMALE(Multi-Material Arbitrary Lagrangian-Eulerian)方法。在拉氏步,流体力学方程组采用中心型间断有限元方法求解。对于混合网格,采用Tipton压力松弛模型更新物理量,用等参坐标法更新物质中心点坐标。界面重构采用一种健壮的MOF(Moment of Fluid)方法。在重映步提出了基于多边形相交的二阶积分守恒重映方法。该方法分为四个部分:多项式重构、多边形相交、积分和后验校正。多边形相交使用"剪裁投影"算法,显著降低了多边形相交算法的复杂度。后验校正是基于MOOD (Multi-dimensional Optimal Order Detection)限制策略,并做了一些改动以适应多介质的计算。数值算例表明,该方法具有二阶的精度和较好的鲁棒性。     

基于阵列流型分离的稳健自适应波束形成

在较少的采样数或期望信号指向存在误差的情况下,传统的自适应波束形成算法的性能会受到较大影响。基于阵列流型分离的思想,提出一种新的自适应波束形成方法,可以在较少的采样数或期望信号指向存在误差时获得较稳健的性能。数值仿真证明了方法的有效性。     

冲压发动机伴随蒸发过程的两相流流场计算

为了计算固液混合式火箭冲压发动机补燃室内带有燃料液滴的流场,用块隐式法计算气相流场,用分离流连续介质模型计算颗粒相的流动与蒸发过程,在控制方程组的联合求荽中,采用特定的算法来体现颗粒相的质量守恒律与尺寸的变化。算例表明,迭代过程能够快速收敛,计算结果定性正确。     

具有复杂边界的火箭冲压发动机流场数值计算

针对实际复杂的补燃结构和计算边界,以块隐式法计算流场,用代数法生成计算网格,通过改变离散方程来处理三维定常Ｎ－Ｓ方程边界条件。算例表明,该方法正确可靠,迭代过程能够快速收敛。     

混合气体流动的一种数值模拟方法

本文建立一个数值模拟完全气体混合流动的理论模型.该模型首先应用混合气体的Euler方程和每种气体组分的质量分数方程来控制流动.为了消除混合网格内气体组分界面附近出现的非物理振荡,我们假定混合气体的每种组分达到了动力学平衡状态然而尚未达到热力学平衡状态.这种思想导致需要另外给定每种气体组分的总能量方程.为使用高分辨格式来求解这组双曲型偏微分方程并且简化对所需要的Jacobi矩阵的推导,混合气体的压力方程也被耦合起来.Godunov型的波传播方法被采用来离散求解所获得的控制方程.从典型算例结果来看,一维问题的数值解与精确解一致,二维问题的数值解与理论分析一致.这说明本文的理论模型是合理的.     

不同喷压比下欠膨胀射流对特征频率激励的响应

为了揭示欠膨胀激励射流的流动机理,以及考察不同喷压比下射流对相同激励的流动响应,采用大涡模拟方法,对喷压比NPR=5.60和9.34的欠膨胀定常射流和激励射流进行了三维数值计算。激励频率为定常射流中固有的轴对称频率fe=14.569kHz,激励形式为在射流喷管入口处施以正弦压力扰动。结果表明,特征频率激励影响射流的声场特征,缩小射流核心区的范围,减少射流近场的激波胞格数目,并影响射流气体与环境空气的混合。同时,激励射流的特征频率转变为激励频率及其高阶倍频,激励射流的主不稳定模态均为轴对称模态。其中,NPR=9.34的欠膨胀射流的主不稳定模态和外加压力扰动的形式相一致,射流与外加激励发生了更加剧烈的流动耦合和响应。这使得在NPR=9.34时,射流核心区长度减小得更多,压力脉动的振幅更大,激励对射流混合的增强作用更加明显。     

固液火箭冲压发动机内两相反应流场数值计算

为了计算固液混合式火箭冲压发动机补燃室内的三维反应流场，用块隐式法求解气相Navier-Stokes方程组，用连续介质模型和k-ε-Ap模型计算颗粒相的湍流流动与蒸发过程，用修正的k-ε-g模型描述燃料的燃烧。为了分析发动机设计参数对反应流场的影响，用不同的条件进行计算，并由此分析了补燃室几何结构和液体燃料切始颗粒直径对燃烧效率的影响。算例表明，计算方法有效可行。数值结果能够反映流场结构、液体燃料的蒸发和两种燃料的燃烧过程。     

非结构网格上Euler方程的区域分裂算法及并行计算

改进了“波阵面”区域分裂算法,并应用于流场区域的划分;对于子区域边界的不光滑现象,为尽量减少通讯消耗,提出了一种边界并行优化策略。利用PVM并行环境,探讨了非结构网格上求解Euler方程的分区并行算法。根据改进的区域分裂算法及优化策略,运用Jameson有限体积法,对二维翼型流场进行了分区并行求解,多区计算的结果与单区计算的结果作了比较,表明了本文研究方法的有效性。