SEARCH · Search NASA

Results for “matched field processing”

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 397 records · Page 22

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Understanding the Angular Momentum Evolution of T Tauri and Herbig Ae/Be Stars

We investigate a sample of six Herbig Ae/Be stars belonging to the Orion OB1 association, as well as 73 low-mass objects, members of the σ Orionis cluster, in order to explore the angular momentum evolution at early stages of evolution, and its possible connection with main-sequence Ap/Bp magnetic stars. Using FIES and HECTOCHELLE spectra, we obtain projected rotational velocities through two independent methods. Individual masses, radii, and ages are computed using evolutionary models, distance, and cluster extinction. Under the assumption that similar physical processes operate in both T Tauri and Herbig Ae/Be stars, we construct snapshots of the protostar’s rotation against mass during the first 10 Myr with the aid of a rotational model that includes a variable disk lifetime, changes in the stellar moment of inertia, a dipolar magnetic field with variable strength, and angular momentum loss through stellar winds powered by accretion. We use these snapshots, as well as the rotational data, to infer a plausible scenario for the angular momentum evolution. We find that magnetic field strengths of a few kilo-Gauss at 3 Myr are required to match the rotational velocities of both groups of stars. Models with masses between 2 M {sub ⊙} and 3 M {sub ⊙} display larger angular momentum values by a factor of ~3, in comparison to stars of similar spectral types on the main sequence. Even though some quantitative estimates on this dramatic decrease with age for Ap/Bp magnetic main-sequence stars are presented, the results obtained for the angular momentum evolution do not explain their low rotation rates.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2016

Properties of the open cluster Tombaugh 1 from high-resolution spectroscopy and uvbyCaHβ photometry

Open clusters can be the key to deepening our knowledge on various issues involving the structure and evolution of the Galactic disk and details of stellar evolution because a cluster’s properties are applicable to all its members. However, the number of open clusters with detailed analysis from high-resolution spectroscopy or precision photometry imposes severe limitations on studies of these objects. To expand the number of open clusters with well-defined chemical abundances and fundamental parameters, we investigate the poorly studied, anticenter open cluster Tombaugh 1. Using precision uvbyCaHβ photometry and high-resolution spectroscopy, we derive the cluster’s reddening, obtain photometric metallicity estimates, and, for the first time, present a detailed abundance analysis of 10 potential cluster stars (nine clump stars and one Cepheid). Using the radial position from the cluster center and multiple color indices, we have isolated a sample of unevolved, probable single-star members of Tombaugh 1. From 51 stars, the cluster reddening is found to be E(b−y) = 0.221 ± 0.006 or E(B−V) = 0.303 ± 0.008, where the errors refer to the internal standard errors of the mean. The weighted photometric metallicity from m{sub 1} and hk is [Fe/H] = −0.10 ± 0.02, while a match to the Victoria-Regina Strömgren isochrones leads to an age of 0.95 ± 0.10 Gyr and an apparent modulus of (m−M) = 13.10 ± 0.10. Radial velocities identify six giants as probable cluster members, and the elemental abundances of Fe, Na, Mg, Al, Si, Ca, Ti, Cr, Ni, Y, Ba, Ce, and Nd have been derived for both the cluster and the field stars. Tombaugh 1 appears to be a typical inner thin disk, intermediate-age open cluster of slightly subsolar metallicity, located just beyond the solar circle, with solar elemental abundance ratios except for the heavy s-process elements, which are a factor of two above solar. Its metallicity is consistent with a steep metallicity gradient in the galactocentric region between 9.5 and 12 kpc. Our study also shows that Cepheid XZ CMa is not a member of Tombaugh 1 and reveals that this Cepheid presents signs of barium enrichment, making it a probable binary star.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

Resolved Molecular Gas Observations of MaNGA Post-starbursts Reveal a Tumultuous Past

