Shallow water equations: Magic and limitations

D Liang from Cambridge University explains the shallow water equations and their applications to the dam-break and other steep-fronted flow modeling. They assume that the horizontal scale of the flow is much greater than the vertical scale, which means the flow is restricted within a thin layer, thus the vertical momentum is insignificant and the pressure distribution is hydrostatic. The left hand sides of the two momentum equations represent the acceleration of the fluid particle in the horizontal plane. If the fluid acceleration is ignored, then the two momentum equations are simplified into the so-called diffusion wave equations. In contrast to the SWEs approach, it is much less convenient to model floods with the Navier-Stokes equations. In conventional computational fluid dynamics (CFD), cumbersome treatments are needed to accurately capture the shape of the free surface. The SWEs are derived using the assumptions of small vertical velocity component, smooth water surface, gradual variation and hydrostatic pressure distribution.

Unsteady transition in an axial-flow turbine. Part 2. Cascade measurements and modeling

Part 1 of this paper reanalyzed previously published measurements from the rotor of a low-speed, single-stage, axial-flow turbine, which highlighted the unsteady nature of the suction surface transition process. Part 2 investigates the significance of the wake jet and the unsteady frequency parameter. Supporting experiments carried out in a linear cascade with varying inlet turbulence are described, together with a simple unsteady transition model explaining the features of seen in the turbine.

Intermittency transport modeling of separated flow transition

An intermittency transport model is proposed for modeling separated-flow transition. The model is based on earlier work on prediction of attached flow bypass transition and is applied for the first time to model transition in a separation bubble at various degrees of free-stream turbulence. The model has been developed so that it takes into account the entrainment of the surrounding fluid. Experimental investigations suggest that it is this phenomena which ultimately determines the extent of the separation bubble. Transition onset is determined via a boundary layer correlation based on momentum thickness at the point of separation. The intermittent flow characteristic of the transition process is modeled via an intermittency transport equation. This accounts for both normal and streamwise variation of intermittency and hence models the entrainment of surrounding flow in a more accurate manner than alternative prescribed intermittency models. The model has been validated against the well established T3L semicircular leading edge flat plate test case for three different degrees of free-stream turbulence characteristic of turbomachinery blade applications.

A FEniCS-based programming framework for modeling turbulent flow by the Reynolds-averaged Navier-Stokes equations

Finding an appropriate turbulence model for a given flow case usually calls for extensive experimentation with both models and numerical solution methods. This work presents the design and implementation of a flexible, programmable software framework for assisting with numerical experiments in computational turbulence. The framework targets Reynolds-averaged Navier-Stokes models, discretized by finite element methods. The novel implementation makes use of Python and the FEniCS package, the combination of which leads to compact and reusable code, where model- and solver-specific code resemble closely the mathematical formulation of equations and algorithms. The presented ideas and programming techniques are also applicable to other fields that involve systems of nonlinear partial differential equations. We demonstrate the framework in two applications and investigate the impact of various linearizations on the convergence properties of nonlinear solvers for a Reynolds-averaged Navier-Stokes model. © 2011 Elsevier Ltd.

CFD modeling of single-phase flow in a packed bed with MRI validation

Computational fluid dynamics (CFD) simulations are becoming increasingly widespread with the advent of more powerful computers and more sophisticated software. The aim of these developments is to facilitate more accurate reactor design and optimization methods compared to traditional lumped-parameter models. However, in order for CFD to be a trusted method, it must be validated using experimental data acquired at sufficiently high spatial resolution. This article validates an in-house CFD code by comparison with flow-field data obtained using magnetic resonance imaging (MRI) for a packed bed with a particle-to-column diameter ratio of 2. Flows characterized by inlet Reynolds numbers, based on particle diameter, of 27, 55, 111, and 216 are considered. The code used employs preconditioning to directly solve for pressure in low-velocity flow regimes. Excellent agreement was found between the MRI and CFD data with relative error between the experimentally determined and numerically predicted flow-fields being in the range of 3-9%. © 2012 American Institute of Chemical Engineers (AIChE).

Turbulence modeling for flows around convex features giving rapid eddy distortion

Reynolds averaged Navier-Stokes model performances in the stagnation and wake regions for turbulent flows with relatively large Lagrangian length scales (generally larger than the scale of geometrical features) approaching small cylinders (both square and circular) is explored. The effective cylinder (or wire) diameter based Reynolds number, ReW ≤ 2.5 × 103. The following turbulence models are considered: a mixing-length; standard Spalart and Allmaras (SA) and streamline curvature (and rotation) corrected SA (SARC); Secundov's νt-92; Secundov et al.'s two equation νt-L; Wolfshtein's k-l model; the Explicit Algebraic Stress Model (EASM) of Abid et al.; the cubic model of Craft et al.; various linear k-ε models including those with wall distance based damping functions; Menter SST, k-ω and Spalding's LVEL model. The use of differential equation distance functions (Poisson and Hamilton-Jacobi equation based) for palliative turbulence modeling purposes is explored. The performance of SA with these distance functions is also considered in the sharp convex geometry region of an airfoil trailing edge. For the cylinder, with ReW ≈ 2.5 × 103 the mixing length and k-l models give strong turbulence production in the wake region. However, in agreement with eddy viscosity estimates, the LVEL and Secundov νt-92 models show relatively little cylinder influence on turbulence. On the other hand, two equation models (as does the one equation SA) suggest the cylinder gives a strong turbulence deficit in the wake region. Also, for SA, an order or magnitude cylinder diameter decrease from ReW = 2500 to 250 surprisingly strengthens the cylinder's disruptive influence. Importantly, results for ReW ≪ 250 are virtually identical to those for ReW = 250 i.e. no matter how small the cylinder/wire its influence does not, as it should, vanish. Similar tests for the Launder-Sharma k-ε, Menter SST and k-ω show, in accordance with physical reality, the cylinder's influence diminishing albeit slowly with size. Results suggest distance functions palliate the SA model's erroneous trait and improve its predictive performance in wire wake regions. Also, results suggest that, along the stagnation line, such functions improve the SA, mixing length, k-l and LVEL results. For the airfoil, with SA, the larger Poisson distance function increases the wake region turbulence levels by just under 5%. © 2007 Elsevier Inc. All rights reserved.