基于分子平均自由程的热流计算壁面网格准则

面向壁面热流Navier-Stokes方程数值计算的应用, 研究了壁面法向网格尺度问题.从气体分子物理学角度, 对壁面热流进行了微观统计分析,指出壁面热流与气体分子平均自由程之间的内在联系是壁面法向网格尺度存在最优下限的原因. 据此建立了一个新的壁面法向网格尺度准则, 即MFP准则.该准则简单、实用, 具有清晰的热力学物理意义, 且只依赖于壁面局部参数. 并以多种对比算例验证了MFP准则的准确性和适用性.      

AN ANALYSIS OF TRANSIENT HEAT CONDUCTION FOR FUNCTIONALLY GRADED SPHERES BASED ON STATE SPACE THEORY1)

基于状态空间理论研究功能梯度圆球的球对称瞬态热传导问题。根据热传导方程和热流密度的定义,取温度场和热流密度为系统的状态向量,通过将圆球分层和在时域内应用差分格式对控制方程进行离散,建立了系统的状态方程,给出了功能梯度圆球瞬态热传导问题的半解析解。算例分析表明：本文解不但结果正确、计算效率高,而且适用于材料参数沿径向任意梯度变化的圆球瞬态热传导分析。     

燃烧室两相流场亚网格燃烧模型的研究

在三维任意曲线坐标系下采用不同的亚网格燃烧模型对环形燃烧室火焰筒气液两相湍流瞬态反应流进行大涡模拟．计算中所采用的数学度模型有：k方程亚网格尺度模型估算亚网格湍流黏性；热通量辐射模型估算辐射换热，分别采用亚网格EBU燃烧模型（E-A model）、亚网格二阶矩输运方程模型（SOM）和亚网格二阶矩代数模型（SOM-A）估算化学反应速率．并在非交错网格系统下气相采用SIMPLE算法和混合差分格式求解，液相采用Lagrange处理，并用PSIC算法对其进行求解．通过实验结果和计算结果的比较，表明在三维任意曲线坐标系下对燃烧室火焰简两相湍流油雾燃烧流场进行大涡模拟，3种不同的亚网格燃烧模型都能真实反映两相湍流化学反应流流动及实际燃烧过程，而采用亚网格二阶矩输运方程模型稍优于其他两种亚网格燃烧模型．     

The effects of control‐volume distributed multi‐point flux approximation (CVD‐MPFA) on upscaling—A study using the CVD‐MPFA schemes

A family of flux‐continuous, locally conservative, finite‐volume schemes has been developed for solving the general geometry‐permeability tensor (petroleum reservoir‐simulation) pressure equation on structured and unstructured grids and are control‐volume distributed (textit Comput. Geo. 1998; 2 :259–290; Comput. Geo. 2002; 6 :433–452). The schemes are applicable to diagonal and full tensor pressure equation with generally discontinuous coefficients and remove the O(1) errors introduced by standard reservoir‐simulation schemes (two‐point flux approximation) when applied to full tensor flow approximation. The family of flux‐continuous schemes is quantified by a quadrature parameterization (Int. J. Numer. Meth. Fluids 2006; 51 :1177–1203). Improved convergence (for two‐ and three‐dimensional formulation) using the quadrature parameterization has been observed for the family of flux‐continuous control‐volume distributed multi‐point flux approximation (CVD‐MPFA) schemes (Ph.D. Thesis, University of Wales, Swansea, U.K., 2007). In this paper family of flux‐continuous (CVD‐MPFA) schemes are used as a part of numerical upscaling procedure for upscaling the fine‐scale grid information (permeability) onto a coarse grid scale. A series of data‐sets (SPE, 2001) are tested where the upscaled permeability tensor is computed on a sequence of grid levels using the same fixed range of quadrature points in each case. The refinement studies presented involve:

• (i) Refinement comparison study: In this study, permeability distribution for cells at each grid level is obtained by upscaling directly from the fine‐scale permeability field as in standard simulation practice.
• (ii) Refinement study with renormalized permeability: In this refinement comparison, the local permeability is upscaled to the next grid level hierarchically, so that permeability values are renormalized to each coarser level. Hence, showing only the effect of increased grid resolution on upscaled permeability, compared with that obtained directly from the fine‐scale solution.
• (iii) Refinement study with invariant permeability distribution: In this study, a classical mathematical convergence test is performed. The same coarse‐scale underlying permeability map is preserved on all grid levels including the fine‐scale reference solution.

The study is carried out for the discretization of the scheme in physical space. The benefit of using specific quadrature points is demonstrated for upscaling in this study and superconvergence is observed. Copyright © 2010 John Wiley & Sons, Ltd.     

