显隐式分析方法在复合材料失效分析中的应用

针对复合材料失效等准静态问题,提出了显隐式联合的求解方法:在加载前期,非线性程度较低时,采用计算效率较高的隐式求解器求解;在材料出现损伤后,非线性程度较高时,改用条件收敛的显式求解器推进计算,中间通过数据传递将材料的应力应变水平及材料损伤状态由隐式求解器传递到显式求解器。该方法在保证计算精度的前提下,兼顾了隐式与显式分析的优点,有效提高了计算效率。     

二维带动力吸气式高超声速飞行器绕流的PNS-NS混合求解

为了提高吸气式高超声速飞行器绕流的求解效率,采用空间推进求解抛物化PNS(parabolizedNavier-Stokes)方程和时间迭代求解Navier-Stokes(N-S)方程的混合计算流体动力学(CFD)方法来求解高超声速飞行器整机绕流.在超声速占主导的流动区域采用空间推进求解抛物化N-S方程的方法,在亚声速和分离区采用时间迭代求解N-S方程的方法.对于求解二维带化学反应的吸气式高超声速飞行器绕流,混合CFD方法和完全时间迭代方法相比,可得到同等准确的数值模拟结果,并且求解效率提高了数倍.      

潜入式喷管喉衬界面间隙优化设计

为了解决固体火箭发动机潜入式喷管喉衬在热载荷与内压联合作用下界面间隙优化问题,数值模拟方法求解了喷管热结构响应.采用流体仿真软件,计算了潜入式喷管流动过程的稳态流场,获取了高温燃气流动参数与固体壁面对流换热系数.并采用三维有限元结构计算平台,编写了非均布壁面压力载荷与非均布对流换热系数子程序,求解了喉衬前后搭接界面间隙...     

滞止式温度传感器动态理论建模研究

根据滞止式温度传感器的换热原理建立了动态理论传热模型,采用数值计算方法对模型方程进行了求解,并对计算结果进行了分析.     

喷管化学反应流的对角化点隐式法求解

采用弱耦合对角化点隐式方法的ＭａｃＣｏｒｍａｃｋ格式,求解喷管的粘性化学反应流动,解决了控制方程组的刚性问题,又大大减少了常规点隐式方法中求逆矩阵的大量工作。对液氢液氧火箭发动机,采用６种组分、８个反应有限速率的化学反应模型和Ｂａｌｄｗｉｎ－Ｌｏｍａｘ代数湍流模型,得到流场参数在喷管中的分布。计算结果表明,采用对角化点隐式方法求解化学反应喷管粘性流动能提高计算效率。     

浸入式贴体网格边界方法

提出一种基于重叠网格的浸入式边界方法,用以模拟复杂的三维黏性跨声速流动。该方法采用一套固定的笛卡儿正交网格,用以主流的求解;采用一套可以移动的贴体等距面网格,用以拟合或者离散物面的作用力,通过空间插值,实现两套网格重叠部分的信息传递。分析了浸入式贴体网格边界方法的优势,介绍了贡献单元的寻找策略和物理通量的插值方法。流场求解采用Spalart-Allmaras带湍流模型的Navier-Stokes方程组,其中对流项采用流通矢量分裂和5阶WENO(weighted essentially non-oscillatory)-Z格式离散,黏性项采用6阶中心差分格式离散,时间项采用龙格库塔显式格式离散。数值验证算例表明:该方法具备高于4阶的空间求解精度,并适用于刚体动网格非定常流场模拟,且无需更新网格形状。等距面贴体网格生成过程简单,避免了繁琐的人工调节过程,与笛卡儿网格结合,可提供足够的壁面附近网格密度,同时有效减少了网格总量需求。     

隐式无网格算法及其应用研究

本文的主要目的在于研究求解Euler方程的隐式无网格算法，并应用于复杂的流场计算。采用无网格算法，计算区域用点云离散代替通常的网格划分；计算点上的空间导数，用当地点云上引入的二次极小曲面逼近。求解的Euler方程用隐式时间后差离散，结合用Roe的近似Riemann解确定通量，并用LU－SGS算法分步计算，数值求解了单翼型或双翼型模拟的复杂绕流。     

吸气式空天飞行器闭环上升制导研究

