拖缆系统直线定常运动仿真计算

推导了一个三因次拖缆系统动力平衡运动方程式 ,采用中心差分方法数值离散拖缆系统动力平衡方程 .根据离散方程编制了拖缆系统直线定常运动程序 ,在设定拖缆尾端自由边界条件及拖点边界条件下 ,数值仿真计算 30m拖缆的阵形 .根据仿真计算讨论了拖缆的材料刚度、拖缆水中重量、拖点速度、水动力系数主要参数等对拖缆阵形的影响 .拖缆系统直线定常运动的仿真计算为拖缆系统设计可提供参考依据 ,同时通过该仿真计算及讨论 ,为拖缆系统非定常运动仿真计算奠定了基础     

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

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

滑移边界条件的收敛性分析及应用

用速度滑移与温度跳跃边界条件代替通常假定的无滑移边界条件,可有效地提高计算流体力学模型对高空滑移流区域流动的预测精度。应用Maxwell滑移边界条件时,通过直接计算速度梯度及温度梯度而得到速度滑移和温度跳跃量的处理方法在网格较密的时候会出现迭代计算发散的问题。理论分析表明,直接计算梯度的方法使边界条件的时间推进过程等价于雅克比迭代过程,因此必须满足相应的收敛性条件。为了消除收敛性条件的限制,给出了一种在任意网格密度下均收敛的边界条件处理方法并通过数值算例验证了该方法的正确性。针对高空高超声速流动,以空天飞机为例,对比了滑移/无滑移边界条件所得结果的差异,分析了滑移效应对飞行器气动特性及热环境的影响。     

基于静叶不同安装角的增压锅炉烟气涡轮多叶片通道数值模拟

基于定常假设,运用Fluent软件,利用混合平面模型传递级间参数,对静叶不同安装角下的烟气涡轮多叶片通道流场进行了数值计算,得到了温度、压力、出口速度和马赫数等参数的变化规律,绘制了静叶不同安装角下各种参数的体积加权平均XY图。     

柴油油滴蒸发和燃烧过程的理论计算

提出了柴油油滴蒸发和燃烧过程理论计算方法.对油滴蒸发考虑了不定常蒸发和定常蒸发2个阶段,提出了一个求解蒸发时,油滴直径瞬时变化的代数方程式.对一台直喷式柴油机进行了实例计算,计算得到的瞬时燃烧放热率和气缸压力,与试验值符合很好.     

以应变能释放率判据为基础的Ⅰ-Ⅱ混合型裂纹相变增韧计算

以应变能释放率判据为基础 ,采用压力敏感准则和权函数法对相变增韧陶瓷Ⅰ -Ⅱ混合型裂纹的增韧效应进行了理论计算 ,分别给出了静止裂纹和定常扩展裂纹相变塑性屏蔽的理论表达式。结果表明 :相变对静止裂纹有负屏蔽效应 ,并随KⅡ/KⅠ 的比值增大而增大 ;对定常扩展裂纹的增韧结果除与材料弹性模量、相变尾区高度和相变体积分数有关外 ,还随KⅡ/KⅠ的比值增大而减小。     

基于传输方程空化模型的定常自然空化流场数值计算

采用基于传输方程的空化模型,应用自主开发的软件对水翼和水下回转体定常自然空化流场进行了数值计算。采用基于压力修正和多块结构化网格的有限体积方法数值求解Favre平均的NS方程,运用k-ε两方程加壁函数的湍流模式封闭雷诺应力。计算结果与已有实验结果符合较好。     

等离子体粒子模拟中的电磁场边界条件

针对等离子体粒子模拟中电磁场边界条件的选取问题 ,对简单吸收边界条件、Lindman边界条件、超吸收边界条件等进行了详细的讨论。结果表明 ,简单吸收边界条件计算量最小、误差最大 ,超吸收边界条件计算量最大、误差最小 ,Lindman边界条件误差较小且其反射系数对入射角的依赖性不大     

高空滑移条件对化学非平衡粘性流动的影响

本文用数值方法计算了不同高度滑移边界条件对钝体绕流流场特性的影响。采用化学非平衡粘性激波层控制方程,边界条件含壁面滑移、激波滑移和无滑移条件。计算结果表明,考虑滑移条件后,激波层较厚,激波处的电子数密度为一有限值。文中给出了计算结果,并与其它文献的结果作了比较。     

鳍对拖式吊舱推进器水动力性能的影响

为考察鳍对拖式吊舱推进器水动力性能的影响,用基于速度势的面元法建立了吊舱推进器水动力性能计算的数学模型,对附鳍的吊舱式推进器的定常及非定常水动力性能进行了预报.螺旋桨和吊舱及鳍之间的相互影响通过迭代计算来处理.采用计算较为简便的关于扰动速度势的基本积分微分方程,并采用双曲面形状的面元以消除面元间的缝隙.在桨叶随边满足压力Kutta条件,用Yanagizawa方法求得物体表面的速度分布.计算了无鳍、单鳍和双鳍拖式吊舱推进器的水动力性能,计算结果表明,鳍会产生推力,使得双鳍时整个推进器的水动力性能最好,单鳍次之,无鳍最低.     

二元翼型跨音速非定常测压实验研究

