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 487 records · Page 27

Evaluation of Potential DAS Array Geometries for the Source Physics Experiment, Phase III

Distributed Acoustic Sensing (DAS) is a rapidly developing technology that can record acoustic wavefields at high sampling rates and with dense spatial spacing difficult to achieve with seismometers. However, the geophysical community has not fully explored DAS survey parameters to characterize the ideal array design. A better quantitative understanding of DAS array behavior prior to SPE Phase III acquisition can help improve the quality of the data recorded by guiding the DAS array design. Here we use array response functions as well as beamforming and backprojection results from forward modelling calculations to assess the performance of varying DAS array geometries to record regional and local sources. A seven-sided polygon DAS array demonstrated improved capabilities for recording regional sources over segmented linear arrays, with potential improvements in recording and locating local sources. These results help reveal DAS array performance as a function of geometry.

58 GEOSCIENCES↗

Millimeter-Wave Imaging Technology Development for Real-Time 2D and 3D Fusion Plasma Diagnostics

Including a one year No-Cost Extension, The University of California, Davis has been funded since 2020 to advance innovative millimeter-wave diagnostics for fusion plasma research. This work includes: (1) Developing and applying cutting-edge System-on-Chip millimeter wave Integrated Circuit (IC) technology for Electron Cyclotron Emission Imaging (ECEI) and Microwave Imaging Reflectometry (MIR) on the DIII-D tokamak and the NSTX-U spherical tokamak. (2) Development of radiation-hardened millimeter-wave chip for diagnostics in harsh environments, utilizing wide bandgap Gallium Nitride (GaN) semiconductors. (3) Integrating forward-modeling synthetic diagnostics to enable advanced, real-time plasma prediction and feedback control. Originally planned as a three-year program (2020–2023), the project was granted a one-year no-cost extension, extending its duration through the end of 2024.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Stochastic Ensemble Generation for Improved Characterization of Representing Geologic Variability in a Reservoir: IBDP Case Study for SMART Initiative

This document is a poster covering the findings from activities on training data generation, specifically geologic ensemble generation. The generated geologic realizations captured the range of possible permeability distributions of the subsurface at the Illinois Basin - Decatur Project (IBDP) site, based on available well log variabilities. The percentages of reservoirs and baffles in the injection zone and a truncation of baffle permeability led to more variance in the simulations. This will be used to build forward modeling, history matching, and optimization workflows. The geologic realizations were also ranked according to dynamic measures of hydraulic diffusivity, and simulations confirm a greater contrast between the reservoir and the baffles during injection.

stochastic ensemble generation↗

Deployment and Evaluation of SciStream on OLCF's Advanced Computing Ecosystem (ACE)

The growing demand for real-time analysis, experimental steering, and decision-making in scientific workflows has created a need for tightly coupled integrations between experimental facilities and high-performance computing (HPC) systems. The Department of Energy’s Integrated Research Infrastructure (IRI) initiative highlights data streaming as a key capability for enabling memory-to-memory data transfers, bypassing the limitations of traditional store-and-forward models. SciStream is a toolkit developed by researchers at Argonne National Laboratory (ANL) to support such streaming by addressing cross-domain security, delegated authentication, and application transparency. We deployed and evaluated SciStream on the Oak Ridge Leadership Computing Facility’s (OLCF) Advanced Computing Ecosystem (ACE) infrastructure, leveraging the Olivine OpenShift cluster and its high-bandwidth Data Streaming Nodes (DSNs) as gateway nodes. Our evaluation included synthetic streaming workloads derived from IRI science workflows, a streaming simulator, and integration with RabbitMQ to handle low-level messaging. This report documents the deployment process, performance evaluation, and challenges encountered, along with opportunities for future improvements.

97 MATHEMATICS AND COMPUTING↗

Performance of synthetic DAS as a function of array geometry

