高速流场气动光学RANS/DSMC混合数值模拟算法研究

提出了适用于高速流场气动光学数值模拟的RANS/DSMC混合算法.通过RANS对全局时均流场进行数值模拟,再对其中局部流场的脉动量采用DSMC进行数值模拟,以Maxwellian速度分布实现RANS宏观量信息向DSMC微观量信息的传递.采用超声速环境下尖劈模型对该混合算法进行校验,对比试验结果论证了算法的正确性.在计算...     

N-S方程数值研究翼型对微型扑翼气动特性的影响

首先基于嵌套网格发展了一套适用于三维扑翼研究的非定常雷诺平均Navier-Stokes(RANS)方程数值模拟方法.为了解决微型扑翼在低马赫数下的收敛问题,使用了预处理方法,湍流模型为BL模型.在该方法的基础上,保持状态参数和扑翼表面形状一定的情况下,分别研究了一系列不同厚度、不同弯度的翼型对于微型扑翼气动特性的影响....     

Scale resolving simulations of a high-temperature turbulent jet in a cold crossflow: Comparison of two approaches

A high-temperature turbulent jet in a cold crossflow is investigated with the help of two scale-resolving simulation approaches. This work aims at improving the methodologies used to predict the thermal footprint of exhaust gases issuing from helicopter engines onto the fuselage. Specific attention is brought to the capability of scale resolving simulations to correctly reproduce flow dynamics and turbulent mixing. Mean flow features, turbulent quantities and temperature fields are compared and validated against wind tunnel test measurements. In addition, the present work highlights the importance of synthetic turbulence injection at pipe inlet to obtain a fair prediction of both flow dynamics and temperature field.     

Investigation of conditional source-term estimation applied to a non-premixed turbulent flame

A turbulent combustion model, Conditional Source-term Estimation (CSE) is applied to a non-premixed turbulent jet methane flame. The conditional chemical source terms are determined on the basis of first order closure and the conditional averaged species concentrations are obtained by inverting an integral equation. The Tikhonov method is implemented for regularisation. Detailed chemistry is tabulated using the trajectory generated low-dimensional manifold method. Radiation due to the gaseous species is included. Reynolds Averaged Navier–Stokes calculations are performed using two different turbulence models. The objectives of the paper are (i) assessment of the impact of the main numerical parameters in CSE and (ii) comparison of the CSE numerical predictions with available experimental data and results from previous simulations for the selected flame. The number of CSE domains and the number of points in each CSE domain are shown to have a significant impact on the results if not selected appropriately. The present CSE calculations always converge to unique and stable predictions. The corrected k–ε model yields mixture fraction profiles in good agreement with the experimental data values for axial locations in the first half of the flame. Farther downstream, the RNG k–ε model performs better. Overall, the current predictions for the mixture fraction are in good agreement with the experimental data. The predicted temperatures using CSE and the k–ε turbulence model with a modified value of C_ε1 = 1.47 are found to be in very good agreement with the experimental data. Further, the current CSE results are of comparable quality with previous simulations using the flamelet model and conditional moment closure. Future work may include further investigation on optimal determination of the regularisation parameter and alternative regularisation techniques, soot modelling within the CSE formulation, and improved formulation of radiation.     

LES of spatially developing turbulent boundary layer over a concave surface

We revisit the problem of a spatially developing turbulent boundary layer over a concave surface. Unlike previous investigations, we simulate the combined effects of streamline curvature as well as curvature-induced pressure gradients on the turbulence. Our focus is on investigating the response of the turbulent boundary layer to the sudden onset of curvature and the destabilising influence of concave surface in the presence of pressure gradients. This is of interest for evaluating the turbulence closure models. At the beginning of the curve, the momentum thickness Reynolds number is 1520 and the ratio of the boundary layer thickness to the radius of curvature is δ₀/R = 0.055. The radial profiles of the mean velocity and turbulence statistics at different locations along the concave surface are presented. Our recently proposed curvature-corrected Reynolds Averaged Navier-Stokes (RANS) model is assessed in an a posteriori sense and the improvements obtained over the base model are reported. From the large Eddy simulation (LES) results, it was found that the maximum influence of concave curvature is on the wall-normal component of the Reynolds stress. The budgets of wall-normal Reynolds stress also confirmed this observation. At the onset of curvature, the effect of adverse pressure gradient is found to be predominant. This decreases the skin friction levels below that in the flat section.     

Nonlinear corrector for Reynolds-averaged Navier-Stokes equations

