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

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

湍流冲击射流流动与传热的数值研究进展

湍流冲击射流由于其冲击表面时具有很高的局部传热率和冲击力,被广泛应用于如表面的加热、电子元件的冷却、纸张的干燥和材料的切割等工程应用和工业过程中.由于其流动的复杂性,也常被作为一种理想的测试实例来评价湍流模型的性能.此外,湍升力射流与地面之间的空气动力作用对V/STOL (垂直或短距离起落)飞机的性能具有很大的影响.长期以来,人们从理论分析、实验测量和数值模拟方面对冲击射流进行了广泛而系统的研究,积累了丰富的资料.本文在分析了湍流冲击射流的数值研究现状的基础上,对近年来有关湍流冲击射流流动与传热的数值研究方面的文献有选择地进行了综述,重点评述了不同湍流模型对冲击射流流动与传热的预测能力,讨论了存在的问题并对该领域今后的研究方向进行了展望.      

A comparison of flow visualisation and wall pressure measurements for a jet impinging on a plane surface

Several studies of jets impinging on a plane surface have already been made. This paper suggests a new approach to studying impingements of jets. Comparisons have been made between visualisation results and wall pressure measurements. It is shown that only one of these two techniques is sufficient for characterising the flow nature near the wall. Visualisations can be sufficient for determining the location of wall pressure maxima and minima. Detachments and reattachments of the flow are thus located, and the main characteristics of a jet impinging on a plane wall can be shown by simpler experiments such as the spreading over method. Received: 15 October 1998/Accepted: 27 July 1999     

Variational boundary integral approach for asymmetric impinging jets of arbitrary two‐dimensional nozzle

We consider asymmetric impinging jets issuing from an arbitrary nozzle. The flow is assumed to be two‐dimensional, inviscid, incompressible, and irrotational. The impinging jet from an arbitrary nozzle has a couple of separated infinite free boundaries, which makes the problem hard to solve. We formulate this problem using the stream function represented with a specific single layer potential. This potential can be extended to the surrounding region of the jet flow, and this extension can be proved to be a bounded function. Using this fact, the formulation yields the boundary integral equations on the entire nozzle and free boundary. In addition, a boundary perturbation produces an extraordinary boundary integral equation for the boundary variation. Based on these variational boundary integral equations, we can provide an efficient algorithm that can treat with the asymmetric impinging jets having arbitrarily shaped nozzles. Particularly, the proposed algorithm uses the infinite computational domain instead of a truncated one. To show the convergence and accuracy of the numerical solution, we compare our solutions with the exact solutions of free jets. Numerical results on diverse impinging jets with nozzles of various shapes are also presented to demonstrate the applicability and reliability of the algorithm.     

On the impact of anisotropic mesh adaptation on computational wind engineering

This paper addresses the critical issue of the accuracy of CFD predictions for wind engineering. Flows around the Silsoe Cube, a high‐rise building (the Jin Mao Tower), and a low‐rise large‐span building (the Pudong International Airport) are computed with the Navier–Stokes solver FENSAP and compared with measurements. Computations are carried out for two wind directions, by solving the steady‐state ensemble‐averaged Navier–Stokes equations with the Spalart–Allmaras one‐equation turbulence model. Pressure coefficients compare well with wind tunnel experiments and the accuracy of the flow solutions is further improved via an automatic mesh adaptation that dynamically places grid points where the flow physics require them, while keeping the number of unknowns and solution time substantially at the same level. Copyright © 2009 John Wiley & Sons, Ltd.     

A simple sliding‐mesh interface procedure and its application to the CFD simulation of a tidal‐stream turbine

An effective way of using computational fluid dynamics (CFD) to simulate flow about a rotating device—for example, a wind or marine turbine—is to embed a rotating region of cells inside a larger, stationary domain, with a sliding interface between. This paper describes a simple but effective method for implementing this as an internal Dirichlet boundary condition, with interfacial values obtained by interpolation from halo nodes. The method is tested in two finite‐volume codes: one using block‐structured meshes and the other unstructured meshes. Validation is performed for flow around simple, isolated, rotating shapes (cylinder, sphere and cube), comparing, where possible, with experiment and the alternative CFD approach of fixed grid with moving walls. Flow variables are shown to vary smoothly across the sliding interface. Simulations of a tidal‐stream turbine, including both rotor and support, are then performed and compared with towing‐tank experiments. Comparison between CFD and experiment is made for thrust and power coefficients as a function of tip‐speed ratio (TSR) using Reynolds‐averaged Navier–Stokes turbulence models and large‐eddy simulation (LES). Performance of most models is good near the optimal TSR, but simulations underestimate mean thrust and power coefficients in off‐design conditions, with the standard k– ? turbulence model performing noticeably worse than shear stress transport k– ω and Reynolds‐stress‐transport closures. LES gave good predictions of mean load coefficients and vital information about wake structures but at substantial computational cost. Grid‐sensitivity studies suggest that Reynolds‐averaged Navier–Stokes models give acceptable predictions of mean power and thrust coefficients on a single device using a mesh of about 4 million cells. Copyright © 2013 John Wiley & Sons, Ltd.     

