SEARCH · Engineering Papers

Results for “model data volumes”

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 919 records · Page 51

All sources NASA NTRS DOE OSTI

DOE OSTI2022

Application of the Redlich-Kister expansion for estimating the density of molten fluoride psuedo-ternary salt systems of nuclear industry interest

The development of next-generation molten salt reactors relies on accurate knowledge of the thermophysical properties of the candidate coolant and fueled molten salts. These thermophysical properties include density, viscosity, thermal conductivity, and heat capacity. Because of difficulties in measuring thermophysical properties of molten salts, there are many gaps in the current state of thermophysical property knowledge of these salts, particularly those that contain actinides or beryllium. Therefore, leveraging modeling techniques to estimate unknown molten salt thermophysical properties and guide future experimental measurements has high value for the nuclear industry. Here, the densities of molten fluoride pseudo-ternary salt systems, which are of interest to the nuclear industry, were estimated using Redlich-Kister expansion and Muggianu interpolation techniques. The pseudo-ternary systems considered for estimation in this study were NaF-LiF-ZrF 4 , LiF-BeF 2 -ZrF 4 , LiF-BeF 2 -ThF 4 , NaF-LiF-BeF 2 , NaF-KF-BeF 2 , NaF-ZrF 4 -UF 4 , and NaF-BeF 2 -UF 4 . This Redlich-Kister estimation approach accounts for nonideal mixing behavior based on pseudo-binary subsystem interaction parameters determined from experimentally measured pseudo-binary system density data sets. The Redlich-Kister estimation was compared with the method of additive molar volumes, which assumes ideal mixing. Additionally, the Redlich-Kister approach was used to determine previously unknown binary and ternary interaction parameters based on experimentally measured density data sets for select pseudo-ternary salt systems. The results of this study show improvement in density estimation using the Redlich-Kister approach for all systems considered compared with estimation by additive molar volumes. Furthermore, this analysis allowed for the estimation of nonideal density behavior in experimentally unstudied ZrF 4 -UF 4 and BeF 2 -UF 4 , as well as the quantification of ternary interaction in NaF-LiF-ZrF 4 , NaF-BeF 2 -UF 4 , and NaF-ZrF 4 -UF 4 .

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Determinants of cognitive and brain resilience to tau pathology: a longitudinal analysis

Mechanisms of resilience against tau pathology in individuals across the Alzheimer’s disease spectrum are insufficiently understood. Longitudinal data are necessary to reveal which factors relate to preserved cognition (i.e. cognitive resilience) and brain structure (i.e. brain resilience) despite abundant tau pathology, and to clarify whether these associations are cross-sectional or longitudinal. We used a longitudinal study design to investigate the role of several demographic, biological and brain structural factors in yielding cognitive and brain resilience to tau pathology as measured with PET. In this multicentre study, we included 366 amyloid-β-positive individuals with mild cognitive impairment or Alzheimer’s disease dementia with baseline 18 F-flortaucipir-PET and longitudinal cognitive assessments. A subset (n = 200) additionally underwent longitudinal structural MRI. We used linear mixed-effects models with global cognition and cortical thickness as dependent variables to investigate determinants of cognitive resilience and brain resilience, respectively. Models assessed whether age, sex, years of education, APOE-ε4 status, intracranial volume (and cortical thickness for cognitive resilience models) modified the association of tau pathology with cognitive decline or cortical thinning. We found that the association between higher baseline tau-PET levels (quantified in a temporal meta-region of interest) and rate of cognitive decline (measured with repeated Mini-Mental State Examination) was adversely modified by older age (Stβ interaction = -0.062, P = 0.032), higher education level (Stβ interaction = -0.072, P = 0.011) and higher intracranial volume (Stβ interaction = -0.07, P = 0.016). Younger age, higher education and greater cortical thickness were associated with better cognitive performance at baseline. Greater cortical thickness was furthermore associated with slower cognitive decline independent of tau burden. Higher education also modified the negative impact of tau-PET on cortical thinning, while older age was associated with higher baseline cortical thickness and slower rate of cortical thinning independent of tau. Our analyses revealed no (cross-sectional or longitudinal) associations for sex and APOE-ε4 status on cognition and cortical thickness. In this longitudinal study of clinically impaired individuals with underlying Alzheimer’s disease neuropathological changes, we identified education as the most robust determinant of both cognitive and brain resilience against tau pathology. The observed interaction with tau burden on cognitive decline suggests that education may be protective against cognitive decline and brain atrophy at lower levels of tau pathology, with a potential depletion of resilience resources with advancing pathology. Finally, we did not find major contributions of sex to brain nor cognitive resilience, suggesting that previous links between sex and resilience might be mainly driven by cross-sectional differences.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2022

