逆风因子法在非正交曲线坐标系中解不可压N-S方程中的应用

所谓的逆风因子法就是分别用二阶精度公式计算出对流项和扩散项的贡献,然后按此两项贡献的比值来确定一个逆风因子,来保证对流扩散方程数值解的精度和稳定性。这种逆风因子法是作者在用有限元素法解不可压N-S方程中提出来的。本文将逆风因子法应用于在非正交曲线坐标系中用有限差分法解不可压N-S方程。      

螺旋桨升力线理论的匹配渐近展开法分析

在固连于螺旋桨桨叶的非惯性坐标系中,对大展弦比叶片的螺旋奖建立了升力线模型。在匹配渐近展开法分析中引入一假想速度位,导出固连坐标系中的控制方程、边界条件和Ｂｅｒｎｏｕｌｌｉ方程。由匹配渐近展开法的分析表明,在螺旋桨的进距比不是太小的情况下,螺旋桨的升力线理论在固连于螺旋桨桨叶的非惯性坐标系中仍然成立。     

机身边界层控制研究

对带壁面抽吸的三维可压缩机身边界层控制进行了计算研究。采用了非正交曲线坐标系边界层方程,通过修改矢量流函数关系式,将壁面法向速度转变成为控制方程的一个可调参数,从而消除了边界条件与方程未知参量之间的非线性耦合;所发展的横向分段推进数值技术,可用于计算横流速度在横向多次反号的复杂边界层流动;设计不同的壁面抽吸系数分布,以研究边界层控制对多种边界层特性的影响。结果表明:壁面抽吸对边界层的控制效果是显著的,合理的抽吸分布能有效地减小机身阻力。     

耦合Navier-Stokes方程的可压缩边界层数值求解

为了在Reynolds-averaged Navier-Stokes(RANS)方程计算中耦合eN方法进行转捩判断,在RANS方程求解过程中耦合求解了可压缩层流边界层方程为判断转捩提供了精确、可靠的边界层信息.利用等熵关系由RANS方程求出的物面压力分布确定边界层外边界的速度分布,进一步确定出边界层外边界.边界层方程的求解使用Keller提出的二阶BOX方法.为了验证方法求解边界层方程的正确性,在低速流动状态下将计算结果和XFOIL的计算得到的边界层解进行了比较,二者吻合较好.     

弯曲进气道三维边界层预测

本文研究了弯曲进气道内三维可压缩湍流边界层的预测方法。采用非正交曲线坐标系中三维边界层积分方程,并引入掺混系数微分方程,用滞后掺混法求解,数值方法为预估-校正法,用三次样条函数解决了离散数据的插值、求导的精度问题。所编程序可在微机上运行。实例计算清楚地揭示了弯曲进气道内强横向流动对边界层发展的影响,计算所得的进气道总压恢复系数与实验值的误差为1～2％。     

平面大头叶栅跨声速流动的一种有效解法

准确数值预估叶片头部表面的压力分布,一直是叶片气膜冷却设计中十分重要而又尚未圆满解决的问题。本文针对这个问题提出了一个高准确度、高效率的计算方法。 在文献[1]的基础上,本文从任意正交曲线坐标系中的基本方程出发,引入Von Misses变换,导出了相应的流线控制方程,并具体地在拟边界层坐标系中计算了RA叶栅流场,结果与实验符合甚好,据作者所知,本文方法是已见到的计算同类叶栅方法中准确度最高、花费机时最少的。     

二维扩压器和进气道内正交贴体曲线坐标的快速生成法

贴体曲线坐标系的形成是用有限差分法或有限元素法在诸如扩压器和进气道这类任意形状的域内数值求解各类场时,所必需要做的第一步工作。 通常,对这种贴体曲线坐标系有若干要求。首先,所形成的曲线坐标系必须是“贴体”的,以便于精确地施加边界条件。其次,曲线坐标系最好是正交的,以使得计算平面内的控制方程变得较为简单（不出现交叉导数项）,同时也使得便于应用双层零方程紊流模型。第三,所生成的网格系统的网格距要能随意控制,以便在域内计算变量变化     

强粘性层流理论及在粘性流计算中的应用

