基于预处理HLLEW格式的全速域数值算法

基于HLLEW(Harten-Lax-Van Leer-Einfeldt-Wada)格式引入预处理技术发展适合求解全速域流场的三维Navier-Stokes求解器.引入低速预处理技术,重新构造HLLEW格式的耗散项,给出预处理后的HLLEW格式,并根据预处理后的雅克比矩阵构造相应的隐式时间推进方程.利用预处理方法求解NACA 4412低速不可压流动与RAE 2822跨声速可压缩流动,并与实验结果及原有方法的计算结果对比.结果表明:预处理HLLEW格式不仅提高低速不可压缩流动的计算效率和精度,也保持了对可压缩流动的处理能力,是一种适用于全速域流场数值模拟的有效方法.     

最小二乘有限元法求解非定常应力的Navier-Stokes方程

为精确求解非定常层流问题,发展一种非定常速度-应力-压力的方法.采用牛顿法对非线性对流项进行线性化处理和预处理共轭梯度法,实现了非定常应力形式Navier-Stokes方程的求解.方腔层流流动比较发现,非定常应力形式比涡量形式与试验结果更加吻合,精度更高.该方法有效地解决亚格子应力项的问题,实现基于最小二乘有限元法的湍流求解.比较方腔湍流流动的试验与仿真结果,证明本文的方法具有可行性,为湍流大涡模拟计算打下基础.     

低速流动的预处理多重网格求解

本文采用预处理的时间推进方法,求解可压形式的Euler/Navier-Stokes方程,数值模拟低马赫数范围内的无粘、粘性流动.为了进一步加速收敛,并使用了多重网格技术.计算表明预处理能够保持可压缩求解方法在高速时的收敛特征,对于低速流动计算可以加速收敛并提高精度.     

JFNK在高温堆扩散计算中的应用

研究了JFNK框架下高温堆中子扩散问题的求解方法。研究结果表明,JFNK方法在求解与源迭代相同形式中子扩散方程时,相对残差下降趋势为逐渐加快并趋于稳定,有利于更高求解精度的实现。使用通量归一化附加方程可以获得更好的JFNK非线性迭代特性,但在算例中其部分牛顿修正方程求解时间偏多,总计算时间高于显式有效增殖系数附加方程法,需要研究更高效的JFNK预处理方法对线性求解环节进行改善。     

基于层流动能湍流模型的数值模拟方法

从层流动能出发提出数值模拟原则;综合考虑自然、旁路和分离流转捩的因素构建实用的层流动能湍流模型,结合预处理和基本求解技术发展出适于转捩流动的数值模拟方法和程序.针对预处理技术,以Weiss-Smith矩阵为基础,考虑湍流粘性的影响;针对基本离散格式和边界条件,结合模型方程进行对角占优强化等特殊处理.最后通过平板边界层和典型翼型,特别是低雷诺数翼型的数值模拟,验证数值方法的有效性和鲁棒性.算例表明本文的方法能够为求解更复杂的流动提供参考.     

任意马赫数非定常流动数值模拟的统一算法

发展适用于从低速到高速任意马赫数非定常流动数值模拟的统一算法.通过引入一个伪时间导数项和一个新的预处理矩阵,得到双时间非定常预处理可压缩Navier-Stokes方程.方程的对流项采用三阶Roe通量近似差分格式离散,粘性项采用二阶中心差分格式离散.基于数值通量的线性化技术,实现伪时间步的隐式ADI-LU格式迭代,进而获得物理时间步的二阶推进精度.重点以低马赫数流动为例,求解了圆柱绕流和NACA0015翼型等速上仰动态失速问题.计算结果表明该统一算法能够较好地模拟低马赫数乃至任意马赫数非定常流动.     

求解非定常不可压N-S方程的预处理方法

应用预处理技术,对不可压非定常N-S方程使用双时间推进法求解.当沿物理时间层推进时,连续性方程和动量方程沿伪时间方向使用隐式线Gauss-Seidel迭代法求解,对流项采用三阶迎风差分法离散.通过对不同Reynolds数、不同深宽比下非定常驱动腔内流动的模拟,数值研究了预处理法计算非定常不可压粘性流动的收敛特性,分析了沿伪时间层的迭代收敛速度对流场Reynolds数的依赖特征.     

一种隐式预处理方法及其在定常和非定常流动数值模拟中的应用

