scholarly journals Use and Scalability of OpenFOAM for Wind Fields and Pollution Dispersion with Building- and Ground-Resolving Topography

Atmosphere ◽  
2021 ◽  
Vol 12 (9) ◽  
pp. 1124
Author(s):  
Daniel Elfverson ◽  
Christian Lejon
Keyword(s):  
Laser Scanning ◽  
Situational Awareness ◽  
Pollutant Dispersion ◽  
Passive Scalar ◽  
Navier Stokes ◽  
Computational Domain ◽  
Wind Fields ◽  
Complex Flow ◽  
Pollution Dispersion ◽  
Urban Geometry

Complex flow and pollutant dispersion simulations in real urban settings were investigated by using computational fluid dynamics (CFD) simulations with the SST k−ω Reynolds-averaged Navier–Stokes (RANS) equation with OpenFOAM. The model was validated with a wind-tunnel experiment using two surface-mounted cubes in tandem, and the flow features were reproduced with the correct qualitative behaviour. The real urban geometry of the Parade Square in Warsaw, Poland was represented with both laser-scanning data for the ground geometry and the CityGML standard to describe the buildings as an example. The Eulerian dispersion of a passive scalar and the flow behaviour could be resolved within minutes over a computational domain with a size of 958 × 758 m2 and a height of 300 m with over 2 M cells due to the good and strong parallel scalability in OpenFOAM. This implies that RANS modelling with parallel computing in OpenFOAM can potentially be used as a tool for situational awareness on a local urban scale; however, entire cities would be too large.

Steady and Unsteady Modeling for Heat Transfer Predictions of High Pressure Turbine Blade Internal Cooling

2012 ◽  
Author(s):  
Rémy Fransen ◽  
Nicolas Gourdain ◽  
Laurent Y. M. Gicquel
Keyword(s):  
Heat Transfer ◽  
Large Scale ◽  
Cooling System ◽  
Transfer Problem ◽  
Operating Conditions ◽  
Navier Stokes ◽  
Wall Region ◽  
Internal Cooling ◽  
Computational Domain ◽  
Complex Flow

This work focuses on numerical simulations of flows in blade internal cooling system. Large Eddy Simulation (LES) and Reynolds-Averaged Navier Stokes (RANS) approaches are compared in a typical blade cooling related problem. The case is a straight rib-roughened channel with high blockage ratio, computed and compared for both a periodic and full spatial domains. The configuration was measured at the Von Karman Institute (VKI) using Particle Image Velocimetry (PIV) in near gas turbine operating conditions. Results show that RANS models used fail to predict the full evolution of the flow within the channels where massive separation and large scale unsteady features are evidenced. In contrast LES succeeds in reproducing these complex flow motions and both mean and fluctuating components are clearly improved in the channels and in the near wall region. Periodic computations are gauged against the spatial computational domain and results on the heat transfer problem are addressed.

scholarly journals Simulation of Turbulent Wake at Mixing of Two Confined Horizontal Flows

2018 ◽  
Vol 2018 ◽  
pp. 1-12 ◽  
Author(s):  
Rok Krpan ◽  
Boštjan Končar
Keyword(s):  
Stokes Equations ◽  
Mixing Layer ◽  
Passive Scalar ◽  
Systematic Investigation ◽  
Navier Stokes ◽  
Computational Domain ◽  
Navier Stokes Equations ◽  
Turbulent Mixing Layer ◽  
Eddy Simulation ◽  
Benchmark Experiment

The development of a turbulent mixing layer at mixing of two horizontal water streams with slightly different densities is studied by the means of numerical simulation. The mixing of such flows can be modelled as the flow of two components, where the concentration of one component in the mixing region is described as a passive scalar. The velocity field remains common over the entire computational domain, where the density and viscosity difference due to the concentration mainly affects the turbulent fluctuations in the mixing region. The numerical simulations are performed with the open source code OpenFOAM using two different approaches for turbulence modelling, Reynolds Averaged Navier Stokes equations (RANS) and Large Eddy Simulation (LES). The simulation results are discussed and compared with the benchmark experiment obtained within the frame of OECD/NEA benchmark test. A good agreement with experimental results is obtained in the case of the single liquid experiment. A high discrepancy between the simulated and the experimental velocity fluctuations in the case of mixing of the flows with the slightly different densities and viscosities triggered a systematic investigation of the modelling approaches that helped us to find out and interpret the main reasons for the disagreement.

