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.

SourceResearch connectionFrom yearThrough year

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

At least 1,369 records · Page 76

All sources NASA NTRS DOE OSTI

NASA NTRSDate not supplied

A Wall-Modeled LES Perspective for the High Lift Common Research Model Using LAVA

A new immersed boundary Wall-Modelled Large Eddy Simulation (WMLES) formulation is developed to study high-lift aerodynamics on the NASA High-Lift Common Research Model (HL-CRM). A sequence of Cartesian Octree grids with sizes ranging from 100 Million through 2.02 Billion grid points is utilized to systematically assess grid-sensitivity and convergence for the in-tunnel (QinetiQ) configuration of the model, and remarkable agreement between the immersed boundary and the curvilinear body-aligned WMLES formulations is reported on grids with comparable resolutions. In the free-air configuration of the model, consistent predictions between the Curvilinear Overset and the Cartesian Octree formulations are reported for angles of attack up to C(L,max) at a=19.57. However, some differences in the onset of stall are seen between the two methods for a>20°: while the curvilinear WMLES experiences wing-root separation with increasing angle of attack (Topology A), the Cartesian Octree formulation shows a different flow topology characterized by boundary layer weakness on the main element, emanating from the pylon-wing attachment (Topology B). In order to obtain further insight into the two-distinct topologies, carefully designed numerical experiments to isolate effects of the model standoff and the tunnel wall-boundary layers are conducted using the immersed boundary WMLES formulation. The increased incidence angle-of-attack on the inboard portion of the wing due to the standoff is shown to be sufficient for triggering a switch from Topology B to Topology A in Cartesian WMLES. The role of the floor boundary layer is further examined in detail by identification of additional corner-flow vorticity generated by the viscous juncture flow interactions between the floor boundary layer and the standoff leading to formation of a strong coherent and persistent vortex on the belly-side of the fuselage. The intensity of this vortex is shown to increase with the thickness of the floor boundary layer. A further increase in the incidence angle of attack near the leading-edge strake caused by the presence of this belly-side vortex is quantified for two-distinct floor boundary layers. Both of the floor boundary layers considered result in the onset of large scale wing-root separation at a=21.47in non-confined (free-air) configurations.

TTT↗

NASA NTRS2015

Investigation of Unsteady Tip Clearance Flow in a Low-Speed One and Half Stage Axial Compressor with LES And PIV

The primary focus of this paper is to investigate how a rotor's unsteady tip clearance flow structure changes in a low speed one and half stage axial compressor when the rotor tip gap size is increased from 0.5 mm (0.49% of rotor tip blade chord, 2% of blade span) to 2.4 mm (2.34% chord, 4% span) at the design condition are investigated. The changes in unsteady tip clearance flow with the 0.62 % tip gap as the flow rate is reduced to near stall condition are also investigated. A Large Eddy Simulation (LES) is applied to calculate the unsteady flow field at these three flow conditions. Detailed Stereoscopic PIV (SPIV) measurements of the current flow fields were also performed at the Johns Hopkins University in a refractive index-matched test facility which renders the compressor blades and casing optically transparent. With this setup, the unsteady velocity field in the entire flow domain, including the flow inside the tip gap, can be measured. Unsteady tip clearance flow fields from LES are compared with the PIV measurements and both LES and PIV results are used to study changes in tip clearance flow structures. The current study shows that the tip clearance vortex is not a single structure as traditionally perceived. The tip clearance vortex is formed by multiple interlaced vorticities. Therefore, the tip clearance vortex is inherently unsteady. The multiple interlaced vortices never roll up to form a single structure. When phased-averaged, the tip clearance vortex appears as a single structure. When flow rate is reduced with the same tip gap, the tip clearance vortex rolls further upstream and the tip clearance vortex moves further radially inward and away from the suction side of the blade. When the tip gap size is increased at the design flow condition, the overall tip clearance vortex becomes stronger and it stays closer to the blade suction side and the vortex core extends all the way to the exit of the blade passage. Measured and calculated unsteady flow fields inside the tip gap agree fairly well. Instantaneous velocity vectors inside the tip gap from both the PIV and LES do show flow separation and reattachment at the entrance of tip gap as some earlier studies suggested. This area at the entrance of tip gap flow (the pressure side of the blade) is confined very close to the rotor tip section. With a small tip gap (0.5mm), the gap flow looks like a simple two-dimensional channel flow with larger velocity near the casing for both flow rates. A small area with a sharp velocity gradient is observed just above the rotor tip. This strong shear layer is turned radially inward when it collides with the incoming flow and forms the core structure of the tip clearance vortex. When tip gap size is increased to 2.4 mm at the design operation, the radial profile of the tip gap flow changes drastically. With the large tip gap, the gap flow looks like a two-dimensional channel flow only near the casing. Near the rotor top section, a bigger region with very large shear and reversed flow is observed.

