Search NASA⌕ Search

SEARCH · Search NASA

Results for “box modeling”

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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Using Knowledge-Guided Machine Learning To Assess Patterns of Areal Change in Waterbodies across the Contiguous United States

Lake and reservoir surface areas are an important proxy for freshwater availability. Advancements in machine learning (ML) techniques and increased accessibility of remote sensing data products have enabled the analysis of waterbody surface area dynamics on broad spatial scales. However, interpreting the ML results remains a challenge. While ML provides important tools for identifying patterns, the resultant models do not include mechanisms. Thus, the “black-box” nature of ML techniques often lacks ecological meaning. Using ML, we characterized temporal patterns in lake and reservoir surface area change from 1984 to 2016 for 103,930 waterbodies in the contiguous United States. We then employed knowledge-guided machine learning (KGML) to classify all waterbodies into seven ecologically interpretable groups representing distinct patterns of surface area change over time. Many waterbodies were classified as having “no change” (43%), whereas the remaining 57% of waterbodies fell into other groups representing both linear and nonlinear patterns. This analysis demonstrates the potential of KGML not only for identifying ecologically relevant patterns of change across time but also for unraveling complex processes that underpin those changes.

58 GEOSCIENCES↗
DOE OSTI2020

A new soil mixing layer model for simulating conservative solute loss from initially saturated soil to surface runoff

Surface water pollution due to solute loss from soil to surface runoff is a serious threat to the environment, and it is of great significance to predict the solute loss. While soil mixing layer theory has been widely used to predict the solute loss, it has been observed that the conventional soil mixing layer models with constant mixing layer depth underestimate solute concentration in surface runoff, and thus underestimate the solute loss after surface runoff starts. This study introduced a new soil mixing layer model with a new concept that the mixing layer depth increases downward with time after surface runoff occurs. This new concept was implemented by adding a mixing front between the soil mixing layer and its underlying soil layer, which is conceptually similar to the wetting front of soil infiltration given by the Green-Ampt equation. This is equivalent to treating the soil mixing layer as a deforming control volume. A mechanistic explanation for the new model was presented. This study also developed a closed form solution of the new model. The new model was evaluated by simulating five laboratory experiments conducted by two groups of researchers. The first two experiments have ponding water, and consider different drainage and runoff conditions at both early and late runoff times. The other three experiments do not have ponding water, and consider different infiltration conditions at the bottom of soil boxes. Simulation results of the five experiments show that the new model resolves the underestimation problem of conventional soil mixing layer models. The new model is subject to a total of six limits such as that the new model only considers conservative solute. Three limits are related to the assumptions used to derive the closed form solution, and they may be alleviated in a future research as briefly discussed at the end of the paper.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

hIPPYlib-MUQ: A Bayesian Inference Software Framework for Integration of Data with Complex Predictive Models under Uncertainty

Bayesian inference provides a systematic framework for integration of data with mathematical models to quantify the uncertainty in the solution of the inverse problem. However, the solution of Bayesian inverse problems governed by complex forward models described by partial differential equations (PDEs) remains prohibitive with black-box Markov chain Monte Carlo (MCMC) methods. We present hIPPYlib-MUQ, an extensible and scalable software framework that contains implementations of state-of-the art algorithms aimed to overcome the challenges of high-dimensional, PDE-constrained Bayesian inverse problems. These algorithms accelerate MCMC sampling by exploiting the geometry and intrinsic low-dimensionality of parameter space via derivative information and low rank approximation. The software integrates two complementary open-source software packages, hIPPYlib and MUQ. hIPPYlib solves PDE-constrained inverse problems using automatically-generated adjoint-based derivatives, but it lacks full Bayesian capabilities. MUQ provides a spectrum of powerful Bayesian inversion models and algorithms, but expects forward models to come equipped with gradients and Hessians to permit large-scale solution. By combining these two complementary libraries, we created a robust, scalable, and efficient software framework that realizes the benefits of each and allows us to tackle complex large-scale Bayesian inverse problems across a broad spectrum of scientific and engineering disciplines. To illustrate the capabilities of hIPPYlib-MUQ, we present a comparison of a number of MCMC methods available in the integrated software on several high-dimensional Bayesian inverse problems. These include problems characterized by both linear and nonlinear PDEs, various noise models, and different parameter dimensions. The results demonstrate that large (~ 50×) speedups over conventional black box and gradient-based MCMC algorithms can be obtained by exploiting Hessian information (from the log-posterior), underscoring the power of the integrated hIPPYlib-MUQ framework.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Comprehensive framework for data-driven model form discovery of the closure laws in thermal-hydraulics codes

