离散事件系统的Petri网图仿真

1.前言 离散事件动态系统DEDS(Discrete Event Dynamic System)的传统分析方法——应用排队论和随机过程理论为离散事件动态系统建立排队网络模型,在很多情况下已不适用。新的理论和研究方法相继提出。其中之一是“皮特里网”(Petri Nets)理论。另一种比较成熟的方法是仿真法,但是,仿真法要以高昂的计算机设备的大量计算时间为代     

用Newton子迭代方法计算前飞旋翼粘性绕流

通过求解Navier-Stokes方程数值模拟了直升机旋翼前飞非定常流场，为了同时保证计算的时间精确性和计算效率，时间推进格式采用了双时间推进方法，在该方法中，子迭代过程由十分高效的LU-SSOR方法完成，且使子迭代过程成为Newton子迭代，空间上应用中心平均的有限体积法进行离散，为了模拟前飞桨时间的相对运动，网格布局采用了运动嵌套网格方法，应用本文方法对-悬停流场进行了数值计算，计算结果与实验吻合较好，尽管缺乏实验数据的验证，对-有升力前飞状态的数值模拟结果是可信的。     

再入飞行器高超声速绕流的一步和二步有限元计算

本文针对再入飞行器的空气动力问题，采用一步和二步有限元格式求解可缩守恒型Ｅｕｌｅｒ方程，并结合不同的人工粘性模型来捕捉激波，给出了高超声速二维钝头绕流的结果。     

钝锥大攻角超声速分离流场的数值模拟

最近,张涵信等人在传统的Beam-Warming隐式、无迭代、空间推进技术的基础上,根据边界层方程的性质,设计了一种可用小步长推进求解抛物化NS(PNS)方程、而不会引起解的“漂移现象”发生的方法。这种方法对轴对称流动的计算是成功的。本文就是将这一思想推广应用于大攻角有周向分离的流场计算。求解的区域为具有薄亚声速层的有粘与无粘干扰的整个激波层内的流场。在对攻角α=0°和α=20°的球钝锥的计算中,关于壁面上的压力、热流率及流场的涡旋结构均得到了满意的结果。文中特别研究了钝锥大攻角绕流的流动分离图象。 为了增强块三对角矩阵的主对角优势,通常在差分方程的左端附加二阶增量项。本文以选取适当小的推进步长的方法来达到增强主对角优势的目的,不需再附加二阶增量项,从而提高了解的精度。     

用欧拉方程计算机翼与翼身组合体绕流

本文用保角转绘加剪切变换生成三维机翼与翼身组合体Ｃ－Ｈ型网格，并在这类网格拓扑上，用有限体积法研制出可供分析机翼与翼身组合体绕流的三维欧拉方程计算程序。本方法的特色是改进了机翼表面网格点的分布，使机翼后缘与网格线一致；采用了当量机身，使机岙表面到对称面光滑过渡；嵌入了边界层粘性修正。改善了计算结果。数值计算表明，对粘性影响较大的超临界绕流，边界层修正效果显著。     

高速风洞支架干扰数值修正研究

提出了一种计算高速风洞支架系统对飞行器模型纵向气动力干扰量的数值计算方法 ,从跨声速全位势积分方程出发 ,编制了适用于飞行器全机模型及其带支架情况下的跨声速绕流计算程序。通过对双垂尾模型和GBM 0 3模型两个算例的计算 ,讨论了尾支撑位置及其几何外形参数对模型气动力的影响 ,并对GBM 0 3模型带短支杆情况下的纵向实验结果进行了修正。表明该方法对于分析研究风洞模型支架干扰问题并进行支架干扰修正是可行的、有效的 ,可以作为选择尾支撑位置及其几何外形参数和对跨声速风洞纵向实验结果进行支架干扰修正的工具。     

一类拟线性椭圆型方程的CHOQUARD-PEKAR问题

讨论了一类拟线性椭圆型方程的CHOQUARD-PEKAR问题在无界区域中的非平凡解的存在性,对于极小问题Iλ和I∞λ,得到了对于每个λ0,存在α∈(0,λ],使得Iα和I∞α可以达到.     

液氢贮箱停放过程中的力热分析

为了寻找到运载火箭长时间停放过程中液氢贮箱的最经济液位，用于制定合理的发射流程以及紧急处置方法，采用计算流体力学(CFD)技术，对某型运载火箭停放期间液氢贮箱的力、热情况进行了仿真计算和分析。计算选用了VOF(Volume-of-fluid)两相流模型以及Lee相变模型，为了提高Lee模型在不同压力情况下对相变过程的模拟精度，采用安托因方程修正了该模型。修正后的模型首先由试验数据校验了其精确性，随后开展的液氢贮箱停放过程仿真结果表明：贮箱的竖直方向与径向均存在温度分层的现象，液相内会形成大的漩涡，该漩涡会使得冷热流体不断进行热交换，并导致贮箱内部的液氢出现气化。贮箱停放期间蒸发率最大值超过2 m 3/h，发生在停放4 h左右；而贮箱液位充填至37 m 3以上或17 m 3以下时蒸发率较低，最小值接近1 m 3/h。