Swirling-strength based large eddy simulation of turbulent flow around single square cylinder at low Reynolds numbers

In view of the fact that large scale vortices play the substantial role of momentum transport in turbulent flows, large eddy simulation(LES) is considered as a better simulation model. However, the sub-grid scale(SGS) models reported so far have not ascertained under what flow conditions the LES can lapse into the direct numerical simulation. To overcome this discrepancy, this paper develops a swirling strength based the SGS model to properly model the turbulence intermittency, with the primary characteristics that when the local swirling strength is zero, the local sub-grid viscosity will be vanished. In this paper, the model is used to investigate the flow characteristics of zero-incident incompressible turbulent flows around a single square cylinder(SC)at a low Reynolds number range Re ∈ [103, 104]. The flow characteristics investigated include the Reynolds number dependence of lift and drag coefficients, the distributions of time-spanwise averaged variables such as the sub-grid viscosity and the logarithm of Kolmogorov micro-scale to the base of 10 at Re = 2 500 and 104, the contours of spanwise and streamwise vorticity components at t = 170. It is revealed that the peak value of sub-grid viscosity ratio and its root mean square(RMS) values grow with the Reynolds number. The dissipation rate of turbulent kinetic energy is larger near the SC solid walls.The instantaneous factor of swirling strength intermittency(FSI) exhibits some laminated structure involved with vortex shedding.     

k-ω SST湍流模式三维激波分离流动修正

“激波?边界层分离”是航空气动领域的典型湍流非平衡流动问题, 准确模拟激波分离对于跨声速飞行器气动性能评估和优化设计具有重要意义. 然而传统涡黏性湍流模式中涡黏性系数的定义方式并不适用于非平衡流动, k-ω SST湍流模式为此引入的Bradshaw假设在应用于三维强逆压梯度和较大分离流动时反而限制了雷诺应力的生成, 导致包括k-ω SST在内的常用涡黏性湍流模式均无法对此类流动进行准确模拟. 同时, 现有的非线性雷诺应力本构关系也并不能有效提高模拟精度. 为此, 针对k-ω SST模式分别提出了基于Bradshaw假设和基于长度尺度的两种激波分离流动修正方法. 前者通过提高Bradshaw常数的方式放宽了对雷诺应力生成的限制, 后者则从湍流长度尺度概念出发, 利用混合长度理论、湍动能生成/耗散之比和一种新定义的长度尺度之比构造了ω方程耗散项修正函数, 提高了模式在三维激波分离流动中的建模长度尺度. 两种方法对ONERA M6机翼跨声速大攻角流动均能得到较雷诺应力模式更好的模拟结果. 进一步的雷诺应力分析表明, 三维激波分离流动中“主雷诺应力分量”的概念不再成立, 各雷诺应力分量大小接近. 网格收敛性分析、对其他攻角状态的验证以及湍流平板边界层壁面律验证进一步确认了所提出的两种修正方法的合理性、有效性和通用性.      

超临界翼型跨音速特性实验研究

针对新设计的超临界翼型,采用风洞实验方法验证和评估了其气动特性。在增压连续式跨音速风洞(NF-6风洞)开展了超临界翼型跨音速特性的实验研究,验证了该翼型设计的压力分布曲线特点。激波位置和波后压力平台区长度表明设计结果和实验结果基本一致,揭示了超临界翼型跨音速的气动特性;阻力发散马赫数达到期望的设计指标,探讨了雷诺数对该翼型气动特性的影响。最后采用升华法实现了翼型表面流动特性的显示。结果表明转捩点约在16%弦长位置。     

Forced oscillations of airfoil flows

