Investigation of the sensitivity of turbulent closures and coupling of hybrid RANS‐LES models for predicting flow fields with separation and reattachment

Hybrid models have found widespread applications for simulation of wall‐bounded flows at high Reynolds numbers. Typically, these models employ Reynolds‐averaged Navier–Stokes (RANS) and large eddy simulation (LES) in the near‐body and off‐body regions, respectively. A number of coupling strategies between the RANS and LES regions have been proposed, tested, and applied in the literature with varying degree of success. Linear eddy‐viscosity models (LEVM) are often used for the closure of turbulent stress tensor in RANS and LES regions. LEVM incorrectly predicts the anisotropy of Reynolds normal stress at the RANS‐LES interface region. To overcome this issue, use of non‐linear eddy‐viscosity models (NLEVM) have started receiving attention. In this study, a generic non‐linear blended modeling framework for performing hybrid simulations is proposed. Flow over the periodic hills is used as the test case for model evaluation. This case is chosen due to complex flow physics with simplified geometry. Analysis of the simulations suggests that the non‐linear hybrid models show a better performance than linear hybrid models. It is also observed that the non‐linear closures are less sensitive to the RANS‐LES coupling and grid resolution. Copyright © 2016 John Wiley & Sons, Ltd.     

Improvement of moving particle semi‐implicit method for simulation of progressive water waves

Precise simulation of the propagation of surface water waves, especially when involving breaking wave, takes a significant place in computational fluid dynamics. Because of the strong nonlinear properties, the treatment of large surface deformation of free surface flow has always been a challenging work in the development of numerical models. In this paper, the moving particle semi‐implicit (MPS) method, an entirely Lagrangian method, is modified to simulate wave motion in a 2‐D numerical wave flume preferably. In terms of consecutive pressure distribution, a new and simple free surface detection criterion is proposed to enhance the free surface recognition in the MPS method. In addition, a revised gradient model is deduced to diminish the effect of nonuniform particle distribution and then to reduce the numerical wave attenuation occurring in the original MPS model. The applicability and stability of the improved MPS method are firstly demonstrated by the calculation of hydrostatic problem. It is revealed that these modifications are effective to suppress the pressure oscillation, weaken the local particle clustering, and boost the stability of numerical algorithm. It is then applied to investigate the propagation of progressive waves on a flat bed and the wave breaking on a mild slope. Comparisons with the analytical solutions and experimental results indicate that the improved MPS model can give better results about the profiles and heights of surface waves in contrast with the previous MPS models. Copyright © 2017 John Wiley & Sons, Ltd.     

波浪破碎的一种混合湍流模拟模式

近岸波浪的变形与破碎,一方面影响水体标识码运输,另一方面对消波护岸具有指导意义.本文提出一种三维混合湍流模拟模式方法,将求解区域分为造波区、波浪传播区和消波区. 造波区采用层流模式,通过基于 Fluent 的二次开发的 UDF 方法在边界进行速度造波.这种方法在给定入口速度的条件标识码据已知波高进行 精准插值,从而控制水的体积分数. 在波浪传播区域,采用大涡标识码行模拟研究,在消波区,采用 RANS 模型并利用多孔介质消波法进行消波. 模标识码 VOF 方法捕捉波浪破碎过程中的自由面变化.本文对波高为 5.5cm 的规则波 (M1)、波高为 13.5cm 的 规则波 (M3)、有效标识码 7.75cm 的 TMA 谱单向不规则波 (U1) 和有效波高为 19cm 的 TMA 谱多向不规则波 (B5) 展开标识码研究,并与前人的相关标识码果作比较,各条件下模拟结果与实验结果吻合.模拟结果说明本文提出的模型能够准确模拟出波浪传播过程中的折射和标识码象,并且能够捕捉的波浪破碎过程中的自由面变化,为三维波浪的传播与破碎的数值模拟提供一种模拟方法.      

超声速燃烧火焰稳定凹腔质量交换特性的数值研究

对超声速冷流条件下用于超燃冲压发动机的凹腔火焰稳定器的质量交换特性进行研 究. 采用混合RANS/LES方法对非定常流场进行数值模拟,考虑了凹腔的长深比和后缘角度两 个关键参数. 计算得到了凹腔剪切层拟序结构的演化过程. 对凹腔压力振荡历程进行幅频分 析,所得到的频率和理论分析结果及一些计算结果十分吻合. 结果表明,凹腔的长深比和后 缘倾角对凹腔质量交换过程都有较大的影响. 凹腔驻留时间随着长深比的增加而增加. 随着 后缘倾角减小,长深比较小的凹腔的驻留时间增加,长深比较大的凹腔则完全相反. 该计算 结果为凹腔设计提供了新认识.     

