Engineering Papers⌕ Search

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.

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

At least 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
NASA NTRS2002

Plasmaspheric Density Troughs: Global IMAGE EUV Observations and Analysis via Global Core Plasma Modeling

To date, the IMAGE EUV camera has observed several plasmaspheric density trough features inside the plasmapause under a wide range of geomagnetic activity. From the perspective of EUV, a density trough feature appears as a channel of diminished pixel counts which spans a width of L-shell (DELTA L) and magnetic local time (MLT) inside the plasmapause. Plasmaspheric density troughs are found to be morphologically complex possessing considerable spatial and temporal variability. We present an analysis of the evolution of trough DELTA L and MLT extent as functions of associated D (sub ST) and K (sub p) history. Trough features range in size from 0.16 less than or equal to DELTA L less than or equal to 1.2 with azimuthal extent from 1500 less than or equal to MLT less than or equal to 1200. All cases of plasmaspheric density troughs studied to date appear to have evolved as a result of the inner edge of the afternoon/evening plasma drainage plume being wrapped around through the nightside plasmasphere. The structure of plasmaspheric density trough features is further probed by analyzing simulated EUV images produced by forward modeling artificially introduced regions of depleted density into both static and dynamic global core plasmaspheric models. Forward modeling suggests that (1) L-shell refilling of density troughs during storm recovery can be modeled as filling from the ionosphere toward the equator (i.e., bottom-up refilling), and (2) that an erosion process is operating within flux tubes beyond the outer L-shell wall of the observed density troughs.

Adrian, M. L.↗
DOE OSTI2021

Bayesian inference of the dense-matter equation of state encapsulating a first-order hadron-quark phase transition from observables of canonical neutron stars

The remarkable progress in recent multimessenger observations of both isolated neutron stars (NSs) and their mergers has provided some of the much needed data to improve our understanding about the equation of state (EOS) of dense neutron-rich matter. Various EOSs with or without some kinds of phase transitions from hadronic to quark matter (QM) have been widely used in many forward modelings of NS properties. Direct comparisons of these predictions with observational data sometimes also using χ 2 minimizations have provided very useful constraints on the model EOSs. Furthermore, it is normally difficult to perform uncertain quantifications and analyze correlations of the EOS model parameters involved in forward modelings especially when the available data are still very limited.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Time-Lapse Integration at FWU: Fluids, Rock Physics, Numerical Model Integration, and Field Data Comparison

We present the current status of time-lapse seismic integration at the Farnsworth (FWU) CO2 WAG (water-alternating-gas) EOR (Enhanced Oil Recovery) project at Ochiltree County, northwest Texas. As a potential carbon sequestration mechanism, CO2 WAG projects will be subject to some degree of monitoring and verification, either as a regulatory requirement or to qualify for economic incentives. In order to evaluate the viability of time-lapse seismic as a monitoring method the Southwest Partnership (SWP) has conducted time-lapse seismic monitoring at FWU using the 3D Vertical Seismic Profiling (VSP) method. The efficacy of seismic time-lapse depends on a number of key factors, which vary widely from one application to another. Most important among these are the thermophysical properties of the original fluid in place and the displacing fluid, followed by the petrophysical properties of the rock matrix, which together determine the effective elastic properties of the rock fluid system. We present systematic analysis of fluid thermodynamics and resulting thermophysical properties, petrophysics and rock frame elastic properties, and elastic property modeling through fluid substitution using data collected at FWU. These analyses will be framed in realistic scenarios presented by the FWU CO2 WAG development. The resulting fluid/rock physics models will be applied to output from the calibrated FWU compositional reservoir simulation model to forward model the time-lapse seismic response. Modeled results are compared with field time-lapse seismic measurements and strategies for numerical model feedback/update are discussed. While mechanical effects are neglected in the work presented here, complementary parallel studies are underway in which laboratory measurements are introduced to introduce stress dependence of matrix elastic moduli.

4D↗
NASA NTRS2018

Active and Passive Radiative Transfer Modeling of the Olympic Mountains Experiment

