Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “predictive mapping”

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
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Bringing 2D Eclipse Mapping out of the Shadows with Leave-one-out Cross Validation

Abstract Eclipse mapping is a technique for inferring 2D brightness maps of transiting exoplanets from the shape of an eclipse light curve. With JWST’s unmatched precision, eclipse mapping is now possible for a large number of exoplanets. However, eclipse mapping has only been applied to two planets, and the nuances of fitting eclipse maps are not yet fully understood. Here, we use Leave-one-out Cross Validation (LOO-CV) to investigate eclipse mapping, with application to a JWST NIRISS/SOSS observation of the ultrahot Jupiter WASP-18b. LOO-CV is a technique that provides insight into the out-of-sample predictive power of models on a data-point-by-data-point basis. We show that constraints on planetary brightness patterns behave as expected, with large-scale variations driven by the phase-curve variation in the light curve and smaller-scale structures constrained by the eclipse ingress and egress. For WASP-18b we show that the need for higher model complexity (smaller-scale features) is driven exclusively by the shape of the eclipse ingress and egress. We use LOO-CV to investigate the relationship between planetary brightness map components when mapping under a positive-flux constraint to better understand the need for complex models. Finally, we use LOO-CV to understand the degeneracy between the competing “hot spot” and “plateau” brightness map models of WASP-18b, showing that the plateau model is driven by the ingress shape and the hot spot model is driven by the egress shape, but preference for neither model is due to outliers or unmodeled signals. Based on this analysis, we make recommendations for the use of LOO-CV in future eclipse-mapping studies.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Automated Classification of Plasma Regions Using 3D Particle Energy Distributions

We investigate the properties of the ion sky maps produced by the Dual Ion Spectrometers (DIS) from the Fast Plasma Investigation (FPI). We have trained a convolutional neural network classifier to predict four regions crossed by the Magnetospheric Multiscale Mission (MMS) on the dayside magnetosphere: solar wind, ion foreshock, magnetosheath, and magnetopause using solely DIS spectrograms. The accuracy of the classifier is urn:x-wiley:21699380:media:jgra56742:jgra56742-math-0001%. We use the classifier to detect mixed plasma regions, in particular to find the bow shock regions. A similar approach can be used to identify the magnetopause crossings and reveal regions prone to magnetic reconnection. Data processing through the trained classifier is fast and efficient and thus can be used for classification for the whole MMS database.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Cholesterol modulates membrane elasticity via unified biophysical laws

Cholesterol and lipid unsaturation underlie a balance of opposing forces that features prominently in adaptive cell responses to diet and environmental cues. These competing factors have resulted in contradictory observations of membrane elasticity across different measurement scales, requiring chemical specificity to explain incompatible structural and elastic effects. Here, we demonstrate that – unlike macroscopic observations – lipid membranes exhibit a unified elastic behavior in the mesoscopic regime between molecular and macroscopic dimensions. Using nuclear spin techniques and computational analysis, we find that mesoscopic bending moduli follow a universal dependence on the lipid packing density regardless of cholesterol content, lipid unsaturation, or temperature. Our observations reveal that compositional complexity can be explained by simple biophysical laws that directly map membrane elasticity to molecular packing associated with biological function, curvature transformations, and protein interactions. The obtained scaling laws closely align with theoretical predictions based on conformational chain entropy and elastic stress fields. These findings provide unique insights into the membrane design rules optimized by nature and unlock predictive capabilities for guiding the functional performance of lipid-based materials in synthetic biology and real-world applications.

