迎风通量差分分裂法对边界层特性的数值模拟

迎风通量差分分裂法作为一种高效率的数值方法,已被广泛地用于模拟各种复杂流场,但用这类方法模拟有边界层存在的超音速粘性流场时,必须保证通量差分分裂为守恒的;否则,将在边界层内引入一个可观的虚假人工扩散项和压力梯度项,歪曲了边界层特性。本文就这个问题进行理论和数值的研究。     

非结构混合网格化学非平衡粘性绕流数值模拟

在非结构混合网格上对化学非平衡粘性绕流进行了数值模拟。控制方程为考虑了化学非平衡效应的二维Navier-Stokes方程,化学动力学模型为7组元、7反应模型。控制方程中的对流项采用VanLeer逆风分裂格式处理,并应用MUSCL方法及Minmod限制器扩展到二阶精度,粘性项用中心差分格式处理。时间推进采用显式5步龙格-库塔方法。为了适应高超声速流场计算,对VanLeer通量分裂方法进行了改进,并引入了化学反应时间步长。对RAMC-II模型的飞行试验流场进行了数值模拟,计算结果与试验测量数据符合较好,并与参考文献中的数值模拟结果吻合。     

用隐式方法和WENO格式计算悬停旋翼跨声速无粘流场

寻找一种能够准确计算以涡为主要特征的复杂流场和克服尾迹耗散问题的数值方法,一直是旋翼空气动力学研究的热点和难点。本文发展了一种基于高阶迎风格式计算悬停旋翼无粘流场的隐式数值方法。无粘通量采用Roe通量差分分裂格式,为提高精度,使用五阶WENO格式进行左右状态插值,并与MUSCL插值进行比较。为提高收敛到定常解的效率,时间推进采用LU-SGS隐式方法。用该方法对一跨声速悬停旋翼无粘流场进行了数值计算,数值结果表明WENO-Roe的激波分辨率高于MUSCL-Roe,体现出了格式精度的提高对计算结果的改善,LU-SGS隐式方法的计算效率比5步Runge-Kutta显式方法的高。     

复杂无粘流场数值模拟的矩形/三角形混合网格技术

建立了一套模拟复杂无粘流场的矩形／三角形混合网格技术，其中三角形仅限于物面附近，发挥非结构网格的几何灵活性，以少量的网格模拟复杂外型；同时在以外的区域采用矩形结构网格，发挥矩形网格计算简单快速的优势，有效地克服全非结构网格计算方法需要较大内存量和较长ＣＰＵ时间的不足．混合网格系统由修正的四分树方法生成．将ＮＮＤ有限差分与ＮＮＤ有限体积格式有机地融合于混合网格计算，消除了全矩形网格模拟曲边界的台阶效应，同时保证了网格间的通量守恒．数值实验表明本方法在模拟复杂无粘流场方面的灵活性和高效性．     

三维高超音速球锥粘性绕流的计算

近二十年内,三维边界层方程的数值计算方法有了很大发展。由于求解三维边界层方程必须先得到无粘流场,或者二者进行迭代,其所需计算机的存储量及计算时间无疑是大量的。为此,发展和完善三维简化N-S方程(亦称广义边界层方程)的数值解法就日益显得重要起来。Rubin等及Helliwell等对抛物化的N-S方程(保留侧向扩散项)提     

应用流体双折射方法计算二维二尖瓣模型的流场

本文运用流体双折射的实验数据与数值差分计算的杂交法对模拟二尖瓣模型的二维稳态流场进行定量计算,得到流场的流函数、速度等流场参量。     

激波风洞内超燃冲压发动机三面压缩进气道流场实验观测

主要进行了超燃冲压发动机三面压缩进气道的实验观测。利用来流马赫数4.5的直通式激波风洞,考察了三组具有不同压缩角度的进气道模型内部的流场情况。实验观测手段为油流法、丝线法和高速纹影,同时,辅以数值模拟以有助于流场细节分析。纹影照片展示了进气道内部以激波边界层相互作用为主要影响因素的流场复杂结构,数值模拟也显示了相近的结果。油流技术与丝线法显示了近壁面处的流动图像,照片中可见激波、分离线、再附线等分界线位置。根据实验结果,可以推测唇口激波与进气道内边界层的相互作用及其引起的壁面分离是影响进气道内流动的主要因素。同时,尝试了利用抽吸方法减弱激波与边界层相互作用诱发的壁面流动分离,并取得一定结果。     

注塑成型充模过程的温度场计算

对注塑成型过程中充填阶段的非等温效应的数值模拟进行了研究 ,模拟了非牛顿流体在任意形状型腔内的流动。这里 ,数学模型是基于非牛顿流体在非等温状态下的广义 Hele-Shaw流动 ,利用有限元 /有限差分混合数值方法求解流场中的压力和温度分布     