Sensor forward models are an important tool for interpreting remote sensing observations of geophysical phenomena. By implementing a three-dimensional framework, we can simulate and analyze observations from various sensors on disparate platforms. To demonstrate our model framework, we simulate observations from the Olympic Mountains Experiment (OLYMPEX). The use of cloud model simulations allows us to understand sensor response to cloud ice, falling snow, and other processes and features, and the application of model tools to observations allows us to quantify precipitation.MIIST 3D Forward ModelThe Multi-Instrument Inverse Solver Testbed(MIIST) uses the Atmospheric Radiative TransferSimulator (ARTS) for solving the vector radiativetransfer (RT) equation in up to three spatialdimensions within a spherical geometry• Gas absorptiono Line-by-line calculationso Fast transmittance tables• Hydrometeor scattering solverso Discrete ordinateo RT4 (Evans, 1D)o Radar Single Scattering (1D or 3D)o Monte Carlo (3D)Scattering TablesHigh-fidelity hydrometeor scatteringtables are necessary for accurateand consistent forward modeling ofmulti-frequency observations• Requires full Stokes matriceso And absorption vector• Randomly oriented particleso Discrete Dipole Approximationo Characteristic Basis Function Method(coming soon)• Horizontally-oriented plateso Invariant Imbedding T-matrix MethodCloud Resolving SimulationsCloud resolving simulations (e.g.,NU-WRF) supply output consistentwith ARTS needs• Atmospheric Informationo Temperatureo Pressure / heighto Water vapor• Hydrometeor Profileso ARTS architecture ripe for explicit binmicrophysics• Examples use Morrison 2M schemeThe Olympic Mountains Experiment (OLYMPEX)Validation for GPM of mid-latitudefrontal systems approaching nearcoastalmountains from the ocean• Large collection of ground-based andairborne sensorso Radarso Radiometerso In situ• Contemporaneous with RADEXo Two sets of radar at same frequenciesRadiometer Simulation (3 km NUWRF, 20151203, 15:00)2018.12.14 7Simulate 166 GHz polarizationdifference• Corresponds to the presence of aligned icecrystals• Look at trends for both simulations andobservations• Simulations can tolerate lower resolutiono Larger domainSimulations from Observations: OLYMPEXSimulate sensor response usinggeophysical retrievals as input• Single frequency radar retrievals• Multiple scattering enhancementapparent at W band• Spatially dependent phenomenonModeling Application: 1D Retrievals03 December 2015• DC-8 and ER-2 flightso Focus on APR-3 (DC-8)• Citationo Stacked microphysics legso Qualitative comparisonso Range of frozen habitso Presence of supercooledliquid cloudsResults• Retrievals match probeso Good qualitative match• Bands of increasedreflectivity correspond tolarge Dm and highaggregate fraction• Significant amounts ofsupercooled liquid water

Adams, Ian S.↗
DOE OSTI2022

A Reanalysis of Public Galactic Bulge Gravitational Microlensing Events from OGLE-III and -IV

Abstract Modern surveys of gravitational microlensing events have progressed to detecting thousands per year, and surveys are capable of probing Galactic structure, stellar evolution, lens populations, black hole physics, and the nature of dark matter. One of the key avenues for doing this is the microlensing Einstein radius crossing time (t E) distribution. However, systematics in individual light curves as well as oversimplistic modeling can lead to biased results. To address this, we developed a model to simultaneously handle the microlensing parallax due to Earth's motion, systematic instrumental effects, and unlensed stellar variability with a Gaussian process model. We used light curves for nearly 10,000 OGLE-III and -IV Milky Way bulge microlensing events and fit each with our model. We also developed a forward model approach to infer the t E distribution by forward modeling from the data rather than using point estimates from individual events. We find that modeling the variability in the baseline removes a source of significant bias in individual events, and the previous analyses overestimated the number of t E > 100 day events due to their oversimplistic model ignoring parallax effects. We use our fits to identify the hundreds filling a regime in the microlensing parameter space that are 50% pure of black holes. Finally, we have released the largest-ever catalog of Markov Chain Monte Carlo parameter estimates for microlensing events.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Efficient high-dimensional variational data assimilation with machine-learned reduced-order models

Abstract. Data assimilation (DA) in geophysical sciences remains the cornerstone of robust forecasts from numerical models. Indeed, DA plays a crucial role in the quality of numerical weather prediction and is a crucial building block that has allowed dramatic improvements in weather forecasting over the past few decades. DA is commonly framed in a variational setting, where one solves an optimization problem within a Bayesian formulation using raw model forecasts as a prior and observations as likelihood. This leads to a DA objective function that needs to be minimized, where the decision variables are the initial conditions specified to the model. In traditional DA, the forward model is numerically and computationally expensive. Here we replace the forward model with a low-dimensional, data-driven, and differentiable emulator. Consequently, gradients of our DA objective function with respect to the decision variables are obtained rapidly via automatic differentiation. We demonstrate our approach by performing an emulator-assisted DA forecast of geopotential height. Our results indicate that emulator-assisted DA is faster than traditional equation-based DA forecasts by 4 orders of magnitude, allowing computations to be performed on a workstation rather than a dedicated high-performance computer. In addition, we describe accuracy benefits of emulator-assisted DA when compared to simply using the emulator for forecasting (i.e., without DA). Our overall formulation is denoted AIEADA (Artificial Intelligence Emulator-Assisted Data Assimilation).