Suitability of the k–ω turbulence model for scramjet flowfield simulations

The suitability of Wilcox's 2006 k– ω turbulence model for scramjet flowfield simulations is demonstrated by validation against five test cases that have flowfields representative of those to be expected in scramjets. The five test cases include a 2D flat plate, an axisymmetric cylinder, a backward‐facing step, the mixing of a pair of coaxial jets and the interaction between a shock wave and turbulent boundary layer. A generally good agreement between the numerical and experimental results is obtained for all test cases. These tests reveal that despite the turbulence model's sensitivity to freestream turbulence properties, the numerically predicted skin friction agrees with experimental data and theoretical correlations to their degree of uncertainty. The tests also confirm the importance of using a y⁺ value of less than 1 in getting accurate surface heat transfer distributions. In the coaxial jets case, the importance of matching the turbulence intensities at the inflow plane in improving the predictions of the turbulent mixing phenomena is also shown. A review of guidelines with regard to the setting up of grids and specification of freestream turbulence properties for turbulent Reynolds‐averaged Navier–Stokes CFD simulations is also included in this paper. Copyright © 2011 John Wiley & Sons, Ltd.     

Aerodynamic interference between high-speed slender bodies

Under supersonic flow conditions, slender bodies in close proximity induce aerodynamic interference effects. This paper aims to quantify the magnitude of the resulting interference loads and to understand the underlying flow-physics mechanisms that cause them. A pair of identical slender bodies are investigated through a series of wind-tunnel experiments and supporting computational fluid dynamics (CFD) predictions. The bodies induce a complex interference flowfield, which tends to be bespoke to each configuration. The flow features include impinging shock and expansion waves, conical shock reflections, strong skewing of the boundary-layer flows and shock diffraction. The effects of axial stagger, lateral separation and the strength of the primary disturbance flow field are evaluated. The interference loads are found to be most sensitive to the initial location of the primary disturbance but are also affected by its strength. In addition, maximum interference loads which equate to an effective incidence of up to 6° are observed. Finally, very good agreement is found between the measurements and the CFD predicted normal force and pitching moment.     

Droplet size and velocity characteristics of water-air impinging jet atomizer

A water-air impinging jets atomizer is investigated in this study, which consists of flow visualization using high speed photography and mean droplet size and velocity distribution measurements of the spray using Phase Doppler Anemometry (PDA). Topological structures and break up details of the generated spray in the far and near fields are presented with and without air jet and for an impinging angle of 90°. Spray angle increases with the water jet velocity, air flow rate and impinging angle. PDA results indicate that droplet size is smallest in the spray center, with minimum value of Sauter mean diameter (SMD) of 50 µm at the air flow rate of Q_m = 13.50 g/min. SMD of droplets increases towards the spray outer region gradually to about 120 µm. The mean droplet velocity component W along the air-jet axis is highest in the spray center and decreases gradually with increasing distance from the spray center. SMD normalized by the air nozzle diameter is found firstly to decrease with gas-to-liquid mass ratio (GLR) and air-to-liquid momentum ratio (ALMR) and then remain almost constant. Its increasing with aerodynamic Weber number indicates an exponential variation. The study sheds light on the performance of water-air impinging jets atomizers providing useful information for future CFD simulation works.     

Research on multiphase flows in Thermo-Energy Engineering Institute of Southeast University

Multiphase flows have received increasing attention over the past decades. This paper describes the research carried out in Thermo-Energy Engineering Institute of Southeast University in recent years, focusing on several common issues associated with multiphase flows in industry, such as: boiling of falling film and complex structure of gas–liquid flow under large difference in temperature, free surface flows involving liquid jets and drop formation, mixing behaviors of gas–liquid–solid three-phase flow, and fluidization characteristics of cylindrical particles. Numerical methods ranging from empirical to CFD models were developed for predictions, and experimental works were essentially conducted for validation and modification. For all cases, simulated results were validated with experiments and good agreements were obtained. Based on the combined modeling and experimental approach, fundamental understanding of multiphase processes in a specific circumstance is achieved under conditions relevant to the actual industrial-scale, such as transport phenomena, flow patterns, fluid dynamics and interactions between phases.     

