叶片安装角槽式轮毂处理的三维数值模拟

采用三维定常不可压Ｎ－Ｓ方程对叶片安装角槽式的轮毂处理的叶栅通道进行了数值模拟,分析了叶片安装角槽式的轮毂处理的扩稳机理,在计算中采用了非交错网格的ＳＩＭＰＬＥ算法,并解决了因使用非交错网格离散方程而带来的压力场振荡的抑制问题,紊流模型采用ｋ－ε两方程模型。     

激波-边界层-分离流相互干扰三维湍流的数值模拟

本文采用数值方法求解时间相关三维可压缩雷诺平均Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程组,模拟激波—边界层—分离流相互干扰三维湍流流动。湍流模型为Ｂａｄｗｉｎ－Ｌｏｍａｘ两层代数模型,改进后用于三维内流问题。采用单元中心有限体积法离散流场控制方程,ＶａｎＬｅｅｒ矢通量格式计算无粘通量,中心差分法计算粘性通量,ＬＵＳＧＳ时间推进格式计算定常流场。本文以二元跨音速扩压器内三流动为算例,数值模拟较强激波—边界层—分离流相互干扰维湍流流动,并与实验结果进行了比较。数值模拟结果,在激波强度、分离点位置和再附点位置等方面,与实验结果吻合较好。     

弯曲管道内湍流流动的数值模拟

为了对弯曲管道内的湍流流动特性进行有效的预测，用显示代数应力模式对－典型的弯曲管道内流动进行了数值模拟和分析，计算结果表明：（1）显式代数应力模式能对弯曲管道内流动的平均量（如速度、压力、摩阻等）进行较为准确的预测；（2）对模式中的ε方程做适当修正后，该模式能较好地模拟出流动曲率效尖对湍流特性的影响，是对雷诺应力模式的一个比较好的近似。此外，该模式在形式上与标准k－ε两方程模式类似，使用方便，具有较强的工程实用价值。     

塞式喷管性能的数值模拟与实验验证

通过特征线法在塞式喷管中的应用，研究了塞式喷管主要结构参数对其性能的影响，同时，研究了塞式喷管的高度特性，最后还对塞式喷管合推力与发动机轴向夹角的高度特性进行了研究和分析。理论研究结果得到了实验的验证。其结果可用于塞式喷管的设计。     

一种适用于大幅值Quasi周期轨道的数值构造方法

针对现有方法在构造大幅值quasi周期轨道时收敛性和通用性不足等问题,提出一种适用于大幅值quasi周期轨道的数值构造方法。首先针对2维不变环面的计算问题,设计了映射间隔约束的环面多步微分修正算法。在此基础上,采用自然参数延拓法和伪弧长延拓法两种不同的思路构造等映射间隔的quasi周期轨道族。接着以quasi-Vertical为对象分析了两种思路各自的特点。最后以quasi-Axial, quasi-DRO等为对象验证了方法对三角平动点及多对复特征值情况的适用性。仿真结果表明,所提出的方法能够简单有效地解决多种类型幅值比超过1.2的大幅值quasi周期轨道的数值计算问题。     

密集型空气雾化流场破碎特征数学模型研究

低速液体射流在高速湍流气体作用下的气液同轴射流雾化流场是瞬态密集型喷雾场,高速湍流中影响和控制雾化的因素很多。针对空气雾化流场液核分布与特征进行分析,开发了模拟液核的随机浸入体模型,结合大涡模拟方法对同轴射流空气雾化喷嘴下游流场进行数值模拟。模拟结果与实验结果对比表明：随机浸入体模型可快速捕捉液核的长度和位置信息。当气液动量比M为3~10000时,能够较准确地预测液核长度,当M>10时,其预测结果远优于唯象模型;该模型能够捕捉回流区、大尺度涡等流场结构。同时,可以准确地预测喷嘴附近的液滴粒径,特别是在气流速度较大（>60m/s）时,液滴平均直径预测误差<10%。     

Fundamental (f) oscillations in a magnetically coupled solar interior-atmosphere system – An analytical approach

