双时间隐式方法求解非定常绕流的相关问题

空间采用三阶精度的WeightedNND(WNND)格式,时间方向则采用具有二阶精度的双时间步隐式方法求解二维非定常Euler方程,数值模拟了NACA0012翼型俯仰强迫振动流场,讨论了内迭代步数、内迭代收敛判据以及物理时间步长对计算结果、计算效率与计算CPU时间的影响,并将计算结果与实验值进行了比较。     

采用不同化学反应源项处理方法的胞格爆轰数值研究

采用一种改进的化学非平衡流解耦方法对Euler反应流方程进行解耦处理,对流项采用五阶WENO格式进行离散,化学反应源项分别采用梯形公式和α-QSS拟稳态逼近两种方法处理,时间步进采用二阶精度的Runge-Kutta方法,对H2/O2/Ar预混气的胞格爆轰进行了二维数值模拟。两种处理方法都得到了合理的三波点结构,计算结果接近,但α-QSS拟稳态逼近处理方法的化学反应较完全,爆轰波传播速度较高;揭示了爆轰波结构的发展变化过程,横波与横波、横波与壁面碰撞引起的二次起爆对爆轰波的稳定发展至关重要;分析了数值胞格结构,胞格长宽比与参考文献基本一致,验证了本计算方法的有效性。     

振动矩形机翼非线性绕流数值研究

发展了可用于模拟包含运动边界的非定常流动的数值方法。该方法采用非结构动网格技术并构造了一种MUSCL类型有限体积格式,求解考虑了动网格效应的三维非定常Euler方程。非结构网格的变形运动由弹簧近似技术实现。用上述方法计算了作正弦振动的矩形机翼非线性绕流,计算结果与实验结果吻合,得到了流体运动与固体运动耦合的非线性流动特性。计算表明,该数值方法具有较好的精度和可靠性,完全可用于气动弹性或其他流固耦合问题的数值模拟。     

航天测控通信系统可靠性分析的Krylov子空间投影算法

基于Markov模型对航天测控通信系统进行可靠性分析的过程中,若系统中测控通信设备数量较多,模型中的状态空间随设备数量呈指数增长,将会导致数值计算困难.提出了一种基于Krylov子空间技术的可靠性分析方法,将大规模问题投影至小规模子空间中,求得问题的近似解.实验结果证明,Krylov子空间方法的计算速度及精度优于Ross方法和前向Euler法(forward Euler method,FEM).     

一种新的基于方向导数的二维自适应网格生成算法

结合形变原理及网格迭代思想,利用方向导数计算控制函数,提出一种新的二维空间自适应网格生成算法。数值实验表明,该算法能较好地适应解函数的空间剧烈变化。与其他自适应算法比较,其主要优点是该算法逻辑简单,避免了解网格偏微分方程,节约了网格计算时间。     

平面运动激波绕障碍物流动的非结构Euler解

用三角形非结构网格及有限体积法求解Euler方程 ,对运动激波跨越不同截面形状的二维障碍物绕流进行了数值模拟。计算结果给出了不同时刻各流场的复杂结构 ,该结构与实验结果一致。计算表明用非结构算法模拟复杂外形绕流具有简单、准确的特点 ,并能避免尖角处非物理结果的产生     

非定常Euler方程数值计算中高精度格式应用

采用五阶精度加权紧致非线性格式(WCNS)和非定常“双时间步”方法求解非定常Euler方程,模拟NACA0012翼型强迫俯仰振动流场,研究了高精度格式应用到非定常计算时“双时间步”方法物理时间步长、子迭代收敛判据、子迭代步数以及物理时间导数离散方法对计算精度和计算效率的影响.     

圆柱壳体内 X 射线热击波的二维力学计算

依据二维流体弹塑性模型及二维断裂处理方法,研究Ｘ射线热击波在圆柱壳体内的传播规律。计算结果揭示了圆柱壳体内热击波的一些二维效应。     

政治工作研究中的计算实验方法初探

政治工作计算实验方法,以系统科学思想为指导,以数学和计算机科学等学科理论为基础,从信息科学角度描述和解释政治工作实践,旨在更深入地观察和理解政治工作实践.政治工作研究中采用计算实验方法是必要的,也是可行的;计算实验方法是描述、揭示政治工作实践活动规律的一种有效方法.     

典型效应弱化的领导因素

“树典型”是我军思想政治工作的一个重要方法。事实证明,典型效应在我军建设中起着重要的作用。但是,近年来典型效应有弱化的趋势。引起弱化的原因是多方面的,其中领导因素占有一定比重。     

SPH的两种并行计算方案及其在自由表面流动模拟中的应用