The two-phase two-fluid model is a basis of many thermal-hydraulics codes used in design, licensing, and safety considerations of nuclear power plants. Thermal-hydraulics codes rely on the closure laws to close the system of conservation equations and describe the interactions between phases. These laws, derived from years of experimental investigations, are semi-empirical correlations that lack generality and have a limited range of applicability. Increase of computational power, availability of new experiments, and development of high-fidelity simulations has increased the number of validation data. The discrepancies between the code predictions and the validation data are a great source of knowledge. Missing physics that are not included in the model but are important for the considered phenomena can be discovered by propagating the information from the experimental results through the model. Furthermore, physics-discovered data-driven model form (P3DM) methodology integrates available integral effect tests and separate effects tests to determine the necessary corrections to the model form of the closure laws. In contrast to existing calibration techniques, the methodology modifies the functional form of the closure laws. Based on the functional form of the correction, the missing physics that were not included in the original model can be discovered. The methodology provides the alternative to the machine learning approach, in which the model is discovered in the form of the intractable black-box relation. In this work, the methodology was applied to the CTF subchannel code to improve the prediction of the two-phase flow phenomena.

42 ENGINEERING↗
DOE OSTI2021

A HPC Theory-Guided Machine Learning Cyberinfrastructure for Communicating Hydrometeorological Data Across Scales

High-resolution predictions of hydrometeorological variables are critical for supporting hydropower generation decisions and flood control at hydroelectric power plants. Traditional climate and hydrologic models rely on the numerical simulation of detailed physical processes. Therefore, running these simulations is time-, labor-, and computation-intensive. Improving the spatial and temporal resolution in these modeling outputs could lead to cubic increases in both the simulation time and computational demands, rendering high-resolution hydrometeorological predictions expensive and impractical. Many past studies apply the super resolution (SR) technique to downscale climate models using deep learners. However, deep learners are deemed “black-boxes,” as their derivation processes from low-resolution outputs to high-resolution outputs are often hidden. Their results are difficult for domain scientists to interpret and validate. Thus, there is a need for an exploratory machine learning approach that can partially integrate domain-specific theory and knowledge into the data-driven mapping process between simulation outputs of different spatial scales. The domain-specific theory and knowledge can be incorporated into the data model through an inductive approach in which process-related environmental variables are used and analyzed as key drivers (i.e., environmental surrogates) to reflect the complex physical processes. Many of these variables, such as land use land cover, soil types, topography, digital elevation, air temperature, and various watershed characteristics, can be directly measured through sensors or remote sensing techniques. Additionally, SR applications that can downscale hydrological and hydrodynamics models to efficiently produce high-resolution (1 m) flood depth grids are still rare. Since the flood depth grid can be used to support critical decisions for flood control operation at hydroelectric power plants, it is crucial to enable an SR-based capability for interpolating high-resolution flood inundation maps.

13 HYDRO ENERGY↗
DOE OSTI2025

Micro2Macro: Origins of Climate Change Uncertainty

Where we are: Global Earth system models (ESMs) are essential tools for seasonal-to decadal environmental predictions, which decision-makers across government and industry require. However, uncertainties originating at the microphysical scale, i.e., in processes occurring on scales smaller than typical ESM grid boxes, remain a major challenge. While resolution continues to improve in our predictive models, we will need to parameterize microphysical processes that contribute the bulk of prediction uncertainty for the foreseeable future. Microphysical uncertainties contribute broadly to remaining limitations on Earth system predictability on seasonal and decadal time scales.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Electroweak box diagrams on the lattice for pion and neutron decay

CKM matrix is unitary by construction in the standard model(SM). The recent analyses on the first row of CKM matrix show ≈3σ tension with unitarity. Nonperturbative calculations of the radiative corrections can reduce the theory uncertainty in CKM matrix elements. Here we compute the electroweak box contribution to the pion and kaon β decays using seven Nf = 2 + 1 + 1 HISQ-Clover lattice with various pion mass and lattice spacing. The continuum and chiral limit is taken using the leading dependence on M π and a, where M π extrapolation is taken to the physical pion mass and SU(3) symmetric mass for pion and kaon box contribution, respectively. Our results are $\square^{VA}_{γW}|_{π}$ = 2.820(28) × 10 –3 and $\square^{VA}_{γW}|_{K}$ = 2.384(17) × 10 –3 .

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Speed vs Accuracy: Effect on Ligand Pose Accuracy of Varying Box Size and Exhaustiveness in AutoDock Vina

