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,783 records · Page 99

Impacts of Solar-Absorbing Aerosol Layers on the Transition of Stratocumulus to Trade Cumulus Clouds

The effects of an initially overlying layer of solar-absorbing aerosol on the transition of stratocumulus to trade cumulus clouds are examined using large-eddy simulations. For lightly drizzling cloud the transition is generally hastened, resulting mainly from increased cloud droplet number concentration (Nc) induced by entrained aerosol. The increased Nc slows sedimentation of cloud droplets and shortens their relaxation time for diffusional growth, both of which accelerate entrainment of overlying air and thereby stratocumulus breakup. However, the decrease in albedo from cloud breakup is more than offset by redistributing cloud water over a greater number of droplets, such that the diurnal-average shortwave forcing at the top of the atmosphere is negative. The negative radiative forcing is enhanced by sizable longwave contributions, which result from the greater cloud breakup and a reduced boundary layer height associated with aerosol heating. A perturbation of moisture instead of aerosol aloft leads to a greater liquid water path and a more gradual transition. Adding absorbing aerosol to that atmosphere results in substantial reductions in liquid water path (LWP) and cloud cover that lead to positive short-wave and negative longwave forcings on average canceling each other. Only for heavily drizzling clouds is the breakup delayed, as inhibition of precipitation overcomes cloud water loss from enhanced entrainment. Considering these simulations as an imperfect proxy for biomass burning plumes influencing Namibian stratocumulus, we expect regional indirect plus semi-direct forcings to be substantially negative to negligible at the top of the atmosphere, with its magnitude sensitive to background and perturbation properties.

Zhou, Xiaoli↗

Validation of a CFD Model Predicting the Effect of High Level Lateral Acceleration Sloshing on the Heat Transfer and Pressure Drop in a Small-Scale Tank in Normal Gravity

A two-phase CFD model is developed to study the effects of sloshing with high level lateral acceleration on the heat transfer and pressure drop in a small scale tank. Computational results are compared to the data provided by a non-isothermal sloshing experiment without phase change conducted by T. Himeno et al. at the University of Tokyo and JAXA in 2011. The results of the current model are, also, compared to CFD predictions reported by Himeno et al. A step change in lateral acceleration was applied in the experiment. Different levels of lateral acceleration amplitude, varying between 0G and 0.5G, were considered. CFD results for interface movement and tank pressure are presented and compared in this paper to the experimental data for the case in which the value of lateral acceleration was set to 0.5G. The effects of initial and boundary conditions and turbulence modeling approach on the tank pressure change during sloshing are discussed in detail. The effect of conjugate heat transfer in the tank wall is also studied to show its important role in determining the tank pressure evolution. The results of the Reynolds Averaged Navier Stokes (RANS) models are compared to the results of the Large Eddy Simulation model (LES) to underscore the importance of correctly capturing the effects of turbulence for high fidelity predictions.

multiphase flow↗

Impact of Wake Dispersion on Axial Compressor Performance

Rotor wake dispersion in a low-speed, one and half stage axial compressor is investigated in detail with a Large Eddy Simulation (LES). The primary focus is to quantify the total pressure recovery due to wake stretching and the total pressure loss from the rotor wake interaction with the stator blade boundary layer. The relative magnitude of the aerodynamic loss due to these two effects is examined at several radial locations. The spacing between the rotor and the stator was varied from 29% to 112% of the rotor axial chord at the mid span to investigate the effects of rotor wake decay before entering the stator passage. The current analysis indicates that the efficiency through the compressor stage is increased about 0.5% when the spacing between the rotor and the stator is decreased from 112% to 29% of the rotor axial chord at mid-span. 22% of the efficiency gain from the narrower axial gap is due to the wake recovery and 63% is due to the stronger unsteady pressure field at the exit of the rotor due to stage interaction. Total pressure loss/recovery across the stator varies significantly in the radial direction for the current compressor, which has a much lower aspect ratio. The total pressure recovery due to wake stretching is larger than the total pressure loss due to the unsteady boundary layer development on the stator blade from 20% to 35% of the span from the hub for 29% spacing and from 35% to 55% of the span for 112% spacing. Above 50% of the span, rotor tip clearance flow affects wake dispersion and the overall wake recovery is less than expected.

Axial Compressor↗

Validation of A CFD Model Predicting the Effect of High Level Lateral Acceleration Sloshing on the Heat Transfer and Pressure Drop in A Small-Scale Tank in Normal Gravity