闭环上升制导是吸气式空天飞行器达到快速优化、实时规划以及自主制导的必然要求。采用GHAME吸气式空天飞行器作为研究对象,提出了闭环上升制导最优控制方案,重点论证了该方案不存在关于节流阀值的奇异问题,并阐述了最优推力控制的Bang-Bang控制。引入有限差分法求解两点边值问题,并通过终端时刻调节算法求解终端时刻自由问题。仿真结果表明,所提出的闭环制导算法是有效、可行的,具有显著的优越性。     

共轴式双旋翼悬停流场和气动力的CFD计算

考虑到共轴式双旋翼流场特征高度复杂、桨叶承受非定常气动载荷的特点,为更好地预测共轴式双旋翼的气动特性,把非结构嵌套网格方法和网格的自适应技术相结合,发展了一套适合于共轴式双旋翼流场数值计算的求解器.在该求解器中,采用非结构嵌套网格方案来模拟桨叶之间存在的相对运动,自适应网格技术用来捕捉尾迹对流场和气动特性的影响,求解惯性坐标系下的非定常N-S主控方程来模拟流场的非定常特性.应用该求解器,首先计算了有试验结果可供对比的一副试验旋翼的诱导速度场分布,在此基础上,计算了共轴式双旋翼的桨尖涡轨迹和拉力分布特性,并与单旋翼的计算结果进行了对比和分析,得到了一些有意义的结论.     

用隐式高精度间断伽辽金方法模拟可压层流和湍流

为了提高隐式高阶间断伽辽金数值方法的稳定性,发展了一种基于解析精确Jacobian矩阵的GMRES隐式方法,用于求解可压缩层流和湍流问题。在GMRES的求解中,无黏通量和黏性通量的Jacobian矩阵采用链式法则解析精确求解,并用于线性系统方程的LU-SGS预处理和GMRES矩阵矢量生成;与此同时,对修正的负Spalart-Allmaras湍流模型的生成源项进行了修正,以避免隐式化求导时出现非物理解。通过典型层流和湍流算例对发展的方法进行了验证,研究结果表明:基于精确Jacobian矩阵的隐式GMRES方法,不仅能够提高隐式高精度间断伽辽金方法计算的稳定性,而且还能够提高计算效率。     

跨声速欧拉方程的一种隐式有限元解法

本文用有限体积Ｇａｌｅｒｋｉｎ法求解了跨声速二维欧拉方程。求通量运用了Ｏｓｈｅｒ格式，通过外插把原来只有一阶精度的逆风格式改进成二阶精度。并采用基于一性化流量的隐式格式和非结构化网格，使程序具有较高的收敛效率。     

一种求解透平叶栅三维流场的高精度TVD格式

本文提出了一种应用三阶精度Chakravarthy-OsherTVD格式求解环形透平叶栅三维定常流场的隐式近似因子分解方法。控制方程为柱座标速度分量表示的任意曲线座标系下的三维欧拉方程。采用Beam-Warming近似分解, 对隐式项进行迎风差分, 得到了时间精度为一阶、空间精度为5点三阶的隐式格式。算例对5点三阶和5点二阶精度的格式进行了比较。表明空间三阶精度隐格式对计算结果没有明显改进, 但收敛残差有较大降低。      

一种新的用于MHD方程逆风格式的Jacobian矩阵分裂方法

针对理想MHD方程,提出了一种新的基于MacCormack算法的雅可比矩阵分裂方法,克服了原有方法稳定性差的问题,并成功地应用于理想MHD方程的求解.控制方程在非结构混合网格上进行空间离散,其中对流项采用本文发展的逆风向量分裂格式,并引入了双曲型磁场散度清除技术,时间推进为显式5步龙格-库塔方法.对MHD激波管流动和带均匀磁场干扰的二维高超声速钝头体绕流流场进行了数值模拟,得到了与参考文献相吻合的数值结果,表明本文发展的数值分裂方法可以有效地捕捉MHD流场的流动特征,并且具有比MacCormack方法更高的稳定性和计算精度.     

捕捉间断的高精度数值方法

