SEARCH · Engineering Papers

Results for “forward modeling”

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 559 records · Page 31

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Combining geometric-optical and spectral invariants theories for modeling canopy fluorescence anisotropy

The spectral invariants theory ( p -theory) has received much attention in the field of quantitative remote sensing over the past few decades and has been adopted for modeling of canopy solar-induced chlorophyll fluorescence (SIF). However, the spectral invariant properties (SIP) in simple analytical formulae have not been applied for modeling canopy fluorescence anisotropy primarily because they are parameterized in terms of leaf total scattering, which precludes the differentiation between forward and backward leaf SIF emissions. In this study, we have developed the canopy-SIP SIF model by combining geometric-optical (GO) theory to account for asymmetric leaf SIF forward and backward emissions at the first-order scattering and by modeling multiple scattering based on the p-theory, thus avoiding the dependence on radiative transfer models. The applicability of the model simulations especially over 3D heterogeneous canopies was improved by incorporating canopy structure through multi-angular clumping index, and by modeling single scattering from the four components of the scene in view according to the GO approach. The results show good consistency with both the state-of-the-art SIF models and multi-angular field SIF observations over grass and chickpea canopies. Further, the coefficient of determination (R²) between the simulated SIF and field measurements was 0.75 (red) and 0.74 (far-red) for chickpea, and 0.65 (both red and far-red) for grass. The average relative error was approximately 3% for 1D homogeneous scenes when comparing the canopy-SIP SIF model simulations to the SCOPE model simulations, and around 4% for the 3D heterogeneous scene when comparing to the LESS model simulations. The results indicate that the proposed approach for separating asymmetric leaf SIF emissions is a robust way to keep a balance between satisfactory simulation accuracy and efficiency. Model simulations suggest that neglecting the leaf SIF asymmetry can lead to an underestimation of canopy red SIF by 6.3% to 42.6% for various leaf biochemical and canopy structural parameters. This study presents a simple but efficient analytical approach for canopy fluorescence modeling, with potential for large-scale canopy fluorescence simulations.

3D heterogeneous structure↗

DOE OSTI2021

Dynamical phase transitions in models of collective neutrino oscillations

Collective neutrino oscillations can potentially play an important role in transporting lepton flavor in astrophysical scenarios where the neutrino density is large, typical examples are the early universe and supernova explosions. It has been argued in the past that simple models of the neutrino Hamiltonian designed to describe forward scattering can support substantial flavor evolution on very short timescales t ≈ log(N) / (G Fρν ), with N the number of neutrinos, G F the Fermi constant and ρ ν the neutrino density. This finding is in tension with results for similar but exactly solvable models for which t ≈ √N/(G Fρν ) instead. In this work we provide a coherent explanation of this tension in terms of dynamical phase transitions (DPT) and study the possible impact that a DPT could have in more realistic models of neutrino oscillations and their mean-field approximation.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2023

Impurity gas monitoring using ultrasonic sensing and neural networks: forward and inverse problems

Ultrasonic sensing is a non-invasive technique for monitoring impurity gas composition in various industrial applications where safety and regulatory compliance are crucial. In this study, ultrasonic sensing and neural networks were used to analyze impurity gases (i.e., air and argon) in helium. An experimental platform was established to acquire ultrasonic data. In the forward problem, an artificial neural network (ANN) model was used to forecast the response and time-of-flight (TOF) based on the excitation, and argon and air concentrations. The inverse problem was solved using a convolutional neural network (CNN) to predict the argon and air concentrations given the ultrasonic response and excitation. The results showed that the ANN accurately predicted the ultrasonic response and the change in TOF with concentration. As the air concentration was increased from 0 to 9.8%, the TOF sensitivity to detect argon decreased by 39.8% and 16.1% from ANN and sound speed theory, respectively. The CNN demonstrated high accuracy in predicting concentrations for inputs in the testing dataset. The application of the trained CNN indicated that it over-predicts air concentration while under-predicting the argon concentration. To improve accuracy, the predicted air and argon concentrations should be corrected by -0.992% and 1.027% bias, respectively.

47 OTHER INSTRUMENTATION↗