58 GEOSCIENCES↗
NASA NTRS2006

Computing Fault Displacements from Surface Deformations

Simplex is a computer program that calculates locations and displacements of subterranean faults from data on Earth-surface deformations. The calculation involves inversion of a forward model (given a point source representing a fault, a forward model calculates the surface deformations) for displacements, and strains caused by a fault located in isotropic, elastic half-space. The inversion involves the use of nonlinear, multiparameter estimation techniques. The input surface-deformation data can be in multiple formats, with absolute or differential positioning. The input data can be derived from multiple sources, including interferometric synthetic-aperture radar, the Global Positioning System, and strain meters. Parameters can be constrained or free. Estimates can be calculated for single or multiple faults. Estimates of parameters are accompanied by reports of their covariances and uncertainties. Simplex has been tested extensively against forward models and against other means of inverting geodetic data and seismic observations. This work

Lyzenga, Gregory↗
DOE OSTI2021

Feasibility of DEIM for retrieving the initial field via dimensionality reduction

When parameter estimation is solved in a high-dimensional space, the dimensionality reduction strategy becomes the primary consideration for alleviating the tremendous computational cost. Here, the discrete empirical interpolation method (DEIM) is explored to retrieve the initial condition (IC) by combining the polynomial chaos (PC) based ensemble Kalman filter (i.e. PC-EnKF), where a non-intrusive PC expansion is considered as a surrogate model in place of the forward model in the prediction step of the ensemble Kalman filter, resulting in fewer forward model integrations but with a comparable accuracy as Monte Carlo-based approaches. The DEIM acts as a hyper-reduction tool to provide the low-dimensional input for the high-dimensional initial field, which can be reconstructed using the information on the sparse interpolation grid points that is adaptively obtained through PC-EnKF data assimilation method. Thus an innovative framework to reconstruct the IC is developed. The detailed procedure at each assimilation iteration includes: the determination of the spatial interpolation points, the estimation of the initial values on the interpolation locations using the optimal observations, and the reconstruction of IC in the full space. The current study uses the reconstruction field of initial conditions of the Navier-Stokes equations as an example to illustrate the efficacy of our method. The experimental results demonstrate the proposed algorithm achieves a satisfactory reconstruction for the initial field. The proposed method helps to extend the applicable area of DEIM in solving inverse problems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

A mechanistic understanding of oxygen isotopic changes in the Western United States at the Last Glacial Maximum

At the Last Glacial Maximum (LGM), records suggest drier conditions in the northwest United States and wetter conditions in the southwest United States relative to present-day as well as widespread changes in the isotopic composition of water. However, the mechanisms responsible for these changes remain ambiguous. Here, we explore differences in western United States hydroclimate between the LGM and preindustrial with a water isotope tracer enabled Earth System Model. We then use proxy forward models to compare simulated and recorded δ 18 O in speleothems. We find that the pattern of hydroclimate response in the western United States at the LGM relates to a combination of 1) increased frequency and southward shifted wintertime extratropical cyclones in the North Pacific, 2) greater rainout of moisture as it moves over the continent, and 3) reduced evaporation in the cooler LGM climate. The simulated lower δ 18 O of precipitation at the LGM relates predominantly to an increase in cool season moisture removal efficiency, with a secondary contribution from relatively more cool season precipitation. Both surface temperatures and North American ice sheets contribute to these hydroclimate changes at the LGM. Comparisons between δ 18 O from proxy forward models and speleothem records in the western United States show general agreement at the LGM, with increasing depletion moving towards the continental interior. Finally, this study highlights the similarities and differences between hydrologic and δ 18 O changes at the LGM and emphasizes the utility of model-proxy comparison for interpretation.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2016

Rough Surfaces: Is the Dark Stuff Just Shadow?: "Who knows what evil lurks in the hearts of men? The shadow knows!"

