SEARCH · Engineering Papers

Results for “Average model”

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 721 records · Page 40

All sources NASA NTRS DOE OSTI

DOE OSTI2022

Towards string order melting of spin-1 particle chains in superconducting transmons using optimal control

Utilizing optimal control to simulate a model Hamiltonian is an emerging strategy that leverages the intrinsic physics of a device with digital quantum simulation methods. Here we evaluate optimal control for probing the nonequilibrium properties of symmetry-protected topological (SPT) states simulated with superconducting hardware. Assuming a tunable transmon architecture, we cast the evolution of these SPT states as a series of one- and two-site pulse optimization problems that are solved in the presence of leakage constraints. From the generated pulses, we classically simulate the time-dependent melting of the perturbed SPT string order across a six-site model with an average state infidelity of 10 -3 . The feasibility of these pulses as well as their efficient application indicate that high-fidelity simulations of string order melting are within reach of current quantum computing systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2022

Probing the Diffuse Ly α Emission on Cosmological Scales: Ly α Emission Intensity Mapping Using the Complete SDSS-IV eBOSS

Based on Sloan Digital Sky Survey Data Release 16, we have detected the large-scale structure of Lyα emission in the universe at redshifts z = 2–3.5 by cross-correlating quasar positions and Lyα emission imprinted in the residual spectra of luminous red galaxies. We apply an analytical model to fit the corresponding Lyα surface brightness profile and multipoles of the redshift-space quasar–Lyα emission cross-correlation function. The model suggests an average cosmic Lyα luminosity density of $6.6^{+3.3}_{-3.1}$ x 10 40 erg s -1 cMpc -3 , a ~2σ detection with a median value about 8–9 times those estimated from deep narrowband surveys of Lyα emitters at similar redshifts. Although the low signal-to-noise ratio prevents us from a significant detection of the Lyα forest–Lyα emission cross-correlation, the measurement is consistent with the prediction of our best-fit model from quasar–Lyα emission cross-correlation within current uncertainties. We rule out the scenario where the Lyα photons mainly originate from quasars. We find that Lyα emission from star-forming galaxies, including contributions from that concentrated around the galaxy centers and that in diffuse Lyα-emitting halos, is able to explain the bulk of the Lyα luminosity density inferred from our measurements. Ongoing and future surveys can further improve the measurements and advance our understanding of the cosmic Lyα emission field.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS1992

Turbomachinery CFD on parallel computers

The role of multistage turbomachinery simulation in the development of propulsion system models is discussed. Particularly, the need for simulations with higher fidelity and faster turnaround time is highlighted. It is shown how such fast simulations can be used in engineering-oriented environments. The use of parallel processing to achieve the required turnaround times is discussed. Current work by several researchers in this area is summarized. Parallel turbomachinery CFD research at the NASA Lewis Research Center is then highlighted. These efforts are focused on implementing the average-passage turbomachinery model on MIMD, distributed memory parallel computers. Performance results are given for inviscid, single blade row and viscous, multistage applications on several parallel computers, including networked workstations.

Blech, Richard A.↗

NASA NTRS1992

Turbomachinery CFD on parallel computers

Blech, R. A.↗

DOE OSTI2023

Plasma Properties in the Earth's Magnetosheath Near the Subsolar Magnetopause: Implications for Geocoronal Density Estimates

Combined in situ ion measurements and remote sensing of energetic neutral atoms are used to determine the geocoronal Hydrogen density at large (∼10 R E ) distances from the Earth. This method for determining the geocoronal density requires global magnetospheric modeling. Observations in the Earth's subsolar magnetosheath from the Magnetospheric Multiscale mission are used to determine the accuracy of using global models to predict the geocoronal density. On average, gas dynamic and magnetohydrodynamic (MHD) models and observations are in reasonable agreement, with differences <25%. In addition, the MHD model subsolar magnetopause is about 0.5 R E sunward of the observed location. However, variations around averages are large (up to a factor of 2), indicating that global models introduce relatively large uncertainties in geocoronal density estimates. Finally, the critical ion flux in the Interstellar Boundary Explorer IBEX‐Hi energy range is often minimally affected by fluctuations of a factor of 2 in the density.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS1984

A Crustal Structure Study of South America