为发展适用于捕捉超声速流场中各种间断的高精度算法,将通量限制的思想引入到紧致格式中,构造了一个传统方法与紧致格式混合组成的通量限制型差分格式.通过在时间方向上利用一阶精度格式计算的一维定常激波,以及在时间方向采用多步Runge-Kutta方法计算的一维非定常激波管问题上的数值试验与二阶精度的TVD格式所计算的结果比较,表明新方法比二阶精度方法在间断的捕捉上具有明显的优势.通过新方法的计算结果与精确解的比较,表明新方法的准度也是非常令人满意的.     

叶轮机转子内流动的三维跨声Navier-Stokes方程组数值计算

应用有限体积法及Baldwin-Lomax代数紊流模型求解叶轮机转子内流动的三维可压缩非定常薄层Navier-Stokes方程组,计算一个绕X轴等转速旋转的转子内流场问题,用当地时间步长方法来加速收敛。采用隐式算法对对流通量进行求解,而扩散通量采用显式算法。采用高阶TVD下的Roe格式计算对流通量,采用中心差分格式计算扩散通量。关于守恒变量Q的非线性方程组,采用牛顿迭代方法求解,线化后的方程用对称的Gauss-Seidle方法求其渐近解。应用该算法模拟了转子NASARotor67流场,与试验做了大量的对比研究,包括性能、等马赫线以及总压和总温等等。在最高效率点附近计算与试验比较接近,但在靠近失速点时计算与试验略有差别。      

高精度差分求解气动方程的几个问题

本文探讨了发展高精度格式的必要性，研究了高精度格式与网格以及熵增条件的关系，并发展一种半离散化的空间为三阶精度的格式。从此半离散化的格式出发，可建立多步显式格式及隐式格式。模型问题的计算表明，该三阶精度的格式具有好的精度，且激波附近基本上没有波动。     

Godunov型显式大时间步长格式研究进展

综述了Godunov型显式大时间步长格式的研究进展。首先介绍了显式大时间步长格式的概念、分类和优势。然后重点阐述了Godunov型显式大时间步长格式的构造方法、高阶精度推广方法、多维问题推广方法和收敛特性、分辨率及计算效率等性能，展示了其在典型问题中的应用和验证。最后给出了Godunov型显式大时间步长格式研究进一步可能的发展方向。     

非均匀网格湍流大涡模拟高精度有限体积解法

为准确预测不可压复杂湍流,提出了一种可用于大涡模拟均匀或非均匀网格上的高精度有限体积法。该方法空间离散采用有限体四阶紧致格式,时间推进采用四阶Runge-Kutta法,压力-速度耦合应用四阶紧致格式的动量插值。通过直接求解顶盖驱动方腔流动证实了该方法具有近四阶的空间精度;并在此基础上,采用动态Smagor-insky亚格子应力模式,成功地实现了充分发展槽道湍流和后台阶湍流流动的大涡模拟计算,所得结果与直接数值模拟结果吻合良好,且采用非均匀网格可在比均匀网格数少的离散系统上得到同样满意的结果。结果表明,该方法是实现高精度湍流大涡数值模拟的一个有效途径。     

非结构网格上欧拉方程的自适应逆风隐式有限元解法

本文发展了一种解二维欧拉方程的隐式逆风有限元格式。空间近似应用的是二阶精度逆风方法；时间近似上根据数值通量函数的线性化处理导出了有限元的隐格式。为减少数值振荡，引入了基于特征变量的限制因子。在经自适应处理的非结构网络上，本文对跨、超、高超声速绕流进行了数值模拟。计算实践表明发展的隐式有限元格式与显格式相比具有较高的计算效率和较强的适应性。     

基于DES的二维和三维空腔流动特性研究

空腔流动广泛存在于航空航天工程中,对其流动特性的研究具有十分重要的工程意义。采用数值方法对比模拟了二维和三维空腔流动,控制方程采用N—S方程,空间离散采用有限体积方法,对流通量计算采用Roe格式,非定常时间离散采用双时间步长方法,湍流粘性计算采用基于SA模型的DES方法。数值计算所得的三维空腔底面压强分布与实验结果一致,所得的二维和三维空腔内的流动结构与相关文献中的分析相吻合。将二维和三维的计算结果进行比较,发现空腔底部压强分布、空腔内的压强振荡、空腔内部流线图及声压级等均有所不同,得出在空腔宽度有限时,横向流动简化了流动结构．削弱了振荡幅值．但使振荡过程复杂化．其三维效应不能忽略．