Upward versus downward non-Boussinesq turbulent fountains

2019 ◽  
Vol 867 ◽  
pp. 374-391 ◽  
Author(s):  
Samuel Vaux ◽  
Rabah Mehaddi ◽  
Olivier Vauquelin ◽  
Fabien Candelier
Keyword(s):  
Steady State ◽  
Numerical Simulations ◽  
Density Ratio ◽  
Large Eddy Simulations ◽  
Linear Scaling ◽  
Power Exponent ◽  
Maximum Height ◽  
Source Density ◽  
Wide Range ◽  
Large Eddy

Turbulent miscible fountains discharged vertically from a round source into quiescent uniform unbounded environments of density $\unicode[STIX]{x1D70C}_{0}$ are investigated numerically using large-eddy simulations. Both upward and downward fountains are considered. The numerical simulations cover a wide range of the density ratio $\unicode[STIX]{x1D70C}_{i}/\unicode[STIX]{x1D70C}_{0}$, where $\unicode[STIX]{x1D70C}_{i}$ is the source density of the released fluid. These simulations are used to evaluate how the initial maximum height $H_{i}$ and the steady state height $H_{ss}$ of the fountains are affected by large density contrasts, i.e. in the general non-Boussinesq case. For both upward and downward non-Boussinesq fountains, the ratio $\unicode[STIX]{x1D706}=H_{i}/H_{ss}$ remains close to $1.45$, as usually observed for Boussinesq fountains. However the Froude (linear) scaling originally introduced by Turner (J. Fluid Mech., vol. 26 (4), 1966, pp. 779–792) for Boussinesq fountains is no longer valid to determine the steady fountain height. The ratio between $H_{ss}$ and the height predicted by the Turner’s relation turns out to be proportional to $(\unicode[STIX]{x1D70C}_{i}/\unicode[STIX]{x1D70C}_{0})^{n}$. Remarkably, it is found that the power exponent $n$ differs following the direction in which the buoyant fluid is released ($n=1/2$ for downward fountains and $n=3/4$ for upward fountains). This new result demonstrates that for non-Boussinesq turbulent fountains the configurations heavy/light and light/heavy are not equivalent. Finally, scalings are proposed for fountains, regardless of the direction (upwards and downwards) and of the density difference (Boussinesq and non-Boussinesq).

Large Eddy Simulations of Staggered Parallel-Plate Flows

2005 ◽  
Author(s):  
M. V. Pham ◽  
F. Plourde ◽  
S. K. Doan
Keyword(s):  
Heat Transfer ◽  
Numerical Simulations ◽  
Heat Transfer Enhancement ◽  
Parallel Plate ◽  
Large Scale ◽  
Primary Objective ◽  
Transfer Enhancement ◽  
Turbulent Structures ◽  
Wide Range ◽  
Large Eddy

Heat transfer enhancement is a subject of major concern in numerous fields of industry and research. Having received undivided attention over the years, it is still studied worldwide. Given the exponential growth of computing power, large-scale numerical simulations are growing steadily more realistic, and it is now possible to obtain accurate time-dependent solutions with far fewer preliminary assumptions about the problems. As a result, an increasingly wide range of physics is now open for exploration. More specifically, it is time to take full advantage of large eddy simulation technique so as to describe heat transfer in staggered parallel-plate flows. In fact, from simple theory through experimental results, it has been demonstrated that surface interruption enhances heat transfer. Staggered parallel-plate geometries are of great potential interest, and yet many numerical works dedicated to them have been tarnished by excessively simple assumptions. That is to say, numerical simulations have generally hypothesized lengthwise periodicity, even though flows are not periodic; moreover, the LES technique has not been employed with sufficient frequency. Actually, our primary objective is to analyze turbulent influence with regard to heat transfers in staggered parallel-plate fin geometries. In order to do so, we have developed a LES code, and numerical results are compared with regard to several grid mesh resolutions. We have focused mainly upon identification of turbulent structures and their role in heat transfer enhancement. Another key point involves the distinct roles of boundary restart and the vortex shedding mechanism on heat transfer and friction factor.

Numerical Modeling of Soot Production in Aero-Engine Combustors Using Large Eddy Simulations

2015 ◽  
Author(s):  
B. Franzelli ◽  
E. Riber ◽  
B. Cuenot ◽  
M. Ihme
Keyword(s):  
Numerical Simulations ◽  
Flame Structure ◽  
Precursor Chemistry ◽  
Turbulent Flames ◽  
Large Eddy Simulations ◽  
Reduced Chemistry ◽  
Tabulated Chemistry ◽  
Large Eddy ◽  
Aero Engine ◽  
Aero Engines

