流体力学分区算法及相关理论问题

分区计算是计算流体力学中处理复杂几何区域、多块物体相互运动和进行并行计算的一类重要方法.对分区算法、相关理论问题和一些数值现象进行了综述.     

分区拼接网格算法数值模拟超声速复杂流场

以分区拼接网格数值模拟方法为研究对象,对二阶Godunov方法在拼接网格中流场的数值模拟进行了研究,发展了适用于Godunov格式的通量守恒型算法,结合二阶Godunov有限体积法离散非定常Euler方程,数值模拟了捆绑火箭及航天飞机的超声速流场,计算结果正确描述了流场中的激波相交、反射等干扰特性.     

一种含运动固壁超声速流动的Descartes网格算法

提出一种Descartes网格算法,用于数值求解含任意复杂及运动固壁的超声速流动问题.采用位标集函数确定和跟踪流-固界面.引入虚网格技术处理流-固边界条件,并沿法向和切向分别进行计算.该算法简单、稳健,可与高阶有限差分格式并用.选取一组一维/二维静止或运动物体绕流算例,验证其有效性.     

超声速后台阶流动/射流相互作用的数值模拟

采用高精度格式求解二维Navier-Stokes方程研究超声速射流与同向超声速后台阶流动相互作用的流场基本结构及规律,分别应用5阶WENO格式、6阶中心差分格式离散对流项和黏性项,时间推进采用3阶Runge-Kutta格式,并应用消息传递接口（message passing interface,MPI）非阻塞式通信实现并行化.分别研究了超声速后台阶流动、超声速射流的基本结构特征,以此讨论和分析超声速后台阶流动/射流相互作用的特征,以及不同来流条件对波系结构、涡结构、剪切层、膨胀扇等的影响,尤其是来流剪切层和射流剪切层的相互作用,形成复杂的波系结构及相互干扰的流动现象.      

基于SIMPLE算法的时间离散格式比较

本文在具有高精度空间离散格式的SIMPLE算法计算代码的基础上研究比较了非稳态计算的时间离散格式.分别采用一阶全隐和二阶全隐格式对非稳态的控制方程进行离散,通过方腔驱动流和圆柱绕流两个经典算例的分析,比较了这两种格式在计算精确性和计算效率等方面的性能。     

微尺度串列双圆柱绕流的分子动力学并行模拟

本文采用基于空间分解算法的分子动力学并行模拟方法,研究了微尺度、低雷诺数（Re=40）下串列等大的双圆柱绕流现象。结果表明：随着间距比L＊／D＊的增加,流动存在3种特征状态：当L＊／D＊〈1．1时,同单一物体的绕钝体流动相似;当1．1〈L＊／0＊〈3．5时,涡脱落现象只在下游圆柱出现,在两圆柱之间有交替附着于下游圆柱的...     

旋转回流的一种多重交错网格计算方法

针对圆柱坐标系下原始变量Navier Stokes方程,在有限控制容积法和压力修正的基础上,引入多重交错网格算法及非线性方程的FAS全近似格式,并对封闭圆柱空腔内的旋转流动进行数值模拟.     

带副翼的翼身组合体绕流的Euler和N-S方程解

将对接分区网格与分区求解算法结合,有效地求解了带副翼偏转的翼身组合体绕流的N S方程.数值方法中选用VanLeer分裂格式离散无粘通量项,采用中心差分格式来离散粘性通量项.分区交界面采用了一种满足通量守恒的内边界耦合条件.数值算例表明该方法是求解带操纵面偏转的翼身组合体绕流的有效方法.     

物体穿越激波非定常流场的统一分区加强解

用数值方法求解非定常可压缩Navier Stokes方程,模拟了物体以超音速穿越与之同向运动激波及逆向运动激波的全过程.为确保数值模拟中的激波和物体的相对运动有足够的网格分辨率,使用了统一分区加强解算法,详细给出了追击穿越和碰撞穿越的流场结构和气动力,并且通过和单区结果的对比,说明统一分区加强解算法是一种处理复杂几何外形和提高网格分辨率的简单有效方法.     

多密度组合体重心研究

利用静力学原理,采用积分算法,求出了多密度物体重心坐标计算公式.并以装水圆柱桶为例,计算出了其重心坐标,给出水流出过程中重心函数变化曲线,求出了重心变化的极值点.     

Nested Cartesian grid method in incompressible viscous fluid flow

