二维非结构网格有限体积格式的壁热现象和人工热流方法

基于Noh人工热流方法的思想,推导出任意多边形网格人工热流的格式,将Noh的结构网格人工热流方法推广到任意多边形非结构网格的拉氏有限体积方法中,抑制或消除了壁热现象.几个数值实验,包括Noh激波反射问题和两块钢金属碰撞问题,验证了方法的可行性.     

一维理想弹塑性流体的数值模拟及壁热现象的抑制

针对一维理想弹塑性流体的Wilkins模型，提出HLLC型的近似黎曼解法器.该解法器引入塑性波并保证波的个数和真实的物理情况一致，其中波速选取由波系的特征分析来确定.整个算法实施简单，无需迭代，为减少强冲击（稀疏）模拟时的壁热误差，设计相应的壁热粘性，有效地抑制了非物理的壁热现象.     

基于变分原理的二维热传导方程差分格式

研究二维热传导方程的差分数值模拟.用变分原理在不规则结构网格上建立热流通量形式的差分格式.将热流通量作为未知函数求泛函极值,并与温度函数联立求解.克服通常九点格式用插值方法计算网格边界上的热传导系数和网格结点上的温度所引入的误差.     

不同精度格式的格子Boltzmann热模型的传热分析

通过在格子Boltzmann (LBM)热模型中添加参数项,使得在对应的宏观传热方程中,消除了一阶非线性误差项,具备二阶精度.通过Rayleigh-Benard对流数值试算,初步探索该二阶精度格式及其对应的一阶精度格式三个热模型的传热特征和适应性,并做出相应对比分析.针对一二阶精度模型在Ra数极高或热传导系数极大时,Nu数的计算与经验值相比出现较大偏差,分析LBM对应宏观热传导方程的截断误差后,在平衡分布函数中引进一个调节因子.通过调节对应宏观传热方程的截断误差项系数,校正Nu数的计算偏差,提高模拟精度,拓展模拟范围,增强了LBM作为一个数值方法在传热中的适应性.     

应用有限差分进行涡轮气热弹耦合计算关键性问题的研究

分析了涡轮发动机涡轮内的气热弹耦合现象,并讨论了气热弹耦合研究方向及基本步骤,为了采用有限差分进行涡轮气热弹耦合数值仿真,详细地讨论了该方法的关键性问题:模型方程、边界条件、差分格式、结构化网格块群、固体应力场的求解以及数值仿真平台的建立.最后通过几个简单的算例初步验证了该方法的可行性和可靠性.     

地上换热对岩土热响应测试的影响分析

建立了岩土热响应测试的三维数值分析模型,分析了不同情况下地上换热对岩土热响应测试结果的影响。结果表明;在相同条件下,埋管深度越小,地上换热引起岩土热响应测试结果的误差就越大;测试环境空气温度的波动会导致岩土热响应测试结果的波动;在实际工程应用中,可通过抑制载流体与地上环境之间一切形式的换热,来减少岩土热响应测试结果的误差。     

热声波数值模拟的虚假振荡研究

采用可压缩流动的SIMPLE算法对一维封闭空腔内由边界突然加热所引起的非稳态热声波进行了数值模拟,对流-扩散项采用了中心差分、一阶迎风差分、QUICK、及MIJSCL等不同格式。计算表明各种格式均存在不同程度的虚假振荡现象,其大小与热声波的强度及离散格式的形式等多种因素有关。这些结果对热声波的进一步研究及高效可靠的对流差分格式的开发具有重要意义。     

碳遮光石英气凝胶传热机制与热性能数值模拟

建立了碳遮光石英气凝胶传热机制及热性能数值模拟方法,在交叉立方阵列导热模型、热辐射传输谱带模型、辐射导热耦合传热模型基础上,采用蒙特卡罗方法与有限体积法数值模拟了气凝胶内的热辐射传输及辐射导热耦合传热,并以表观导热系数描述气凝胶传热性能.以某石英气凝胶为例,定量模拟了热性能、各种传热方式的作用及温度依赖性,分析了应用Rosseland扩散近似引起的误差.     

湍流边界层近壁区多个相干结构的数值模拟

从流动稳定性理论中的一般共振三波概念出发,提出一种湍流边界层近壁区多个相干结构的理论模型,采用高精度差分格式和Fourier谱展开相结合的方法,求解三维不可压Navier-Stokes方程,直接数值模拟近壁湍流多个相干结构的演化问题.并将得到的湍流边界层近壁区多个相干结构的数值演化特性与实验观察到的特性进行了比较.     

抛物方程显隐修正迎风区域分裂差分法

给出抛物方程一种有效的区域分裂差分格式,提高了计算效率.对一阶项采用二阶迎风差分格式,内边界点和各子区域分别采用显隐差分格式.在较弱的稳定性条件下,得到离散l2模误差估计结果.最后给出具体的数值算例,以验证方法的实用性.     