DOE OSTI2025

A thermodynamically consistent discretization of 1D thermal-fluid models using their metriplectic 4-bracket structure

Thermodynamically consistent models in continuum physics, i.e. models which satisfy the first and second laws of thermodynamics, may be expressed using the metriplectic formalism. In this work, we leverage the structures underlying this modeling formalism to preserve thermodynamic consistency in discretizations of a fluid model. The procedure relies (1) on ensuring that the spatial semi-discretization retains certain symmetries and degeneracies of the Poisson and metriplectic 4-brackets, and (2) on the use of an appropriate energy conserving time-stepping method. Here, the minimally simple yet nontrivial example of a one-dimensional thermal-fluid model is treated. It is found that preservation of the requisite symmetries and degeneracies of the 4-bracket is relatively simple to ensure in Galerkin spatial discretizations, suggesting a path forward for thermodynamically consistent discretizations of more complex fluid models using more specialized Galerkin methods.

Hamiltonian structure↗

DOE OSTI2022

Study of very forward energy and its correlation with particle production at midrapidity in pp and p-Pb collisions at the LHC

The energy deposited at very forward rapidities (very forward energy) is a powerful tool for characterising proton fragmentation in pp and p-Pb collisions. The correlation of very forward energy with particle production at midrapidity provides direct insights into the initial stages and the subsequent evolution of the collision. Furthermore, the correlation with the production of particles with large transverse momenta at midrapidity provides information complementary to the measurements of the underlying event, which are usually interpreted in the framework of models implementing centrality-dependent multiple parton interactions. Results about very forward energy, measured by the ALICE zero degree calorimeters (ZDCs), and its dependence on the activity measured at midrapidity in pp collisions at $\sqrt{s}$ = 13 TeV and in p-Pb collisions at $\sqrt{s{NN}}$ = 8.16 TeV are discussed. The measurements performed in pp collisions are compared with the expectations of three hadronic interaction event generators: PYTHIA 6 (Perugia 2011 tune), PYTHIA 8 (Monash tune), and EPOS LHC. These results provide new constraints on the validity of models in describing the beam remnants at very forward rapidities, where perturbative QCD cannot be used.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2024

Using a 1D Radially Symmetric Coda Envelope Model for Robust Moment Magnitude ( M w ) Estimation in Iraq’s Tectonically Diverse Zones

Robust estimation of moment magnitude (M w ) can be challenging for Iraq due to the strong lateral heterogeneity across diverse tectonic zones. We aim to improve moment magnitude estimation by investigating the reliability of using a 1D coda envelope model in diverse tectonic zones of different lateral effects and offer a way forward for reliable estimates of M w for small events that are difficult to waveform model. Iraq comprises two main tectonic zones: (1) the Outer platform, consisting of the northwestern Zagros fold-thrust belt and the Mesopotamian foredeep, and (2) the Inner Arabian platform which is overlain by the Iraqi desert. A simple 1D coda envelope model was used because coda waves have a low sensitivity to the source and path heterogeneity. Three separate coda calibrations were conducted to investigate the robustness of a single 1D calibration to fit the country: Whole-region calibration, Zagros calibration, and Mesopotamia calibration. In the whole-region calibration, we used stations from both the Zagros and Mesopotamia zones. In the two other calibration models, we used only stations that were in those particular zones. Ground-truth reference spectra derived from the coda spectral ratio method were used to constrain high-frequency site terms. There was no drastic difference when comparing the moment magnitudes calculated from the waveform modeling and the three calibration models. The results show that the 1D coda envelope model is a reliable method even for a region with diverse tectonic zones. Therefore, we recommend using the whole-region calibration model for moment magnitude estimation that provides more complete path coverage and avoids biases introduced by path correction failures. The proposed calibration is a fundamental step in updating the comprehensive earthquake catalog and probabilistic seismic hazard assessments for Iraq.

58 GEOSCIENCES↗

DOE OSTI2022

Retail Rate Projections for Long-Term Electricity System Models