Remote observations of the surfaces of airless planetary objects are fundamental to inferring the physical structure and compositional makeup of the surface material. A number of forward models have been developed to reproduce the photometric behavior of these surfaces, based on specific, assumed structural properties such as macroscopic roughness and associated shadowing. Most work of this type is applied to geometric albedos, which are affected by complicated effects near zero phase angle that represent only a tiny fraction of the net energy reflected by the object. Other applications include parameter fits to resolved portions of some planetary surface as viewed over a range of geometries. The spherical albedo of the entire object (when it can be determined) captures the net energy balance of the particle more robustly than the geometric albedo. In most treatments involving spherical albedos, spherical albedos and particle phase functions are often treated as if they are independent, neglecting the effects of roughness. In this paper we take a different approach. We note that whatever function captures the phase angle dependence of the brightness of a realistic rough, shadowed, flat surface element relative to that of a smooth granular surface of the same material, it is manifested directly in both the integral phase function and the spherical albedo of the object. We suggest that, where broad phase angle coverage is possible, spherical albedos may be easily corrected for the effects of shadowing using observed (or assumed) phase functions, and then modeled more robustly using smooth-surface regolith radiative transfer models without further imposed (forward-modeled) shadowing corrections. Our approach attributes observed "power law" phase functions of various slope (and "linear" ranges of magnitude-vs.-phase angle) to shadowing, as have others, and goes on to suggest that regolith-model-based inferences of composition based on shadow-uncorrected spherical albedos overestimate the amount of absorbing material contained in the regolith.

geometric albedos↗
NASA NTRS1996

Research concerning the net flux of radiation in the atmosphere of Jupiter

The plan of the NFR (Net Flux of Radiation) team is for the data from the two solar channels (B and E) of NFR to be reduced with the goal of determining the solar heating rate. In order to determine the solar heating rate from the NFR measurements, effects due to the instrument's spatial and spectral response functions, to the temperature variation of the instrument (and associated drift of calibration), to the setting sun, and to the rotation of the probe (initially at a rate comparable to the NFR sampling frequency), all must be well modelled. In the past year, a forward modeling routine was created to simulate NFR data return in the B and E channels. The effects of varying parameters describing the atmospheric model (such as cloud location and thickness) and the descent profile (such as rotation rate) were investigated and an inversion routine was developed. For the forward modeling, existing radiative transfer codes were used to determine intensity fields within the Jovian atmosphere. A routine was developed to determine instantaneous instrument response by integrating the intensity field over the instrument response functions. A second routine was developed to determine the actual output of the NFR by integrating along an arbitrary descent trajectory. Near the top of the atmosphere, the upflux data alone are used to constrain the cloud structure of he atmosphere. To accomplish this, models are used to describe the variation in up flux between consecutive measurements in terms of variations of cloud opacity and variations in known parameters such as the solar zenith angle. This allows us to develop a zero-order model of cloud structure. Lower in the atmosphere, at levels where there is little or no azimuthal structure to the net flux measurements, both the up flux and net flux are used to derive layer transmission and reflection functions, which then determine layer opacity and single scattering albedo. A preliminary analysis of the data began in December 1995. In these data we could see the rapid oscillations expected at the beginning of the data due to probe rotation and the sun passing through the edge of the field of view. In addition, the time when this oscillation stopped was clearly visible. This sets the rough optical depth above the probe at this time.

Tomasko, M. G.↗
NASA NTRS1989

Comparison of measured induced velocities with results from a closed-form finite state wake model in forward flight

An unsteady, three-dimensional induced-flow model is applied to compute the induced-flow distribution of a rotor in forward flight, and numerical results are compared against LDV measurements for both time-averaged and unsteady induced flow at the disk. The former flow shows good agreement with measured data, except just behind the pylon at the lowest advanced ratio and near the blade tips for rectangular blades at high advance ratios. The method performs as well or better than other codes that have been applied to these data, takes less computing time, and is better suited to aeroelastic analysis. Results with only four harmonics and 15 state variables converge to all fundamental characteristics of the time-averaged flow.

Peters, David A.↗
NASA NTRS2019

Modeling Mg II During Solar Flares. I. Partial Frequency Redistribution, Opacity, and Coronal Irradiation