对高雷诺数流动计算，为了解决部分粘性项甚至全部惯性粘性项落入误差与流动物理尺度的关系问题，本文提出强粘性流动理论。强粘性流中至少有一个粘性项与惯性项同量阶，理论包含了物理尺度各向极限、经典边界层和多层边界层理论为其特例，给出了从经典边界层向物理尺度各向相同极限演化的尺度规律和粘性惯性诸项变化的量阶关系，阐明了粘性与惯性力强相互作用将在剪切层的法向以流向同时“激发”涉尺度结构。对粘性流计算，利用强粘性剪切流尺度律重新标度NS格式的修正微分方程，给出临界网格尺度与流动物理尺度和差分格式精度的关系，得到部分粘性项落入误差和计算结果为非物理数值粘性解的二个判据。并以流场中的边界层、驻点和分离点领域计算为例说明理论的应用，对强粘性剪切流计算、证实部分粘性项甚至全部粘性惯性项落入误差的问题值得重视。     

轴对称弯道流动的N—S方程数值解

利用涡量—流函数方法求解离心压气机进气弯道中的轴对称粘性流场。本文推导出了q—τ湍流模型方程在正交曲线坐标系下的表达式,并用以预测湍流流场。文中在涡量方程及q—τ方程中引进非定常项,采用时间推进法求解。通过对二维直通道及复杂边界的轴对称弯道内的层流、湍流流动的计算,获得了令人满意的结果。      

基于轴对称比拟的高超声速复杂外形气动热预测方法研究

轴对称比拟法是高超声速飞行器气动热计算的一种有效方法。针对轴对称比拟法应用时公式推导繁琐、变量迭代复杂、计算量大的缺陷,直接在笛卡尔坐标系下采用三维线性方程拟合物面方程及流场变量,推导出相应的流线尺度因子计算公式,扩展了驻点区热流密度的计算方法,提出了一种驻点区下游流线推进格式。以改进的轴对称比拟法为基础,将边界层外无粘流场数值方法与边界层内气动热工程算法进行耦合,发展了一套适用于三维复杂外形飞行器的气动热计算方法。通过对球头钝锥和双椭球算例进行验证,结果表明:方法计算效率较高,适用范围较广,热流计算结果和实验数据吻合良好。     

旋转对三维边界层的影响研究

总结了目前已知的一些旋转的基本现象。在对三维旋转不可压边界层方程进行认真地分析的基础上,详细、深入地讨论了旋转产生的离心惯性力、哥氏惯性力以及法向压力梯度不为零(?p/?x2≠ 0 )对三维边界层计算结果的影响,得到了一些具有一定参考价值的结论。      

N-S方程数值模拟方法研究及边界层中相干结构的非线性演化

建立了一套三维不可压流动的直接数值模拟方法。该算法在x及y向构造了基于非等间距网格的紧致有限差分格式和非线性项的迎风紧致型格式,在z向采用Fourier谱方法。并研究、提出了三维、非定常流体运动下游边界的无反射出流条件。另外,本文还构建了多涡结构的初始扰动理论模型,利用直接数值模拟的方法研究了边界层中多涡结构的非线性演化特点。结果表明多涡结构非线性演化中的许多特点,如雷诺应力分布、高剪切层的形成、喷射和扫掠等与壁湍流中相干结构的发展规律和现象十分相似。同时,通过对边界层中多涡结构的非线性演化问题的直接数值模拟,验证了该算法具有计算精度高,稳定性好,收敛速度快,出流边界影响小等优点。     

边界层流动的非平行稳定性研究

采用极为有效的抛物化稳定性方程（PSE）方法研究边界民支的非平行稳定性。并利用抛物线坐标变换、有限远延拓边界条件、及流向和法向全差分的数值方法，对不同频率的二维Tollmien-Schichting波的非平行稳定尾进行了计算和分析，计算结果与Orr-Sommerfeld方程（OES）的解以及谱配置方法的PSE结果作了详细的比较，得到了满意的结果。     

旋转叶片上的三维附面层

在附壁的非正交曲线坐标系下,见图1。张量形式的流体运动方程见参考文献。 我们假设:(1)流体呈定常运动;(2)忽略彻体力;(3)TSL条件成立。在基本方程中,引入逆变物理分量代替逆变分量,并根据假设进行简化和量级分析,得到一组普遍适用的附面层方程:     