compressor↗

NASA NTRS2015

Investigation of Unsteady Tip Clearance Flow in a Low-Speed One and Half Stage Axial Compressor with LES And PIV

unsteady↗

DOE OSTI2016

Surface waves affect frontogenesis

Here, this paper provides a detailed analysis of momentum, angular momentum, vorticity, and energy budgets of a submesoscale front undergoing frontogenesis driven by an upper–ocean, submesoscale eddy field in a Large Eddy Simulation (LES). The LES solves the wave–averaged, or Craik–Leibovich, equations in order to account for the Stokes forces that result from interactions between nonbreaking surface waves and currents, and resolves both submesoscale eddies and boundary layer turbulence down to 4.9 m × 4.9 m × 1.25 m grid scales. It is found that submesoscale frontogenesis differs from traditional frontogenesis theory due to four effects: Stokes forces, momentum and kinetic energy transfer from submesoscale eddies to frontal secondary circulations, resolved turbulent stresses, and unbalanced torque. In the energy, momentum, angular momentum, and vorticity budgets for the frontal overturning circulation, the Stokes shear force is a leading–order contributor, typically either the second or third largest source of frontal overturning. These effects violate hydrostatic and thermal wind balances during submesoscale frontogenesis. The effect of the Stokes shear force becomes stronger with increasing alignment of the front and Stokes shear and with a nondimensional scaling. The Stokes shear force and momentum transfer from submesoscale eddies significantly energize the frontal secondary circulation along with the buoyancy.

58 GEOSCIENCES↗

DOE OSTI2016

Characterization of cumulus cloud fields using trajectories in the center of gravity versus water mass phase space: 1. Cloud tracking and phase space description

We study the evolution of warm convective cloud fields using large eddy simulations of continental and trade cumulus. Individual clouds are tracked a posteriori from formation to dissipation using a 3D cloud tracking algorithm and results are presented in the phase- space of center of gravity altitude versus cloud liquid water mass (CvM space). Here, the CvM space is shown to contain rich information on cloud field characteristics, cloud morphology, and common cloud development pathways, together facilitating a comprehensive understanding of the cloud field. In this part we show how the meteorological (thermodynamic) conditions that determine the cloud properties are projected on the CvM phase space and how changes in the initial conditions affect the clouds' trajectories in this space. This part sets the stage for a detailed microphysical analysis that will be shown in part II.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2016

Estimation of Cloud Fraction Profile in Shallow Convection Using a Scanning Cloud Radar

Large spatial heterogeneities in shallow convection result in uncertainties in estimations of domain-averaged cloud fraction profiles (CFP). This issue is addressed using large eddy simulations of shallow convection over land coupled with a radar simulator. Results indicate that zenith profiling observations are inadequate to provide reliable CFP estimates. Use of Scanning Cloud Radar (SCR), performing a sequence of cross-wind horizon-to-horizon scans, is not straightforward due to the strong dependence of radar sensitivity to target distance. An objective method for estimating domain-averaged CFP is proposed that uses observed statistics of SCR hydrometeor detection with height to estimate optimum sampling regions. This method shows good agreement with the model CFP. Results indicate that CFP estimates require more than 35 min of SCR scans to converge on the model domain average. Lastly, the proposed technique is expected to improve our ability to compare model output with cloud radar observations in shallow cumulus cloud conditions.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2017