Hybrid finite‐volume finite‐difference scheme for the solution of Boussinesq equations

A hybrid scheme composed of finite‐volume and finite‐difference methods is introduced for the solution of the Boussinesq equations. While the finite‐volume method with a Riemann solver is applied to the conservative part of the equations, the higher‐order Boussinesq terms are discretized using the finite‐difference scheme. Fourth‐order accuracy in space for the finite‐volume solution is achieved using the MUSCL‐TVD scheme. Within this, four limiters have been tested, of which van‐Leer limiter is found to be the most suitable. The Adams–Basforth third‐order predictor and Adams–Moulton fourth‐order corrector methods are used to obtain fourth‐order accuracy in time. A recently introduced surface gradient technique is employed for the treatment of the bottom slope. A new model ‘HYWAVE’, based on this hybrid solution, has been applied to a number of wave propagation examples, most of which are taken from previous studies. Examples include sinusoidal waves and bi‐chromatic wave propagation in deep water, sinusoidal wave propagation in shallow water and sinusoidal wave propagation from deep to shallow water demonstrating the linear shoaling properties of the model. Finally, sinusoidal wave propagation over a bar is simulated. The results are in good agreement with the theoretical expectations and published experimental results. Copyright © 2005 John Wiley & Sons, Ltd.     

Derivation,implementation, and initial testing of a compressible wall‐layer model

Traditional wall functions have been used successfully for decades to decrease the computational cost for obtaining solutions to incompressible flows with equilibrium turbulence. However, these traditional, analytic wall functions are poorly suited for more complex flows. The present work describes an alternative approach named the ‘diffusion model’. The diffusion model is a subgrid model developed by Blottner and Bond that solves a system of ODEs in the near‐wall region instead of assuming an analytic profile. The diffusion model has previously been shown to reproduce profiles of various turbulence models through the log layer with fixed outer boundary conditions. This paper documents the implementation and the testing of the diffusion model fully coupled with a 3‐D Reynolds‐averaged Navier–Stokes (RANS) algorithm, using the Spalart–Allmaras model. The results indicate that the coupled algorithm is valid when the interface between diffusion model and the 3‐D RANS algorithm is below the outer edge of the log layer. Although the current work focuses on using steady 3‐D RANS outside of the log layer, modeling assumptions introduced in the current work could be used to derive a diffusion model for coupling with a large eddy simulation region. Published in 2010 by John Wiley & Sons, Ltd.     

SPH computation of plunging waves using a 2‐D sub‐particle scale (SPS) turbulence model

The paper presents a 2‐D large eddy simulation (LES) modelling approach to investigate the properties of the plunging waves. The numerical model is based on the smoothed particle hydrodynamics (SPH) method. SPH is a mesh‐free Lagrangian particle approach which is capable of tracking the free surfaces of large deformation in an easy and accurate way. The Smagorinsky model is used as the turbulence model due to its simplicity and effectiveness. The proposed 2‐D SPH–LES model is applied to a cnoidal wave breaking and plunging over a mild slope. The computations are in good agreement with the documented data. Especially the computed turbulence quantities under the breaking waves agree better with the experiments as compared with the numerical results obtained by using the k–ε model. The sensitivity analyses of the SPH–LES computations indicate that both the turbulence model and the spatial resolution play an important role in the model predictions and the contributions from the sub‐particle scale (SPS) turbulence decrease with the particle size refinement. Copyright © 2006 John Wiley & Sons, Ltd.     

High‐order detached eddy simulation,zonal LES and URANS of cavity and labyrinth seal flows

Hybrid numerical large eddy simulation (NLES), detached eddy simulation (DES) and URANS methods are assessed on a cavity and a labyrinth seal geometry. A high sixth‐order discretization scheme is used and is validated using the test case of a two‐dimensional vortex. The hybrid approach adopts a new blending function. For the URANS simulations, the flow within the cavity remains steady, and the results show significant variation between models. Surprisingly, low levels of resolved turbulence are observed in the cavity for the DES simulation, and the cavity shear layer remains two dimensional. The hybrid RANS–NLES approach does not suffer from this trait. For the labyrinth seal, both the URANS and DES approaches give low levels of resolved turbulence. The zonal Hamilton–Jacobi approach on the other had given significantly more resolved content. Both DES and hybrid RANS–NLES give good agreement with the experimentally measured velocity profiles. Again, there is significant variation between the URANS models, and swirl velocities are overpredicted. Copyright © 2013 John Wiley & Sons, Ltd.     

A multi‐energy‐level lattice Boltzmann model for two‐dimensional wave equation

