计算流体力学中的验证与确认

计算流体力学(CFD)在航空航天等诸多领域的应用越来越广泛.特别是近年来, CFD在实际飞行器的设计中扮演着越来越重要的角色, 许多设计参数直接来源于CFD的计算结果.由此, 飞行器设计师对CFD提供结果的可信度提出了更高的要求.验证(verification)与确认(validation)是评价数值解精度和可信度的主要手段.本文综述了国内外开展CFD验证与确认研究的进展.在引言中论述了开展CFD验证与确认的重要性和必要性, 简述了国内外CFD验证与确认研究的历史和发展现状.第2节中讨论了CFD验证与确认的一些基本概念,以及这些概念定义的形成过程, 并指出了进行CFD验证与确认的基本步骤.第3节和第4节分别讨论了CFD验证与确认的方法, 如CFD验证中的精确解比较方法, 制造解比较方法, 网格收敛性研究;CFD确认中的层次结构, 流动分类法, 确认实验指南.在第5节中我们列举了几个CFD验证与确认的应用实例.最后, 对我国开展CFD验证与确认研究工作提出了若干建议, 包括:(1)开展流动分类法研究,(2)推行软件质量工程方法,(3)开展规范精细的实验, 建立国内的网络数据库.      

HyperFLOW软件非结构网格亚跨声速湍流模拟的验证与确认

计算流体力学(computational fluid dynamics, CFD)数值模拟在航空航天等领域发挥越来越重要的作用,然而CFD数值模拟结果的可信度仍然需要通过不断地验证与确认来提高.本文给出了从制造解精度测试、简单到复杂外形湍流模拟网格收敛性研究等三个方面开展CFD软件验证与确认的方法,并对自主研发的CFD软件平台HyperFLOW在非结构网格上模拟亚跨声速湍流问题的能力进行了验证与确认.首先通过基于Euler方程和标量扩散方程的制造解精度测试,分别验证了HyperFLOW在非结构网格上对Euler方程和黏性项的求解精度,结果表明其能够在任意非结构网格上达到设计的二阶精度.其次,通过NASA Turbulence Modeling Resource中的湍流平板、二维翼型近尾迹流动、二维Bump等几个典型的亚声速湍流算例的网格收敛性研究,量化考察了数值结果的观测精度阶和网格收敛性指数,并与国外知名CFD解算器CFL3D,FUN3D的计算结果进行了对比,验证了HyperFLOW对简单湍流问题的模拟能力,且具有良好的网格收敛性和计算精度(阶).最后,通过NASA Common Research Model标模定升力系数的网格收敛性研究和升阻极曲线预测,验证了软件在复杂外形亚跨声速湍流流动数值模拟中也具有良好的可信度.     

HyperFLOW软件非结构网格亚跨声速湍流模拟的验证与确认

计算流体力学(computational fluid dynamics,CFD)数值模拟在航空航天等领域发挥越来越重要的作用,然而CFD数值模拟结果的可信度仍然需要通过不断地验证与确认来提高.本文给出了从制造解精度测试、简单到复杂外形湍流模拟网格收敛性研究等三个方面开展CFD软件验证与确认的方法,并对自主研发的CFD软件平台HyperFLOW在非结构网格上模拟亚跨声速湍流问题的能力进行了验证与确认.首先通过基于Euler方程和标量扩散方程的制造解精度测试,分别验证了HyperFLOW在非结构网格上对Euler方程和黏性项的求解精度,结果表明其能够在任意非结构网格上达到设计的二阶精度. 其次,通过NASATurbulence Modeling Resource中的湍流平板、二维翼型近尾迹流动、二维Bump等几个典型的亚声速湍流算例的网格收敛性研究,量化考察了数值结果的观测精度阶和网格收敛性指数,并与国外知名CFD解算器CFL3D,FUN3D的计算结果进行了对比,验证了HyperFLOW对简单湍流问题的模拟能力,且具有良好的网格收敛性和计算精度(阶). 最后,通过NASA CommonResearchModel标模定升力系数的网格收敛性研究和升阻极曲线预测,验证了软件在复杂外形亚跨声速湍流流动数值模拟中也具有良好的可信度.      