Evidence of galaxy assembly bias in SDSS DR7 galaxy samples from count statistics

We present observational constraints on the galaxy–halo connection, focusing particularly on galaxy assembly bias from a novel combination of counts-in-cylinders statistics, P(N CIC ), with the standard measurements of the projected two-point correlation function w p (r p ), and number density n gal of galaxies. We measure n gal , w p (r p ), and P(N CIC ) for volume-limited, luminosity-threshold samples of galaxies selected from SDSS DR7, and use them to constrain halo occupation distribution (HOD) models, including a model in which galaxy occupation depends upon a secondary halo property, namely halo concentration. We detect significant positive central assembly bias for the M r < -20.0 and M r <-19.5 samples. Central galaxies preferentially reside within haloes of high concentration at fixed mass. Positive central assembly bias is also favoured in the M r < -20.5 and M r < -19.0 samples. We find no evidence of central assembly bias in the M r < -21.0 sample. We observe only a marginal preference for negative satellite assembly bias in the M r < -20.0 and M r < -19.0 samples, and non-zero satellite assembly bias is not indicated in other samples. Our findings underscore the necessity of accounting for galaxy assembly bias when interpreting galaxy survey data, and demonstrate the potential of count statistics in extracting information from the spatial distribution of galaxies, which could be applied to both galaxy–halo connection studies and cosmological analyses.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2024

Physics and chemistry from parsimonious representations: image analysis via invariant variational autoencoders

Electron, optical, and scanning probe microscopy methods are generating ever increasing volume of image data containing information on atomic and mesoscale structures and functionalities. This necessitates the development of the machine learning methods for discovery of physical and chemical phenomena from the data, such as manifestations of symmetry breaking phenomena in electron and scanning tunneling microscopy images, or variability of the nanoparticles. Variational autoencoders (VAEs) are emerging as a powerful paradigm for the unsupervised data analysis, allowing to disentangle the factors of variability and discover optimal parsimonious representation. Here, we summarize recent developments in VAEs, covering the basic principles and intuition behind the VAEs. The invariant VAEs are introduced as an approach to accommodate scale and translation invariances present in imaging data and separate known factors of variations from the ones to be discovered. We further describe the opportunities enabled by the control over VAE architecture, including conditional, semi-supervised, and joint VAEs. Several case studies of VAE applications for toy models and experimental datasets in Scanning Transmission Electron Microscopy are discussed, emphasizing the deep connection between VAE and basic physical principles. Python codes and datasets discussed in this article are available at https://github.com/saimani5/VAE-tutorials and can be used by researchers as an application guide when applying these to their own datasets.

36 MATERIALS SCIENCE↗

DOE OSTI2026

Low energy neutron light output characterization of EJ301D and deuterated stilbene with a comparison of light output characterization methods

The neutron-induced light yield of a 2.54 cm diameter by 2.54 cm long right circular cylinder of EJ301D and a (5.08 cm)3 custom made cube of deuterated trans-stilbene-d12 (d-stilbene) were measured over incident neutron energies from 300 keV to 2.2 MeV and 200 keV to 2.4 MeV, respectively. The measurements were performed using a time-of-flight experiment with a Cf-252 source and an approximately 1.5 m flight path. We compare three light output spectrum full energy deposition edge estimation methods: (1) simulating the neutron energy spectrum edge and fitting it to the light output spectrum, (2) using the inflection point of the light output spectrum edge (derivative method, a.k.a. Kornilov’s method), and (3) using an empirical model fit to the edge of the light output spectrum. Both the derivative and equation fit methods do not account for physical processes such as multiple neutron scattering in the detectors. They instead rely on assumptions about the linear shape continuum shape of the light output spectrum and the direct correlation between the location of the spectrum’s inflection point and maximum energy deposition. These assumptions were found to introduce bias into those methods when tested against simulated spectra with known edge locations. When tested against measured spectra the derivative method was found to differ from the simulation fit by greater than 30% at low energies with large discontinuities for adjacent data points above 800 keV incident neutron energy. The empirical equation fitting method was found to also exhibit bias of a similar magnitude, but with significantly more continuous behavior, especially with the lower count data of the smaller volumed EJ301D scintillator. Experimental light output yield for this neutron energy range is reported using the simulated spectrum fitting method because it includes physics neglected by the other methods, and did not exhibit the bias observed in the other methods

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2023

High-bandwidth density silicon photonic resonators for energy-efficient optical interconnects