Solar fundamental (f) acoustic mode oscillations are investigated analytically in a magnetohydrodynamic (MHD) model. The model consists of three layers in planar geometry, representing the solar interior, the magnetic atmosphere, and a transitional layer sandwiched between them. Since we focus on the fundamental mode here, we assume the plasma is incompressible. A horizontal, canopy-like, magnetic field is introduced to the atmosphere, in which degenerated slow MHD waves can exist. The global (f-mode) oscillations can couple to local atmospheric Alfvén waves, resulting, e.g., in a frequency shift of the oscillations. The dispersion relation of the global oscillation mode is derived, and is solved analytically for the thin-transitional layer approximation and for the weak-field approximation. Analytical formulae are also provided for the frequency shifts due to the presence of a thin transitional layer and a weak atmospheric magnetic field. The analytical results generally indicate that, compared to the fundamental value ($\omega =\sqrt{\mathit{gk}}$), the mode frequency is reduced by the presence of an atmosphere by a few per cent. A thin transitional layer reduces the eigen-frequencies further by about an additional hundred microhertz. Finally, a weak atmospheric magnetic field can slightly, by a few percent, increase the frequency of the eigen-mode. Stronger magnetic fields, however, can increase the f-mode frequency by even up to ten per cent, which cannot be seen in observed data. The presence of a magnetic atmosphere in the three-layer model also introduces non-permitted propagation windows in the frequency spectrum; here, f-mode oscillations cannot exist with certain values of the harmonic degree. The eigen-frequencies can be sensitive to the background physical parameters, such as an atmospheric density scale-height or the rate of the plasma density drop at the photosphere. Such information, if ever observed with high-resolution instrumentation and inverted, could help to gain further insight into solar magnetic structures by means of solar magneto-seismology, and could provide further insight into the role of magnetism in solar oscillations.     

Kink mode m = 1 in magnetic tube with discontinuous magnetic field

In this paper, we study conditions of realization and stability of kink modes with azimuthal wave numbers $m=\pm 1$ in cylindrical plasma flex with twisted magnetic field and homogeneous current along its axis. We assume permanent axial magnetic field both inside and outside the flex, surrounded by currentless plasma. Azimuthal magnetic field decreases inversely proportional to the distance from the boundary beyond the flex. We derived dispersion equations for stable and unstable modes in approximation of “thin” plasma flex. The analysis of equations has been provided for the case of discontinuous axial magnetic field on flex’s boundary. Conditions of propagation of wave modes have been defined. It was shown, that unstable modes can be implemented in certain interval of longitudinal wavenumbers. Results can be applied for the interpretation of solar magnetic tubes behavior, using measurements, provided by spacecrafts.     

NNW-TopViz流场可视分析系统

流场可视化技术采用图形图像直观地表现CFD数值模拟的计算结果,使用户能够方便地对这些数据进行分析、比较和研究。然而,CFD数值模拟的流动复杂,其产生的流场数据规模巨大、数据类型复杂、特征提取困难,传统的串行可视化软件效率低、交互手段单一,难以满足数据分析的需求。国家数值风洞（NNW）工程研制了一套流场数据处理可视化软件系统（NNW-TopViz,简称TopViz）,具有对流场数据处理与特征提取、几何图形绘制等可视化与交互功能。根据可视分析效率需求,TopViz实现了线程并行,在多核计算环境下有效提高了可视化计算和交互效率;针对流场特征提取困难、常规方法效率低的问题,TopViz实现了基于卷积神经网络的流场旋涡特征提取方法,提升了特征提取准确率和效率;为提高软件交互效率并提供便捷的交互方式和体验,基于头戴式显示设备和体感控制器构建沉浸式虚拟显示与交互平台,TopViz实现了手势和眼球凝视2种交互方法,提供沉浸式环境下多视图、多角度流场探测方式。     

翼伞系统纵向飞行性能分析

翼伞系统的飞行性能不仅取决于翼伞本身的气动特性,而且与安装角、伞绳长度、回收物阻力特征、翼载荷等系统参数密切相关。文章应用拉格朗日方程建立翼伞系统的纵向飞行力学模型,对翼伞系统进行飞行力学数值仿真,深入分析了系统参数以及开伞状态对翼伞系统纵向飞行性能的影响规律。结果表明：只有安装角在0°～20°时,翼伞系统才能达到稳定的滑翔状态,且安装角在4°～6°时对应两个稳定的滑翔状态,具体由开伞姿态和速度决定;伞绳特征长度的增加使系统的静稳定性增加;回收物的阻力特征增加6m2,翼伞系统的稳定滑翔角增加15°左右,而迎角减小不到1°;翼伞飞行速度随着翼载荷的增加而增加,其平方与回收物质量成正比。上述结论可为翼伞系统的工程实际应用提供参考。