SEARCH · Engineering Papers
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.
HyRAM is a prototype software toolkit that integrates data and methods relevant to assessing the safety of hydrogen fueling and storage infrastructure. The HyRAM toolkit integrates deterministic and probabilistic models for quantifying accident scenarios, predicting physical effects, and characterizing the impact of hydrogen hazards, including thermal effects from jet fires and thermal pressure effects from deflagration.
In remote sensing systems, the capabilities of the system are constrained by the complex interactions between size, weight, and power (SWAP) of potential designs. In electro-optical (EO) systems, examples of these critical parameters include the system’s sensitivity and resolution. Those parameters can be increased by ever larger optical apertures and focal planes but at the cost of more SWAP. Multi-image super resolution (MISR) techniques allow resolution to be enhanced via computation rather than more sophisticated optical hardware. These algorithms combine multiple images together into a single, higher resolution image, trading temporal resolution and computation for spatial resolution. Fielded MISR techniques, such as Drizzle, can require several hundred images to create a single super resolved image, implying reduced temporal resolution, increased data acquisition load, and limiting mission applications. Iterative techniques, such as model-based image reconstruction and compressive sensing, have been shown to create super resolved images using fewer images than Drizzle. They do this by posing an optimization problem that balances accuracy between a highly accurate physical model and an image model. In the case of super resolution, the physical model is defined by the relation between low resolution input images and the desired high resolution output image. The image model encodes some assumptions about the super resolved image. These assumptions are meant to suppress reconstruction artifacts that arise due to deterministic physical model error, stochastic measurement noise, and potential undersampling. In practice, the performance of iterative methods are limited by imaging models compatible with optimization. Deep learning-based methods can effectively learn image models of arbitrary complexity, but lack the theoretical explainability and robustness of iterative techniques. Consensus equilibrium (CE) generalizes the iterative techniques beyond optimization, enabling blackbox algorithms such as traditional and neural image denoisers to be used as the image model. CE-based approaches retain much of the explainability and robustness of iterative techniques while allowing the expressiveness of machine learning image models to be used. Additionally, by unrolling iterations of CE with an embedded image denoiser, the image denoiser can be further trained and specialized to the specific application with potentially higher quality reconstructions. Under this project, we demonstrated the feasibility of training an unrolled neural network based upon CE. While we didn’t train one, we showed that the CE process is differentiable and its gradient can be tractably computed. We also explored the usage of a variants of CE akin to generative neural works. Most importantly, we applied the CE framework to a number of problems including non-blind deconvolution, upsampling, single-image super resolution, MISR, event-based sensing, and saturated deconvolution. Our MISR prototype creates high quality reconstructions with an order of magnitude fewer images than previous approaches and, critically, produces these reconstructions fast enough for practical usage.
Machine learning technology spans many areas and today plays a significant role in addressing a wide range of problems in critical domains,i.e., healthcare, autonomous driving, finance, manufacturing, cybersecurity,etc. Metamorphic testing (MT) is considered a simple but very powerful approach in testing such computationally complex systems for which either an oracle is not available or is available but difficult to apply. Conventional metamorphic testing techniques have certain limitations in verifying deep learning-based models (i.e., convolutional neural networks (CNNs)) that have a stochastic nature (because of randomly initializing the network weights) in their training. In this article, we attempt to address this problem by using a statistical metamorphic testing (SMT) technique that does not require software testers to worry about fixing the random seeds (to get deterministic results) to verify the metamorphic relations (MRs). We propose seven MRs combined with different statistical methods to statistically verify whether the program under test adheres to the relation(s) specified in the MR(s). We further use mutation testing techniques to show the usefulness of the proposed approach in the healthcare space and test two CNN-based deep learning models (used for pneumonia detection among patients). The empirical results show that our proposed approach uncovers 85.71% of the implementation faults in the classifiers under test (CUT). Furthermore, we also propose an MRs minimization algorithm for the CUT, thus saving computational costs and organizational testing resources.
HyRAM is a software toolkit that integrates data and methods relevant to assessing the safety of hydrogen fueling and storage infrastructure. The HyRAM toolkit integrates deterministic and probabilistic models for quantifying accident scenarios, predicting physical effects, and characterizing the impact of hydrogen hazards (thermal effects from jet fires, thermal pressure effects from deflagrations) on people and structures. HyRAM incorporates generic probabilities for equipment failures for nine types of components, and probabilistic models for the impact of heat flux on humans and structures, with computationally and experimentally validated models of hydrogen release and flame physics. Version 1.0.0.280 can be used to quantify the likelihood and thermal consequences associated with gaseous hydrogen releases from user-defined hydrogen installations.
Multiphysics reactor analysis has become one of the leading activities of most nuclear engineering research groups. Coupled neutronics and thermal/hydraulics (T/H) simulations have traditionally been performed by adopting deterministic transport tools and reactor system T/H codes. The main reason for avoiding higher fidelity numerical methods like continuous-energy Monte Carlo analysis has been the high computational cost of coupled, full-scale simulations. Nowadays, high fidelity multiphysics simulations represent a hot research topic, with several applications involving advanced nuclear systems. Liquid fuel Molten Salt Reactors (MSRs) feature peculiar physical phenomena that are not present in solid fuel reactors. Some of these phenomena are not correctly resolved by legacy reactor physics tools, and require specific treatments. This work focuses in particular on the analysis of the effect of the delayed neutron precursors drift and the impact of liquid fluid thermal expansion and compressibility in postulated criticality excursion transients in fast-spectrum MSRs. A new modeling tool based on a Monte Carlo approach for neutronics (Serpent code) coupled to a Computational Fluid Dynamic (CFD) analysis (OpenFOAM) was developed to obtain accurate modeling and simulation of these phenomena. The development of the Serpent Monte Carlo code [1] at VTT in the past few years has made this code particularly attractive for performing multiphysics studies. The availability of a flexible multiphysics interface, on-the-fly Doppler broadening and density-based collision rejection, and Delta-tracking, proved of paramount importance for the computationally efficient performance of coupled Monte Carlo/CFD calculations [2]. OpenFOAM is a C++ multiphysics toolkit featured by automatic matrix construction and solution capabilities for scalar and vector partial differential equations (PDEs), adopting several possible differencing and interpolation schemes, based on the Finite-Volume method. The ease of implementing the constitutive equations of different physical phenomena makes OpenFOAM particularly suitable for the modeling of complex coupled phenomena. In the past, several research activities proposed the adoption of both Serpent and OpenFOAM for the modeling and simulation of MSRs and the associated fuel cycles. In this work, these tools have been internally coupled for the steady-state and transient analysis of fast-spectrum systems, accounting for the effect of delayed neutron precursors transport and the salt thermal expansion and compressibility.
Cancer treatment using radiotherapy techniques needs an accurate prediction of the dose delivered to the patient at each position. Treatment planning systems (TPS) allow the optimization of the dose delivered to the tumor volume while minimizing the dose absorbed by the surrounding healthy tissue. Commercial treatment planning systems use empirical deterministic algorithms that do not correctly estimate the dose in the surroundings of heterogeneities. This lack of accuracy is taken into account as bigger security margins in the treatment, reducing the potential effect of the treatment. However, Monte Carlo (MC) simulation techniques have proved to predict the dose distribution more accurately and also representing experimental data better. Nevertheless, the clinical application of these techniques has been limited by the long calculation times they require. As high performance computers decrease in cost MC simulations become more feasible to be used in TPS. The present work improves a simplified model of the MLC-Elekta Precise linear accelerator (linac or unit from now on) previously developed and correctly validated by the group. A multileaf collimator (MLC) has been modelled according to the manufacturer specifications and added to the previous model to allow higher conformity in the studies. The full geometry has been entered to the MCNP6 Monte Carlo code to simulate the transport of particles through the linac. The geometry also includes a cubic water phantom for dosimetry studies and validation with experimental data. The experimental data employed for the dose distribution comparison has been measured at the linac from the linac Hospital Universitari Clinic de Valencia (Spain). Additionally, particles are stored in a phase space surface before their entrance in the water to be analyzed and for performance optimization. Given the results presented, we conclude that this first approach to validate a model of the Elekta Precise linac from the Hospital Cl nic Universitari de Valancia represents well the depth-dose distribution, but fails at representing the dose profiles. It has been shown that they predict too much penumbra, probably caused by the filtration between the leaves. In order to solve this problems, further work is needed to improve the simulated model. This work will mainly focus on the multileaf collimator. Also, the codes developed by the group to analyze phase space files produced by MCNP6 have proven highly useful. As work progress and the codes are debugged, they will be published so that the community can take advantage of them.
In this work, we develop a probabilistic estimator for the voltage-to-current map arising in electrical impedance tomography. This novel so-called partially reflecting random walk on spheres estimator enables Monte Carlo methods to compute the voltage-to-current map in an embarrassingly parallel manner, which is an important issue with regard to the corresponding inverse problem. Our method uses the well-known random walk on spheres algorithm inside subdomains where the diffusion coefficient is constant and employs replacement techniques motivated by finite difference discretization to deal with both mixed boundary conditions and interface transmission conditions. We analyze the global bias and the variance of the new estimator both theoretically and experimentally. Subsequently, the variance of the new estimator is considerably reduced via a novel control variate conditional sampling technique which yields a highly efficient hybrid forward solver coupling probabilistic and deterministic algorithms.
Scientific computations are expected to be increasingly distributed across wide-area networks, and Message Passing Interface (MPI) has been shown to scale to support their communications over long distances. Application-level measurements of MPI operations reflect the connection Round-Trip Time (RTT) and loss rate, and machine learningmethods have been previously developed to estimate them under deterministic periodic losses. In this paper, we consider more complex, random losses with uniform, Poisson and Gaussian distributions. We study five disparate machine leaning methods, with linear and non-linear, and smooth and non-smooth properties, to estimate RTT and loss rate over 10 Gbps connections with 0–366 ms RTT. The diversity and complexity of these estimators combined with the randomness of losses and TCP’s non-linear response together rule out the selection of a single best among them; instead, we fuse them to retain their design diversity. Overall, the results show that accurate estimates can be generated at low loss ratesbut become inaccurate at loss rates 10% and higher, thereby illustrating both their strengths and limitations.
The increasing penetration of renewable energy sources and the retirement of conventional generation units have decreased system inertia, making power systems more vulnerable to resilience and stability issues. To address this problem, this paper proposes a novel approach using grid-supportive loads (GSLs) to provide a fast and concise primary frequency response and a deep deterministic policy gradient agent-based secondary controller to restore the system frequency to the nominal value. The proposed method is evaluated on the single-area and multi-area test systems. The simulation results demonstrate that using GSLs enhances the power system's stability and resilience. Compared to conventional controllers, the frequency nadir is improved with GSLs. Additionally, the proposed method effectively enhances resilience even with high penetration. These findings indicate that the proposed approach can improve the resilience and stability of power systems and provide a promising solution for future power systems. The results of this study underscore the importance of utilizing innovative approaches to enhance the stability and resilience of power systems in the context of high penetration of renewable energy sources and the retirement of conventional generation.
Bayesian inverse modeling techniques are computationally expensive because many forward simulations are needed when sampling the posterior distribution of the parameters. Here, in this paper, we combine the implicit sampling method and generalized polynomial chaos expansion (gPCE) to significantly reduce the computational cost of performing Bayesian inverse modeling. There are three steps in this approach: (1) find the maximizer of the likelihood function using deterministic approaches; (2) construct a gPCE-based surrogate model using the results from a limited number of forward simulations; and (3) efficiently sample the posterior distribution of the parameters using implicit sampling method. The cost of constructing the gPCE-based surrogate model is further decreased by using sparse Bayesian learning to reduce the number of gPCE coefficients that have to be determined. We demonstrate the approach for a synthetic ponded infiltration experiment simulated with TOUGH2. The surrogate model is highly accurate with mean relative error that is < 0.035% in predicting saturation and < 0.25% in predicting the likelihood function. The posterior distribution of the parameters obtained using our proposed technique is nearly indistinguishable from the results obtained from either an implicit sampling method or a Markov chain Monte Carlo method utilizing the full mode.
Calculation of the shutdown dose rate is crucial for safe fusion reactor operations. The Rigorous-two-step (R2S) method is a method that requires connected neutron transport, activation, and gamma transport. Shift has integrated variance reduction with a deterministic solver Denovo, supports multiple geometry formats, and is scalable. These features make it an attractive transport solver choice for an R2S workflow. An R2S workflow for the Shift Monte Carlo code is developed and compared to the existing Oak Ridge National Laboratory Shutdown Dose Rate Code Suite (ORCS) workflow. Also, a Python framework for integrating two R2S workflows is developed to mix and match each step in the R2S workflow for improved collaboration and verification experience. In this study, results show that the Shift-Denovo R2S workflow and the ORCS workflow calculate the shutdown dose rate of the ITER Shutdown Dose Rate benchmark problem with an average relative error of 2.285%.
Reliable transfer processes that enable manipulation of two-dimensional (2D) materials, e.g., transition metal dichalcogenides (TMDCs) and MXenes, from one substrate to another has been a necessity for successful device fabrication. Here, with both mechanical exfoliation and chemical vapor deposition (CVD) widely used, a versatile, clean, deterministic, and yet simple transfer technique is highly needed. To address such need, we developed a transfer method that takes advantage of wettability contrast between interfaces without the use of sacrificial layers or chemical processes. More importantly, a setup was developed to carry out this transfer method with high sample selectivity and fine control of the position and orientation of transferred TMDC crystals, a feature required for fabrication of the devices based on vertical 2D heterostructures. Using both exfoliated and CVD grown materials and subsequent atomic force microscopy (AFM), photoluminescence (PL), confocal Raman and tip enhanced Raman spectroscopy (TERS) characterization, we ascertained the quality of interfaces resulting from the transfer process while preserving excellent 2D material integrity. PL and TERS maps revealed nanometer-scale heterogeneities in the interfaces of fabricated heterostructures, which should enable further perfection of the transfer technique. TERS/TEPL information were employed to identify areas suitable for nanodevice fabrication, making the reported transfer and characterization methods ideal for making high quality assembly of 2D heterostructure more accessible, which should facilitate exploration of vertical 2D heterostructures for applications in electronics, batteries, solar cells, and twistronics.
Here, we develop a multistage, stochastic mixed-integer model to support biofuel supply chain expansion under evolving uncertainties. By utilizing the block-separable recourse property, we reformulate the multistage program in an equivalent two-stage program and solve it using an enhanced nested decomposition method with maximal non-dominated cuts. We conduct extensive numerical experiments and demonstrate the application of the model and algorithm in a case study based on the South Carolina settings. The value of multistage stochastic programming method is also explored by comparing the model solution with the counterparts of an expected value based deterministic model and a two-stage stochastic model.
In this paper, we propose a new heuristic symbolic tool for unveiling chaotic and stochastic dynamics: the permutation spectrum test. Several numerical examples allow us to confirm the usefulness of the introduced methodology. Indeed, we show that it is robust in situations in which other techniques fail (intermittent chaos, hyperchaotic dynamics, stochastic linear and nonlinear correlated dynamics, and deterministic non-chaotic noise-driven dynamics). We illustrate the applicability and reliability of this pragmatic method by examining real complex time series from diverse scientific fields. Taking into account that the proposed test has the advantages of being conceptually simple and computationally fast, we think that it can be of practical utility as an alternative test for determinism.
Throughout core life, the reactor vessel is continually bombarded with neutron flux embrittling the vessel and limiting its useful life. As a result, it is imperative to adequately predict how much flux has passed through the vessel during all stages of its life, as a means of ensuring it has adequate ductility to preclude a brittle fracture even under design accident conditions. At present many of these predictions are made either by purely deterministic or purely stochastic means, each offering their own advantages and disadvantages. This paper introduces a coupled method for calculating fluence through the VERA framework, which enables modeling incremental accumulations of fluence while simultaneously accounting for several factors as they evolve through core life. The intention of this method is to simplify the calculation process and provide better resource utilization in obtaining the required accuracy.
We present a randomized dynamical decoupling (DD) protocol that can substantially improve the performance of any given deterministic DD scheme for suppressing coherent noise by using no more than two additional pulses. Our construction is implemented by probabilistically applying sequences of pulses, which, when combined, effectively eliminate the error terms that scale linearly with the system-environment coupling strength. As a result, we show that a randomized protocol using a few pulses can outperform deterministic DD protocols that require considerably more pulses. Furthermore, we prove that the randomized protocol provides an improvement compared to deterministic DD sequences that aim to reduce the error in the system’s Hilbert space, such as Uhrig DD, which had been previously regarded to be optimal. To rigorously evaluate the performance, we introduce new analytical methods suitable for analyzing higher-order DD protocols that might be of independent interest. Here, we also present numerical simulations confirming the significant advantage of using randomized protocols compared to widely used deterministic protocols.
Introducing risk-informed safety analyses for external events can bring sensible benefits to the nuclear power industry easing the regulatory burden introduced in the aftermath of the Fukushima accident. In the framework of the US-DOE Light-Water Reactor Sustainability program, Risk-Informed Margin Characterization project (LWRS-RISMC), research activities were conducted at INL for developing and applying a toolkit for performing risk-informed external hazards analysis. The “EEVE” (External EVEnts) toolkit, a chain of INL and state-of-the-art codes, was set-up and applied for simulating earthquake-induced transients with internal flooding in a PWR. This was done using methodology set-up for evaluating new data and determining critical areas that would benefit most from advanced analysis methods. The methodology is agnostic and it can be applied to other external hazards than just seismic and flooding (e.g., high winds, intense precipitation, etc.). The first step was the deterministic calculations for evaluating the propagation of a set of earthquakes and for determining the effects on a fire suppression system in an auxiliary building. These deterministic calculations were performed using the LS-DYNA and the OPENSEE structural mechanics codes. Probability failures and event trees of a generic 3-loops Westinghouse PWR were derived from a generic SAPHIRE model. Seismically induced Loss of Offsite Power, Station Blackout, and fire suppression piping failure events were considered. SAPHIRE code allowed a pre-screening of all possible branches of the event tree and a selection of those with the largest Core Damage Frequency increase for further analysis using dynamic PRA methods. INL's EMRALD code performed dynamic PRA calculations of the identified sequences, integrating risk analysis with on-line deterministic safety analysis results (seismically-induced internal flooding calculations by NEUTRINO code and system analysis by RELAP5-3D code). RELAP5-3D Best-Estimate calculations coupled with RAVEN statistical code were used for deriving limit surfaces as function of relevant transient parameters. This data was used by EMRALD for determining the final core status (core damage/core safe). Finally, the coupling of seismic, flooding, system thermal-hydraulic, uncertainty quantification and probabilistic risk analysis for an advanced risk-informed external hazards analysis was demonstrated. Results showed the sensible decrease in the level of conservatism achieved using the developed risk-informed methodology, which allowed also improving the evaluation of the failure modes and failure probabilities of different components.
This study illustrates recent computational advances in the application of the Sjöstrand (area), Rossi, and Feynman methods to estimate the effective multiplication factor of a subcritical system driven by an external neutron source. The methodologies introduced in this study have been validated with the experimental results from the KUKA facility of Japan by Monte Carlo (MCNP6 and MCNPX) and deterministic (ERANOS, VARIANT, and PARTISN) codes. When the assembly is driven by a pulsed neutron source generated by a particle accelerator and delayed neutrons are at equilibrium, the Sjöstrand method becomes extremely fast if the integral of the reaction rate from a single pulse is split into two parts. These two integrals distinguish between the neutron counts during and after the pulse period. To conclude, when the facility is driven by a spontaneous fission neutron source, the timestamps of the detector neutron counts can be obtained up to the nanosecond precision using MCNP6, which allows obtaining the Rossi and Feynman distributions.