scholarly journals A Code for Simulating Heat Transfer in Turbulent Channel Flow

Mathematics ◽  
2021 ◽  
Vol 9 (7) ◽  
pp. 756
Author(s):  
Federico Lluesma-Rodríguez ◽  
Francisco Álcantara-Ávila ◽  
María Jezabel Pérez-Quiles ◽  
Sergio Hoyas
Keyword(s):  
Heat Transfer ◽  
Stokes Equations ◽  
Three Dimensional ◽  
Turbulent Channel Flow ◽  
Passive Scalar ◽  
Direct Numerical Simulations ◽  
Time Dependent ◽  
Navier Stokes ◽  
Thermal Flow ◽  
Navier Stokes Equations

One numerical method was designed to solve the time-dependent, three-dimensional, incompressible Navier–Stokes equations in turbulent thermal channel flows. Its originality lies in the use of several well-known methods to discretize the problem and its parallel nature. Vorticy-Laplacian of velocity formulation has been used, so pressure has been removed from the system. Heat is modeled as a passive scalar. Any other quantity modeled as passive scalar can be very easily studied, including several of them at the same time. These methods have been successfully used for extensive direct numerical simulations of passive thermal flow for several boundary conditions.

On LES Methods Applied to Seal Geometries

2012 ◽  
Author(s):  
James Tyacke ◽  
Richard Jefferson-Loveday ◽  
Paul Tucker
Keyword(s):  
Maximum Error ◽  
Labyrinth Seal ◽  
Navier Stokes ◽  
Final Solution ◽  
Complex Flow ◽  
Eddy Simulation ◽  
Wide Range ◽  
Large Eddy ◽  
Rans Models ◽  
Flow Through

Nine Large Eddy Simulation (LES) methods are used to simulate flow through two labyrinth seal geometries and are compared with a wide range of Reynolds-Averaged Navier-Stokes (RANS) solutions. These involve one-equation, two-equation and Reynolds Stress RANS models. Also applied are linear and nonlinear pure LES models, hybrid RANS-Numerical-LES (RANS-NLES) and Numerical-LES (NLES). RANS is found to have a maximum error and a scatter of 20%. A similar level of scatter is also found among the same turbulence model implemented in different codes. In a design context, this makes RANS unusable as a final solution. Results show that LES and RANS-NLES is capable of accurately predicting flow behaviour of two seals with a scatter of less than 5%. The complex flow physics gives rise to both laminar and turbulent zones making most LES models inappropriate. Nonetheless, this is found to have minimal tangible results impact. In accord with experimental observations, the ability of LES to find multiple solutions due to solution non-uniqueness is also observed.

ANALYSIS OF A HETEROGENEOUS DOMAIN DECOMPOSITION FOR COMPRESSIBLE VISCOUS FLOW

2001 ◽  
Vol 11 (04) ◽  
pp. 565-599 ◽  
Author(s):  
CRISTIAN A. COCLICI ◽  
WOLFGANG L. WENDLAND
Keyword(s):  
Domain Decomposition ◽  
Stokes Equations ◽  
Domain Decomposition Method ◽  
Near Field ◽  
Outer Region ◽  
Navier Stokes ◽  
Far Field ◽  
Computational Domain ◽  
Linearized Euler Equations ◽  
Naca0012 Airfoil

We analyze a nonoverlapping domain decomposition method for the treatment of two-dimensional compressible viscous flows around airfoils. Since at some distance to the given profile the inertial forces are strongly dominant, there the viscosity effects are neglected and the flow is assumed to be inviscid. Accordingly, we consider a decomposition of the original flow field into a bounded computational domain (near field) and a complementary outer region (far field). The compressible Navier–Stokes equations are used close to the profile and are coupled with the linearized Euler equations in the far field by appropriate transmission conditions, according to the physical properties and the mathematical type of the corresponding partial differential equations. We present some results of flow around the NACA0012 airfoil and develop an a posteriori analysis of the approximate solution, showing that conservation of mass, momentum and energy are asymptotically attained with the linear model in the far field.