A two-phase CFD model is developed to study the effects of sloshing with high level lateral acceleration on the heat transfer and pressure drop in a small scale tank. Computational results are compared to the data provided by a non-isothermal sloshing experiment without phase change conducted by T. Himeno et al. at the University of Tokyo and JAXA in 2011 [1]. The results of the current model are, also, compared to CFD predictions reported by Himeno et al. [2]. A step change in lateral acceleration was applied in the experiment. Different levels of lateral acceleration amplitude, varying between 0G and 0.5G, were considered. CFD results for interface movement and tank pressure are presented and compared in this paper to the experimental data for the case in which the value of lateral acceleration was set to 0.5G. The effects of initial and boundary conditions and turbulence modeling approach on the tank pressure change during sloshing are discussed in detail. The effect of conjugate heat transfer in the tank wall is also studied to show its important role in determining the tank pressure evolution. The results of the Reynolds Averaged Navier Stokes (RANS) models are compared to the results of the Large Eddy Simulation model (LES) to underscore the importance of correctly capturing the effects of turbulence for high fidelity predictions.

multiphase flow↗

Development of Dynamic Sub-Grid Models for Variational Multiscale Methods

A dynamic Variational Multiscale Method (Hughes et al. 1998) is developed by leveraging the Germano procedure from classical Large-eddy Simulations (LES). The similarity between the classical and variational approaches is analyzed in the context of incompressible flow. This analysis leads to a consistent modeling approach for both incompressible and compressible flows, the latter being demonstrated in a priori testing for low-speed attached and separated boundary layers. Similar to the classical LES procedure from which it is derived, the variational dynamic procedure does not guarantee a positive semi-definite coefficient in the general case. However, reproducing the behavior of the classical LES dynamic approach is seen as a necessary first step to develop a VMM that automatically adjusts to the local resolution and flow physics.

Variational multiscale↗

Sidewall Effects on Exact Reynolds-Stress Budgets in an Impinging Shock Wave/Boundary Layer Interaction

Large-eddy simulations are performed using wall-resolved mesh for a Mach 2.29 impinging shock wave/boundary-layer interaction. Flow conditions are based on an experiment and therefore entire span was simulated, including the two sidewalls. Mean flow comparison with the experimental data showed that the interaction was larger in the simulation. Time-series analysis of a rake of pressure probes immediately downstream of the mean reflected shock position showed a peak in weighted power spectral density occurred about $St_{Lint}=0.01$, owing to a larger interaction length. Budgets of Reynolds-stress transport calculated across the span and along the corner bisector showed high degree of anisotropy. Merging of the secondary flows and separation along the corner gives rise to unstablecounter rotating vortices, which straddle the corner and grow in size. This also leads to a development of new behavior in the viscous sublayer along the corner bisector, where the pressure strain andmolecular diffusion mechanisms become prominent.

Vyas, Manan A.↗

Application of Lattice Boltzmann and Navier-Stokes Methods to NASA's Wall Mounted Hump

Lattice Boltzmann (LB) based Large Eddy Simulation (LES), Reynolds-averaged Navier-Stokes (RANS) as well as hybrid RANS/LES methods within the Launch Ascent and Vehicle Aerodynamics (LAVA) solver framework are applied to NASA's wall-mounted hump. Computational results are compared with experiments performed by Greenblatt et al. A detailed comparison between the accuracy and resolution requirements of the two approaches for turbulence resolving simulations, as well as the suitability of different grid paradigms (body-fitted curvilinear and block structured Cartesian) are presented. This test case is part of NASA's Revolutionary Computational Aerosciences (RCA) sub-project which addresses the technical challenge of predicting flow separation and reattachment accurately. Improvements in predictive accuracy by as much as 90% are demonstrated using LB as well as hybrid RANS/LES approaches compared to state-of-the-art steady state RANS simulations.

Lattice Boltzmann↗

Bowles-Tatnall Wake Vortex Encounter Hazard Metric