Electricity prices reported in most electricity-system planning studies leave out many price components and do not translate into retail rates, making it difficult to interpret how projected electricity system changes will impact costs to consumers. Full transmission costs are left out of many studies; distribution and administration costs are similarly excluded or highly simplified. Here, we present a detailed bottom-up accounting method for projecting future retail electricity rates in the United States. Making the simplifying assumption that each state is served by an investor-owned utility (IOU), we translate projected generation and transmission capacity and costs from the Regional Energy Deployment System (ReEDS) capacity-expansion model into IOU balance sheet expenditures, accounting for depreciation, taxes, and the breakdown between operating and capitalized (rate-based) expenses. Distribution, administration, and intra-regional transmission costs are projected forward based on empirical trends over the past decade. Modeled bottom-up electricity rates are compared to historical rates from 2010-2019, and the sensitivity of modeled rates to a range of financing and modeling assumptions is explored. Distribution and administration rate components account for roughly 40% (4.4 ¢/kWh) of the projected national-average retail rate over 2020-2050 under central assumptions.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

DOE OSTI2026

Pedestal origin and extrapolation of high-density small edge-localised-modes peak parallel energy fluence in ITER and SPARC

Experimental analysis and simulations with the BOUT++ code show that small edge-localised modes (ELMs) in reactor-relevant high-density regimes originate in a region close to the separatrix and only marginally perturb the pedestal structure. The measured divertor peak parallel energy fluence (ε ∥,peak ) for a database of small ELM scenarios in DIII-D and ASDEX Upgrade can be reproduced, within 40 % accuracy on average, if an ad hoc modification of the Eich peak parallel ELM energy fluence model is applied to account for the small ELM pedestal birth location. This allows for first-order extrapolation of small-ELM divertor ε ∥,peak to ITER and SPARC, resulting in values that satisfy the nominal melting threshold of tungsten monoblocks of 12 MJ m −2 . The findings reported in this study, both via modelling and direct measurements, constitute a step forward in assessing small ELMs in high edge-collisionality scenarios as a viable plasma regime for the operation of next-generation fusion machines.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Effect of realistic out-of-plane dopant potentials on the superfluid density of overdoped cuprates

Recent experimental papers on hole-doped overdoped cuprates have argued that a series of observations showing unexpected behavior in the superconducting state imply the breakdown of the quasiparticle-based Landau–BCS paradigm in that doping range. In contrast, some of the present authors have argued that a phenomenological “dirty d-wave” theoretical analysis explains essentially all aspects of thermodynamic and transport properties in the superconducting state, provided the unusual effects of weak, out-of-plane dopant impurities are properly accounted for. Here we attempt to place this theory on a more quantitative basis by performing ab initio calculations of dopant impurity potentials for LSCO and Tl-2201. These potentials are more complex than the pointlike impurity models considered previously, and require calculation of forward scattering corrections to transport properties. Including realistic, ARPES-derived band structures, Fermi liquid renormalizations, and vertex corrections, we show that the theory can explain semiquantitatively the unusual superfluid density measurements of the two most studied overdoped materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2024

Higgs to b b ¯ from vector boson fusion for high-scale physics

Vector boson fusion is arguably the most direct collider probe of electroweak symmetry breaking. Typically, the signature includes two forward/backward jets with low transverse momenta with a scale that is set by the mass of the vector boson. For this reason, an upper cut is used when searching for vector boson fusion processes in the Standard Model. Alternatively, the upper cut on the forward jets can be removed and the high-momentum exchange region of vector boson fusion can be studied. This phase space region has sensitivity to new physics via higher dimensional operators and form factors. In this work, we study the high-momentum region of the vector boson fusion channel where the Higgs decays to b b ¯ . We show that, depending on the form of new physics, the limits on the new physics scale range from 0.5 to 1.8 TeV. Published by the American Physical Society 2024

Astronomy & Astrophysics↗

DOE OSTI2021

Regulatory Implications of Embedded Grid Energy Storage

