时域推进与边界元混合瞬态热传导问题

对时间变量在小区间内应用多项式插值,对空间变量采用边界元求解,提出一种新的分析瞬态热传导问题的计算格式。与现行的利用含时间变量的基本解所得到的边界元计算格式相比较,具有基本解简单易得、系数矩阵容易计算以及无重复积分等优点;与向前差分消去时间因子的边界元计算格式相比较,则在算法稳定性、计算精度以及时间步长的选取等方面具有明显的优势,并且很容易推广到求解非线性问题。算例表明,采用时域推进与边界元混合方法分析瞬态热传导问题有效可行。     

考虑辐射效应的多孔材料传热有限元算法

应用有限元法求解一维半透明多孔介质辐射对流传热问题。通过对辐射传递方程与边界条件中积分项的离散 ,推导出有限元求解列式。对于多孔介质传热控制方程和气体传热控制方程的求解 ,采用时间积分技术对两相温度统一进行求解。而对于边界方程处理 ,则通过传热方程与边界方程的循环求解 ,求出任意时刻温度场的瞬态解。给出数值算例 ,得出了放置在通风口、同时受到伴随辐射作用的半透明多孔介质传热的瞬态解 ,讨论了部分参数对瞬态温度场以及换热效果的影响      

求解非定常N-S方程的一种自适应时间步长的罚函数有限元方法

 本文提出了一种求解速度-压力变量形式的非定常不可压N-S方程的数值方法。其中采用罚函数的Galerkin有限元进行空间离散,改进了梯形法的时间积分格式并结合自适应时间步长技术使计算效率最优。文中对圆柱突然启动到Re=100和200时的流场进行了数值模拟,所得到的Karman涡街的诱发过程是满意的。     

化学非平衡解耦算法精度分析及其改进

模拟化学非平衡流问题涉及到流动方程和化学反应方程两部分的求解,流动方程组中包含所有组元的偏微分方程,导致求解变量成数量级增加,点隐算法和全隐算法在求解源项刚性问题时,通常会涉及到矩阵求逆运算,这些因素带来的巨大计算量限制了隐式算法在复杂工程问题中的应用。1993年刘君提出了化学非平衡流解耦算法,将控制方程组分解为流动和化学反应两部分,流动方程的求解采用冻结流假设来描述流体微团沿流线的运动过程,化学反应方程的求解描述流体微团在随体坐标系下发生绝热、定容的爆炸过程。结合解耦算法和有限体积法的特点对这种解耦算法进行改进,提出的优化算法不需要求解组元变量所对应的偏微分方程组,只求解由5个基本变量构成的,形式上与量热完全气体近似的偏微分方程组,通过对比计算结果发现优化算法可以显著地提高计算效率。同时,将精细积分方法应用于化学非平衡流问题的求解中,通过与传统的VODE方法和α-QSS方法对比发现,精细积分方法的鲁棒性更优、精度对时间步长不敏感,适当的选取时间步长,可以充分发挥精细积分方法的优势。     

基于离散伴随方程的三维雷达散射截面几何敏感度计算

针对有限差分法计算雷达散射截面（RCS）梯度效率低，采用高精度雷达散射截面评估时计算代价高的问题，提出了一种基于麦克斯韦积分方程离散伴随方程的RCS梯度高效计算方法。基于伴随方程的梯度计算可以通过一次雷达散射截面求解、一次伴随方程求解获得RCS关于所有设计变量的梯度。其中麦克斯韦积分方程离散伴随方程的形式与原方程基本一致，可以采用矩量法（MOM）及多层快速多极子算法（MLFMA）求解。伴随方程求解计算量与直接雷达散射截面评估基本一致，存储量在直接雷达散射截面评估的基础上增加不明显。通过双椎体模型、导弹模型对基于矩量法、多层快速多极子算法的伴随梯度进行验证，证明了基于伴随方法的RCS梯度计算可以实现复杂外形中RCS梯度的高效、高精度求解，为基于梯度的高精度气动/隐身一体化优化提供了基础。     

改进的电热除冰系统仿真

