Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “imaging simulators”

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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

High fidelity simulations of contaminant dispersion in an urban environment with comparison to magnetic resonance imaging measurements

The dispersion of a contaminant in an urban environment has the potential to impact a large population of people. In this work, a complex urban canopy flow based on the Oklahoma City downtown business district circa 2003 is studied using Magnetic Resonance Imaging (MRI) and high-fidelity Large Eddy Simulations (LES). MRI is a novel experimental technique that can provide high-resolution measurements in four dimensions (three spatial and temporal) for lab scale models. The experiments and simulations use the same geometry and boundary conditions providing a one-to-one comparison of the two methods. Results are presented on the time-averaged velocity and concentration fields, the temporal dynamics of the concentration plumes for a transient release, and a novel Cloud Identification Algorithm that can separate plumes produced by periodic contaminant releases used for ensemble averaging over many releases. The MRI and LES datasets both include millions of measurement voxels and the comparisons highlight the complex 3D nature of the flow including strong vertical velocities in spanwise street canyons and flow acceleration in streamwise street canyons. The concentration fields are qualitatively similar albeit the LES shows larger dispersion. A quantitative analysis with performance measures compares the datasets pointwise and demonstrates that the two 3D datasets are similar with respect to many measures including a fractional bias of 0.02 (ideal=0.0), correlation coefficient of 0.87 (ideal = 1.0), and the fraction points within a factor of 2 is 0.98 (ideal = 1.0). Plume analysis compares the arrival and residence time of contaminant and is found to vary significantly with location within the urban environment with arrival times between 0 and 1.25 and differences within the contaminant cloud less than 10% at most locations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Bayesian merging of numerical modeling and remote sensing for saltwater intrusion quantification in the Vietnamese Mekong Delta

Saltwater intrusion has become one of the most concerning issues in the Vietnamese Mekong Delta (VMD) due to its increasing impacts on agriculture and food security of Vietnam. Reliable estimation of salinity plays a crucial role to mitigate the impacts of saltwater intrusion. Here, this study developed a hybrid technique that merges satellite imagery with numerical simulations to improve the estimation of salinity in the VMD. The salinity derived from Landsat images and by numerical simulations was fused using the Bayesian inference technique. The results indicate that our technique significantly reduces the uncertainties and improves the accuracy of salinity estimates. The Nash–Sutcliffe coefficient is 0.74, which is much higher than that of numerical simulation (0.63) and Landsat estimation (0.6). The correlation coefficient between the ensemble and measured salinity is relatively high (0.88). The variance of the ensemble salinity errors (5.0 ppt 2 ) is lower than that of Landsat estimation (10.4 ppt 2 ) and numerical simulations (9.6 ppt 2 ). The proposed approach shows a great potential to combine multiple data sources of a variable of interest to improve its accuracy and reliability wherever these data are available.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Comparison of Polarized Radiative Transfer Codes Used by the EHT Collaboration

Interpretation of resolved polarized images of black holes by the Event Horizon Telescope (EHT) requires predictions of the polarized emission observable by an Earth-based instrument for a particular model of the black hole accretion system. Such predictions are generated by general relativistic radiative transfer (GRRT) codes, which integrate the equations of polarized radiative transfer in curved spacetime. A selection of ray-tracing GRRT codes used within the EHT Collaboration is evaluated for accuracy and consistency in producing a selection of test images, demonstrating that the various methods and implementations of radiative transfer calculations are highly consistent. When imaging an analytic accretion model, we find that all codes produce images similar within a pixel-wise normalized mean squared error (NMSE) of 0.012 in the worst case. When imaging a snapshot from a cell-based magnetohydrodynamic simulation, we find all test images to be similar within NMSEs of 0.02, 0.04, 0.04, and 0.12 in Stokes I , Q , U , and V , respectively. We additionally find the values of several image metrics relevant to published EHT results to be in agreement to much better precision than measurement uncertainties.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Validation of the imaging neutral particle analyzer in nearly MHD quiescent plasmas using injected beam ions on DIII-D