Variations in the average crustal thickness of the South American platform were determined by measuring fundamental mode Rayleigh wave dispersion between nine pairs of South American WWSSN seismograph stations and inverting the measurements to obtain average shear wave velocity-depth models across the platform. Additional models were derived from the Rayleigh wave dispersion data previously measured over eastern South America. Both sets of models were interpreted in terms of average crustal thickness and average crustal shear wave velocity. The results obtained were depicted in a contour map of crustal thickness. Average values for seismic velocities in the crust and upper mantle were calculated and compared with maps of regional gravity models, MAGSAT magnetic anomalies, and major tectonic features. The correlations and relationships discovered are discussed.

Renbarger, K. S.↗

DOE OSTI2024

ExaWind: Open‐source CFD for hybrid‐RANS/LES geometry‐resolved wind turbine simulations in atmospheric flows

Abstract Predictive high‐fidelity modeling of wind turbines with computational fluid dynamics, wherein turbine geometry is resolved in an atmospheric boundary layer, is important to understanding complex flow accounting for design strategies and operational phenomena such as blade erosion, pitch‐control, stall/vortex‐induced vibrations, and aftermarket add‐ons. The biggest challenge with high‐fidelity modeling is the realization of numerical algorithms that can capture the relevant physics in detail through effective use of high‐performance computing. For modern supercomputers, that means relying on GPUs for acceleration. In this paper, we present ExaWind, a GPU‐enabled open‐source incompressible‐flow hybrid‐computational fluid dynamics framework, comprising the near‐body unstructured grid solver Nalu‐Wind, and the off‐body block‐structured‐grid solver AMR‐Wind, which are coupled using the Topology Independent Overset Grid Assembler. Turbine simulations employ either a pure Reynolds‐averaged Navier–Stokes turbulence model or hybrid turbulence modeling wherein Reynolds‐averaged Navier–Stokes is used for near‐body flow and large eddy simulation is used for off‐body flow. Being two‐way coupled through overset grids, the two solvers enable simulation of flows across a huge range of length scales, for example, 10 orders of magnitude going from O(μm) boundary layers along the blades to O(10 km) across a wind farm. In this paper, we describe the numerical algorithms for geometry‐resolved turbine simulations in atmospheric boundary layers using ExaWind. We present verification studies using canonical flow problems. Validation studies are presented using megawatt‐scale turbines established in literature. Additionally presented are demonstration simulations of a small wind farm under atmospheric inflow with different stability states.

17 WIND ENERGY↗

NASA NTRS2006

The Effect of Cumulus Cloud Field Anisotropy on Domain-Averaged Solar Fluxes and Atmospheric Heating Rates

Cumulus clouds can become tilted or elongated in the presence of wind shear. Nevertheless, most studies of the interaction of cumulus clouds and radiation have assumed these clouds to be isotropic. This paper describes an investigation of the effect of fair-weather cumulus cloud field anisotropy on domain-averaged solar fluxes and atmospheric heating rate profiles. A stochastic field generation algorithm was used to produce twenty three-dimensional liquid water content fields based on the statistical properties of cloud scenes from a large eddy simulation. Progressively greater degrees of x-z plane tilting and horizontal stretching were imposed on each of these scenes, so that an ensemble of scenes was produced for each level of distortion. The resulting scenes were used as input to a three-dimensional Monte Carlo radiative transfer model. Domain-average transmission, reflection, and absorption of broadband solar radiation were computed for each scene along with the average heating rate profile. Both tilt and horizontal stretching were found to significantly affect calculated fluxes, with the amount and sign of flux differences depending strongly on sun position relative to cloud distortion geometry. The mechanisms by which anisotropy interacts with solar fluxes were investigated by comparisons to independent pixel approximation and tilted independent pixel approximation computations for the same scenes. Cumulus anisotropy was found to most strongly impact solar radiative transfer by changing the effective cloud fraction, i.e., the cloud fraction when the field is projected on a surface perpendicular to the direction of the incident solar beam.

Hinkelman, Laura M.↗

DOE OSTI2021

In-Situ Calibrated Modeling of Residual Stresses Induced in Machining under Various Cooling and Lubricating Environments

Although many functional characteristics, such as fatigue life and damage resistance depend on residual stresses, there are currently no industrially viable ‘Digital Process Twin’ models (DPTs) capable of efficiently and quickly predicting machining-induced stresses. By leveraging advances in ultra-high-speed in-situ experimental characterization of machining and finishing processes under plane strain (orthogonal/2D) conditions, we have developed a set of physics-based semi-analytical models to predict residual stress evolution in light of the extreme gradients of stress, strain and temperature, which are unique to these thermo-mechanical processes. Initial validation trials of this novel paradigm were carried out in Ti-6Al4V and AISI 4340 alloy steel. A variety dry, cryogenically cooled and oil lubricated conditions were evaluated to determine the model’s ability to capture the tribological changes induced due to lubrication and cooling. The preliminarily calibrated and validated model exhibited an average correlation of better than 20% between the predicted stresses and experimental data, with calculation times of less than a second. Based on such fast-acting DPTs, the authors envision future capabilities in pro-active surface engineering of advanced structural components (e.g., turbine blades).