Post-starburst (PSB) galaxies have recently and rapidly quenched their star formation; thus, they are an important way to understand how galaxies transition from star-forming late types to quiescent early types. The recent discovery of large cold gas reservoirs in PSB galaxies calls into question the theory that galaxies must lose their gas to become quiescent. Optical Integral Field Spectroscopy (IFS) surveys have revealed two classes of PSB galaxies: central PSB (cPSB) galaxies with central quenching regions and ring PSB (rPSB) galaxies with quenching in their outskirts. We analyze a sample of 13 nearby (z < 0.1) PSB galaxies with spatially resolved optical IFS data from the Mapping Nearby Galaxies at Apache Point Observatory survey and matched resolution Atacama Large Millimeter/submillimeter Array observations of 12 CO(1–0). Disturbed stellar kinematics in 7/13 of our PSB galaxies and centrally concentrated molecular gas is consistent with a recent merger for most of our sample. In galaxies without merger evidence, alternate processes may funnel gas inward and suppress star formation, which may include outflows, stellar bars, and minor mergers or interactions. The star formation efficiencies of the PSB regions in nearly half our galaxies are suppressed while the gas fractions are consistent with star-forming galaxies. Active galactic nucleus (AGN) feedback may drive this stabilization, and we observe AGN-consistent emission in the centers of 5/13 galaxies. Finally, our cPSB and rPSB galaxies have similar properties except the ionized and molecular gas in cPSB galaxies is more disturbed. Overall, the molecular gas in our PSB galaxies tends to be compact and highly disturbed, resulting in concentrated gas reservoirs unable to form stars efficiently.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2024

Boundary of the Distribution of Solar Wind Proton Beta versus Temperature Anisotropy

Abstract The frequency distribution of solar wind protons, measured in the vicinity of Earth’s orbit, is customarily plotted in ( β ∥ , T ⊥ / T ∥ ) phase space. Here, T ⊥ / T ∥ is the ratio of perpendicular and parallel temperatures, and β ∥ = 8 π nT ∥ / B 2 is the ratio of parallel thermal energy to background magnetic field energy, the so-called “parallel beta,” with ⊥ and ∥ denoting directions with respect to the ambient magnetic field. Such a frequency distribution, plotted as a two-dimensional histogram, forms a peculiar rhombic shape defined with an outer boundary in the said phase space. Past studies reveal that the threshold conditions for temperature anisotropy–driven plasma instability partially account for the boundary on the high- β ∥ side. The low- β ∥ side remains largely unexplained despite some efforts. Work by Vafin et al. recently showed that certain contours of collisional relaxation frequency, ν pp , when parameterized by T ⊥ / T ∥ and β ∥ , could match the overall shape of the left-hand boundary, thus suggesting that the collisional relaxation process might be closely related to the formation of the left-hand boundary. The present paper extends the analysis by Vafin et al. and carries out the dynamical computation of the collisional relaxation process for an ensemble of initial proton states with varying degrees of anisotropic temperatures. The final states of the relaxed protons are shown to closely match the observed boundary to the left of the ( β ∥ , T ⊥ / T ∥ ) phase space. When coupled with a similar set of calculations for the ensemble in the collective instability regime, it is found that the combined collisional/collective effects provide the baseline explanation for the observation.

Astronomy & Astrophysics↗

NASA NTRS2001

High Altitude Measurements of Radiance at High Spectral and Spatial Resolution for SIMBIOS Sensor Calibration, Validation, and Intercomparisons

The successful combination of data from different ocean color sensors depends on the correct interpretation of signal from each of these sensors. Ideally, the sensor measured signals are calibrated to geophysical units of spectral radiance, and sensor artifacts are removed and corrected. The calibration process resamples the signal into a common radiometric data space so that subsequent ocean color algorithms that are applied to the data are based on physical processes and are inherently sensor independent. The objective of this project is to calibrate and validate the on-orbit radiometric characteristics of Sea-Viewing Wide Field-of-view Sensor (SeaWiFS) with underflights of NASA's calibrated Airborne Visible/Infrared Imaging Spectrometer (AVIRIS). This objective is feasible because AVIRIS measures the same spectral range as SeaWIFS at higher spectral resolution. In addition to satellite sensor underflights, the AVIRIS project has supported comparison and analysis of the radiometric calibration standards used for AVIRIS and SeaWIFS. To date, both the OCTS and SeaWIFS satellite sensors have been underflown by AVIRIS with matching spectral, spatial, geometric, radiometric, and temporal domains. The calibration and validation objective of this project is pursued for the following reasons: (1) Calibration is essential for the quantitative use of SeaWIFS and other SIMBIOS (Sensor Intercomparison and Merger for Biological and Interdisciplinary Oceanic Studies) sensor data; (2) Calibration in the laboratory of spaceborne sensors is challenging; (3) Satellite sensors are subjected aging on the ground and to trauma during launch; (4) The Earth orbit environment is significantly different than the laboratory calibration environment; (5) Through years of effort AVIRIS has been demonstrated to be well calibrated; and (6) AVIRIS can match the spectral and spatial observation characteristics near the top of the atmosphere at the time of SeaWIFS measurements.