基于POD和代理模型的静气动弹性分析方法

静气动弹性问题考虑弹性结构与定常气动力间的相互耦合作用,对飞行器的性能和安全具有显著的影响.在现代飞行器设计阶段,计算流体力学(CFD)/计算结构力学(CSD)直接耦合方法是精确考察静气动弹性影响的重要手段.然而,基于CFD技术的气动力仿真手段在耦合过程中计算量大且耗时长,难以满足设计阶段的需求.因此,为了兼顾计算精度与效率,文章采用本征正交分解(POD)和Kriging代理模型相结合的模型降阶方法,替代CFD求解过程并耦合有限元分析(FEA)方法,建立了高效、准确的静气动弹性分析框架.相较于传统的以模态法为主的静气动弹性分析方法,该方法能够解决更为复杂的静气动弹性问题以及提供静气动弹性变形过程中的气动分布载荷.针对典型三维跨声速HIRENASD机翼模型开展的马赫数、迎角变化的算例验证表明:由建立的静气动弹性分析方法与CFD/CSD直接耦合方法计算得到机翼翼梢处的静变形量间的相对误差在5%以内;同时该方法预测静平衡位置处的气动分布载荷的误差在5%以内,静气动弹性分析的计算效率至少提升了6倍.     

CFD结合降阶模型预测阵风响应

传统的阵风响应主要在频域内进行分析,气动载荷基于线性方法计算,不能考虑黏性和跨声速流动影响. 飞机设计需考虑不同频率和不同形状阵风的响应,基于CFD的阵风响应预测由于计算工况太多,工作量巨大. 本文发展了一种CFD结合非定常气动力ARMA(autoregressive-moving-averagemodel)降阶模型的阵风响应分析方法,CFD只要针对给定频率和形状的一种阵风响应进行计算,对获得的气动力时间历程运用线性最小二乘法参数辨识ARMA降阶模型的系数,则对任意频率和形状的阵风,代入降阶模型即可确定该阵风的响应,大大提高了计算效率. 为效验发展的方法,先计算NACA0012翼型在低马赫数0.11的阵风响应,通过对比CFD、ARMA降阶模型及早期发展的不可压阵风响应预测方法的结果,验证了方法的有效性. 再对比CFD、ARMA在跨声速马赫数0.8的阵风响应预测结果,证实所发展的方法对跨声速阵风响应预测亦是有效的.      

孤立两叶螺旋桨风洞试验准定常数值模拟

为研究某型两叶螺旋桨的气动性能,对某型两叶螺旋桨进行了风洞试验及计算流体动力学(CFD)数值模拟。通过求解多重参考坐标系(MRF)模型下的准定常雷诺平均Navier-Stokes方程,计算了直径为960mm的某型两叶木质螺旋桨静态和动态的气动特性。并在西北工业大学NF-3风洞的三元试验段获得了该螺旋桨桨叶的拉力、扭矩、功率等气动性能数据。螺旋桨风洞试验与CFD数值模拟对比验证表明:基于MRF模型的数值模拟结果与风洞试验结果具有一致性,两者的螺旋桨拉力的偏差控制在5%以内,扭矩的偏差控制在10%以内,拉力功率比的偏差控制在6.08%以内,证明了CFD数值模拟对螺旋桨气动性能预测的准确性。本文可为通用电动飞机螺旋桨的设计与模拟验证提供参考。     

覆冰导线气动力特性的规律性研究