Abstract Structure‐based virtual high‐throughput screening involves docking chemical libraries to targets of interest. A parameter pertinent to the accuracy of the resulting pose is the root mean square deviation (RMSD) from a known crystallographic structure, i. e., the ‘docking power’. Here, using a popular algorithm, Autodock Vina, as a model program, we evaluate the effects of varying two common docking parameters: the box size (the size of docking search space) and the exhaustiveness of the global search (the number of independent runs starting from random ligand conformations) on the RMSD from the PDBbind v2017 refined dataset of experimental protein‐ligand complexes. Although it is clear that exhaustiveness is an important parameter, there is wide variation in the values used, with variation between 1 and >100. We, therefore, evaluated a combination of cubic boxes of different sizes and five exhaustiveness values (1, 8, 25, 50, 75, 100) within the range of those commonly adopted. The results show that the default exhaustiveness value of 8 performs well overall for most box sizes. In contrast, for all box sizes, but particularly for large boxes, an exhaustiveness value of 1 led to significantly higher median RMSD (mRMSD) values. The docking power was slightly improved with an exhaustiveness of 25, but the mRMSD changes little with values higher than 25. Therefore, although low exhaustiveness is computationally faster, the results are more likely to be far from reality, and, conversely, values >25 led to little improvement at the expense of computational resources. Overall, we recommend users to use at least the default exhaustiveness value of 8 for virtual screening calculations.

Agarwal, Rupesh↗
DOE OSTI2021

Quantum dynamics of a nucleon in the Fermi accelerator

Highlights: • In nuclear physics, the Fermi accelerator is a coherent open quantum system. • There is not a true equilibrium for a nucleon‘s mean-energy in the Fermi accelerator. • Occupation of higher energy eigen-states remains even after the mean-energy plateaus. • This work paves the way for work on quantum tunnelling in the Fermi accelerator. The quantum dynamics of a particle in a one-dimensional box with an oscillating wall (the Fermi accelerator) is investigated. The model is applied to the motion of a single nucleon in the mean-field potential of a heavy atomic nucleus whose surface vibrates. By directly solving the time-dependent Schrödinger equation, both the state of the particle and its mean-energy are studied. The effects of the frequency of the wall oscillation on the nucleon’s energy are addressed. Its energy oscillates in phase with the moving wall for all frequencies, showing no chaotic behaviour. There is a large initial peak of the nucleon’s energy as the particle adjusts to the sudden change in the size of the box and a varying relaxation time as it plateaus towards lower energy and a partial equilibrium. Small oscillations in energy continue, since there cannot be a true equilibrium while the wall is moving. The quantum coherence between the different parts of the nucleon’s wave-function in real space is very much preserved. This research lays the foundation for future investigations into quantum tunnelling in the Fermi accelerator.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

PORFLOW Implementation of Special Waste from Models for Slit and Engineered Trenches in the E-Area Low Level Waste Facility Performance Assessment

The implementation of special waste forms (SWFs) in the deterministic PORFLOW Slit Trench (ST) and Engineered Trench (ET) models to be used in the next revision of the E-Area Low Level Waste Facility’s (ELLWFs) Performance Assessment (PA) is outlined in this report. Four SWF implementation methods will be used: effective K d , delayed release, solubility-controlled/diffusion-controlled release, and complex SWF model updates. In addition, the implementation of models that address the presence of tall used equipment storage boxes in ST08-10 is considered a special waste form and will be described.

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

Advancing Enhanced Raman Spectroscopy to Detect and Measure NOW Pheromones (CPRB 2022 proposal)

