Search NASA⌕ Search

SEARCH · Search NASA

Results for “Large eddy simulation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,477 records · Page 82
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Numerical investigation of initial condition effects on Rayleigh-Taylor instability with acceleration reversals

The influence of initial conditions on miscible incompressible baroclinically driven Rayleigh-Taylor instability undergoing nonuniform acceleration is explored computationally using an implicit large eddy simulation (ILES) technique. We consider here the particular case of evolution during multiple reversals of acceleration direction, where the flow is alternately statically stable or unstable. In the unstable phase, the flow is driven by the baroclinic release of potential energy, whereas in the stable phase, work is done against the density stratification with the energy exchange taking place by wavelike mechanisms. These dynamics are fundamentally different; here, we track the evolution of volume-averaged turbulent statistics that are most sensitive to changes in the distribution of spectral power and bandwidth of the initial conditions as the flow alternates between dynamical regimes due to acceleration reversal.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Self-contained filtered density function

The filtered density function (FDF) closure is extended to a “self-contained” format to include the subgrid-scale (SGS) statistics of all of the hydro-thermo-chemical variables in turbulent flows. These are the thermodynamic pressure, the specific internal energy, the velocity vector, and the composition field. In this format, the model is comprehensive and facilitates large-eddy simulation (LES) of flows at both low and high compressibility levels. A transport equation is developed for the joint pressure-energy-velocity-composition filtered mass density function (PEVC-FMDF). In this equation, the effect of convection appears in closed form. The coupling of the hydrodynamics and thermochemistry is modeled via a set of stochastic differential equation for each of the transport variables. This yields a self-contained SGS closure. We demonstrated how LES is conducted of a turbulent shear flow with transport of a passive scalar. Finally, the consistency of the PEVC-FMDF formulation is established, and its overall predictive capability is appraised via comparison with direct numerical simulation (DNS) data.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Inviscid criterion for decomposing scales

The proper scale decomposition in flows with significant density variations is not as straightforward as in incompressible flows, with many possible ways to define a “length scale.” A choice can be made according to the so-called inviscid criterion [Aluie, Physica D 24, 54 (2013)]. It is a kinematic requirement that a scale decomposition yield negligible viscous effects at large enough length scales. It has been proved [Aluie, Physica D 24, 54 (2013)] recently that a Favre decomposition satisfies the inviscid criterion, which is necessary to unravel inertial-range dynamics and the cascade. We present numerical demonstrations of those results. We also show that two other commonly used decompositions can violate the inviscid criterion and, therefore, are not suitable to study inertial-range dynamics in variable-density and compressible turbulence. Our results have practical modeling implication in showing that viscous terms in Large Eddy Simulations do not need to be modeled and can be neglected.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Partially averaged Navier-Stokes closure modeling for variable-density turbulent flow

We report this work extends the framework of the partially averaged Navier-Stokes (PANS) equations to variable-density flow, i.e., multimaterial and/or compressible mixing problems with density variations and production of turbulence kinetic energy by both shear and buoyancy mechanisms. The proposed methodology is utilized to derive the PANS BHR-LEVM closure. This includes a priori testing to analyze and develop guidelines toward the efficient selection of the parameters controlling the physical resolution and, consequently, the range of resolved scales of PANS. Two archetypal test-cases involving transient turbulence, hydrodynamic instabilities, and coherent structures are used to illustrate the accuracy and potential of the method: the Taylor-Green vortex at Reynolds number Re = 3000, and the Rayleigh-Taylor flow at Atwood number 0.5 and (Re) max ≈ 500. These representative problems, for which turbulence is generated by shear and buoyancy processes, constitute the initial validation space of the new model, and their results are comprehensively discussed in two subsequent studies. The computations indicate that PANS can accurately predict the selected flow problems, resolving only a fraction of the scales of large-eddy simulation and direct numerical simulation strategies. The results also reiterate that the physical resolution of the PANS model must guarantee that the key instabilities and coherent structures of the flow are resolved. The remaining scales can be modeled through an adequate turbulence scale-dependent closure.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Dimensionless parameters for cloudy Rayleigh-Bénard convection: Supersaturation, Damköhler, and Nusselt numbers

