海军工程学院学报1985年总目次

第一期(总第三十期)求测不规则刚体质量矩阵的方法······················································……朱石坚(1)壁面湍流换热全断面统一温度分布公式···························     

五阶HWCNS在低速复杂流场中的应用

采用五阶精度显式混合加权紧致非线性格式求解雷诺平均NS方程;利用多块对接结构网格技术,对30P-30N多段翼型进行网格收敛性研究。在不考虑转捩的情况下,采用SA一方程湍流模型研究混合加权紧致非线性格式与二阶精度MUSCL格式对该翼型压力分布和典型站位速度型的影响,并与实验结果进行对比分析。采用混合加权紧致非线性格式和SA一方程湍流模型模拟梯形翼高升力构型低速复杂流场,通过对总体气动特性和压力分布的分析,探讨五阶精度显式混合加权紧致非线性格式在低速复杂外形流动中的应用能力。结果表明,对30P-30N三段翼型,采用全湍流模拟方法可以得到较好的压力分布;对梯形翼高升力构型,在附着流和边界层小分离情况下混合加权紧致非线性格式有较好的模拟能力。     

分解炉内三维湍流流场数值模拟

针对某大型高性能分解炉的实际尺寸,用Realizable(带旋流修正的)k-ε模型,采用SIMPLE算法,二阶迎风差分格式,模拟了炉内三维湍流流场.计算所得的速度、压力分布与该炉热工标定的结果以及其实际运行的性能符合较好.     

火灾中受热铜导线断面形貌的研究

在火灾调查工作中,初步估计局部火灾现场的温度是非常必要的.利用扫描电子显微镜(SEM)对不同受热温度下铜导线断面形貌特征进行了研究,并通过金相显微镜对导线微观组织进行研究,得出断面形貌特征与受热温度之间的关系;提出了一种能够根据断面形貌特征判断导线所处火场温度的方法,从而为火灾调查工作提供技术参考.     

超声速欠膨胀冲击射流数值模拟

超声速欠膨胀冲击射流有着重要的实际应用价值,如S/TOVL飞行器、火箭发射、除尘等。其流场结构复杂,包含间断激波、反射激波、马赫盘、滞止泡以及冲击平面传热不同于亚声速冲击射流的特点。为了分析超声速欠膨胀冲击射流流场和传热,采用有限体积法,结合k-l湍流模型以及二阶精度的TVD格式进行数值模拟:对比实验和k-ε湍流模型的努赛尔数,得出k-l湍流模型在传热问题中更具优势;对比阴影图和计算密度云图以及对比冲击平面压力系数的实验值和计算值,验证了k-l湍流模型模拟超声速欠膨胀冲击射流流场的合理性;采用k-l湍流模型研究3种冲击高度(3D,6D,10D),3种压比(2.0,3.4,4.4),3种喷管总温(493 K,591 K,580 K)下,冲击平面温度分布。数值研究结果对分析超声速欠膨胀冲击平面的烧蚀有一定的指导意义。     

碳化硼表面金属化初步研究

为了降低碳化硼(B,C)与金属钎焊连接的温度,采用真空熔结Al基钎料的方法,在1000℃时对其进行了表面金属化处理,SEM和EDAX观察分析B4C与钎料结合断面的形貌和元素面分布,发现B4C与钎料界面发生了互扩散,界面结合致密,形成冶金结合.     

尾流温度场对其速度场的影响与加热减阻

本文利用LDV—02型激光测速仪,在循环水槽里测量了加热和非加热流线型细长旋成体尾流的速度分布和湍流度,并把Hama&Peterson的半经验理论(它适用于非加热体的类似尾流)修正为本文的半经验分析,较好地描述了尾流温度场对其速度场的影响,还从测得的尾流速度数据推断出加热对水下运动物体阻力系数的影响。     

环管型燃烧室喷雾燃烧DES模拟

为研究燃气轮机燃烧室内复杂的两相三维湍流燃烧问题,采用基于Realizable k-ε模型的分离涡模拟(DES)方法对环管型燃烧室进行喷雾燃烧数值模拟.数值模拟结果表明:DES模拟在燃烧区发挥出大涡模拟(LES)的作用,能获取结构复杂的流场、温度场和涡量场,其对流场内涡结构的捕捉和分辨能力明显强于传统的RANS(雷诺时均N-S方法)模拟.与实验数据对比发现,DES模拟所得燃烧产物含量和燃烧室出口温度分布更符合.因此,DES方法可用于复杂形状燃烧室的燃烧数值模拟,能够避免RANS模型对湍流脉动进行平均处理导致耗散过大、抹平涡结构的缺点,从而可获取更准确的计算结果.     