Insect sex pheromones are chemical compounds that insects release to attract their partners over distances of hundreds of meters or even kilometers, in complete darkness and without any audible signals. Use of synthetic forms of key compounds have in some cases become an essential component of monitoring and/or managing key pests of agricultural crops, including navel orangeworm (Amyelois transitella) (NOW) in California tree nuts. These pheromone-based strategies can include monitoring, mating disruption, mass trapping, attract-and-kill and push-pull. There are currently multiple commercially available mating disruption products available for NOW and recent studies have demonstrated that they can be effectively used to reduce crop damage by this pest. Just how mating disruption works is not fully established and likely varies across products and target species. For instance, the extent to which synthetic pheromones compete with natural pheromone is not well understood, or in the case of monitoring, how efficaciously the insect follows the diffusing plumes, especially across large blocks and at plot borders. Furthermore, there may be specific conditions under which poor or impeded diffusion of synthetic pheromone diminishes the disruption effect, which could result in some males effectively locating females for reproduction. At the same time, many commercially available synthetic lures are also available, and growers have been effectively using them to track population development – although the attractive range of these lures remains unclear. That is, while pheromone lures can attract many moths, the relationship between trap capture and local populations, much less crop damage, remains unclear. Regardless of the emission source (aerosol puffers, meso emitters, lures etc.), we currently lack the ability to fully understand how these synthetic pheromone compounds diffuse away from their point-source of emission – and subsequently how this might affect the efficacy of mating disruption and/or the accuracy of monitoring efforts. Raman spectroscopy (RS) may offer a means to detect synthetic pheromone because of its analytical properties and recent technological advancements. In combination with nanostructured probes, RS has a demonstrated capability to detect volatiles at extremely low concentrations. While other methods for detecting pheromones have proven to be successful, such analytical methods like mass spectrometry (MS) gas chromatography (GC) require sampling in the fields and are not real-time. Other methods such as EAG/ EAD are also invasive (use antennas) but are of course valuable. Opto-electronic, micro and nano sensors just recently enticed interest. We fall in this latter category, proposing the use of optical spectroscopy (Raman) in a fiber aided by nanostructures to provide a chemical fingerprint without transductions and in real-time (i.e. no sampling). To date, we have been able to generate Raman signatures for the main chemicals in synthetic pheromones (a first for science, to the best of our knowledge), observe their trends at different concentrations, and also provide a clear distinction of such signatures in mixed samples via Principal Component Analysis (PCA) – all of which strengthens our conclusions on the ability of RS to detect emissions in more realistic conditions. We also established processes to enhance the detection sensitivity via preconcentration with Solid Phase Micro Extraction (SPME) fibers and Surface Enhanced Raman Spectroscopy (SERS) to increase the overall optical signals. Furthermore, we explored limitations and requirements for detecting from dispensers, lures, or live moths, including Dept. of Energy and CA State regulatory approvals. Finally, we setup some of our first “orchard in a box” experiments to measure pheromone diffusion along with developing appropriate models.

42 ENGINEERING↗
DOE OSTI2024

Reduced Geometry for the NB5 Velocity Selector Shielding Box

This report details the reduced NB5 velocity selector shielding box that has been isolated from the complete HFIR CGH model in order to speed up detailed calculations of the velocity selector shielding. This shielding model is intended for outsourcing shielding optimization calculations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Fast, High-fidelity Lyα Forests with Convolutional Neural Networks

Full-physics cosmological simulations are powerful tools for studying the formation and evolution of structure in the universe but require extreme computational resources. Here, we train a convolutional neural network to use a cheaper N-body-only simulation to reconstruct the baryon hydrodynamic variables (density, temperature, and velocity) on scales relevant to the Lyα forest, using data from Nyx simulations. We show that our method enables rapid estimation of these fields at a resolution of ~20 kpc, and captures the statistics of the Lyα forest with much greater accuracy than existing approximations. Because our model is fully convolutional, we can train on smaller simulation boxes and deploy on much larger ones, enabling substantial computational savings. Furthermore, as our method produces an approximation for the hydrodynamic fields instead of Lyα flux directly, it is not limited to a particular choice of ionizing background or mean transmitted flux.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Do Minihalos Affect Cosmic Reionization?

The role of unresolved structures (“minihalos”) in determining the consumption of ionizing photons during cosmic reionization remains an unsolved problem in modeling cosmic reionization, despite recent extensive studies with small-box high-resolution simulations by Park et al. and Chan et al., because the small-box studies are not able to fully sample all environments. In this paper these simulations are combined with large-box simulations from the “Cosmic Reionization On Computers” (CROC) project, allowing one to account for the full range of environments and to produce an estimate for the number of recombinations per hydrogen atom that are missed in large-scale simulations like CROC or Thesan. I find that recombinations in unresolved minihalos are unimportant compared to recombinations produced in large-scale cosmic structures and inside more massive, fully resolved halos. Since both the Park et al. and Chan et al. studies have severe limitations, the conclusions of this paper may need to be verified with more representative sets of small-box high-resolution simulations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

A field guide for aging passerine nestlings using growth data and predictive modeling