36 MATERIALS SCIENCE↗

DOE OSTI2026

A Multi-Region SEIR Model Incorporating Inter-County Mobility and Time-Dependent Transmission Dynamics: Application to COVID-19 Disease Outbreak Data in North Carolina.

Classical infectious disease compartmental models typically do not incorporate spatial heterogeneity or mobility. We develop a multi-region susceptible-exposed-infected-recovered (SEIR) model in which disease dynamics are coupled via inter-region mobility and the transmission rate is both region and time dependent. We calibrate the model using rolling averages of daily COVID-19 data in all 100 North Carolina counties. Mobility parameters are prescribed using daily inter-county commuter data. The number of transmission rate parameters is substantially reduced by hypothesizing that the dynamics correlate with county-level population density. Parameter estimation is carried out using several objective functions with error terms at different scales. An additive combination of least squares error at the county-level and the state-level, along with a quadratic transmission rate polynomial, yields the lowest overall error at both spatial scales. The calibrated model is used to simulate regional effects of perturbing disease transmission rates in adjacent counties and to illustrate effects of the state’s mobility infrastructure on disease dynamics and spread for a new disease outbreak.

COVID-19 modeling↗

DOE OSTI2021

Floating platform effects on power generation in spar and semisubmersible wind turbines

The design and financing of commercial-scale floating offshore wind projects require a better understanding of how power generation differs between newer floating turbines and well-established fixed-bottom turbines. In floating turbines, platform mobility causes additional rotor motion that can change the time-averaged power generation. In this work, OpenFAST simulations examine the power generated by the National Renewable Energy Laboratory's 5-MW reference turbine mounted on the OC3-UMaine spar and OC4-DeepCWind semisubmersible floating platforms, subjected to extreme irregular waves and below-rated turbulent inflow wind from large-eddy simulations of a neutral atmospheric boundary layer. For these below-rated conditions, average power generation in floating turbines is most affected by two types of turbine displacements: an average rotor pitch angle that reduces power, caused by platform pitch; and rotor motion upwind-downwind that increases power, caused by platform surge and pitch. The relative balance between these two effects determines whether a floating platform causes power gains or losses compared to a fixed-bottom turbine; for example, the spar creates modest (3.1%–4.5%) power gains, whereas the semisubmersible creates insignificant (0.1%–0.2%) power gains for the simulated conditions. Furthermore, platform surge and pitch motions must be analyzed concurrently to fully capture power generation in floating turbines, which is not yet universal practice. Finally, a simple analytical model for predicting average power in floating turbines under below-rated wind speeds is proposed, incorporating effects from both the time-averaged pitch displacement and the dynamic upwind-downwind displacements.

17 WIND ENERGY↗

DOE OSTI2025

Revealing Defect-Seeded and Interfacial Generation Mechanisms of Photoinduced Coherent Phonons with 4D Ultrafast Electron Microscopy

Ultrafast photoexcitation of coherent phonons is driven by an impulsive, collective displacement of constituent atoms from their average equilibrium lattice positions [1]. Models describing the generation of coherent acoustic modes typically invoke the creation of an anisotropic strain profile arising from the relatively instantaneous absorption of an ultrafast laser pulse [2]. If the skin depth is shallow relative to the specimen thickness, a steep tensile strain gradient perpendicular to the surface ($\frac{∂ε}{∂z}$) results. Initial relaxation occurs via rapid contraction of the surface layers followed by subsequent coherent oscillations of the lattice and launch of a train of coherent elastic strain waves [i.e., coherent acoustic phonons (CAPs)]. Macroscopically, responses are generally well-described by constitutive relations as gleaned from data gathered using ultrasonic methods or ultrafast spectroscopies. Furthermore, at the atomic to nanoscale level, individual lattice discontinuities and their impact on CAP behaviors can be modeled using multiscale methods [3,4]. Further, average unit-cell level responses on ultrafast timescales can be probed using femtosecond electron and X-ray scattering [5,6].

