Computations of strongly swirling flows with second‐moment closures

The present study is concerned with simulating turbulent, strongly swirling flows by eddy viscosity model and Reynolds stress transport model variants adopting linear and quadratic form of the pressure–strain models. Flows with different inlet swirl numbers, 2.25 and 0.85, were investigated. Detailed comparisons of the predicted results and measurements were presented to assess the merits of model variants. For the swirl number 2.25 case, due to the inherent capability of the Reynolds stress models to capture the strong swirl and turbulence interaction, both the linear and quadratic form of the pressure–strain models predict the flow adequately. In strong contrast, the k–ϵ model predicts an excessively diffusive flow fields. For the swirl number 0.85 case, both the k–ϵ and Reynolds stress model with linear pressure–strain process, show an excessive diffusive transport of the flow fields. The quadratic pressure–strain model, on the other hand, mimics the correct flow development with the recirculating region being correctly predicted. Copyright © 1999 John Wiley & Sons, Ltd.     

可压缩多介质粘性流体和湍流的大涡模拟

在可压缩多介质粘性流体动力学高精度计算方法MVPPM(multi-viscous-fluid piecewise parabolicmethod)基础上,引入Smagorinsky和Vreman亚格子湍流模型,采用大涡数值模拟方法求解可压缩粘性流体NS(Navier-Stokes)方程,给出适用于可压缩多介质流体界面不稳定性发展演化至湍流阶段的计算方法和二维计算程序MVFT(multi-viscosity-fluid and turbulence)。在2种亚格子湍流模型下计算了LANL(Los Ala-mos National Laboratory)激波管单气柱RM不稳定性实验,分析了气柱的形状、流场速度以及涡的特征,通过与LANL实验和计算结果的比较可知,Vreman模型略优于Smagorinsky模型,MVFT方法和计算程序可用于对界面不稳定性发展演化至湍流阶段的数值模拟。     

Residual‐based artificial viscosity for simulation of turbulent compressible flow using adaptive finite element methods

In this paper, we present a finite element method with a residual‐based artificial viscosity for simulation of turbulent compressible flow, with adaptive mesh refinement based on a posteriori error estimation with sensitivity information from an associated dual problem. The artificial viscosity acts as a numerical stabilization, as shock capturing, and as turbulence capturing for large eddy simulation of turbulent flow. The adaptive method resolves parts of the flow indicated by the a posteriori error estimates but leaves shocks and turbulence under‐resolved in a large eddy simulation. The method is tested for examples in 2D and 3D and is validated against experimental data. Copyright © 2012 John Wiley & Sons, Ltd.     

Subgrid-scale modelling at low mesh reynolds number

Subgrid-scale models are derived for large-eddy simulations in the limit of low mesh Reynolds number, or, equivalently, resolution approaching that required for full resolution of the simulated turbulent flow. The models are constructed from standard forms of the dissipation spectrum in a manner analogous to that used to derive the classical Smagorinsky-Lilly model from the inertial range spectrum. Practical methods for computing the subgrid-scale eddy viscosity are described, together with examples of the effects of using such models in a real simulation.     

高雷诺数湍流风场大涡模拟的并行直接求解方法

在环境流体力学中,风场是风沙流、风雪流等自然环境特性问题研究的动力源和基础. 通常采用壁湍流模型进行风场大涡模拟(large eddy simulation, LES)计算,但受到计算规模的限制使得 高雷诺数风场的模拟计算难以实现. 并行计算技术是解决大规模高雷诺数风场大涡模拟的关键技术之一. 在不可压湍流风场的LES模拟中,压力泊松方程的并行计算技术是进行规模并行计算的困难点. 根据风场流动模拟计算的特点,采用水平网格等距而垂直于地面网格非等距,在解决规模并行计算中求解压力泊松方程的难点问题时,利用FFT解耦三维泊松方程使其变为垂向的一维三对角方程, 并利用可并行的三对角方程PDD求解技术,可建立三维泊松方程的直接并行求解技术. 结合其它容易并行的动量方程计算,本文建立风场LES模拟的并行直接求解方法(parallel direct method-LES, PDM-LES). 在超级计算机上对新方法进行并行效率测试,并行计算效率达到90${\%}$. 新的方法可用于进行湍流风场大涡模拟的大规模并行计算. 计算结果表明,湍流风场瞬时速度分布近壁面存在条带状的拟序结构,平均场的速度分布符合速度对数律特性,风场湍流特性基本合理.      