建立了电热除冰系统矩形微段的二维简化物理模型,采用焓法模型建立相变导热问题的微分控制方程,用控制容积法对该模型的控制方程及其边界、初始条件进行离散,使用块修正技术对原有的离散方程组的求解方法进行了改进,得到各个时刻的温度分布及冰的融解情况.计算结果与文献实验结果吻合,验证了改进的算法和程序的正确性.并计算分析了时间步长的选择、网格的划分对计算结果的影响.      

改进的电热除冰系统的仿真方法

建立了电热除冰系统矩形微段的二维简化物理模型，采用焓法模型建立相变导热问题的微分控制方程，用控制容积法对该模型的控制方程及其边界、初始条件进行离散，使用块修正技术对原有的离散方程组的求解方法进行了改进，得到各个时刻的温度分布及冰的融解情况。计算结果与国外文献结果吻合，验证了改进的算法和程序的正确性。并计算分析了时间步长的选择、网格的划分对计算结果的影响。     

金属蜂窝夹芯面板有效导热系数的数值计算

对具有漫灰体内壁的金属蜂窝夹芯面板,忽略蜂窝内腔的空气导热,综合考虑蜂窝结构的热传导和热辐射等热传递过程,利用有限元方法求解了周期性分布的蜂窝单胞稳态热传导控制方程,蜂窝内壁的边界条件是由净热量法得到的热辐射换热积分方程,由胞元的温度场分布数据及Fourier定律得到了蜂窝结构的有效导热系数。与现有文献相比,采用了较少近似的模型及较高精度的离散方法,计算结果表明计算模型及方法是可靠有效的。     

基于多物理场耦合的双脉冲发动机点火过程数值模拟

为研究双脉冲固体火箭发动机Ⅱ脉冲点火瞬态过程,发展一套多物理场耦合求解器。流体控制方程基于有限体积法求解,时间推进采用双时间步LU-SGS(Lower Upper Symmetric Guass Seidel)方法;固体推进剂表面温度基于耦合传热方法计算;结构动力学运动方程基于有限元方法离散,采用经典的Newmark格式进行时间推进,流固耦合采用松耦合算法,并通过算例验证求解器的可靠性。计算结果表明:该求解器能够数值模拟Ⅱ脉冲启动过程中的点火药气体冲击、燃气非定常流动及金属膜片机械响应过程,获得金属膜片的破裂时间和压强;且随着点火质量流率增加,推进剂装药首次点燃时间和金属膜片破裂时间变短,膜片破裂压强降低;金属膜片破裂时间和压强不仅与作用在其表面的压强载荷大小相关,而且与压强载荷加载的过程相关;金属膜片厚度越薄,膜片破裂时间越短,膜片轴向位移越大,膜片破裂压强越低。     

Navier-Stokes方程两层稳定有限元算法分析

分析了定常Navier-Stokes方程的两种两层稳定有限元算法。它们将局部Gauss积分稳定化技术和两层算法的思想充分结合,采用低次等阶有限元P1-P1或Q1-Q1对N-S方程进行数值求解。误差分析和数值算例结果表明,当粗、细网格尺度H=O（h1/2）时,它们与在细网格上的单层有限元算法具有相同的收敛速度,而两层算法却节省了大量的计算时间。相比之下,Simple算法具有更高的计算效率。而且进一步发现Oseen算法能够对小粘性系数N-S方程进行有效求解。     

N-S方程有限分析数值研究

 本文将有限分析方法用于曲线座标系上紊流N-S方程的数值计算。分别计算了单列和串列叶栅内部流场,计算中采用了k-ε紊流模型和壁面函数。计算结果与试验结果相比较,吻合程度令人满意。有限分析方法在网格单元上对N-S方程进行线化处理,以解析边界条件作为约束,得出解析解,在解析解基础上构造离散代数方程。有限分析方法的最大特点是可以适应对流速度大小和方向,自动调整格式系数,因而具有数值扩散小和稳定性高等优点。     

机翼和翼身融合体跨音速绕流有限元算法的研究