Validation of CFD‐CSD coupling interface methodology using commercial codes

Coupling interface between computational fluid dynamics (CFD) and computational structural dynamics (CSD) is required to provide exchange of information for the simulation of fluid–structure interaction (FSI) phenomena. Accuracy and consistency of information exchanged through coupling interface between the independent CFD and CSD solvers plays a central role in the simulation and prediction of FSI phenomenon, like flutter. In this paper validation of an implemented coupling interface methodology is presented for subsonic, transonic and near supersonic mach regime. The test case chosen for this purpose is the flutter of AGARD445.6 standard I‐wing weakened model configuration for subsonic to near transonic flow regime. Gambit^® and Fluent^® are used for CFD grid generation and solution of fluid dynamic equations, respectively. CSD modeling and simulation are provided by numerical time integration of modal dynamic equations derived through the finite element modeling in ANSYS^® environment. Copyright © 2009 John Wiley & Sons, Ltd.     

Experiments on free and impinging supersonic microjets

The fluid dynamics of microflows has recently commanded considerable attention because of their potential applications. Until now, with a few exceptions, most of the studies have been limited to low speed flows. This experimental study examines supersonic microjets of 100–1,000 μm in size with exit velocities in the range of 300–500 m/s. Such microjets are presently being used to actively control larger supersonic impinging jets, which occur in STOVL (short takeoff and vertical landing) aircraft, cavity flows, and flow separation. Flow properties of free as well as impinging supersonic microjets have been experimentally investigated over a range of geometric and flow parameters. The flowfield is visualized using a micro-schlieren system with a high magnification. These schlieren images clearly show the characteristic shock cell structure typically observed in larger supersonic jets. Quantitative measurements of the jet decay and spreading rates as well as shock cell spacing are obtained using micro-pitot probe surveys. In general, the mean flow features of free microjets are similar to larger supersonic jets operating at higher Reynolds numbers. However, some differences are also observed, most likely due to pronounced viscous effects associated with jets at these small scales. Limited studies of impinging microjets were also conducted. They reveal that, similar to the behavior of free microjets, the flow structure of impinging microjets strongly resembles that of larger supersonic impinging jets.     

Experimental study of an impinging jet with different swirl rates

A stereo PIV technique using advanced pre- and post-processing algorithms is implemented for the experimental study of the local structure of turbulent swirling impinging jets. The main emphasis of the present work is the analysis of the influence of swirl rate on the flow structure. During measurements, the Reynolds number was 8900, the nozzle-to-plate distance was equal to three nozzle diameters and the swirl rate was varied from 0 to 1.0. For the studied flows, spatial distributions of the mean velocity and statistical moments (including triple moments) of turbulent pulsations were measured.

The influence of the PIV finite spatial resolution on the measured dissipation rate and velocity moments was analyzed and compared with theoretical predictions. For this purpose, a special series of 2D PIV measurements was carried out with vector spacing up to several Kolmogorov lengthscales.

All terms of the axial mean momentum and the turbulent kinetic energy budget equations were obtained for the cross-section located one nozzle diameter from the impinging plate. For the TKE budget, the dissipation term was directly calculated from the instantaneous velocity fields, thereby allowing the pressure diffusion term to be found as a residual one. It was found that the magnitude of pressure diffusion decreased with the growth of the swirl rate. In general, the studied swirling impinging jets had a greater spread rate and a more rapid decay in absolute velocity when compared to the non-swirling jet.     

Determination of heat transfer correlations for plate-fin-and-tube heat exchangers

This paper presents a numerical method for determining heat transfer coefficients in cross-flow heat exchangers with extended heat exchange surfaces. Coefficients in the correlations defining heat transfer on the liquid- and air-side were determined using a non-linear regression method. Correlation coefficients were determined from the condition that the sum of squared fluid temperature differences at the heat exchanger outlet, obtained by measurements and those calculated, achieved minimum. Minimum of the sum of the squares was found using the Levenberg-Marquardt method. The outlet temperature of the fluid leaving the heat exchanger was calculated using the mathematical model describing the heat transfer in the heat exchanger. Since the conditions at the liquid-side and those at the air-side are identified simultaneously, the derived correlations are valid in a wide range of flow rate changes of the air and liquid. This is especially important for partial loads of the exchanger, when the heat transfer rate is lower than the nominal load. The correlation for the average heat transfer coefficient on the air-side based on the experimental data was compared with the correlation obtained from numerical simulation of 3D fluid and heat flow, performed by means of the commercially available CFD code. The numerical predictions are in good agreement with the experimental data.     

