多块网格计算跨音压气机带叶尖间隙流场

开发了三维数值模拟程序研究轴流跨音转子叶尖间隙流动,应用高雷诺数k-ε湍流模型加壁面函数的方法,计算了轴流跨音转子NASA Rotor37在设计转速下的流场.叶尖间隙采用分区的H型网格和主流区连续对接耦合计算,没有用间隙模型,也没有考虑Vena收缩效应而减小间隙量.在用有限体积法对Navier-Stokes方程和湍流方程进行空间离散的过程中采用了交错网格的方法将N-S方程与湍流方程紧密地耦合在一起,从而提高了计算精度.计算结果和实验数据进行了详细的比较和分析.结果表明,中部叶展具有与实验结果非常一致的流场特征,根、尖区流场则因涡粘假设和激波问题的存在而使流动细节与实验结果略有偏差.     

轴向斜槽式轮毂处理的三维流场分析

首次采用网格分区的方法,对带有轴向斜槽式轮毂处理的压气机叶栅通道流场进行了数值模拟,分析了轴向斜槽式轮毂处理的扩稳机理以及它对压气机性能的影响的原因所在,进一步探讨轴向斜槽的几何参数的变化对其扩稳效果的影响.在计算中采用了非交错网格的SIMPLE算法,求解时均N-S方程,紊流模型采用基于RNG(重整化群)方法的K-ε紊流模型.      

相变储热/辐射器式热沉的传热特性研究

对一种用于舱外活动的相变储热/辐射器式热沉进行了实验研究和数值模拟.设计了相变储热器的实验件和实验方案,进行了储热器的冷却实验.建立了具有导热与对流换热耦合边界的相变储热器数学模型,模型中考虑了自然对流及空穴的影响.编制了相变储热器的性能分析计算程序.通过实验值与计算值的比较分析,验证了程序的可靠性.该模型及计算程序为相变储热器的设计和优化提供了有利的分析工具.     

有限体积法中面积分离散格式的精度分析

讨论了应用顶点中心格式对欧拉方程进行有限体积离散时的精度,详细推导了顶点中心格式在平滑网格和扭曲网格上的截断误差,同时与中心格式和顶点格式进行了比较.格式精度分析采用泰勒级数展开的方法.分析结果表明:顶点中心格式和顶点格式类似,总能达到至少是一阶的精度,并且可以通过加密网格的办法减少误差,而且顶点中心格式更适合于分区计算.通过圆弧亚音速绕流和NACA0012翼型绕流流场的数值计算,验证了分析的结论.     

一种微波及毫米波滤波器的精确设计方法

介绍了一种微波及毫米波滤波器精确设计方法:由滤波器的原型参数,导出滤波器的广义阻抗变换系数;采用模匹配法或耦合积分方程技术等数值算法精确计算各个耦合结的S参数,进而转换成广义阻抗变换系数;通过插值,获得满足原型要求的滤波器各个耦合结的结构尺寸以及谐振腔长度,将结构参数代入模匹配程序、HFSS或CST软件进行分析,分析结果与预期值一致。设计了26GHz矩形波导带通滤波器,模匹配法和HFSS以及CST软件的仿真结果、测试数据与设计要求一致,说明该种设计方法是有效的。     

迎风格式在二维非结构网格中的应用

将原来在结构网格中提出的AUSM+ (AdvectionUpstreamSplittingMethod)格式 ,推广应用于二维非结构网格中 .利用格林 高斯公式对控制体内的变量进行线性重构 ,获得空间二阶精度 ,并采用Barth型限制器以抑制数值解的振荡 .为提高计算效率 ,采用了自适应多重网格法 .最后给出了NACA0 0 1 2翼型和平板激波反射的几个Euler方程的算例 ,并进行了分析比较 .     

基于多点约束的大展弦比机翼静气动弹性计算

基于雷诺平均N-S方程和多块结构网格技术计算大展弦比机翼的气动力.机翼采用梁模型结构,利用有限元方法计算结构变形.建立基于多点约束(MPC,Multi-Point Constrain)的气动、结构数据的双向传递方法.采用MPC数据插值方法快捷地实现自主研发的计算流体力学(CFD,Computational Fluid Dynamics)程序和NASTRAN计算结构力学(CSD,Computational Structural Dynamics)软件的耦合计算,发展出可以考虑结构非线性的静气动弹性CFD/CSD耦合计算方法,该方法比只考虑单方向变形的柔度法更精确.开展了某大展弦比机翼的静气动弹性数值计算分析,结果显示只考虑机翼纵向变形而忽略展向变形对计算结果会有一定影响,而该算例几何非线性影响较小.     