多维迎风方法在非结构/混合网格热流计算中的应用研究

非结构/混合网格具有极强的几何灵活性,在复杂外形飞行器的气动力特性数值模拟中已得到广泛应用,但目前还难以准确地预测气动热环境。本文从非结构/混合网格热流计算的三个需求出发,选取了多维迎风方法,并与其他方法进行了对比研究。以二维圆柱高超声速绕流这一Benchmark典型问题为例,对比研究了多维迎风方法和几种广泛使用的无粘通量格式(Roe格式、Van Leer格式和AUSMDV格式)对混合网格热流计算精度的影响。结果表明,多维迎风方法在热流计算精度、鲁棒性以及收敛性方面表现良好。最后,将多维迎风方法应用于常规混合网格上的圆柱和钝双锥绕流问题,均得到了较好的热流计算结果,为非结构/混合网格热流计算在复杂高超飞行器中的应用奠定了基础。     

An adaptive multilevel wavelet framework for scale‐selective WENO reconstruction schemes

We put forth a dynamic computing framework for scale‐selective adaptation of weighted essential nonoscillatory (WENO) schemes for the simulation of hyperbolic conservation laws exhibiting strong discontinuities. A multilevel wavelet‐based multiresolution procedure, embedded in a conservative finite volume formulation, is used for a twofold purpose. (i) a dynamic grid adaptation of the solution field for redistributing grid points optimally (in some sense) according to the underlying flow structures, and (ii) a dynamic minimization of the in built artificial dissipation of WENO schemes. Taking advantage of the structure detection properties of this multiresolution algorithm, the nonlinear weights of the conventional WENO implementation are selectively modified to ensure lower dissipation in smoother areas. This modification is implemented through a linear transition from the fifth‐order upwind stencil at the coarsest regions of the adaptive grid to a fully nonlinear fifth‐order WENO scheme at areas of high irregularity. Therefore, our computing algorithm consists of a dynamic grid adaptation strategy, a scale‐selective state reconstruction, a conservative flux calculation, and a total variation diminishing Runge‐Kutta scheme for time advancement. Results are presented for canonical examples drawn from the inviscid Burgers, shallow water, Euler, and magnetohydrodynamic equations. Our findings represent a novel direction for providing a scale‐selective dissipation process without a compromise on shock capturing behavior for conservation laws, which would be a strong contender for dynamic implicit large eddy simulation approaches.     

THE NUMERICAL STUDY OF LATTICE BOLTZMANN METHOD BASED ON DIFFERENT GRID STRUCTURE

传统的格子波尔兹曼方法(lattice-Boltzmann method, LBM)通常基于标准均匀网格, 这主要取决于速度的空 间离散格式.均匀网格结构的特点, 使LBM在处理具有复杂边界的问题时遇到较大的困难, 从而限制了它的应用.另外, 对于较为复杂的流动, 其流场存在流动变化剧烈和平缓的区域, 在流动变化剧烈的区域, 往往需要足够的网格点才能更好地捕捉到流场信息, 而均匀网格会使得网格数量过多, 这会增加计算量, 但网格数量过少又无法获得必要的流场信息, 使LBM的计算效率降低.为了解决上述问题, 用不同的网格结构, 以顶盖驱动的腔体内流、柱体绕流和翼型绕流为例, 探讨了提高LBM算法的计算效率和适用性问题.     

THE OPTIMAL SOLUTION FOR THE TWO-DIMENSIONAL CONVECTION-DIFFUSION MODEL IN BATCHING PRODUCT PIPELINE

将一种简化的自适应网格应用于顺序输送二维混油模型的计算中,网格的疏密程度由速度和浓度梯度控制且用动网格取代传统的固定网格。混油模型中采用混掺长度理论计算速度场和扩散系数,并考虑了油品黏度变化对混油的影响。提高传统数值解法的精度,减少计算量,使之适应于长输管线顺序输送的混油数值计算。通过在计算机上的数值模拟实验,证明了该方法的有效性,分析了雷诺数、输送顺序对混油的影响并得到混油随输送距离的变化规律。     

Conduction–radiation interaction in 3D irregular enclosures using the finite volume method

Interaction of conduction and radiation in 3D enclosures is carried out with a gray participating media. Application of block structured grid is shown with the finite volume method (FVM). Radiation modeling is performed with the FVM and is coupled with an ‘in-house’ code to solve the set of transport equations. The detailed numerical results are presented for a cubical and a cylindrical enclosure as these results are not available in the literature. The numerical simulation for the cylindrical enclosure is performed using a block-structured ‘O’ grid. Two additional geometries are considered in order to show the applicability of the present work. Results of temperature, radiative heat flux and total heat flux distributions are presented for different optical thicknesses, scattering albedoes, emissivities and conduction–radiation parameters. The 3D results are validated with the available 2D results or results with pure radiation problems as limiting cases.     

