Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “sparse data”

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 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

On efficient Hessian computation using the edge pushing algorithm in Julia

Evaluating the Hessian matrix of second-order derivatives at a sequence of points can be costly when applying second-order methods for nonlinear optimization. In this work, we discuss our experiences implementing the recently proposed Edge Pushing (EP) method in Julia as an experimental replacement for the current colouring-based methods used by JuMP, an open-source algebraic modelling language. We propose an alternative data structure for sparse Hessians to avoid the use of hash tables and analyse the space and time complexity of EP method. In our benchmarks, we find that EP is very competitive in terms of both preprocessing time and Hessian evaluation time. As a result, we identify cases where EP closes the performance gap between JuMP's previous implementation and the implementation in AMPL, a commercial software package with similar functionality.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Measurement of the mean excitation energy of liquid argon

The mean excitation energy (I-value) of liquid argon is a critical input for energy estimation in neutrino oscillation experiments. It is measured to be 199 +4 -5 eV using the range of 402.2 MeV protons from the Fermilab Linac. This compares to the author's recent evaluation of (197 ± 7) eV based on a combination of an oscillator strength distribution analysis, gaseous argon range measurements, sparse stopping power data on solid argon, and an extrapolation of data on the effect of phase from other substances. Using all sources of information, we recommend a value of (198 ± 4) eV for liquid argon, which is significantly higher than 188 eV, from ICRU-37's gaseous argon evaluation, commonly used in Monte Carlo codes such as Geant4.

47 OTHER INSTRUMENTATION↗
DOE OSTI2023

The evolution of HCO + in molecular clouds using a novel chemical post-processing algorithm

Modelling the chemistry of molecular clouds is critical to accurately simulating their evolution. To reduce computational cost, 3D simulations generally restrict their chemistry to species with strong heating and cooling effects. Time-dependent information about the evolution of other species is therefore often neglected. We address this gap by post-processing tracer particles in the SILCC-Zoom molecular cloud simulations. Using a chemical network of 39 species and 301 reactions (including freeze-out of CO and H 2 O) and a novel algorithm to reconstruct a density grid from sparse tracer particle data, we produce time-dependent density distributions for various species. We focus upon the evolution of HCO + , which is a critical formation reactant of CO but is not typically modelled on the fly. We find that ∼ 90 per cent of the HCO + content of the cold molecular gas forms in situ around n HCO + ∼ 10 3 –10 4 cm −3 , over a time-scale of approximately 1 Myr. The remaining ∼ 10 per cent forms at high extinction sites, with minimal turbulent mixing out into the less dense gas. We further show that the dominant HCO + formation pathway is dependent on the visual extinction, with the reaction H 3 + + CO contributing 90 per cent of the total HCO + production above A V, 3D = 3. We produce the very first maps of the HCO + column density, N(HCO + ), and show that it reaches values as high as 10 15 cm −2 . We find that 50 per cent of the HCO + mass is located within AV ∼ 10–30 in a density range of 10 3.5 –10 4.5 cm −3 . Our maps of N(HCO + ) are shown to be in good agreement with recent observations of the W49A star-forming region.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Connect the Dots: In Situ 4-D Seismic Monitoring of CO 2 Storage With Spatio-Temporal CNNs

4-D seismic imaging has been widely used in CO 2 sequestration projects to monitor the fluid flow in the volumetric subsurface region that is not sampled by wells. Ideally, real-time monitoring and near-future forecasting would provide site operators with great insights to understand the dynamics of the subsurface reservoir and assess any potential risks. However, due to obstacles such as high deployment cost, availability of acquisition equipment, exclusion zones around surface structures, only very sparse seismic imaging data can be obtained during monitoring. That leads to an unavoidable and growing knowledge gap over time. The operator needs to understand the fluid flow throughout the project lifetime and the seismic data are only available at a limited number of times. This is insufficient for understanding reservoir behavior. To overcome those challenges, we have developed spatio-temporal neural-network-based models that can produce high-fidelity interpolated or extrapolated images effectively and efficiently. Specifically, our models are built on an autoencoder, and incorporate the long short-term memory (LSTM) structure with a new loss function regularized by optical flow. We validate the performance of our models using real 4-D post-stack seismic imaging data acquired at the Sleipner CO 2 sequestration field. We employ two different strategies in evaluating our models. Numerically, we compare our models with different baseline approaches using classic pixel-based metrics. We also conduct a blind survey and collect a total of 20 responses from domain experts to evaluate the quality of data generated by our models. Finally, via both numerical and expert evaluation, we conclude that our models can produce high-quality 2-D/3-D seismic imaging data at a reasonable cost, offering the possibility of real-time monitoring or even near-future forecasting of the CO 2 storage reservoir.