弹性变形对柔性机翼气动特性影响分析

针对大展弦比全复合材料机翼的非线性静气弹响应行为,采用CFD/CSM(Computational Fluid Dynamics/Computational Structural Mechanics)弱耦合方法,使用三维N-S方程和结构力学方程以及动网格生成技术和压强插值技术,通过高精度的数值模拟,求解了大展弦比柔性机翼在结构几何非线性变形状态下的非线性静气弹响应问题.算例分别对典型的大展弦比复合材料前掠机翼和后掠机翼进行了求解,计算结果表明,与前掠翼相比,后掠翼的升力系数明显降低,升阻比降低,严重偏离了刚性机翼的设计点,在柔性机翼的气动设计中必须加以考虑.      

二维非结构网格的生成及其Euler方程解

 采用推进阵面法生成三角形非结构网格,用多重网格技术求解泊桑方程以取得结构背景网格.提出了一种高效的边界剖分法,并用它来生成初始阵面.另外,还提出了一种新颖的数据结构——数组链接表,用它来存贮网格数据.用有限体积法在非结构网格中求解了二维可压缩Euler方程,给出了单段翼型和多段翼型的算例,并与实验结果进行了比较和分析.     

水洞收缩段流场的数值模拟及优选

介绍了一种水洞非对称收缩段的数值模拟方法和程序,计算采用有限体积预估修正格式求解非定常欧拉方程.结合具体水洞设计对几种典型收缩曲线的三维流场进行了计算和分析,给出了沿轴向的压力分布,出口截面速度分布及分布的均匀性,选出最佳曲线,为水洞的设计提供依据.结果表明:本文的计算方法和程序对于水洞或风洞的三元收缩段,是有价值的分析和辅助设计手段.     

使用GAO-YONG湍流方程组 对顶盖驱动方腔流的计算

运用SIMPLER方法、QUICK格式及交错网格技术求解GAO-YONG湍流方程组,对二维顶盖驱动方腔流进行数值模拟,得到方腔内的涡量分布和方腔中心水平剖线上的速度变化曲线,并与DNS(Directly Numerical Simulation)结果进行定性对比,有较好的一致性.结果表明,GAO-YONG湍流方程组中与湍流多尺度现象相对应的机械能方程对此湍流方程组能够合理地模拟出顶盖驱动方腔流中复杂的湍流粘性分布以及湍流能量逆转现象起到了关键作用.研究进一步证明了GAO-YONG湍流方程组能够得到复杂流动的真实粘性场.     

使用GAO-YONG湍流方程组对扩压器流动的计算

采用基于同位网格的SIMPLE方法求解GAO-YONG不可压湍流方程组,对二维扩压器流动进行了数值模拟.通过界面速度动量插值法解决了压力锯齿波问题,并且采用了正交贴体网格和二阶QUICK格式离散对流项以提高计算精度.与实验结果以及BL(Baldw-Lomax)模型计算结果的比较表明,不需要任何经验系数及壁面函数的GAO-YONG不可压湍流模型方程组能够对有压力梯度的湍流流动做出很好的预测.计算发现,在机械能方程中引入平均流压力梯度的作用,对GAO-YONG湍流方程组正确模拟逆压力梯度流动起到了关键作用.      

A method of estimating GPS instrumental biases with a convolution algorithm

This paper presents a method of deriving the instrumental differential code biases (DCBs) of GPS satellites and dual frequency receivers. Considering that the total electron content (TEC) varies smoothly over a small area, one ionospheric pierce point (IPP) and four more nearby IPPs were selected to build an equation with a convolution algorithm. In addition, unknown DCB parameters were arranged into a set of equations with GPS observations in a day unit by assuming that DCBs do not vary within a day. Then, the DCBs of satellites and receivers were determined by solving the equation set with the least-squares fitting technique. The performance of this method is examined by applying it to 361?days in 2014 using the observation data from 1311 GPS Earth Observation Network (GEONET) receivers. The result was crosswise-compared with the DCB estimated by the mesh method and the IONEX products from the Center for Orbit Determination in Europe (CODE). The DCB values derived by this method agree with those of the mesh method and the CODE products, with biases of 0.091?ns and 0.321?ns, respectively. The convolution method's accuracy and stability were quite good and showed improvements over the mesh method.     