Green, Robert O.↗

DOE OSTI2019

Importance of Ambipolar Electric Field in Driving Ion Loss From Mars: Results From a Multifluid MHD Model With the Electron Pressure Equation Included

The multi-fluid (MF) magnetohydrodynamic model of Mars is improved by solving an additional electron pressure equation. Through the electron pressure equation, the electron temperature is calculated based on the effects from various electron-related heating and cooling processes (e.g., photo-electron heating, electron-neutral collision, and electron-ion collision), and thus, the improved model can calculate the electron temperature and the electron pressure force terms self-consistently. Model results of a typical case using the MF with electron pressure equation included model are compared in detail to identical cases using the MF and multi-species models to identify the effect of the improved physics. Here, we find that when the electron pressure equation is included, the general interaction patterns are similar to those with no electron pressure equation. However, the MF with electron pressure equation included model predicts that the electron temperature is much larger than the ion temperature in the ionosphere, consistent with both Viking and Mars Atmosphere and Volatile EvolutioN (MAVEN) observations. Using our numerical model, we also examined in detail the relative importance of different forces in the plasma interaction region. All three models are also applied to a MAVEN event study using identical input conditions; overall, the improved model matches best with MAVEN observations. All of the simulation cases are examined in terms of the total ion loss, and the results show that the inclusion of the electron pressure equation increases the escape rates by 50–110% in total mass, depending on solar condition and strong crustal field orientation, clearly demonstrating the importance of the ambipolar electric field in facilitating ion escape.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2016

SIFT-based dense pixel tracking on 0.35 T cine-MR images acquired during image-guided radiation therapy with application to gating optimization

Purpose: To first demonstrate the viability of applying an image processing technique for tracking regions on low-contrast cine-MR images acquired during image-guided radiation therapy, and then outline a scheme that uses tracking data for optimizing gating results in a patient-specific manner. Methods: A first-generation MR-IGRT system—treating patients since January 2014—integrates a 0.35 T MR scanner into an annular gantry consisting of three independent Co-60 sources. Obtaining adequate frame rates for capturing relevant patient motion across large fields-of-view currently requires coarse in-plane spatial resolution. This study initially (1) investigate the feasibility of rapidly tracking dense pixel correspondences across single, sagittal plane images (with both moderate signal-to-noise and spatial resolution) using a matching objective for highly descriptive vectors called scale-invariant feature transform (SIFT) descriptors associated to all pixels that describe intensity gradients in local regions around each pixel. To more accurately track features, (2) harmonic analysis was then applied to all pixel trajectories within a region-of-interest across a short training period. In particular, the procedure adjusts the motion of outlying trajectories whose relative spectral power within a frequency bandwidth consistent with respiration (or another form of periodic motion) does not exceed a threshold value that is manually specified following the training period. To evaluate the tracking reliability after applying this correction, conventional metrics—including Dice similarity coefficients (DSCs), mean tracking errors (MTEs), and Hausdorff distances (HD)—were used to compare target segmentations obtained via tracking to manually delineated segmentations. Upon confirming the viability of this descriptor-based procedure for reliably tracking features, the study (3) outlines a scheme for optimizing gating parameters—including relative target position and a tolerable margin about this position—derived from a probability density function that is constructed using tracking results obtained just prior to treatment. Results: The feasibility of applying the matching objective for SIFT descriptors toward pixel-by-pixel tracking on cine-MR acquisitions was first retrospectively demonstrated for 19 treatments (spanning various sites). Both with and without motion correction based on harmonic analysis, sub-pixel MTEs were obtained. A mean DSC value spanning all patients of 0.916 ± 0.001 was obtained without motion correction, with DSC values exceeding 0.85 for all patients considered. While most patients show accurate tracking without motion correction, harmonic analysis does yield substantial gain in accuracy (defined using HDs) for three particularly challenging subjects. An application of tracking toward a gating optimization procedure was then demonstrated that should allow a physician to balance beam-on time and tissue sparing in a patient-specific manner by tuning several intuitive parameters. Conclusions: Tracking results show high fidelity in assessing intrafractional motion observed on cine-MR acquisitions. Incorporating harmonic analysis during a training period improves the robustness of the tracking for challenging targets. The concomitant gating optimization procedure should allow for physicians to quantitatively assess gating effectiveness quickly just prior to treatment in a patient-specific manner.