Electricity is unique among commodities in that its supply chain was developed without a storage component. Every other resource commodity has the ability to store excess quantities built into its supply chain – in the form of granaries, warehouses, reservoirs, etc. This embedded storage creates a buffer for mismatches between supply and demand, stabilizing prices and protecting customers. Recent technological advances in cost-competitive energy storage technologies that are scalable and flexible have made the concept of embedded electrical storage feasible, but several regulatory questions remain. This paper summarizes energy regulatory structures in the U.S. and the implications that they would have for embedded storage on the electric grid. Regulatory challenges for embedded storage include the lack of underlying standards, inclusion in planning processes, ownership models, compensation, and metrics. Two possible pathways forward are explored: an incremental one including regulatory guidance regarding the evaluation of embedded storage in existing planning processes, and a complex one involving the revision of existing reliability standards or the development of a new standard specific to embedded storage.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2022

Adapting Technology Learning Curves for Prospective Techno-Economic and Life Cycle Assessments of Emerging Carbon Capture and Utilization Pathways

Comparisons of emerging carbon capture and utilization (CCU) technologies with equivalent incumbent technologies are necessary to support technology developers and to help policy-makers design appropriate long-term incentives to mitigate climate change through the deployment of CCU. In particular, early-stage CCU technologies must prove their economic viability and environmental reduction potential compared to already-deployed technologies. These comparisons can be misleading, as emerging technologies typically experience a drastic increase in performance and decrease in cost and greenhouse gas emissions as they develop from research to mass-market deployment due to various forms of learning. These changes complicate the interpretation of early techno-economic assessments (TEAs) and life cycle assessments (LCAs) of emerging CCU technologies. The effects of learning over time or cumulative production themselves can be quantitatively described using technology learning curves (TLCs). While learning curve approaches have been developed for various technologies, a harmonized methodology for using TLCs in TEA and LCA for CCU in particular is required. To address this, we describe a methodology that incorporates TLCs into TEA and LCA to forecast the environmental and economic performance of emerging CCU technologies. This methodology is based on both an evaluation of the state of the art of learning curve assessment and a literature review of TLC approaches developed in various manufacturing and energy generation sectors. Additionally, we demonstrate how to implement this methodology using a case study on a CO 2 mineralization pathway. Finally, commentary is provided on how researchers, technology developers, and LCA and TEA practitioners can advance the use of TLCs to allow for consistent, high-resolution modeling of technological learning for CCU going forward and enable holistic assessments and fairer comparisons with other climate technologies.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

DOE OSTI2025

Application of physics-informed neural networks (PINNs) solution to coupled thermal and hydraulic processes in silty sands

Abstract The accurate modeling of water and heat transport in soils is crucial for both geo-environmental and geothermal engineering. Traditional modeling methods are problematic because they require well-defined boundaries and initial conditions. Recently, physics-informed neural networks (PINNs), which incorporate partial differential equations (PDEs) to solve forward and inverse problems, have attracted increasing attention in machine learning research. In this study, we applied PINNs to tackle hydraulic and thermal transport coupling forward problems in silty sands. A fully connected deep neural network was utilized for training. This neural network model leverages automatic differentiation to apply the governing equations as constraints, based on the mathematical approximations established by the neural network itself. We conducted forward problems and compared the solutions derived from PINNs with those from Finite Element Method (FEM) simulations. The forward problem results demonstrate the PINNs model’s capability in predicting hydraulic transport, heat transport, and thermal–hydraulic coupling in silty sands under various boundary conditions. The PINNs exhibited great performance in simulating the thermal–hydraulic coupling problem. The accuracy of the PINNs solutions shows its potential for simulation in geotechnical engineering.

Feng, Yuan↗

DOE OSTI2020

Code for Calculating Regional Seismic Travel Time

The RSTT software computes predictions of the travel time of seismic energy traveling from a source to a receiver through 2.5D models of the seismic velocity distribution within the Earth. The two primary applications for the RSTT library are tomographic inversion studies and seismic event location calculations. In tomographic inversions studies, a seismologist begins with number of source-receiver travel time observations and an initial starting model of the velocity distribution within the Earth. A forward travel time calculator, such as the RSTT library, is used to compute predictions of each observed travel time and all of the residuals (observed minus predicted travel time) are calculated. The Earth model is then modified in some systematic way with the goal of minimizing the residuals. The Earth model obtained in this way is assumed to be a better model than the starting model if it has lower residuals. The other major application for the RSTT library is seismic event location. Given an Earth model, an initial estimate of the location of a seismic event, and some number of observations of seismic travel time thought to have originated from that event, location codes systematically modify the estimate of the location of the event with the goal of minimizing the difference between the observed and predicted travel times. The second application, seismic event location, is routinely implemented by the military as part of its effort to monitor the Earth for nuclear tests conducted by foreign countries.