NUMERICAL MODELLING OF VELOCITY AND TEMPERATURE DISTRIBUTIONS OF THE BUOYANCY CONVECTION EFFECT IN KNbO3 MELT

Numerical modelling of velocity and temperature fields in high-temperature KNbO3 melt of a loop-shaped Pt wire heater is carried out by using the commercial com putational code ANSYS for the mathematical solution of the governing equations.Based on the experimental boundary conditions and the Boussinesq approximation,the numerical modelling of a steady and two-dimensional model is applied to study the process under consideration of the buoyancy-driven convection condition. The result is compared with the previous experimental and theoretical data obtained in our laboratory, and the former is in agreement with the latter. Thus a theoretical guide for reasonable growth conditions is provided by studying in depth the real fluid flow effects in the crystal growth from the melt.     

通气位置对潜射航行体流体动力特性影响分析

基于均质平衡流理论,通过求解混合介质的RANS方程、SST湍流输运方程和各相之间的质量输运方程,开展了通气位置对潜射航行体流体动力特性影响的三维数值模拟研究,对比分析了不同通气位置条件下空泡形态特性、表面压力分布以及阻力变化特性.结果表明:当两个通气口间距增加到一定值后,在两个通气口之间区域空泡发生断裂,在断裂区域迎流面和背流面表面压力值升高,且在空泡断裂闭合位置出现压力峰值;在相同通气量条件下,随着通气口间距的增大,压差阻力系数和粘性阻力系数均呈减小趋势变化.     

随机振动条件火箭贮箱增压系统工作过程仿真

为研究随机振动对贮箱增压系统动态特性和稳定性的影响,引入管路元件壳体坐标实现外激振动的加载,以谐波叠加法模拟飞行环境的随机振动过程,对管路和管路元件分别建立拉格朗日坐标下的一维流动方程和常微分方程组以模拟增压系统动态工作过程,对某贮箱增压系统有/无外界激励振动条件下的工作过程进行了全程计算.计算结果表明:振动环境下减压器、稳压器等运动元件的瞬时流量变化在增压系统中产生压力波动,某些情况下,会导致减压器等出现具有一个主要频率的振荡过程,影响减压器工作稳定性;振动过程中贮箱气枕压力基本不受影响.无振动条件下的计算结果与试验数据基本吻合.提出的模型和方法在贮箱增压系统工作过程仿真中显示了较好的有效性.     

基于N-S方程串并行计算的机翼优化设计

对典型的大展弦比和中等展弦比机翼在0.7~0.9马赫数区间进行了串行和并行计算的数值校验,证实了以N-S (Navier-Stokes)方程为主控方程的串行和并行流场求解器的正确性,并讨论了并行效率和加速比.结合Powell算法,讨论了在确定的机翼平面形状和翼型的条件下,以升阻比最大为目标的三维机翼截面翼型最大厚度与扭角的优化设计.算例结果表明,厚度的非线性分布和负的扭角会改善机翼流场的流动状态,使机翼的升阻比得到提高,优化设计方法是可行的.      

不同厚度翼型动态失速涡运动数值研究

在低马赫数下, 对三种不同相对厚度的NACA系列基本厚度翼型在俯仰振荡运动中的动态失速现象进行了数值研究.数值模拟时采用双时间法和LU-SGS隐式解法相结合的模式求解了非惯性坐标系下的非定常纳维-司托克斯(NS)方程组, 空间离散采用Roe格式,并结合刚性动网格生成技术.采用全湍流计算,通过引入BL(Baldwin-Lomax)模型计入湍流影响.计算出的气动力迟滞曲线与实验结果变化趋势符合较好,表明了数值方法的有效性,同时通过比较和分析不同厚度翼型在轻失速和深失速下的流场结构,发现绕翼型的失速旋涡产生和发展规律是明显不同的.