In steady-state Rayleigh-Bénard convection, heat is transported by turbulent thermal convection from the bottom, hot surface to the top, cold surface, leading to a height-independent sensible heat flux. When water vapor is present and cloud formation occurs, there is also an additional latent heat flux. Heat transport in cloudy Rayleigh-Bénard convection depends on turbulent flow as well as the microphysical state of the clouds: specifically, whether substantial supersaturations exist and whether cloud liquid water is removed through sedimentation/precipitation. In this article we bridge between the Rayleigh-Bénard convection literature and the atmospheric literature. We express the governing equations for cloudy convection in dimensionless form, thereby explicitly identifying the governing parameters relevant to the cloudy case, including Schmidt, Damköhler, supersaturation, and sedimentation numbers. We further connect to the atmospheric literature by obtaining a Nusselt number (dimensionless heat flux) for a cloud-convection system, directly from the conservation equations for temperature and water vapor. This flux has the same form as that identified by Zhang et al. [L. Zhang, K. L. Chong, and K.-Q. Xia, J. Fluid Mech. 874, 1041 (2019)] for convection with water vapor, but is extended to the cloudy case. For equal thermal and water vapor diffusivities, the flux corresponds to the widely used atmospheric quantities equivalent temperature and moist static energy. Using large eddy simulation (LES) of an idealized cloudy Rayleigh-Bénard convection system with fixed boundary conditions, we find that the equivalent heat flux (Nusselt number) is only weakly dependent on the microphysical details of the system, such as liquid water mixing ratio and cloud droplet number concentration. Finally, from the results, we show the vertical profiles of sensible and latent heat fluxes depend on the liquid water content, whereas the equivalent heat flux remains a constant throughout the height of the chamber.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Analysis of real-fluid thermodynamic effects on turbulent statistics in transcritical channel flows

Wall-bounded turbulence at high-pressure transcritical conditions with intense density fluctuations are encountered in many technical applications. In this study, we analyze the turbulent energy transport in transcritical channel flows specifically focusing on dissipation rate, turbulent kinetic energy budgets, heat fluxes, and momentum-fluctuation statistics; results from this analysis are used to guide the development of turbulent scaling laws. Furthermore, we find that the dissipation rate of turbulent kinetic energy is dominated by the enstrophy in the logarithmic layer, and the fluctuating viscosity results in the reduced tilting of the vortical structures and the attenuation of streamwise vorticity in the near-wall layer; the fluctuating viscosity attenuates the dissipation rate by reducing the shear strain and the enstrophy production. Local equilibrium of the turbulent kinetic energy exists in the logarithmic layer. We show that the real-fluid thermodynamic effects significantly change the turbulent heat flux correlated with the sweep and the ejection events; the density changes alter the turbulent transport and result in noticeable magnitudes of density-fluctuation-related momentum-fluctuation statistics. From these results, scaling laws for the turbulent length scales and turbulent kinetic energy budgets are proposed, thereby contributing to improvement of the wall models in large-eddy simulations and Reynolds-averaged Navier-Stokes (RANS) simulations.

42 ENGINEERING↗
DOE OSTI2024

Mechanisms and models of the turbulent boundary layers at transcritical conditions

Computational models for high-pressure transcritical turbulence in wall-modeled large-eddy simulation typically rely on wall-function models to overcome the need for resolving the boundary layer structures. However, the mechanisms and models of turbulent boundary layer at transcritical conditions remain poorly understood since the near-wall flow and heat transfer are significantly affected by the enhanced fluctuations and steep gradients in thermodynamic properties. Here, to address this issue, we study the mechanisms of transcritical turbulent boundary layers and wall-attached models at transcritical conditions. It is shown that the real-fluid variable-property effects are associated with the coupling between the near-wall cycle and the wall-normal coherent motions in the outer layer, resulting in the amplification of turbulent energy in the outer layer and noticeable energy transfer between the log-layer and the outer layer; hence, turbulence in the log-layer is modulated by the outer layer. Based on this underlying physical principle, we propose the characteristic velocity and length scales for the attached eddy at transcritical conditions, and extend the attached eddy model to transcritical turbulent boundary layers by introducing a mixed scaling that incorporates both inner and outer scalings. We show that the new characteristic scales and extended attached-eddy model perform well in characterizing the structures in transcritical turbulent boundary layers.

Li, Fangbo↗
DOE OSTI2016

Wind Farm Flow Modeling using an Input-Output Reduced-Order Model

