小波变换数值均匀化方法求解振荡系数的方程

将小波数值均匀化方法用于求解二维快速振荡系数的抛物型方程.小波数值均匀化方法基于小波的多分辨分析,在细尺度空间上得到原方程的半离散方程,然后利用小波映射得到粗尺度空间的数值均匀化方程,即可在粗尺度上求解原方程,有效地减小了计算量.计算结果表明:小波数值均匀化方法与精细剖分的有限差分法相比较,既大大地节省了计算时间又获得了较好的精度.     

某型涡扇发动机起动和加速过程数值模拟

借助国内外涡扇发动机过渡态数值模拟研究的最新成果,开展了某型涡扇发动机起动加速过程的数值模拟研究。首先参照风机和泵类机械低转速部件特性相似理论,结合指数外插法,成功地发展了航空燃气涡轮发动机低转速部件特性扩展方法和计算程序;其次,基于燃气涡轮发动机部件匹配原理,写出描述部件内动态流动过程的三个基本方程：动量守恒方程（功率平衡方程）、流量连续方程和能量守恒方程;最后发展了具有一定计算精度的涡扇发动机起动加速过程的数值计算模型和程序,对某型涡扇发动机的起动和加速过程进行了数值分析,并与试验数据进行了对比分析。结果表明,给出的数值计算结果与试验数据具有较好的一致性,表明发展的模拟涡扇发动机起动加速过程的数学模型和程序是合理的。     

CFD技术在动态压力校准中的应用研究

采用计算流体力学(CFD)流场可视化的手段,用数值模拟的方法研究了动态压力校准装置.研究工作主要集中于激波管动态压力校准装置和正弦动态压力校准装置的数值模拟,控制方程采用雷诺平均的N-S方程描述,方程的封闭采用精度相对较高的高雷诺数湍流模型,用有限体积法对控制方程进行空间离散和时间离散.并通过实验验证了控制方程和计算方法的正确性.     

电弧喷射推力器流动区域电磁场数值解法

为了有效地进行电弧喷射推力器流动区域全场(流场、电场、磁场和化学反应)耦合数值模拟并揭示其电磁场特征，建立了适当的电磁场模型，对电磁场控制方程的数值解法进行了详细研究。模型基于推力器稳态工作特征，控制方程由麦克斯韦方程组简化得到，采用有限控制容积积分方法离散椭圆型控制方程，采用9种不同的迭代方法求解离散方程。给出了不同迭代方法的收敛速度、数值稳定性和最终精度。研究表明，Gauss—Seidel逐线超松弛迭代法是推力器流动区域电磁场离散方程的一种快速有效的数值解法。     

基于控制理论和NS方程的气动设计方法研究

研究了基于控制理论和NS方程的气动设计方法,针对给定的目标函数表达形式,应用该设计理论在计算坐标下详细推导了相应的共轭方程及边界条件具体表达形式,以及梯度方程求解表达式,通过合理的数学变换,导出了共轭方程在笛卡尔坐标系下的直观表达形式,发展了有效的共轭方程数值求解方法,通过流动控制方程数值求解、共轭方程数值求解、目标函数对设计变量的梯度求解和优化算法等方面的有效结合,研究与发展了一种新的气动设计方法,以二维机翼气动设计为例,成功进行了亚、跨音速情形下的相关设计算例研究,研究结果表明应用控制理论和NS方程的气动设计方法在设计理论、适用性以及时间花费等方面都有着很好的特色和优点,且设计结果也更为可靠.     

二维舵面绕流的N-S方程数值计算研究

本文基于结构化网格,应用Baldwin-Lomax湍流模型进行绕舵面流动的二维N-S方程数值模拟,将部分计算结果与风洞试验结果进行了对比分析研究,分析了计算方法的可靠性和计算结果的合理性.为开展三维舵面绕流的N-S方程数值模拟和铰链力矩数值计算研究奠定了基础.     

低速叶型气动反问题设计方法

 低马赫数不可压流动中声速与流速大小差别巨大,采用基于可压缩流动控制方程的计算格式求解流场时,由于数值黏性的污染,解的精度低且收敛性差,通常可使用时间预处理技术来解决这一问题。在基于控制理论的优化方法中,共轭方程的Jacobian矩阵和流动方程的系数矩阵相似,因而在低流动马赫数下,求解共轭方程存在着与求解流动方程相同的数值污染和数值刚性问题。首先推导了带有预处理的Roe格式,然后发展了适合全速度流动的共轭方程求解方法,最后选取翼型和叶栅两个典型算例进行了验证。计算结果表明所发展的方法可很好地用于低马赫数时的气动反问题设计。     

平底圆筒底部垂直撞水的实验及数值计算

实验采用 4种不同厚度底板的圆筒模型。以不同下落高度撞水,测量底板径向压力分布和应变响应。数值计算中联立求解底板弹性动力响应方程和流场方程。实验和数值分析表明底板的初始几何缺陷和弹性变形对冲击压力具有较大影响     

二维引射—混合器流场的数值研究与验证

本文运用一维引射特性方程和不可压缩流动Navier-Stokes方程,对二维引射—混合流场进行了数值研究。首先用一维引射特性方程估算理想的引射气流流量;然后用N-S方程对混合管内的气流混合流动进行计算,用计算得到的混合管速度分布和壁面粘性剪切力来修正一维引射特性方程中的经验常数,逐次逼近真实的引射流量比。计算结果表明,本文的数值方法可以有效地预测引射流量比和混合流场。     

