三维机翼可压缩层湍流边界层计算

 从三维时间相依可压缩边界层动量积分方程和平均流动能积分方程出发,将二维可压缩层、湍流边界层积分方程算法推广到求解有限翼展后掠机翼的边界层流动。利用四步Runge-Kutta时间步进方案数值求解积分方程组,并利用当地时间步长加速迭代收敛。文中分析了数值方法的稳定性与收敛性,并考查了横向耗散项对计算结果的影响。算例表明,能获得令人满意的三维机翼定常可压缩层、湍流边界层的计算结果。     

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

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

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

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

跨声速翼型绕流的Euler／边界层方程干扰数值解

本文利用Euler方程和可压缩湍流边界层积分方程研究绕跨声速翼型的有粘与无粘强干扰流动。应用有限差分法在贴体的网格上求解时间相关的Euler方程,以剪功积分方法求解翼面贴附和分离湍流边界层流动,并引入一个松弛方程描述剪应力对上游湍流历程的延迟响应。有粘/无粘干扰采用表面源模型。计算结果表明,对翼面存在强干扰流动情况,获得了与实验值基本吻合的结果。     

网格湍流CAARC模型风洞实验

对大气湍流边界层的真实模拟是风工程风洞模拟实验所要满足的基本条件之一，平均速度剖面，湍流度剖面，积分尺度和风谱是反映大气湍流边界层流动的四个最基本的因素。本文首先研究了不同孔隙率的均匀网格湍流场的流动特性，得到湍流度，积分尺度以及风谱的变化规律，然后在变化基中某一个因素而保持另外三个因素不变的条件下，研究了处于均匀网格湍流场中的CAARC模型风荷载响应的变化，从而确定每个因素的独立变化对模型所受风荷的影响。实验结果表明湍流度和积分尺度对结构物所受风载都有很重要的影响，在风洞大气湍流边界层模拟中应该予以充分的考虑，否则将会引起风荷载实验结果的偏差。     

边界层中湍动能和耗散能最大的尺度分量特征研究

阐述了利用离散正交小波理论分析 1 -D湍流测试信号的方法,使用热线风速仪,在 3个雷诺数下,分别测量了光滑壁面湍流边界层的速度和加速度信号。运用所述方法分析了边界层测试信号,给出了最大湍动能密度和最大耗散能密度及其对应尺度沿 y+ 的分布曲线,并且给出该尺度分量的发生频率沿 y+ 分布曲线。分析可见 :两种最大能量密度的无量纲尺度沿 y+分布具有雷诺相似性,其发生频率随着雷诺数增加而增加,但是最大湍动能密度尺度的频率沿 y+基本是一常数;无量纲的最大湍动能密度沿 y+在边界层小部分区域具有雷诺相似性,而无量纲的最大耗散能密度沿 y+ 在整个边界层内具有雷诺相似性。     

湍流的耗散及弥散相互作用理论

 湍流具有涡团散裂、耗散和弥散的特性。根据Prigogine倡导的耗散结构理论,推导了表征耗散与弥散相互作用的新的湍流控制方程组。其特点是:用稳定性分析得到湍流动能产生项,再根据广义熵增原理推出并列存在的,分别适用于强弱涡量的二个湍流动量方程。运用该理论已成功地计算了一些典型的湍流问题:湍流边界层中的马蹄涡拟序结构、钝体涡尾区的湍流能量逆转、湍流涡团散裂弛豫及各向异性分布。文中给出了部分算例。     

低雷诺数翼型和跨声速翼型边界层计算

本文利用边界层动量积分方程和平均流动能积分方程兼容计算了翼面的层流和湍流边界层流动。文中采用一种e~9型转捩判别公式预测翼面存在层流分离气泡情况的转捩位置。并引入对剪应力系数C_r的滞后方程以体现湍流中雷诺应力的历程效应。为避免在翼面流动分离时边界层方程的奇异,采用反解法计算。算例表明,计算与测量结果吻合良好。     

几种跨声速边界层积分方法的比较

本文研究了几种基于不同类典型速度型有代表性的层流和湍流边界层积分方法，通过对这些方法的理论基础的考察以及对计算结果和实验结果的对比，讨论了它们的特点，评价了它们的适用性。     