Failure of Taylor's hypothesis in the atmospheric surface layer and its correction for eddy-covariance measurements

Taylors’ frozen turbulence hypothesis suggests that all turbulent eddies are advected by the mean streamwise velocity, without changes in their properties. This hypothesis has been widely invoked to compute Reynolds’ averaging using temporal turbulence data measured at a single point in space. However, in the atmospheric surface layer, the exact relationship between convection velocity and wavenumber k has not been fully revealed since previous observations were limited by either their spatial resolution or by the sampling length. Using Distributed Temperature Sensing (DTS), acquiring turbulent temperature fluctuations at high temporal and spatial frequencies, we computed convection velocities across wavenumbers using a phase spectrum method. We found that convection velocity decreases as k –1/3 at the higher wavenumbers of the inertial subrange instead of being independent of wavenumber as suggested by Taylor's hypothesis. We further corroborated this result using large eddy simulations. Applying Taylor's hypothesis thus systematically underestimates turbulent spectrum in the inertial subrange. As a result, a correction is proposed for point-based eddy-covariance measurements, which can improve surface energy budget closure and estimates of CO 2 fluxes.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2018

The Challenge of Identifying Controls on Cloud Properties and Precipitation Onset for Cumulus Congestus Sampled During MC3E

Here, the controls on precipitation onset and the transition from shallow cumulus to congestus are explored using a suite of 16 large–eddy simulations based on the 25 May 2011 event from the Midlatitude Continental Convective Clouds Experiment (MC3E). The thermodynamic variables in the model are relaxed at various timescales to observationally constrained temperature and moisture profiles in order to better reproduce the observed behavior of precipitation onset and total precipitation. Three of the simulations stand out as best matching the precipitation observations and also perform well for independent comparisons of cloud fraction, precipitation area fraction, and evolution of cloud top occurrence. All three simulations exhibit a destabilization over time, which leads to a transition to deeper clouds, but the evolution of traditional stability metrics by themselves is not able to explain differences in the simulations. Conditionally sampled cloud properties (in particular, mean cloud buoyancy), however, do elicit differences among the simulations. The inability of environmental profiles alone to discern subtle differences among the simulations and the usefulness of conditionally sampled model quantities argue for hybrid observational/modeling approaches. These combined approaches enable a more complete physical understanding of cloud systems by combining observational sampling of time–varying three–dimensional meteorological quantities and cloud properties, along with detailed representation of cloud microphysical and dynamical processes from numerical models.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Composited structure of non-precipitating shallow cumulus clouds

The normalized distributions of thermodynamic and dynamical variables both within and outside shallow clouds are investigated through a composite algorithm using large-eddy simulations of oceanic and continental cases. The normalized magnitude is maximum near the cloud centre and decreases outwards. While relative humidity (RH) and cloud liquid water ( q l ) decrease smoothly to match the environment, the vertical velocity, virtual potential temperature (θ v ), and potential temperature (θ) perturbations have more complicated behaviour towards the cloud boundary. Below the inversion layer, θ' v becomes negative before the vertical velocity has turned from an updraft to a subsiding shell outside the cloud, indicating the presence of a transition zone where the updraft is negatively buoyant. Due to the downdraft outside the cloud and enhanced horizontal turbulent mixing across the edge, the normalized turbulent kinetic energy (TKE) and horizontal turbulent kinetic energy (HTKE) decrease more slowly from the cloud centre outwards than the thermodynamic variables. The distributions all present asymmetric structures in response to the vertical wind shear, with more negatively buoyant air, stronger downdrafts, and larger TKE on the downshear side. We discuss several implications of the distributions for theoretical models and parameterizations. Positive buoyancy near the cloud base is mostly due to the virtual effect of water vapour, emphasizing the role of moisture in triggering. The mean vertical velocity is found to be approximately half the maximum vertical velocity within each cloud, providing a constraint to achieve possible power-law distributions for some models. Finally, the normalized distributions for different variables are used to estimate the vertical heat and moisture fluxes within clouds. The results suggest that distributions near the cloud edge and variability of maximum perturbations need careful treatment. The fluxes are underestimated in the inversion layer because cloud-top downdrafts cannot be captured well.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