从全速位方程出发研究了机翼和翼身融合体跨音速绕流的有限元算法。给出适用于有限元法和复杂外形的块结构网格生成技术和使单元系数矩阵积分数目减少2/3的方法;并发展了一种内存量少的总体系数矩阵存贮方式。为了提高有限元离散方程组解的精度和加速收敛,提出α—GMRES算法代替传统全速位方程有限元解法中所习用收敛慢的线松弛迭代解法。通过机翼和翼身融合体实例考核,效果是好的,为飞行器先进气动布局的数值计算提供了一种新的有效算法。     

单列叶栅紊流流场有限分析数值研究

本文将有限分析方法用于曲线座标系上紊流 N-S方程数值计算,研究了高雷诺数时叶栅粘性紊流流场。有限分析方法在网格单元上对非线性偏微分程进行线化处理,在解析边界条件下求出线化偏微分方程的解析解,以此解构造离散方程。有限分析方法能够根据对流的方向和大小自动改变格式系数,并具有数值扩散小、精度高和稳定性好等优点。本文以 k-ε紊流模型模化紊流,以壁面函数方法处理近壁区流动参数。      

各向异性Chaboche黏塑性本构方程隐式应力积分算法

针对各向异性Chaboche黏塑性本构方程,利用有效参数控制法推导了该方程的隐式积分算法,该算法把所有未知量都用同一个控制参数表示,并形成控制方程,从而把本构方程应力积分的问题转化成一个非线性控制方程的求解问题.该积分算法程序,与通用有限元软件Marc结合,通过单个单元对单向拉伸、循环、松弛的响应的实验数据进行了数值模拟,通过与显式算法对比,该隐式算法的稳定性和收敛性有所提高.构件的对比分析也表明隐式算法比显式算法具有较高的效率.      

梁/液体耦合振动的边界积分对称列式法

 利用格林积分导出的液体内部扰动压力函数作为测试函数,通过变分得到一组具有对称系数的代数方程,将该方程同梁的有限元振动方程联立后,即可求解梁在液体中的耦合振动问题。本方法具有未知数目少,数据存储和计算效率高的特点。     

跟踪问题最优控制律精细积分

构成有限时间最优跟踪系统的控制律需要求解 Riccati微分方程及外部控制输入向量满足的微分方程 ,前者是非线性矩阵微分方程 ,后者是变系数线性微分方程。在结构力学与最优控制的模拟理论基础上所发展的精细积分方法借鉴了计算结构力学中的算法 ,可以精确有效地求解这些微分方程。这种方法的特点之一在于步长幅度变化较大时 ,Riccati微分方程的数值解仍可以保持很高的精度 ,并且变系数线性微分方程的求解亦可纳入其体系而不必用通常的差分方法。本文介绍了用精细积分方法求解这些方程的过程 ,并给出了数值算例。     

线性GMRES算法在求解全速位方程中的应用

从全速位方程出发,利用压强极小积分得到有限元方程组。采用适用于有限元法的块结构网格生成技术,给出求解有限元总体方程组的线性GMRES算法。该算法比线松弛迭代法的收敛速度快得多。对M6机翼跨音速有升力情况进行了计算,其结果与实验值符合较好。     

轮盘振动特性计算

本文提出了一种计算轮盘振动特性的数值解法。采用薄板和厚板两种理论,将问题化成一阶微分方程组进行数值计算。编制了计算机程序,得出了令人满意的结果。由于没有涉及到复杂的数学问题,计算机程序大为简化,计算时间少。      

求解平面翼型亚,跨声速绕流的一个新方法

本文提出了求解平面翼型亚、跨声速绕流的一个新方法。引入流函数和Ｖｏｎ Ｍｉｓｅｓ变换后，亚、跨声速绕机翼无旋流动的基本方程组被化为以流线纵坐标ｙ为未知量的单个二阶偏微分方程－流线控制方程。并通过变换将物理平面上的无限域变为计算平面上有限的矩形域，而后在计算平面采用有限差分线松弛迭代法求解。作为算例，计算了对称翼型ＮＡＣＡ００１２－３４和非对称翼型ＮＡＣＡ４４１２的亚、跨声速有攻角绕流，所得数值结果