Wake vortex spacing standards constrict the terminal area throughput and impose severe constraints on the overall capacity and efficiency of the National Airspace System. For more than two decades starting in the early 1990s, the National Aeronautics and Space Administration conducted extensive research on characterizing the formation and evolution of aircraft wakes. This multidisciplinary work included comprehensive field experiments (Pruis et al. 2016), flight tests (Vicroy et al. 1998), and wind tunnel tests (Rossow 1994; Chow et al. 1997). Parametric studies using large eddy simulations (Proctor 1998; Proctor et al. 2006) were conducted in order to develop fast-time models for the prediction of wake transport and decay (Ahmad et al. 2016). Substantial effort was spent on the formulation of acceptable vortex hazard metrics (Tatnall 1995; Hinton and Tatnall 1997). Several wake encounter severity metrics have been suggested in the past, which include the wake circulation strength, vortex-induced rolling moment coefficient (Clv), bank angle, and the roll control ratio (Tatnall 1995; Hinton and Tatnall 1997; Van der Geest 2012). The vortex-induced rolling moment coefficient introduced by Bowles and Tatnall (Tatnall 1995; Gloudemans et al. 2016) has been used extensively for risk and safety analysis of newly proposed air traffic management concepts and procedures. The original method of Bowles and Tatnall assumed a constant wing loading (the wing lift-curve slope, CL is constant), which resulted in an overestimation of the vortexinduced rolling moment coefficient. Bowles (2014) suggested a correction to the original method that provides more accurate values of Clv and which is also consistent with the underlying physics of the problem. The overestimation of Clv in the original method can be corrected by assuming an elliptical lift distribution. Figure 1.1 illustrates the correction in Clv achieved by the modified method.

Joel Malissa↗

The Cumulus and Stratocumulus CloudSat-CALIPSO Dataset (CASCCAD)