The growth of artificial intelligence applications demands ever larger and more complex deep learning models, dominating today's—and tomorrow's—data center and high-performance computing systems. While traditional electronics are failing to keep pace with application demands, silicon photonic (SiPh) interconnects have emerged as a necessary technology to support these systems. SiPh-driven wavelength-division multiplexing (WDM) offers a particularly promising path toward supporting incredibly high-aggregate link bandwidth in a compact and efficient form factor. One of the basic building blocks of these integrated WDM interconnects is the SiPh resonator. Their inherent wavelength selectivity and compact footprint allow for efficient data transmission multiplexed across dozens of carrier wavelengths. Used as add-drop (AD) filters, SiPh resonators are critical to constructing integrated tunable wavelength-selective optical circuit switches as well as for demultiplexing the different carrier wavelengths toward independent wavelength-insensitive photodiodes in a dense wavelength-division multiplexing receiver. Resonators in the all-pass (AP) configuration are widespread as well, allowing for wavelength-selective modulation to drive aggregate link bandwidths far beyond the individual channel data rate. Unlike SiPh Mach–Zehnder modulators (MZM), resonant modulators can be driven using low, complementary metal-oxide-semiconductor drive voltages, allowing for tight co-integration between photonic integrated circuits, fabricated with larger process node technologies, and electronic integrated circuits, designed to exploit the advantages of the latest node. To push toward practical peta-scale interconnects, a comprehensive review of SiPh resonators is required, addressing bottlenecks and design constraints at both the architecture and device levels. We first describe the predominant integrated link architectures and identify their limits. We then discuss the device-level design considerations that can be made for both AD and AP configuration resonators to overcome the system level limits with novel resonator device designs. Analytical models and numerical simulation of resonators are validated by experimental measurement of devices fabricated in a commercial 300-mm foundry, showing a clear path toward volume manufacturing. The demonstrated resonant modulators and filters support the feasibility of increasing the aggregate bandwidth of resonator-driven SiPh interconnects into the peta-scale regime.

42 ENGINEERING↗

DOE OSTI2020

LLNL Kimberlina 1.2 NUFT Simulations June 2018 (v2)