Assessment of the SST-IDDES with a shear-layer-adapted subgrid length scale for attached and separated flows

A shear-layer-adapted subgrid length scale is applied to the improved delayed detached eddy simulation using the shear stress transport background model (SST-IDDES). The aim is to assess the combination of the wall-modeled LES (WMLES) branch of the SST-IDDES with the new length scale in computing attached flows, as well as to assess the effect of the new length scale when it is applied to the SST-IDDES for mitigating the “grey area” issue through initiating a dramatic drop of eddy viscosity in the initial region of a free shear layer. The assessment is conducted through simulations of a turbulent boundary layer, a fully developed channel flow, a near-sonic turbulent jet and a backward-facing step flow. The results provide strong evidence for the conclusion that the SST-IDDES combined with the new length scale performs the same as the original SST-IDDES when its WMLES branch is applied to compute the resolved parts of an attached flow, and the combination helps mitigate the “grey area” issue of the SST-IDDES and accurately represent the K-H instability in the initial region of a free shear layer. In addition, the superiority is particularly remarkable for the simulations with coarse grids.     

剪切湍流的涡黏张量模式

对剪切湍流提出了涡黏系数为四阶张量的涡黏张量模式。引入近代数学中Ｍｏｏｒｅ－Ｐｅｎｒｏｓｅ广义逆矩阵的研究结果，给出了构造涡黏张量各分量的计算公式。用平面后台阶流动验算了剪切湍流的涡黏张量模式，比ＲＳＭ和ｋ－ε模式更接近实验结果。提出了剪切湍流涡黏张量模式的应用设想。     

End-to-end differentiable learning of turbulence models from indirect observations

The emerging push of the differentiable programming paradigm in scientific computing is conducive to training deep learning turbulence models using indirect observations. This paper demonstrates the viability of this approach and presents an end-to-end differentiable framework for training deep neural networks to learn eddy viscosity models from indirect observations derived from the velocity and pressure fields. The framework consists of a Reynolds-averaged Navier–Stokes(RANS) solver and a neuralnetwork-represented turbulence model, each accompanied by its derivative computations. For computing the sensitivities of the indirect observations to the Reynolds stress field, we use the continuous adjoint equations for the RANS equations, while the gradient of the neural network is obtained via its built-in automatic differentiation capability. We demonstrate the ability of this approach to learn the true underlying turbulence closure when one exists by training models using synthetic velocity data from linear and nonlinear closures. We also train a linear eddy viscosity model using synthetic velocity measurements from direct numerical simulations of the Navier–Stokes equations for which no true underlying linear closure exists. The trained deep-neural-network turbulence model showed predictive capability on similar flows.     

Balance equation of the generalised sub-grid scale (SGS) turbulent kinetic energy in a new tensorial dynamic mixed SGS model

A new dynamic model is proposed in which the eddy viscosity is defined as a symmetric second rank tensor, proportional to the product of a turbulent length scale with an ellipsoid of turbulent velocity scales. The employed definition of the eddy viscosity allows to remove the local balance assumption of the SGS turbulent kinetic energy formulated in all the dynamic Smagorinsky-type SGS models. Furthermore, because of the tensorial structure of the eddy viscosity the alignment assumption between the principal axes of the SGS turbulent stress tensor and the resolved strain-rate tensor is equally removed, an assumption which is employed in the scalar eddy viscosity SGS models. The proposed model is tested for a turbulent channel flow. Comparison with the results obtained with other dynamic SGS models (Dynamic Smagorinsky Model, Dynamic Mixed Model and Dynamic K-equation Model) shows that the tensorial definition of the eddy viscosity and the removal of the local balance assumption of the SGS turbulent kinetic energy considerably improves the agreement between results obtained with Large Eddy simulation (LES) and Direct Numerical Simulations (DNS), respectevely. Received August 26, 1999     