高精度ENO格式在喷管流动模拟中的应用

本文分析了ＥＮＯ格式的特点,并应用于Ｅｕｌｅｒ方程和全ＮＳ方程的迁移项和压力项,模拟了各种喷管流动。首先计算了ＪＰＬ喷管流动,取得了和实验一致的结果,分析了喷管内产生弱激波的原因;其次计算了非定常二维喷管流动,给出了相应的计算结果,并分析了壁面附近流场的拓扑结构;最后计算了摆动喷管的流动及推力矢量计算,取得了推力随摆角变化的规律,和实验结果表明基本一致。     

应用GAO-YONG湍流模式数值模拟三维激波/湍流边界层干扰

应用GAO-YONG可压缩湍流模式数值模拟了三维激波/湍流边界层干扰算例之一——单鳍流动.攻角20°,来流马赫数2.93,雷诺数9.8×105.对流项和扩散项分别采用Roe格式和二阶中心差分格式计算.Runge-Kutta显示时间推进方法求解了半离散的控制方程.包括壁面压力分布,边界层内流动偏移角等在内的计算值与试验数据进行了比较.准确地预测出了三维激波/湍流边界层干扰流场的主要流动特性——λ波结构,主分离涡核,膨胀区,滑移线等.计算与Alvi等提出的单鳍流动的理论模型符合很好,得到了平板表面压力以及分离线、再附线等在单鳍流动中所独有的半圆锥特性.      

胶接连接件应力分析新模型与新方法

在工程梁理论基础上根据力法变分原理得到薄壁结构胶接接头应力分布的解析闭合解,克服依照传统剪滞理论所得之结果不能满足边界条件的缺点,指出使胶接接头失效的真正原因是接头端部的剥离正应力,从而为胶接接头设计提供正确依据。     

基于非Favre过滤的动力学非线性亚格子模型

采用直接过滤的Navier-Stocks（NS）方程组作为可压缩湍流大涡模拟控制方程组。方程组中因过滤产生的高阶相关项用Taylor级数展开近似表达,并仅保留级数的一阶导数项。因为截断级数产生的误差相当于丢失了模型的部分耗散作用,采用一种动力学非线性亚格子模型来补偿丢失的耗散影响。用模型模拟了超音速湍流平板边界层及绕双椭球高超声速湍流,并将所得结果与经验公式和计算结果进行了对比。对比表明模型可以有效的模拟可压缩湍流流场。     

喷注壁面楔状流层流边界层速度与温度相似解

建立了壁面有喷注的楔状流层流边界层冷却数学模型,求解经相似变换得到的描述无量纲流函数和无量纲温度的常微分方程,获得了楔状流层流边界层无量纲速度和温度的相似解,给出了考虑壁面喷注边界条件的形式简洁、物理意义明确的无量纲流函数拟合式;用Runge-Kutta法求解无量纲速度与温度的常微分方程,获得了壁面有喷注的楔状流层流边界层无量纲速度和温度随楔形角、吹风参数、冷却介质温度的变化规律.通过计算0°,18°和36°楔形角的楔状流层流边界冷却速度与温度分布得到以下结论:楔状流速度边界层和温度边界层厚度都随楔形角的增大而变小,随着吹风参数的增大而增大;当冷却介质温度越低、吹风参数越小、楔形角越大时,靠近壁面处楔状流层流边界层内温度梯度越大.      

半波GLM方法封闭的湍流方程对转捩过程的数值模拟

ＴＲＡＮＳＩＴＩＯＮＳＩＭＵＬＡＴＩＯＮＵＳＩＮＧＧＬＭ－ＥＭＢＡＳＥＤＣＬＯＳＥＤＥＱＵＡＴＩＯＮＳＯＦＴＵＲＢＵＬＥＮＣＥＧａｏＧｅ（４ｔｈＤｅｐｔ．ＢｅｉｊｉｎｇＵｎｉｖｅｒｓｉｔｙｏｆＡｅｒｏｎａｕｔｉｃｓａｎｄＡｓｔｒｏｎａｕｔｉｃｓ，Ｂｅ...