This dataset contains the output 6,000, 3-dimensional reactive multi-phase flow and transport aquifer simulations of brine and CO2 leakage into a protective aquiver in California’s San Joaquin Valley and input data files detailing the geologic mesh, aquifer physical properties and CO2 and brine injection rates. This data set was generated as an ongoing effort with the US DOE National Risk Assessment Partnership (NRAP) to evaluate the effectiveness of monitoring techniques to detect brine and CO2 leakage from legacy wells into underground sources of drinking water overlaying a CO2 storage reservoir. Each simulation contains a unique set of input parameters, generated stochastically. The outputs consist of these upper three geologic layers (from top): the Etchegoin, Macoma-Chanac, Santa Margarita-McLure formations. These simulations span the several distances (1, 3 and 6 km or wells W31-0.2, W31-0.5 and W31-1.0, respectively) from the CO2 injector, initiated from bottom hole pressure and saturation to calculate wellbore leakage from the storage reservoir, with low and high regional groundwater gradients and wellbore leakage into 5 leaky nodes. The dataset includes 1,000 unique simulations for each distance, which each contain a unique aquifer heterogeneity, aquifer and caprock permeability, and two model generations are included with a high permeability (prod07) and hybrid permeability (prod09). The range of permeability distributions is listed in Table 1. Each model generation consists of 3,000 simulations. Included in the dataset are the leakage rates determined from 2D wellbore models which utilize the pressure and CO2 saturation from LBL's reservoir simulations, NUFT mesh files with distributed lithology, NUFT rocktab files which describe the material properties for the geologic layers and the NUFT input files and post-processed output 'ntab' files. Each ntab file contains spatial (rows) and temporal (columns) model output tables for each model cell, the locations (x,y,z) and dimensions for each cells (dx, dy, dz). Table 1. Permeability distribution ranges for prod07 and prod09 model generations Geologic Layer: Permeability Range (log10 m^2) prod07 prod09 Etchegoin -12.92 to -10.92 -13.70 to -11.44 Macoma-Chanac -12.72 to -10.72 -13.50 to -11.24 Santa Margarita-McLure -12.70 to -10.70 -13.48 to -11.22 The input files used to generate the model include which are included in the dataset are: Time series of CO2 leakage input into the model (ex: Q_brn.W31-0.2.sim1000.layers123.tab) Time series of CO2 leakage input into the model (ex: Q_CO2.W31-0.2.sim1000.layers123.tab) Physical properties of the aquifer materials detailing the aquifer porosity, solid density, partitioning coefficients, permeabilities and van-Genuchten parameters detailed in a NUFT rocktab file: (ex: sim1000.usnt.rocktab) Numerical mesh and geologic data assigned to each model cell detailed in a NUFT genmsh format (ex: sim1000.mesh_k16.prod07.trans.genmsh) The primary output parameters are: pH (use absolute value) Change in TDS (mg/kg) Change in Pressure (Pa) Change CO2 gas saturation (fraction range 0.0-1.0) for example, the directory /p/lscratchh/mansoor1/nrap/kimberlina/prod09/mainfiles/sim1000/W31- 0.2 contains: sim1000.W31-0.2.trans.pH.red.ntab sim1000.W31-0.2.no_bg.trans.TDS.red.ntab sim1000.W31-0.2.usnt.P.deltabg.red.ntab sim1000.W31-0.2.usnt.CO2_sat.deltabg.red.ntab Each row in the NTAB files consist of model output per numerical grid cell. Each output file contains 33 columns (variables), including the information of numerical records, geologic location and sizes and the simulated parameter values over time. The first 13 variables are about numerical records and relative geologic information for a simulation grid: 1. index: simulation index 2. i: the ith grid of x-axis 3. j: the ith grid of y-axis 4. k: the ith grid of z-axis 5. element_ref: element reference 6. nuft_ind: nuft index 7. x: grid location in the x axis direction 8. y: grid location in the y axis direction 9. z: grid location in the z axis direction 10. dx: grid length in the x axis direction 11. dy: grid length in the y axis direction 12. dz: grid length in the z axis direction 13. volume: volume of the simulation grid The remainder (14, 15, 16...) variables are the simulated parameter values over time, take Pressure as an example, are: 14. 0.0y: initial pressure per cell. 15. 10.0y: simulated pressure at the end of the 10th year. 16. 20.0y: simulated pressure at the end of the 20th year. ... (repeated for every 10 years until 200 years)... The model extends 10,000 m, 5,000 m and 1,411 m in the x,y and z dimensions, respectively. The mesh consists of 164,832 cells with mesh dimensions of 101 x 51 x 32 (nx, ny, nz), with cell dimensions ranging from 100 m laterally (along x and y-axis) and model layers are as designated in the z-axis: Layer 1: atmosphere (1e-30 m thick) Layer 2: upper caprock (10 m thick) Layers 3-13: Etchegoin (536.23 m thck) Layers 14-27: Macoma-Chanac (679.04 m thick) Layers 28-32: Santa Margarita-McLure (185.94 m thick) The wellbore is placed along node i=51, j=26, and extends vertically along 5 nodes from the top to the bottom of the model. Special instructions when extracting files: Each Gzip archive (ex: prod07.sim1000-sim00099.tar.gz) contains 100 simulations. Gzip archives should be transferred into base directories (ie. In Linux: mkdir prod07; mv prod07.*.tar.gz prod07/.) before extracting, or files will be overwritten. Each sub-simulation tree should have the following file structure pattern (using the linux 'tree' command): |-- prod07 | |-- sim0001 | |-- W31-0.2 | | |-- Q_brn.W31-0.2.sim0001.layers123.tab | | |-- Q_co2.W31-0.2.sim0001.layers123.tab | | |-- sim0001.W31-0.2.no_bg.trans.TDS.red.ntab | | |-- sim0001.W31-0.2.trans.pH.red.ntab | | |-- sim0001.W31-0.2.usnt.CO2_sat.deltabg.red.ntab | | |-- sim0001.W31-0.2.usnt.P.deltabg.red.ntab | |-- W31-0.5 | | |-- Q_brn.W31-0.5.sim0001.layers123.tab | | |-- Q_co2.W31-0.5.sim0001.layers123.tab | | |-- sim0001.W31-0.5.no_bg.trans.TDS.red.ntab | | |-- sim0001.W31-0.5.trans.pH.red.ntab | | |-- sim0001.W31-0.5.usnt.CO2_sat.deltabg.red.ntab | | |-- sim0001.W31-0.5.usnt.P.deltabg.red.ntab | |-- W31-1.0 | | |-- Q_brn.W31-1.0.sim0001.layers123.tab | | |-- Q_co2.W31-1.0.sim0001.layers123.tab | | |-- sim0001.W31-1.0.no_bg.trans.TDS.red.ntab | | |-- sim0001.W31-1.0.trans.pH.red.ntab | | |-- sim0001.W31-1.0.usnt.CO2_sat.deltabg.red.ntab | | |-- sim0001.W31-1.0.usnt.P.deltabg.red.ntab | |-- sim0001.mesh_k16.prod07.trans.genmsh Disclaimer This document was prepared as an account of work sponsored by an agency of the United States government. Neither the United States government nor Lawrence Livermore National Security, LLC, nor any of their employees makes any warranty, expressed or implied, or assumes any legal liability or responsibility for the accuracy, completeness, or usefulness of any information, apparatus, product, or process disclosed, or represents that its use would not infringe privately owned rights. Reference herein to any specific commercial product, process, or service by trade name, trademark, manufacturer, or otherwise does not necessarily constitute or imply its endorsement, recommendation, or favoring by the United States government or Lawrence Livermore National Security, LLC. The views and opinions of authors expressed herein do not necessarily state or reflect those of the United States government or Lawrence Livermore National Security, LLC, and shall not be used for advertising or product endorsement purposes. Lawrence Livermore National Laboratory is operated by Lawrence Livermore National Security, LLC, for the U.S. Department of Energy, National Nuclear Security Administration under Contract DE-AC52-07NA27344. This report was reviewed and released as LLNL-MI-753464.