The imaging neutral particle analyzer (INPA) is a scintillator-based diagnostic that provides energy and radially resolved measurements of confined fast ions. To verify operation, understanding, and modeling of the diagnostic performance, distributions of deuterium beam ions in low density, nearly MHD-quiescent plasmas are measured on the DIII-D tokamak using the INPA. Here, the distribution is altered by changing the pitch-angle scattering rate with electron cyclotron heating (ECH). Experimental images are mostly reproduced with simulation of core fast ions charge exchanging with injected beam and edge cold neutrals in classical cases. Image comparisons between the experiment and simulation show agreement with less than 25% discrepancies near the beam injection energy of 50 keV. Other strong passive signals, observed in a localized region on the image, are suggested to come from charge exchange between cold neutrals and fast ions around the plasma boundary.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

HAIMOS Ensemble Forecasts for Intra-day and Day- Ahead GHI, DNI and Ramps

The objective of this research is to develop a hybrid physics-based/data-driven forecast model to improve direct normal and global horizontal irradiance (DNI and GHI) prediction for horizons ranging from 1 to 72 hours. Project objectives also address key gaps in state-of-the-art solar forecasting: accurate probabilistic solar forecasts and the forecasting of large irradiance ramps (ramp onset and magnitude). The proposed model ensembles Numerical Weather Prediction (NWP) forecasts, determinist physics-based algorithms, and new-generation cloud cover products (high-resolution rapid refresh satellite images and Large Eddy Simulations). The result is the Hybrid Adaptive Input Model Objective Selection (HAIMOS) ensemble model. HAIMOS blends state of the art machine learning methodologies with physics-based models for cloud cover and cloud optical depth forecasts. The technical activities followed a two-pronged strategy. First, the preprocessing of data, the selection of inputs to the nonlinear approximators, the type of approximator and objective functions, and post-processing ensembling techniques included in HAIMOS were all optimized adaptively to find the best model for a specific goal (reduce DNI/GHI forecast error, improve the prediction of ramp onset, etc.). Second, a large effort was put in improving cloud identification and the forecast of cloud cover and cloud optical depth. To this end, new-generation cloud parametrization products were developed in this work. These include improved algorithms to assist in cloud identification, cloud classification and cloud parametrization from satellite images – three key factors in the accuracy of 1 to 6-hours irradiance forecasts and prediction of ramp onset. Furthermore, we also included cloud information extracted from high resolution rapid refresh satellite images (GOES-16) and Large Eddy Simulations (LES). LES was used to model the atmosphere in detail over locations of interest and produce cloud optical depth forecasts. Once these data streams were validated, they were used as input data to the HAIMOS forecast. The model was developed using data from several climatologically distinct locations with potential for high solar penetration. In the last year of the project, we conducted a validation campaign according to the guidelines stipulated by the Topic Area 1 project as described in the FOA. This effort brings, for the first time, proven machine-learning methodologies for generating state-of-the-art solar forecasts interweaved with detailed physics-based models for cloud detection, and cloud optical depth forecasts. HAIMOS will generate accurate irradiance probabilistic forecast to assist in reducing solar generation prediction error. Globally optimized solar forecast models are more likely to impact solar energy stakeholders. The goal of this project was to increase the state-of-the-art forecast skill from their present values of 10 to 35%. At the end of the project, we achieved between 30% and 50% forecast skill across a wide range of horizons for both GHI and DNI.

14 SOLAR ENERGY↗
DOE OSTI2024

Analytical noise bias correction for precise weak lensing shear inference