BALLARD, SANFORD↗

DOE OSTI2020

Regional Seismic Travel Time v.3.2.0

BALLARD, SANFORD↗

DOE OSTI2023

Verification and Benchmarking of High Fidelity Physics Peat Smoldering Model

Peat fires are a major contributor to greenhouse gas emissions. The estimates of these emissions currently contain major uncertainties, due to the difficulty of determining the mass of peat burned in a fire. To address these uncertainties, we develop a computational physics-based peat smoldering model, which will be leveraged for high-fidelity quantitative estimates of peat fire emissions relevant to climate change. We present the verification of the 2-D axisymmetric model, a first step towards developing a full 3-D model. Verification includes the solution verification against a literature model for the 0-D smoldering case and verification of the heat transfer problem in 1-D and 2-D. Also presented is the effect of reaction mechanism on the smoldering model, for which we found a relatively simple three-step reaction mechanism is able to capture key behavior. These verification results provide the foundation for moving forward with validation against experimental data of the 2-D model.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2022

PyApprox: Enabling efficient model analysis

PyApprox is a Python-based one-stop-shop for probabilistic analysis of scientific numerical models. Easy to use and extendable tools are provided for constructing surrogates, sensitivity analysis, Bayesian inference, experimental design, and forward uncertainty quantification. The algorithms implemented represent the most popular methods for model analysis developed over the past two decades, including recent advances in multi-fidelity approaches that use multiple model discretizations and/or simplified physics to significantly reduce the computational cost of various types of analyses. Simple interfaces are provided for the most commonly-used algorithms to limit a user’s need to tune the various hyper-parameters of each algorithm. However, more advanced work flows that require customization of hyper-parameters is also supported. An extensive set of Benchmarks from the literature is also provided to facilitate the easy comparison of different algorithms for a wide range of model analyses. This paper introduces PyApprox and its various features, and presents results demonstrating the utility of PyApprox on a benchmark problem modeling the advection of a tracer in ground water.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2021

Transport Modeling of Kinetically Limited Microscale Extraction Systems: Droplet and Supported Liquid Membrane Separations

Microfluidic droplet extraction devices can be used to efficiently measure solute phase transfer kinetics in liquid–liquid extraction systems. In this work, a droplet extraction apparatus was constructed using off-the-shelf and custom 3D printed parts for less than 110 USD. The interfacial kinetics of the transfer of Pu(IV) between an aqueous nitric acid phase and a 5 v/v% tridodecylamine (TDA) in toluene phase were characterized by determining the extent of Pu(IV) transfer between phases at different flow rates. Forward interfacial rate constants calculated from this data using an analytical model are 5.6 × 10 –3 mm s –1 under extraction (aqueous-to-organic solute transfer) conditions and 3.3 × 10 –3 mm s –1 under stripping (organic-to-aqueous solute transfer) conditions. The analytical model is valid when the rate of solute transfer across the interface between the two liquid phases has pseudo-first-order reversible reaction kinetics and is slow relative to the rate of solute transport to the interface. The change in the outlet concentration with residence time determined using the analytical model was compared with the results of a numerical model accounting for advection/diffusion to the interface in addition to kinetically limited interfacial solute transfer. These calculations revealed that for systems with forward interfacial rate constants greater than approximately 10 –2 mm s –1 , the analytical model must be supplemented with a numerical transport model accounting for an additional advection/diffusion mass transfer resistance. This threshold corresponds to a Damköhler number of one. The measured interfacial kinetic rate constants for Pu(IV)/TDA were used in a supported liquid membrane (SLM) numerical model, which showed that the Pu(IV) phase transfer kinetics are fast enough for use in a practical separation.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