In this work, the local grid refinement procedure is focused by using a nested Cartesian grid formulation. The method is developed for simulating unsteady viscous incompressible flows with complex immersed boundaries. A finite-volume formulation based on globally second-order accurate central-difference schemes is adopted here in conjunction with a two-step fractional-step procedure. The key aspects that needed to be considered in developing such a nested grid solver are proper imposition of interface conditions on the nested-block boundaries, and accurate discretization of the governing equations in cells that are with block-interface as a control-surface. The interpolation procedure adopted in the study allows systematic development of a discretization scheme that preserves global second-order spatial accuracy of the underlying solver, and as a result high efficiency/accuracy nested grid discretization method is developed. Herein the proposed nested grid method has been widely tested through effective simulation of four different classes of unsteady incompressible viscous flows, thereby demonstrating its performance in the solution of various complex flow–structure interactions. The numerical examples include a lid-driven cavity flow and Pearson vortex problems, flow past a circular cylinder symmetrically installed in a channel, flow past an elliptic cylinder at an angle of attack, and flow past two tandem circular cylinders of unequal diameters. For the numerical simulations of flows past bluff bodies an immersed boundary (IB) method has been implemented in which the solid object is represented by a distributed body force in the Navier–Stokes equations. The main advantages of the implemented immersed boundary method are that the simulations could be performed on a regular Cartesian grid and applied to multiple nested-block (Cartesian) structured grids without any difficulty. Through the numerical experiments the strength of the solver in effectively/accurately simulating various complex flows past different forms of immersed boundaries is extensively demonstrated, in which the nested Cartesian grid method was suitably combined together with the fractional-step algorithm to speed up the solution procedure.     

An immersed boundary method based on discrete stream function formulation for two- and three-dimensional incompressible flows

An immersed boundary method is proposed in the framework of discrete stream function formulation for incompressible flows. In order to impose the non-slip boundary condition, the forcing term is determined implicitly by solving a linear system. The number of unknowns of the linear system is the same as that of the Lagrangian points representing the body surface. Thus the extra cost in force calculation is negligible if compared with that in the basic flow solver. In order to handle three-dimensional flows at moderate Reynolds numbers, a parallelized flow solver based on the present method is developed using the domain decomposition strategy. To verify the accuracy of the immersed-boundary method proposed in this work, flow problems of different complexity (decaying vortices, flows over stationary and oscillating cylinders and a stationary sphere, and flow over low-aspect-ratio flat-plate) are simulated and the results are in good agreement with the experimental or computational data in previously published literatures.     

基于OpenFOAM的投影法磁流体求解器开发与验证

在开源计算流体力学C++工具包OpenFOAM环境下开发了低磁雷诺数条件下的磁流体求解器,并进行了验证。采用投影算法求解动量方程和压力泊松方程;采用非结构网格同位相容守恒算法求解电势泊松方程、感应电流和洛伦兹力;采用边界耦合方法求解流固耦合电势场。通过对均匀磁场下导电方管和导电圆管内的完全发展磁流体层流的数值模拟和解析解的对比,对求解器进行了验证。进一步对非均匀强磁场作用下导电方管和导电圆管内完全发展磁流体层流进行了数值模拟,并与ALEX实验结果进行了比较。数值解和实验结果吻合良好。所开发的求解器可用于复杂结构强磁场作用下磁流体的数值模拟研究。     

Filtered density function simulator on unstructured meshes

A new computational filtered density function (FDF) methodology is developed for large eddy simulation (LES) of turbulent reacting flows. This methodology is based on a Lagrangian Monte Carlo (MC) FDF solver constructed on a domain portrayed by an unstructured mesh. The base filtered transport equations on this mesh are solved by a finite-volume (FV) method. The consistency of the hybrid FV-MC solver and the realizability of the simulated results are demonstrated via LES of a temporally developing mixing layer. The overall performance of the model is appraised by comparison with direct numerical simulation (DNS) data. The algorithmic implementation in the commercial software ANSYS-FLUENT facilitates future FDF-LES of turbulent combustion in complex configurations.     

Modelling atmospheric flows with adaptive moving meshes

An anelastic atmospheric flow solver has been developed that combines semi-implicit non-oscillatory forward-in-time numerics with a solution-adaptive mesh capability. A key feature of the solver is the unification of a mesh adaptation apparatus, based on moving mesh partial differential equations (PDEs), with the rigorous formulation of the governing anelastic PDEs in generalised time-dependent curvilinear coordinates. The solver development includes an enhancement of the flux-form multidimensional positive definite advection transport algorithm (MPDATA) — employed in the integration of the underlying anelastic PDEs — that ensures full compatibility with mass continuity under moving meshes. In addition, to satisfy the geometric conservation law (GCL) tensor identity under general moving meshes, a diagnostic approach is proposed based on the treatment of the GCL as an elliptic problem. The benefits of the solution-adaptive moving mesh technique for the simulation of multiscale atmospheric flows are demonstrated. The developed solver is verified for two idealised flow problems with distinct levels of complexity: passive scalar advection in a prescribed deformational flow, and the life cycle of a large-scale atmospheric baroclinic wave instability showing fine-scale phenomena of fronts and internal gravity waves.     