Noise bias is a significant source of systematic error in weak gravitational lensing measurements that must be corrected to satisfy the stringent standards of modern imaging surveys in the era of precision cosmology. This paper reviews the analytical noise bias correction method and provides analytical derivations demonstrating that we can recover shear to its second order using the ‘renoising’ noise bias correction approach introduced by METACALIBRATION. We implement this analytical noise bias correction within the AnaCal shear estimation framework and propose several enhancements to the noise bias correction algorithm. We evaluate the improved AnaCal using simulations designed to replicate Rubin Legacy Survey of Space and Time (LSST) imaging data. These simulations feature semi-realistic galaxies and stars, complete with representative distributions of magnitudes and Galactic spatial density. We conduct tests under various observational challenges, including cosmic rays, defective CCD columns, bright star saturation, bleed trails, and spatially variable point spread functions. Our results indicate a multiplicative bias in weak lensing shear recovery of less than a few tenths of a per cent, meeting LSST Dark Energy Science Collaboration requirements without requiring calibration from external image simulations. Additionally, our algorithm achieves rapid processing, handling one galaxy in less than a millisecond.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

X-ray free-electron laser based dark-field X-ray microscopy: a simulation-based study

Dark-field X-ray microscopy (DFXM) is a nondestructive full-field imaging technique providing three-dimensional mapping of microstructure and local strain fields in deeply embedded crystalline elements. This is achieved by placing an objective lens in the diffracted beam, giving a magnified projection image. So far, the method has been applied with a time resolution of milliseconds to hours. In this work, the feasibility of DFXM at the picosecond time scale using an X-ray free-electron laser source and a pump–probe scheme is considered. Thermomechanical strain-wave simulations are combined with geometrical optics and wavefront propagation optics to simulate DFXM images of phonon dynamics in a diamond single crystal. Using the specifications of the XCS instrument at the Linac Coherent Light Source as an example results in simulated DFXM images clearly showing the propagation of a strain wave.

47 OTHER INSTRUMENTATION↗
DOE OSTI2022

Snowmass2021 Computational Frontier White Paper: Cosmological Simulations and Modeling

Powerful new observational facilities will come online over the next decade, enabling a number of discovery opportunities in the "Cosmic Frontier", which targets understanding of the physics of the early universe, dark matter and dark energy, and cosmological probes of fundamental physics, such as neutrino masses and modifications of Einstein gravity. Synergies between different experiments will be leveraged to present new classes of cosmic probes as well as to minimize systematic biases present in individual surveys. Success of this observational program requires actively pairing it with a well-matched state-of-the-art simulation and modeling effort. Next-generation cosmological modeling will increasingly focus on physically rich simulations able to model outputs of sky surveys spanning multiple wavebands. These simulations will have unprecedented resolution, volume coverage, and must deliver guaranteed high-fidelity results for individual surveys as well as for the cross-correlations across different surveys. The needed advances are as follows: (1) Development of scientifically rich and broadly-scoped simulations, which capture the relevant physics and correlations between probes (2) Accurate translation of simulation results into realistic image or spectral data to be directly compared with observations (3) Improved emulators and/or data-driven methods serving as surrogates for expensive simulations, constructed from a finite set of full-physics simulations (4) Detailed and transparent verification and validation programs for both simulations and analysis tools. (Abridged)

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2024

Using AI for Wave-front Estimation with the Rubin Observatory Active Optics System

Abstract The Vera C. Rubin Observatory will, over a period of 10 yr, repeatedly survey the southern sky. To ensure that images generated by Rubin meet the quality requirements for precision science, the observatory will use an active-optics system (AOS) to correct for alignment and mirror surface perturbations introduced by gravity and temperature gradients in the optical system. To accomplish this, Rubin will use out-of-focus images from sensors located at the edge of the focal plane to learn and correct for perturbations to the wave front. We have designed and integrated a deep-learning (DL) model for wave-front estimation into the AOS pipeline. In this paper, we compare the performance of this DL approach to Rubin’s baseline algorithm when applied to images from two different simulations of the Rubin optical system. We show the DL approach is faster and more accurate, achieving the atmospheric error floor both for high-quality images and low-quality images with heavy blending and vignetting. Compared to the baseline algorithm, the DL model is 40× faster, the median error 2× better under ideal conditions, 5× better in the presence of vignetting by the Rubin camera, and 14× better in the presence of blending in crowded fields. In addition, the DL model surpasses the required optical quality in simulations of the AOS closed loop. This system promises to increase the survey area useful for precision science by up to 8%. We discuss how this system might be deployed when commissioning and operating Rubin.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Sample environment effects on synchrotron-measured temperature profiles in an approximant of optical floating zone crystal growth