61 RADIATION PROTECTION AND DOSIMETRY↗

NASA NTRS2009

Near-UV Sources in the Hubble Ultra Deep Field: The Catalog

The catalog from the first high resolution U-band image of the Hubble Ultra Deep Field, taken with Hubble s Wide Field Planetary Camera 2 through the F300W filter, is presented. We detect 96 U-band objects and compare and combine this catalog with a Great Observatories Origins Deep Survey (GOODS) B-selected catalog that provides B, V, i, and z photometry, spectral types, and photometric redshifts. We have also obtained Far-Ultraviolet (FUV, 1614 Angstroms) data with Hubble s Advanced Camera for Surveys Solar Blind Channel (ACS/SBC) and with Galaxy Evolution Explorer (GALEX). We detected 31 sources with ACS/SBC, 28 with GALEX/FUV, and 45 with GALEX/NUV. The methods of observations, image processing, object identification, catalog preparation, and catalog matching are presented.

Gardner, Jonathan P.↗

NASA NTRS1989

IRAS observations of dust heating and energy balance in the Rho Ophiuchi dark cloud

The equilibrium process dust emission in the Rho Ophiuchi dark cloud is studied. The luminosity of the cloud is found to closely match the luminosity of the clouds's known embedded and external radiation sources. There is no evidence for a large population of undetected low-luminosity sources within the cloud and unknown external heating is also only a minor source of energy. Most of the cloud's luminosity is emitted in the mid-to-far-IR. Dust temperature maps indicate that the dust is not hot enough to heat the gas to observed temperatures. A simple cloud model with a radiation field composed of flux HD 147889, S1, and Sco OB2 associations predicts the observed IRAS 60 to 100 micron in-band flux ratios for a mean cloud density n(H2) = 1400. Flattened 12 and 25 micron observations show much extended emission in these bands, suggesting stochastic heating of very small grains or large molecules.

Greene, Thomas P.↗

NASA NTRS1991

A dynamic analysis of the radiation excitation from the activation of a current collecting system in space

Current collecting systems moving in the ionosphere will induce electromagnetic wave radiation. The commonly used static analysis is incapable of studying the situation when such systems undergo transient processes. A dynamic analysis has been developed, and the radiation excitation processes are studied. This dynamic analysis is applied to study the temporal wave radiation from the activation of current collecting systems in space. The global scale electrodynamic interactions between a space-station-like structure and the ionospheric plasma are studied. The temporal evolution and spatial propagation of the electric wave field after the activation are described. The wave excitations by tethered systems are also studied. The dependencies of the temporal Alfven wave and lower hybrid wave radiation on the activation time and the space system structure are discussed. It is shown that the characteristics of wave radiation are determined by the matching of two sets of characteristic frequencies, and a rapid change in the current collection can give rise to substantial transient radiation interference. The limitations of the static and linear analysis are examined, and the condition under which the static assumption is valid is obtained.

Wang, J.↗

DOE OSTI2022

High harmonic generation in mixed XUV and NIR fields at a free-electron laser

Abstract We present the results of an experiment investigating the generation of high-order harmonics by a femtosecond near-infrared (NIR) laser pulse in the presence of an extreme ultraviolet (XUV) field provided by a free-electron laser (FEL), a process referred to as XUV-assisted high-order harmonic generation (HHG). Our experimental findings show that the XUV field can lead to a small enhancement in the harmonic yield when the XUV and NIR pulses overlap in time, while a strong decrease of the HHG yield and a red shift of the HHG spectrum is observed when the XUV precedes the NIR pulse. The latter observations are in qualitative agreement with model calculations that consider the effect of a decreased number of neutral emitters but are at odds with the predicted effect of the correspondingly increased ionization fraction on the phase matching. Our study demonstrates the technical feasibility of XUV-assisted HHG experiments at FELs, which may provide new avenues to investigate correlation-driven electron dynamics as well as novel ways to study and control propagation effects and phase matching in HHG.

47 OTHER INSTRUMENTATION↗

DOE OSTI2010

Improved semileptonic form factor calculations in lattice QCD

