一种求解任意形状导体目标内谐振时表面电流的有效方法

用矩量法求解闭合导体目标的电场积分方程或磁场积分方程将得到不正确的表面电流。提出了一种确定正确表面电流的有效方法，即正确的表面电流由非谐振模电流和归一化的谐振模电流与一未知复因子的乘积组成，通过令闭合导体内给定点的总电场为零可得该未知复因子。计算了一无限长理想导体圆柱谐振时的表面电流，所得结果与解析解一致，进而证明了本方法的有效性和准确性。     

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

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

改进节点积分的无单元Galerkin法及其在流动问题中的应用

无单元Galerkin法需要在背景网格上积分,计算量大,且在求解对流占优问题时会出现非物理的数值伪振荡现象。为此,基于局部Taylor展开思想,采用节点处的局部Taylor展开计算积分,建立了局部Taylor展开积分无单元Galerkin法。该方法同时解决了标准的无单元Galerkin法计算量大和对流占优时会出现数值伪振荡的问题。一维定常对流扩散方程和二维Burgers方程的求解说明了该方法的有效性,且计算效率远高于无单元Galerkin法。     

奇异积分方程方法在接触压力分析中的应用

总结回顾了奇异积分方程的数值解法,对于第二类奇异积分方程使用分段连续函数法进行了求解;对于两个弹性体接触问题,通过接触体之间的滑移函数和间隙函数,建立求解接触压力的奇异积分方程.分别针对圆柱体与弹性半空间体、抛物线型压头与弹性半空间体、榫头与榫槽3类接触问题,确立奇异积分方程的具体表达式,而后使用分段连续函数方法进行求解,获得接触面上的接触压力.最后将计算所得的接触压力分别与理论解和有限元解进行了对比.对于圆柱体与弹性半空间体接触问题,奇异积分方程法的最大接触压力与理论解和有限元解的相对误差分别为0.3%和0.5%;对于榫头与榫槽接触问题,奇异积分方程方法计算所得的最大接触压力与有限元解的相对误差为1.8%,验证了奇异积分方程方法的有效性.      

弯曲网格上的间断有限元湍流数值解法研究

采用间断有限元方法对雷诺平均Navier-Stokes(RANS)方程进行了数值求解,对Spalart-Allmaras单方程湍流模型进行了部分修正,使得求解器更加鲁棒。构造了分段高次多项式来逼近真实物面,同时物面附近采用多层弯曲网格来避免网格交叉,此外提出了一种快速计算积分点的曲面物面距的方法。采用混合网格对NACA0012翼型以及RAE翼型进行了数值计算,并与实验数据以及前人数据进行了对比。计算结果表明,采用物面弯曲网格结合修正的湍流模型方法在相对稀疏的网格上就能得到比较好的数值解。     

可压缩紊流分离边界层的计算

本文介绍了一种求解二维可压缩紊流分离边界层的正反积分方法。该方法以边界层动量积分方程和平均流动能积分方程作为基本方程,采用适用于贴附和分离紊流的Swafford速度型解析式,对比计算了四种代数涡流粘性模型对沿流向各点的耗散积分及诸边界层特征参数的影响。计算结果和存在激波边界层较强干扰的跨音速翼型边界层的测量数据吻合。该方法可扩展用于计算翼型的尾迹特性。     

非紧致气动噪声半空间传播基本解的边界积分方法

为了获得低马赫数流动诱发的非紧致气动噪声在半空间内传播的基本解,结合复等效源方法和边界元方法建立了半空间精确格林函数的边界积分方程,当半空间边界为刚度型阻抗边界时可避免奇异性积分。同时基于等效源方法提出一种半空间二维非紧致圆柱声散射模型,推导了静止介质中声散射基本解的理论表达式。对静止介质中的二维圆柱声散射,数值解在研究的频率与观察点处与理论解一致。采用数值方法计算二维圆柱绕流诱发的半空间声场基本解,结果显示半空间边界强化了声源的声辐射能力,小于马赫数0.2的介质运动对声传播的影响可以忽略。     

可渗透面对流FW-H方程伪声产生机制和抑制方法

可渗透面对流FW-H方程及其积分解在远场气动噪声预测中广泛应用，但当涡波通过可渗透面时会产生伪声传播。本文旨在清晰阐明这种伪声的产生机制，并提出有效抑制方法。耦合分析可渗透面对流FW-H方程的右端厚度源和载荷源，发现厚度源的物质导数和载荷源的散度操作可使积分面源自动过滤伪声源。伪声的产生源于对流FWH方程的求解采用了Farassat提出的分部积分公式，因部分积分项被忽略，导致伪声源自动过滤功能失效。从厚度源和载荷源中抽取含有涡波扰动的项，在对流波动方程求解过程中保留相应的物质导数和散度操作，以抑制涡波伪声传播。数值测试算例验证了伪声产生机制和抑制方法的正确性。     

全速位积分方程跨音速绕流的数值解

从全速位方程出发,利用Green公式将其化为激波捕获积分方程和激波装配积分方程,然后离散进行数值解。流场出现激波时,对激波捕获积分方程应用上风技术捕捉到激波,然后应用激波装配技术计算,得到了满意的结果。经算例考核,该方法具有计算区域小,收敛快和CPU时间较少等优点。     

准则方程形式的探讨