采用数值模拟结合风洞试验验证的方法对影响覆冰导线气动力特性的主要因素进行规律性研究,建立了覆冰导线周围的空气流场二维CFD模型.对比湖北省中山口覆冰分裂导线的数值计算结果与风洞试验数据,验证CFD数值模拟方法的准确度和可靠性.在此基础上将CFD数值模拟方法推广至更多工况下的覆冰导线气动力系数数值模拟中,考察了覆冰导线气动力系数随覆冰形状、覆冰厚度、来流风速、导线分裂数目的影响规律.结果表明,该方法是一种实用、有效、准确的气动力系数数值模拟方法,所得气动力系数的变化规律对覆冰导线空气动力特性及舞动研究具有一定的指导意义.     

不同攻角下2:1矩形截面棱柱绕流的大涡模拟研究

本文研究了高雷诺数下(Re=24000),宽高比为2:1的矩形截面棱柱在不同风攻角下的绕流问题.基于大涡模拟理论和Smagorinsky亚格子应力模型,运用开源CFD软件OpenFOAM对Navier-Stokes控制方程进行求解,对均匀来流作用下的矩形截面棱柱绕流问题进行了三维数值模拟.通过对平均分量、脉动分量的比较,验证了数值模拟结果的准确性与可靠性,并给出了阻力系数、升力系数、斯特罗哈数随攻角的变化规律以及平均流场和湍流流场的流场特征.基于本文研究结果,通过开源CFD可视化软件ParaView给出了不同攻角下的流型变化,结合平均压力系数的分布,详细解释了不同流型的特征,这些结果验证了基于风洞实验所提出的假设流型,相关结论可为结构风工程设计提供参考依据.本文所采用的软件均为开源软件,其开源特性具有良好的科研价值.     

一种适用于超音速边界层的湍流转捩模式

建立一种合理反映扰动模态和可压缩性影响的新型k-ω-γ转捩模式.其主要特点为:(1)假设脉动动能k 由湍流脉动部分和非湍流脉动部分组成,且在模化后者时采用了稳定性分析的结果; (2)在间歇因子γ方程的源项中,构造了具有``自动判断转捩起始位置'功能的函数; (3)通过构造新型的物面法向长度尺度,保证了模式中所有的表达式均由当地变量构成,可以方便地应用于现代CFD程序之中. 该模式在亚音速、超音速和高超音速条件下的边界层流动中进行了验证. 计算结果表明,该模式可应用于较宽马赫数范围内的自然转捩以及旁路转捩过程,所具有的捕捉流动转捩的性能优于国际上的现有模式.      

基于制造解的非结构二阶有限体积离散格式的精度测试与验证

随着计算机技术的飞速进步,计算流体力学得到迅猛发展,数值计算虽能够快速得到离散结果,但是数值结果的正确性与精度则需要通过严谨的方法来进行验证和确认.制造解方法和网格收敛性研究作为验证与确认的重要手段已经广泛应用于计算流体力学代码验证、精度分析、边界条件验证等方面.本文在实现标量制造解和分量制造解方法的基础上,通过将制造解方法精度测试结果与经典精确解(二维无黏等熵涡)精度测试结果进行对比,进一步证实了制造解精度测试方法的有效性,并将两种制造解方法应用于非结构网格二阶精度有限体积离散格式的精度测试与验证,对各种常用的梯度重构方法、对流通量格式、扩散通量格式进行了网格收敛性精度测试.结果显示,基于Green-Gauss公式的梯度重构方法在不规则网格上会出现精度降阶的情况,导致流动模拟精度严重下降,而基于最小二乘(least squares)的梯度重构方法对网格是否规则并不敏感.对流通量格式的精度测试显示,所测试的各种对流通量格式均能达到二阶精度,且各方法精度几乎相同;而扩散通量离散中界面梯度求解方法的选择对流动模拟精度有显著影响.     

CFD in research for the petrochemical industry

Existing computer programs for turbulent flow simulations have been tested on their capacity to solve practical problems in engineering research for the petrochemical industry. This paper presents results of six validation studies which are all related to experimental work carried out in the Shell Group laboratories. On the basis of these studies needs for further improvement and validation of both numerical schemes and turbulence models have been formulated.The present as well as future relevance of Computational Fluid Dynamics (CFD) in engineering research is discussed in a wider context than covered by the validation studies. This analysis formed the basis for the identification of a number of CFD areas which have great potential in engineering application but need substantial further research and development to make them suitable for this purpose.     