4-D seismic imaging↗
DOE OSTI2017

Direct coupling of tomography and ptychography

We present a generalization of the ptychographic phase problem for recovering refractive properties of a three-dimensional object in a tomography setting. Our approach, which ignores the lateral overlapping probe requirements in existing ptychography algorithms, can enable the reconstruction of objects using highly flexible acquisition patterns and pave the way for sparse and rapid data collection with lower radiation exposure.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

High-Resolution 2D Lg and Pg Attenuation Models in the Basin and Range Region with Implications for Frequency-Dependent Q

Here, we develop high–resolution, laterally varying attenuation models for the regional crustal phases of Pg and Lg in the area surrounding the Basin and Range Province in the western United States. The models are part of the characterization effort for the Source Physics Experiment (SPE), a series of chemical explosions at the Nevada National Security Site designed to improve our understanding of explosion source phenomenology. To aid in SPE modeling efforts, we focus on improving our ability to accurately predict amplitudes in a set of narrow frequency bands ranging from 0.5 to 16.0 Hz. To explore constraints at higher frequencies where data become more sparse, we test the robustness of the empirically observed power–law relationship between quality factor Q and frequency (Q=Q 0 f γ ). Our methodology uses a staged approach to consider attenuation, physics–based source terms, site terms, and geometrical spreading contributions to amplitude measurements. Tomographic inversion results indicate that the frequency dependence is a reasonable assumption as attenuation varies laterally for this region through all frequency bands considered. Our 2D Pg and Lg attenuation models correlate with underlying physiographic provinces, with the highest Q located in the Sierra Nevada Mountains and the Colorado plateau. Compared to a best–fitting 1D model for the region, the 2D model provides an 81% variance reduction overall for Lg residuals and a 75% reduction for Pg . These detailed attenuation maps at high frequencies will facilitate further study of local and regional distance P/S amplitude discriminants that are typically used to distinguish between earthquakes and underground explosions.

58 GEOSCIENCES↗
DOE OSTI2017

Streamline-Based Time-Lapse-Seismic-Data Integration Incorporating Pressure and Saturation Effects

In this paper, we present an efficient history-matching technique that simultaneously integrates 4D repeat seismic surveys with well-production data. This approach is particularly well-suited for the calibration of the reservoir properties of high-resolution geologic models because the seismic data are a really dense but sparse in time, whereas the production data are finely sampled in time but spatially averaged. The joint history matching is performed by use of streamline-based sensitivities derived from either finite-difference or streamline-based flow simulation. For the most part, earlier approaches have focused on the role of saturation changes, but the effects of pressure have largely been ignored. Here, we present a streamline-based semianalytic approach for computing model-parameter sensitivities, accounting for both pressure and saturation effects. The novelty of the method lies in the semianalytic sensitivity computations, making it computationally efficient for high-resolution geologic models. The approach is implemented by use of a finite-difference simulator incorporating the detailed physics. Its efficacy is demonstrated by use of both synthetic and field applications. Lastly, for both the synthetic and the field cases, the advantages of incorporating the time-lapse variations are clear, seen through the improved estimation of the permeability distribution, the pressure profile, the evolution of the fluid saturation, and the swept volumes.

58 GEOSCIENCES↗
DOE OSTI2016

Comparison of continuous and discrete-time data-based modeling for hypoelliptic systems

In this work, we compare two approaches to the predictive modeling of dynamical systems from partial observations at discrete times. The first is continuous in time, where one uses data to infer a model in the form of stochastic differential equations, which are then discretized for numerical solution. The second is discrete in time, where one directly infers a discrete-time model in the form of a nonlinear autoregression moving average model. The comparison is performed in a special case where the observations are known to have been obtained from a hypoelliptic stochastic differential equation. We show that the discrete-time approach has better predictive skills, especially when the data are relatively sparse in time. Finally , we discuss open questions as well as the broader significance of the results.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2018

Final Technical Report for Collaborative Research: Developing and Implementing Ocean-Atmosphere Reanalyses for Climate Applications (OARCA)