压电热弹性体开口壳的响应分析

根据广义的Hamilton变分原理推导出了压电热弹性体非齐次的Hamilton正则方程。结合压电热弹性体平衡方程和热平衡方程,成功地导出了压电热弹性体机、电、热耦合问题的齐次状态方程。将非齐次方程转化为齐次方程不仅使问题变得大为简化,同时也减少了数值计算的工作量。数值实例研究了温度载荷和力载荷作用下压电热弹性材料四边简支层合开口壳的响应问题。     

三维机翼的静气动弹性特性数值模拟研究

发展了一种计算流体动力学（CFD）和计算结构动力学（CSD）的耦合计算方法,对三维机翼的静气动弹性进行了数值模拟研究。采用三维欧拉方程为控制方程基于直角网格计算气动力,并耦合结构静平衡方程进行静气动弹性数值模拟,设计了CFD/CSD耦合计算的数据交换的方式。以M6机翼作为算例,进行了机翼静气动弹性的数值模拟,计算结果表明所发展的三维机翼静气动弹性数值模拟方法是合理可行的,并可作为机翼结构优化设计和考虑结构弹性变形影响的气动外形优化设计的基础。     

高Ra数自然对流问题的张力样条数值模拟

本文将基于压力泊松方程技术的原参数Ｎ－Ｓ方程的张力样条数值解法推广应用于自然对流问题的数值模拟，提出了自然对流问题一个新的数值模拟方法，对一方内高Ｒａ数自然对流问题的数值模拟表明，结果是令人满意的。     

考虑气动弹性影响的机翼气动外形设计研究

采用三维Euler方程为控制方程,耦合静气动弹性平衡方程,进行机翼静气动弹性数值模拟;在机翼静气动弹性分析的基础上,结合Takanashi余量修正方法对三维大展弦比机翼进行气动外形反设计,以确定机翼的型架外形.以某型支线飞机的大展弦比机翼为算例,进行了静气动弹性数值模拟和机翼型架外形设计研究,设计结果表明发展的机翼静气动弹性数值模拟和型架外形设计方法是合理可行的.     

样条Lax-Wendroff格式和样条跳蛙格式

1.引言 用三次样条求解粘性流体力学问题的数值方法已有了相当的发展。为了探讨用该方法求解无粘流动问题的可能性,本文尝试把差分方法中一些典型的计算格式扩展到样条方法,提出了样条Lax-Wendroff格式和样条跳蛙格式。以一维模型方程为例,按照Warming等的方法导出了它们的修正方程,分析了格式的稳定条件及数值弥散和耗散效应。然后扩展到对一般非线性双曲方程求解的研究。最后,讨论了某些数值解例。     

Euler及N-S方程的二维四边形网格自适应算法

描述了一种基于贴体网格的Euler和N-S方程自适应算法。采用中心有限体积法,对Euler及N-S方程进行数值求解,N-S方程的计算中加入了基于SA湍流模型的DES模型。在流场中,运用了网格自适应算法,提高了数值计算对激波、流动分离、涡等特性的捕捉和分辨能力。采用上述方法,数值分析了单段翼型的绕流问题。文末给出的2个二维问题的对比算例表明本文方法有效。     

计算翼型跨音速有势绕流的压强极小积分有限元法

本文采用压强极小积分有限元法计算了NACA 0012翼型的跨音速流场。出发方程为全速位方程。用人工密度的方式引入人工粘性。计算网格是通过解析/数值变换的方法自动生成的。采用线松弛法和人工时间相关法求解有限元方程。给出了有关的数值结果。     

不等熵一维定常管流数值解法

不等熵一维定常管流控制方程的数值解法,已有许多讨论,但大多局限于某些特定情况,特别是控制方程在奇点(M=1)附近的计算,还没有一个简便的方法.本文引进新变量,把方程变为非奇异的,提出了对各种不等熵流都适用的数值解法.由于引进了喉部条件,保证了数值解的唯一性.经过对模型方程的计算精度分析,证明该数值解除了在奇点附近稍有误差外,其它地方与精确解一致.本文最后给出了该方法在无喷管发动机内弹道计算中的应用.     

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

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

基于黏性伴随方法的跨声速机翼气动优化设计

 进行了基于黏性伴随方法和Navier Stokes方程的跨声速机翼气动优化设计研究。分别推导了适用于三维跨声速机翼气动反设计和减阻设计的黏性伴随方程、边界条件和梯度求解表达式,并研究了伴随方程的数值求解方法。通过将网格生成、流场计算、黏性伴随方程数值求解、梯度求解和拟牛顿优化算法等几方面的有效结合,发展了一种跨声速机翼气动优化设计方法。为了提高计算效率,将多重网格方法应用到方程的数值求解中来加速收敛。跨声速机翼反设计和减阻设计算例验证了本文所发展的方法的正确性。采用本文的方法进行优化设计,一般通过20~30次迭代就能得到满意的结果。     

昆虫悬停飞行的动稳定性:理论分析与数值模拟

首先,从昆虫(身体和拍动翅)的动力学方程和N-S方程出发,在一定假设下,将运动方程简化为6自由度刚体的方程,并用线化理论给出了反映昆虫悬停飞行纵向动稳定性的理论解.然后,用完全的动力学方程和N-S方程偶合的数值解对上述理论的简化假设进行了检验.