Flow measurement around a model ship with propeller and rudder

For the design of hull forms with better resistance and propulsive performance, it is essential to understand flow characteristics, such as wave and wake development, around a ship. Experimental data detailing the local flow characteristics are invaluable for the validation of the physical and numerical modeling of computational fluid dynamics (CFD) codes, which are recently gaining attention as efficient tools for hull form evaluation. This paper describes velocity and wave profiles measured in the towing tank for the KRISO 138,000 m³ LNG carrier model with propeller and rudder. The effects of propeller and rudder on the wake and wave profiles in the stern region are clearly identified. The results contained in this paper can provide an opportunity to explore integrated flow phenomena around a model ship in the self-propelled condition, and can be added to the International Towing Tank Conference benchmark data for CFD validation as the previous KCS and KVLCC cases.     

CFD modeling of hydrodynamics and mass transfer of Rhodamine B in annular reactor

The hydrodynamics and mass transfer are the two crucial issues in annular reactors. An accurate prediction of these issues is required for design, optimization and scale-up applications. The present study deals with the modeling and simulation of flow through an annular reactor at different hydrodynamic conditions using computational fluid dynamics (CFD) to investigate the hydrodynamics and mass transfer. CFD modeling was utilized to predict velocity distribution, average velocity and average mass transfer coefficient in the annular geometry. The results of CFD simulations were compared with the mathematically derived equations and already developed correlations for validation purposes. CFD modeling was found suitable for predicting hydrodynamics and mass transfer for annular geometry under laminar flow conditions. It was observed that CFD also provides local values of the parameters of interest in addition to the average values for the simulated geometry.     

Verification and validation of impinging round jet simulations using an adaptive FEM

This paper illustrates the use of an adaptive finite element method as a means of achieving verification of codes and simulations of impinging round jets, that is obtaining numerical predictions with controlled accuracy. Validation of these grid‐independent solution is then performed by comparing predictions to measurements. We adopt the standard and accepted definitions of verification and validation (Technical Report AIAA‐G‐077‐1998, American Institute of Aeronautics and Astronautics, 1998; Verification and Validation in Computational Science and Engineering. Hermosa Publishers: Albuquerque, NM, 1998). Mesh adaptation is used to perform the systematic and rigorous grid refinement studies required for both verification and validation in CFD. This ensures that discrepancies observed between predictions and measurements are due to deficiencies in the mathematical model of the flow. Issues in verification and validation are discussed. The paper presents an example of code verification by the method of manufactured solution. Examples of successful and unsuccessful validation for laminar and turbulent impinging jets show that agreement with experiments is achieved only with a good mathematical model of the flow physics combined with accurate numerical solution of the differential equations. The paper emphasizes good CFD practice to systematically achieve verification so that validation studies are always performed on solid grounds. Copyright © 2004 John Wiley & Sons, Ltd.     

CFD analysis of the erosion of a light gas stratification by means of a hot air jet in the MiniPanda facility

Following the Fukushima-Daiichi accident, many countries decided to strengthen the safety systems of their nuclear power plants in order to increase the capabilities of managing severe nuclear accidents. To achieve this goal detailed analysis of postulated design and beyond-design-basis accidents is essential. The main tools available for the analysis of these complex scenarios are Advanced Lumped Parameters (LP) and Computational Fluid Dynamics (CFD) codes. Currently, the most limiting factor in the application of these tools is their validation. This activity is performed in support of the development and validation of a CFD model, built in the numerical environment of OpenFOAM, for hydrogen behavior in the containment of Light Water Reactors. For validation purpose, numerical simulations of the erosion of a light gas stratification by means of a hot air jet have been performed, and results of the simulations were analysed and compared against experimental results obtained by Ritterath (2012) in the MiniPanda facility. Two different scenarios have been considered, each characterized by a different value of the mass flow rate of the jet used to erode the stratification.     