As an important step toward a coupled data assimilation system for generating reanalysis fields needed to assess climate model projections, the Ocean Atmosphere Coupled Reanalysis for Climate Applications (OARCA) project assesses and improves the longest reanalyses currently available of the atmosphere and ocean: the 20th Century Reanalysis Project (20CR) and the Simple Ocean Data Assimilation with sparse observational input (SODAsi) system, respectively. In this project, we make off-line but coordinated improvements in the 20CR and SODAsi datasets, with improvements in one feeding into improvements of the other through an iterative generation of new versions. These datasets now span from the 19th to 21st centuries. We then study the extreme weather and variability from days to decades of the resulting datasets. A total of 24 publications have been produced in this project.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Efficient Learning of Accurate Surrogates for Simulations of Complex Systems

Machine learning methods are increasingly deployed to construct surrogate models for complex physical systems at a reduced computational cost. However, the predictive capability of these surrogates degrades in the presence of noisy, sparse, or dynamic data. We introduce an online learning method empowered by optimizer-driven sampling that has two advantages over current approaches: it ensures that all local extrema (including endpoints) of the model response surface are included in the training data and it employs a continuous validation and update process in which surrogates undergo retraining when their performance falls below a validity threshold. We find, using benchmark functions, that optimizer-directed sampling generally outperforms traditional sampling methods in terms of accuracy around local extrema even when the scoring metric is biased towards assessing overall accuracy. Finally, the application to dense nuclear matter demonstrates that highly accurate surrogates for a nuclear equation of state model can be reliably auto-generated from expensive calculations using few model evaluations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Evaluation of Hypotheses for the Fate of Historical Discharges from Hanford Site BC Cribs and Trenches

A preliminary review of declassified reports and historical documents from the 1950s and 1960s indicated they may contain valuable information regarding potential impacts to groundwater from liquid waste discharged to the 216-BC Cribs and Trenches during the late 1950s. Transport of contaminants in the vadose zone beneath the BC Cribs area is simulated by a vadose zone model developed for the Composite Analysis (DOE/RL-2019-52, Composite Analysis for Low-Level Waste Disposal in the Hanford Site Central Plateau ). Results from this model indicate there has been no substantial breakthrough of contaminants (particularly technetium-99) to groundwater from the vadose zone up to the present day. However, there is uncertainty in the anisotropy and other model input parameters as well as acknowledged sparseness of sampling data for this area. In this study, historical documents and other information were reviewed to evaluate the possibility that contaminant breakthrough to groundwater occurred during or shortly after operation of the BC Cribs sites, and whether refinements of the vadose zone model parametrization should be considered.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2020

Reduced Models for Chemical Kinetics derived from Parallel Ensemble Simulations of Stirred Reactors

Detailed chemical kinetic models are essential to improve the predictive capability of computational fluid dynamics (CFD) simulations and to impact the development of next generation combustion technologies. However, such detailed kinetic models are beyond the reach of CFD codes due to the computational cost of transporting additional chemical degrees of freedom. Reduced order models for combustion kinetics exploit the presence of lower dimensional manifolds in the state space and require far fewer number of scalars to be transported as part of CFD simulations. In this paper, we describe the software infrastructure developed to allow large ensembles of stirred reactor configurations with detailed chemical kinetics. The canonical reactor simulations are incorporated in a workflow with sparse grid and data reduction methodologies aimed at development of reduced order models. We present results from analyzing the results of the stirred reactor calculations for syngas mixture using principal component analysis (PCA) and neural network methodologies. The results show that nonlinear autoencoder model can surpass PCA in representing the original dataset with far fewer degrees of freedom.

Zhang, Pei↗
DOE OSTI2016

SPECTRUM AND ANISOTROPY OF TURBULENCE FROM MULTI-FREQUENCY MEASUREMENT OF SYNCHROTRON POLARIZATION