本文描述了在南航ＮＨ－１风洞中,对ＮＡＣＡ００１２翼型绕其１／４弦长处轴作俯仰正弦振动时,跨音速非定常压力测量的实验。实验目的在于研究在跨音速风洞中二元翼型非定常压力实验中的技术问题,包括实验装置的设计,数据采集和处理等,并初步研究了Ｍ数、攻角、振幅等参数对非定常压力分布的影响。     

多尺度方法在强非线性系统中应用的一个注释

多尺度法在求解非线性系统时，目前所得到的结果多集中在弱非线性系统上。论文对多尺度法进行了改进。对强非线性杜芬系统进行了变换，得到了强非线性杜芬系统存在3个解的结论，其中大的和小的是稳定的，中间的是不稳定的，物理上实现哪个稳定运动取决于初始条件。     

液管道系统动力响应的研究进展

综述了自70年代以来关于输液管道系统的动态响应问题的研究进展情况.根据管路结构的特点,将问题分为两端支撑直管、悬臂管、不稳定流直管、圆柱壳管和曲管5个部分分别进行了阐述.并在此基础上,概括了应用于该领域的一些计算方法.最后,在总结研究成果的同时分析了这一领域的研究趋向.     

球锥体高超声速绕流的俯仰阻尼导数的数值计算

采用Navier Stokes方程描述球锥体外形非定常振动流场 ,在Etkin理论下给出俯仰阻尼导数的计算公式。对定常流场的计算采用ADI形式的NND格式 ,对非定常流场的计算采用四步Runge Kutta方法 ,并引入变系数残值光顺技术加速收敛。将俯仰阻尼导数计算结果与实验及工程计算进行了比较 ,并数值研究了其随振动中心的变化规律。     

势流理论的一种数值方法在空泡形态中的应用

基于势流理论的假设,采取常规边界元法对具有自由面的定常平板和轴对称空泡流动数值方法进行了研究,运用Riabourichinsky镜像空化模型,建立了一种空泡自由面求解的优化迭代方法,以自由流线长度为参数求解模拟了定常无重工况下流经二维平板和轴对称圆盘的空化流动.数值计算结果与解析解的相对误差不超过3%,证实了轴对称体空泡外形的近似回转椭球体假说在空化器附近是不成立的.计算结果与同类文献数据的对比显示了边界元法进行空化流场计算具有着显著的优越性,所采用的迭代方法精度高、收敛速度快,所得结论对深入了解空化特性具有参考价值.     

三维超声速开式空腔振荡特性研究

采用时间3阶、空间5阶的高精度、高分辨率格式,对稳态来流条件下超声速开式空腔的非定常振荡特性进行数值研究。通过二维和三维空腔压力变化、声压级、振荡频率等参数的计算和分析表明:空腔内存在典型的压力振荡,振荡有近似的周期性,且变化规律较复杂,三维空腔对主振荡频率影响较小,但振荡强度有较大的不同。     

非结构动网格分区并行计算方法

针对包含运动边界的非定常流动问题,基于弹簧近似法设计了一套非结构动网格分区并行计算流程,引入了METIS软件进行非结构网格二次剖分,研究了对接边界通信类型及其处理方法,基于消息传递并行编程模式实现了非结构动网格求解器的并行化。测试结果表明:该求解器具有较高的并行效率,网格变形模块对非结构动网格并行计算性能具有重要影响。     

水下应用栅格翼动态展开参数预示方法

基于水下应用栅格翼动态展开过程动力学模型,根据展开特征角度定常水动力方法,拟合获得展开全程受到的流体力矩,并引入考虑相对速度影响的修正因子,形成水下应用栅格翼动态展开过程参数工程预示方法,以典型展开时序点航行体运动参数为设计输入,对栅格翼展开过程运动参数进行预示。通过与水下应用栅格翼非定常流场仿真计算数据以及水下航行体弹射试验数据对比,验证了上述预示方法的正确性及工程适用性,为水下应用栅格翼方案设计的优化及展开不同步性分析提供设计参考。     

钝体高超声速气动加热与结构热传递耦合的数值计算

气动加热与结构热传递耦合问题在航天和工程应用领域非常重要。分别采用松耦合与紧耦合方法,数值模拟了高超声速二维圆管绕流的流场与结构传热耦合的非定常过程。在紧耦合方法中,流场部分采用基于Navier Stokes方程的有限体积法,将AUSM+格式与时间方向的显式多步Runge Kutta法结合;结构传热部分采用基于二维热传导方程的Galerkin有限元法。流场与结构区通过交界面的热流和温度边界条件实现耦合。计算结果分别与实验、文献做了对比,结构内部温度变化关系以及壁面的热流分布均较好地吻合。两种耦合方式的计算结果对比表明,对于流场特征时间远小于结构传热特征时间的问题,松耦合方法计算效率高,精度与紧耦合方法接近。     

An M/M/1 queue with delayed feedback

We present some results for M/M/1 queues with finite capacities with delayed feedback. The delay in the feedback to an M/M/1 queue is modelled as another M-server queue with a finite capacity. The steady state probabilities for the two dimensional Markov process {N(t), M(t)} are solved when N(t) = queue length at server 1 at t and M(t) = queue length at server 2 at t. It is shown that a matrix operation can be performed to obtain the steady state probabilities. The eigenvalues of the operator and its eigenvectors are found. The problem is solved by fitting boundary conditions to the general solution and by normalizing. A sample problem is run to show that the solution methods can be programmed and meaningful results obtained numerically.