Wind turbines in a wind farm operate individually to maximize their own power regardless of the impact of aerodynamic interactions on neighboring turbines. There is the potential to increase power and reduce overall structural loads by properly coordinating turbines. To perform control design and analysis, a model needs to be of low computational cost, but retains the necessary dynamics seen in high-fidelity models. The objective of this work is to obtain a reduced-order model that represents the full-order flow computed using a high-fidelity model. A variety of methods, including proper orthogonal decomposition and dynamic mode decomposition, can be used to extract the dominant flow structures and obtain a reduced-order model. In this paper, we combine proper orthogonal decomposition with a system identification technique to produce an input-output reduced-order model. This technique is used to construct a reduced-order model of the flow within a two-turbine array computed using a large-eddy simulation.

17 WIND ENERGY↗
DOE OSTI2017

A Validation Study of the Compressible Rayleigh–Taylor Instability Comparing the Ares and Miranda Codes

In this paper, the compressible Rayleigh–Taylor (RT) instability is studied by performing a suite of large eddy simulations (LES) using the Miranda and Ares codes. A grid convergence study is carried out for each of these computational methods, and the convergence properties of integral mixing diagnostics and late-time spectra are established. A comparison between the methods is made using the data from the highest resolution simulations in order to validate the Ares hydro scheme. We find that the integral mixing measures, which capture the global properties of the RT instability, show good agreement between the two codes at this resolution. The late-time turbulent kinetic energy and mass fraction spectra roughly follow a Kolmogorov spectrum, and drop off as k approaches the Nyquist wave number of each simulation. The spectra from the highest resolution Miranda simulation follow a Kolmogorov spectrum for longer than the corresponding spectra from the Ares simulation, and have a more abrupt drop off at high wave numbers. The growth rate is determined to be between around 0.03 and 0.05 at late times; however, it has not fully converged by the end of the simulation. Finally, we study the transition from direct numerical simulation (DNS) to LES. The highest resolution simulations become LES at around t/τ ≃ 1.5. Finally, to have a fully resolved DNS through the end of our simulations, the grid spacing must be 3.6 (3.1) times finer than our highest resolution mesh when using Miranda (Ares).

42 ENGINEERING↗
DOE OSTI2018

Using Machine Learning to Analyze Factors Determining Cycle-to-Cycle Variation in a Spark-Ignited Gasoline Engine

In this work, we have applied a machine learning (ML) technique to provide insights into the causes of cycle-to-cycle variation (CCV) in a gasoline spark-ignited (SI) engine. The analysis was performed on a set of large eddy simulation (LES) calculations of a single cylinder of a four-cylinder port-fueled SI engine. The operating condition was stoichiometric, without significant knock, at a load of 16 bar brake mean effective pressure (BMEP), at an engine speed of 2500 revolutions per minute. A total of 123 cycles was simulated. Of these, 49 were run in sequence, while 74 were run in parallel. For the parallel approach, each cycle is initialized with its own synthetic turbulent field to generate CCV, as part of another work performed by us. In the current work, we used 3D information from all 123 cycles to compute flame topology and pre-ignition flow-field metrics. We then evaluated correlations between these metrics, and peak cylinder pressure (PCP) employing an ML technique called random forest. The computed metrics form the inputs to the random forest model, and PCP is the output. This model captures the effect of all inputs, as well as interactions between them owing to its decision-tree structure. The goal of this work is to demonstrate (as a first step) that ML models can implicitly learn complex relationships between pre-ignition flow-fields, flame shapes, and the eventual outcome of the cycle (whether a cycle will be a high or a low cycle).

42 ENGINEERING↗
DOE OSTI2019

Oxy-Combustion Modeling for Direct-Fired Supercritical CO 2 Power Cycles

We report that supercritical CO 2 power cycles for fossil energy power generation will likely employ oxy-combustion at very high pressures, possibly exceeding 300 bar. At these high pressures, a direct fired oxy-combustor is more likely to behave like a rocket engine than any type of conventional gas turbine combustor. Issues such as injector design, wall heat transfer, and combustion dynamics may play a challenging role in combustor design. Computational fluid dynamics modeling will not only be useful, but may be a necessity in the combustor design process. To accurately model turbulent reacting flows, combustion submodels appropriate for the conditions of interest as defined by the turbulent time and length scales as well as chemical kinetic time scales are necessary. Here, this paper presents a comparison of various turbulence–chemistry interaction (TCI) modeling approaches on a canonical, single injector, direct-fired sCO 2 combustor. Large eddy simulation is used to model the turbulent combustion process with varying levels of injector oxygen concentration while comparing the effect of the combustion submodel on CO emissions and flame shape. While experimental data are not yet available to validate the simulations, the sensitivity of CO production and flame shape can be studied as a function of combustion modeling approach and oxygen concentration in an effort to better understand how to approach combustion modeling at these unique conditions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2019