The scope of this paper is to present a nonlinear error estimation and correction for Navier-Stokes and Reynolds-averaged Navier-Stokes equations. This nonlinear corrector enables better solution or functional output predictions at fixed mesh complexity and can be considered in a mesh adaptation process. After solving the problem at hand, a corrected solution is obtained by solving again the problem with an added source term. This source term is deduced from the evaluation of the residual of the numerical solution interpolated on the h/2 mesh. To avoid the generation of the h/2 mesh (which is prohibitive for realistic applications), the residual at each vertex is computed by local refinement only in the neighborhood of the considered vertex. One of the main feature of this approach is that it automatically takes into account all the properties of the considered numerical method. The numerical examples point out that it successfully improves solution predictions and yields a sharp estimate of the numerical error. Moreover, we demonstrate the superiority of the nonlinear corrector with respect to linear corrector that can be found in the literature.     

A divergence-conforming hybridized discontinuous Galerkin method for the incompressible Reynolds-averaged Navier-Stokes equations

We present a hybridized discontinuous Galerkin (HDG) method for the incompressible Reynolds-averaged Navier-Stokes equations coupled with the Spalart-Allmaras one-equation turbulence model. The method extends upon an HDG method recently introduced by Rhebergen and Wells for the incompressible Navier-Stokes equations. With a special choice of velocity and pressure spaces for both element and trace degrees of freedom (DOFs), the method returns pointwise divergence-free mean velocity fields and properly balances momentum and energy. We further examine the use of different polynomial degrees and meshes to see how the order of the scalar eddy viscosity affects the convergence of the mean velocity and pressure fields, specifically for the method of manufactured solutions. As is standard with HDG methods, static condensation can be employed to remove the element DOFs and thus dramatically reduce the global number of DOFs. Numerical results illustrate the effectiveness of the proposed methodology.     

Zonal PANS: evaluation of different treatments of the RANS–LES interface

The partially Reynolds-averaged Navier–Stokes (PANS) model can be used to simulate turbulent flows either as RANS, large eddy simulation (LES) or DNS. Its main parameter is f_k whose physical meaning is the ratio of the modelled to the total turbulent kinetic energy. In RANS f_k = 1, in DNS f_k = 0 and in LES f_k takes values between 0 and 1. Three different ways of prescribing f_k are evaluated for decaying grid turbulence and fully developed channel flow: f_k = 0.4, f_k = k^3/2 _tot/? and, from its definition, f_k = k/k_tot where k_tot is the sum of the modelled, k, and resolved, k_res, turbulent kinetic energy. It is found that the f_k = 0.4 gives the best results. In Girimaji and Wallin, a method was proposed to include the effect of the gradient of f_k. This approach is used at RANS– LES interface in the present study. Four different interface models are evaluated in fully developed channel flow and embedded LES of channel flow: in both cases, PANS is used as a zonal model with f_k = 1 in the unsteady RANS (URANS) region and f_k = 0.4 in the LES region. In fully developed channel flow, the RANS– LES interface is parallel to the wall (horizontal) and in embedded LES, it is parallel to the inlet (vertical). The importance of the location of the horizontal interface in fully developed channel flow is also investigated. It is found that the location – and the choice of the treatment at the interface – may be critical at low Reynolds number or if the interface is placed too close to the wall. The reason is that the modelled turbulent shear stress at the interface is large and hence the relative strength of the resolved turbulence is small. In RANS, the turbulent viscosity – and consequently also the modelled Reynolds shear stress – is only weakly dependent on Reynolds number. It is found in the present work that it also applies in the URANS region.     

Simulation of vortex shedding behind a bluff body flame stabilizer using a hybrid U-RANS/PDF method

The present study aims at the investigation of the effects of turbulence-chemistry interaction on combustion instabilities using a probability density function(PDF) method.The instantaneous quantities in the flow field were decomposed into the Favre-averaged variables and the stochastic fluctuations,which were calculated by unsteady Reynolds averaged Navier-Stokes(U-RANS) equations and the PDF model,respectively.A joint fluctuating velocityfrequency-composition PDF was used.The governing equations are solved by a consistent hybrid finite volume/MonteCarlo algorithm on triangular unstructured meshes.A nonreacting flow behind a triangular-shaped bluff body flame stabilizer in a rectilinear combustor was simulated by the present method.The results demonstrate the capability of the present method to capture the large-scale coherent structures.The triple decomposition was performed,by dividing the coherent Favre-averaged velocity into time-averaged value and periodical coherent part,to analyze the coherent and incoherent contributions to Reynolds stresses.A simple modification to the coefficients in the turbulent frequency model will help to improve the simulation results.Unsteady flow fields were depicted by streamlines and vorticity contours.Moreover,the association between turbulence production and vorticity saddle points is illustrated.