绕近地面二维平板非定常不可压缩位势流的精确解

Yih给出了极端地效条件下二维水平平板垂直下落时有势绕流的首阶近似解;给出了同一问题在平板同时还有水平速度情况下的近似解。本文利用椭圆函数将平板和地面围成的区域保角变换至矩形内部,然后在变换平面上求解,得到了有地面效应时水平平板铅垂运动所诱导的势流精确解,以及平板同时还有水平速度情况下的“准定常”连续绕流近似解。     

随机涡法及对钝体后气流结构的模拟

一、问题与方法 本文应用随机涡法,对三棱柱在无界均匀不可压流场中的绕流问题进行了计算。图1为计算简图。      

非结构网格上欧拉方程的自适应逆风隐式有限元解法

本文发展了一种解二维欧拉方程的隐式逆风有限元格式。空间近似应用的是二阶精度逆风方法；时间近似上根据数值通量函数的线性化处理导出了有限元的隐格式。为减少数值振荡，引入了基于特征变量的限制因子。在经自适应处理的非结构网络上，本文对跨、超、高超声速绕流进行了数值模拟。计算实践表明发展的隐式有限元格式与显格式相比具有较高的计算效率和较强的适应性。     

超音速钝体绕流数值计算

本文提出了关于超音速钝体绕流数值求解的一种简便方法——非定常直线法。该方法在一般直线法的基础上利用具有时间项的双曲型方程组,将原求解一阶拟线性混合型方程组的边值问题,化为求解拟线性双曲型方程组的初边值问题。本方法的特点是对初始值要求不高,可克服一般直线法所遇到的奇异性问题,同非定常差分法比较,本方法的算法简单,花费机时少,计算结果也较为满意。用本方法给出了物形有曲率间断的超音速钝体绕流数值计算结果。     

二维可压无粘流的自适应流量修正有限元解

从非定常Ｅｕｌｅｒ方程组出发,采用非结构三角形网格,研究二维无粘流动问题的自适应流量修正有限元解。将ＦＥＭ－ＦＣＴ与网格自适应相结合计算了ＮＡＣＡ００１２翼型的三个算例。计算结果与实验及其他计算结果比较,亚音速和超音速绕流符合很好;跨音速绕流情况存在一些差别。     

翼型和叶栅前后缘局部流场的新精化解法

本文提出一种利用前后驻点自动定位来求前后缘局部流场精化解的方法 ,即在全流场解的基础上 ,把局部流场解分为两部分的叠加 ,其中之一是驻点邻区的近似解析解 ,它吸收了原解中变化急剧的部分 ,使剩余部分成为变化比较平缓、具有更好的解析性态的分布 ,因而易于求准确。算例表明 ,本文方法为提高亚声绕统计算精度 ,特别是带攻角绕流工况 ,提供了有效的手段 ,本文思想也可同变分差分法结合应用     

亚声速圆柱绕流的摄动Galerkin杂交解

本文采用摄动伽辽金杂交方法，应用计算机代数与符号运算求解亚声速圆柱绕流。对给定的马赫数，获得了速度势的近似解表达式，并与Ｊａｎｚｅｎ－Ｒａｙｌｅ－ｉｇｈ方法进行了对比。计算结果表明本方法精度高，计算工作量小，解的表达式简洁。     

高超声速验证飞行器助推分离段流场数值研究

为了模拟有相对运动的多体非定常绕流，发展了基于弹簧近似的非结构动网格方法及耦合动网格的Euler解算器，作为方法验证，模拟了俯爷振动的NACA0012翼型绕流，计算结果与实验及文献结果非常接近。采用二维近似外形，对高超声速验证习行器助推分离段流场进行了数值研究，得到了不同时刻由于多体相对运动形成的干扰流场结构以及分离过程的气动力参数。     

壁面喷流可控环量机翼的数值计算

本文提出一种适用于可控环量机翼研究的势流/粘性流的数值解法,包括钝体分离绕流的一种势流解,以及改进的Spalding统一边界层理论的解。通过对具有壁面切向喷流可控环量圆柱绕流和类椭圆翼绕流所作出的数值计算表明,在其计算模型确定后,就不必依赖实验数据来模拟实际问题,计算结果与实验数据符合较好。     

圆锥低超声速有攻角绕流非对称涡流的数值模拟与非线性分岔的研究

本文用数值方法研究了圆锥低超声速有攻角绕流的对称和非对称定常解，扰动响应以及在更大角时出现的准周期解问题。     

突然起动圆柱的非定常绕流计算

本文用分区方法定量模拟高Re数情况下突然起动圆柱的初期对称流动。绕圆柱的二维不可压流场可分为无粘区、附着粘性区和分离区。无粘区有解析解，用非定常层流边界层方程解附着粘性区，用离散涡方法处理随时间变化的分离区。试用了混合的Rankine和Lamb涡模型及二次离散涡。给出了Re=300、9500两个算例，基本符合有关实验结果。     

跨声速伽辽金有限元法所用上风技术的改进

伽辽金有限元法是解跨声速绕流的有效方法,对超临界流动,为了自动捕捉到激波,通常采用了上风技术,但导致了气流通过激波时质量守恒条件得不到满足。本文通过Bateman变分原理和加罚泛函找到了强制满足质量守恒约束条件的附加项,大大改善了用伽辽金有限元法解超临界绕流的计算结果。     