The effect of background flow oscillations on a transonic airfoil (NACA 0012) flow was investigated experimentally. The oscillations were generated by means of a rotating plate placed downstream of the airfoil. Owing to the expansion and compression waves generated at the plate, the flow over the airfoil flow was drastically disturbed. This resulted in the presence of high intensity oscillations of a shock wave and a separation bubble on the suction surface of the airfoil. For relatively large values of the airfoil angle of attack, weak shock waves (transonic sound waves) were periodically shed upstream of the airfoil.This work was supported by Commission of the European Communities (Communit's Action for Cooperation in Science and Technology with Central and Eastern European Countries).The authors wish to thank Mr P. Koperski for his effective assistance in taking the photographs.     

Buffeting in transonic flow prediction using time‐dependent turbulence model

In transonic flow conditions, the shock wave/turbulent boundary layer interaction and flow separations on wing upper surface induce flow instabilities, ‘buffet’, and then the buffeting (structure vibrations). This phenomenon can greatly influence the aerodynamic performance. These flow excitations are self‐sustained and lead to a surface effort due to pressure fluctuations. They can produce enough energy to excite the structure. The objective of the present work is to predict this unsteady phenomenon correctly by using unsteady Navier–Stokes‐averaged equations with a time‐dependent turbulence model based on the suitable (k–ε) turbulent eddy viscosity model. The model used is based on the turbulent viscosity concept where the turbulent viscosity coefficient (Cμ) is related to local deformation and rotation rates. To validate this model, flow over a flat plate at Mach number of 0.6 is first computed, then the flow around a NACA0012 airfoil. The comparison with the analytical and experimental results shows a good agreement. The ONERA OAT15A transonic airfoil was chosen to describe buffeting phenomena. Numerical simulations are done by using a Navier–Stokes SUPG (streamline upwind Petrov–Galerkin) finite‐element solver. Computational results show the ability of the present model to predict physical phenomena of the flow oscillations. The unsteady shock wave/boundary layer interaction is described. Copyright © 2005 John Wiley & Sons, Ltd.     

High-order strand grid methods for shock turbulence interaction

ABSTRACT

In this work, we examine the flux correction method for three-dimensional transonic turbulent flows on strand grids. Building upon previous work, we treat flux derivatives along strands with high-order summation-by-parts operators and penalty-based boundary conditions. A finite-volume like limiting strategy is implemented in the flux correction algorithm in order to sharply capture shocks. To achieve turbulence closure in the Reynolds-Averaged Navier–Stokes equations, a robust version of the Spalart–Allmaras turbulence model is employed that accommodates negative values of the turbulence working variable. Validation studies are considered which demonstrate the flux correction method achieves a high degree of accuracy for turbulent shock interaction flows.     

The application of the gradient-based adjoint multi-point optimization of single and double shock control bumps for transonic airfoils

A shock control bump (SCB) is a flow control method that uses local small deformations in a flexible wing surface to considerably reduce the strength of shock waves and the resulting wave drag in transonic flows. Most of the reported research is devoted to optimization in a single flow condition. Here, we have used a multi-point adjoint optimization scheme to optimize shape and location of the SCB. Practically, this introduces transonic airfoils equipped with the SCB that are simultaneously optimized for different off-design transonic flight conditions. Here, we use this optimization algorithm to enhance and optimize the performance of SCBs in two benchmark airfoils, i.e., RAE-2822 and NACA-64-A010, over a wide range of off-design Mach numbers. All results are compared with the usual single-point optimization. We use numerical simulation of the turbulent viscous flow and a gradient-based adjoint algorithm to find the optimum location and shape of the SCB. We show that the application of SCBs may increase the aerodynamic performance of an RAE-2822 airfoil by 21.9 and by 22.8 % for a NACA-64-A010 airfoil compared to the no-bump design in a particular flight condition. We have also investigated the simultaneous usage of two bumps for the upper and the lower surfaces of the airfoil. This has resulted in a 26.1 % improvement for the RAE-2822 compared to the clean airfoil in one flight condition.     

A SOLVER USING NEWTON‘S METHOD FOR UNSTEADY VISCOUS FLOWS

A solver is developed for time-accurate computations of viscous flows based on the conception of Newton‘s method. A set of pseudo-time derivatives are added into governing equations and the discretized system is solved using GMRES algorithm. Due to some special properties of GMRES algorithm, the solution procedure for unsteady flows could be regarded as a kind of Newton iteration. The physical-time derivatives of governing equations are discretized using two different approaches, I.e., 3-point Euler backward, and Crank-Nicolson formulas, both with 2nd-order accuracy in time but with different truncation errors. The turbulent eddy viscosity is calculated by using a version of Spalart~Allmaras one-equation model modified by authors for turbulent flows. Two cases of unsteady viscous flow are investigated to validate and assess the solver, I.e., low Reynolds number flow around a row of cylinders and transonic bi-circular-arc airfoil flow featuring the vortex shedding and shock buffeting problems, respectively. Meanwhile, comparisons between the two schemes of timederivative discretizations are carefully made. It is illustrated that the developed unsteady flow solver shows a considerable efficiency and the Crank-Nicolson scheme gives better results compared with Euler method.     

Fourier time spectral method for subsonic and transonic flows

The time accuracy of the exponentially accu-rate Fourier time spectral method (TSM) is examined and compared with a conventional 2nd-order backward differ-ence formula (BDF) method for periodic unsteady flows. In particular, detailed error analysis based on numerical com-putations is performed on the accuracy of resolving the local pressure coefficient and global integrated force coefficients for smooth subsonic and non-smooth transonic flows with moving shock waves on a pitching airfoil. For smooth sub-sonic flows, the Fourier TSM method offers a significant accuracy advantage over the BDF method for the predic-tion of both the local pressure coefficient and integrated force coefficients. For transonic flows where the motion of the discontinuous shock wave contributes significant higher-order harmonic contents to the local pressure fluctuations, a sufficient number of modes must be included before the Fourier TSM provides an advantage over the BDF method. The Fourier TSM, however, still offers better accuracy than the BDF method for integrated force coefficients even for transonic flows. A problem of non-symmetric solutions for symmetric periodic flows due to the use of odd numbers of intervals is uncovered and analyzed. A frequency-searching method is proposed for problems where the frequency is not known a priori. The method is tested on the vortex shedding problem of the flow over a circular cylinder.     

Numerical simulation of viscous transonic flows over an airfoil

High Reynolds number viscous transonic flow is described based on an interaction of the potential outer flow with the boundary layer and wake. Following the procedure of Lighthill (1958), the solutions in these domains are matched to each other through boundary conditions. The solution to the complete problem is obtained iteratively through successive computations of the flows in the outer and inner domains. Both old and new algorithms are used for the iteration process and subsequent problem solution. Results are given for all the airfoils from the Experimental Data Base for Computer Program Assessment (AGARD-AR-138, 1979). A comparison of these results with experimental data shows the degree of agreement between these unbounded airfoil flow simulations and real transonic flow over the central part of a straight wing.     

Experiences with two detached eddy simulation approaches through use of a high-order finite volume solver

In the present study, two advanced detached eddy simulation (DES) approaches, shear-layer-adapted delayed DES and zonal DES in mode II, which are known to help transition from RANS to LES mode, are employed in various flow problems in conjunction with a high-order finite volume solver. The numerical scheme, being only applicable on structured grids, has low-dissipation and low-dispersion features. Such features benefit mostly in the LES mode, minimizing the interference of numerical diffusion with subgrid eddy viscosity. First, corresponding subgrid models are validated via decaying homogeneous turbulence benchmark case. Then, a channel flow problem is chosen to examine these models in attached flow situations. Finally, flow around an airfoil at low Reynolds number is solved using the shear-layer-adapted delayed DES approach only, in an aim to obtain trailing-edge noise spectrum at an observer location. Despite some log-layer mismatch over turbulent boundary layers, which is typical of most DES methods, the combined application of high-resolution numerical method and advanced DES approaches, which are implemented on a stabilized Spalart-Allmaras turbulence model, shows merit in resolution of turbulence in regions of interest.     

Frequency lock-in phenomenon for oscillating airfoils in buffeting flows

Navier-Stokes based computer simulations are conducted to determine the aerodynamic flow field response that is observed for a NACA0012 airfoil that undergoes prescribed harmonic oscillation in transonic buffeting flows, and also in pre-buffet flow conditions. Shock buffet is the term for the self-sustained shock oscillations that are observed for certain combinations of Mach number and steady mean flow angle of attack even in the absence of structural motion. The shock buffet frequencies are typically on the order of the elastic structural frequencies, and therefore may be a contributor to transonic aeroelastic response phenomena, including limit-cycle oscillations. Numerical simulations indicate that the pre-shock-buffet flow natural frequency increases with mean angle of attack, while the flow damping decreases and approaches zero at the onset of buffet. Airfoil harmonic heave motions are prescribed to study the interaction between the flow fields induced by the shock buffet and airfoil motion, respectively. At pre-shock-buffet conditions the flow response is predominantly at the airfoil motion frequency, with some smaller response at multiplies of this frequency. At shock buffet conditions, a key effect of prescribed airfoil motions on the buffeting flow is to create the possibility of a lock-in phenomenon, in which the shock buffet frequency is synchronized to the prescribed airfoil motion frequency for certain combinations of airfoil motion frequencies and amplitudes. Aerodynamic gain-phase models for the lock-in region, as well as for the pre-shock-buffet conditions are suggested, and also a possible relationship between the lock-in mechanism and limit-cycle oscillation is discussed.     

A two‐scale low‐Reynolds number turbulence model

In this study, a two‐scale low‐Reynolds number turbulence model is proposed. The Kolmogorov turbulence time scale, based on fluid kinematic viscosity and the dissipation rate of turbulent kinetic energy (ν, ε), is adopted to address the viscous effects and the rapid increasing of dissipation rate in the near‐wall region. As a wall is approached, the turbulence time scale transits smoothly from a turbulent kinetic energy based (k, ε) scale to a (ν, ε) scale. The damping functions of the low‐Reynolds number models can thus be simplified and the near‐wall turbulence characteristics, such as the ε distribution, are correctly reproduced. The proposed two‐scale low‐Reynolds number turbulence model is first examined in detail by predicting a two‐dimensional channel flow, and then it is applied to predict a backward‐facing step flow. Numerical results are compared with the direct numerical simulation (DNS) budgets, experimental data and the model results of Chien, and Lam and Bremhorst respectively. It is proved that the proposed two‐scale model indeed improves the predictions of the turbulent flows considered. Copyright © 2000 John Wiley & Sons, Ltd.     

Performance Improvement of a Supercritical Airfoil by a Multi-point Optimized Shock Control Channel

A shock control channel (SCC) is a flow control method introduced here to control the shock wave/boundarylayer interaction (SWBLI) in order to reduce the resulting wave drag in transonic flows. An SCC transfers an appropriate amount of mass and momentum from downstream of the shock wave location to its upstream to decrease the pressure gradient across the shock wave and as a result the shock-wave strength is reduced. Here, a multi-point optimization method under a constant-lift-coefficient constraint is used to find the optimum design of the SCC. This flow control method is implemented on a RAE-2822 supercritical airfoil for a wide range of off-design transonic Mach numbers. The RANS flow equations are solved using Roe’s averages scheme and a gradient-based adjoint algorithm is used to find the optimum location and shape of the SCC. The solver is validated against experimental works studying effect of cavities in transonic airfoils. It is shown that the application of an SCC improves the average aerodynamic efficiency in a range of off-design conditions by 13.2% in comparison with the original airfoil. The SCC is shown to be an effective tool also for higher angle of attack at transonic flows. We have also studied the SWBLI and how the optimization algorithm makes the flow wave structure and interactions of the shock wave with the boundary layer favorable.     

Semi-hyperbolic Patches of Solutions to the Two-dimensional Euler Equations

We construct semi-hyperbolic patches of solutions, in which one family out of two nonlinear families of characteristics starts on sonic curves and ends on transonic shock waves, to the two-dimensional Euler equations. This type of solution appears in the transonic flow over an airfoil and Guderley reflection, and is common in the numerical solutions of Riemann problems.     

LES of incompressible heat and fluid flows past a square cylinder at high Reynolds numbers

This paper presents large eddy simulation (LES) results of incompressible heat and fluid flows around a square cylinder (SC) at zero incident angle at high Reynolds numbers (Re) in the range from 1.25×10⁵ to 3.5×10⁵. LES results are obtained on the basis of swirling strength based sub-grid model, and a higher order upwind scheme developed with respect to the Taylor expansion. It was found that, for the zero incident SC wake flows at a Reynolds number in the range {Re₅ = Re/10⁵ ∈ [1.25, 3.5]}, the Strouhal number equals to 0.1079, completely independent of the Reynolds number; the coefficient of drag is around 1.835 with an uncertainty of about 1.9%, almost non-sensitive to the Re. When Re is beyond 3.0×10⁵, the time-averaged peak value of sub-grid viscosity is over 340, implying that the role of sub-grid model is crucial in some regions where vortex motion is active and vortex interaction is intense. The time–spanwise (t-z) averaged sub-grid viscosity ratio profiles and the profiles of fluctuations of the sub-grid viscosity ratio and velocity components at four locations downstream of the SC are presented. The fields of the t-z averaged sub-grid viscosity ratio, and the instantaneous fields of streamwise and spanwise vorticities are also reported and discussed. The predicted mean Nusselt number is compared with empirical correlations, revealing that swirling strength based LES has its potential in predicting natural and industrial flows.     

Comparative study of implicit and subgrid‐scale model large‐eddy simulation techniques for low‐Reynolds number airfoil applications

A computational study of a high‐fidelity, implicit large‐eddy simulation (ILES) technique with and without the use of the dynamic Smagorinsky subgrid‐scale (SGS) model is conducted to examine the contributions of the SGS model on solutions of transitional flow over the SD7003 airfoil section. ILES without an SGS model has been shown in the past to produce comparable and sometimes favorable results to traditional SGS‐based large‐eddy simulation (LES) when applied to canonical turbulent flows. This paper evaluates the necessity of the SGS model for low‐Reynolds number airfoil applications to affirm the use of ILES without SGS‐modeling for a broader class of problems such as those pertaining to micro air vehicles and low‐pressure turbines. It is determined that the addition of the dynamic Smagorinsky model does not significantly affect the time‐mean flow or statistical quantities measured around the airfoil section for the spatial resolutions and Reynolds numbers examined in this study. Additionally, the robustness and reduced computational cost of ILES without the SGS model demonstrates the attractiveness of ILES as an alternative to traditional LES. Published 2012. This article is a US Government work and is in the public domain in the USA.     

Time‐dependent turbulent cavitating flow computations with interfacial transport and filter‐based models

Turbulent cavitating flow computations need to address both cavitation and turbulence modelling issues. A recently developed interfacial dynamics‐based cavitation model (IDCM) incorporates the interfacial transport into the computational modelling of cavitation dynamics. For time‐dependent flows, it is known that the engineering turbulence closure such as the original k–ε model often over‐predicts the eddy viscosity values reducing the unsteadiness. A recently proposed filter‐based modification has shown that it can effectively modulate the eddy viscosity, rendering better simulation capabilities for time‐dependent flow computations in term of the unsteady characteristics. In the present study, the IDCM along with the filter‐based k–ε turbulence model is adopted to simulate 2‐D cavitating flows over the Clark‐Y airfoil. The chord Reynolds number is Re=7.0 × 10⁵. Two angles‐of‐attack of 5 and 8° associated with several cavitation numbers covering different flow regimes are conducted. The simulation results are assessed with the experimental data including lift, drag and velocity profiles. The interplay between cavitation and turbulence models reveals substantial differences in time‐dependent flow results even though the time‐averaged characteristics are similar. Copyright © 2005 John Wiley & Sons, Ltd.     

Turbulence model and numerical scheme assessment for buffet computations

The prediction of shock‐induced oscillations over transonic rigid airfoils is important for a better understanding of the buffeting phenomenon. The unsteady resolution of the Navier–Stokes equations is performed with various transport‐equation turbulence models in which corrections are added for non‐equilibrium flows. The lack of numerical efficiency due to the CFL stability condition is circumvented by the use of a wall law approach and a dual time stepping method. Moreover, various numerical schemes are used to try and be independent of the numerical discretization. Comparisons are made with the experimental results obtained for the supercritical RA16SC1 airfoil. They show the interest in using the SST correction or realizability conditions to get correct predictions of the frequency, amplitude and pressure fluctuations over the airfoil. Copyright © 2004 John Wiley & Sons, Ltd.