COMPUTATIONAL STUDIES OF IMPINGING JETS USING K-ϵ TURBULENCE MODELS

This paper reports numerical modelling of impinging jet flows using Rodi and Malin corrections to the k–ϵ turbulence model, carried out using the PHOENICS finite volume code. Axisymmetric calculations were performed on single round free jets and impinging jets and the effects of pressure ratio, height and nozzle exit velocity profile were investigated numerically. It was found that both the Rodi and Malin corrections tend to improve the prediction of the hydrodynamic field of free and impinging jets but still leave significant errors in the predicted wall jet growth. These numerical experiments suggest that conditions before impingement significantly affect radial wall jet development, primarily by changing the wall jet's initial thickness.     

Predictions of axisymmetric and two-dimensional impinging turbulent jets

The k − turbulence model and a version of a second-moment closure, modified to include the effect of pressure reflections from a solid surface, have been used as the basis of predictions of the flow that results from the orthogonal impingement of circular and two-dimensional (2-D) jets on a flat surface. Comparison of model predictions has been made with velocity measurements obtained in the stagnation and wall jet regions of the impinging flows. Results, in general, confirm the superiority of the Reynolds stress transport equation model for predicting mean and fluctuating velocities within the latter regions of such flows. In particular, modifications to the second-moment closure to account for the influence of the surface in distorting the fluctuating pressure field away from the wall successfully predict the damping of normal-to-wall velocity fluctuations throughout the impinging flows. In contrast, results derived from the eddy-viscosity-based approach do not, in general, accurately reproduce experimental observations.     

A Coupled Level Set and Volume of Fluid method for automotive exterior water management applications

Motivated by the need for practical, high fidelity, simulation of water over surface features of road vehicles a Coupled Level Set Volume of Fluid (CLSVOF) method has been implemented into a general purpose CFD code. It has been implemented such that it can be used with unstructured and non-orthogonal meshes. The interface reconstruction step needed for CLSVOF has been implemented using an iterative ‘clipping and capping’ algorithm for arbitrary cell shapes and a re-initialisation algorithm suitable for unstructured meshes is also presented. Successful verification tests of interface capturing on orthogonal and tetrahedral meshes are presented. Two macroscopic contact angle models have been implemented and the method is seen to give very good agreement with experimental data for a droplet impinging on a flat plate for both orthogonal and non-orthogonal meshes. Finally the flow of a droplet over a round edged channel is simulated in order to demonstrate the ability of the method developed to simulate surface flows over the sort of curved geometry that makes the use of a non-orthogonal grid desirable.     

Numerical study on the performance of nozzle flow for supersonic chemical oxygen–iodine lasers

Laser performance is greatly dependent on its operating conditions due to the strong coupling among multi- physics such as gas-dynamics, chemical reaction kinetics and optics in the mixing nozzle of COIL. In this paper, 3D CFD technology is used to simulate the mixing and reactive flow of subsonic cross jet scheme at different conditions. Results obtained show that the jet penetration depth plays a dominant role in the spatial distribution of small signal gains. In the case of over-penetration, unsteady flow structures are induced by impinging between the opposing jets. The optimum spatial distribution of the chemical performance cannot be obtained even if the full penetration condition is achieved through the subsonic transverse jet mixing scheme in the COIL nozzle flow.     

A dissipation‐free numerical method to solve one‐dimensional hyperbolic flow equations

In this paper, a numerical method to capture the shock wave propagation in 1‐dimensional fluid flow problems with 0 numerical dissipation is presented. Instead of using a traditional discrete grid, the new numerical method is built on a range‐discrete grid, which is obtained by a direct subdivision of values around the shock area. The range discrete grid consists of 2 types: continuous points and shock points. Numerical solution is achieved by tracking characteristics and shocks for the movements of continuous and shock points, respectively. Shocks can be generated or eliminated when triggering entropy conditions in a marking step. The method is conservative and total variation diminishing. We apply this new method to several examples, including solving Burgers equation for aerodynamics, Buckley‐Leverett equation for fractional flow in porous media, and the classical traffic flow. The solutions were verified against analytical solutions under simple conditions. Comparisons with several other traditional methods showed that the new method achieves a higher accuracy in capturing the shock while using much less grid number. The new method can serve as a fast tool to assess the shock wave propagation in various flow problems with good accuracy.