求解多维欧拉方程的二阶非结构网格混合旋转Riemann求解器

将基于旋转近似Riemann求解器的二阶精度迎风型有限体积方法推广到非结构网格,采用基于网格中心的有限体积法,梯度的计算采用基于节点的方法引入更多的控制体模板,限制器的构造采用与非结构化网格相适应的形式.在求解Riemann问题时,沿具有一定物理意义的两个迎风方向,即控制体界面两侧速度差矢量方向及与之正交的方向.能够完全消除基于Riemann求解器的通量差分裂格式存在的激波不稳定或"红斑"现象.为减小计算量,采用HLL和Roe FDS混合旋转格式.     

Multidimensional HLLE Riemann solver: Application to Euler and magnetohydrodynamic flows

In this work we present a general strategy for constructing multidimensional HLLE Riemann solvers, with particular attention paid to detailing the two-dimensional HLLE Riemann solver. This is accomplished by introducing a constant resolved state between the states being considered, which introduces sufficient dissipation for systems of conservation laws. Closed form expressions for the resolved fluxes are also provided to facilitate numerical implementation. The Riemann solver is proved to be positively conservative for the density variable; the positivity of the pressure variable has been demonstrated for Euler flows when the divergence in the fluid velocities is suitably restricted so as to prevent the formation of cavitation in the flow.We also focus on the construction of multidimensionally upwinded electric fields for divergence-free magnetohydrodynamical (MHD) flows. A robust and efficient second order accurate numerical scheme for two and three-dimensional Euler and MHD flows is presented. The scheme is built on the current multidimensional Riemann solver and has been implemented in the author’s RIEMANN code. The number of zones updated per second by this scheme on a modern processor is shown to be cost-competitive with schemes that are based on a one-dimensional Riemann solver. However, the present scheme permits larger timesteps.Accuracy analysis for multidimensional Euler and MHD problems shows that the scheme meets its design accuracy. Several stringent test problems involving Euler and MHD flows are also presented and the scheme is shown to perform robustly on all of them.     

High-order finite-volume methods for the shallow-water equations on the sphere

This paper presents a third-order and fourth-order finite-volume method for solving the shallow-water equations on a non-orthogonal equiangular cubed-sphere grid. Such a grid is built upon an inflated cube placed inside a sphere and provides an almost uniform grid point distribution. The numerical schemes are based on a high-order variant of the Monotone Upstream-centered Schemes for Conservation Laws (MUSCL) pioneered by van Leer. In each cell the reconstructed left and right states are either obtained via a dimension-split piecewise-parabolic method or a piecewise-cubic reconstruction. The reconstructed states then serve as input to an approximate Riemann solver that determines the numerical fluxes at two Gaussian quadrature points along the cell boundary. The use of multiple quadrature points renders the resulting flux high-order. Three types of approximate Riemann solvers are compared, including the widely used solver of Rusanov, the solver of Roe and the new AUSM⁺-up solver of Liou that has been designed for low-Mach number flows. Spatial discretizations are paired with either a third-order or fourth-order total-variation-diminishing Runge–Kutta timestepping scheme to match the order of the spatial discretization. The numerical schemes are evaluated with several standard shallow-water test cases that emphasize accuracy and conservation properties. These tests show that the AUSM⁺-up flux provides the best overall accuracy, followed closely by the Roe solver. The Rusanov flux, with its simplicity, provides significantly larger errors by comparison. A brief discussion on extending the method to arbitrary order-of-accuracy is included.     

A high-order cell-centered Lagrangian scheme for compressible fluid flows in two-dimensional cylindrical geometry

The goal of this paper is to present high-order cell-centered schemes for solving the equations of Lagrangian gas dynamics written in cylindrical geometry. A node-based discretization of the numerical fluxes is obtained through the computation of the time rate of change of the cell volume. It allows to derive finite volume numerical schemes that are compatible with the geometric conservation law (GCL). Two discretizations of the momentum equations are proposed depending on the form of the discrete gradient operator. The first one corresponds to the control volume scheme while the second one corresponds to the so-called area-weighted scheme. Both formulations share the same discretization for the total energy equation. In both schemes, fluxes are computed using the same nodal solver which can be viewed as a two-dimensional extension of an approximate Riemann solver. The control volume scheme is conservative for momentum, total energy and satisfies a local entropy inequality in its first-order semi-discrete form. However, it does not preserve spherical symmetry. On the other hand, the area-weighted scheme is conservative for total energy and preserves spherical symmetry for one-dimensional spherical flow on equi-angular polar grid. The two-dimensional high-order extensions of these two schemes are constructed employing the generalized Riemann problem (GRP) in the acoustic approximation. Many numerical tests are presented in order to assess these new schemes. The results obtained for various representative configurations of one and two-dimensional compressible fluid flows show the robustness and the accuracy of our new schemes.     