Distributed Acoustic Sensing (DAS) can record acoustic wavefields at high sampling rates and with dense spatial resolution difficult to achieve with seismometers. Using optical scattering induced by cable deformation, DAS can record strain fields with spatial resolution of a few meters. However, many experiments utilizing DAS have relied on unused, dark telecommunication fibers. As a result, the geophysical community has not fully explored DAS survey parameters to characterize the ideal array design. This limits our understanding of guiding principles in array design to deploy DAS effectively and efficiently in the field. A better quantitative understanding of DAS array behavior can improve the quality of the data recorded by guiding the DAS array design. Here we use steered response functions, which account for DAS fiber’s directional sensitivity, as well as beamforming and back-projection results from forward modelling calculations to assess the performance of varying DAS array geometries to record regional and local sources. A regular heptagon DAS array demonstrated improved capabilities for recording regional sources over other polygonal arrays, with potential improvements in recording and locating local sources. These results help reveal DAS array performance as a function of geometry and can guide future DAS deployments.

58 GEOSCIENCES↗

Machine Learning Predictions of Transition Probabilities in Atomic Spectra

Forward modeling of optical spectra with absolute radiometric intensities requires knowledge of the individual transition probabilities for every transition in the spectrum. In many cases, these transition probabilities, or Einstein A-coefficients, quickly become practically impossible to obtain through either theoretical or experimental methods. Complicated electronic orbitals with higher order effects will reduce the accuracy of theoretical models. Experimental measurements can be prohibitively expensive and are rarely comprehensive due to physical constraints and sheer volume of required measurements. Due to these limitations, spectral predictions for many element transitions are not attainable. In this work, we investigate the efficacy of using machine learning models, specifically fully connected neural networks (FCNN), to predict Einstein A-coefficients using data from the NIST Atomic Spectra Database. For simple elements where closed form quantum calculations are possible, the data-driven modeling workflow performs well but can still have lower precision than theoretical calculations. For more complicated nuclei, deep learning emerged more comparable to theoretical predictions, such as Hartree–Fock. Unlike experiment or theory, the deep learning approach scales favorably with the number of transitions in a spectrum, especially if the transition probabilities are distributed across a wide range of values. It is also capable of being trained on both theoretical and experimental values simultaneously. In addition, the model performance improves when training on multiple elements prior to testing. The scalability of the machine learning approach makes it a potentially promising technique for estimating transition probabilities in previously inaccessible regions of the spectral and thermal domains on a significantly reduced timeline.

74 ATOMIC AND MOLECULAR PHYSICS↗

DIAS: A Data-Informed Active Subspace Regularization Framework for Inverse Problems

This paper presents a regularization framework that aims to improve the fidelity of Tikhonov inverse solutions. At the heart of the framework is the data-informed regularization idea that only data-uninformed parameters need to be regularized, while the data-informed parameters, on which data and forward model are integrated, should remain untouched. We propose to employ the active subspace method to determine the data-informativeness of a parameter. The resulting framework is thus called a data-informed (DI) active subspace (DIAS) regularization. Four proposed DIAS variants are rigorously analyzed, shown to be robust with the regularization parameter and capable of avoiding polluting solution features informed by the data. They are thus well suited for problems with small or reasonably small noise corruptions in the data. Furthermore, the DIAS approaches can effectively reuse any Tikhonov regularization codes/libraries. Though they are readily applicable for nonlinear inverse problems, we focus on linear problems in this paper in order to gain insights into the framework. Various numerical results for linear inverse problems are presented to verify theoretical findings and to demonstrate advantages of the DIAS framework over the Tikhonov, truncated SVD, and the TSVD-based DI approaches.

Nguyen, Hai↗

Monitoring Offshore CO 2 Sequestration Using Marine CSEM Methods; Constraints Inferred from Field- and Laboratory-Based Gas Hydrate Studies

