导热-辐射耦合传热的多尺度分析和数值模型

复合材料高温传热特性的准确预测对飞行器热防护结构的设计有重要意义,相关模型也是国家数值风洞工程中多相多介质计算模型的重要组成部分。针对周期性结构复合材料的高温传热问题,利用多尺度渐进分析方法,对包含导热方程和辐射传输方程的导热-辐射耦合传热模型开展了研究。建立了表征单元尺度模型及宏观平均导热-辐射耦合传热模型,获得了材料宏观等效导热系数与表征单元模型间的关系,发现宏观等效辐射吸收和散射等系数可通过表征单元内的体积平均求取。结合有限容积方法与格子Boltzmann方法,建立了复合材料导热-辐射耦合传热多尺度数值模型。采用二维常物性材料传热过程的模拟验证了多尺度模型的有效性,结果表明所建立的多尺度模型能够对温度场给出准确高效的计算结果。该方法有助于为复合材料传热特性的预测提供数值手段。     

不规则地形条件下双向DMFT电波传播特性算法研究

 在传统中心差分离散混合傅里叶变换(DMFT)抛物方程模型基础上,针对传统算法计算动态阻抗边界条件时的数值振荡及忽略后向传播与散射造成计算误差增大的问题,提出基于前-后向混合差分方法的递归双向DMFT(TW-DMFT)模型,以一阶前-后向差分方程代替传统的中心差分方程拟合阻抗边界条件,以递归双向传播模型改进传统单向预测模型.仿真分析了双向DMFT模型计算不同边界及媒介条件下单刃峰和实际地形的电波传播分布特性,与传统算法及几何绕射理论(GTD)进行比较,新算法提高稳定性及精度的同时弥补了传统算法的不足,更具有普适性.     

航空发动机故障诊断方程的求解及应用

用换算系数方法建立了基线模型，分析了测量参数偏差的计算方法，重点介绍了发动机故障诊断系统中故障方程的求解方法，并应用发动机试车数据对此进行了求解检验，结果表明故障诊断方程的求解方法可行。     

基于MAP估计和广义高斯分布的SAR图像边缘比率检测方法

 合成孔径雷达(SAR)图像的低信噪比和乘性相干斑噪声给SAR图像的边缘检测带来了极大的困难.通过引入广义高斯(GG)分布作为局部均值功率的先验分布模型,给出了局部均值功率在最大后验概率(MAP)意义下的最优估计,进而提出一种新的SAR图像边缘比率检测算子.利用以梅林变换为基础的对数累积量(MoLC)方法估计GG分布的参数,在此基础上给出一种局部均值功率MAP估计和GG分布参数估计的联合迭代求解方法.利用SAR实测数据对本文提出的边缘检测算子进行仿真验证,并将其与平均比率(RoA)算子和指数加权均值比(ROEWA)算子进行了对比,结果表明该算子可以有效克服相干斑噪声的影响,边缘定位准确且虚假边缘明显减少.     

某型涡轴发动机等压差活门建模分析

考虑某型涡轴发动机燃油流量调节器通道截面的局部能量损失,推导出流经计量油针的燃油流量方程,在此基础上详细分析等压差活门的功用、结构组成及其工作原理,以连续方程和力平衡方程为基础,采用线性化处理方法,建立等压差活门的数学模型。通过分析等压差活门系统结构,得到等压差活门系统的传递函数,进而对其稳定性进行分析,计算单位阶跃输入的稳态误差。结果表明:以流量方程和力平衡方程为基础建立等压差活门数学模型的方法可行,分析出等压差活门稳定工作条件和稳态误差的主要影响因素,并且为燃油调节器数学模型的建立以及仿真计算奠定了基础。     

基于二维厚板模型的波箔片轴承静特性

建立平箔片的二维厚板有限元模型,运用有限单元法和有限差分法耦合求解Reynolds方程和气膜厚度方程,研究了在两个工作转速下气体波箔片轴承在中截面和边缘处最小气膜厚度随轴承承载力变化规律.通过数值仿真对该模型、一维梁模型、二维薄壳模型和文献实验数据进行对比分析,结果表明:在轴承中截面处,3个模型的最小气膜厚度仿真结果都与实验结果符合得很好,但在轴承边缘处,由于二维厚板模型考虑了平箔片的剪切效应,因此其最小气膜厚度比二维薄壳模型的结果更接近实验值,而一维梁模型只考虑轴承圆周方向,因此不能体现气膜厚度沿轴承长度方向的变化规律.通过研究,为分析箔片轴承动力学特性奠定了理论基础.      

一种MEMS陀螺仪零偏建模与估计方法