Prediction of unsteady,separated boundary layer over a blunt body for laminar,turbulent, and transitional flow

The focus of this paper is to study the ability of unsteady RANS‐based CFD to predict separation over a blunt body for a wide range of Reynolds numbers particularly the ability to capture laminar‐to‐turbulent transition. A perfect test case to demonstrate this point is the cylinder‐in‐crossflow for which a comparison between experimental results from the open literature and a series of unsteady simulations is made. Reynolds number based on cylinder diameter is varied from 10⁴ to 10⁷ (subcritical through supercritical flow). Two methods are used to account for the turbulence in the simulations: currently available eddy–viscosity models, including standard and realizable forms of the k–ε model; and a newly developed eddy–viscosity model capable of resolving boundary layer transition, which is absolutely necessary for the type and range of flow under consideration. The new model does not require user input or ‘empirical’ fixes to force transition. For the first time in the open literature, three distinct flow regimes and the drag crisis due to the downstream shift of the separation point are predicted using an eddy–viscosity based model with transition effects. Discrepancies between experimental and computational results are discussed, and difficulties for CFD prediction are highlighted. Copyright © 2004 John Wiley & Sons, Ltd.     

Turbulence modeling for the axially rotating pipe from the viewpoint of analytical closures

A single-point model eddy viscosity model of rotation effects on the turbulent flow in an axially rotating pipe is developed based on two-point closure theories. Rotation is known to impede energy transfer in turbulence; this fact is reflected in the present model through a reduced eddy viscosity, leading to laminarization of the mean velocity profile and return to a laminar friction law in the rapid rotation limit. This model is compared with other proposals including linear redistribution effects through the rapid pressure-strain correlation, Richardson number modification of the eddy viscosity in a model of non-rotating turbulence, and the reduction of turbulence through the suppression of near-wall production mechanisms. PACS 47.27.Eq, 47.32.-y     

Hybrid RANS/LES of flow and heat transfer in round impinging jets

Fluid flow and convective heat transfer predictions are presented of round impinging jets for several combinations of nozzle-plate distances H/D = 2, 6 and 13.5 (where D is the nozzle diameter) and Reynolds numbers Re = 5000, 23,000 and 70,000 with the newest version of the k-ω model of Wilcox (2008) and three hybrid RANS/LES models. In the RANS mode of the hybrid RANS/LES models, the k-ω model is recovered. Three formulations are considered to activate the LES mode. The first model is similar to the hybrid models of Davidson and Peng (2003) and Kok et al. (2004). The turbulent length scale is replaced by the grid size in the destruction term of the k-equation and in the definition of the RANS eddy viscosity. As grid size, a maximum measure of the hexahedral grid cell is used. The second model has the same k-equation, but the eddy viscosity is the minimum of the k-ω eddy viscosity and the Smagorinsky eddy viscosity, following a proposal by Batten et al. (2004). The Smagorinsky eddy viscosity is formed with the cube root of the cell volume. The third model has, again, the same k-equation, but has an eddy viscosity which is an intermediate between the eddy viscosities of the first and second models. This is reached by using the cube root of the cell volume in the eddy viscosity formula of the first model.The simulation results are compared with experimental data for the high Reynolds number cases Re = 23,000 and Re = 70,000 and LES data for the low-Reynolds number case Re = 5000. The Reynolds numbers are defined with the nozzle diameter and the bulk velocity at nozzle outlet. At low nozzle-plate distance (the impingement plate is in the core of the jet), turbulent kinetic energy is overpredicted by RANS in the stagnation flow region. This leads to overprediction of the heat transfer rate along the impingement plate in the impact zone. At high nozzle-plate distance (the impingement plate is in the mixed-out region of the jet), the turbulence mixing is underpredicted by RANS in the shear layer of the jet which gives a too high length of the jet core. This also results in overprediction of the heat transfer rate in the impingement zone caused by too big temperature gradients at impingement.All hybrid RANS/LES models are able to correct the heat transfer overprediction of the RANS model. For good predictions at low nozzle-plate distance, it is necessary to sufficiently resolve the formation and development of the near-wall vortices in the jet impingement region. At high nozzle-plate distance, the essence is to capture the evolution and breakup of the flow unsteadiness in the shear layer of the jet, so that accurate mean and fluctuating velocity profiles are obtained in the impingement region. Although the models have a quite different theoretical justification and generate a quite different eddy viscosity in some flow regions, their overall results are very comparable. The reason is that in zones that are crucial for the results, the models behave similarly.     