Kumarage, Teshani [Virginia Polytechnic Inst. and ↗
DOE OSTI2024

Accuracy of kinetic equilibrium reconstruction of NSTX and NSTX-U plasmas and its impact on the transport and stability analysis

An accurate magnetohydrodynamic (MHD) equilibrium reconstruction is an essential starting point for stability and transport plasma analysis. Herein this work describes an approach for obtaining kinetic equilibrium reconstructions using the OMFIT framework, which has been applied for the first time to spherical tokamak data from NSTX and NSTX-U. The EFIT equilibrium solver is integrated with experimental data analysis procedures and subsequent TRANSP transport simulations to enhance the accuracy of the reconstruction, in particular, at the edge region, by adding constraints on the total pressure and current density profiles, based on the transport code solution. The accuracy of the equilibrium reconstruction depends on the uncertainty and number of constraints, as well as the choice of basis functions to represent the pressure and current density profiles. Improved fidelity of the equilibrium reconstruction is demonstrated by reducing the variability of the magnetic axis and boundary locations from several centimeters, for reconstructions based on magnetic and experimental pressure constraints, to only several millimeters, for kinetic reconstructions based on transport code constraints, when different representations of basis functions were tested. The variability of the safety factor on axis was reduced ten times in the same sensitivity study. The accuracy of the equilibrium reconstruction and subsequent mapping of the experimental kinetic profile data have a significant impact on the trapped gyro Landau fluid and linear CGYRO turbulence simulations, which predict different spectra of unstable modes and turbulent fluxes for cases with different numbers of constraints in the equilibrium reconstruction. Conversely, the stability analysis performed using the GATO code shows plasmas that are stable to n = 1 MHD modes in both equilibria using magnetic and experimental pressure constraints as well as the transport code constrained equilibrium. However, a scan of parameters away from these conditions shows considerable deviation in the threshold of unstable modes between these reconstructions. Therefore, for reliable plasma analysis and use in turbulence and stability calculations, a high-fidelity equilibrium reconstruction with accurate kinetic constraints based on transport code solutions is necessary.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

HumoNet: A Framework for Realistic Modeling and Simulation of Human Mobility Network

Understanding, analyzing, and predicting human mobility and dynamics are valuable to solving pressing problems, developing effective plans, and prescribing timely remedies. As a computational approach, realistic human mobility simulations allow us to understand, analyze, and predict complex systems, including human societies. Accurate simulations rely on (1) the model that captures interactions and behaviors of myriad entities in our society and (2) the mapping of model instances to real-world entities. Taking this into account, this paper introduces the Human Mobility Network simulation framework (HumoNet), an integrated patterns of life (POL) simulation framework that leverages real-world data layers including transportation networks, points of interest, populations, popularity, and human trajectories. HumoNet is a data informed model in which agents are equipped with activities, locomotion, and planning capabilities. To simulate realistic kinematic maneuvers of individuals in transportation networks, HumoNet harnesses a microscopic traffic simulator that provides interaction among vehicles and traffic objects. In this paper, we describe the framework, outline our methodologies, and discuss the data processing and challenges of each data layer. Through experiments, we demonstrate that our simulations capture key features of human mobility by comparing them to the literature and real data using standard measures of human mobility (i.e., the radius of gyration, number of locations visited, level of exploration) and metrics scoring (i.e., Jensen-Shannon divergence). We envision that the synthetic data produced by HumoNet will serve as a benchmark for analyzing epidemics, deploying EV charging networks, and validating AI/ML tasks such as location prediction.

Kim, Joon-Seok↗
DOE OSTI2022

Asteroid Measurements at Millimeter Wavelengths with the South Pole Telescope

Abstract We present the first measurements of asteroids in millimeter wavelength data from the South Pole Telescope (SPT), which is used primarily to study the cosmic microwave background (CMB). We analyze maps of two ∼270 deg2 sky regions near the ecliptic plane, each observed with the SPTpol camera ∼100 times over 1 month. We subtract the mean of all maps of a given field, removing static sky signal, and then average the mean-subtracted maps at known asteroid locations. We detect three asteroids—(324) Bamberga, (13) Egeria, and (22) Kalliope—with signal-to-noise ratios (S/N) of 11.2, 10.4, and 6.1, respectively, at 2.0 mm (150 GHz); we also detect (324) Bamberga with an S/N of 4.1 at 3.2 mm (95 GHz). We place constraints on these asteroids’ effective emissivities, brightness temperatures, and light-curve modulation amplitude. Our flux density measurements of (324) Bamberga and (13) Egeria roughly agree with predictions, while our measurements of (22) Kalliope suggest lower flux, corresponding to effective emissivities of 0.64 ± 0.11 at 2.0 and < 0.47 at 3.2 mm. We predict the asteroids detectable in other SPT data sets and find good agreement with detections of (772) Tanete and (1093) Freda in recent data from the SPT-3G camera, which has ∼10× the mapping speed of SPTpol. This work is the first focused analysis of asteroids in data from CMB surveys, and it demonstrates we can repurpose historic and future data sets for asteroid studies. Future SPT measurements can help constrain the distribution of surface properties over a larger asteroid population.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Reliable extrapolation of deep neural operators informed by physics or sparse observations

Deep neural operators can learn nonlinear mappings between infinite-dimensional function spaces via deep neural networks. As promising surrogate solvers of partial differential equations (PDEs) for real-time prediction, deep neural operators such as deep operator networks (DeepONets) provide a new simulation paradigm in science and engineering. Pure data-driven neural operators and deep learning models, in general, are usually limited to interpolation scenarios, where new predictions utilize inputs within the support of the training set. However, in the inference stage of real-world applications, the input may lie outside the support, i.e., extrapolation is required, which may result to large errors and unavoidable failure of deep learning models. Here, we address this challenge of extrapolation for deep neural operators. First, we systematically investigate the extrapolation behavior of DeepONets by quantifying the extrapolation complexity, via the 2-Wasserstein distance between two function spaces and propose a new strategy of bias–variance trade-off for extrapolation with respect to model capacity. Subsequently, we develop a complete workflow, including extrapolation determination, and we propose five reliable learning methods that guarantee a safe prediction under extrapolation by requiring additional information—the governing PDEs of the system or sparse new observations. The proposed methods are based on either fine-tuning a pre-trained DeepONet or multifidelity learning. We demonstrate the effectiveness of the proposed framework for various types of parametric PDEs. Furthermore, our systematic comparisons provide practical guidelines for selecting a proper extrapolation method depending on the available information, desired accuracy, and required inference speed.

42 ENGINEERING↗
DOE OSTI2026

Phase diagram and spectroscopic signatures of a supersolid in the quantum ising magnet K 2 Co(SeO 3 ) 2

Supersolid phases are quantum-entangled states of matter exhibiting the dual characteristics of superfluidity and solidity. Theory predicts that hard-core bosons on a triangular lattice can form such phases at half filling and near complete filling. Leveraging an exact mapping between bosons and spin-$\frac{1}{2}$ degrees of freedom, here we show that these phases are realized in the triangular-lattice antiferromagnet K 2 Co(SeO 3 ) 2 . At zero field, neutron diffraction reveals the development of quasi-two-dimensional $\sqrt3$ x $\sqrt3$ magnetic order with Z 3 translational symmetry breaking (solidity), though with reduced amplitude indicating strong quantum fluctuations. These fluctuations manifest as equidistant bands of continuum neutron scattering, where the lowest-energy mode is gapless at K ($\frac{1}{3}$ $\frac{1}{3}$), consistent with broken U(1) spin rotational symmetry (superfluidity). For c-axis-oriented magnetic fields near saturation, we find a second phase consistent with a high-field supersolid. These two supersolids are separated by a pronounced 1/3 magnetization plateau phase that supports coherent spin waves, from which we determine the underlying spin Hamiltonian.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Reaction dynamics of S( 3 P) with 1,3-butadiene and isoprene: crossed-beam scattering, low-temperature flow experiments, and high-level electronic structure calculations

Sulfur atoms serve as key players in diverse chemical processes, from astrochemistry at very low temperature to combustion at high temperature. Building upon our prior findings, showing cyclization to thiophenes following the reaction of ground-state sulfur atoms with dienes, we here extend this investigation to include many additional reaction products, guided by detailed theoretical predictions. The outcomes highlight the complex formation of products during intersystem crossing (ISC) to the singlet surfaces. Here, we employed crossed-beam velocity map imaging and high-level ab initio methods to explore the reaction of S( 3 P) with 1,3-butadiene and isoprene under single-collision conditions and in low-temperature flows. For the butadiene reaction, our experimental results show the formation of thiophene via H 2 loss, a 2H-thiophenyl radical through H loss, and thioketene through ethene loss at a slightly higher collision energy compared to previous observations. Complementary Chirped-Pulse Fourier-Transform mmWave spectroscopy (CP-FTmmW) measurements in a uniform flow confirmed the formation of thioketene in the reaction at 20 K. For the isoprene reaction, we observed analogous products along with the 2H-thiophenyl radical arising from methyl loss and C 3 H 4 S (loss of ethene or H 2 + acetylene). CP-FTmmW detected the formation of thioformaldehyde via loss of 1,3-butadiene, again in the 20 K flow. Coupled-cluster calculations on the pathways found by the automated kinetic workflow code KinBot support these findings and indicate ISC to the singlet surface, leading to the generation of various long-lived intermediates, including 5-membered heterocycles.

Lang, Jinxin↗
DOE OSTI2022

Treating random sequential addition via the replica method

While many physical processes are non-equilibrium in nature, the theory and modeling of such phenomena lag behind theoretical treatments of equilibrium systems. The diversity of powerful theoretical tools available to describe equilibrium systems has inspired strategies that map non-equilibrium systems onto equivalent equilibrium analogs so that interrogation with standard statistical mechanical approaches is possible. In this work, we revisit the mapping from the non-equilibrium random sequential addition process onto an equilibrium multi-component mixture via the replica method, allowing for theoretical predictions of non-equilibrium structural quantities. We validate the above approach by comparing the theoretical predictions to numerical simulations of random sequential addition.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Modelling the Milky Way – I. Method and first results fitting the thick disc and halo with DES-Y3 data

ABSTRACT We present a technique to fit the stellar components of the Galaxy by comparing Hess Diagrams (HDs) generated from trilegal models to real data. We apply this technique, which we call mwfitting, to photometric data from the first 3 yr of the Dark Energy Survey (DES). After removing regions containing known resolved stellar systems such as globular clusters, dwarf galaxies, nearby galaxies, the Large Magellanic Cloud, and the Sagittarius Stream, our main sample spans a total area of ∼2300 deg2. We further explore a smaller subset (∼1300 deg2) that excludes all regions with known stellar streams and stellar overdensities. Validation tests on synthetic data possessing similar properties to the DES data show that the method is able to recover input parameters with a precision better than 3 per cent. We fit the DES data with an exponential thick disc model and an oblate double power-law halo model. We find that the best-fitting thick disc model has radial and vertical scale heights of 2.67 ± 0.09 kpc and 925 ± 40 pc, respectively. The stellar halo is fit with a broken power-law density profile with an oblateness of 0.75 ± 0.01, an inner index of 1.82 ± 0.08, an outer index of 4.14 ± 0.05, and a break at 18.52 ± 0.27 kpc from the Galactic centre. Several previously discovered stellar overdensities are recovered in the residual stellar density map, showing the reliability of mwfitting in determining the Galactic components. Simulations made with the best-fitting parameters are a promising way to predict Milky Way star counts for surveys such as the LSST and Euclid.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Covariance matrices for variance-suppressed simulations

ABSTRACT Cosmological N-body simulations provide numerical predictions of the structure of the Universe against which to compare data from ongoing and future surveys, but the growing volume of the Universe mapped by surveys requires correspondingly lower statistical uncertainties in simulations, usually achieved by increasing simulation sizes at the expense of computational power. It was recently proposed to reduce simulation variance without incurring additional computational costs by adopting fixed-amplitude initial conditions. This method has been demonstrated not to introduce bias in various statistics, including the two-point statistics of galaxy samples typically used for extracting cosmological parameters from galaxy redshift survey data, but requires us to revisit current methods for estimating covariance matrices of clustering statistics for simulations. In this work, we find that it is not trivial to construct covariance matrices analytically for fixed-amplitude simulations, but we demonstrate that ezmock (Effective Zel’dovich approximation mock catalogue), the most efficient method for constructing mock catalogues with accurate two- and three-point statistics, provides reasonable covariance matrix estimates for such simulations. We further examine how the variance suppression obtained by amplitude-fixing depends on three-point clustering, small-scale clustering, and galaxy bias, and propose intuitive explanations for the effects we observe based on the ezmock bias model.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Phase-field model of alloy solidification far from chemical equilibrium at the solid-liquid interface

We further develop a recently introduced phase-field model of far-from-equilibrium alloy solidification under additive manufacturing conditions [K. Ji et al., Phys. Rev. Lett. 130, 026203 (2023)]. This model utilizes enhanced solute diffusivity within the spatially diffuse interface region to quantitatively capture solute trapping with a larger interface width, thereby making simulations on experimentally relevant length and timescales computationally feasible. The main developments presented here include testing the robustness of different variational formulations, extending the model to concentrated alloys by incorporating solid and liquid free energies from thermodynamic databases, as illustrated for hypoeutectic Al-Ag alloys with CALPHAD, extending convergence tests as a function of interface width to 3D, and carrying out simulations in both 2D and 3D to examine existing theories of microstructure development. Our results indicate that the simplest variational formulation that interpolates the bulk free-energy density between its solid and liquid forms is the most robust. Remarkably, for hypoeutectic Al-Ag alloys, this formulation yields a high-velocity nonequilibrium phase diagram that is independent of interface width, thereby demonstrating that the framework of enhanced solute diffusivity can be nontrivially extended to concentrated alloys. Other variational formulations have restricted ranges of materials or processing parameters that can be reliably modeled. We use 2D simulations to construct high-velocity microstructure selection maps for dilute Al-Cu alloys. The results validate the important role of latent heat rejection at the interface and extend the limited predictions of linear stability analysis [A. Karma and A. Sarkissian, Phys. Rev. E 47, 513 (1993)] and sharp-interface 1D simulations to fully nonlinear regimes. Furthermore, 3D simulations, carried out using a computationally tractable axisymmetric cellular/dendritic interface shape, demonstrate a good convergence similar to that observed in 2D as a function of interface width. Full 3D simulations, in turn, reveal that the standard theory of absolute stability is a good predictor of the upper critical velocity beyond which steady-state growth becomes unstable, despite the different morphological manifestations of this instability in 2D and 3D.

36 MATERIALS SCIENCE↗
DOE OSTI2022

ENIGMA: Ecosystems and Networks Integrated with Genes and Molecular Assemblies. Final Report

ENIGMA research is performed at a very large scale involving collaborations between multiple research groups to address two major driving concepts: For any microbe, it is possible to discover efficiently the genetic determinants of adaptation to life in dynamic environments in contact with other living members of its ecology and reveal hidden functional phenomena that are only possible by using a systems biology approach; and for any environment, it is possible to predict the temporal changes in geochemistry and accompanying system behavior, for some length of time, given biogeochemical and physical inputs. These driving concepts also map directly to grand challenges the ENIGMA team has outlined: “Define the levels of biological organization most relevant to scaling from single cells to ecosystems and global cycles; capture how that organization varies in time and space; and identify critical interactions that dictate rates of carbon, nutrient, and energy transformation in different environments.”

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

QCD Equation of State and Phase Diagram from Holographic Black Holes

By using the AdS/CFT correspondence, we construct an Einstein-Maxwell-Dilaton model to map the thermodynamics of strongly interacting matter. The holographic model, constrained to reproduce the lattice QCD equation of state at zero baryon chemical potential, predicts a critical end point and a first order phase transition line. We also obtain the equation of state of the model for a large region of the phase diagram. We characterize the crossover transition region by two lines, one associated to the inflection of the second order baryon susceptibility, and another one associated to the minimum of the square of the speed of sound along trajectories of constant entropy per baryon number. We observe that both lines merge at the critical point.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

A Hydro-particle-mesh Code for Efficient and Rapid Simulations of the Intracluster Medium

We introduce the cosmological HYPER code based on an innovative hydro-particle-mesh (HPM) algorithm for efficient and rapid simulations of gas and dark matter. For the HPM algorithm, we update the approach of Gnedin & Hui to expand the scope of its application from the lower-density intergalactic medium (IGM) to the higher-density intracluster medium (ICM). While the original algorithm tracks only one effective particle species, the updated version separately tracks the gas and dark matter particles, as they do not exactly trace each other on small scales. For the approximate hydrodynamics solver, the pressure term in the gas equations of motion is calculated using robust physical models. In particular, we use a dark matter halo model, ICM pressure profile, and IGM temperature–density relation, all of which can be systematically varied for parameter-space studies. We show that the HYPER simulation results are in good agreement with the halo model expectations for the density, temperature, and pressure radial profiles. Simulated galaxy cluster scaling relations for Sunyaev–Zel’dovich (SZ) and X-ray observables are also in good agreement with mean predictions, with scatter comparable to that found in hydrodynamic simulations. HYPER also produces lightcone catalogs of dark matter halos and full-sky tomographic maps of the lensing convergence, SZ effect, and X-ray emission. These simulation products are useful for testing data analysis pipelines, generating training data for machine learning, understanding selection and systematic effects, and for interpreting astrophysical and cosmological constraints.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

High–Resolution Maps of Near–Surface Permafrost for Three Watersheds on the Seward Peninsula, Alaska Derived From Machine Learning

Permafrost soils are a critical component of the global carbon cycle and are locally important because they regulate the hydrologic flux from uplands to rivers. Furthermore, degradation of permafrost soils causes land surface subsidence, damaging infrastructure that is crucial for local communities. Regional and hemispherical maps of permafrost are too coarse to resolve distributions at a scale relevant to assessments of infrastructure stability or to illuminate geomorphic impacts of permafrost thaw. Here we train machine learning models to generate meter–scale maps of near–surface permafrost for three watersheds in the discontinuous permafrost region. The models were trained using ground truth determinations of near–surface permafrost presence from measurements of soil temperature and electrical resistivity. We trained three classifiers: extremely randomized trees (ERTr), support vector machines (SVM), and an artificial neural network (ANN). Model uncertainty was determined using k–fold cross validation, and the modeled extents of near–surface permafrost were compared to the observed extents at each site. At–a–site near–surface permafrost distributions predicted by the ERTr produced the highest accuracy (70%–90%). However, the transferability of the ERTr to the sites outside of the training data set was poor, with accuracies ranging from 50% to 77%. The SVM and ANN models had lower accuracies for at–a–site prediction (70%–83%), yet they had greater accuracy when transferred to the non–training site (62%–78%). These models demonstrate the potential for integrating high–resolution spatial data and machine learning models to develop maps of near–surface permafrost extent at resolutions fine enough to assess infrastructure vulnerability and landscape morphology influenced by permafrost thaw.

54 ENVIRONMENTAL SCIENCES↗