The Interface Region Imaging Spectrograph has routinely observed the flaring Mg II near-ultraviolet (NUV) spectrum, offering excellent diagnostic potential and a window into the location of energy deposition. A number of studies have forward-modeled both the general properties of these lines and specific flare observations. Generally these have forward-modeled radiation via post-processing of snapshots from hydrodynamic flare simulations through radiation transfer codes. There has, however, not been a study of how the physics included in these radiation transport codes affects the solution. A baseline setup for forward-modeling Mg II in flares is presented and contrasted with approaches that add or remove complexity. It is shown for Mg II that (1) partial frequency distribution (PRD) is still required during flare simulations despite the increased densities; (2) using full angle-dependent PRD affects the solution but takes significantly longer to process a snapshot; (3) including Mg I in non-LTE (NLTE) results in negligible differences to the Mg II lines but does affect the NUV quasi-continuum; (4) only hydrogen and Mg II need to be included in NLTE; (5) ideally the nonequilibrium hydrogen populations, with nonthermal collisional rates, should be used rather than the statistical equilibrium populations; (6) an atom consisting of only the ground state, h and k upper levels, and continuum level is insufficient to model the resonance lines; and (7) irradiation from a hot, dense flaring transition region can affect the formation of Mg II. We discuss modifications to the RH code allowing straightforward inclusion of the transition region and coronal irradiation in flares.

Kerr, Graham S.↗
DOE OSTI2024

Sequential Kalman tuning of the t -preconditioned Crank-Nicolson algorithm: efficient, adaptive and gradient-free inference for Bayesian inverse problems

Ensemble Kalman Inversion (EKI) has been proposed as an efficient method for the approximate solution of Bayesian inverse problems with expensive forward models. However, when applied to the Bayesian inverse problem EKI is only exact in the regime of Gaussian target measures and linear forward models. Here, in this work we propose embedding EKI and Flow Annealed Kalman Inversion, its normalizing flow (NF) preconditioned variant, within a Bayesian annealing scheme as part of an adaptive implementation of the t-preconditioned Crank-Nicolson (tpCN) sampler. The tpCN sampler differs from standard pCN in that its proposal is reversible with respect to the multivariate t-distribution. The more flexible tail behaviour allows for better adaptation to sampling from non-Gaussian targets. Within our Sequential Kalman Tuning (SKT) adaptation scheme, EKI is used to initialize and precondition the tpCN sampler for each annealed target. The subsequent tpCN iterations ensure particles are correctly distributed according to each annealed target, avoiding the accumulation of errors that would otherwise impact EKI. We demonstrate the performance of SKT for tpCN on three challenging numerical benchmarks, showing significant improvements in the rate of convergence compared to adaptation within standard SMC with importance weighted resampling at each temperature level, and compared to similar adaptive implementations of standard pCN. The SKT scheme applied to tpCN offers an efficient, practical solution for solving the Bayesian inverse problem when gradients of the forward model are not available. Code implementing the SKT schemes for tpCN is available at https://github.com/RichardGrumitt/KalmanMC.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1973

Acoustic characteristics of a large-scale augmentor wing model at forward speed

The augmentor wing concept is being studied as one means of attaining short takeoff and landing (STOL) performance in turbofan powered aircraft. Because of the stringent noise requirements for STOL operation, the acoustics of the augmentor wing are undergoing extensive research. The results of a wind tunnel investigation of a large-scale swept augmentor model at forward speed are presented. The augmentor was not acoustically treated, although the compressor supplying the high pressure primary air was treated to allow the measurement of only the augmentor noise. Installing the augmentor flap and shroud on the slot primary nozzle caused the acoustic dependence on jet velocity to change from eighth power to sixth power. Deflecting the augmentor at constant power increased the perceived noise level in the forward quadrant. The effect of airspeed was small. A small aft shift in perceived noise directivity was experienced with no significant change in sound power. Sealing the lower augmentor slot at a flap deflection of 70 deg reduced the perceived noise level in the aft quadrant. The seal prevented noise from propagating through the slot.

Falarski, M. D.↗
DOE OSTI2024

Differentiable stochastic halo occupation distribution

ABSTRACT In this work, we demonstrate how differentiable stochastic sampling techniques developed in the context of deep reinforcement learning can be used to perform efficient parameter inference over stochastic, simulation-based, forward models. As a particular example, we focus on the problem of estimating parameters of halo occupation distribution (HOD) models that are used to connect galaxies with their dark matter haloes. Using a combination of continuous relaxation and gradient re-parametrization techniques, we can obtain well-defined gradients with respect to HOD parameters through discrete galaxy catalogue realizations. Having access to these gradients allows us to leverage efficient sampling schemes, such as Hamiltonian Monte Carlo, and greatly speed up parameter inference. We demonstrate our technique on a mock galaxy catalogue generated from the Bolshoi simulation using a standard HOD model and find near-identical posteriors as standard Markov chain Monte Carlo techniques with an increase of ∼8× in convergence efficiency. Our differentiable HOD model also has broad applications in full forward model approaches to cosmic structure and cosmological analysis.

79 ASTRONOMY AND ASTROPHYSICS↗