Offshore geological sequestration of CO 2 offers a viable approach for reducing greenhouse gas emissions into the atmosphere. Strategies include injection of CO 2 into the deep-ocean or ocean-floor sediments, whereby depending on pressure–temperature conditions, CO 2 can be trapped physically, gravitationally, or converted to CO 2 hydrate. Energy-driven research continues to also advance CO 2 -for-CH 4 replacement strategies in the gas hydrate stability zone (GHSZ), producing methane for natural gas needs while sequestering CO 2 . In all cases, safe storage of CO 2 requires reliable monitoring of the targeted CO 2 injection sites and the integrity of the repository over time, including possible leakage. Electromagnetic technologies used for oil and gas exploration, sensitive to electrical conductivity, have long been considered an optimal monitoring method, as CO 2 , similar to hydrocarbons, typically exhibits lower conductivity than the surrounding medium. We apply 3D controlled-source electromagnetic (CSEM) forward modeling code to simulate an evolving CO 2 reservoir in deep-ocean sediments, demonstrating sufficient sensitivity and resolution of CSEM data to detect reservoir changes even before sophisticated inversion of data. Laboratory measurements place further constraints on evaluating certain systems within the GHSZ; notably, CO 2 hydrate is measurably weaker than methane hydrate, and >1 order of magnitude more conductive, properties that may affect site selection, stability, and modeling considerations.

58 GEOSCIENCES↗

First Images of the Protoplanetary Disk around PDS 201

Scattered light imaging has revealed nearly a dozen circumstellar disks around young Herbig Ae/Be stars—enabling studies of structures in the upper disk layers as potential signs of ongoing planet formation. We present the first images of the disk around the variable Herbig Ae star PDS 201 (V* V351 Ori) and an analysis of the images and spectral energy distribution through 3D Monte Carlo radiative transfer simulations and forward modeling. The disk is detected in three data sets with the Large Binocular Telescope Interferometer/Large Binocular Telescope mid-infrared camera at the Large Binocular Telescope, including direct observations in the Ks and L′ filters, and an L′ observation with the 360° vector apodizing phase plate coronagraph. The scattered light disk extends to a very large radius of ∼250 au, which places it among the largest of such disks. Exterior to the disk, we establish detection limits on substellar companions down to ∼5 M {sub Jup} at ≳1.″5 (≳500 au), assuming the Baraffe et al. models. The images show a radial gap extending to ∼0.″4 (∼140 au at a distance of 340 pc) that is also evident in the spectral energy distribution. The large gap is a possible signpost of multiple high-mass giant planets at orbital distances (∼60–100 au) that are unusually massive and widely separated compared to those of planet populations previously inferred from protoplanetary disk substructures.

79 ASTRONOMY AND ASTROPHYSICS↗

Detection and Bulk Properties of the HR 8799 Planets with High-resolution Spectroscopy

Using the Keck Planet Imager and Characterizer, we obtained high-resolution (R ~ 35,000) K-band spectra of the four planets orbiting HR 8799. We clearly detected H{sub 2}O and CO in the atmospheres of HR 8799 c, d, and e, and tentatively detected a combination of CO and H{sub 2}O in b. These are the most challenging directly imaged exoplanets that have been observed at high spectral resolution to date when considering both their angular separations and flux ratios. We developed a forward-modeling framework that allows us to jointly fit the spectra of the planets and the diffracted starlight simultaneously in a likelihood-based approach and obtained posterior probabilities on their effective temperatures, surface gravities, radial velocities, and spins. We measured vsin (i) values of 10.1{sub -2.7}{sup +2.8} km s{sup -1} for HR 8799 d and 15.0{sub -2.6}{sup +2.3} km s{sup -1} for HR 8799 e, and placed an upper limit of <14 km s{sup -1} of HR 8799 c. Under two different assumptions of their obliquities, we found tentative evidence that rotation velocity is anticorrelated with companion mass, which could indicate that magnetic braking with a circumplanetary disk at early times is less efficient at spinning down lower-mass planets.

74 ATOMIC AND MOLECULAR PHYSICS↗

Angular Correlation Function Estimators Accounting for Contamination from Probabilistic Distance Measurements