An optimized spectral difference scheme for CAA problems

In the implementation of spectral difference (SD) method, the conserved variables at the flux points are calculated from the solution points using extrapolation or interpolation schemes. The errors incurred in using extrapolation and interpolation would result in instability. On the other hand, the difference between the left and right conserved variables at the edge interface will introduce dissipation to the SD method when applying a Riemann solver to compute the flux at the element interface. In this paper, an optimization of the extrapolation and interpolation schemes for the fourth order SD method on quadrilateral element is carried out in the wavenumber space through minimizing their dispersion error over a selected band of wavenumbers. The optimized coefficients of the extrapolation and interpolation are presented. And the dispersion error of the original and optimized schemes is plotted and compared. An improvement of the dispersion error over the resolvable wavenumber range of SD method is obtained. The stability of the optimized fourth order SD scheme is analyzed. It is found that the stability of the 4th order scheme with Chebyshev–Gauss–Lobatto flux points, which is originally weakly unstable, has been improved through the optimization. The weak instability is eliminated completely if an additional second order filter is applied on selected flux points. One and two dimensional linear wave propagation analyses are carried out for the optimized scheme. It is found that in the resolvable wavenumber range the new SD scheme is less dispersive and less dissipative than the original scheme, and the new scheme is less anisotropic for 2D wave propagation. The optimized SD solver is validated with four computational aeroacoustics (CAA) workshop benchmark problems. The numerical results with optimized schemes agree much better with the analytical data than those with the original schemes.     

On the role of Riemann solvers in Discontinuous Galerkin methods for magnetohydrodynamics

It has been claimed that the particular numerical flux used in Runge–Kutta Discontinuous Galerkin (RKDG) methods does not have a significant effect on the results of high-order simulations. We investigate this claim for the case of compressible ideal magnetohydrodynamics (MHD). We also address the role of limiting in RKDG methods.For smooth nonlinear solutions, we find that the use of a more accurate Riemann solver in third-order simulations results in lower errors and more rapid convergence. However, in the corresponding fourth-order simulations we find that varying the Riemann solver has a negligible effect on the solutions.In the vicinity of discontinuities, we find that high-order RKDG methods behave in a similar manner to the second-order method due to the use of a piecewise linear limiter. Thus, for solutions dominated by discontinuities, the choice of Riemann solver in a high-order method has similar significance to that in a second-order method. Our analysis of second-order methods indicates that the choice of Riemann solver is highly significant, with the more accurate Riemann solvers having the lowest computational effort required to obtain a given accuracy. This allows the error in fourth-order simulations of a discontinuous solution to be mitigated through the use of a more accurate Riemann solver.We demonstrate the minmod limiter is unsuitable for use in a high-order RKDG method. It tends to restrict the polynomial order of the trial space, and hence the order of accuracy of the method, even when this is not needed to maintain the TVD property of the scheme.     

A new class of gas-kinetic relaxation schemes for the compressible Euler equations

Starting from the gas-kinetic model, a new class of relaxation schemes for the Euler equations is presented. In contrast to the Riemann solver, these schemes provide a multidimensional dynamical gas evolution model, which combines both Lax-Wendroff and kinetic flux vector splitting schemes, and their coupling is based on the fact that a nonequilibrium state will evolve into an equilibrium state along with the increase of entropy. The numerical fluxes are constructed without getting into the details of the particle collisions. The results for many well-defined test cases are presented to indicate the robustness and accuracy of the current scheme.     

HLLC-type Riemann solver for the Baer–Nunziato equations of compressible two-phase flow

We first construct an approximate Riemann solver of the HLLC-type for the Baer–Nunziato equations of compressible two-phase flow for the “subsonic” wave configuration. The solver is fully nonlinear. It is also complete, that is, it contains all the characteristic fields present in the exact solution of the Riemann problem. In particular, stationary contact waves are resolved exactly. We then implement and test a new upwind variant of the path-conservative approach; such schemes are suitable for solving numerically nonconservative systems. Finally, we use locally the new HLLC solver for the Baer–Nunziato equations in the framework of finite volume, discontinuous Galerkin finite element and path-conservative schemes. We systematically assess the solver on a series of carefully chosen test problems.     

A positive MUSCL-Hancock scheme for ideal magnetohydrodynamics

We present a highly robust second order accurate scheme for the Euler equations and the ideal MHD equations. The scheme is of predictor–corrector type, with a MUSCL scheme following as a special case. The crucial ingredients are an entropy stable approximate Riemann solver and a new spatial reconstruction that ensures positivity of mass density and pressure. For multidimensional MHD, a new discrete form of the Powell source terms is vital to ensure the stability properties. The numerical examples show that the scheme has superior stability compared to standard schemes, while maintaining accuracy. In particular, the method can handle very low values of pressure (i.e. low plasma β$\beta$ or high Mach numbers) and low mass densities.