We consider turbulent synchrotron-emitting media that also exhibit Faraday rotation and provide a statistical description of synchrotron polarization fluctuations. In particular, we consider these fluctuations as a function of the spatial separation of the direction of the measurements and as a function of wavelength for the same line of sight. On the basis of our general analytical approach, we introduce several measures that can be used to obtain the spectral slopes and correlation scales of both the underlying magnetic turbulence responsible for emission and the spectrum of the Faraday rotation fluctuations. We show the synergetic nature of these measures and discuss how the study can be performed using sparsely sampled interferometric data. We also discuss how additional characteristics of turbulence can be obtained, including the turbulence anisotropy and the three-dimensional direction of the mean magnetic field. In addition, we consider the cases when the synchrotron emission and Faraday rotation regions are spatially separated. Appealing to our earlier study, we explain that our new results are applicable to a wide range of spectral indexes of relativistic electrons responsible for synchrotron emission. We expect wide application of our techniques, both with existing synchrotron data sets and with big forthcoming data sets from LOFAR and SKA.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

SUBSTRUCTURE WITHIN THE SSA22 PROTOCLUSTER AT z ≈ 3.09

We present the results of a densely sampled spectroscopic survey of the SSA22 protocluster at z ≈ 3.09. Our sample with Keck/LRIS spectroscopy includes 106 Ly α emitters (LAEs) and 40 Lyman break galaxies (LBGs) at z = 3.05–3.12. These galaxies are contained within the 9′ × 9′ region in which the protocluster was discovered, which also hosts the maximum galaxy overdensity in the SSA22 region. The redshift histogram of our spectroscopic sample reveals two distinct peaks, at z = 3.069 (blue; 43 galaxies) and z = 3.095 (red; 103 galaxies). Furthermore, objects in the blue and red peaks are segregated on the sky, with galaxies in the blue peak concentrating toward the western half of the field. These results suggest that the blue and red redshift peaks represent two distinct structures in physical space. Although the double-peaked redshift histogram is traced in the same manner by LBGs and LAEs, and brighter and fainter galaxies, we find that 9 out of 10 X-ray AGNs in SSA22, and all 7 spectroscopically confirmed giant Ly α “blobs,” reside in the red peak. We combine our data set with sparsely sampled spectroscopy from the literature over a significantly wider area, finding preliminary evidence that the double-peaked structure in redshift space extends beyond the region of our dense spectroscopic sampling. In order to fully characterize the three-dimensional structure, dynamics, and evolution of large-scale structure in the SSA22 overdensity, we require the measurement of large samples of LAE and LBG redshifts over a significantly wider area, as well as detailed comparisons with cosmological simulations of massive cluster formation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

A production-tagged aerosol module for Earth system models, OsloAero5.3 – extensions and updates for CAM5.3-Oslo

Abstract. We document model updates and present and discuss modeling and validation results from a further developed production-tagged aerosol module, OsloAero5.3, for use in Earth system models. The aerosol module has in this study been implemented and applied in CAM5.3-Oslo. This model is based on CAM5.3–CESM1.2 and its own predecessor model version CAM4-Oslo. OsloAero5.3 has improved treatment of emissions, aerosol chemistry, particle life cycle, and aerosol–cloud interactions compared to its predecessor OsloAero4.0 in CAM4-Oslo. The main new features consist of improved aerosol sources; the module now explicitly accounts for aerosol particle nucleation and secondary organic aerosol production, with new emissions schemes also for sea salt, dimethyl sulfide (DMS), and marine primary organics. Mineral dust emissions are updated as well, adopting the formulation of CESM1.2. The improved model representation of aerosol–cloud interactions now resolves heterogeneous ice nucleation based on black carbon (BC) and mineral dust calculated by the model and treats the activation of cloud condensation nuclei (CCN) as in CAM5.3. Compared to OsloAero4.0 in CAM4-Oslo, the black carbon (BC) mass concentrations are less excessive aloft, with a better fit to observations. Near-surface mass concentrations of BC and sea salt aerosols are also less biased, while sulfate and mineral dust are slightly more biased. Although appearing quite similar for CAM5.3-Oslo and CAM4-Oslo, the validation results for organic matter (OM) are inconclusive, since both of the respective versions of OsloAero are equipped with a limited number of OM tracers for the sake of computational efficiency. Any information about the assumed mass ratios of OM to organic carbon (OC) for different types of OM sources is lost in the transport module. Assuming that observed OC concentrations scaled by 1.4 are representative for the modeled OM concentrations, CAM5.3-Oslo with OsloAero5.3 is slightly inferior for the very sparsely available observation data. Comparing clear-sky column-integrated optical properties with data from ground-based remote sensing, we find a negative bias in optical depth globally; however, it is not as strong as in CAM4-Oslo, but has positive biases in some areas typically dominated by mineral dust emissions. Aerosol absorption has a larger negative bias than the optical depth globally. This is reflected in a lower positive bias in areas where mineral dust is the main contributor to absorption. Globally, the low bias in absorption is smaller than in CAM4-Oslo. The Ångström parameter exhibits small biases both globally and regionally, suggesting that the aerosol particle sizes are reasonably well represented. Cloud-top droplet number concentrations over oceans are generally underestimated compared to satellite retrievals, but seem to be overestimated downwind of major emissions of dust and biomass burning sources. Finally, we find small changes in direct radiative forcing at the top of the atmosphere, while the cloud radiative forcing due to anthropogenic aerosols is now more negative than in CAM4-Oslo, being on the strong side compared to the multi-model estimate in IPCC AR5. Although not all validation results in this study show improvement for the present CAM5.3-Oslo version, the extended and updated aerosol module OsloAero5.3 is more advanced and applicable than its predecessor OsloAero4.0, as it includes new parameterizations that more readily facilitate sensitivity and process studies and use in climate and Earth system model studies in general.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Tensor Extraction of Latent Features (TELF)