Even though the growth of crystals using optical floating zone furnaces has had an immense scientific impact, the implementation of this method remains more of an art than a science due to the difficulty of obtaining quantitative information about the sample thermal profile during crystal growth. Building on recent work demonstrating that in situ synchrotron studies can be used to map sample rod temperatures during heating, investigations were carried out to better understand how the sample environment affects the sample temperature profile. Through a combination of experimental studies and modeling efforts, it is shown that the environment in the furnace can strongly influence the sample temperature at the lamp focus, the steepness of the vertical temperature gradient, and the timescale required for sample heating and cooling - effects which can combine to produce a strong history-dependence to sample temperature profiles. It is demonstrated that the furnace effects can be effectively captured in thermal models, allowing both the steady-state and time-dependent behavior of the sample to be accurately reproduced with predictive models, and providing a launching point for improved furnace designs that can more readily deliver desired thermal profiles.

36 MATERIALS SCIENCE↗
DOE OSTI2026

High-resolution bandpass x-ray imaging with crystal reflectors: Overcoming geometric aberrations

The imaging problem of a specular reflector is revisited. Retaining terms through the second order in the reflector surface expansion, we derive the form of the aberration-limiting aperture for arbitrary magnification, assuming no bandwidth limitations. A permissible relative aperture size of the reflector is limited by a set relative aberration tolerance and scales with the tangent of the central glancing angle of incidence. These limiting aberrations become practically insignificant near backscattering. The results extend to x-ray diffracting crystals in symmetric Bragg geometry shaped as an ellipsoid of revolution. This geometry permits polychromatic imaging for hard x-rays over a bandwidth defined by the accepted range of Bragg angles, thereby suppressing aberrations of higher orders. We assess ellipsoidal crystal imagers using ray tracing simulations for two high-magnification designs with Bragg angles far from and close to backscattering. Finally, in both cases, the ellipsoidal crystals produce images of higher quality compared to those formed by equivalent toroidal crystal imagers.

Bragg reflection↗
DOE OSTI2021

Going Forward with the Nancy Grace Roman Space Telescope Transient Survey: Validation of Precision Forward-modeling Photometry for Undersampled Imaging

The Nancy Grace Roman Space Telescope (Roman) is an observatory for both wide-field observations and coronagraphy that is scheduled for launch in the mid-2020s. Part of the planned survey is a deep, cadenced field or fields that enable cosmological measurements with type Ia supernovae (SNe Ia). With a pixel scale of 0$_.^"$11, the Wide Field Instrument will be undersampled, presenting a difficulty for precisely subtracting the galaxy light underneath the SNe. We use simulated data to validate the ability of a forward-model code (such codes are frequently also called “scene-modeling” codes) to perform precision supernova photometry for the Roman SN survey. Our simulation includes over 760,000 image cutouts around SNe Ia or host galaxies (~10% of a full-scale survey). To have a realistic 2D distribution of underlying galaxy light, we use the VELA simulated high-resolution images of galaxies. We run each set of cutouts through our forward-modeling code which automatically measures time-dependent SN fluxes. Given our assumed inputs of a perfect model of the instrument point-spread functions and calibration, we find biases at the millimagnitude level from this method in four red filters (Y106, J129, H158, and F184), easily meeting the 0.5% Roman inter-filter calibration requirement for a cutting-edge measurement of cosmological parameters using SNe Ia. Simulated data in the bluer Z087 filter shows larger ~ 2–3 mmg biases, also meeting this requirement, but with more room for improvement. Our forward-model code has been released on Zenodo.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Integrating Experiments and Simulations to Reveal Anisotropic Growth Mechanisms and Interfaces of a One-Dimensional Zeolite