高超声速激波湍流边界层干扰直接数值模拟研究

高超声速激波与湍流边界层干扰会导致飞行器表面出现局部热流峰值,严重影响飞行器气动性能和飞行安全. 针对高马赫数激波干扰问题,以往数值研究多采用雷诺平均方法,而在直接数值模拟方面的相关工作较为少见. 开展高超声速激波与湍流边界层干扰的直接数值模拟研究,有助于进一步提升对其复杂流动机理认识和理解,同时也将为现有湍流模型和亚格子应力模型的改进提供理论依据. 采用直接数值模拟方法对来流马赫数6.0,34°压缩拐角内激波与湍流边界层的干扰问题进行了研究. 基于雷诺应力各向异性张量,分析了高超声速湍流边界层在压缩拐角内的演化特性. 通过对湍动能输运方程的逐项分析,系统地研究了可压缩效应对湍动能及其输运的影响机制. 采用动态模态分解方法,探讨了干扰流场的非定常运动历程. 研究结果表明,随着湍流边界层往下游发展,近壁湍流的雷诺应力状态由两组元轴对称状态逐渐演化为两组元状态,外层区域则由轴对称膨胀趋近于各向同性. 干扰流场内存在强内在压缩性效应(声效应),其对湍动能输运的影响主要体现在压力--膨胀项,而对膨胀--耗散项影响较小. 高超声速下压缩拐角内的非定常运动仍存在以分离泡膨胀/收缩为特征的低频振荡特性,其物理机制与分离泡剪切层密切相关.      

一种激波稳定的对流-压力通量分裂格式

针对欧拉方程三种流行的对流-压力通量分裂方法(Liou-Steffen,Zha-Bilgen和Toro-Vázquez)进行特征分析,进而提出一种新的对流-压力通量分裂格式。采用Zha-Bilgen分裂方法将欧拉方程的通量分裂成对流项和压力项两部分,使用TV格式来计算这两部分的数值通量。利用压力比构造激波探测函数,并且在强激波附近的亚声速区域增加TV格式的剪切粘性来克服数值模拟中的激波不稳定性。数值算例的计算结果表明,新的对流-压力通量分裂格式不仅保留了原始TV格式精确分辨接触间断的优点,而且具有更好的鲁棒性,在数值模拟多维强激波问题时不会出现不稳定现象。因此,该格式是一种精确并且具有强鲁棒性的数值方法,可以广泛地应用于可压缩流体的数值计算中。     

An extended study of the hydrodynamics of gravity-driven film flow of power-law fluids

A new similarity transformation has been devised for extensive studies of accelerating non-Newtonian film flow. The partial differential equations governing the hydrodynamics of the flow of a power-law fluid down along an inclined plane surface are transformed into a set of two ordinary differential equations by means of the dimensionless velocity component approach. Although the analysis is applicable for any angle of inclination (0<π/2), the resulting one-parameter problem involves only the power-law index n. Nevertheless, physically essential quantities, like the velocity components and the skin-friction coefficient, do depend on and relevant relationships are deduced between the vertical and inclined cases. Accurate numerical similarity solutions are provided for n in the range from 0.1 to 2.0. The present method enables solutions to be obtained also for highly pseudo-plastic films, i.e. for n below 0.5. The mass flow rate entrained into the momentum boundary layer from the inviscid freestream is expressed in terms of a dimensionless mass flux parameter Φ, which depends on the dimensionless boundary layer thickness and the velocity components at the edge of the viscous boundary layer. Φ, which is thus an integral part of the similarity solution, turns out to decrease monotonically with n. This parameter is of particular relevance in the determination of the streamwise position at which the entire freestream has been entrained and viscous stresses prevail all the way to the free surface of the film. A short-cut method to facilitate rapid and yet accurate estimates of the mass flux parameter is developed to this end.     

Transient flow rate behaviour in an external natural convection boundary layer

A theoretical approach is proposed to investigate the transient dynamic behaviour of a free convection boundary layer-type flow. The set of continuity, momentum and energy equations are solved with the classical Boussinesq approximation using the Karman–Pohlhausen integral method. Applying a step variation of the uniform heat flux on a vertical wall, the boundary layer thickness and velocity profiles within the viscous layer, streamline patterns and volumetric flow rate are evaluated as a function of time. In addition, corresponding fully analytical asymptotic solutions are derived to be readily used in engineering applications.     

Calculation of Reynolds stresses in turbulent supersonic flows