THE PROGRESS OF THE OVERLAPPING GRID TECHNIQUES

从网格装配和插值计算两个主要方面对现有的重叠网格方法进行了综述。首先,从挖洞方法和建立嵌入网格关系环节的寻点技术出发归纳和介绍了网格装配方法;其次,介绍了数值迭代过程中的插值计算方法,并特别讨论了插值守恒性以及插值计算精度等问题;另外,对重叠网格方法的并行计算和应用成果也作了介绍;最后,通过总结认为重叠网格方法在改进网格装配方法、改善插值和并行计算效率等方面仍需进一步研究。     

A VARIATIONAL INEQUALITY APPROACH FOR NON-STEADY SEEPAGE FLOW THROUGH THREE-DIMENSIONAL FRACTURE NETWORK

为了求解裂隙岩体有自由面非稳定渗流问题,将Darcy定律延拓至整个研究区域,使得潜在溢出边界条件满足Signorini型边界条件,建立了三维裂隙网络非稳定渗流问题的抛物型变分不等式（parabolic variational inequality,PVI）提法,并证明其与偏微分方程（partial differential equation,PDE）提法的等价性,从而将自由面上的流量条件以及潜在溢出边界上的互补条件转化成自然边界条件,降低该问题求解难度。同时给出了基于PVI提法的有限元数值求解方法,通过与交叉裂隙模型理论解的对比分析,证明了该方法的正确性。最后将该方法对含复杂三维裂隙网络的边坡进行非稳定渗流分析,计算结果表明该方法对于复杂裂隙网络求解具有较强的可靠性和适应性。     

STUDY ON ANALYSIS METHOD OF “THREE BOX” IN HEAT AND MASS TRANSFER PROCESS OF RESERVOIR MEDIA

油藏多孔介质孔隙组成及结构变化多样,一些特性参数很难全部获得,精确描述和分析困难;另外,多孔介质内渗流过程水力条件和作用机理复杂,存在热流固耦合作用,目前的一些分析方法和研究模型具有一定的局限性.提出了油藏多孔介质的表征单元体（representative elementary volume,REV）描述表征方法;基于表征单元体建立了多孔介质的黑箱模型、灰箱模型和白箱模型,据此提出了多孔介质的“黑箱→灰箱→白箱”分析过程.基于黑箱模型和灰箱模型推导了REV导热系数计算公式、给出了REV热质传递过程的热平衡方程.结合中国油藏热采情况,对多孔介质导热系数变化规律和蒸汽驱热质传递特性进行了分析,得到了一些有意义的结果.该工作为多孔介质热质传递过程分析提供了新思路和新方法.     

VISCOUS FLOWFIELD BASED ON DISCONTINUOUS BOUNDARY ELEMENT METHOD AND VORTEX METHOD

基于非协调边界元方法和涡方法的联合应用, 模拟了二维和三维黏性不可压缩流场. 计算中利用离散涡元对漩涡的产生、凝聚和输送过程进行模拟, 并将整体计算域分解为采用涡泡模拟的内部区域和用涡列模拟的数字边界层区域. 计算域中涡量场的拉伸和对流由Lagrangian涡方法模拟, 用随机走步模拟涡量场的扩散. 内部区域涡元涡量场速度由广义Biot-Savart公式计算, 势流场速度则采用非协调边界元方法计算. 非协调边界元将所有节点均取在光滑边界处, 从而避免了法向速度的不连续现象; 而对于系数矩阵不对称的大型边界元方程组,引入了非常高效的预处理循环型广义极小残余(the generalized minimum residual, GMRES)迭代算法, 使得边界元法的优势得到了充分发挥, 同时, 在内部涡元势流场计算中对近边界点采用了正则化算法, 该算法将奇异积分转化为沿单元围道上一系列线积分, 消除了势流计算中速度及速度梯度的奇异性. 二维、三维流场算例证明了所用方法的正确性, 也验证了该算法可以大幅度提高模拟精度和效率.     

OVERLAPPING GRIDS FOR THE SIMULATION OF PARTICLE INTERACTION AT THE MICRO SCALE