Modelling the non-linear interaction of wind and tide: Its influence on current profiles

A single-point model in the vertical is used to examine the coupling between tidal currents and wind-driven flows in shallow near-coastal regions. Calculations using both a linear slip and a no-slip condition at the sea bed clearly show that coupling between tidal and wind-driven currents cannot occur in a linear model with a time-independent eddy viscosity. However with a physically more realistic time-varying viscosity related to the flow field, coupling does occur, the magnitude of this non-linear interaction depending upon the change in eddy viscosity over a tidal cycle and the intensity of shear in the vertical. A point model in the vertical with flow induced by an oscillatory pressure gradient and an additional constant wind stress is used to examine the influence of viscosity parametrization and water depth upon this coupling. The solution in the vertical is accomplished using both a functional approach and a finite difference method. Some conclusions as to the relative merits of these approaches, particularly the use of a transformed grid in the case of high-shear surface and bed boundary layers, are made in the paper.     

A dynamic subgrid-scale tensorial Eddy viscosity model

In the Navier-Stokes equations the removal of the turbulent fluctuating velocities with a frequency above a certain fixed threshold, employed in the Large Eddy Simulation (LES), causes the appearance of a turbulent stress tensor that requires a number of closure assumptions. In this paper insufficiencies are demonstrated for those closure models which are based on a scalar eddy viscosity coefficient. A new model, based on a tensorial eddy viscosity, is therefore proposed; it employs the Germano identity [1] and allows dynamical evaluation of the single required input coefficient. The tensorial expression for the eddy viscosity is deduced by removing the widely used scalar assumption of the high-frequency viscous dissipation and replacing it by its tensorial counterpart arising in the balance of the Reynolds stress tensor. The numerical simulations performed for a lid driven cavity flow show that the proposed model allows to overcome the drawbacks encountered by the scalar eddy viscosity models. Received November 25, 1997     

Numerical problems in simulating tidal flows with a frictional-velocity-dependent eddy viscosity and the influence of stratification

The paper deals with the accurate determination of tidal current profiles in both homogeneous and stratified regions when a no-slip condition is used at the seabed with a flow-dependent eddy viscosity related to the depth-mean current or the bed frictional velocity. Calculations show that it is essential to accurately resolve the high-shear region which occurs at the bed and across the pycnocline/thermocline in the case of stratified flow. A computationally accurate and economic method of resolving these regions is demonstrated using the Galerkin method with a set of basis functions designed to accurately reproduce the high-shear layers which occur in these regions. With a flow-independent eddy viscosity a stability analysis can be readily performed and an unconditionally stable algorithm developed. However, with a flow-dependent viscosity, in particular a viscosity computed from the frictional velocity, a non-linear numerical instability can occur. A method of maintaining numerical stability in this case is also described. The importance of near-bed resolution to the computed value of the frictional velocity is demonstrated and its influence on the total tidal velocity profile is illustrated by a number of idealized calculations using various eddy viscosity formulations. The influence of stratification on the computed tidal profiles is shown in the latter part of the paper.     

A dynamic framework for functional parameterisations of the eddy viscosity coefficient in two-dimensional turbulence

