Engineering Papers⌕ Search

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.

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

At least 289 records · Page 16

A hybrid Penman-Monteith and machine learning model for simulating evapotranspiration and its components

Integrating physical processes with machine learning has advanced evapotranspiration (ET) simulation, yet most hybrid models fail to partition total ET into its components: soil evaporation (E) and vegetation transpiration (T). This study introduces Residual Neural Network–Penman–Monteith (RNN-PM), a novel hybrid dual-source ET model designed to overcome this limitation. The model synergizes the physically-based Penman–Monteith framework with three specialized residual neural networks trained to estimate key conductance parameters (canopy conductance, soil surface conductance, and aerodynamic conductance). Furthermore this explicit parameterization allows for the direct partitioning of total ET. Validation at National Ecological Observatory Network (NEON) flux sites using high-frequency partitioned E and T shows that RNN-PM reliably reproduces ET and the transpiration fraction (T/ET). For ET, the model achieves an average Kling–Gupta efficiency (KGE) of 0.89 and a root-mean-square error (RMSE) of 0.55 mm/day; for T/ET, the KGE is 0.87 with an RMSE of 0.06. Furthermore, RNN-PM demonstrates robust generalization, accurately simulating ET and its components well beyond the initial training dataset, even under extreme climatic conditions. This study extended the analysis by comparing the RNN-PM model with seven established dual-source ET models. The results indicate that RNN-PM outperforms both conventional machine learning models and purely physical process-based models in simulating ET components in most cases. Among the purely physical process-based dual-source models, those based on surface temperature decomposition showed improved performance as the leaf area index (LAI) decreased when evaluated against high-frequency ET component datasets. In contrast, the performance of conductance-based dual-source models declined with decreasing LAI. Although purely machine learning-based models can produce relatively accurate simulations of ET components, they often exhibit limited generalization capability, an issue that the RNN-PM model effectively overcomes. Ultimately, the RNN-PM model represents a significant advance in simulating ET components, offering a novel and scalable approach for improving the representation of land–atmosphere interactions in Earth system models.

54 ENVIRONMENTAL SCIENCES↗

A Comparison of Power Take-Off Architectures for Wave-Powered Reverse Osmosis Desalination of Seawater with Co-Production of Electricity

Several power take-off (PTO) architectures for wave-powered reverse osmosis (RO) desalination of seawater are introduced and compared based on the annual average freshwater production and the size of the components, which strongly relate to the costs of the system. The set of architectures compared includes a novel series-type PTO architecture not previously considered. These seawater hydraulic PTO architectures are composed of a WEC-driven pump, an RO module, an intake charge pump driven by an electric motor, and a hydraulic motor driving an electric generator for electric power production. This study is performed using an efficient two-way coupled steady-state model for the average performance of the system in a given sea state, including freshwater permeate production, electric power production, and electric power consumption. A multi-objective design problem is formulated for the purposes of this comparative study, with the objectives of maximizing annual freshwater production, minimizing the displacement of the WEC-driven pump, and minimizing the installed RO membrane area. This establishes a framework for comparison in the absence of a mature techno-economic model. The requirement that the system produces enough electric power to meet its consumption is applied as a constraint on the operation of the system. The oscillating wave surge converter Oyster 1 is assumed as the WEC. Weights on performance of the system in a given sea state are based on historical data from Humboldt Bay, CA. This study finds that (1) architectures in a series configuration allow for a reduction in the WEC-driven pump size of 59–92% compared to prior work, (2) varying the displacement of the WEC-driven pump between sea conditions does not provide any significant advantage in performance, and (3) varying the active RO membrane area between sea condition offers improvements between 7% and 41% in each design objective.

Energy & Fuels↗

GIGA-Lens: Fast Bayesian Inference for Strong Gravitational Lens Modeling

We present GIGA-Lens: a gradient-informed, GPU-accelerated Bayesian framework for modeling strong gravitational lensing systems, implemented in TensorFlow and JAX. The three components, optimization using multistart gradient descent, posterior covariance estimation with variational inference, and sampling via Hamiltonian Monte Carlo, all take advantage of gradient information through automatic differentiation and massive parallelization on graphics processing units (GPUs). We test our pipeline on a large set of simulated systems and demonstrate in detail its high level of performance. The average time to model a single system on four Nvidia A100 GPUs is 105 s. The robustness, speed, and scalability offered by this framework make it possible to model the large number of strong lenses found in current surveys and present a very promising prospect for the modeling of ${ \mathcal O }({10}^{5})$ lensing systems expected to be discovered in the era of the Vera C. Rubin Observatory, Euclid, and the Nancy Grace Roman Space Telescope.

79 ASTRONOMY AND ASTROPHYSICS↗

Assessing China’s efforts to pursue the 1.5 degrees C warming limit

Given the increasing interest in keeping global warming below 1.5 ?, a key question is what this would mean for China’s emission pathway, energy restructuring and decarbonization. By conducting a multi-model study, we find that the 1.5 ?-consistent goal would require China to reduce its carbon emissions and energy consumption by over 90% and 39%, respectively, compared to the No Policy case. Negative emission technologies play an important role in achieving near-zero emissions, with captured carbon averagely accounting for 23% of the total reductions in 2050. Our multi-model comparisons reveal large differences in necessary emission reductions across sectors, while what’s consistent is that the power sector is required to achieve full decarbonization by 2050. The cross-model averages indicate that China’s accumulated policy costs may amount to 2.75-5.73% of its GDP by 2050, given the 1.5 ? warming limit.

