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 523 records · Page 29

A Highly Accurate Technique for the Treatment of Flow Equations at the Polar Axis in Cylindrical Coordinates using Series Expansions

Numerical methods for solving the flow equations in cylindrical or spherical coordinates should be able to capture the behavior of the exact solution near the regions where the particular form of the governing equations is singular. In this work we focus on the treatment of these numerical singularities for finite-differences methods by reinterpreting the regularity conditions developed in the context of pseudo-spectral methods. A generally applicable numerical method for treating the singularities present at the polar axis, when nonaxisymmetric flows are solved in cylindrical, coordinates using highly accurate finite differences schemes (e.g., Pade schemes) on non-staggered grids, is presented. Governing equations for the flow at the polar axis are derived using series expansions near r=0. The only information needed to calculate the coefficients in these equations are the values of the flow variables and their radial derivatives at the previous iteration (or time) level. These derivatives, which are multi-valued at the polar axis, are calculated without dropping the accuracy of the numerical method using a mapping of the flow domain from (0,R)*(0,2pi) to (-R,R)*(0,pi), where R is the radius of the computational domain. This allows the radial derivatives to be evaluated using high-order differencing schemes (e.g., compact schemes) at points located on the polar axis. The proposed technique is illustrated by results from simulations of laminar-forced jets and turbulent compressible jets using large eddy simulation (LES) methods. In term of the general robustness of the numerical method and smoothness of the solution close to the polar axis, the present results compare very favorably to similar calculations in which the equations are solved in Cartesian coordinates at the polar axis, or in which the singularity is removed by employing a staggered mesh in the radial direction without a mesh point at r=0, following the method proposed recently by Mohseni and Colonius (1). Extension of the method described here for incompressible flows or for any other set of equations that are solved on a non-staggered mesh in cylindrical or spherical coordinates with finite-differences schemes of various level of accuracy is immediate.

Constantinescu, George S.↗

Terminal Area Simulation System User's Guide - Version 10.0

The Terminal Area Simulation System (TASS) is a three-dimensional, time-dependent, large eddy simulation model that has been developed for studies of wake vortex and weather hazards to aviation, along with other atmospheric turbulence, and cloud-scale weather phenomenology. This document describes the source code for TASS version 10.0 and provides users with needed documentation to run the model. The source code is programed in Fortran language and is formulated to take advantage of vector and efficient multi-processor scaling for execution on massively-parallel supercomputer clusters. The code contains different initialization modules allowing the study of aircraft wake vortex interaction with the atmosphere and ground, atmospheric turbulence, atmospheric boundary layers, precipitating convective clouds, hail storms, gust fronts, microburst windshear, supercell and mesoscale convective systems, tornadic storms, and ring vortices. The model is able to operate in either two- or three-dimensions with equations numerically formulated on a Cartesian grid. The primary output from the TASS is time-dependent domain fields generated by the prognostic equations and diagnosed variables. This document will enable a user to understand the general logic of TASS, and will show how to configure and initialize the model domain. Also described are the formats of the input and output files, as well as the parameters that control the input and output.

Switzer, George F.↗

Modeling and Analysis of Realistic Fire Scenarios in Spacecraft

An accidental fire inside a spacecraft is an unlikely, but very real emergency situation that can easily have dire consequences. While much has been learned over the past 25+ years of dedicated research on flame behavior in microgravity, a quantitative understanding of the initiation, spread, detection and extinguishment of a realistic fire aboard a spacecraft is lacking. Virtually all combustion experiments in microgravity have been small-scale, by necessity (hardware limitations in ground-based facilities and safety concerns in space-based facilities). Large-scale, realistic fire experiments are unlikely for the foreseeable future (unlike in terrestrial situations). Therefore, NASA will have to rely on scale modeling, extrapolation of small-scale experiments and detailed numerical modeling to provide the data necessary for vehicle and safety system design. This paper presents the results of parallel efforts to better model the initiation, spread, detection and extinguishment of fires aboard spacecraft. The first is a detailed numerical model using the freely available Fire Dynamics Simulator (FDS). FDS is a CFD code that numerically solves a large eddy simulation form of the Navier-Stokes equations. FDS provides a detailed treatment of the smoke and energy transport from a fire. The simulations provide a wealth of information, but are computationally intensive and not suitable for parametric studies where the detailed treatment of the mass and energy transport are unnecessary. The second path extends a model previously documented at ICES meetings that attempted to predict maximum survivable fires aboard space-craft. This one-dimensional model implies the heat and mass transfer as well as toxic species production from a fire. These simplifications result in a code that is faster and more suitable for parametric studies (having already been used to help in the hatch design of the Multi-Purpose Crew Vehicle, MPCV).