Abstract Background Accurate nestling age is valuable for studies on nesting strategies, productivity, and impacts on reproductive success. Most aging guides consist of descriptions and photographs that are time consuming to read and subjective to interpret. The Western Bluebird ( Sialia mexicana ) is a secondary cavity-nesting passerine that nests in coniferous and open deciduous forests. Nest box programs for cavity-nesting species have provided suitable nesting locations and opportunities for data collection on nestling growth and development. Methods We developed models for predicting the age of Western Bluebird nestlings from morphometric measurements using model training and validation. These were developed for mass, tarsus, and two different culmen measurements. Results Our models were accurate to within less than a day, and each model worked best for a specific age range. The mass and tarsus models can be used to estimate the ages of Western Bluebird nestlings 0–10 days old and were accurate to within 0.5 days for mass and 0.7 days for tarsus. The culmen models can be used to estimate ages of nestlings 0–15 days old and were also accurate to within less than a day. The daily mean, minimum, and maximum values of each morphometric measurement are provided and can be used in the field for accurate nestling age estimations in real time. Conclusions The model training and validation procedures used here demonstrate that this method can create aging models that are highly accurate. The methods can be applied to any passerine species provided sufficient nestling morphometric data are available.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Tractable minor-free generalization of planar zero-field Ising models

In this work, we present a new family of zero-field Ising models over N binary variables/spins obtained by consecutive 'gluing' of planar and O(1)-sized components and subsets of at most three vertices into a tree. The polynomial time algorithm of the dynamic programming type for solving exact inference (computing partition function) and exact sampling (generating i.i.d. samples) consists of sequential application of an efficient (for planar) or brute-force (for O(1)-sized) inference and sampling to the components as a black box. To illustrate the utility of the new family of tractable graphical models, we first build a polynomial algorithm for inference and sampling of zero-field Ising models over K 33 -minor-free topologies and over K 5 -minor-free topologies—both of which are extensions of the planar zero-field Ising models—which are neither genus- nor treewidth-bounded. Second, we empirically demonstrate an improvement in the approximation quality of the NP-hard problem of inference over the square-grid Ising model in a node-dependent nonzero 'magnetic' field.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Improving Convection Trigger Functions in Deep Convective Parameterization Schemes Using Machine Learning

Abstract Deficiencies in convection trigger functions, used in deep convection parameterizations in General Circulation Models (GCMs), have critical impacts on climate simulations. A novel convection trigger function is developed using the machine learning (ML) classification model XGBoost. The large‐scale environmental information associated with convective events is obtained from the long‐term constrained variational analysis forcing data from the Atmospheric Radiation Measurement (ARM) program at its Southern Great Plains (SGP) and Manaus (MAO) sites representing, respectively, continental mid‐latitude and tropical convection. The ML trigger is separately trained and evaluated per site, and jointly trained and evaluated at both sites as a unified trigger. The performance of the ML trigger is compared with four convective trigger functions commonly used in GCMs: dilute convective available potential energy (CAPE), undilute CAPE, dilute dynamic CAPE (dCAPE), and undilute dCAPE. The ML trigger substantially outperforms the four CAPE‐based triggers in terms of the F 1 score metric, widely used to estimate the performance of ML methods. The site‐specific ML trigger functions can achieve, respectively, 91% and 93% F 1 scores at SGP and MAO. The unified trigger also has a 91% F 1 score, with virtually no degradation from the site‐specific training, suggesting the potential of a global ML trigger function. The ML trigger alleviates a GCM deficiency regarding the overprediction of convection occurrence, offering a promising improvement to the simulation of the diurnal cycle of precipitation. Furthermore, to overcome the black box issue of the ML methods, insights derived from the ML model are discussed, which may be leveraged to improve traditional CAPE‐based triggers.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

In with the old, in with the new: machine learning for time to event biomedical research

The predictive modeling literature for biomedical applications is dominated by biostatistical methods for survival analysis, and more recently some out of the box machine learning approaches. In this article, we show a presentation of a machine learning method appropriate for time-to-event modeling in the area of prostate cancer long-term disease progression. Using XGBoost adapted to long-term disease progression, we developed a predictive model for 118 788 patients with localized prostate cancer at diagnosis from the Department of Veterans Affairs (VA). Our model accounted for patient censoring. Harrell’s c-index for our model using only features available at the time of diagnosis was 0.757 95% confidence interval [0.756, 0.757]. Our results show that machine learning methods like XGBoost can be adapted to use accelerated failure time (AFT) with censoring to model long-term risk of disease progression. Furthermore, the long median survival justifies and requires censoring. Overall, we show that an existing machine learning approach can be used for AFT outcome modeling in prostate cancer, and more generally for other chronic diseases with long observation times.

96 KNOWLEDGE MANAGEMENT AND PRESERVATION↗