This paper puts forth a dynamic framework for investigating the subgrid scale physics of decaying two-dimensional turbulence utilising a modular approach with eddy viscosities in various functional forms. The derivation of the low-pass spatially filtered implementation of the Navier–Stokes equations is given by using the vorticity-streamfunction formulation. Two different implementations of the viscosity kernels based on the representation of the eddy viscosity terms are proposed and tested by solving a canonical two-dimensional decaying turbulence problem. It is seen that the implementation of the eddy viscosity formulation plays a distinct role in the dissipative behaviour of the different viscosity kernels. Among eddy viscosity kernels tested, we found that the Leith eddy viscosity formulation yields superior results with higher correlation coefficients.     

An eddy viscosity model with near‐wall modifications

An extended version of the isotropic k–ε model is proposed that accounts for the distinct effects of low‐Reynolds number (LRN) and wall proximity. It incorporates a near‐wall correction term to amplify the level of dissipation in nonequilibrium flow regions, thus reducing the kinetic energy and length scale magnitudes to improve prediction of adverse pressure gradient flows, involving flow separation and reattachment. The eddy viscosity formulation maintains the positivity of normal Reynolds stresses and the Schwarz' inequality for turbulent shear stresses. The model coefficients/functions preserve the anisotropic characteristics of turbulence. The model is validated against a few flow cases, yielding predictions in good agreement with the direct numerical simulation (DNS) and experimental data. Comparisons indicate that the present model is a significant improvement over the standard eddy viscosity formulation. Copyright © 2005 John Wiley & Sons, Ltd.     

湍流边界层复涡黏模式的实验研究

在开口式循环水槽底部湍流边界层外区引入周期性扰动。利用X型热膜探针在扰动下游进行测量。用实验的方法研究了周期性大尺度结构下壁湍流涡黏模式中涡黏系数的形式，结果发现和周期性扰动对应的变形率及与之对应的雷诺应力间存在着相位差。这是目前许多最终导致涡黏系数的湍流模式理论都没有考虑到的一个重要因素。     

Godunov‐type adaptive grid model of wave–current interaction at cuspate beaches

This paper presents a second‐order accurate Godunov‐type numerical scheme for depth‐ and period‐averaged wave–current interaction. A flux Jacobian is derived for the wave conservation equations and its eigensystem determined, enabling Roe's approximate Riemann solver to be used to evaluate convective fluxes. Dynamically adaptive quadtree grids are used to focus on local hydrodynamic features, where sharp gradients occur in the flow variables. Adaptation criteria based on depth‐averaged vorticity, wave‐height gradient, wave steepness and the magnitude of velocity gradients are found to produce accurate solutions for nearshore circulation at a half‐sinusoidal beach. However, the simultaneous combination of two or more separate criteria produces numerical instability and interference unless all criteria are satisfied for mesh depletion. Simulations of wave–current interaction at a multi‐cusped beach match laboratory data from the United Kingdom Coastal Research Facility (UKCRF). A parameter study demonstrates the sensitivity of nearshore flow patterns to changes in relative cusp height, angle of wave incidence, bed roughness, offshore wave height and assumed turbulent eddy viscosity. Only a small deviation from normal wave incidence is required to initiate a meandering longshore current. Nearshore circulation patterns are highly dependent on the offshore wave height. Reduction of the assumed eddy viscosity parameter causes the primary circulation cells for normally incident waves to increase in strength whilst producing rip‐like currents cutting diagonally across the surf zone. Copyright © 2004 John Wiley & Sons, Ltd.     

Compound wall treatment with low‐Re turbulence model

A generalized treatment for the wall boundary conditions relating to turbulent flows is developed that blends the integration to a solid wall with wall functions. The blending function ensures a smooth transition between the viscous and turbulent regions. An improved low Reynolds number k?ε model is coupled with the proposed compound wall treatment to determine the turbulence field. The eddy viscosity formulation maintains the positivity of normal Reynolds stresses and Schwarz' inequality for turbulent shear stresses. The model coefficients/functions preserve the anisotropic characteristics of turbulence. Computations with fine and coarse meshes of a few flow cases yield appreciably good agreement with the direct numerical simulation and experimental data. The method is recommended for computing the complex flows where computational grids cannot satisfy a priori the prerequisites of viscous/turbulence regions. Copyright © 2011 John Wiley & Sons, Ltd.