Low clouds continue to contribute greatly to the uncertainty in cloud feedback estimates. Depending on whether a region is dominated by cumulus (Cu) or stratocumulus (Sc) clouds, the interannual low-cloud feedback is somewhat different in both spaceborne and large-eddy simulation studies. Therefore, simulating the correct amount and variation of the Cu and Sc cloud distributions could be crucial to predict future cloud feedbacks. Here we document spatial distributions and profiles of Sc and Cu clouds derived from Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observations (CALIPSO) and CloudSat measurements. For this purpose, we create a new dataset called the Cumulus And Stratocumulus CloudSat-CALIPSO Dataset (CASCCAD), which identifies Sc, broken Sc, Cu under Sc, Cu with stratiform outflow and Cu. To separate the Cu from Sc, we design an original method based on the cloud height, horizontal extent, vertical variability and horizontal continuity, which is separately applied to both CALIPSO and combined CloudSat–CALIPSO observations. First, the choice of parameters used in the discrimination algorithm is investigated and validated in selected Cu, Sc and Sc–Cu transition case studies. Then, the global statistics are compared against those from existing passive- and active-sensor satellite observations. Our results indicate that the cloud optical thickness – as used in passive-sensor observations – is not a sufficient parameter to discriminate Cu from Sc clouds, in agreement with previous literature. Using clustering-derived datasets shows better results although one cannot completely separate cloud types with such an approach. On the contrary, classifying Cu and Sc clouds and the transition between them based on their geometrical shape and spatial heterogeneity leads to spatial distributions consistent with prior knowledge of these clouds, from ground-based, ship-based and field campaigns. Furthermore, we show that our method improves existing Sc–Cu classifications by using additional information on cloud height and vertical cloud fraction variation. Finally, the CASCCAD datasets provide a basis to evaluate shallow convection and stratocumulus clouds on a global scale in climate models and potentially improve our understanding of low-level cloud feedbacks. The CASCCAD dataset (Cesana, 2019, https://doi.org/10.5281/zenodo.2667637) is available on the Goddard Institute for Space Studies (GISS) website at https://data.giss.nasa.gov/clouds/casccad/ (last access: 5 November 2019) and on the zenodo website at https://zenodo.org/record/2667637 (last access: 5 November 2019).

Cesana, Gregory V.↗

Adaptive Embedded LES of the NASA Hump

A scheme for adaptive embedded LES is proposed which automatically determines boundaries for LES regions in a hybrid LES-RANS computation, with the goal of minimizing the LES part of the computation for maximum accuracy with minimum cost. The model-invariant hybrid formulation enables this scheme through greater flexibility in the placement of RANS-LES transitions. An adaptive embedded large-eddy simulation is carried out for the NASA hump test case and adaptive meshing is added to show how additional adaptive features may be controlled by the adaptive hybrid scheme.

Woodruff, Stephen↗

Wall-modeled LES of the NASA Juncture Flow Experiment

We perform wall-modeled large eddy simulation (WMLES) of the NASA Juncture Flow experiment (Kegerise & Neuhart, NASA/TM–2019–20286) using an equilibrium wall model and unstructured finite volume solver, to assess its predictive ability for complex geometries. The flow condition simulated corresponds to 5 degrees angle of attack with a freestream Mach number of 0.189, and a Reynolds number based on the mean crank chord of 2.4 million. All simulations neglect the effects of wind tunnel walls, sting and mast present in the experiment. To assess the performance of WMLES on simpler flows with very coarse grids, we simulate flow in a turbulent channel at a friction Reynolds number, Reτ≈2000, and flow through a square duct at anReτ≈1000. The duct flow results using the equilibrium wall model indicate that the skin friction is inaccurate near the corner for grid topologies and resolutions typically used in WMLES, having implications for more complicated flows such as the juncture flow. For the juncture flow simulation, two different approaches are investigated. The first approach consists of a truncated-domain simulation wherein the inflow plane of the computational domain is placed at≈0.2chord length of the wing, and we prescribe the mean flow from a separate Reynolds-averaged Navier-Stokes (RANS) solution along with synthetic turbulence to initiate realistic unsteadiness in the domain. The second approach involves simulating the entire geometry with trips to trigger transition to turbulence. The truncated- and full-domain simulations contain about 62 and 90 million cells, respectively; with 8-10 points per boundary layer thickness in the wall-parallel directions, 14-20 points in the wall-normal direction and a near-wall viscous spacing (∆n+1) of≈100. Preliminary results are encouraging overall in terms of the prediction of wall pressure, wall skin friction, velocity and stresses; but indicate that further work is required to improve the predictions in the separation bubble and wing-fuselage corner regions. The grid resolution used in this study is still fairly coarse, and the results should be interpreted as work in progress.

Prahladh S Iyer↗

TPSAS-NF1676L-35395-DND

Reynolds-averaged and hybrid Reynolds-averaged / large eddy simulations are performed for a flush-wall hypervelocity fuel injector that has been studied experimentally as part of the Enhanced Injection and Mixing Project at NASA Langley.

Robert A Baurle↗

Physics Informed Broadband Noise Source Identification and Prediction of an Ideally Twisted Rotor

This work aims at a physics based approach for the prediction and isolation of broadband noise emanating from different regions of an ideally twisted rotor. A preliminary prediction was conducted using a lattice-Boltzmann method–very-large-eddy simulation (LBM-VLES) implemented within the software suite, PowerFLOW. Regions of particular interest to broadband noise were investigated by calculating one-third octave sound pressure levels of the unsteady pressure fluctuations acting on the rotor. These regions were then isolated and treated as individual Ffowcs Williams-Hawkings (FW-H) surfaces for three run conditions at a much finer spatial resolution to delineate the broadband noise from these separate regions from the total acoustic spectra. These predictions were then compared to experimentally acquired data as well as to semi-empirical prediction methods to highlight the acoustic contributions of various broadband noise generation mechanisms as well as to exemplify and explain shortcomings in the semi-empirical methodology.

Christopher S Thurman↗

Computational Prediction of Broadband Noise from a Representative Unmanned Aerial System Rotor

This work details the application of a lattice-Boltzmann method–very-large-eddy simulation (LBM-VLES) employed by the software suite, PowerFLOW. This LBM-VLES simulation predicted the aeroacoustic noise emanating from a representative, small unmanned aircraft system rotor, namely, the DJI-9450 in a hover condition. Predicted total aerodynamic loading as well as 2D aerodynamic loading along discrete spanwise sections of a rotor blade were compared to lower fidelity predictions and experimental results acquired in the Structural Acoustic Loads and Transmission anechoic chamber facility at the NASA Langley Research Center. The total acoustic spectra were decomposed into tonal and broadband components, which showed that broadband noise was a dominant contributor above 1 kHz for this rotor. These data were then compared to experimentally acquired data, showing good agreement up to approximately 11 kHz. Above 11 kHz, however, a grid sensitivity study showed dependency of the highest resolvable frequency on the spatial resolution of the computational domain, explaining the roll off in predicted data. Individual broadband noise sources were further investigated by calculating one-third octave sound pressure levels of the unsteady pressure fluctuations acting on the rotor, providing evidence that blade self-noise was the prominent noise source. Using these results, blade wake interaction noise was seen to be negligible for this particular rotor, which was further validated by calculating blade vortex miss distances and comparing to theory.

Christopher S. Thurman↗

Physics-Informed Broadband Noise Source Identification and Prediction of an Ideally Twisted Rotor

This work aims to provide a physics-based approach for the prediction and isolation of broadband noise emanating from different regions of an ideally twisted rotor. A preliminary prediction was conducted using a lattice-Boltzmann method–very-large-eddy simulation (LBM-VLES) implemented within the software suite, PowerFLOW. Regions of particular interest to broadband noise were investigated by calculating one-third octave sound pressure levels of the unsteady pressure fluctuations acting on the rotor. These regions were then treated as individual Ffowcs Williams and Hawkings (FW-H) surfaces and simulated for three run conditions at a finer spatial resolution to identify the broadband noise from these separate regions and isolate it from the total acoustic spectra. These predictions were then compared to experimentally acquired data and to semiempirical prediction methods to highlight the acoustic contributions of various broadband noise generation mechanisms as well as to exemplify and explain shortcomings in the semiempirical methodology.

Christopher S Thurman↗

Multiphysics Computational Analysis of a Perforated Plate Cooling Flow

A multiphysics analysis of a film cooling flow was conducted using a computational fluid dynamics (CFD) solver in conjunction with a conjugate heat transfer (CHT) solver. The film cooling problem consisted of a square nozzle producing a heated subsonic jet over a perforated plate with three sections of 45 cooling holes each. Results from the multiphysics analysis were compared to a CFD-only approach and experimental results. Reynolds-averaged Navier-Stokes (RANS) computations were employed first, with and without the CHT solver. The multiphysics approach provided much closer agreement with experimentally measured plate temperatures than the CFD-only approach. Comparisons of computed velocities and temperatures were also made with experimental data. Hybrid RANS-large eddy simulation calculations employing the Delayed Detached Eddy Simulation (DDES) approach were also conducted within the multiphysics framework. The DDES approach was unable to resolve turbulent structures in the near wall regions. Hence, the RANS part of the DDES model was used in near wall regions, and produced surface temperatures along with velocity and temperature profiles close to the wall that were very similar to the RANS results. The DDES approach did resolve structures in the outer portion of the boundary layer especially near the end of the plate where the boundary layer was thicker. The key result was the improved surface temperatures produced by the multiphysics approach employing CHT over the fluid-only approach.

Computational Fluid Dynamics↗

Improving the Representation of Subtropical Boundary Layer Clouds in the NASA GEOS Model with the Eddy-Diffusivity/Mass-Flux Parameterization

A systematic underestimation of subtropical planetary boundary layer (PBL) stratocu-mulus clouds by the GEOS model has been significantly improved by a new Eddy-Diffusivity/Mass-Flux (EDMF) parameterization. The EDMF parameterization representsthe subgrid-scale transport in the dry and moist parts of the PBL in a unified mannerand it combines an adjusted eddy-diffusivity PBL scheme from GEOS with a stochasticmulti-plume mass-flux model. The new EDMF version of the GEOS model is first com-pared against the CONTROL version in a single-column model (SCM) framework for twobenchmark cases representing subtropical stratocumulus and shallow cumulus clouds, andvalidated against large eddy simulations. Global simulations are performed and comparedagainst observations and reanalysis data. The results show that the EDMF version of theGEOS model produces more realistic subtropical PBL clouds. Noteworthy, the EDMF im-provements first detected in the SCM framework translate into similar improvements of theglobal GEOS model.

Kay Suselj↗

Turbulence, Transition, and Numerical Method Technologies

The close collaboration in the validation experimentation effort was excellent, the complementary research efforts related tofundamental numerical model development, and the activities related to effective HPC utilization on near term architectures that are coming down the pipeline were exactly what NASA should be doing. Even though the TQR panel cannot say that the success criteria were fully met, we do feel that the research funded by TTT/RCA was successful and useful in determining the next suite of CFD validation tests and the problematic flow conditions that continue to be difficult to predict using our current suite of computational tools. Given the TQR Panel’s review criteria: 1. The deliverables meet the technical requirements; 2.The appropriate technical approaches were followed in producing the deliverables; and 3. All technical caveats and concerns have been identified and the deliverables: a.Development of more accurate physics-based methods (e.g., higher moment closure); b. Large eddy simulation (LES); c. Advanced numerical methods; d. Transition prediction and modeling; e.Validation experiments; f. Multidisciplinary analysis and design (high fidelity). TQR panel feels that the technical processes were all completed and that the review criteria and all the deliverables were met. Thus, the Milestone TACP-2016-001 has been completed.

Mujeeb R Malik↗