The baseline numerical procedure of interest in this study combines flux vector splitting, flux difference splitting and an explicit treatment of the diffusion terms of the flow equations. The viscous terms are treated explicitly to preserve the wave propagation properties of the Euler fluxes and permit splitting. The experience with this scheme has been limited to laminar or, at best, ‘eddy viscosity’ flows. In this paper the applicability of the scheme is extended to include the calculation of turbulent Reynolds stresses in supersonic flows. The schemes and our implementation are discussed. Both laminar and turbulence subsets of the Reynolds/Favre-averaged equations are tested, with a discussion of relative performance. The test problem for turbulence consists of a zero-pressure-gradient supersonic boundary layer as well as a supersonic boundary layer experiencing the combined effects of adverse pressure gradient, bulk compression and a concave streamline curvature. Excellent agreement with experimental measurements is observed for most of the quantities compared, which suggests that the numerical procedures presented in this paper are potentially very useful.     

Numerical investigation of condensing steam flow in boundary layers

The paper describes a numerical method for the prediction of condensing steam flow within compressible boundary layers. The method is based on a simple stream function technique, which enables straightforward integration of the nucleation and droplet growth equations in a Lagrangian frame of reference. Calculations show how viscous dissipation and reduced expansion rate within a typical boundary layer influence nucleation and growth, leading to droplet radii and size distributions that differ substantially from those predicted in inviscid flow. The impact of condensation on temperature and velocity profiles, and the implications for thermodynamic loss are also considered.     

Numerical solution of the problem of the mixing of the boundary layers shed from the trailing edge of a wing

During the mixing of viscous incompressible flows with different velocities, in the vicinity of a trailing edge an interaction region with a three-layer structure is formed, similar to that in the case of symmetric shedding with equal velocities. The boundary layers developing on the upper and lower sides of the airfoil form a viscous mixing layer, or vortex sheet, which separates the flows downstream of the trailing edge. The boundary value problem corresponding to the flow in the viscous sublayer in the vicinity of the trailing edge of a flat plate is solved for high Reynolds numbers using an efficient numerical method for solving the equations of asymptotic interaction theory.     

Comment on the paper “Transient MHD free convective flow past an infinite vertical plate embedded in a porous medium with viscous dissipation,Siva Reddy Sheri,R. Srinivasa Raju,Meccanica, DOI 10.1007/s11012-015-0285-y”

A numerical investigation of transient magnetohydrodynamic free convection flow past an infinite vertical plate embedded in a porous medium with viscous dissipation is presented in the above paper. The governing differential equations are transformed into a set of non-linear coupled partial differential equations and are solved numerically using the finite element method. Numerical results for the velocity, temperature and concentration profiles within the boundary layer are presented and discussed.     

Computational modeling of heat transfer over an unsteady stretching surface embedded in a porous medium

The present paper deals with the study of heat transfer characteristics in the laminar boundary layer flow of an incompressible viscous fluid over an unsteady stretching sheet which is placed in a porous medium in the presence of viscous dissipation and internal absorption or generation. Similarity transformations are used to convert the governing time dependent nonlinear boundary layer equations into a system of non-linear ordinary differential equations containing Prandtl number, Eckert number, heat source/sink parameter, porous parameter and unsteadiness parameter with appropriate boundary conditions. These equations are solved numerically by applying shooting method using Runge-Kutta-Fehlberg method. Comparison of numerical results is made with the earlier published results under limiting cases. The effects of the parameters which determine the velocity and temperature fields are discussed in detail.     

ANALYSIS/DESIGN CALCULATIONS OF TRANSONIC FLOWS

SUMMARY

Analysis/design calculations of transonic flow are discussed and several improvements are made. The nonisentropic potential method is used to calculate the inviscid transonic flow analysis problem instead of the traditional potential method. An inverse integral 3D boundary layer method is used to calculate the boundary layer in the viscous transonic flow analysis problem. The viscous/inviscid interaction calculations are carried out by a semi-inverse coupling scheme. In design problem calculations, an improved residual-correction method is used. Three individual methods are combined in a global algorithm and computing code. The improvements speed up the convergence, increase applicability and computational efficiency. Some numerical results are given to illustrate that the present method provides an effective engineering tool of high accuracy and efficiency in three dimensional transonic analysis and design situations.     

Flow and heat transfer over hyperbolic stretching sheets

The boundary layer flow and heat transfer analysis of an incompressible viscous fluid for a hyperbolically stretching sheet is presented. The analytical and numerical results are obtained by a series expansion method and a local non-similarity (LNS) method, respectively. The analytical and numerical results for the skin friction and the Nusselt number are calculated and compared with each other. The significant observation is that the momentum and the thermal boundary layer thickness decrease as the distance from the leading edge increases. The well-known solution of linear stretching is found as the leading order solution for the hyperbolic stretching.