Optimization of Kaplan Hydroturbine at Very Low Head With Rim-Driven Generator

The objective of the paper is to study the design and optimization of Kaplan hydroturbines for a very low head (less than 3 m), with a particular emphasis on the use of rim-drive electrical generators. The work is based on an experimental setup and computational fluid dynamics (CFD) analysis of various design parameters for maximum output power and efficiency. Two designs are presented in this paper. One is a 90-cm (35-in.) diameter vertical-oriented Kaplan hydroturbine system as an intended product capable of generating over 50 kW. The other is a smaller, 7.6-cm (3-in.) diameter horizontal-oriented system for prototyping and laboratory verification. Both are analyzed through CFD based on large eddy simulation (LES) of transient turbulence. Specific design for the runner and the stator, intake tube shape, as well as guide vanes upstream of the turbine was studied to get the most from the available head. The intent is to use 3D-printing manufacturing techniques, which may offer original design opportunities as well as the possibility of turbine and water conduit design customization as a function of the head and flow available from a specific site. Based on the CFD analysis, the 7.6-cm diameter system achieved the highest power output and the maximum efficiency at the rotational speed range of 1500–2000 rpm, while for the experimental testing, the optimum rotational speed range was 1000–1500 rpm. Because of the mismatch between CFD and experimental results, the CFD results were correlated due to the presence of air and friction; moreover, error and uncertainty analysis were presented for both methods. For the 90-cm case, the optimum performance was found at a rotational speed around 350 rpm according to the CFD results. Finally, investigating the shape of the intake tube of the hydroturbine setup can significantly increase the power output and the efficiency of the system.

Energy & Fuels↗
DOE OSTI2020

Sweeping Jet Film Cooling at High Blowing Ratio on a Turbine Vane

Abstract Experimental and numerical investigations were performed to study the effects of high blowing ratios and high freestream turbulence on sweeping jet film cooling. Experiments were conducted on a nozzle guide vane suction surface in a low-speed linear cascade at a range of blowing ratios of 0.5–3.5 and freestream turbulence of 0.6% and 14.3%. Infrared thermography was used to estimate the adiabatic cooling effectiveness. Thermal field and boundary layer measurements were conducted at a cross-plane at x/D = 12 downstream of the hole exit. Results were compared with a baseline 777-shaped hole and showed that the sweeping jet hole has improved cooling effectiveness at high blowing ratios (M > 1). The thermal field data revealed that the coolant separates from the surface at high blowing ratios for the 777-shaped hole while the coolant remains attached for the sweeping jet hole. Boundary layer measurements further confirmed that due to the sweeping motion of the jet, the effective jet momentum of the sweeping jet hole remains much lower than that of a 777-shaped hole. Thus, the coolant remains closer to the wall even at high blowing ratios. Large Eddy simulations (LES) were performed for both the sweeping jet and 777-shaped hole to evaluate the interaction between the coolant jet and the freestream in the near hole regions. Results showed that 777-shaped hole has a strong jetting at high blowing ratio that originates inside the hole breakout edges thus causing the jet to blow-off from the surface. In contrast, the sweeping jet does not show this behavior due to its internal geometry and flapping motion of the jet.

Engineering↗
DOE OSTI2021

Application of Scalar Filtered Density Function to Turbulent Flows Under Supercritical Condition

Abstract The scalar filtered density function (FDF) methodology is extended and employed for large eddy simulation (LES) of turbulent flows under supercritical condition. To describe real fluid behavior, the extended methodology incorporates the generalized heat and mass diffusion models along with real fluid thermodynamic relations which are derived using the cubic Peng–Robinson equation of state. These models are implemented within the stochastic differential equations comprising the scalar FDF transport. Simulations are conducted of a temporally developing mixing layer under supercritical condition and the results are assessed by comparing with data generated by direct numerical simulation (DNS) of the same layer. The consistency of the proposed FDF methodology is assessed. The LES-FDF predictions are shown to agree favorably with the DNS data and exhibit several key features pertaining to supercritical turbulent flows.

Energy & Fuels↗
DOE OSTI2021

Experimental Study on Flow Behavior and Heat Transfer Enhancement With Distinct Dimpled Gas Turbine Blade Internal Cooling Channel