Numerical simulations are regarded as an essential tool for improving the design of combustion systems since they can provide information that is complementary to experiments. However, although numerical simulations have already been successfully applied to the prediction of temperature and species concentration in turbulent flames, the production of soot is far from being conclusive due to the complexity of the processes involved in soot production. In this context, first Large Eddy Simulations (LES) of soot production in turbulent flames are reported in the literature in laboratory-scale configurations, thereby confirming the feasibility of the approach. However numerous modeling and numerical issues have not been completely solved. Moreover, validation of the models through comparisons with measurements in realistic complex flows typical of aero-engines is still rare. This work therefore proposes to evaluate the LES approach for the prediction of soot production in an experimental swirl-stabilized non-premixed ethylene/air aero-engine combustor, for which soot and flame data are available. Two simulations are carried out using a two-equation soot model to compare the performance of a hybrid chemical description (reduced chemistry for the flame structure/tabulated chemistry for soot precursor chemistry) to a classical full tabulation method. Discrepancies of soot concentration between the two LES calculations will be analyzed and the sensitivity to the chemical models will be investigated.

scholarly journals Alternating direction implicit type preconditioners for the steady state inhomogeneous Vlasov equation

2017 ◽  
Vol 83 (1) ◽  
Author(s):  
Markus Gasteiger ◽  
Lukas Einkemmer ◽  
Alexander Ostermann ◽  
David Tskhakaya
Keyword(s):  
Steady State ◽  
Numerical Simulations ◽  
Vlasov Equation ◽  
Numerical Solutions ◽  
Splitting Method ◽  
Time Integration ◽  
Computational Effort ◽  
Alternating Direction Implicit ◽  
Alternating Direction ◽  
Wide Range

The purpose of the current work is to find numerical solutions of the steady state inhomogeneous Vlasov equation. This problem has a wide range of applications in the kinetic simulation of non-thermal plasmas. However, the direct application of either time stepping schemes or iterative methods (such as Krylov-based methods such as the generalized minimal residual method (GMRES) or relaxation schemes) is computationally expensive. In the former case the slowest time scale in the system forces us to perform a long time integration while in the latter case a large number of iterations is required. In this paper we propose a preconditioner based on an alternating direction implicit type splitting method. This preconditioner is then combined with both GMRES and Richardson iteration. The resulting numerical schemes scale almost ideally (i.e. the computational effort is proportional to the number of grid points). Numerical simulations conducted show that this can result in a speed-up of close to two orders of magnitude (even for intermediate grid sizes) with respect to the unpreconditioned case. In addition, we discuss the characteristics of these numerical methods and show the results for a number of numerical simulations.

Effect of Blowing Ratio in the Near-Stagnation Region of a Three-Row Leading Edge Film Cooling Geometry Using Large Eddy Simulations

2009 ◽  
Author(s):  
Sai Shrinivas Sreedharan ◽  
Danesh K. Tafti
Keyword(s):  
Film Cooling ◽  
Density Ratio ◽  
Leading Edge ◽  
Large Eddy Simulations ◽  
Stagnation Region ◽  
Stagnation Line ◽  
Blowing Ratio ◽  
Large Eddy ◽  
Adiabatic Effectiveness ◽  
Compound Angle

Computational studies are carried out using Large Eddy Simulations (LES) to investigate the effect of coolant to mainstream blowing ratio in a leading edge region of a film cooled vane. The three row leading edge vane geometry is modeled as a symmetric semi-cylinder with a flat afterbody. One row of coolant holes is located along the stagnation line and the other two rows of coolant holes are located at ±21.3° from the stagnation line. The coolant is injected at 45° to the vane surface with 90° compound angle injection. The coolant to mainstream density ratio is set to unity and the freestream Reynolds number based on leading edge diameter is 32000. Blowing ratios (B.R.) of 0.5, 1.0, 1.5, and 2.0 are investigated. It is found that the stagnation cooling jets penetrate much further into the mainstream, both in the normal and lateral directions, than the off-stagnation jets for all blowing ratios. Jet dilution is characterized by turbulent diffusion and entrainment. The strength of both mechanisms increases with blowing ratio. The adiabatic effectiveness in the stagnation region initially increases with blowing ratio but then generally decreases as the blowing ratio increases further. Immediately downstream of off-stagnation injection, the adiabatic effectiveness is highest at B.R. = 0.5. However, further downstream the larger mass of coolant injected at higher blowing ratios, in spite of the larger jet penetration and dilution, increases the effectiveness with blowing ratio.