后向台阶湍流边界层分离、再附及发展

应用激光测速仪对后向台阶湍流边界层分离、再附及发展进行了测量,得到时均流速、流向湍流度、平坦系数和偏斜系数等沿程分布。依据时均流速分布特性,得到分离区边界及再附点位置ｘ＿ｒ＝５．２ｈ＿ｓ（ｈ＿ｓ为台阶高度）,再附后湍流边界层恢复速度型相似的位置ｘ＿Ｔ＝１３．５ｈ＿ｓ。在详细分析湍流边界层区域时均及紊动特性的基础上,与光滑平板湍流边界层进行了比较。结果表明,除流速分布对数公式中积分常数较大外,无其它差别。此外,湍流边界层区存在尾流律,尾流参数π＝０．２９。     

二维非定常跨音速时间线化积分方程的推导

1.引言 求解非定常跨音速流动的主要困难是非线性问题。对于微幅翼型振动问题(同时引起激波的微幅振动)可做时间线化简化处理,得到时间线化微分方程定解问题。时间线化积分方程是由时间线化微分方程推导出来的。     

动量矩方程的剪切积分解析表示式

本文采用两层涡粘性模型和速度剖面的指数律,导出了可压缩动量矩方程的剪切积分解析表示式,从而大大简化了边界层计算。和其它理论结果相比,本结果更接近于实验数据。     

二项分布参数整体推断方法

针对工程中常见的二态数据,提出了二项分布参数整体推断方法,该方法可以很好地确定二项分布参数与状态变量之间的数学关系式,给出了整体推断的点估计和置信限、置信区间估计公式。它能够将不同状态下的二态数据作为一个整体来进行统计推断,全面开发利用了不同状态下二态数据之间相互提供的横向信息,从而使其可利用的信息量远远大于只能分别在各自状态下对二态数据进行统计分析的传统方法。所以,在精度相同的情况下,可以节省大量试样;而在试样数一定的条件下,又能提高预测精度。最后给出了应用实例。     

二维复合广义J积分

本文推导出了含复合型裂纹的变厚度板在承受体力和面力、存在非均匀温度场的情况下的能量差率表达式。并由此定义了复合型广义J积分 J=J₁cosφ+J₂sinφ J₁=lim（1/t）{∫Wtdx₂-∫W（t/x₁）d A-∫S_i（u_i/x₁）tds -∫（W/T） （T/x₁）td A-∫Bi（u_i/x₁）tdA} J₂=lim1/t{-∫Wtdx₁-∫W（t/x₂）dA-∫si（u_i/x₂）tds -∫（W/T）（T/x₂）tdA-∫Bi（u_i/x₂）tdA} 在此定义的基础上,我们证明了广义J积分的与路径无关性,并利用其路径无关性和裂纹尖端的应力场奇异性得到了广义J积分与线弹...     

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

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

用四坐标数控铣床加工整体叶轮

介绍一种形成包络面(等距面)的实用方法,提出一种用四坐标数控铣床加工整体叶轮的程序计算方法。还介绍了刀具选择、叶轮装夹和干涉检查等内容。根据本文所述方法用四坐标数控铣床在国内首次加工出导弹涡喷发动机整体叶轮,各项指标均达到设计标准。     

Navier-Stokes方程积分表达式的有限元解法

本文采用定常不可压Navier-Stokes（NS）方程的积分表达式方法计算粘性流动问题,用线性三角形单元划分计算区域,建立了适当的迭代求解过程。本文改进了压力计算方法,用全变量松驰加快了收敛速度,用矩阵分块计算减小了内存限制,并可以通过对矩阵特性的分析更好地控制迭代过程的稳定性。     

旋转变厚度圆盘的样条积分方程法

一、前言 对于旋转变厚度圆盘的应力分析,一般可用有限元法或差分法,但相对于边界元来说精度较差,若要提高精度就必须将单元分细,随之而来的是方程阶数的增加和计算时间的增长。为了克服上述缺点,本文提出了一种新的样条积分方程法,直接利用弹性力学二维问题的环基本解建立积分方程并用样条函数进行求解。这样只需划分少量单元,就可获得精度很高的数值解,可大大节约计算机内存和计算时间。