An efficient multi-scale Poisson solver for the incompressible Navier–Stokes equations with immersed boundaries

The iterative-multi-scale-finite-volume (IMSFV) procedure is applied as an efficient solver for the pressure Poisson equation arising in numerical methods for the simulation of incompressible flows with the immersed-interface method (IIM). Motivated by the requirements of the specific IIM implementation, a modified version of the IMSFV algorithm is presented to allow the solution of problems, in which the varying coefficient of the elliptic equation (e.g. the permeability of the medium in the context of the simulation of flows in porous media) varies over several orders of magnitude or even becomes zero within the integration domain. Furthermore, a strategy is proposed to incorporate the iterative procedure needed by the IIM to converge out constraints at immersed boundaries into the iterative IMSFV cycle. No significant deterioration of performance of the IMSFV method is observed with respect to cases, in which no iterative improvement of the boundary conditions is considered.     

A new solver for granular avalanche simulation: Indoor experiment verification and field scale case study

A new solver based on the high-resolution scheme with novel treatments of source terms and interface capture for the Savage-Hutter model is developed to simulate granular avalanche flows. The capability to simulate flow spread and deposit processes is verified through indoor experiments of a two-dimensional granular avalanche. Parameter studies show that reduction in bed friction enhances runout efficiency, and that lower earth pressure restraints enlarge the deposit spread. The April 9, 2000, Yigong avalanche in Tibet, China, is simulated as a case study by this new solver. The predicted results, including evolution process, deposit spread, and hazard impacts, generally agree with site observations. It is concluded that the new solver for the Savage-Hutter equation provides a comprehensive software platform for granular avalanche simulation at both experimental and field scales. In particular, the solver can be a valuable tool for providing necessary information for hazard forecasts, disaster mitigation, and countermeasure decisions in mountainous areas.     

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.     

A high-order multi-dimensional HLL-Riemann solver for non-linear Euler equations

This article presents a numerical model that enables to solve on unstructured triangular meshes and with a high-order of accuracy, a multi-dimensional Riemann problem that appears when solving hyperbolic problems.For this purpose, we use a MUSCL-like procedure in a “cell-vertex” finite-volume framework. In the first part of this procedure, we devise a four-state bi-dimensional HLL solver (HLL-2D). This solver is based upon the Riemann problem generated at the centre of gravity of a triangular cell, from surrounding cell-averages. A new three-wave model makes it possible to solve this problem, approximately. A first-order version of the bi-dimensional Riemann solver is then generated for discretizing the full compressible Euler equations.In the second part of the MUSCL procedure, we develop a polynomial reconstruction that uses all the surrounding numerical data of a given point, to give at best third-order accuracy. The resulting over determined system is solved by using a least-square methodology. To enforce monotonicity conditions into the polynomial interpolation, we develop a simplified central WENO (CWENO) procedure.Numerical tests and comparisons with competing numerical methods enable to identify the salient features of the whole model.     

Development of Lattice Boltzmann Flux Solver for Simulation of Incompressible Flows

A lattice Boltzmann flux solver (LBFS) is presented in this work for simulation of incompressible viscous and inviscid flows. The new solver is based on Chapman-Enskog expansion analysis, which is the bridge to link Navier-Stokes (N-S) equations and lattice Boltzmann equation (LBE). The macroscopic differential equations are discretized by the finite volume method, where the flux at the cell interface is evaluated by local reconstruction of lattice Boltzmann solution from macroscopic flow variables at cell centers. The new solver removes the drawbacks of conventional lattice Boltzmann method such as limitation to uniform mesh, tie-up of mesh spacing and time interval, limitation to viscous flows. LBFS is validated by its application to simulate the viscous decaying vortex flow, the driven cavity flow, the viscous flow past a circular cylinder, and the inviscid flow past a circular cylinder. The obtained numerical results compare very well with available data in the literature, which show that LBFS has the second order of accuracy in space, and can be well applied to viscous and inviscid flow problems with non-uniform mesh and curved boundary.