Experimental measurements and computational modeling of the flow field in an idealized human oropharynx

The velocity field in the central sagittal plane of an idealized representation of the human oropharynx (HOP) during steady inspiration, simulating oral inhalation through an inhaler mouthpiece, was measured experimentally using endoscopic particle image velocimetry (PIV). Measurements were made at three flow rates: 15, 30, and 90 L/min, which correspond to a wide range of physiological conditions. Extensive tests were performed to verify the veracity of the PIV data. The flow was also modeled computationally using Reynolds-averaged Navier–Stokes (RANS) computational fluid dynamics (CFD) methods. The PIV data clearly indicate the complex nature of HOP flow, with three-dimensionality and several regions of separation and recirculation evident. Comparison of the experimental and computational results shows that, although the RANS CFD reproduces the basic features of the flow, it does not adequately capture the increased viscous effects at lower Reynolds numbers. The results demonstrate the need for more development and validation of CFD modeling, in particular RANS methods, in these flows.     

Steam condensing flow modeling in turbine channels

The numerical method for modeling of the transonic steam flows with homogeneous and/or heterogeneous condensation has been presented. The experiments carried out for the Laval nozzles, for 2-D turbine cascades and for a 3-D flow in real turbine were selected to validate an in-house CFD code adjusted to the calculations of the steam condensing flows in complicated geometries. The sensitivity of the condensation model and difficulties in the validation process of the CFD code have been discussed. These difficulties limit the possibilities of verification and improvement of the condensation theory based on the existing experimental data.     

2D numerical simulation of passive autocatalytic recombiner for hydrogen mitigation

Resolving hydrogen related safety issues, pertaining to nuclear reactor safety has been an important area of research world over for the past decade. The studies on hydrogen transport behavior and development of hydrogen mitigation systems are still being pursued actively in various research labs, including Bhabha Atomic Research Centre (BARC), in India. The passive autocatalytic recombiner (PAR) is one of such hydrogen mitigating device consisting of catalyst surfaces arranged in an open-ended enclosure. In the plate type recombiner design sheets made of stainless steel and coated with platinum catalyst material are arranged in parallel inside a flow channel. The catalyst elements are exposed to a constant flow of a mixture of air, hydrogen and steam, a catalytic reaction occurs spontaneously at the catalyst surfaces and the heat of reaction produces natural convection flow through the enclosure. Numerical simulation and experiments are required for an in-depth knowledge of such plate type PAR. Specific finite volume based in-house 2D computational fluid dynamics (CFD) code has been developed to model and analyse the working of these recombiners and has been used to simulate one literature quoted experiment. The validation results were in good agreement against literature quoted German REKO experiments. Parametric study has been performed for particular recombiner geometry for various inlet conditions. Salient features of the simplified CFD model developed at BARC and results of the present model calculations are presented in this paper.     

Development of co-current air–water flow in a vertical pipe

Measurements of the cross-sectional distribution of the gas fraction and bubble size distributions were conducted in a vertical pipe with an inner diameter of 51.2 mm and a length of about 3 m for air/water bubbly and slug flow regimes. The use of a wire-mesh sensor obtained a high resolution of the gas fraction data in space as well as in time. From this data, time averaged values for the two-dimensional gas fraction profiles were decomposed into a large number of bubble size classes. This allowed the extraction of the radial gas fraction profiles for a given range of bubble sizes as well as data for local bubble size distributions. The structure of the flow can be characterized by such data. The measurements were performed for up to 10 different inlet lengths and for about 100 combinations of gas and liquid volume flow rates. The data is very useful for the development and validation of meso-scale models to account for the forces acting on a bubble in a shear liquid flow and models for bubble coalescence and break-up. Such models are necessary for the validation of CFD codes for the simulation of bubbly flows.