将Choi-Merkle矩阵预处理方法与LU-SGS隐式方法、双时间法以及多重网格方法结合,发展适用于绕飞行器定常和非定常粘性流动的高效隐式预处理计算方法和程序.介绍一种针对定常和非定常流动的LU-SGS隐式预处理方法的统一表述方法.在不改变流动解的前提下,对Navier-Stokes方程的伪时间导数项实施Choi-Merkle矩阵预处理,从而改善可压缩控制方程在低速情况下的系统刚性,使基于LU-SGS时间推进格式的数值模拟方法同时适用于从极低马赫数到可压缩范围内的数值模拟.对Jameson中心格式的人工粘性进行相应的修改,以提高低速流动的计算精度.翼型、机翼以及翼身组合体绕流的数值模拟研究表明,隐式预处理方法获得了很高的计算效率,可使马赫数0.1左右的低速流动计算时间减少50%以上;通过对现有可压缩计算程序进行小量改动,便可使其均匀覆盖整个低速流动范围,提高CFD程序在飞行器绕流数值模拟中的实用性.     

完全变光滑长度SPH法及其实现

提出完全变光滑长度SPH法及其算法实现.方程组基于对称形式核函数近似,SPH密度演化方程与变光滑长度方程隐式关联;在Springel提出的全守恒SPH方程组基础上,通过将分散核近似形式改进为对称核近似形式得到SPH动量方程和能量方程.采用迭代求解密度演化方程和变光滑长度方程,显式求解SPH动量方程和能量方程,增加的计算量相对很少.给出三个1D激波管算例和2D Sedov算例验证方法的有效性.数值结果表明,算法保持动量和能量的守恒律,解决了传统SPH法中由于变光滑长度影响带来的计算误差,且在模拟2DSedov问题时能得到比Springel方法更准确的压强峰值位置和中心压强值.特别适合于模拟爆炸与冲击、大变形大扭曲等密度梯度和光滑长度梯度剧烈变化的问题.     

叶栅内定常及非定常粘性流动数值模拟

本文给出了一个模拟叶栅内准三维定常和非定常粘性流动的数值方法。对于定常流动,采用ＴＶＤ Ｌａｘ－Ｗｅｎｄｒｏｆｆ格式和代数湍流模型求解雷诺平均Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程,使用当地时间步长和多网格技术使计算加速收敛到定常状态;对于非定常流动,使用双时间步长和全隐式离散,采用与求解定常流动相似的多网格方法求解隐式离散方程。文中给出了ＶＫＩ透平叶栅内的定常流结果和１．５级透平叶栅内的非定常数值结果。     

YLSG 107高升力翼型主动流动控制的数值模拟

数值模拟零质量射流与YLSG 107翼型绕流的干扰流场,探讨零质量射流在高升力翼型失速控制中的控制效果、控制特性及控制机理.数值模拟以积分形式雷诺平均Navier-Stokes(N-S)方程为控制方程,采用格心有限体积法进行求解.通过在喷口上施加非定常边界吹/吸边界条件模拟射流对翼型绕流的干扰.采用与风洞实验相同的来流状态和控制参数进行数值模拟,得到与实验相吻合的结果.为进一步研究控制特性和控制规律、提出改进的实验方案,研究不同动量系数、不同射流偏角对控制效果的影响,并对法向射流和近切向射流进行较深入的比较.研究表明,先前的风洞实验对应的射流动量系数(0.000 014)偏小是控制效果不显著的重要原因之一,必须达到0.001以上才有明显控制效果(射流动量系数为0.005时可使该翼型失速迎角增大2°,最大升力提高8.7%);近切向射流在失速控制方面明显优于法向射流.     

耦合广义非线性薛定谔方程的相互作用表象龙格库塔算法及其误差分析

本文通过表象变换, 将耦合广义非线性薛定谔方程 (C-GNLSE) 变换成相互作用表象中的向量方程, 再利用向量形式的4阶龙格-库塔迭代格式, 建立了一种在频域内求解C-GNLSE的同步更新迭代算法. 通过将该向量形式的相互作用表象中的4阶龙格-库塔 (V-JH-RK4IP) 算法应用于高双折射光子晶体光纤中超连续谱产生的数值模拟, 验证了算法的有效性, 通过与现有其他典型算法的比较, 表明以V-JH-RK4IP算法求解C-GNLSE具有最高的计算精度和计算效率. 关键词： 耦合广义非线性薛定谔方程(C-GNLSE) 相互作用表象 4阶龙格-库塔算法 超连续谱产生     

A robust implicit multigrid method for RANS equations with two-equation turbulence models