二维凹槽过渡流的DSMC方法模拟

采用DSMC (DirectSimulationMonte -Carlo)方法模拟二维凹槽过渡流动 ,给出了在不同Knudsen数、不同展弦比、不同壁温条件下凹槽流的速度和温度分布 ,并对DSMC方法在模拟全速度域流场时存在的局限性进行了讨论。     

全附体潜艇模型回转运动流场数值模拟

为研究潜艇的回转运动,在RANS方程中引入旋转坐标系科氏力和艇体惯性离心力,并采用分块结构化网格技术,对全附体潜艇模型回转运动的流场进行了CFD数值模拟。首先,推导了旋转坐标系下流场计算的控制方程,给出了科氏力和艇体惯性离心力的表达式,并将其引入到动量方程中;然后,采用RNGk-ε湍流模型以及用户自定义函数方法设置了计算边界条件,数值模拟了全附体潜艇模型(DRAPA SUBOFF)的操舵定常回转运动流场;最后,分析了所计算流场的速度、压力和涡的分布及形态,分析表明计算结果符合其规律性。     

DARPA2潜艇模型定常流动 黏性流场和水动力计算

为了进一步验证SA湍流模式模拟分离流动问题的能力,采用SA(C_υ=30)湍流模型,数值求解定常不可压缩的RANS方程,计算DARPA2潜艇模型定常流动黏性流场和水动力;速度-压力耦合采用SIMPLE方法处理,对流项采用二阶迎风格式,隐式迭代求解方程组。得到的数值计算结果与实验值符合较好。潜艇不同部分对水动力的贡献及受攻角的影响是不同的,升力系数随攻角的增加,艉附体的作用明显增强,而指挥台的作用不是很明显;指挥台和艉附体分别存在时,水动力系数随攻角变化的趋势与全附体时是一致的,曲线不存在奇异性变化。     

脉冲防暴水炮非淹没气体射流的RANS和LES模拟

为提高脉冲防暴水炮刺激剂的雾化水平,必须提高射流的出口速度和改善雾化带形状,而这些因素与管内非淹没气体射流的湍流流场息息相关。为准确描述湍流分布和水柱加速过程,对脉冲防暴水炮的非淹没气体射流进行了RANS和LES模拟。结果表示LES模型较RANS的标准k-ε模型能更为精确地描述湍流分布和管内非淹没气体射流速度的各向异性。但两者对于水柱加速过程的描述基本一致,只是由于LES模型的气液参混较为充分,气液速度相差较小,而可能导致两模型气体在喷口膨胀规律的异同。LES模型气液界面复杂,是用于研究管内气液流型的理想方法。该方法为三维模拟打下了较好的基础,所得结果可以作为管外射流研究的初始条件。     

水下航行器不同燃料燃烧性能的仿真研究

为了解热动力水下航行器不同燃料的燃烧性能,利用FLUENT软件,分别对OTTO-Ⅱ单组元燃料和HAP+OTTO-Ⅱ+H2O的三组元燃料在旋转燃烧室内的两相湍流燃烧进行了数值仿真.湍流计算采用标准模型κ-ε、气相燃烧采用ED模型、液相采用离散液滴模型.得到两种燃料在燃烧室内的温度分布,化学反应速率分布、燃烧产物CO2和CO的质量分数分布以及两种燃料液滴的运动轨迹,并对计算结果进行了比较和分析.结果表明,三组元燃料的燃烧性能要优于单组元燃料的燃烧性能,但由于三组元液滴蒸发较慢,对燃烧稳定性造成影响.     

非均匀湍流场中螺旋桨噪声研究

针对螺旋桨在湍流中的非定常噪声问题,在均匀湍流下螺旋桨宽频噪声相关性法的基础上,推导了非均匀湍流下螺旋桨宽频噪声的公式。在均匀湍流的研究中使用了十桨叶螺旋桨模型,并将计算结果与试验数据进行了对比,发现吻合较好。对于非均匀湍流中的问题,研究发现湍流的不均匀程度不会对宽频谱波峰的位置产生影响,噪声谱波峰的位置仍然取决于叶频,湍流的不均匀程度只会改变叶频处的波峰高度。最后研究了均匀湍流和非均匀湍流中侧斜角对螺旋桨噪声的影响,发现桨叶侧斜角的存在会降低螺旋桨在湍流中受到的宽频噪声。     

纵向涡发生器对棒状燃料元件温度场的影响分析