A two-fluid single-column model of turbulent shallow convection, Part III: Results and parameter sensitivity

Two-fluid modelling has recently emerged as a promising approach to representing cumulus convection in weather and climate models. This study applies the two-fluid model described in Part II Thuburn et al., 2021b to a shallow cumulus convection case study over land (ARM). Large eddy simulation data is used to tune the majority of the closures that determine the properties of entrained and detrained air. The two-fluid model is generally able to reproduce the profiles of the mean and turbulent quantities over all stages of the diurnal cycle. As such, the initiation of shallow convection and the evolution of the cloud layer are well captured. The robustness of the two-fluid model is further verified using a steady-state test case (BOMEX), in which the cloud properties are also well modelled.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Numerical investigation of the influence of shear and thermal stratification on the wind turbine tip‐vortex stability

Summary The interaction between wind turbine wakes and atmospheric turbulence is characterised by complex dynamics. In this study, two major components of the atmospheric boundary layer dynamics have been isolated, namely, the mean velocity profile shear and the thermal stratification, to examine their impact on the near‐wake development by undertaking a series of highly resolved large‐eddy simulations. Subsequently, instantaneous flow fields are extracted from the simulations and used to conduct Fourier analysis and proper orthogonal decomposition (POD) and compute the mean kinetic energy fluxes by different POD modes to better understand the tip‐vortex instability mechanisms. Our findings indicate that shear can significantly affect the breakup of the wind turbine tip‐vortices and the shape and stable length of the wake, whereas thermal stratification seems to only have limited contribution to the spatial development of the near‐wake field. Finally, our analysis shows that the applied perturbation frequency determines the tip‐vortex breakup location as it controls the onset of the mutual inductance instability.

17 WIND ENERGY↗

DOE OSTI2023

A dynamic model of wind turbine yaw for active farm control

This paper presents a graph-based dynamic yaw model to predict the dynamic response of the hub-height velocities and the power of a wind farm to a change in yaw. The model builds on previous work where the turbines define the nodes of the graph and the edges represent the interactions between turbines. Advances associated with the dynamic yaw model include a novel analytical description of the deformation of wind turbine wakes under yaw to represent the velocity deficits and a more accurate representation of the interturbine travel time of wakes. The accuracy of the model is improved by coupling it with time- and space-dependent estimates of the wind farm inflow based on real-time data from the wind farm. The model is validated both statically and dynamically using large-eddy simulations. An application of the model is presented that incorporates the model into an optimal control loop to control the farm power output.

17 WIND ENERGY↗

DOE OSTI2025

Investigation of Main Bearing Fatigue Estimate Sensitivity to Synthetic Turbulence Models Using a Novel Drivetrain Model Implemented in OpenFAST