Tensor ELF is a user-friendly parallel tensor decomposition Python toolbox that includes a suite of machine learning algorithms for CPU and GPU architectures for the analysis of sparse and dense data including utility tools for pre-processing and post-processing.

Eren, Maksim↗
DOE OSTI2022

Guidelines for Publicly Archiving Terrestrial Model Data to Enhance Usability, Intercomparison, and Synthesis

Scientific communities are increasingly publishing data to evaluate, accredit, and build on published research. However, guidelines for curating data for publication are sparse for model-related research, limiting the usability of archived simulation data. In particular, there are no established guidelines for archiving data related to terrestrial models that simulate land processes and their coupled interactions with climate. Terrestrial modelers have a unique set of challenges when publishing data due to the diversity of scientific domains, research questions, and the types and scales of simulations. Researchers in the U.S. Department of Energy’s (DOE) projects use a variety of multiscale models to advance robust predictions of terrestrial and subsurface ecosystem processes. Here, we synthesize archiving needs for data associated with different DOE models, and provide guidelines for publishing terrestrial model data components following FAIR (Findable, Accessible, Interoperable, Reusable) principles. The guidelines recommend archiving model inputs and testing data used in final simulation runs along with associated codes, workflow scripts, and metadata in public repositories. Researchers should consider archiving model outputs if they are within the storage limits of the repository. We also provide considerations for how to bundle files into different data publications with citable digital object identifiers. Finally, we identify repository features and tools that would enable storage and reuse of model data. Given the diversity of DOE terrestrial models, these guidelines are transferable to other model types and will enable efficient reuse of simulation data for purposes such as model intercomparisons, initialization, benchmarking, synthesis, and comparisons with field observations.

58 GEOSCIENCES↗
DOE OSTI2022

Experimental Steps toward a Density Law for Chlorine-Crediting Criticality Models of Aqueous Plutonium solutions

Criticality accident prevention is an essential safety consideration for all operations involving fissionable and fissile materials. Predictive criticality calculations using advanced neutron transport codes, such as Monte Carlo N-Particle Transport (MCNP), are invaluable tools for designing and implementing operational limits. While indispensable, these tools are limited by the quality and accuracy of the inputs that the user provides to define the modeled system. Parameters such as atomic composition, shape, and material density must be accurately defined to obtain a meaningful result. In the case of material property details for a plutonium model, accurate material characterization data is sometimes sparse. Thus limitations in calculation accuracy can be reduced by improving our knowledge of material properties in the targeted fissile systems. One of the biggest remaining challenges to accurately defining fissile systems is the description of aqueous fissile solutions. Even simple properties, such as density, are not well-known for fissile solutions relevant to nuclear energy and security. Described here are initial efforts undertaken to improve criticality calculation inputs for fissile plutonium chloride solutions in water. This effort is focused on experimentally determining accurate solution characteristics for ternary plutonium chloride/hydrochloric acid/water systems, by measurement of water activity and solution density. The effect of inputting experimental densities for these solutions into MCNP criticality calculations is compared to the traditional approach of modeling an idealized (and fictitious) plutonium-water mixture. Expansion of these efforts to a working density law for aqueous plutonium chloride solutions is also discussed.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