Fluid Flow Structure in Arterial Bypass Anastomosis

2005 ◽  
Vol 127 (4) ◽  
pp. 611-618 ◽  
Author(s):  
C. M. Su ◽  
D. Lee ◽  
R. Tran-Son-Tay ◽  
W. Shyy
Keyword(s):  
Fluid Flow ◽  
Flow Structure ◽  
Bypass Graft ◽  
Flow Characteristics ◽  
Navier Stokes ◽  
Complex Flow ◽  
Arterial Bypass ◽  
Flow Recirculation ◽  
Area Reduction ◽  
Complex Flow Structure

The fluid flow through a stenosed artery and its bypass graft in an anastomosis can substantially influence the outcome of bypass surgery. To help improve our understanding of this and related issues, the steady Navier-Stokes flows are computed in an idealized arterial bypass system with partially occluded host artery. Both the residual flow issued from the stenosis—which is potentially important at an earlier stage after grafting—and the complex flow structure induced by the bypass graft are investigated. Seven geometric models, including symmetric and asymmetric stenoses in the host artery, and two major aspects of the bypass system, namely, the effects of area reduction and stenosis asymmetry, are considered. By analyzing the flow characteristics in these configurations, it is found that (1) substantial area reduction leads to flow recirculation in both upstream and downstream of the stenosis and in the host artery near the toe, while diminishes the recirculation zone in the bypass graft near the bifurcation junction, (2) the asymmetry and position of the stenosis can affect the location and size of these recirculation zones, and (3) the curvature of the bypass graft can modify the fluid flow structure in the entire bypass system.

Case Study of Nanjing Inner Ring on Air Pollution Dispersion from Roadway Tunnel Portals

2012 ◽  
Vol 610-613 ◽  
pp. 1895-1900 ◽  
Author(s):  
Shu Jiang Miao ◽  
Da Fang Fu
Keyword(s):  
Air Pollution ◽  
Pollutant Dispersion ◽  
Air Pollutant ◽  
Urban Traffic ◽  
Lagrangian Model ◽  
Pollution Dispersion ◽  
Piecewise Function ◽  
Air Pollutant Dispersion ◽  
Air Pollution Dispersion

The tunnel module of a rather simple Lagrangian model GRAL (Grazer Langrange model) has been chosen to study air pollutant dispersion around tunnel portals in Nanjing inner ring. Two points have been made to popularize GRAL3.5TM (the tunnel module of a Lagrangian model GRAL; the update was in May 2003) and assure it more suitable for the actual situations in Nanjing. One is to derive a piecewise function of the intermediate parameter ‘stiffness’. Another is to take Romberg NOx-NO2 scheme into account. After these 2 works on GRAL3.5TM, NO2 dispersion from portals of all the 6 tunnels in Nanjing inner ring has been simulated. The importance of limiting urban traffic volume to control air quality around tunnel portals and roadways has been emphasized.

Reducing Undesired Wave Reflection at Domain Boundaries in 3D Finite Volume-Based Flow Simulations via Forcing Zones

2020 ◽  
Vol 64 (01) ◽  
pp. 23-47
Author(s):  
Robinson Peric ◽  
Moustafa Abdel-Maksoud
Keyword(s):  
Finite Volume ◽  
Wave Reflection ◽  
Three Dimensional ◽  
Navier Stokes ◽  
Special Focus ◽  
Computational Domain ◽  
Wave Reflections ◽  
Flow Simulations ◽  
Large Eddy ◽  
Domain Boundaries