A lattice Boltzmann model for two‐dimensional wave equation is presented. In this model, we used higher‐order moment method, multi‐scale technique and Chapman–Enskog expansion, and multi‐energy‐level to obtain wave equation and energy conservation equation. As numerical examples, the interference and diffraction of wave are simulated. The numerical results show this model can be used to simulate two‐dimensional wave propagation. Copyright © 2009 John Wiley & Sons, Ltd.     

海堤越浪的数值模拟

基于RANS方程和两方程湍流模型,采用有限体积法,将人射波波场作为人工的分布源项加人动量方程,提出了适用于VOF方法的源造波一消波技术。通过对行波及驻波的计算,分别考察了数值波浪水槽前端及末端消波段的有效性。在本文建立的数值波浪水槽内对规则波在海堤上爬高和越浪过程进行了数值模拟,并将计算结果与现有实验结果进行了比较。验证计算结果表明,数值模拟结果较好地复演了海堤越浪过程。为了研究模型尺度对越浪量的影响,文中设计了一组满足重力相似但具有不同几何比尺的数值实验模型。系列数值实验结果表明,若按重力相似换算越浪量,计算结果与实验预报值间的偏差随模型比尺的增大和堤前波浪破碎强度的增强而增大,建议在进行越浪物理模型实验时需进一步考虑模型比尺对原型预报值的影响。     

Computation of turbulent free‐surface flows around modern ships

This paper presents the calculated results for three classes of typical modern ships in modelling of ship‐generated waves. Simulations of turbulent free‐surface flows around ships are performed in a numerical water tank, based on the FINFLO‐RANS SHIP solver developed at Helsinki University of Technology. The Reynolds‐averaged Navier–Stokes (RANS) equations with the artificial compressibility and the non‐linear free‐surface boundary conditions are discretized by means of a cell‐centred finite‐volume scheme. The convergence performance is improved with the multigrid method. A free surface is tracked using a moving mesh technology, in which the non‐linear free‐surface boundary conditions are given on the actual location of the free surface. Test cases recommended are a container ship, a US Navy combatant and a tanker. The calculated results are compared with the experimental data available in the literature in terms of the wave profiles, wave pattern, and turbulent flow fields for two turbulence models, Chien's low Reynolds number k–εmodel and Baldwin–Lomax's model. Furthermore, the convergence performance, the grid refinement study and the effect of turbulence models on the waves have been investigated. Additionally, comparison of two types of the dynamic free‐surface boundary conditions is made. Copyright © 2003 John Wiley& Sons, Ltd.     

A lattice Boltzmann model for two‐dimensional sound wave in the small perturbation compressible flows

A multi‐entropy‐level lattice Boltzmann model for two‐dimensional sound wave equation in the small perturbation flows is presented. In this model, we used higher‐order moment method, multi‐scale technique and the Chapman–Enskog expansion, and multi‐entropy‐level to obtain sound wave equation with isentropic equation. As numerical examples, the Doppler effects in the sound wave propagation, the sound scattering from circular cylinder are simulated. The numerical results show that this model can be used to simulate sound wave propagation. Copyright © 2010 John Wiley & Sons, Ltd.     

Numerical simulations of wave interactions with vertical wave barriers using the SPH method

This paper focuses on the fluid boundary separation problem of the conventional dynamic solid boundary treatment (DSBT) and proposes a modified DSBT (MDSBT). Classic 2D free dam break flows and 3D dam break flows against a rectangular box are used to assess the performance of this MDSBT in free surface flow and violent fluid–structure interaction, respectively. Another test, water column oscillations in a U‐tube, is specially designed to reveal the applicability of dealing with two types of particular boundaries: the wet–dry solid boundary and the large‐curvature solid boundary. A comparison between the numerical results and the experimental data shows that the MDSBT is capable of eliminating the fluid boundary separation, improving the accuracy of the solid boundary pressure calculations and preventing the unphysical penetration of fluid particles. Using a 2D SPH numerical wave tank with MDSBT, the interactions between regular waves and a simplified vertical wave barrier are simulated. The numerical results reveal that the maximum horizontal force occurs at the endpoint of the vertical board, and with the enlargement of the relative submerged board length, the maximum moment grows linearly; furthermore, the relative average mass transportation under the breakwater initially increases to 11.14 per wave strike but is later reduced. The numerical simulation of a full‐scale 3D wave barrier with two vertical boards shows that the wave and structure interactions in the practical project are far more complicated than in the simplified 2D models. The SPH model using the MDSBT is capable of providing a reference for engineering designs. Copyright © 2014 John Wiley & Sons, Ltd.     

Numerical prediction of turbulent flows using Reynolds‐averaged Navier–Stokes and large‐eddy simulation with uncertain inflow conditions