ABSTRACT A coupled medium‐fidelity drivetrain model is developed and implemented in OpenFAST for a 10‐MW land‐based reference turbine. The implementation is verified against a fully coupled multibody wind turbine model, including a detailed drivetrain. The new model can simultaneously and accurately estimate main bearing loads and represent elastic bending of the drivetrain. It has low computational cost and is useful for early design phases, sensitivity analyses and complex systems like wind farms (where computational expense must be expended elsewhere). Here, the model is implemented for a monopile offshore wind turbine and used to investigate the sensitivity of main bearing basic rating life to different synthetic turbulence models. Large‐eddy simulations (LES) targeting stable, neutral, and unstable atmospheric conditions at below‐, near‐ and above‐rated wind speeds are used as a reference. The turbulence models recommended by the International Electrotechnical Commission, the Mann spectral tensor model, and the Kaimal spectral model with exponential coherence are fitted to the LES data. Additionally, a constrained turbulence generator, PyConTurb (short for Python Constrained Turbulence ), based on LES data, is applied in the aero‐hydro‐servo‐elastic simulations. Taking PyConTurb as the baseline, the Kaimal model significantly underestimates fatigue of the downwind main bearing, with between 10% and 40% less damage. The Mann model also underestimates the downwind main bearing fatigue by up to 30%. The upwind main bearing damage is driven by mean loads, and differences between models are less significant, although the trends are similar. Reasons for these discrepancies are investigated and attributed to differences in spatial and temporal variations among the turbulence models.

17 WIND ENERGY↗

DOE OSTI2021

Reynolds stress tensor measurements using magnetic resonance velocimetry: expansion of the dynamic measurement range and analysis of systematic measurement errors

This study presents magnetic resonance velocimetry (MRV) Reynolds Stress measurements in a periodic hill channel with a hill Reynolds number of Re = 29,500. The velocity encoding scheme is based on the ICOSA6 method with six icosahedral encoding directions and multiple encoding values are measured to increase the dynamic range. The full Reynolds stress tensor is obtained from a voxel-wise three-dimensional Gaussian fit using the magnitude data of all acquisitions. The MRV results are compared to a wall-resolved large eddy simulation and laser Doppler velocimetry measurements conducted in the same channel. It is shown that the MRV Reynolds stress data have excellent precision and agree qualitatively with the reference data. However, there are apparent systematic deviations. One of the most prominent error contributions is the signal attenuation caused by higher orders of motion, which leads to an overestimation of the turbulence level. Another fundamental error is identified in the assumption that the turbulence is Gaussian distributed. With the presented reconstruction technique, the MRV data are fitted to a statistical model, and depending on the examined flow setup, the Gaussian model can lead to considerable errors. Possible ways of how to reduce all identified errors are presented. In summary, this technique enables Reynolds stress tensor measurements in complex internal flows with high dynamic range and excellent precision. However, several issues need to be resolved to make the turbulence quantification more accurate.

42 ENGINEERING↗

DOE OSTI2023

Parameterization and Explicit Modeling of Cloud Microphysics: Approaches, Challenges, and Future Directions

Cloud microphysical processes occur at the smallest end of scales among cloud-related processes and thus must be parameterized not only in large-scale global circulation models (GCMs) but also in various higher-resolution limited-area models such as cloud-resolving models (CRMs) and large-eddy simulation (LES) models. Instead of giving a comprehensive review of existing microphysical parameterizations that have been developed over the years, this study concentrates purposely on several topics that we believe are understudied but hold great potential for further advancing bulk microphysics parameterizations: multi-moment bulk microphysics parameterizations and the role of the spectral shape of hydrometeor size distributions; discrete vs “continuous” representation of hydrometeor types; turbulence-microphysics interactions including turbulent entrainment-mixing processes and stochastic condensation; theoretical foundations for the mathematical expressions used to describe hydrometeor size distributions and hydrometeor morphology; and approaches for developing bulk microphysics parameterizations. Also presented are the spectral bin scheme and particle-based scheme (especially, super-droplet method) for representing explicit microphysics. Their advantages and disadvantages are elucidated for constructing cloud models with detailed microphysics that are essential to developing processes understanding and bulk microphysics parameterizations. Particle-resolved direct numerical simulation (DNS) models are described as an emerging technique to investigate turbulence-microphysics interactions at the most fundamental level by tracking individual particles and resolving the smallest turbulent eddies in turbulent clouds. Outstanding challenges and future research directions are explored as well.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Main bearing response in a waked 15-MW floating wind turbine in below-rated conditions