This article reviews different types of forcing zones (sponge layers, damping zones, relaxation zones, etc.) as used in finite volume-based flow simulations to reduce undesired wave reflections at domain boundaries, with special focus on the case of strongly reflecting bodies subjected to long-crested incidence waves. Limitations and possible sources of errors are discussed. A novel forcing-zone arrangement is presented and validated via three-dimensional (3D) flow simulations. Furthermore, a recently published theory for predicting the forcing-zone behavior was investigated with regard to its relevance for practical 3D hydrodynamics problems. It was found that the theory can be used to optimally tune the case-dependent parameters of the forcing zones before running the simulations. 1. Introduction Wave reflections at the boundaries of the computational domain can cause significant errors in flow simulations, and must therefore be reduced. In contrast to boundary element codes, where much progress in this respect has been made decades ago (see e.g., Clement 1996; Grilli &Horillo 1997), for finite volume-based flow solvers, there are many unresolved questions, especially:How to reliably reduce reflections and disturbances from the domain boundaries?How to predict the amount of undesired wave reflection before running the simulation? This work aims to provide further insight to these questions for flow simulations based on Navier-Stokes-type equations (Reynolds-averaged Navier-Stokes, Euler equations, Large Eddy Simulations, etc.), when using forcing zones to reduce undesired reflections. The term "forcing zones" is used here to describe approaches that gradually force the solution in the vicinity of the boundary towards some reference solution, as described in Section 2; some examples are absorbing layers, sponge layers, damping zones, relaxation zones, or the Euler overlay method (Mayer et al. 1998; Park et al. 1999; Chen et al. 2006; Choi &Yoon 2009; Jacobsen et al. 2012; Kimet al. 2012; Schmitt & Elsaesser 2015; Perić & Abdel-Maksoud 2016a; Vukčević et al. 2016).

Computation of Discrete-Hole Film Cooling by a Near-Wall Second-Moment Turbulence Closure

1999 ◽  
Author(s):  
Hamn-Ching Chen ◽  
Gengsheng Wei ◽  
Je-Chin Han
Keyword(s):  
Film Cooling ◽  
Moment Closure ◽  
Navier Stokes ◽  
Closure Model ◽  
Complex Flow ◽  
Second Moment ◽  
Finite Analytic Method ◽  
Second Moment Closure ◽  
Effectiveness Prediction ◽  
Near Wall

Abstract A multiblock Favre-Averaged Navier-Stokes (FANS) method has been developed in conjunction with a chimera domain decomposition technique for investigation of flat surface, discrete-hole film cooling performance. The finite-analytic method solves the FANS equations in conjunction with a near-wall second-order Reynolds stress (second-moment) closure model and a two-layer k-ε model. Comparisons of flow fields and turbulence quantities with experimental data clearly demonstrate the capability of the near-wall second-moment closure model for accurate resolution of the complex flow interaction bewteen the coolant jet and the mainstream. The near-wall second-moment anisotropic model provides better agreement in adiabatic film effectiveness prediction than the two-layer k-ε model.

Two Dimensional Simulations of Enhanced Heat Transfer in an Intermittently Grooved Channel

2000 ◽  
Author(s):  
M. Greiner ◽  
P. F. Fischer ◽  
H. M. Tufo
Keyword(s):  
Heat Transfer ◽  
Local Heat Transfer ◽  
Spectral Element ◽  
Local Heat ◽  
Navier Stokes ◽  
Computational Domain ◽  
Two Dimensional ◽  
Number Range ◽  
Element Technique ◽  
Outflow Boundary

Abstract Two-dimensional Navier-Stokes simulations of heat and momentum transport in an intermittently grooved passage are performed using the spectral element technique for the Reynolds number range 600 ≤ Re ≤ 1800. The computational domain has seven contiguous transverse grooves cut symmetrically into opposite walls, followed by a flat section with the same length. Periodic inflow/outflow boundary conditions are employed. The development and decay of unsteady flow is observed in the grooved and flat sections, respectively. The axial variation of the unsteady component of velocity is compared to the local heat transfer, shear stress and pressure gradient. The results suggest that intermittently grooved passages may offer even higher heat transfer for a given pumping power than the levels observed in fully grooved passages.

Sign in / Sign up

Export Citation Format

Share Document