With the advent of surveys containing millions to billions of galaxies, it is imperative to develop analysis techniques that utilize the available statistical power. In galaxy clustering, even small sample contamination arising from distance uncertainties can lead to large artifacts, which the standard estimator for two-point correlation functions does not account for. We first introduce a formalism, termed decontamination, that corrects for sample contamination by utilizing the observed cross-correlations in the contaminated samples; this corrects any correlation function estimator for contamination. Using this formalism, we present a new estimator that uses the standard estimator to measure correlation functions in the contaminated samples but then corrects for contamination. We also introduce a weighted estimator that assigns each galaxy a weight in each redshift bin based on its probability of being in that bin. We demonstrate that these estimators effectively recover the true correlation functions and their covariance matrices. Our estimators can correct for sample contamination caused by misclassification between object types as well as photometric redshifts; they should be particularly helpful for studies of galaxy evolution and baryonic acoustic oscillations, where forward modeling the clustering signal using the contaminated redshift distribution is undesirable.

79 ASTRONOMY AND ASTROPHYSICS↗

The HST See Change Program. I. Survey Design, Pipeline, and Supernova Discoveries*

The See Change survey was designed to make z > 1 cosmological measurements by efficiently discovering high-redshift Type Ia supernovae (SNe Ia) and improving cluster mass measurements through weak lensing. This survey observed twelve galaxy clusters with the Hubble Space Telescope (HST) spanning the redshift range z = 1.13–1.75, discovering 57 likely transients and 27 likely SNe Ia at z ~ 0.8–2.3. As in similar previous surveys, this proved to be a highly efficient use of HST for supernova observations; the See Change survey additionally tested the feasibility of maintaining, or further increasing, the efficiency at yet higher redshifts, where we have less detailed information on the expected cluster masses and star formation rates. Furthermore, we find that the resulting number of SNe Ia per orbit is a factor of ~8 higher than for a field search, and 45% of our orbits contained an active SN Ia within 22 rest-frame days of peak, with one of the clusters by itself yielding 6 of the SNe Ia. We present the survey design, pipeline, and supernova discoveries. Novel features include fully blinded supernova searches, the first random forest candidate classifier for undersampled IR data (with a 50% detection threshold within 0.05 mag of human searchers), real-time forward-modeling photometry of candidates, and semi-automated photometric classifications and follow-up forecasts. We also describe the spectroscopic follow-up, instrumental in measuring host galaxy redshifts. The cosmology analysis of our sample will be presented in a companion paper.

79 ASTRONOMY AND ASTROPHYSICS↗

AGNs on the Move: A Search for Off-nuclear AGNs from Recoiling Supermassive Black Holes and Ongoing Galaxy Mergers with the Zwicky Transient Facility

A supermassive black hole (SMBH) ejected from the potential well of its host galaxy via gravitational wave recoil carries important information about the mass ratio and spin alignment of the pre-merger SMBH binary. Such a recoiling SMBH may be detectable as an active galactic nucleus (AGN) broad-line region offset by up to 10 kpc from a disturbed host galaxy. We describe a novel methodology using forward modeling with The Tractor to search for such offset AGNs in a sample of 5493 optically variable AGNs detected with the Zwicky Transient Facility (ZTF). We present the discovery of nine AGNs that may be spatially offset from their host galaxies and are candidates for recoiling SMBHs. Five of these offset AGNs exhibit double-peaked broad Balmer lines, which may have arisen from unobscured accretion disk emission, and four show radio emission indicative of a relativistic jet. The fraction of double-peaked emitters in our spatially offset AGN sample is significantly larger than the 16% double-peaked emitter fraction observed for ZTF AGNs overall. In our sample of variable AGNs we also identified 52 merging galaxies, including a new spectroscopically confirmed dual AGN. Finally, we detected the dramatic rebrightening of SDSS 1133, a previously discovered variable object and recoiling SMBH candidate, in ZTF. The flare was accompanied by the reemergence of strong P Cygni line features, indicating that SDSS 1133 may be an outbursting luminous blue variable star.

79 ASTRONOMY AND ASTROPHYSICS↗

Variability-selected Intermediate-mass Black Hole Candidates in Dwarf Galaxies from ZTF and WISE