Large Eddy Simulations of Flow and Heat Transfer in Rotating Ribbed Duct Flows

2005 ◽  
Vol 127 (5) ◽  
pp. 486-498 ◽  
Author(s):  
Mayank Tyagi ◽  
Sumanta Acharya
Keyword(s):  
Heat Transfer ◽  
Temperature Field ◽  
Large Scale ◽  
Density Ratio ◽  
Low Frequency ◽  
Turbulent Energy ◽  
Large Eddy Simulations ◽  
Temperature Fields ◽  
Periodic Domain ◽  
Large Eddy

Large eddy simulations are performed in a periodic domain of a rotating square duct with normal rib turbulators. Both the Coriolis force as well as the centrifugal buoyancy forces are included in this study. A direct approach is presented for the unsteady calculation of the nondimensional temperature field in the periodic domain. The calculations are performed at a Reynolds number (Re) of 12,500, a rotation number (Ro) of 0.12, and an inlet coolant-to-wall density ratio Δρ/ρ of 0.13. The predicted time and space-averaged Nusselt numbers are shown to compare satisfactorily with the published experimental data. Time sequences of the vorticity components and the temperature fields are presented to understand the flow physics and the unsteady heat transfer behavior. Large scale coherent structures are seen to play an important role in the mixing and heat transfer. The temperature field appears to contain a low frequency mode that extends beyond a single inter-rib geometric module, and indicates the necessity of using at least two inter-rib modules for streamwise periodicity to be satisfied. Proper orthogonal decomposition (POD) of the flowfield indicates a low dimensionality of this system with almost 99% of turbulent energy in the first 80 POD modes.

Large Eddy Simulations of Flow and Heat Transfer in Rotating Ribbed Duct Flows

2004 ◽  
Author(s):  
Mayank Tyagi ◽  
Sumanta Acharya
Keyword(s):  
Heat Transfer ◽  
Temperature Field ◽  
Large Scale ◽  
Density Ratio ◽  
Low Frequency ◽  
Turbulent Energy ◽  
Large Eddy Simulations ◽  
Temperature Fields ◽  
Periodic Domain ◽  
Large Eddy

Large eddy simulations are performed in a periodic domain of a rotating square duct with normal rib turbulators. Both the Coriolis force as well as the centrifugal buoyancy force are included in this study. A direct approach is presented for the unsteady calculation of the non-dimensional temperature field in the periodic domain. The calculations are performed at a Reynolds number (Re) of 12, 500, a Rotation number (Ro) of 0.12 and an inlet coolant-to-wall density ratio (Δρ/ρ) of 0.13. The time-averaged Nusselt numbers compare satisfactorily with the data of Wagner et al. (J. Turbomachinery, Vol. 114, pp. 847–857). Time-sequences of the vorticity components and the temperature fields are presented to understand the flow physics and the unsteady heat transfer processes. Large scale coherent structures are seen to play an important role in the mixing and heat transfer. The temperature field appears to contain a low frequency mode that extends beyond a single inter-rib geometric module, and indicates the necessity of using at least two inter-rib modules for streamwise periodicity to be satisfied. Proper orthogonal decomposition (POD) of 200 snapshots indicates a low dimensionality of this system with almost 99% of turbulent energy in the first 80 POD modes.

Numerical Simulation of Free-Surface Turbulence

1994 ◽  
Vol 47 (6S) ◽  
pp. S163-S165
Author(s):  
Douglas G. Dommermuth ◽  
Rebecca C. Y. Mui
Keyword(s):  
Numerical Simulation ◽  
Numerical Modeling ◽  
Free Surface ◽  
Numerical Simulations ◽  
Free Surface Flows ◽  
Direct Numerical Simulations ◽  
Large Eddy Simulations ◽  
Surface Flows ◽  
Large Eddy ◽  
Free Surface Turbulence

Direct numerical simulations and large-eddy simulations of turbulent free-surface flows are currently being performed to investigate the roughening of the surface, and the scattering, radiation, and dissipation of waves by turbulence. The numerical simulation of turbulent free-surface flows is briefly reviewed. The numerical, modeling, and hardware issues are discussed.

scholarly journals A Closer Look at Boundary Layer Inversion in Large-Eddy Simulations and Bulk Models: Buoyancy-Driven Case