aquifer↗

DOE OSTI2020

Measurements of inclusive and differential cross-sections of combined $ t\overline{t}\gamma$ and tWγ production in the eμ channel at 13 TeV with the ATLAS detector

Inclusive and differential cross-sections for the production of top quarks in association with a photon are measured with proton-proton collision data corresponding to an integrated luminosity of 139 fb - 1 . The data were collected by the ATLAS detector at the LHC during Run 2 between 2015 and 2018 at a centre-of-mass energy of 13 TeV. The measurements are performed in a fiducial volume defined at parton level. Events with exactly one photon, one electron and one muon of opposite sign, and at least two jets, of which at least one is b -tagged, are selected. The fiducial cross-section is measured to be $ {39.6}_{-2.3}^{+2.7} $ fb. Differential cross-sections as functions of several observables are compared with state-of-the-art Monte Carlo simulations and next-to-leading-order theoretical calculations. These include cross-sections as functions of photon kinematic variables, angular variables related to the photon and the leptons, and angular separations between the two leptons in the event. All measurements are in agreement with the predictions from the Standard Model.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2020

Land Model Testbed: Accelerating Development, Benchmarking and Analysis of Land Surface Models

A Land Model Testbed (LMT), designed to provide a computational framework for systematically assessing model fidelity and supporting rapid development of complex multiscale models, offers a general-purpose workflow for conducting large ensemble simulations of multiple land surface models, post-processing large volumes of model output, and evaluating model results. It leverages existing tools for launching model simulations and the International Land Model Benchmarking (ILAMB) package for assessing model fidelity through comparison with best-available observational datasets. Increased complexity and proliferation of uncertain parameters in process representations in land surface models has driven the need for frequent and intensive testing and evaluating of models to quantify uncertainties and optimize parameters such that results are consistent with observations. The LMT described here meets these needs by providing tools to run thousands of ensemble simulations simultaneously and post-process their output files, by automating execution of an enhanced version of ILAMB with site-specific benchmarks and multivariate functional relationships, and by offering ensemble diagnostics and a customizable dashboard for displaying model performance metrics and associated graphics. We envision the LMT capabilities will serve as a foundational computational resource for a proposed user facility focused on terrestrial multiscale model--data integration.

Sreepathi, Sarat↗

DOE OSTI2026

An Integrated Modeling Framework for Sediment Dynamics During Urban Flooding: Application to Hurricane Harvey in Houston

Floodwater can mobilize and redistribute large volumes of sediment from upland to downstream urban areas, threatening infrastructure, water quality, and ecosystem health. However, existing modeling approaches often fail to capture sediment dynamics in urban floodplains due to the lack of integration between upland hydrological processes and riverine sediment transport. This study presents the first integrated modeling framework that couples the Energy Exascale Earth System Model (E3SM) land component, which simulates runoff and hillslope erosion, with TELEMAC-GAIA, a two-dimensional hydrodynamic and sediment transport model. This framework enables the fully distributed, process-based simulation of high-resolution (as fine as 30 m) sediment dynamics from hillslopes to floodplains. Applied to a highly urbanized watershed in Houston during Hurricane Harvey, this framework reproduced observed water levels at 16 USGS gauges (median R 2 = 0.83 and KGE = 0.78), key sediment dynamics such as sediment transport and deposition processes, and reproduced spatial deposition patterns consistent with LiDAR-derived data. Based on the simulation, we estimate 8.0 million m 3 of event-scale sediment deposition, including 5.7 million m 3 trapped in the flood-control reservoirs and 2.3 million m 3 deposited along major channels and floodplains. Using a representative unit removal cost, this corresponds to an estimated dredging cost of $581 million for total deposition. These results provide a first-order, physically based quantification of Harvey-scale sediment impacts. This study provides a valuable tool for the holistic analysis of sediment dynamics triggered by extreme urban flooding, supporting flood-resilience planning. More broadly, it highlights the importance of integrating physically based hydrological processes for urban flooding and sediment research.