We investigate the computational efficiency of two stochastic based alternatives to the sequential propagator method used in lattice QCD calculations of heavy-light semileptonic form factors. In the first method, we replace the sequential propagator, which couples the calculation of two of the three propagators required for the calculation, with a stochastic propagator so that the calculations of all three propagators are independent. This method is more flexible than the sequential propagator method but introduces stochastic noise. We study the noise to determine when this method becomes competitive with the sequential propagator method, and find that for any practical calculation it is competitive with or superior to the sequential propagator method. We also examine a second stochastic method, the so-called 'one-end trick', concluding it is relatively inefficient in this context. The investigation is carried out on two gauge field ensembles, using the nonperturbatively improved Wilson-Sheikholeslami-Wohlert action with N{sub f}=2 mass-degenerate sea quarks. The two ensembles have similar lattice spacings but different sea-quark masses. We use the first stochastic method to extract O(a)-improved, matched lattice results for the semileptonic form factors on the ensemble with lighter sea quarks, extracting f{sub +}(0).

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2020

N-jettiness beam functions at N3LO

We present the first complete calculation for the quark and gluon N -jettiness ( $$ {\mathcal{T}}_N $$ T N ) beam functions at next-to-next-to-next-to-leading order (N 3 LO) in perturbative QCD. Our calculation is based on an expansion of the differential Higgs boson and Drell-Yan production cross sections about their collinear limit. This method allows us to employ cutting edge techniques for the computation of cross sections to extract the universal building blocks in question. The class of functions appearing in the matching coefficents for all channels includes iterated integrals with non-rational kernels, thus going beyond the one of harmonic polylogarithms. Our results are a key step in extending the $$ {\mathcal{T}}_N $$ T N subtraction methods to N 3 LO, and to resum $$ {\mathcal{T}}_N $$ T N distributions at N 3 LL' accuracy both for quark as well as for gluon initiated processes.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2013

The Northwest Geysers EGS Demonstration Project, California. Pre-stimulation Modeling and Interpretation of the Stimulation

The Northwest Geysers Enhanced Geothermal System (EGS) demonstration project aims to create an EGS by directly and systematically injecting cool water at relatively low pressure into a known High Temperature (280–400 °C) Zone (HTZ) located under the conventional (240 °C) geothermal steam reservoir at The Geysers geothermal field in California. Here we report that , the results of coupled thermal, hydraulic, and mechanical (THM) analyses made using a model developed as part of the pre-stimulation phase of the EGS demonstration project is presented. The model simulations were conducted in order to investigate injection strategies and the resulting effects of cold-water injection upon the EGS system; in particular to predict the extent of the stimulation zone for a given injection schedule. The actual injection began on October 6, 2011, and in this paper a comparison of pre-stimulation model predictions with micro-earthquake (MEQ) monitoring data over the first few months of a one-year injection program is presented. The results show that, by using a calibrated THM model based on historic injection and MEQ data at a nearby well, the predicted extent of the stimulation zone (defined as a zone of high MEQ density around the injection well) compares well with observed seismicity. The modeling indicates that the MEQ events are related to shear reactivation of preexisting fractures, which is triggered by the combined effects of injection-induced cooling around the injection well and small changes in steam pressure as far as half a kilometer away from the injection well. Pressure-monitoring data at adjacent wells and satellite-based ground-surface deformation data were also used to validate and further calibrate reservoir-scale hydraulic and mechanical model properties. The pressure signature monitored from the start of the injection was particularly useful for a precise back-calculation of reservoir porosity. Ultimately, the first few months of reservoir pressure and surface deformation data were useful for estimating the reservoir-rock permeability and elastic modulus. Finally, although the extent of the calculated stimulation zone matches the field observations over the first few months of injection, the observed surface deformations and MEQ evolution showed more heterogeneous behavior as a result of more complex geology, including minor faults and fracture zones that are important for consideration in the analysis of energy production and the long-term evolution of the EGS system.

58 GEOSCIENCES↗

DOE OSTI2011

Phase-matched emission from an optically thin medium following one-photon pulse excitation: Energy considerations

Scully and coworkers [M. O. Scully, E. S. Fry, C. H. R. Oii, and K. Wodkiewicz, Phys. Rev. Lett. 96, 010501 (2006)] demonstrated that there is directional, phase-matched emission following the excitation of an ensemble of atoms by a single-photon pulse. While the phase-matched emission intensity is proportional to the the number of atoms, for optically thin samples the total energy emitted in the phase-matched direction is much less than that radiated in other directions. Moreover, even for optically thin samples, it is necessary to take into account effects related to cooperative decay if energy is to be conserved in the overall emission process. An analytic calculation is presented to show explicitly how cooperative decay reduces the incoherent emission and restores energy conservation in this low-density limit.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2018