Abstract The study presents the investigation on heat transfer distribution along a gas turbine blade internal cooling channel. Six different cases were considered in this study, using the smooth surface channel as a baseline. Three different dimples depth-to-diameter ratios with 0.1, 0.25, and 0.50 were considered. Different combinations of partial spherical and leaf dimples were also studied with the Reynolds numbers of 6,000, 20,000, 30,000, 40,000, and 50,000. In addition to the experimental investigation, the numerical study was conducted using a large eddy simulation (LES) to validate the data. It was found that the highest depth-to-diameter ratio showed the highest heat transfer rate. However, the deepest dimpled channel shows the highest pressure drop, which affects the overall thermal performance of the cooling channel. The results showed that the leaf dimpled surface is the best cooling channel based on the highest Reynolds number's heat transfer enhancement and friction factor. However, at the lowest Reynolds number, partial spherical dimples with a 0.25 depth-to-diameter ratio showed the highest thermal performance.

Energy & Fuels↗
DOE OSTI2022

Evaluation of Heat Transfer Enhancement on Rotational Gas Turbine Blade Internal Cooling Channel With Dimpled Surface

Abstract The present investigation represents the rotational effect on gas turbine blade internal cooling with a uniform heat flux of 2000 W/m2 at the bottom wall. The experiment was conducted with three different rpms, such as 300 rpm, 600 rpm, and 900 rpm, with Reynolds number (Re) ranging from 6000 to 50,000 with a two-pass cooling channel. The numerical investigation was conducted with the large eddy simulation (LES) technique to understand the rotational flow behavior of the cooling channel. Four distinct arrangements of dimpled cooling channel surfaces were considered with two different dimple shapes, i.e., partial spherical and leaf. It is found that the rotation effect, dimple arrangement, and design have significant influences on heat transfer. Results indicated that the partial spherical 1-row dimpled surface experienced the highest heat transfer coefficient and pressure drop. In contrast, the leaf-shaped dimpled cooling channel experienced the highest thermal efficiency.

Energy & Fuels↗
DOE OSTI2022

Exploring the Effects of Guide Vane on Dimpled Cooling Channel of Gas Turbine Blade

Abstract The purpose of the current study is to find out the heat transfer and pressure drop phenomena of the cooling channel with dimples and guide vane and compare the results with the no guide vane dimpled cooling channel. The first leg of the cooling channel is 490 mm, and the second leg is 460 mm. The two legs relate to the 180 deg turn region. The guide vane was inserted at the bend region of the dimpled cooling channel. The study was conducted with two different guide vanes geometry at two different orientations, i.e., U-guide vane with protrusion and depression orientation and curve-guide vane with protrusion and depression orientations both experimentally and numerically. The numerical study was performed with the large eddy simulation method. The result shows that for both stationary and rotational motion, the U-guide vane with depression experiences the highest thermal performance. The friction factor is comparatively higher for curve-guide vane with protrusion under stationary motion. However, under rotation, the curve-protrusion guide encounters the highest friction factor, which is higher compared to the no guide vane cooling channel.

Energy & Fuels↗
DOE OSTI2024

Validation Assessment of Turbulent Reacting Flow Model Using The Area-Validation Metric on Medium-Scale Methanol Pool Fire Results

Accident analysis and ensuring power plant safety are pivotal in the nuclear energy sector. Significant strides have been achieved over the past few decades regarding fire protection and safety, primarily centered on design and regulatory compliance. Yet, after the Fukushima accident a decade ago, the imperative to enhance measures against fire, internal flooding, and power loss has intensified. Hence, a comprehensive, multilayered protection strategy against severe accidents is needed. Consequently, gaining a deeper insight into pool fires and their behavior through extensive validated data can greatly aid in improving these measures using advanced validation techniques. A model validation study was performed at Sandia National Laboratories in which a 30-cm diameter methanol pool fire was modeled using the SIERRA/Fuego turbulent reacting flow code. This validation study used a standard validation experiment to compare model results against, and conclusions have been published. The fire was modeled with a Large Eddy Simulation (LES) turbulence model with subgrid turbulent kinetic energy closure. Combustion was modeled using a strained laminar flamelet library approach. Radiative heat transfer was accounted for with a model utilizing the gray-gas approximation. In the present study, additional validation analysis is performed using the area validation metric (AVM). These activities are done on multiple datasets involving different variables and temporal/spatial ranges and intervals. In conclusion, the results provide insight into the use of the area validation metric on such temporally varying datasets and the importance of physics-aware use of the metric for proper analysis.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