This paper describes the use of overlapping grids for the calculation of flow around single and multipleparticle configurations at the micro scale. The basic equations for calculation are those for conservation of mass and momentum which are solved using a common Finite-Volume formulation. The hydrodynamic particle-particle and particle-wall interaction can be calculated by using an overlapping or Chimera grid scheme. With the grid structuring procedure it is possible to use simple and structured grids around the particles and the overall main grid geometry. The particle grids are lapped over the main grid such that they can move independently after each time step without remeshing the whole geometry. The paper gives results for the validation of the code developed for general test cases, for a rotating ellipsoid in simple shear flow, the flow around particles attached to a wall, the motion of a particle in the vicinity of a wall and some results for the flow through a packed bed configuration.     

A sliding interface method for unsteady unstructured flow simulations

The objective of this work is to develop a sliding interface method for simulations involving relative grid motion that is fast and efficient and involves no grid deformation, remeshing, or hole cutting. The method is implemented into a parallel, node‐centred finite volume, unstructured viscous flow solver. The rotational motion is accomplished by rigidly rotating the subdomain representing the moving component. At the subdomain interface boundary, the faces along the interface are extruded into the adjacent subdomain to create new volume elements forming a one‐cell overlap. These new volume elements are used to compute a flux across the subdomain interface. An interface flux is computed independently for each subdomain. The values of the solution variables and other quantities for the nodes created by the extrusion process are determined by linear interpolation. The extrusion is done so that the interpolation will maintain information as localized as possible. The grid on the interface surface is arbitrary. The boundary between the two subdomains is completely independent from one another; meaning that they do not have to connect in a one‐to‐one manner and no symmetry or pattern restrictions are placed on the grid. A variety of numerical simulations were performed on model problems and large‐scale applications to examine conservation of the interface flux. Overall solution errors were found to be comparable to that for fully connected and fully conservative simulations. Excellent agreement is obtained with theoretical results and results from other solution methodologies. Copyright © 2006 John Wiley & Sons, Ltd.     

A new approach in modeling of constant temperature boundary condition in thermal lattice‐Boltzmann method

The lattice‐Boltzmann method is being applied to a diversity of fluid flow and heat transfer problems nowadays. Because of its microscale nature, strict attention should be paid when introducing macroscopic inputs to the model. One of the challenging issues dealing with macroscale and microscale treatment is the implementation of boundary conditions. In this regard constant‐temperature boundaries are frequently used in energy transfer problems. Such boundaries are simply modeled in Navier–Stokes based solvers, but they are not so harnessed in lattice‐Boltzmann models. One of the problems is that the calculated tangential heat flux is not zero along such boundaries in most of the previous models. In the present paper, a model has been developed, which has the capability of controlling tangential heat flux along the constant‐temperature boundaries. It aims to set the heat flux nearly zero along the boundary in midplane grid schemes. Copyright © 2012 John Wiley & Sons, Ltd.     

Choice of optimal blunting radius of axisymmetric body with account for surface radiation

We consider the problem of determining the optimal blunting with respect to heat transfer of an axisymmetric body in supersonic gas flow at zero angle of attack with account for body surface radiation. Results are presented of a calculation of the optimal blunting radius of a cone with half-angle 10° for various values of the ratio of convective heat flux at the stagnation point to the radiative heat flux. It is shown that for small values of this ratio the small bluntings are optimal.     

Application of Effective Flux Boundary Conditions to Two-Phase Upscaling in Porous Media

A new algorithm is introduced for upscaling relative permeabilities, and tested in simulations of two-dimensional reservoir displacement processes. The algorithm is similar to existing algorithms for computing upscaled relative permeabilities from subgrid simulations, but uses new boundary conditions for the pressure field. The new 'effective flux boundary conditions' were introduced in a previous paper and provide a more accurate estimate of flux through high permeability channels. The algorithm was tested in conjunction with uniform grid coarsening and upscaled absolute permeabilities for a broad range of coarsenings. The permeability fields were highly heteroge-neous and layered, and were obtained from synthetic data and from conditioned realizations of actual oil reservoirs. The algorithm was tested for a wide variety of grid aspect ratios, and for both viscous-and gravity-dominated flow. Typical fine grids were of the order of 100×100 cells; the coarsest scaled-up grids were on the order of 5×5 cells. The quality of scale up was evaluated by comparing oil cut curves for the fine and coarse grid simulations. We consistently obtained excellent agreement, even at the coarsest levels of scale up.     

Free-molecular gas flow in plane channels and grids

The problem of the flow of a rarefied gas in plane channels and grids in the free-molecular regime is considered. The parameters of the grid (channel) cells and the over-all geometrical dimensions of the grid are considered smaller than the mean free path of the molecules. The grid profiles are considered thin, the profile material is a perfect heat conductor. For the case of diffuse reflection of the molecules from the surface of the grid (channel), a calculation is made of the probability of passage of a molecule and the aerodynamic characteristics of the grid are determined.