应用积分方程求解换热问题,不仅数学运算较简单而且计算结果具有很好的直观性,容易看出理论解的物理意义,其计算结果使人们一目了然。下面我们利用积分方程,在第一类边界条件求解Pr=0.001～7500的流体绕流平壁层流换热问题。      

一种无条件稳定的时域推进法与有限元混合分析瞬态热传导问题

本文将一种新的无条件稳定的时域推进法与有限元结合,用于分析瞬态热传导问题。根据微-积分型热传导控制方程,对时间变量在小区闻内插值积分,变初-边值问题为一系列离散时刻的边界值问题,再应用有限元法求解之。由于推进中的每一时刻的解,都严格满足原给定的初始条件,就消除了累积误差的影响。此外,这种时域推进法对时间变量的数值积分也优于数值微分。因而,与现行的基于微分型控制方程,采用时-空有限元模型同时离散,再逐步求解的直接积分算法相比,计算精度可大大提高。不难预料,当求解较长时间后的瞬态值时,本算法的优越性会更加明显。     

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

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

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

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

月球探测器发动机羽流辐射热流计算

火箭发动机羽流的辐射的计算需要求解气体介质的辐射传递方程(RTE).计算模型中基于单线组(SLG)模型考虑气体组分的光谱特性,并采用离散传递法(DTM)求解RTE,对波长积分得到辐射热流密度.通过对圆柱高温CO₂气团辐射热流计算来验证计算程序,然后在流场计算的基础上,计算了月球探测器主发动机和姿控发动机羽流对各部位和仪器的辐射热流密度,根据对不同工况计算结果的分析为探测器的设计提供参考.      

导弹运动方程的逆解析解

针对在采用反馈线性化理论设计导弹飞行控制系统时控制器算法对导弹逆模型的要求 ,推导了导弹运动方程的逆解析解 ,包括导弹质心运动的运动学方程及其动力学方程的逆解析解、姿态运动的运动学方程及其动力学方程的逆解析解 ,详细给出了逆模型的求解过程 ,为一类飞行控制系统运用全局线性化设计奠定了基础。最后 ,结合工程应用进行了仿真计算 ,仿真结果表明逆解析解的推导是正确的     

超燃冲压发动机流量匹配的临界流量法

为了使用基于积分方程的超燃冲压发动机性能计算模型模拟跨燃、亚燃工作模态,用临界流量法建立用于表征燃烧室热壅塞程度的残差方程,并以隔离段出口静压为独立变量,通过牛顿迭代法迭代求解残差方程,模拟了燃烧室内存在临界截面（即喉道）但不出现热壅塞的实际跨燃、亚燃工作状态。结果表明：该方法能够确定有喉道的超燃冲压发动机燃烧室一维参数分布,并具有计算精度高、计算速度快、收敛性好等优点,可以为基于积分方程的超燃冲压发动机性能计算模型求解提供一定的参考。     

计算机翼跨声速绕流的全速位积分方程数值解

本文利用格林定理将全速位方程转换为全速位积分方程，应用高斯定理对沿激波面的显式积分进行简化。将全速位积分方程的解表示成为面元项和场元项之和，然后离散进行数值求解。本文将激波捕捉方法和激波装配技术结合起来，对机翼跨声速绕流进行数值实验，计算结果令人满意。     

轴对称喷管内外流场与结构温度场耦合计算

利用有限体积的对流项二次迎风插值格式和重整化群(RNG)k-ε湍流模型,二层增强型壁面函数,同时利用球形谐波法考虑热辐射的影响,以灰气体加权模型(WSGGM)确定气体介质的辐射性质,求解N-S方程、热传导方程、考虑吸收-发射性气体介质的辐射传输方程。采用流固耦合的流动与换热模型,流场与结构温度场互为边界条件交换数据,实现了流场解算与温度场解算的耦合数值分析。采用此计算模型对静止的轴对称收-扩喷管进行了数值模拟,计算得到的喷管壁温与试验数据吻合良好。在此基础上,进一步开展了巡航状态下轴对称收-扩喷管内外流场与结构温度场的耦合数值研究。     

飞机稳态尾旋平衡点的综合求解方法研究

为求解飞机稳态尾旋平衡点,推导了飞机稳态尾旋的动力学方程,并针对该方程给出了三种求解方法:作图法、数值求解法以及在前两种方法的基础上所提出的一种综合求解方法。综合求解方法首先使用作图法得出方程的所有近似解,并用稳态尾旋判据判断解是否为稳定解,然后将稳定近似解作为数值求解法的初值进行精确求解。运用综合求解方法对算例飞机进行了稳态尾旋平衡点的计算,算法收敛迅速,结果精确,表明该方法具有很好的工程应用价值。     

改进虚拟流体方法在气-水界面捕捉中的应用

提出了一种改进的Ghost方法,通过求解双介质的Riemann问题,将全部解参量外推得到虚拟流体状态参量,使单一的气-水间断问题转化为两个单介质的Riemann问题。流场控制方程求解采用二阶NND格式,LevelSet方程用五阶WENO格式离散,实现了对气-水界面稳定而准确的捕捉。激波与气-水界面作用的一维和二维算例验证了程序的有效性。对气液两相均运用可压模型,计算了单个水滴被马赫数为2的空气激波扫过后的二维变形过程。