life support↗

Summary of the 5th Propulsion Aerodynamics Workshop Nozzle Test Case: Heated Nozzle Exhaust Passing Over A Film-Cooled Plate

This paper summarizes findings from the fifth AIAA Propulsion Aerodynamics Workshop nozzle test case. The experimental configuration examined by workshop participants was a subsonic nozzle with a square exit blowing a heated exhaust over a film-cooled plate. Computational fluid dynamics solutions were obtained and compared to experimental measurements of flow velocities and temperatures, as well as plate surface temperatures. The heated nozzle operated with a Mach 0.3 exit flow at a static temperature ratio of 2.7. Cooling air blowing ratios of 0, 1, and 2 were considered. Computational meshes for the flow domain were provided for participants. Workshop participants from eight separate organizations represented government, industry, and academia. A variety of flow solutions were obtained: Most of the flow solutions employed a Reynolds-averaged Navier-Stokes (RANS) approach, but scale-resolving simulations were used in some cases, including wall-modeled large-eddy simulation (LES), and a few hybrid RANS-LES approaches. One Lattice-Boltzmann solver was employed. For heat transfer to the test article, many participants used a conjugate heat transfer approach. Effects of mesh sensitivity, flow solution approach, and wall heat transfer approach are considered. In general, a fully three-dimensional conjugate heat transfer approach enabled better prediction of surface temperatures than simpler wall temperature boundary treatments. Also, the scale resolving approaches provided better prediction of static temperatures in the flow immediately above the plate consisting of the hot jet exhaust boundary layer interacting with the cooling film. Comparisons of plate surface friction drag and heat transfer obtained from the computations are also presented.

nozzle↗

Summary of the 5th Propulsion Aerodynamics Workshop Nozzle Test Case: Heated Nozzle Exhaust Passing Over A Film-Cooled Plate

This paper summarizes findings from the fifth AIAA Propulsion Aerodynamics Workshop nozzle test case. The experimental configuration examined by workshop participants was a subsonic nozzle with a square exit blowing a heated exhaust over a film-cooled plate. Computational fluid dynamics solutions were obtained and compared to experimental measurements of flow velocities and temperatures, as well as plate surface temperatures. The heated nozzle operated with a Mach 0.3 exit flow at a static temperature ratio of 2.7. Cooling air blowing ratios of 0, 1, and 2 were considered. Computational meshes for the flow domain were provided for participants. Workshop participants from eight separate organizations represented government, industry, and academia. A variety of flow solutions were obtained: Most of the flow solutions employed a Reynolds-averaged Navier-Stokes (RANS) approach, but scale-resolving simulations were used in some cases, including wall-modeled large-eddy simulation (LES), and a few hybrid RANS-LES approaches. One Lattice-Boltzmann solver was employed. For heat transfer to the test article, many participants used a conjugate heat transfer approach. Effects of mesh sensitivity, flow solution approach, and wall heat transfer approach are considered. In general, a fully three-dimensional conjugate heat transfer approach enabled better prediction of surface temperatures than simpler wall temperature boundary treatments. Also, the scale resolving approaches provided better prediction of static temperatures in the flow immediately above the plate consisting of the hot jet exhaust boundary layer interacting with the cooling film. Comparisons of plate surface friction drag and heat transfer obtained from the computations are also presented.

nozzle↗

Identifying Meteorological Influences on Marine Low Cloud Mesoscale Morphology Using Satellite Classifications

Marine low cloud mesoscale morphology in the southeastern Pacific Ocean is analyzed using a large dataset of machine-learning generated classifications spanning three years. Meteorological variables and cloud properties are composited 10by mesoscale cloud type, showing distinct meteorological regimes of marine low cloud organization from the tropics to the midlatitudes. The presentation of mesoscale cellular convection, with respect to geographic distribution, boundary layer structure, and large-scale environmental conditions, agrees with prior knowledge. Two tropical and subtropical cumuliform boundary layer regimes, suppressed cumulus and clustered cumulus, are studied in detail. The patterns in precipitation, circulation, column water vapor, and cloudiness are consistent with the representation of marine shallow mesoscale convective 15 self-aggregation by large eddy simulations of the boundary layer. Although they occur under similar large-scale conditions, the suppressed and clustered low cloud types are found to be well-separated by variables associated with low-level mesoscale circulation, with surface wind divergence being the clearest discriminator between them, whether reanalysis or satellite observations are used. Clustered regimes are associated with surface convergence and suppressed regimes are associated with surface divergence.