While it is difficult to observe the first black hole seeds in the early universe, we can study intermediate-mass black holes (IMBHs) in local dwarf galaxies for clues about their origins. In this paper we present a sample of variability-selected active galactic nuclei (AGN) in dwarf galaxies using optical photometry from the Zwicky Transient Facility (ZTF) and forward-modeled mid-IR photometry of time-resolved Wide-field Infrared Survey Explorer (WISE) co-added images. We found that 44 out of 25,714 dwarf galaxies had optically variable AGN candidates and 148 out of 79,879 dwarf galaxies had mid-IR variable AGN candidates, corresponding to active fractions of 0.17% ± 0.03% and 0.19% ± 0.02%, respectively. We found that spectroscopic approaches to AGN identification would have missed 81% of our ZTF IMBH candidates and 69% of our WISE IMBH candidates. Only nine candidates have been detected previously in radio, X-ray, and variability searches for dwarf galaxy AGN. The ZTF and WISE dwarf galaxy AGN with broad Balmer lines have virial masses of 10 5 M ⊙ < M BH < 10 7 M ⊙ , but for the rest of the sample, BH masses predicted from host galaxy mass range between 10 5.2 M ⊙ < M BH < 10 7.25 M ⊙ . We found that only 5 of 152 previously reported variability-selected AGN candidates from the Palomar Transient Factory in common with our parent sample were variable in ZTF. We also determined a nuclear supernova fraction of 0.05% ± 0.01% yr -1 for dwarf galaxies in ZTF. Our ZTF and WISE IMBH candidates show the promise of variability searches for the discovery of otherwise hidden low-mass AGN.

79 ASTRONOMY AND ASTROPHYSICS↗

BASS. XXXVI. Constraining the Local Supermassive Black Hole–Halo Connection with BASS DR2 AGNs

We investigate the connection between supermassive black holes (SMBHs) and their host dark matter halos in the local universe using the clustering statistics and luminosity function of active galactic nuclei (AGNs) from the Swift/BAT AGN Spectroscopic Survey (BASS DR2). By forward-modeling AGN activity into snapshot halo catalogs from N -body simulations, we test a scenario in which SMBH mass correlates with dark matter (sub)halo mass for fixed stellar mass. We compare this to a model absent of this correlation, where stellar mass alone determines the SMBH mass. We find that while both simple models are able to largely reproduce the abundance and overall clustering of AGNs, the model in which black hole mass is tightly correlated with halo mass is preferred by the data by 1.8 σ . When including an independent measurement on the black hole mass–halo mass correlation, this model is preferred by 4.6 σ . We show that the clustering trends with black hole mass can further break the degeneracies between the two scenarios and that our preferred model reproduces the measured clustering differences on one-halo scales between large and small black hole masses. These results indicate that the halo binding energy is fundamentally connected to the growth of SMBHs.

79 ASTRONOMY AND ASTROPHYSICS↗

A Machine Learning Framework to Predict Images of Edge-on Protoplanetary Disks

The physical structure and properties of protoplanetary disks are typically derived from spatially resolved disk images. Edge-on disks in particular provide an important view point on the vertical structure and degree of settling of disks. Such analyses rely on radiative transfer (RT) calculations that are generally computationally intensive due to the high optical depth of disks. Here we present a machine learning framework that has the potential to dramatically speed up the forward modeling process by approximating the results of RT calculations. This framework, trained on an initial set of RT calculations, utilizes an autoencoder neural network to enable the generation of synthetic scattered light images of edge-on disks directly from a set of physical parameters. We demonstrate that this framework generates synthetic images 2–3 orders of magnitude faster than using RT calculations. These machine learning-generated images appear to approximate the RT images well, in particular preserving their size and shape. We also find a strong correlation between the latent space representations of the generated disk images and several of their associated physical parameters. Finally, we discuss potential changes to the framework, such as methods to further improve the image quality, extending the framework to multiple wavelengths, and inverting the process to infer physical parameters from observed images. Overall, these new tools have the potential to enable a more efficient and uniform analysis of edge-on disk properties and the initial conditions of planet formation.

79 ASTRONOMY AND ASTROPHYSICS↗

Revisiting matrix-based inversion of scanning mobility particle sizer (SMPS) and humidified tandem differential mobility analyzer (HTDMA) data