Hurricane Harvey↗

DOE OSTI2025

Estimating Soil Thermal Inertia Profiles From the Passive Equilibration of a Temperature Probe

Knowledge of the distribution of soil thermal properties is important for understanding subsurface hydrological and biogeochemical processes. This study describes and evaluates quick thermal profiling (QTP), a new measurement technique aimed at providing rapid, depth-resolved measurements of soil thermal inertia at numerous locations across the landscape. A cylindrical probe with temperature sensors at multiple depths is quickly inserted into the ground, and soil thermal inertia is estimated from how quickly the probe temperature equilibrates with the soil. To this end, a finite volume heat transfer model is used to generate temperature equilibration time series across combinations of controlling factors, and a gridded search inversion approach is applied to infer soil thermal inertia. Field tests in the Arctic indicate that QTP measurements have a minimum uncertainty of 0.14 J m −2 K −1 s −1/2 and covary with dual-probe heat pulse thermal analyzer measurements (concordance correlation coefficient = 0.56) with a root-mean-square error of 0.40 J m −2 K −1 s −1/2 . Besides demonstrating the value of QTP for estimating thermal inertia, this study identifies various sources of measurement uncertainty, particularly probe-soil contact resistance and frictional heating. Further, analysis of soil samples indicates that thermal inertia can be used to estimate thermal conductivity and dry bulk density in the studied area, although such inferences are highly site-specific. Overall, the QTP method holds promise to generate thermal inertia data products and to complement other characterization approaches for advancing understanding of soil properties across far more locations than is currently possible.