Johannes Mohrmann↗

Simulation of turbulent shear flows at Stanford and NASA-Ames - What can we do and what have we learned?

The capabilities and limitations of large eddy simulation (LES) and full turbulence simulation (FTS) are outlined. It is pointed out that LES, although limited at the present time by the need for periodic boundary conditions, produces large-scale flow behavior in general agreement with experiments. What is more, FTS computations produce small-scale behavior that is consistent with available experiments. The importance of the development work being done on the National Aerodynamic Simulator is emphasized. Studies at present are limited to situations in which periodic boundary conditions can be applied on boundaries of the computational domain where the flow is turbulent.

Reynolds, W. C.↗

Recent progress on coarse graining simulations

We focus on coarse graining simulations based on the primary conservation equations, effectively codesigned physics and algorithms, and low-Mach-number corrected (LMC) hydrodynamics. Simulation methods involve LANL’s x-Radiation-Adaptive-Grid-Eulerian Large-Eddy Simulation, Besnard-Harlow-Rauenzahn (BHR) Reynolds-Averaged Navier-Stokes (RANS) approach, and Dynamic BHR – a paradigm bridging RANS and LES. A relevant question addressed relates to whether 3D RANS and RANS/LES hybrids – the industry standards for aerospace and automotive research, are presently relevant for practical variable-density applications involving shocked and accelerated interface instabilities. Furthermore, recent simulations of the GaTECH inclined mixing-layer shock-tube and NIF ICF-capsule experiments are used to demonstrate issues, challenges, and potential for 3D coarse grained LMC simulation strategies for robustly simulating complex transitional and coupled hydrodynamics-multiphysics with coarser resolution. Present LES readiness to provide accurate predictions at scale is demonstrated – whereas 3D RANS and RANS/LES bridging do not appear impactful in this context.

42 ENGINEERING↗

Idealized simulations of wind farm interactions with intermittent turbulence in stable boundary layer conditions

Stable atmospheric boundary layer conditions typically correspond to weak turbulence levels, but intermittent periods of elevated turbulence can occur during otherwise quiescent conditions. The interaction between intermittent turbulence and wind turbines is not well understood because of sparse observations, as well as the difficulty in realistically resolving small-scale turbulence during strongly stable conditions with numerical simulations. In this study, an explicit filtering and reconstruction approach for large-eddy simulation (LES) is used to simulate weakly and strongly stable conditions, with surface cooling rates of −0.2 and −2.0 K h −1 , respectively. This approach can sustain resolved background turbulence at relatively coarse grid spacing and stronger stratification compared to conventional closures, permitting more realistic intermittent stable boundary layer (SBL) turbulence. The idealized LES capability of the Weather Research and Forecasting model is employed with turbine rotors parameterized using generalized actuator disks to examine (1) how the presence of turbine wakes affects SBL evolution and (2) the effect of intermittent turbulence on power production and wake recovery. Wakes increase mixing and deepen the SBL, with a stronger effect under strongly stable conditions, primarily because the SBL is shallower and closer to the top of the wind turbine rotor layer. Intermittent turbulence does not have a significant impact on mean power generation and wake recovery because the relevant intermittent turbulent structures in this study only affect the bottom half of the rotor disk. Power production is, however, more variable during periods of elevated turbulence, demonstrating the impact of SBL intermittency. This study uses an idealized configuration, focusing on LES model performance and physical understanding, with the goal of informing future simulations of the conditions observed during the American Wake Experiment.

Energy - Wind↗

Investigating the interactions between wakes and floating wind turbines using FAST.Farm