针对MEMS陀螺仪零偏随时间变化的问题,提出了一种MEMS陀螺仪的零偏建模与估计方法.通过分析静态条件下陀螺仪零偏变化的影响因素,建立了陀螺仪零偏和温度之间的差分方程模型.采用系统辨识的方法,给出了系统模型的参数辨识流程.基于最小二乘法,实现了模型参数的辨识.设计了测试验证试验,基于陀螺仪的零偏样本特性建立了三阶差分方...     

TRIP2．0软件湍流模块的开发与确认

基于TRIP2.0_SOLVER数值模拟软件,开发了湍流模型方法库,目前该CFD软件中集成的湍流模型主要包括B-L代数湍流模型、SA一方程模型和SST两方程模型.为了考核不同的湍流模型工程适用性,本文采用对接网格技术,通过有限体积法数值离散三维任意坐标系下的RANS方程组,应用LU分解、MUSCL差分格式和低雷诺数SA和SST两种湍流模型,数值模拟了二维NLR-7301两段翼型和三维DLR-F6翼身组合体的绕流流场,计算与试验比较的内容包括了表面的压力系数分布、典型站位的速度型和气动特性曲线等内容.通过计算与试验的比较,在本文的计算范围内,采用两种湍流模型均可以得到与试验结果相吻合的压力分布和升力曲线,但在边界层内的速度型、粘性阻力和力矩特性等方面,不同的湍流模型具有明显的差异.     

相变调温服相变传热研究(英文)

为了研究相变调温服的热防护性能,介绍了人体热调节模型,利用焓法建立了相变传热数学模型。在整个区域建立统一的能量方程,采用控制容积有限差分法对方程进行隐式离散,用追赶法求解方程组得到各节点的焓值,求出相变材料内表面的温度。通过相变调温服人体热调节模型,计算出人体的动态温度分布及出汗量。计算结果表明,核心温度的变化与试验结果相符,出汗量相差不大,所建相变传热模型合理可靠。通过比较人体穿着相变调温服和未穿该服装时的动态温度分布及出汗量,可以看到穿着相变调温服能显著地减轻人体的热负荷。     

Burgers方程高精度差分格式分析

以非线性扩散波动方程Burgers方程为模型,通过调节方程右端粘性项的大小,在分析了方程的性质的基础上.探讨方程左端项和右端项差分离散格式匹配问题,以及高精度格式稳定性影响和计算格式精度与粘性系数和计算网格之间的关系等问题.研究结果将对方程左右端项离散时应如何采取合适的差分格式有一定的指导意义.     

圆筒燃烧室紊流流场数值研究

计算方法 本方法不需解压力修正方程,而是通过速度散度余量修正迭代法使连续方程直接得到满足。紊流模型采用k-ε模型,壁面处理采用无滑流固壁条件,避免复杂的壁面函数计算。对燃烧室内火焰稳定器,无论是圆盘或锥形都采用简单的取除格子法,从而大大简化计算程序。同样算例本方法所用计算语句只是SIMPLE算法的1/4～1/3。      

三维高超声速喷流干扰流场的数值模拟

采用ＬＵ隐式方法求解了有限体积法离散的完全Ｎ－Ｓ方程,数值模拟了三维高超声速喷流干扰流场。计算中采用了Ｊａｍｅｓｏｎ的当地极值递减（ＬＥＤ）格式和Ｂａｌｄｗｉｎ－Ｌｏｍａｘ代数湍流模型,与实验结果相比较,数值方法对附体流动计算较准确,对分离区附近流场计算有一定误差。     

大迎角非定常气动力建模方法研究

以战斗机俯仰机动为例,建立了飞机大迎角非定常气动力模型,包括非线性代数模型、Fourier函数分析模型、状态空间模型、差分方程模型以及模糊逻辑模型.并用SDM标模大振幅俯仰振荡非定常气动力风洞试验数据对模型参数进行了辨识,验证了模型的有效性.从模型物理意义、参数辨识难易程度及模型通用性和精确度等方面对几种非定常气动力模型进行了比较研究.     

箔片结构库伦摩擦效应对径向箔片轴承特性的影响

基于有限元法建立波箔型气体径向轴承箔片结构的库伦摩擦模型,通过改变平箔片与波纹箔片之间以及波纹箔片与轴承座之间的摩擦因数,对比分析了在各种载荷分布条件下波纹箔片库伦摩擦模型与文献中线性弹簧模型的刚度特性,研究了库伦摩擦效应对波纹箔片刚度特性的影响规律.在此基础上,运用有限单元法和有限差分法求解雷诺方程和气膜厚度方程,研究了在两个工作转速下气体波箔片轴承中截面处最小气膜厚度随轴承承载力的变化规律以及承载力随偏心率的变化规律.通过数值仿真对该模型、文献中线性弹簧模型和刚性表面气体轴承进行对比分析,并把气膜厚度分布与文献结果进行了对比.结果表明:箔片的库伦摩擦力在一定程度上增大了波纹箔片的刚度,并且随着摩擦因数的增大其刚度以及两端固定的波纹箔片个数也增加,使得箔片轴承表面变“刚”,因此轴承静特性更趋于刚性表面轴承,此外当轴承承载力一定时,箔片摩擦因数越大轴承的最小气膜厚度越小.      