1.5 degrees, China, Paris agreement, ENGAGE↗

Modeling laser-wakefield accelerators using the time-averaged ponderomotive approximation in a Lorentz boosted frame

Future, high-fidelity simulations of multi-GeV-class laser Wakefield accelerators (LWFAs) will need to model the propagation of high-intensity laser drivers over meter-scale plasmas with high spatial and temporal resolutions, thus requiring high amounts of computational resources. Various techniques have been devised over the years to reduce the computational cost of such simulations, including the time-averaged ponderomotive approximation, and the use of the Lorentz boosted frame technique. In this paper we discuss the combination of these two computational techniques, highlighting the resulting significant reduction in the computational cost of LWFA simulations and the limitations of this approach. The combination of the two techniques can potentially become essential for the modeling of a multi-TeV, LWFA-based collider.

Laser Wakefield Acceleration↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

Non-intrusive reduced-order modeling for dynamical systems with spatially localized features

This work presents a non-intrusive reduced-order modeling framework for dynamical systems with spatially localized features characterized by slow singular value decay. The proposed approach builds upon two existing methodologies for reduced and full-order non-intrusive modeling, namely Operator Inference (OpInf) and sparse Full-Order Model (sFOM) inference. We decompose the domain into two complementary subdomains that exhibit fast and slow singular value decay. The dynamics of the subdomain exhibiting slow singular value decay are learned with sFOM while the dynamics with intrinsically low dimensionality on the complementary subdomain are learned with OpInf. The resulting, coupled OpInf-sFOM formulation leverages the computational efficiency of OpInf and the high resolution of sFOM, and thus enables fast non-intrusive predictions for conditions beyond those sampled in the training data set. A novel regularization technique with a closed-form solution based on the Gershgorin disk theorem is introduced to promote stable sFOM and OpInf models. We also provide a data-driven indicator for subdomain selection and ensure solution smoothness over the interface via a post-processing interpolation step. We evaluate the efficiency of the approach in terms of offline and online speedup through a quantitative, parametric computational cost analysis. We demonstrate the coupled OpInf-sFOM formulation for two test cases: a one-dimensional Burgers’ model for which accurate predictions beyond the span of the training snapshots are presented, and a two-dimensional parametric model for the Pine Island Glacier ice thickness dynamics, for which the OpInf-sFOM model achieves an average prediction error on the order of 1% with an online speedup factor of approximately 8$\times$ compared to the numerical simulation.

42 ENGINEERING↗

Modeling Yield, Biogenic Emissions, and Carbon Sequestration in Southeastern Cropping Systems With Winter Carinata

Sustainable aviation fuel (SAF) production from lipids is a technologically mature approach for replacing conventional fossil fuel use in the aviation sector, and there is increasing demand for such feedstocks. The oilseed Brassica carinata (known as Ethiopian mustard or simply carinata) is a promising SAF feedstock that can be grown as a supplemental cash crop over the winter fallow season of various annual crop rotations in the Southeast US, avoiding land use changes and potentially achieving some of the soil carbon sequestration and ecosystem service benefits of winter cover crops. However, carinata may require more intensive management than traditional cover crops, potentially leading to additional soil greenhouse gas (GHG) emissions through increased carbon losses from soil tillage and nitrous oxide (N 2 O) emissions from nitrogen fertilizer application. In this work, the 2017 version of the process-based DayCent ecosystem model was used to establish initial expectations for the total regional SAF production potential and associated soil GHG emissions when carinata is integrated as a winter crop into the existing crop rotations across its current suitability range in southern Alabama, southern Georgia, and northern Florida. Using data from academic and industry carinata field trials in the region, DayCent was calibrated to reproduce carinata yield, nitrogen response, harvest index, and biomass carbon-to-nitrogen ratio. The resulting model was then used to simulate the integration of carinata every third winter across all 2.1 Mha of actively cultivated cropland in the study area. The model predicted regional average yields of 2.9–3.0 Mg carinata seed per hectare depending on crop management assumptions. That results in the production of more than two million Mg of carinata seed annually across the study area, enough to supply approximately one billion liters of SAF. Conventional management of carinata led to only modest increases in soil carbon storage that were largely offset by additional N2O emissions. Climate-smart management via adopting no-till carinata establishment or using poultry litter as a nitrogen source resulted in a substantial net soil GHG sink (0.23–0.31 Mg CO 2 e ha -1 y -1 , or 0.24–0.32 Mg CO2e per Mg of seed produced) at the farms where carinata is cultivated.

54 ENVIRONMENTAL SCIENCES↗

Accuracy of the time-averaged ponderomotive approximation for laser-plasma accelerator modeling

Reliable modeling of laser-plasma accelerators, where a short and intense laser pulse propagates in an underdense plasma over long distances, is a computationally challenging task. This is due to the great disparity among the scales involved in the modeling, ranging from the micrometer scale of the laser wavelength to, for instance, the meter scale of the laser-plasma interaction length for a multi-GeV-class laser-plasma accelerator. To reduce such imbalance, the time-averaged ponderomotive approximation may be used, where the plasma particle dynamics is analytically averaged over the laser frequency, and only spatiotemporal scales associated with the laser envelope are retained in the calculations, resulting in significant computational savings. Here, we characterize the accuracy and robustness of the time-averaged ponderomotive approximation for a range of laser parameters of interest for present and future laser-plasma accelerators, and we show that the error introduced by the averaging process is small in all relevant cases.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