As floating offshore wind progresses to commercial maturity, wake and array effects across a farm of floating offshore wind turbines (FOWTs) will become increasingly important. While wakes of land-based and bottom-fixed offshore wind turbines have been extensively studied, only recently has this topic become relevant for floating turbines. This work presents an investigation of the mutual interaction between the motions of floating wind turbines and wakes using FAST.Farm. While FAST.Farm has been extensively validated across a wide range of conditions, it has never been validated for FOWT applications. Hence, in the first part of this work, we validate FAST.Farm by comparing simulations of a single FOWT against high-fidelity results from large-eddy simulations available in the literature. The validation is based on wake meandering, mean wake deflection, and velocity deficit at different downstream locations. This validation showed that the original axisymmetric (polar) wake model of FAST.Farm overpredicts the vertical wake deflection induced by shaft tilt and floater pitch, while the new curled wake model is capable of properly capturing the vertical wake deflection. In the second part, we use FAST.Farm to analyze a small three-unit array of FOWTs with a spacing of 7 diameters across a wide range of environmental conditions. The same National Renewable Energy Laboratory 5 MW reference wind turbine atop the OC4-DeepCwind semisubmersible is adopted for the three FOWTs and for the validation against high-fidelity simulations. To assess the effect of the floating substructure, we compare the power production, tower-base moments, and blade-root moments obtained for the floating turbines with the results obtained in a fixed-bottom configuration. The main differences introduced by the floating substructure are the motions induced by the waves, the change in the natural frequencies of the tower caused by differences in the boundary condition at its base, and the larger vertical deflection of the wake deficit due to the mean pitch of the platform. The impact of these differences, as well as other minor effects, are analyzed in detail.

16 TIDAL AND WAVE POWER↗

High speed turbulent reacting flows: DNS and LES

Work on understanding the mechanisms of mixing and reaction in high speed turbulent reacting flows was continued. Efforts, in particular, were concentrated on taking advantage of modern computational methods to simulate high speed turbulent flows. In doing so, two methodologies were used: large eddy simulations (LES) and direct numerical simulations (DNS). In the work related with LES the objective is to study the behavior of the probability density functions (pdfs) of scalar properties within the subgrid in reacting turbulent flows. The data base obtained by DNS for a detailed study of the pdf characteristics within the subgrid was used. Simulations are performed for flows under various initializations to include the effects of compressibility on mixing and chemical reactions. In the work related with DNS, a two-dimensional temporally developing high speed mixing layer under the influence of a second-order non-equilibrium chemical reaction of the type A + B yields products + heat was considered. Simulations were performed with different magnitudes of the convective Mach numbers and with different chemical kinetic parameters for the purpose of examining the isolated effects of the compressibility and the heat released by the chemical reactions on the structure of the layer. A full compressible code was developed and utilized, so that the coupling between mixing and chemical reactions is captured in a realistic manner.

Givi, Peyman↗

N-S Simulations of Crow-Type Instabilities in Vortex Wake

Three-dimensional Navier-Stokes simulations of the Crow instability of wake vortices are conducted using large-eddy simulations without background turbulence. Sinusoidal displacement has been specified as the initial perturbation for the vortex system. The results have shown that the minimum Crow instability wavelength is about one vortex spacing shorter than predicted by Crow's linear stability theory. The planar- standing-wave-angle value and the amplitude indifference behavior agree with Crow's analysis. Simulations with periodicity in the axial direction have indicated minor influence of axial flow on the Crow instability.

Zheng, Z. C.↗

Meteorology and Wake Vortex Influence on American Airlines FL-587 Accident

The atmospheric environment surrounding the crash of American Airlines Flight 587 is investigated. Examined are evidence for any unusual atmospheric conditions and the potential for encounters with aircraft wake vortices. Computer simulations are carried out with two different vortex prediction models and a Large Eddy Simulation model. Wind models are proposed for studying aircraft and pilot response to the wake vortex encounter.

Proctor, Fred H.↗

Exposing and Reducing Biases of Simulating Mixed-Phase Clouds in the Convection-Permitting E3SM Atmosphere Model: Lessons From an Arctic Cold-Air Outbreak

Mixed-phase clouds modulate the water and energy cycles of high-latitude regions, yet their liquid-ice phase partitioning has long been poorly simulated in climate models. Here, simulations of Arctic mixed-phase clouds by the Simple Cloud-Resolving E3SM Atmosphere Model (SCREAM) are assessed against large-eddy simulations, satellite data, and ground-based observations during the Cold-Air Outbreaks in the Marine Boundary Layer Experiment field campaign. SCREAM simulates nearly completely frozen clouds, which is attributed largely to the unreasonably strong Wegener–Bergeron–Findeisen (WBF) process that converts liquid to ice excessively and partly to the early over-abundant ice production at cold temperatures from a temperature-deterministic deposition ice nucleation scheme. Assuming no subgrid variation for the WBF process in the original formulation particularly conflicts with the instantaneous saturation adjustment assumption in the condensation scheme that assumes subgrid variability, leading to exaggerated WBF process rates. A proposed simple physically-based improvement on the treatment of subgrid cloud overlap substantially increases supercooled liquid water content and notably improves cloud-top phase partitioning, aligning better with observations. Improvement of supercooled liquid water content also converges with increasing horizontal resolution. The deposition ice nucleation scheme is found responsible for a falsely-produced ice cloud aloft that is not observed, biasing the simulated cloud radiative effects and top-of-atmosphere radiative fluxes. This study identifies key deficiencies in cloud parameterizations that continue to challenge convection-permitting models.