针对光滑粒子动力学(SPH)主要计算量是近邻粒子搜索这一特点,提出了一种基于粒子分解的SPH并行计算方案。利用该方案可以方便的将任意串行SPH代码并行计算,而且每一个时间步内的信息传递量只和粒子总数有关,而和粒子的分布无关,因而特别适合于自由表面流动等大变形问题的并行数值模拟。对一个粒子总数为40万的三维溃坝问题的模拟结果表明,此方案能达到的最大加速比约为16,这一结果可能比空间分解方案(不考虑动态负载均衡)更优。     

SPH并行计算方案及其在自由表面流动模拟中的应用

针对光滑粒子动力学主要计算量是近邻粒子搜索这一特点,提出了一种基于粒子分解的光滑粒子动力学并行计算方案。利用该方案可以方便地将任意串行光滑粒子动力学代码并行计算,而且每一个时间步内的信息传递量只和粒子总数有关,而和粒子的分布无关,因而特别适合于自由表面流动等大变形问题的并行数值模拟。对一个粒子总数为40万的三维溃坝问题的模拟结果表明:此方案能达到的最大加速比约为16,这一结果可能比空间分解方案(不考虑动态负载均衡)更优。     

带横向射流的绕后台阶超音速外流场的数值模拟

将一股射流喷入超音速外流是飞行器飞行控制的常用方法。本文通过应用MacCormack 显格式和Ba(?)dwin-Lomax 代数湍流模型求解二维可压缩的RANS 方程,对带横向射流的绕后台阶的超音速外流场进行了数值模拟,并给出了流场的速度向量图、等压线图和等马赫线图。与无喷射的情形相比,台阶底部压力增加了一倍左右。     

WCNS格式在热化学非平衡流数值模拟中的应用研究

高精度格式因其在较低的计算资源消耗下仍能保持足够的精度而在精细流动结构的模拟中具有独特优势。首次将显式5阶加权紧致非线性格式(WCNS-E-5)引入二维高温非平衡流计算。进行了自由流速度5.0~5.7km/s的高焓风洞和高空条件下的圆柱非平衡流场数值模拟,得到了正确的流场参数分布,壁面压力和壁面热流等结果与实验值吻合较好。计算结果表明WCNS-E-5格式较2阶的MUSCL格式具有更好的热流网格收敛特性。分析并针对算例条件初步解决了高阶格式在刚性化学反应流模拟中的鲁棒性问题,为进一步采用高阶格式开展流场结构更加复杂的热化学非平衡流数值研究奠定了良好基础。     

地表目标复合散射的强迫激励法分析

研究了强迫激励方法在电磁散射计算中的应用 ,此方法将地表目标散射分析化为两步 :首先分析地面的散射 ,再将地面的散射作为激励作用于目标 ,从而得到地面背景下目标的电磁散射信号。应用强迫激励方法计算了地表背景下车体的时域电磁散射 ,并与试验作了对比 ,两者比较吻合     

坦克稳定器和电气故障诊断系统

现役59式坦克稳定器和电气系统,由于研制时测试性设计不周,在使用中一旦发生故障,如何迅速排除,在部队一直是个老大难问题.使用本故障诊断系统,能够指导维修人员快速找到发生故障的部件或外部电路;可以在显示器上用颜色作相应变化的电路图,一步步演示出故障压缩过程,直至最后找出故障部位.本文将介绍该系统的功能、技术方案、技术关键及使用特点.     

激波灭火技术的初步研究

讨论了激波灭火技术的基本原理， 给出可行的研究路线。采用新型的计算效率较高的ＥＮＯ差分格式， 从非定常ＥＵＬＥＲ方程出发， 对近地爆炸波产生的流场进行数值模拟， 分析后认为控制爆炸波后风场是关键。     

SAR图像中道路的检测

提出一种两步算法用于从合成孔径雷达 (SAR)图像中无监督地提取线性特征 ,特别是提取公路网。第一步利用局部线特征检测算子提取线基元 ,为此提出一种改进的均值比率 (RoA)算子 ,应用SAR图像的内在统计特性 ,保证较高的准确率 ;第二步引入公路结构特征的先验信息 ,对提取的线基元进行组织 ,获得比较宏观的线性结构。文中的方法通过ERS -1图像加以验证。     

非定常超音速流动数值计算方法比较

为适用于非定常流场计算 ,在保证时、空二阶精度的前提下 ,通过引入变系数残值光顺技术 ,提高了四步Runge Kutta方法的计算效率 ;并针对球锥简谐振动过程 ,数值比较Runge Kutta多步法和“双时间法”的计算精度。