高阶精度间断Galerkin有限元方法研究

在二维非结构网格上,对高阶精度间断Galerkin有限元方法求解跨音速欧拉方程进行研究。运用间断有限元理论,采用施密特正交化多项式基函数对流场解进行近似描述,使用HLLC近似黎曼解方法计算网格单元边界处的数值通量,积分项通过高斯积分求解,时间推进采用经典四步Runge-Kutta方法,并引入斜率限制器,抑制数值振荡。对NACA0012翼型跨音速无粘流动进行数值模拟,数值结果表明：间断Galerkin有限元方法具有较高的精度,较小的数值耗散和较强的激波捕捉能力。     

三维水下爆炸气泡的数值模拟

在理想不可压势流假定下,采用高阶有限元离散边界积分方程求解水流场。假定泡内气体遵循等熵关系。用Ｅｕｌｅｒ－ｌａｇｒａｎｇｅ方法模拟气泡的演化。计算结果与球泡“精确解”及实验值比较,表明数学模型合理且三维数值方法精度高。计算给出气泡迁移、射流等重要动力学行为和物体受力分析     

非结构网格上新型的NND有限元格式

张涵信发展的ＮＮＤ差分格式是由中心差分格式,二阶迎风格式和一阶迎风格式混合组成的杂交型格式。众所周知,和中心差分格式相对应的是Ｇａｌｅｒｋｉｎ有限元格式。通过对中心型有限元格式加修正项的方法本文成功地构造出二阶迎风型有限元格式和一阶迎风型有限元格式。     

3-D NUMERICAL SIMULATION OF UNDERWATER EXPLOSION BUBBLE

３－ＤＮＵＭＥＲＩＣＡＬＳＩＭＵＬＡＴＩＯＮＯＦＵＮＤＥＲＷＡＴＥＲＥＸＰＬＯＳＩＯＮＢＵＢＢＬＥＬｕＣｈｕａｎｊｉｎｇ（ＳｈａｎｇｈａｉＪｉａｏＴｏｎｇＵｎｉｖｅｒｓｉｔｙ，Ｓｈａｎｇｈａｉ，Ｃｈｉｎａ，２０００３０）Ａｂｓｔｒａｃｔ：Ｔｈｅｐｒｏ...     

高超声速内收缩进气道分步优化设计方法

提出了基准流场与唇口平面形状分步优化的高超声速内收缩进气道设计方法。基准流场以反射激波不均匀性最小和总压恢复最大进行多目标优化设计,使用结合Tayler-Maccoll方程的有旋特征线方法(MOC)进行流场计算,获得双拐点母线内收缩锥基准流场。进气道唇口形状以沿流线积分(Streamline Integral Method, SIM)获得的进气道无黏阻力最小为目标进行优化设计,获得类椭圆形唇口平面形状。针对优化设计结果进行数值模拟,与传统直母线基准流场相比,双拐点母线基准流场反射激波后流动不均匀性下降40%左右,总压损失减少35%左右,总体性能提升明显。类椭圆唇口进气道在设计点的单位质量流量无黏阻力相较于圆形唇口降低6%,具有良好的压缩特性和气动效率,能够减弱进气系统对飞行器气动性能的不利影响。研究结果表明该方法是一种高效且实用的高超声速内收缩进气道设计方法。     

亚临界情况下翼型的非线性压缩性修正

本文从有限基本解方法出发,介绍一种计算亚音速流中非线性压缩性效应的方法。采用逐步逼近的办法,可以把有限基本解方法推广到求解非线性的无粘可压流体,从而为非线性方程的求解提供了新的途径。 编制了FORTRAN语言程序。与精确解和其他理论值的比较表明,结果比传统的有限差分法好,尤其是在翼型头部。     

肋角度对流量系数影响的数值模拟研究

利用有限体积法对三维不可压缩的N-S方程进行离散,对同时带肋和气膜孔出流的内流冷却通道在45°,60°,90°,120°肋的结构下进行了数值模拟。网格划分采用非结构化网格,湍流模型为Realizablek-ε模型,近壁处湍流采用壁面函数,采用SIMPLE算法求解速度与压力的耦合。计算获得了同时带肋和气膜孔出流的内流冷却通道在不同肋角度时的三维流场分布和气膜孔的流量系数值。结果表明同时带肋和气膜孔出流的通道流场结构比较复杂,肋在通道中诱发的二次流改变了气体在进入气膜孔时偏转角度,从而影响了流量系数值。通道中气膜孔的流量系数值随肋角度的增大而呈增大的趋势。计算结果与实验结果进行了比较,两者规律符合良好。      

采用近似边界求解壁板颤振的CBS有限元算法

为了简化Euler方程求解壁板颤振问题,将近似边界条件成功地推广到CBS(characteristic-based split)有限元方法求解壁板颤振领域,把运动壁面的影响转化为Euler方程的边界条件,避免了对动网格的处理.通过对近似边界方法、动网格方法以及3阶活塞理论的数值模拟结果进行对比发现:首先,加载近似边界的CBS有限元方法操作非常简便,且模拟结果与动网格方法完全吻合;其次,Euler方程可以准确捕捉到壁板前后两端不光滑尖点所引起的压力突跳,而3阶活塞理论则无法做到;最后,采用近似边界的CBS有限元方法在求解无黏可压缩流动中细长体小变形的流-固耦合问题上具有潜力和优势.