Geosciences↗

Exploring Unstructured Mesh Adaptation for Hybrid Reynolds-Averaged Navier–Stokes/Large Eddy Simulation

Mesh adaptation methods for the Reynolds-averaged Navier–Stokes (RANS) equations are rapidly maturing and beginning to impact the design of aerospace vehicles. RANS turbulence modeling improvements have slowed and may stagnate. Wall-modeled large eddy simulation (LES) and hybrid RANS/LES (HRLES) may provide an improved modeling capability but require specialized expertise to construct appropriate meshes and are considered be too ex-pensive for routine practical use. The realization of the CFD Vision 2030 Study includes improving geometry linkage, mesh generation/adaptation, and turbulence modeling/resolving methods for automated management of errors and uncertainties of physics-based, predictive modeling that can set the stage for ensuring a vehicle is in compliance with a regulation or specification (i.e., certification or qualification by analysis). An exploration of mesh adaptation for HRLES is performed to document synergies and challenges between mesh adaptation and HRLES. Vortex breakdown over a delta wing is examined to show the improvement of HRLES over RANS turbulence modeling approaches. A high lift configuration is shown to demonstrate complex geometry capability. Research and development opportunities are identified to advocate for continuing investments that may allow HRLES to enter routine practical use as a tool for aerospace vehicle analysis and design.

Michael A Park↗

Numerical Simulation of Wake Vortices Measured During the Idaho Falls and Memphis Field Programs

A numerical large-eddy simulation model is under modification and testing for application to aircraft wake vortices. The model, having a meteorological framework, permits the interaction of wake vortices with environments characterized by crosswind shear, stratification, and humidity. As part of the validation process, model results are compared with measured field data from the 1990 Idaho Falls and the 1994-1995 Memphis field experiments. Cases are selected that represent different aircraft and a cross section of meteorological environments. Also included is one case with wake vortex generation in ground effect. The model simulations are initialized with the appropriate meteorological conditions and a post roll-up vortex system. No ambient turbulence is assumed in our initial set of experiments, although turbulence can be self generated by the interaction of the model wakes with the ground and environment.

Proctor, Fred H.↗

A Model Intercomparison Study to Investigate Mixing Characteristics in Non‐Precipitating Stratocumulus Clouds

Recent aircraft observations of marine stratocumulus clouds consistently showed that cloud microphysical relationships vary with altitude, indicating inhomogeneous mixing characteristics near cloud top and homogeneous mixing characteristics in mid-levels of clouds. Here, we conduct model intercomparison of an idealized, non-precipitating stratocumulus cloud to evaluate model consistency and examine whether simulations can reproduce the observed mixing characteristics. The results show that eleven large-eddy simulations with various dynamics and microphysics schemes show good agreement on the thermodynamical, microphysical, and dynamical properties of the stratocumulus-topped boundary layer in a steady state. The inter-model spread in steady-state liquid water path is significantly reduced compared to previous model intercomparison studies. This improvement might be due to better models and more consistent initial conditions than those used decades ago. In addition, most simulations, including a low-dimensional simulation, capture inhomogeneous mixing characteristics near the cloud top and homogeneous mixing characteristics inside the cloud. Moreover, simulations using Lagrangian microphysics schemes agree better with the observed mixing characteristics compared with those using the bin microphysics schemes. Since most simulations do not fully resolve the entrainment process, the apparent mixing characteristics arise from the variations in the resolved cloud properties. Our results support the vertical circulation mixing hypothesis, which suggests that homogeneous mixing characteristics in mid-levels of clouds are due to the vertical circulation of entrainment-affected and diluted parcels from the cloud top moved to lower levels.

54 ENVIRONMENTAL SCIENCES↗

Numerical simulation of wall-bounded turbulent shear flows

Developments in three dimensional, time dependent numerical simulation of turbulent flows bounded by a wall are reviewed. Both direct and large eddy simulation techniques are considered within the same computational framework. The computational spatial grid requirements as dictated by the known structure of turbulent boundary layers are presented. The numerical methods currently in use are reviewed and some of the features of these algorithms, including spatial differencing and accuracy, time advancement, and data management are discussed. A selection of the results of the recent calculations of turbulent channel flow, including the effects of system rotation and transpiration on the flow are included.

Moin, P.↗