An Improved Representation of Vegetative Filter Strips in SWAT

Abstract. Vegetative filter strips (VFS) are popular conservation practices installed at the edges of agricultural fields to reduce losses of pollutants from agricultural areas to receiving waterbodies. The recent interest in using VFS areas as multifunctional landscapes necessitates an improved depiction of VFS in simulation models. This study is aimed to enhance the physical representation of VFS in the Soil and Water Assessment Tool (SWAT) to improve the representation of ecohydrologic processes and land management practices in VFS areas. The proposed enhancement enables routing of water, sediment, and nutrients from the source area through the VFS area and makes the infiltrated water and nutrients available for filter crop uptake. The improvements are implemented in SWAT by modifying input files through Matlab scripts and by changing SWAT subroutines to enable routing. The model improvements are tested with three paired watershed studies with and without edge-of-field VFS in central Iowa. The improved model estimated 46% runoff reduction, 91% sediment reduction, 83% total phosphorus reduction, and 54% nitrate reduction with VFS, and these estimates closely matched the measured VFS reductions. The improved model was able to represent increased infiltration, soil moisture, and denitrification in the VFS area, along with accurately capturing crop yields in the source and VFS areas. Overall, the results indicate improved physical representation of VFS in SWAT. Keywords: Conservation measures, Flexible buffer strips, Multifunctional buffers, SWAT model, Vegetative filter strips.

Agriculture↗

NASA NTRSDate not supplied

Generating a 4D Global CH(4) Product by Assimilating TROPOMI column CH(4) in NASA’s GEOS GCM

Examination of temporal and spatial CH4 variability is crucial for better understanding the human and natural processes driving climate change and ultimately designing mitigation strategies. Here we present an analysis framework that uses NASA’s GEOS General Circulation Model (GCM) to construct a high-resolution, time varying picture of atmospheric CH4 consistent with measurements from a variety of platforms, both in situ and remotely sensed. The resulting time varying atmospheric CH4 product can (i) support interpretation of high-resolution point source detection approaches, (ii) provide reanalysis fields for CH4 and other greenhouse gases, and (iii) supply boundary conditions for regional models. Our approach starts with a set of CH4 emissions from various inventories that have been adjusted to match the global annual growth rate over recent decades. These emissions are transported by the GEOS GCM, which in turn is constrained by meteorology from NASA’s Modern-Era Retrospective analysis for Research and Applications, Version 2 (MERRA-2) product. The simulated atmospheric field is compared with CH4 measurements, such as those from the TROPOspheric Monitoring Instrument (TROPOMI), and adjustments calculated following a Bayesian protocol. The accuracy of the resultant optimal atmospheric CH4 field can be demonstrated by its improved agreement (compared to a direct simulation of the CH4 inventories) with a host of independent CH4 measurements, such as those from the Total Carbon Column Observation Network (TCCON) and in situ observations from surface and airborne platforms.

Nikolay V. Balashov↗

DOE OSTI2026

Latent space dynamics identification for interface tracking with application to shock-induced pore collapse

Capturing sharp, evolving interfaces remains a central challenge in reduced-order modeling, especially when data is limited and the system exhibits localized nonlinearities or discontinuities. Here, we propose LaSDI-IT (Latent Space Dynamics Identification for Interface Tracking), a data-driven framework that combines low-dimensional latent dynamics learning with explicit interface-aware encoding to enable accurate and efficient modeling of physical systems involving moving material boundaries. At the core of LaSDI-IT is a revised autoencoder architecture that jointly reconstructs the physical field and an indicator function representing material regions or phases, allowing the model to track complex interface evolution without requiring detailed physical models or mesh adaptation. The latent dynamics are learned through linear regression in the encoded space and generalized across parameter regimes using Gaussian process interpolation with greedy sampling. We demonstrate LaSDI-IT on the problem of shock-induced pore collapse in high explosives, a process characterized by sharp temperature gradients and dynamically deforming pore geometries. The method achieves relative prediction errors below 9% across the parameter space, accurately recovers key quantities of interest such as pore area and hot spot formation, and matches the performance of dense training with only half the data. This latent dynamics prediction was 10 6 times faster than the conventional high-fidelity simulation, proving its utility for multi-query applications. These results highlight LaSDI-IT as a general, data-efficient framework for modeling discontinuity-rich systems in computational physics, with potential applications in multiphase flows, fracture mechanics, and phase change problems.

Gaussian process↗