飞机大攻角空间机动气动力建模研究

 准确的气动力模型是大攻角机动仿真与控制等飞行动力学问题研究的基础。针对纵向运动建立的非线性非定常气动力微分方程模型扩展应用于一般空间运动形式,在分析大攻角非定常流动机理和试验数据的基础上,发展了一种多自由度非线性非定常气动力数学模型,并探讨了模型辨识和参数估计问题。用全机构形的大攻角动态风洞试验数据对模型进行的验证表明所发展的气动力模型能够描述大攻角气动力主要的非线性和非定常特性。     

时间推进积分方程法求解非定常跨声速流动

本文提出了一种求解非定常跨声速流动的新方法——时间推进积分方程法,此法克服了时间线化积分方程法的限制,能较好地模拟激波的运动。本文首先用一维波(?)问题——模型问题阐明此法的基本思想,然后将它应用于二维非定常跨声速流动中。本文还首次引入拟速度位的概念,使时间推进积分方程式得到简化,尾涡条件和Kutta条件更易处理。数值计算表明时间推进积分方程法是合理可靠的。     

高升阻比翼型的设计

使用解析形状函数法和参数化翼型表示方法,将求解绕翼型流场的N-S方程解与优化方法相结合,可设计出新的高升阻比翼型.计算结果表明,设计的新翼型满足设计要求,且获得比原始翼型高得多的升阻比,翼型设计方法是实用的.参数化翼型表示法比解析形状函数法对给定的初始翼型离散数据的准确度影响更小,更适合用于新翼型的设计.     

用虚拟黏性法构造Navier-Stokes方程黏性项的隐式格式

提出的虚拟黏性法是一种关于Navier-Stokes方程黏性项(黏性应力项和传热项)的隐式新方法。通过引入虚拟时间和虚拟黏性项,将隐式格式的构造大大简化,从而避免了大型复杂隐式差分方程组的常规求解。在虚拟时间推进过程中所需求解的方程组的系数矩阵是一个三对角矩阵,它具有计算简单且计算量小的优点。用模型方程和Navier-Stokes方程进行了数值仿真,研究各种参数对计算的影响,并在精度和效率上与显式方法进行比较,证实了算法的正确性和优势。最后分析了本文格式的适用范围。     

基于Reynolds平均的高超声速二方程湍流模型

 在Reynolds平均的框架下推导了可压缩湍流Reynolds应力方程和湍动能方程。根据一定的假设和尺度分析简化并封闭了所推导的湍动能方程。在均匀湍流假设下,湍动能耗散率可分解成可压缩性耗散和旋度耗散,并对其中的可压缩性耗散进行了封闭;同时认为旋度耗散不受可压缩性影响,直接引用传统的Reynolds平均不可压缩湍动能耗散率模型方程。由此构造了适用于高马赫数的二方程可压缩湍流模型。应用所发展的模型计算了高超声速平板绕流,并与若干现有模型的计算、实验与半经验公式的计算结果进行了对比,验证了所发展的模型。在此基础上,通过对压缩拐角的高超声速湍流的数值模拟,对所发展的湍流模型,以及若干现有模型进行了对比,显示了不同湍流模型及可压缩性修正在计算壁面压力分布和热流分布上的特点,说明了湍流模型可压缩修正的必要性,得到了所发展模型的计算结果最接近实验结果的结论。     

一种半球谐振陀螺谐振子动力学建模方法

准确完备的半球谐振陀螺（HRG）谐振子动力学模型是陀螺误差分析的基础。为建立半球壳谐振子动力学模型,基于弹性力学薄壳理论,提出了一种谐振子动力学建模方法。首先,在薄壳的弹性力学几何方程基础上,推导了半球壳谐振子的变形几何方程。其次,在提高受力分析计算精度的基础上,推导了半球壳谐振子的物理方程。然后,分析了谐振子中面的受力平衡关系,推导了谐振子的平衡微分方程。最后,基于以上对整个谐振子的动力学分析,建立了谐振子动力学方程。根据谐振子的不同外载荷形式,利用布勃诺夫-伽辽金法求解得到谐振子2阶谐振状态动力学模型,并得到了谐振子比例系数和2阶谐振频率的表达式。通过对比验证可以看出,参数计算值与实测数据结果一致,证明了所建立的动力学模型的准确性。