Zeolites are nanoporous crystalline materials critical for diverse industrial applications, yet their growth mechanisms are poorly understood. Here, this study presents a novel integrated framework combining experimental synthesis, high-resolution imaging, coarse-grained molecular dynamics simulations, and computer vision to uncover the mechanisms of growth of SSZ-24, a 1D channel zeolite. We demonstrate how synthesis conditions, such as temperature and reactant concentration, govern crystal anisotropy and surface roughness with growth dynamics differing markedly by crystallographic orientation. Along the channels, growth involves minimal energy barriers and rapid nucleation, resulting in rough surfaces. In contrast, growth perpendicular to the channels requires cooperative molecular organization and is highly sensitive to thermodynamic and kinetic conditions, yielding smooth anisotropic surfaces under low driving forces. By simulating transmission electron microscopy (TEM) images, we bridge molecular-scale simulations with experimental observations, identifying distinct growth mechanisms along different crystal planes. This work offers molecular-level insights into zeolite crystallization, advancing the rational design of nanoporous materials. The integration of cross-disciplinary methodologies establishes a transformative framework for optimizing zeolite synthesis, with implications for broader classes of materials.

Bertolazzo, Andressa A. [Univ. of Utah, Salt Lake ↗
DOE OSTI2025

Multimode turbulent flow measurements using magnetic resonance imaging- and laser-based techniques and computational fluid dynamics simulations

We studied the flow field characteristics of a turbulent flow over a regularized cube array with a perpendicular injection flow through the floor between the second and third cubical elements, representing the complex flow interactions of a 3D jet and the wake flows behind cubical obstacles. Four different experimental measurements were performed: two magnetic resonance imaging-based measurements for three-dimensional three-component velocity (MRV) and concentration (MRC) and two laser-based techniques, particle image velocimetry (PIV) and planar laser-induced fluorescence (PLIF), for two-dimensional two-component velocity and concentration measurement, respectively. The mainstream Reynolds number is Re = 15 000⁠, based on the primary inlet velocity U m and channel height D H ⁠, whereas the injector Reynolds number is Re j = 3400⁠, based on the injector velocity U j and the injector's exit width D j ⁠. Numerical simulations were performed for the studied flow configuration of turbulent flow over a regularized cube array using Reynolds-averaged Navier–Stokes (RANS) and large-eddy simulation (LES) approaches. Results obtained from experimental measurements—including MRV, MRC, PIV, and PLIF—as well as RANS and LES simulations are discussed and compared along several horizontal and vertical planes of the studied configuration. In addition, 3D turbulent flow structures, such as leading-edge vortex, horseshoe vortex, and jet shear-layer vortex, and the isosurfaces of scalar concentration successfully revealed by MRV and MRC techniques were found to be in very good agreement with those 3D features extracted from RANS and LES simulations. In conclusion, the high-resolution experimental and numerical database obtained from this study could be useful for validation and verification of numerical codes.

Computational fluid dynamics↗
DOE OSTI2021

Seeding the explosion of a high-current-density conductor in a controlled manner through the addition of micron-scale surface defects

Inhomogeneities in a current-carrying conductor promote non-uniform heating and expansion through the complex feedback between current density, electrical resistivity, Ohmic heating, temperature, and hydrodynamics. Three-dimensional-magnetohydrodynamic (3D-MHD) simulations suggest that μm-scale resistive inclusions or voids seed local overheating and through hydrodynamic explosion generate continuously growing craters which become several times larger than the initial perturbation. The ejected mass is the genesis of an electrothermally driven plasma filament which develops at lower current than plasmas on uniform surfaces adjacent to the defect. This result suggests that 1D or even 2D treatments are largely inadequate for detailed prediction of plasma formation. To test computational predictions, z-pinch experiments driven to 1 MA studied ultra-high-purity aluminum rods which were then machined to include pairs of quasi-hemispherical voids or “engineered defects (ED)” on the current-carrying surface. ED are the dominant current-density perturbation and reproducibly drive local overheating which can be compared with 3D-MHD simulation. Data from high-resolution-gated imagers of visible surface emissions confirm many simulation predictions, including the surface topography of local overheating, and the propensity for neighboring ED to prematurely source plasmas which then connect to form a plasma filament. Results also provide conditional support of theory which suggests heating similarity; that is, heating is independent of ED size for geometrically scaled ED.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Exploring the Great Pyramid: Detector Technical Design Report with Stand-Alone Monte Carlo Simulations

Cosmic-ray muon imaging has been used to non-destructively examine the Pyramids of Khufu and Khafre on the Giza Plateau; the EGP project will continue this line of research by undertaking a full tomographic scan of the former and by doing so will increase the sensitivity of the technique by upwards of two orders of magnitude. For this purpose, a muon telescope using triangular (vernier) detector cells far outperforms one with rectangular cells, providing an angular and positional RMS improvement by a factor of 4 for the same cost per unit area. A refinement algorithm was developed to handle tracks that produce secondaries. The triangular detector yields a positional resolution for a back-projection to the center of the pyramid of less than 20 cm, ensuring that a feature large enough to be of significance will still be seen by the telescope. This method is shown to be able to accurately and precisely reconstruct tracks of muons that pass through the King s Chamber, Queen s Chamber, and Grand Gallery.

43 PARTICLE ACCELERATORS↗
DOE OSTI2026

Collision Tracking in OpenMC: Methods and Applications in Neutron Noise, Neutron Imaging, Time-of-Flight, and Multiplicity Counting

We present the development and application of a collision tracking feature within the OpenMC Monte Carlo particle transport code, designed for diverse applications such as neutron spectroscopy, scatter camera system, neutron noise, and multiplicity counting simulations. This feature enables the tracking of individual particle collisions, with potential applications in nuclear nonproliferation, reactor physics, and nuclear security. Additionally, the feature holds potential for the calibration of neutron detectors, specifically in converting light output into energy deposited within the detectors. The implementation consists of a set of filters—such as reaction type, energy, cell, and material—that constrain the set of collisions that are tracked, extensions to the Python API to enable simple input specification, and support for writing either OpenMC’s native HDF5-based format or the Monte Carlo particle list format. This feature was added to the official OpenMC release in version 0.15.3. In this work, the feature will be applied to showcase scenarios such as time-of-flight simulations, scatter-camera imaging for neutron source localization, neutron-noise analysis to extract integral kinetic parameters such as the prompt decay constant α, and multiplicity counting to estimate the mass of special nuclear materials. Ultimately, this feature aims to expand the application scope of open-source Monte Carlo particle transport codes such as OpenMC.

Monte Carlo code↗
DOE OSTI2024

A Strong Gravitational Lens Is Worth a Thousand Dark Matter Halos: Inference on Small-scale Structure Using Sequential Methods

Strong gravitational lenses are a singular probe of the Universe’s small-scale structure—they are sensitive to the gravitational effects of low-mass (<10 10 M ⊙ ) halos even without a luminous counterpart. Recent strong-lensing analyses of dark matter structure rely on simulation-based inference (SBI). Modern SBI methods, which leverage neural networks as density estimators, have shown promise in extracting the halo-population signal. However, it is unclear whether the constraints from these models are limited by the methodology or the data. In this study, we introduce an accelerator-optimized simulation pipeline that can generate lens images with realistic subhalo populations in milliseconds. Leveraging this simulator, we identify the main limitation of our fiducial SBI analysis: training set size. We then adopt a sequential neural posterior estimation (SNPE) approach, allowing us to refine the training distribution to align with the observed data. Using only one-fifth as many mock Hubble Space Telescope images, SNPE matches the constraints on the low-mass halo population produced by our best nonsequential model. Our experiments suggest that an over 3 order-of-magnitude increase in training set size and GPU hours would be required to achieve an equivalent result without sequential methods. While the full potential of the existing lens sample remains to be explored, the notable improvement in constraining power enabled by our sequential approach highlights that current constraints are limited primarily by methodology and not the data itself. Moreover, our results emphasize the need to treat training set generation and model optimization as interconnected stages of any cosmological analysis using SBI.

79 ASTRONOMY AND ASTROPHYSICS↗