Abstract. Tikhonov regularization is a tool for reducing noise amplification during data inversion. This work introduces RegularizationTools.jl, a general-purpose software package for applying Tikhonov regularization to data. The package implements well-established numerical algorithms and is suitable for systems of up to ∼ 1000 equations. Included is an abstraction to systematically categorize specific inversion configurations and their associated hyperparameters. A generic interface translates arbitrary linear forward models defined by a computer function into the corresponding design matrix. This obviates the need to explicitly write out and discretize the Fredholm integral equation, thus facilitating fast prototyping of new regularization schemes associated with measurement techniques. Example applications include the inversion involving data from scanning mobility particle sizers (SMPSs) and humidified tandem differential mobility analyzers (HTDMAs). Inversion of SMPS size distributions reported in this work builds upon the freely available software DifferentialMobilityAnalyzers.jl. The speed of inversion is improved by a factor of ∼ 200, now requiring between 2 and 5 ms per SMPS scan when using 120 size bins. Previously reported occasional failure to converge to a valid solution is reduced by switching from the L-curve method to generalized cross-validation as the metric to search for the optimal regularization parameter. Higher-order inversions resulting in smooth, denoised reconstructions of size distributions are now included in DifferentialMobilityAnalyzers.jl. This work also demonstrates that an SMPS-style matrix-based inversion can be applied to find the growth factor frequency distribution from raw HTDMA data while also accounting for multiply charged particles. The outcome of the aerosol-related inversion methods is showcased by inverting multi-week SMPS and HTDMA datasets from ground-based observations, including SMPS data obtained at Bodega Marine Laboratory during the CalWater 2/ACAPEX campaign and co-located SMPS and HTDMA data collected at the US Department of Energy observatory located at the Southern Great Plains site in Oklahoma, USA. Results show that the proposed approaches are suitable for unsupervised, nonparametric inversion of large-scale datasets as well as inversion in real time during data acquisition on low-cost reduced-instruction-set architectures used in single-board computers. The included software implementation of Tikhonov regularization is freely available, general, and domain-independent and thus can be applied to many other inverse problems arising in atmospheric measurement techniques and beyond.

54 ENVIRONMENTAL SCIENCES↗

Extending water vapor measurement capability of photon-limited differential absorption lidars through simultaneous denoising and inversion

Abstract. The micropulse differential absorption lidar (MPD) was developed at Montana State University (MSU) and the National Center for Atmospheric Research (NCAR) to perform range-resolved water vapor (WV) measurements using low-power lasers and photon-counting detectors. The MPD has proven to produce accurate WV measurements up to 6 km altitude. However, the MPD's ability to produce accurate higher-altitude WV measurements is impeded by the current standard differential absorption lidar (DIAL) retrieval methods. These methods are built upon a fundamental methodology that algebraically solves for the WV using the MPD forward models and noisy observations, which exacerbates any random noise in the lidar observations. The work in this paper introduces the adapted Poisson total variation (PTV) specifically for the MPD instrument. PTV was originally developed for a ground-based high spectral resolution lidar, and this paper reports on the adaptations that were required in order to apply PTV on MPD WV observations. The adapted PTV method, coined PTV-MPD, extends the maximum altitude of the MPD from 6 to 8 km and substantially increases the accuracy of the WV retrievals starting above 2 km. PTV-MPD achieves the improvement by simultaneously denoising the MPD noisy observations and inferring the WV by separating the random noise from the non-random WV. An analysis with 130 radiosonde (RS) comparisons shows that the relative root-mean-square difference (RRMSE) of WV measurements between RS and PTV-MPD exceeds 100 % between 6 and 8 km, whereas the RRMSE between RS and the standard method exceeds 100 % near 3 km. In addition, we show that by employing PTV-MPD, the MPD is able to extend its useful range of WV estimates beyond that of the ARM Southern Great Plains Raman lidar (RRMSE exceeding 100 % between 3 and 4 km); the Raman lidar has a power-aperture product 500 times greater than that of the MPD.

54 ENVIRONMENTAL SCIENCES↗