运用CFD方法,模拟计算了在反应堆棒状燃料元件包壳表面加装纵向涡产生器(LVG)后,燃料棒、气隙、包壳和冷却剂的温度场和冷却剂的流场。计算时采用分块结构网格对带LVG的计算区域划分网格,使用SSTk-ω模型模拟冷却剂的湍流流动,对于控制方程和湍流方程的离散采用二阶格式以保证计算精度。计算结果表明:加装LVG后,包壳和冷却剂间的换热有所改善,燃料棒、气隙和包壳的最高温度有所降低。在LVG附近,燃料棒、气隙、包壳和冷却剂温度均会出现一个峰值,然后急剧下降。冷却剂流动的湍动能突然增强,有效地增强了冷却剂的横向流动,经过LVG后,湍动能逐渐减弱,横向流动速度也逐渐减小。     

低速大雷诺数混合飞艇气动性能分析

针对混合飞艇体积巨大同时气动外形复杂使得现有条件的风洞试验很难精确测量其气动性能的问题,开展了适用于混合飞艇气动性能分析的计算流体力学(CFD)的数值分析方法研究。考虑混合飞艇低速大雷诺数的特点,将变分多尺度方法 (VMS)与动态Smagorinsky大涡模拟(LES)模型相结合,提出了组合的VMSLES湍流模型。将基于RANS方法和LES方法的其他三种湍流模型相对比,利用雷诺数相近、实验数据丰富的6:1长椭球飞艇对不同的湍流模型进行了对比验证。结果显示LES方法预测结果与实验结果吻合较好,优于RANS方法,并能显示更多流动细节,而组合的VMS-LES模型能够更精确地捕获实验研究中观察到的二次涡。利用组合的VMS-LES模型对有翼HAV与多囊瓣HAV进行了气动性能分析,并研究了不同部件对飞艇气动特性的影响。结果表明,由于尾翼表面产生的一次涡与二次涡相互作用,尾翼在增加气动升力的同时也增加了阻力。     

几个寿命分布类封闭性的统一证明

利用累积失效函数给出失效率递增(减)类,平均失效率递增(减)类,新比旧好(差)类,按失效率新比旧好(差)类和按平均失效率新比旧好(差)类的等价刻画。通过与指数分布的对比,给出了描述这些分布类的等价条件。利用特定函数类的复合封闭性,给出以上分布类在加速寿命系统和协同系统中的封闭条件的统一证明。     

单个燃料元件热工水力三维数值模拟

针对在以往反应堆热工水力分析中采用一雏模型计算燃料元件径向传热及采用固定热流密度模拟包壳外壁热流量的情况,建立了运行条件下燃料元件及其周围冷却剂温度场和流场的三堆数值计算模型.考虑燃料元件功率的轴向分布和燃料芯块热传导率随温度的变化,运用标准k-e模型和壁面函数法模拟了冷却荆湍流流动和近壁区域粘性流动.经计算发现燃料芯...     

脊状表面湍流边界层高阶统计量特性

通过对不同风速下V型脊状表面及光洁平板表面湍流边界层内偏斜系数和平坦系数的测试,对比分析了脊状表面边界层内偏斜系数和平坦系数的分布规律。试验在一小型专用风洞中开展,流场测试中使用恒温式IFA300智能型流动分析仪,测试模型则采用有机玻璃材质的矩形平板结构,且试验中模型表面脊状结构的方向与流向一致。最终研究结果表明,脊状表面湍流边界层内偏斜系数和平坦系数分布规律与平板表面基本一致,但脊状结构的存在降低了距壁面无量纲高度y′10区域(包括整个粘性底层和过渡区的一部分)的偏斜系数和平坦系数,而对边界层中过渡区以外区域则影响不明显。由此可以推断,脊状结构主要影响边界层流场的近壁区。     

旋转燃烧室内燃气湍流流动的数值分析

为了解水下航行器旋转燃烧室内燃气的流动特性,利用FLUENT软件,对旋转燃烧室内的气相流场进行了数值模拟。湍流计算采用RNGκ-ε模型、气相燃烧采用ED模型、液相采用离散液滴模型。得到燃烧室不同轴向截面位置的速度变化曲线和燃烧室中心纵截面的湍流粘性系数分布。结果表明,在燃烧室中间燃气形成范围较大的中心回流区,气体切向速度分布呈现Rankine涡结构,径向速度较小且为向心向。由于出口位置的偏心性和空间突缩,燃烧室后部的流场呈现明显的非轴对称分布,而且燃气速度变化剧烈。湍流粘性系数在回流区和出口上游处较大,而在壁面附近都小于0.01。