A stabilized Nitsche‐type extended embedding mesh approach for 3D low‐ and high‐Reynolds‐number flows

This paper presents a stabilized extended finite element method (XFEM) based fluid formulation to embed arbitrary fluid patches into a fixed background fluid mesh. The new approach is highly beneficial when it comes to computational grid generation for complex domains, as it allows locally increased resolutions independent from size and structure of the background mesh. Motivating applications for such a domain decomposition technique are complex fluid‐structure interaction problems, where an additional boundary layer mesh is used to accurately capture the flow around the structure. The objective of this work is to provide an accurate and robust XFEM‐based coupling for low‐ as well as high‐Reynolds‐number flows. Our formulation is built from the following essential ingredients: Coupling conditions on the embedded interface are imposed weakly using Nitsche's method supported by extra terms to guarantee mass conservation and to control the convective mass transport across the interface for transient viscous‐dominated and convection‐dominated flows. Residual‐based fluid stabilizations in the interior of the fluid subdomains and accompanying face‐oriented fluid and ghost‐penalty stabilizations in the interface zone stabilize the formulation in the entire fluid domain. A detailed numerical study of our stabilized embedded fluid formulation, including an investigation of variants of Nitsche's method for viscous flows, shows optimal error convergence for viscous‐dominated and convection‐dominated flow problems independent of the interface position. Challenging two‐dimensional and three‐dimensional numerical examples highlight the robustness of our approach in all flow regimes: benchmark computations for laminar flow around a cylinder, a turbulent driven cavity flow at Re = 10000 and the flow interacting with a three‐dimensional flexible wall. Copyright © 2016 John Wiley & Sons, Ltd.     

Low Reynolds number k–ε model for near‐wall flow

A wall‐distance free k–ε turbulence model is developed that accounts for the near‐wall and low Reynolds number effects emanating from the physical requirements. The model coefficients/functions depend non‐linearly on both the strain rate and vorticity invariants. Included diffusion terms and modified C_ε(1,2) coefficients amplify the level of dissipation in non‐equilibrium 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 model is validated against a few flow cases, yielding predictions in good agreement with the direct numerical simulation (DNS) and experimental data. Copyright © 2005 John Wiley & Sons, Ltd.     

A γ‐model BGK scheme for compressible multifluids

We present a γ‐model BGK scheme for the numerical simulation of compressible multifluids. The scheme is based on the incorporation of a conservative γ‐model scheme given in (J. Comput. Phys. 1996; 125 :150–160) into the gas kinetic BGK scheme (J. Comput. Phys. 1993; 109 :53–66, J. Comput. Phys. 1994; 114 :9–17), and is simple to implement. Several numerical examples presented in this paper validate the scheme in the application of compressible multimaterial flows. Copyright © 2004 John Wiley & Sons, Ltd.     

Stability and consistency of nonhydrostatic free‐surface models using the semi‐implicit θ‐method

The θ‐method is a popular semi‐implicit finite‐difference method for simulating free‐surface flows. Problem stiffness, arising because of the presence of both fast and slow timescale processes, is easily handled by the θ‐method. In most ocean, coastal, and estuary modeling applications, stiffness is caused by fast surface gravity wave timescales imposed on slower timescales of baroclinic variability. The method is well known to be unconditionally stable for shallow water (hydrostatic) models when , where θ is the implicitness parameter. In this paper, we demonstrate that the method is also unconditionally stable for nonhydrostatic models, when for both pressure projection and pressure correction methods. However, the methods result in artificial damping of the barotropic mode. In addition to investigating stability, we also estimate the form of artificial damping induced by both the free surface and nonhydrostatic pressure solution methods. Finally, this analysis may be used to estimate the damping or growth associated with a particular wavenumber and the overall order of accuracy of the discretization. Copyright © 2012 John Wiley & Sons, Ltd.     

Implicit meanflow–multigrid algorithms for Reynolds stress model computation of 3‐D anisotropy‐driven and compressible flows

The present paper investigates the multigrid (MG) acceleration of compressible Reynolds‐averaged Navier–Stokes computations using Reynolds‐stress model 7‐equation turbulence closures, as well as lower‐level 2‐equation models. The basic single‐grid SG algorithm combines upwind‐biased discretization with a subiterative local‐dual‐time‐stepping time‐integration procedure. MG acceleration, using characteristic MG restriction and prolongation operators, is applied on meanflow variables only (MF–MG), turbulence variables being simply injected onto coarser grids. A previously developed non‐time‐consistent (for steady flows) full‐approximation‐multigrid (s–MG) is assessed for 3‐D anisotropy‐driven and/or separated flows, which are dominated by the convergence of turbulence variables. Even for these difficult test cases CPU‐speed‐ups r_CPUSUP∈[3, 5] are obtained. Alternative, potentially time‐consistent approaches (unsteady u–MG), where MG acceleration is applied at each subiteration, are also examined, using different subiterative strategies, MG cycles, and turbulence models. For 2‐D shock wave/turbulent boundary layer interaction, the fastest s–MG approach, with a V(2, 0) sawtooth cycle, systematically yields CPU‐speed‐ups of 5±½, quasi‐independent of the particular turbulence closure used. Copyright © 2008 John Wiley & Sons, Ltd.     

A non‐linear k–ε model with realizability for prediction of flows around bluff bodies

The incompressible flow around bluff bodies (a square cylinder and a cube) is investigated numerically using turbulence models. A non‐linear k–ε model, which can take into account the anisotropy of turbulence with less CPU time and computer memory then RSM or LES, is adopted as a turbulence model. In tuning of the model coefficients of the non‐linear terms are adjusted through the examination of previous experimental studies in simple shear flows. For the tuning of the coefficient in the eddy viscosity (=C_μ), the realizability constraints are derived in three types of basic 2D flow patterns, namely, a simple shear flow, flow around a saddle and a focal point. C_μ is then determined as a function of the strain and rotation parameters to satisfy the realizability. The turbulence model is first applied to a 2D flow around a square cylinder and the model performance for unsteady flows is examined focussing on the period and the amplitude of the flow oscillation induced by Karman vortex shedding. The applicability of the model to 3D flows is examined through the computation of the flow around a surface‐mounted cubic obstacle. The numerical results show that the present model performs satisfactorily to reproduce complex turbulent flows around bluff bodies. Copyright © 2003 John Wiley & Sons, Ltd.