The design of a new, truly robust multigrid framework for the solution of steady-state Reynolds-Averaged Navier–Stokes (RANS) equations with two-equation turbulence models is presented. While the mean-flow equations and the turbulence model equations are advanced in time in a loosely-coupled manner, their multigrid cycling is strongly coupled (FC-MG). Thanks to the loosely-coupled approach, the unconditionally positive-convergent implicit time-integration scheme for two-equation turbulence models (UPC) is used. An improvement to the basic UPC scheme convergence characteristics is developed and its extension within the multigrid method is proposed. The resulting novel FC-MG-UPC algorithm is nearly free of artificial stabilizing techniques, leading to increased multigrid efficiency. To demonstrate the robustness of the proposed algorithm, it is applied to linear and non-linear two-equation turbulence models. Numerical experiments are conducted, simulating separated flow about the NACA4412 airfoil, transonic flow about the RAE2822 airfoil and internal flow through a plane asymmetric diffuser. Results obtained from numerical simulations demonstrate the strong consistency and case-independence of the method.     

Conservative form of interpolated differential operator scheme for compressible and incompressible fluid dynamics

The proposed scheme, which is a conservative form of the interpolated differential operator scheme (IDO-CF), can provide high accurate solutions for both compressible and incompressible fluid equations. Spatial discretizations with fourth-order accuracy are derived from interpolation functions locally constructed by both cell-integrated values and point values. These values are coupled and time-integrated by solving fluid equations in the flux forms for the cell-integrated values and in the derivative forms for the point values. The IDO-CF scheme exactly conserves mass, momentum, and energy, retaining the high resolution more than the non-conservative form of the IDO scheme. A direct numerical simulation of turbulence is carried out with comparable accuracy to that of spectral methods. Benchmark tests of Riemann problems and lid-driven cavity flows show that the IDO-CF scheme is immensely promising in compressible and incompressible fluid dynamics studies.     

The effect of turbulence damping on acoustic wave propagation in tubes

The attenuation of sound due to the interaction between a low Mach number turbulent boundary layer and acoustic waves can be significant at low frequencies or in narrow tubes. In a recent publication by the present authors the acoustics of charge air coolers for passenger cars has been identified as an interesting application where turbulence attenuation can be of importance. Favourable low-frequency damping has been observed that could be used for control of the in-duct sound that is created by the engine gas exchange process. Analytical frequency-dependent models for the eddy viscosity that controls the momentum and thermal boundary layers are available but are restricted to thin acoustic boundary layers. For cases with cross-sections of a few millimetres a model based on thin acoustic boundary layers will not be applicable in the frequency range of interest.In the present paper a frequency-dependent axis-symmetric numerical model for interaction between turbulence and acoustic waves is proposed. A finite element scheme is used to formulate the time harmonic linearized convective equations for conservation of mass, momentum and energy into one coupled system of equations. The turbulence is introduced with a linear model for the eddy viscosity that is added to the shear viscosity. The proposed model is validated by comparison with experimental data from the literature.     

An Asymptotic-Preserving all-speed scheme for the Euler and Navier–Stokes equations

We present an Asymptotic-Preserving ‘all-speed’ scheme for the simulation of compressible flows valid at all Mach-numbers ranging from very small to order unity. The scheme is based on a semi-implicit discretization which treats the acoustic part implicitly and the convective and diffusive parts explicitly. This discretization, which is the key to the Asymptotic-Preserving property, provides a consistent approximation of both the hyperbolic compressible regime and the elliptic incompressible regime. The divergence-free condition on the velocity in the incompressible regime is respected, and an the pressure is computed via an elliptic equation resulting from a suitable combination of the momentum and energy equations. The implicit treatment of the acoustic part allows the time-step to be independent of the Mach number. The scheme is conservative and applies to steady or unsteady flows and to general equations of state. One and two-dimensional numerical results provide a validation of the Asymptotic-Preserving ‘all-speed’ properties.     

Iterative nonlinear beam propagation using Hamiltonian ray tracing and Wigner distribution function

We present an iterative method for simulating beam propagation in nonlinear media using Hamiltonian ray tracing. The Wigner distribution function of the input beam is computed at the entrance plane and is used as the initial condition for solving the Hamiltonian equations. Examples are given for the study of periodic self-focusing, spatial solitons, and Gaussian-Schell model in Kerr-effect media. Simulation results show good agreement with the split-step beam propagation method. The main advantage of ray tracing, even in the nonlinear case, is that ray diagrams are intuitive and easy to interpret in terms of traditional optical engineering terms, such as aberrations, ray-intercept plots, etc.     

Binary quantum turbulence arising from countersuperflow instability in two-component Bose-Einstein condensates

We theoretically study the development of quantum turbulence from two counter-propagating superfluids of miscible Bose-Einstein condensates by numerically solving the coupled Gross-Pitaevskii equations. When the relative velocity exceeds a critical value, the countersuperflow becomes unstable and quantized vortices are nucleated, which leads to isotropic quantum turbulence consisting of two superflows. It is shown that the binary turbulence can be realized experimentally in a trapped system.