Increased wind turbine size raises unknowns related to structural flexibility. Moreover, moving to deeper waters, component reliability becomes more critical. This work investigates main bearing response dependence on drivetrain flexibility and wake impingement in a two-turbine wind farm. A 15-MW floating direct-drive turbine is considered. Large eddy simulations (LES) are employed to model neutral, stable and unstable atmospheric conditions at below-rated mean wind speed, while the engineering codes OpenFAST and FAST.Farm simulate turbine and wake behavior. Results indicate significant sensitivities in fatigue estimates to lateral distance between the upstream and downstream turbine. The trends are most substantial in stable conditions, where the waked downwind main bearing sees twice the fatigue damage estimates of the upstream turbine for one position and 50% for another. Main bearing fatigue sensitivity to drivetrain flexibility is minor, while properly including generator rotor inertia loads is important for the axial forces of the locating (axially fixed) bearing, especially in stable conditions.

17 WIND ENERGY↗

DOE OSTI2021

Experimental study and demonstration of pilot-scale, dry feed, oxy-coal combustion under pressure

In this paper, we discuss how pressurized oxy-combustion is a promising technology for low-carbon, fossil fuel utilization. It has the potential of improved efficiency and economics compared with conventional atmospheric pressure oxy-combustion. Washington University in St. Louis has proposed a new pressurized oxy-combustion process, namely Staged Pressurized Oxy-Combustion, which has the potential to improve further the plant efficiency, operational flexibility, and economics. This process burns pulverized coal in a pressurized, oxy-combustion environment, which has not been demonstrated in a pilot-scale system before. To address this gap, a 100 kWth pressurized oxy-combustion facility was designed and constructed. This facility has a unique burner and furnace design featuring a co-axial, low-mixing flow field, which is drastically different from conventional coal-fired boiler designs where strong mixing is sought by introducing swirl or recirculating flows. This work aims to present the first pilot-scale experimental results of a dry-feed, pressurized oxy-combustion system. The tests focused on exploring flame stability and shape, char burnout, and fine particulate matter formation. Testing results suggest that the burner has excellent flame stability. The flame shape is consistent with the design philosophy and agrees with large eddy simulations. Importantly, complete char combustion can be achieved with an oxygen mole fraction in the flue gas of only 0.8%, as opposed to a required value of ~3% for conventional atmospheric pressure air-fired or oxyfuel combustion, which reduces the costs of both oxygen generation upstream and oxygen removal downstream. The testing results show promise for dry-feed, pressurized oxy-combustion, and the new burner and furnace design.

01 COAL, LIGNITE, AND PEAT↗

DOE OSTI2017

A minimally-dissipative low-Mach number solver for complex reacting flows in Open$\mathrm{FOAM}$

Large eddy simulation (LES) has become the de-facto computational tool for modeling complex reacting flows, especially in gas turbine applications. However, readily usable general-purpose LES codes for complex geometries are typically academic or proprietary/commercial in nature. In this report the objective is to develop and disseminate an open source LES tool for low-Mach number turbulent combustion using the OpenFOAM framework. In particular, a collocated-mesh approach suited for unstructured grid formulation is provided. Unlike other fluid dynamics models, LES accuracy is intricately linked to so-called primary and secondary conservation properties of the numerical discretization schemes. This implies that although the solver only evolves equations for mass, momentum, and energy, the implied discrete equation for kinetic energy (square of velocity) should be minimally-dissipative. Here, a specific spatial and temporal discretization is imposed such that this kinetic energy dissipation is minimized. The method is demonstrated using manufactured solutions approach on regular and skewed meshes for a canonical flow and a lab-scale turbulent flow problem.

42 ENGINEERING↗