Numerous comparisons between Reynolds‐averaged Navier–Stokes (RANS) and large‐eddy simulation (LES) modeling have already been performed for a large variety of turbulent flows in the context of fully deterministic flows, that is, with fixed flow and model parameters. More recently, RANS and LES have been separately assessed in conjunction with stochastic flow and/or model parameters. The present paper performs a comparison of the RANS k ? ε model and the LES dynamic Smagorinsky model for turbulent flow in a pipe geometry subject to uncertain inflow conditions. The influence of the experimental uncertainties on the computed flow is analyzed using a non‐intrusive polynomial chaos approach for two flow configurations (with or without swirl). Measured quantities including an estimation of the measurement error are then compared with the statistical representation (mean value and variance) of their RANS and LES numerical approximations in order to check whether experiment/simulation discrepancies can be explained within the uncertainty inherent to the studied configuration. The statistics of the RANS prediction are found in poor agreement with experimental results when the flow is characterized by a strong swirl, whereas the computationally more expensive LES prediction remains statistically well inside the measurement intervals for the key flow quantities.Copyright © 2012 John Wiley & Sons, Ltd.     

Simulation of dam‐ and dyke‐break hydrodynamics on dynamically adaptive quadtree grids

Flooding due to the failure of a dam or dyke has potentially disastrous consequences. This paper presents a Godunov‐type finite volume solver of the shallow water equations based on dynamically adaptive quadtree grids. The Harten, Lax and van Leer approximate Riemann solver with the Contact wave restored (HLLC) scheme is used to evaluate interface fluxes in both wet‐ and dry‐bed applications. The numerical model is validated against results from alternative numerical models for idealized circular and rectangular dam breaks. Close agreement is achieved with experimental measurements from the CADAM dam break test and data from a laboratory dyke break undertaken at Delft University of Technology. Copyright © 2004 John Wiley Sons, Ltd.     

Validation of a 2-D semi-coupled numerical model for fluid–structure–seabed interaction

A 2-D semi-coupled model PORO-WSSI 2D (also be referred as FSSI-CAS 2D) for the Fluid-Structure-Seabed Interaction (FSSI) has been developed by employing RANS equations for wave motion in fluid domain, VARANS equations for porous flow in porous structures; and taking the dynamic Biot's equations (known as "u − p" approximation) for soil as the governing equations. The finite difference two-step projection method and the forward time difference method are adopted to solve the RANS, VARANS equations; and the finite element method is adopted to solve the "u − p" approximation. A data exchange port is developed to couple the RANS, VARANS equations and the dynamic Biot's equations together. The analytical solution proposed by Hsu and Jeng (1994) and some experiments conducted in wave flume or geotechnical centrifuge in which various waves involved are used to validate the developed semi-coupled numerical model. The sandy bed involved in these experiments is poro-elastic or poro-elastoplastic. The inclusion of the interaction between fluid, marine structures and poro-elastoplastic seabed foundation is a special point and highlight in this paper, which is essentially different with other previous coupled models The excellent agreement between the numerical results and the experiment data indicates that the developed coupled model is highly reliablefor the FSSI problem.     

A numerical model for aerated‐water wave breaking

This work presents a numerical model designed for the simulation of water‐wave impacts on a structure when aeration of the liquid phase is considered. The model is based on a multifluid Navier–Stokes approach in which all fluids are assumed compressible. The numerical method is based on a finite volume algorithm in space and a second order Runge–Kutta method in time. A validation of this model is performed. It shows a good accuracy for acoustic and shock wave propagation in a bubbly liquid and for wave breaking. Copyright © 2011 John Wiley & Sons, Ltd.     

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.     

A semi‐implicit finite element model for non‐hydrostatic (dispersive) surface waves

The objective of this research is to develop a model that will adequately simulate the dynamics of tsunami propagating across the continental shelf. In practical terms, a large spatial domain with high resolution is required so that source areas and runup areas are adequately resolved. Hence efficiency of the model is a major issue. The three‐dimensional Reynolds averaged Navier–Stokes equations are depth‐averaged to yield a set of equations that are similar to the shallow water equations but retain the non‐hydrostatic pressure terms. This approach differs from the development of the Boussinesq equations where pressure is eliminated in favour of high‐order velocity and geometry terms. The model gives good results for several test problems including an oscillating basin, propagation of a solitary wave, and a wave transformation over a bar. The hydrostatic and non‐hydrostatic versions of the model are compared for a large‐scale problem where a fault rupture generates a tsunami on the New Zealand continental shelf. The model efficiency is also very good and execution times are about a factor of 1.8 to 5 slower than the standard shallow water model, depending on problem size. Moreover, there are at least two methods to increase model accuracy when warranted: choosing a more optimal vertical interpolation function, and dividing the problem into layers. Copyright © 2005 John Wiley & Sons, Ltd.