2015 ◽  
Vol 72 (2) ◽  
pp. 728-749 ◽  
Author(s):  
Pierre Gentine ◽  
Gilles Bellon ◽  
Chiel C. van Heerwaarden
Keyword(s):  
Boundary Layer ◽  
Convective Boundary Layer ◽  
Weather Prediction ◽  
Inversion Layer ◽  
Large Eddy Simulations ◽  
Buoyancy Flux ◽  
Convective Boundary ◽  
Strong Inversion ◽  
Wide Range ◽  
Large Eddy

Abstract The inversion layer (IL) of a clear-sky, buoyancy-driven convective boundary layer is investigated using large-eddy simulations covering a wide range of convective Richardson numbers. A new model of the IL is suggested and tested. The model performs better than previous first-order models of the entrainment and provides physical insights into the main controls of the mixed-layer and IL growths. A consistent prognostic equation of the IL growth is derived, with explicit dependence on the position of the minimum buoyancy flux, convective Richardson number, and relative stratification across the inversion G. The IL model expresses the interrelationship between the position and magnitude of the minimum buoyancy flux and inversion-layer depth. These relationships emphasize why zero-order jump models of the convective boundary layer perform well under a strong inversion and show that these models miss the additional parameter G to fully characterize the entrainment process under a weak inversion. Additionally, the position of the minimum buoyancy flux within the new IL model is shown to be a key component of convective boundary layer entrainment. The new IL model is sufficiently simple to be used in numerical weather prediction or general circulation models as a way to resolve the IL in a low-vertical-resolution model.

scholarly journals The Aerodynamics of the Curled Wake: A Simplified Model in View of Flow Control

2018 ◽  
Author(s):  
Luis A. Martínez-Tossas ◽  
Jennifer Annoni ◽  
Paul A. Fleming ◽  
Matthew J. Churchfield
Keyword(s):  
Boundary Layer ◽  
Steady State ◽  
Wind Farm ◽  
Stokes Equations ◽  
Large Eddy Simulations ◽  
Simplified Model ◽  
Navier Stokes ◽  
Navier Stokes Equations ◽  
Actuator Disk ◽  
Large Eddy

Abstract. When a wind turbine is yawed, the shape of the wake changes and a curled wake profile is generated. The curled wake has drawn a lot of interest because of its aerodynamic complexity and applicability to wind farm controls. The main mechanism for the creation of the curled wake has been identified in the literature as a collection of vortices that are shed from the rotor plane when the turbine is yawed. This work extends that idea by using aerodynamic concepts to develop a control-oriented model for the curled wake based on approximations to the Navier-Stokes equations. The model is tested and compared to large-eddy simulations using actuator disk and line models. The model is able to capture the curling mechanism for a turbine under uniform inflow and in the case of a neutral atmospheric boundary layer. The model is then tested inside the FLOw Redirection and Induction in Steady State framework and provides excellent agreement with power predictions for cases with two and three turbines in a row.

Large Eddy Simulations of Discrete Hole Film Cooling With Plenum Inflow Orientation Effects

2012 ◽  
Vol 135 (1) ◽  
Author(s):  
Sumanta Acharya ◽  
David Houston Leedom
Keyword(s):  
Film Cooling ◽  
Flow Behavior ◽  
Density Ratio ◽  
Large Eddy Simulations ◽  
Delivery Tube ◽  
Cooling Effectiveness ◽  
Freestream Velocity ◽  
Plenum Chamber ◽  
Discharge Coefficients ◽  
Large Eddy

Large eddy simulations of film cooling from a discrete hole inclined 35 deg and fed by a plenum chamber are performed at a density ratio of 2 and blowing ratios from 0.5 to 2.0. Cylindrical holes at a length to diameter ratio of 1.75 and 3.5 are simulated issuing into a crossflow at a Reynolds number of approximately 16,000 based on freestream velocity and hole diameter. In addition to the baseline case of vertical inflow into the plenum, flow orientation into the plenum chamber parallel to and perpendicular to the mainstream flow are investigated. The predicted results are validated with reported measurements of the flow field and surface adiabatic effectiveness. Results show that the longer delivery tubes (L/D = 3.5) have higher cooling effectiveness except in the very near field of the coolant hole. The flow orientation in the plenum is demonstrated to have a significant effect on cooling effectiveness and on flow behavior in the delivery tube and downstream of the hole. The perpendicular plenum inflow exhibits the lowest cooling effectiveness, the lowest discharge coefficients, asymmetric jetting behavior, swirl, and a low-velocity core at the exit of the delivery tube. The parallel plenum flow orientation is shown to exhibit the highest cooling effectiveness and discharge coefficients.

Sign in / Sign up

Export Citation Format

Share Document