Lamb, J. R. [Lawrence Berkeley National Laboratory↗

DOE OSTI2022

Adaptive Sampling for In Situ Cloud Probe (Final Report)

Clouds play a leading role in the Earth's global energy and solar radiation balance and hydrological cycle. Improving cloud models requires detailed information on the cloud microphysical properties, such as droplet size distribution and number density, liquid water content and cloud composition (droplets, ice particles), which can only be provided by aerial in situ measurements. However, for many atmospheric measurement instruments, the lack of flexibility in selecting the operational mode during operation can lead to uncertainties in sampling and measurement characteristics under continuously varying atmospheric conditions. This SBIR project is developing an advanced, compact optical imaging technology for in situ characterization of cloud hydrometeors. The development involves a deep modification of the existing Mesa Photonics’ Cloud Droplet Measurement System (CDMS) in order to implement real-time automatic adaptive sampling based on the acquired in situ data and environmental parameters. The new system, CDMS-2, implements two measurement modes: side-scatter imaging for smaller hydrometeors and direct bright-field-illumination imaging for larger hydrometeors in a significantly larger sample volume. The system measures the droplet size distribution (DSD) and number density with an added capability of discriminating between liquid water and ice hydrometeors (based on polarization-resolved side-scatter imaging). The instrument will implement automatic switching or alternating between the regular side-scatter imaging mode and sparse/large hydrometeor mode (based on the acquired data). Other adaptive sampling capabilities include variable sample volume and dynamic range (based on the measured DSD). The preferred deployment platforms are uncrewed aircraft systems (UAS) and tethered balloon/kite systems (TBS). The Phase I project achieved (or exceeded) the goals listed in the Work Plan. A CDMS-2 laboratory prototype implementing the polarization-resolved side-scatter imaging mode and direct bright-field-illumination imaging mode was designed and built. Additional capabilities included the variable illumination pulse energy and sample volume. The smallest detectable droplet diameter was improved to 3–4 μm (from the nominal 10 μm value specified for the original CDMS). Discrimination between water droplets and ice particles was experimentally demonstrated. The Phase I prototype was extensively tested and calibrated in the laboratory and also tested in the Pi Cloud Chamber at Michigan Technological University (MTU). The two intensive experimental campaigns at MTU provided unique opportunities of testing the CDMS-2 laboratory prototype under realistic warm and mixed-phase cloud conditions (stable for long periods of time), testing different sampling modes and intercomparing the CDMS-2 prototype to other co-located cloud characterization instruments. The Phase I project successfully demonstrated the feasibility of the proposed technology and identified the engineering challenges of designing a field deployable prototype instrument in Phase II. The Phase I study provides a solid basis for development, characterization and field-testing of the proposed advanced cloud probe with adaptive sampling in Phase II followed by commercialization of the technology in Phase III.

47 OTHER INSTRUMENTATION↗

DOE OSTI2022

Volumetric Rendering on Wavelet-Based Adaptive Grid

Numerical modeling of physical phenomena frequently involves processes across a wide range of spatial and temporal scales. In the last two decades, the advancements in wavelet-based numerical methodologies to solve partial differential equations, combined with the unique properties of wavelet analysis to resolve localized structures of the solution on dynamically adaptive computational meshes, make it feasible to perform large-scale numerical simulations of a variety of physical systems on a dynamically adaptive computational mesh that changes both in space and time. Volumetric visualization of the solution is an essential part of scientific computing, yet the existing volumetric visualization techniques do not take full advantage of multi-resolution wavelet analysis and are not fully tailored for visualization of a compressed solution on the wavelet-based adaptive computational mesh. Our objective is to explore the alternatives for the visualization of time-dependent data on space-time varying adaptive mesh using volume rendering while capitalizing on the available sparse data representation. Two alternative formulations are explored. The first one is based on volumetric ray casting of multi-scale datasets in wavelet space. Rather than working with the wavelets at the finest possible resolution, a partial inverse wavelet transform is performed as a preprocessing step to obtain scaling functions on a uniform grid at a user-prescribed resolution. As a result, a solution in physical space is represented by a superposition of scaling functions on a coarse regular grid and wavelets on an adaptive mesh. An efficient and accurate ray casting algorithm is based just on these coarse scaling functions. Additional details are added during the ray tracing by taking an appropriate number of wavelets into account based on support overlap with the interpolation point, wavelet coefficient magnitude, and other characteristics, such as opacity accumulation (front to back ordering) and deviation from frontal viewing direction. The second approach is based on complementing of wavelet-based adaptive mesh to the traditional Adaptive Mesh Refinement (AMR) mesh. Both algorithms are illustrated and compared to the existing volume visualization software for Rayleigh-Benard thermal convection and electron density data sets in terms of rendering time and visual quality for different data compression of both wavelet-based and AMR adaptive meshes.

Vezolainen, Alexei V.↗

DOE OSTI2021

DELVE-ing into the Jet: A Thin Stellar Stream on a Retrograde Orbit at 30 kpc

We perform a detailed photometric and astrometric analysis of stars in the Jet stream using data from the first data release of the DECam Local Volume Exploration Survey DR1 and Gaia EDR3. We discover that the stream extends over ∼ 29° on the sky (increasing the known length by 18°), which is comparable to the kinematically cold Phoenix, ATLAS, and GD-1 streams. Using blue horizontal branch stars, we resolve a distance gradient along the Jet stream of 0.2 kpc deg$^{−1}$, with distances ranging from D $_{⊙}$ ∼ 27–34 kpc. We use natural splines to simultaneously fit the stream track, width, and intensity to quantitatively characterize density variations in the Jet stream, including a large gap, and identify substructure off the main track of the stream. Furthermore, we report the first measurement of the proper motion of the Jet stream and find that it is well aligned with the stream track, suggesting the stream has likely not been significantly perturbed perpendicular to the line of sight. Finally, we fit the stream with a dynamical model and find that it is on a retrograde orbit, and is well fit by a gravitational potential including the Milky Way and Large Magellanic Cloud. These results indicate the Jet stream is an excellent candidate for future studies with deeper photometry, astrometry, and spectroscopy to study the potential of the Milky Way and probe perturbations from baryonic and dark matter substructure.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

Neural network for 3D inertial confinement fusion shell reconstruction from single radiographs

In inertial confinement fusion (ICF), x-ray radiography is a critical diagnostic for measuring implosion dynamics, which contain rich three-dimensional (3D) information. Traditional methods for reconstructing 3D volumes from 2D radiographs, such as filtered backprojection, require radiographs from at least two different angles or lines of sight (LOS). In ICF experiments, the space for diagnostics is limited, and cameras that can operate on fast timescales are expensive to implement, limiting the number of projections that can be acquired. To improve the imaging quality as a result of this limitation, convolutional neural networks (CNNs) have recently been shown to be capable of producing 3D models from visible light images or medical x-ray images rendered by volumetric computed tomography. Here, we propose a CNN to reconstruct 3D ICF spherical shells from single radiographs. We also examine the sensitivity of the 3D reconstruction to different illumination models using preprocessing techniques such as pseudo-flatfielding. To resolve the issue of the lack of 3D supervision, we show that training the CNN utilizing synthetic radiographs produced by known simulation methods allows for reconstruction of experimental data as long as the experimental data are similar to the synthetic data. We also show that the CNN allows for 3D reconstruction of shells that possess low mode asymmetries. Further comparisons of the 3D reconstructions with direct multiple LOS measurements are justified.

3D Reconstruction↗

DOE OSTI2023

Smart Preprocessing & Robust Integration Emulator

To achieve the desired particle size of biomass feedstocks during preprocessing for trouble-free handling and conversion to produce biofuels and bioproducts, the raw materials must undergo a crucial milling process. The particle size of biomass plays a critical role in subsequent biofuel manufacturing, where a larger area-to-volume ratio facilitates efficient synthesis while balancing the impact of moisture on biomass storage. To optimize biofuel production efficiency and overcome these challenges, it is imperative to accurately predict the particle size distribution (PSD) of the biomass in the design of efficient preprocessing systems. The population balance model (PBM), upon empirical calibration and validation, can provide rapid prediction of post-milling PSD of granular biomass. However, PSD has limitations related to mass conservation and the absence of moisture considerations. To overcome these drawbacks, a deep learning model called the enhanced deep neural operator (DNO+) is implemented in the code. This model not only retains the capabilities of the PBM in handling complex mapping functions but also incorporates additional factors influencing the system. By considering various experimental conditions such as sieve size and moisture content, the trained DNO+ model can effectively predict the PSD after milling for any given feed PSD. To further reduce the reliance on experimental data, the PBM is integrated into the DNO+ model, resulting in a physics-informed DNO+ (PIDNO+). The PIDNO+ model addresses the non-conservation of quality exhibited by the PBM while inheriting the advantages of the DNO+ model in considering multiple influencing factors. Moreover, the PIDNO+ model significantly reduces the amount of data required for model training. Both deep learning models, i.e., DNO+ and PIDNO+, are excellent in predictive performance, offering swift and accurate machine learning-based predictions. The use of this code that contains these models will assist in guiding the proper milling equipment selection and operational conditions to achieve the desired biomass particle sizes, ensuring the efficiency of subsequent biofuel and bioproduct production processes.

Xia, Yidong [Idaho National Laboratory (INL), Idah↗

DOE OSTI2021

Total X-Ray Emission from the LMC Observed with HaloSat

We constrain the global parameters of the Large Magellanic Cloud (LMC) by modeling the first soft X-ray (0.4–6.8 keV) observations of the entire LMC produced in a single pointing with moderate spectral resolution by the HaloSat CubeSat. These data are best fit with the sum of two thermal plasma components in collisional ionization equilibrium and a power law. We find cool (0.210 ± 0.014 keV) and hot (0.89 ± 0.14 keV) components. The total X-ray luminosity of the LMC is (1.08 ± 0.14) × 10{sup 39} erg s{sup −1}. X-ray binaries make up a large fraction of the emission with a luminosity of (6.0 ± 0.8) × 10{sup 38} erg s{sup −1}, followed by cool gas from superbubbles, supernovae, and diffuse emission with a luminosity of (3.0 ± 0.3) × 10{sup 38} erg s{sup −1}. The hot gas from star formation contributes the smallest fraction, with a luminosity of (1.9 ± 0.5) × 10{sup 38} erg s{sup −1}. We estimate the total volume of the cool gas to be between (0.2–1.2) × 10{sup 10} pc{sup 3} and the hot gas to be between (1.0–5) × 10{sup 7} pc{sup 3} for filling factors of f = 1 and 0.2. These volumes result in a total thermal energy for the cool gas between (1.4–3) × 10{sup 54} erg for electron densities of 0.017–0.04 cm{sup −3}, and a thermal energy for the hot gas between (1.7–4) × 10{sup 53} erg for electron densities of 0.12–0.3 cm{sup −3}. This yields cooling timescales for the cool and hot gas of (1.5–3) × 10{sup 8} yr and (1.8–6) × 10{sup 7} yr, respectively.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2021

Offline Arterial Signal Timing Optimization for Closely Spaced Intersections

The stop-and-go driving behavior at a busy arterial causes a significant amount of fuel waste and emissions that could be reduced. For closely spaced intersections, the traditional bottom-up signal timing approach could fail under high volumes when the queues reach the upstream intersections. A well-coordinated arterial should manage the queue length to prevent the control from failing. Traditionally, coordination among signalized intersections for an arterial means finding the good offsets for each intersection. The control of the intersections along an arterial can also be coordinated by optimizing splits. Splits impact the capacity of an intersection. Similar to the highway, the capacity drop could cause growing queues. This work is based on the Virtual Phase-Link (VPL) model, a street traffic model designed for online traffic model predictive control, to obtain a top-down offline arterial signal timing. The VPL-based model can guarantee capacity consistency in intersections along an arterial. We, therefore, adopted the VPL-based model and developed an offline signal timing optimization approach. The proposed timing derived from the VPL-based offline signal timing optimization showed very good results in simulation. We also collected field experiment data, which demonstrated overall energy reductions and speed improvements on some sections of the study arterial.

47 OTHER INSTRUMENTATION↗