Flannigan, David J. [University of Minnesota, Minn↗

DOE OSTI2022

Impact of particle size and particle-flow-wall coupling on pressurized oxy-combustion in the down-fired burner

Concerns over climate change have led to numerous efforts in developing low-carbon energy technologies. Pressurized oxy-combustion (POC) is a promising candidate to reduce carbon emissions in power generation. Designing an effective burner plays a vital role in developing new coal combustion technologies. Because of the high pressure, the volume fraction of coal particles at the fuel inlet of the burner of a pressurized oxy-combustor is close or even higher than the maximum limit that commercial CFD codes (e.g., ANSYS FLUENT) can handle. At this high particle volume fraction, the interactions among particles, fluid flow, and wall need to be re-evaluated. The present computational work is the first step of a systematic analysis of the particle influence, like releasing method, releasing locating, and particle size, in a pilot-scale POC combustor, developed at Washington University in St. Louis (WUSTL). In a POC process, pulverized coal is burned under elevated pressure and O2-CO2 environment. Specifically, a 15-bar, 100 kWth, POC combustor is modeled employing ANSYS FLUENT, using Reynolds-averaged Navier-Stokes (RANS) modeling. It is revealed that for this pilot-scale, pressurized burner, velocity profiles of the near-wall region in this POC facility exhibit some discrepancy against the near-wall turbulent flow velocity profile. In order to investigate the particle influence in the near-wall region, particle releasing location will be investigated. Numerical simulation results exhibit the coupling effect of turbulence flow and particles in this case. The particle size also demonstrates a great effect on particle trajectory, then further has an impact on flame stability and temperature.

Li, Lei↗

NASA NTRS2007

Spectroscopic Detection of COClF in the Tropical and Mid-Latitude Lower Stratosphere

We report retrievals of COClF (carbonyl chlorofluoride) based on atmospheric chemistry experiment (ACE) solar occultation spectra recorded at tropical and mid-latitudes during 2004-2005. The COClF molecule is a temporary reservoir of both chlorine and fluorine and has not been measured previously by remote sensing. A maximum COClF mixing ratio of 99.7+/-48.0 pptv (10(exp -12) per unit volume, 1 sigma) is measured at 28km for tropical and subtropical occultations (latitudes below 20deg in both hemispheres) with lower mixing ratios at both higher and lower altitudes. Northern hemisphere mid-latitude mixing ratios (30-50degN) resulted in an average profile with a peak mixing ratio of 51.7+/-32.1 pptv, 1 sigma, at 27 km, also decreasing above and below that altitude. We compare the measured average profiles with the one reported set of in situ lower stratospheric mid-latitude measurements from 1986 and 1987, a previous two-dimensional (2-D) model calculation for 1987 and 1993, and a 2-D-model prediction for 2004. The measured average tropical profile is in close agreement with the model prediction; the northern mid-latitude profile is also consistent, although the peak in the measured profile occurs at a higher altitude (2.5-4.5km offset) than in the model prediction. Seasonal average 2-D-model predictions of the COClF stratospheric distribution for 2004 are also reported.

Rinsland, Curtis P.↗

DOE OSTI2020

Integration of Total-Sky Imager Data with a Physics-Based Smart Persistence Model for Intra-Hour Forecasting of Solar Radiation

Short-term solar forecasting models based solely on global horizontal irradiance (GHI) measurements are often unable to discriminate the forecasting of the factors affecting GHI from those that can be precisely computed by atmospheric models. Our previous study introduced a Physics-based Smart Persistence model for Intra-hour forecasting of solar radiation (PSPI) that decomposed the forecasting of GHI into the computation of extraterrestrial solar radiation and solar zenith angle and the forecasting of cloud albedo and cloud fraction. The extraterrestrial solar radiation and solar zenith angle were accurately computed by the Solar Position Algorithm (SPA) developed at the National Renewable Energy Laboratory (NREL). A cloud retrieval technique was used to estimate cloud albedo and cloud fraction from surface-based observations of GHI. With the assumption of persistent cloud structures, the cloud albedo and cloud fraction were predicted for future time steps using a two-stream approximation and a 5-minute exponential weighted moving average, respectively. The model evaluation indicated the estimation and forecast of cloud fraction mostly contributed to the uncertainty of the PSPI though it overcame the persistence and smart persistence models in all forecast time horizons between 5 and 60 minutes. This study aims to enhance the PSPI by ingesting surface-based observations of cloud fraction from a total sky imager (TSI). The estimation and forecast of cloud albedo is correspondingly improved by utilizing the cloud fraction observations and thus leads to more accurate GHI forecast. Various time-series analysis methods are also investigated on the forecasting of cloud fraction and cloud albedo for further improving the GHI forecast. These improvements are valuable for many applications, such as forecasting energy use for buildings, grid operations, and ultimately bringing down the cost of solar energy.

14 SOLAR ENERGY↗

NASA NTRS2018

Cloud Impacts on Photochemistry: Building a Climatology of Photolysis Rates from the Atmospheric Tomography Mission

Measurements from actinic flux spectroradiometers on board the NASA DC-8 during the Atmospheric Tomography (ATom) mission provide an extensive set of statistics on how clouds alter photolysis rates (J values) throughout the remote Pacific and Atlantic Ocean basins. J values control tropospheric ozone and methane abundances, and thus clouds have been included for more than three decades in tropospheric chemistry modeling. ATom made four profiling circumnavigations of the troposphere capturing each of the seasons during 2016–2018. This work examines J values from the Pacific Ocean flights of the first deployment, but publishes the complete Atom-1 data set (29 July to 23 August 2016). We compare the observed J values (every 3 s along flight track) with those calculated by nine global chemistry– climate/transport models (globally gridded, hourly, for a mid-August day). To compare these disparate data sets, we build a commensurate statistical picture of the impact of clouds on J values using the ratio of J -cloudy (standard, sometimes cloudy conditions) to J -clear (artificially cleared of clouds). The range of modeled cloud effects is inconsistently large but they fall into two distinct classes: (1) models with large cloud effects showing mostly enhanced J values aloft and or diminished at the surface and (2) models with small effects having nearly clear-sky J values much of the time. The ATom-1 measurements generally favor large cloud effects but are not precise or robust enough to point out the best cloud-modeling approach. The models here have resolutions of 50–200 km and thus reduce the occurrence of clear sky when averaging over grid cells. In situ measurements also average scattered sunlight over a mixed cloud field, but only out to scales of tens of kilometers. A primary uncertainty remains in the role of clouds in chemistry, in particular, how models average over cloud fields, and how such averages can simulate measurements.

actinic flux spectroradiometers↗

DOE OSTI2020

Hydrologic Model Files (PRMS) for Historical Conditions in the East River Watershed, Colorado between 1987-2019

The data package contains model input and output file for the East River watershed (85 km2) located in the headwaters of the Colorado River. The model is the U.S. Geological Survey (USGS) Precipitation-Modeling Runoff System (PRMS) for the historical, or baseline, condition (1987-2019) to assess the efficiency of summer monsoon rains to generate streamflow. The finite difference grid resolution is 100 m with elevations resampled from the USGS National Elevation Dataset. LANDFIRE raster data is used to derive parameters of dominant cover type, summer and winter cover density, canopy interception characteristics for snow and rain and transmission coefficient for shortwave radiation. Climate forcing uses minimum and maximum daily temperature lapse rates defined by the two Snow Telemetry (SNOTEL) stations located in the vicinity of the East River adjusted for aspect. Schofield SNOTEL snowfall is spatially distributed using LiDAR derived snow depth observations from the Airborne Snow Observatory (ASO) flown 4 April 2016. Snow depths are converted to snow water equivalent based on ground surveys and density modeling. Rainfall is spatially distributed using the monthly Parameter-elevation Relationships on Independent Slopes Model 30-year monthly averages. Simulated solar radiation is calibrated to match weather station observations. Model verification of snow water equivalent (SWE) accumulation and ablation relies on the 2018 and 2019 ASO maps and weather station snow depth. The data package contains model input files in text format, model documentation in pfd and a readme.txt file to help users understand the file structure to run the program.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS1989

The mean ozone profile and its temperature sensitivity in the upper stratosphere and lower mesosphere - An analysis of LIMS observations

Multiple simultaneous LIMS (Limb Infrared Monitor of the Stratosphere) observations of both O3 and temperature T in the upper stratosphere and lower mesosphere are analyzed to derive O3-T correlations for comparisons with model calculations. The zonally averaged O3 profile is compared with results from a simplified photochemical model that assumes O3 to be in a photochemical steady state. The model profile is systematically lower than the observed profile. Key parameters are identified in which changes can systematically increase the O3 profile in both the stratosphere and the mesosphere. The LIMS-derived values for the sensitivity of O3 to changes in T are compared with equilibrium model calculations which include the T-driven opacity feedback effect on photodissociation rate constants. The theoretical O3 response to temperature perturbations is investigated. It is shown how the O3-T sensitivity coefficient is affected by zonal and vertical advection terms as well as by the photochemical coupling between O3 and T.

Froidevaux, Lucien↗