Search NASA⌕ Search

SEARCH · Search NASA

Results for “Monte Carlo Analysis”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Kraft black liquor concentration with graphene oxide membranes: Process simulations and technoeconomic analysis

Abstract Black liquor (BL) dewatering by multi‐effect evaporation in the kraft papermaking process is highly energy‐intensive. It was previously shown that graphene oxide (GO) nanofiltration membranes can remove lignin, other organics, and inorganic salts from BL while exhibiting stability in caustic BL conditions. Here, we design and simulate several candidate dewatering processes and evaluate their technoeconomic characteristics. All processes concentrate BL from 15 to 30 wt% solutes while producing aqueous permeate. Two process options were analyzed—option A including “last‐mile” permeate treatment to reduce solutes to 0.2 wt%, and option B excluding this treatment and producing a 3–4 wt% solutes stream. These processes were simulated in custom‐built ASPEN Plus flowsheets interfaced with Microsoft Excel and MATLAB. All processes deliver large (>40%) energy savings. Detailed technoeconomic analysis showed that option A processes are profitable in mills equipped with condensing turbines, but unprofitable with only purchased fuel savings. Option B processes are profitable in both situations, but require the caustic permeate to be utilized in other kraft process units. They are also profitable with electricity generation when operated at smaller scales matching the requirements of other process units. Monte‐Carlo sensitivity analysis shows that Option A can yield median 20‐year NPVs up to ~$10MM and Option B up to ~$25MM. Overall, GO membrane‐based BL dewatering is economically promising, assuming successful slipstream piloting and scale‐up campaigns. It would have immediate sustainability benefits from large energy savings, and broader implications for biorefinery processes due to the ability to fractionate biomass feedstock components under harsh conditions.

Wang, Zhongzhen↗
DOE OSTI2022

A Hydro-based MCMC Analysis of SNR 0509–67.5: Revealing the Explosion Properties from Fluid Discontinuities Alone

Using Hubble Space Telescope/ACS Hα images of SNR 0509–67.5 taken ~10 yr apart, we measure the forward shock (FS) proper motions (PMs) at 231 rim locations. The average shock radius and velocity are 3.66 ± 0.036 pc and 6315 ± 310 km s –1 . Hydrodynamic simulations, recast as similarity solutions, provide models for the supernova remnant’s expansion into a uniform ambient medium. These are coupled to a Markov chain Monte Carlo (MCMC) analysis to determine explosion parameters, constrained by the FS measurements. For our baseline model, the MCMC posteriors yield an age of 315.5 ± 1.8 yr, a dynamical explosion center at 5 h 09 m 31. s 16, $-67^\circ 31^{\prime} 17\buildrel{\prime\prime}\over{.} 1$ and ambient medium densities at each azimuth ranging over 3.7–8.0 × 10 –25 g cm –3 . The age uncertainty can be an order of magnitude larger when considering other models, or subsets of the data. We detect stellar PMs corresponding to speeds in the Large Magellanic Cloud ≥ 770 km s –1 . Five stars in the remnant show measurable PMs but none are moving radially from the dynamical center. There are four stars 1."4 from the center, including three faint, previously unidentified ones. Using coronal [Fe XIV ] λ5303 emission as a proxy for the reverse shock location, we constrain the explosion energy (for a compression factor of 4) to a value of E = (1.30 ± 0.41) × 10 51 erg for the first time from shock kinematics alone. Higher compression factors (7 or more) are strongly disfavored based on multiple criteria, arguing for inefficient particle acceleration in the Balmer shocks of SNR 0509–67.5.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

A Fundamental Economic Assessment of Recovering Rare Earth Elements and Critical Minerals from Acid Mine Drainage Using a Network Sourcing Strategy

In recent years, acid mine drainage (AMD) has emerged as a promising unconventional source of rare earth elements (REEs) and other critical minerals (CMs) such as cobalt and manganese. In this regard, AMD provides a natural heap leaching effect that extracts and concentrates REE/CM from the host strata creating a partially enriched feedstock suitable for downstream extraction, separation, and recovery. While several prior studies have described processes and approaches for the valorization of AMD, very few have described the supply chain and infrastructure requirements as well as the associated economic assessment. To that end, this paper provides a fundamental economic assessment of REE/CM recovery from AMD using a network sourcing strategy in addition to a robust, flexible feedstock separations and refining facility. The methodology of this paper follows that of a typical techno-economic analysis with capital and operating costs estimated using AACE Class IV (FEL-2) guidelines. To demonstrate the range of possible outcomes, four pricing scenarios were modeled including contemporary prices (September, 2021) as well as the minimum and maximum prices over the last decade. In addition, five production scenarios were considered reflecting variations in the product suite, ranging from full elemental separation to magnet REE and CM production only (i.e., Pr, Nd, Tb, Dy, Y, Sc, Co, and Mn). The results of this analysis show that, with the exception of the minimum price scenario, all operational configurations have positive economic indicators with rates of return varying from 25% to 32% for the contemporary price scenario. The optimal configuration was determined to be production of Co, Mn, and all REEs except for mischmetal, which is not recovered. Sensitivity analysis and Monte Carlo simulation show that capital cost and HCl consumption are the two major factors influencing rate of return, thus indicating opportunities for future technology development and cost optimization. Implications of the study and a cooperative profit-sharing model for sourcing are also described.

42 ENGINEERING↗
DOE OSTI2025

Developing A Custom-Built Metal Aerosol Processing Chamber: Analysis of Aerosol Coagulation at Low Humidities

We have developed an intermediate size (906 L) aerosol processing chamber, and this work reports on the design and initial characterization of dry aerosol experiments. Specifically, we are determining wall-loss and coagulation correction factors using the observed size distribution measurements for surrogates of common aerosol classes: sodium chloride, sucrose, and biomass burning aerosol smoke. Results show that, on average, sodium chloride, sucrose, and smoke wall-loss rates converge to similar values on relatively short time scales (< 1 h). The fitted coagulation correction factor, W C -1 , for smoke particles (1.23 ± 0.312), indicates that on average they adhere to each other more than sodium chloride (0.969 ± 0.524) and sucrose (1.16 ± 1.38). The relative uncertainty is high for the coagulation correction, but it is consistent with our Monte Carlo error analysis. This study lays the foundation for future experiments at elevated humidity and supersaturation conditions to characterize the influence of particle shape on coagulation and cloud parameters.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

TIC 168789840: A Sextuply Eclipsing Sextuple Star System

We report the discovery of a sextuply eclipsing sextuple star system from TESS data, TIC 168789840, also known as TYC 7037-89-1, the first known sextuple system consisting of three eclipsing binaries. The target was observed in Sectors 4 and 5 during Cycle 1, with lightcurves extracted from TESS Full Frame Image data. It was also previously observed by the WASP survey and ASAS-SN. The system consists of three gravitationally bound eclipsing binaries in a hierarchical structure of an inner quadruple system with an outer binary subsystem. Follow-up observations from several different observatories were conducted as a means of determining additional parameters. The system was resolved by speckle interferometry with a 0.″42 separation between the inner quadruple and outer binary, inferring an estimated outer period of ∼2 kyr. It was determined that the fainter of the two resolved components is an 8.217 day eclipsing binary, which orbits the inner quadruple that contains two eclipsing binaries with periods of 1.570 days and 1.306 days. Markov Chain Monte Carlo (MCMC) analysis of the stellar parameters has shown that the three binaries of TIC 168789840 are “triplets,” as each binary is quite similar to the others in terms of mass, radius, and T {sub eff}. As a consequence of its rare composition, structure, and orientation, this object can provide important new insight into the formation, dynamics, and evolution of multiple star systems. Future observations could reveal if the intermediate and outer orbital planes are all aligned with the planes of the three inner eclipsing binaries.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2022

Skewing the CMB×LSS: a fast method for bispectrum analysis

Upcoming cosmic microwave background (CMB) lensing measurements and tomographic galaxy surveys are expected to provide us with high-precision data sets in the coming years, thus paving the way for fruitful cross-correlation analyses. Here, in this paper we study the information content of the weighted skew-spectrum, a nearly-optimal estimator of the angular bispectrum amplitude, as a means to extract non-Gaussian information on both bias and cosmological parameters from the bispectra of galaxies cross-correlated with CMB lensing, while gaining significantly on speed. Our results show that for the combination of the Planck satellite and the Dark Energy Spectroscopic Instrument (DESI), the difference in the constraints on bias and cosmological parameters from the skew-spectrum and the bispectrum is at most 17%. We further compare and find agreement between our theoretical skew-spectra and those estimated from N-body simulations, for which it is important to include gravitational non-linearities beyond perturbation theory and the post-Born effect for CMB lensing. We define an algorithm to apply the skew-spectrum estimator to the data and, as a preliminary step, we use the skew-spectra to constrain bias parameters and the amplitude of shot noise from the simulations through a Markov chain Monte Carlo likelihood analysis, finding that it may be possible to reach percent-level estimates for the linear bias parameter b 1 .

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Diffuse back-illumination temperature imaging (DBI-TI), a novel soot thermometry technique

To meet stringent emissions regulations on soot emissions, it is critical to further advance the fundamental understanding of in-cylinder soot formation and oxidation processes. Among several optical techniques for soot quantification, diffuse back-illumination extinction imaging (DBI-EI) has recently gained traction mainly due to its ability to compensate for beam steering, which if not addressed, can cause unacceptably high measurement uncertainty. Until now, DBI-EI has only been used to measure the amount of soot along the line of sight, and here, we extend the capabilities of a DBI-EI setup to also measure in-cylinder soot temperature. This proof of concept of diffuse back-illumination temperature imaging (DBI-TI) as a soot thermometry technique is presented by implementing DBI-TI in a single cylinder, heavy-duty, optical diesel engine to provide 2-D line-of-sight integrated soot temperature maps. The potential of DBI-TI to be an accurate thermometry technique for use in optical engines is analyzed. The achievable accuracy is due in part to simultaneous measurement of the soot extinction, which circumvents the uncertainty in dispersion coefficients that depend on the optical properties of soot and the wavelength of light utilized. Analysis shows that DBI-TI provides temperature estimates that are closer to the mass-averaged soot temperature when compared to other thermometry techniques that are more sensitive to soot temperature closer to the detector. Furthermore, uncertainty analysis and Monte Carlo (MC) simulations provide estimates of the temperature measurement errors associated with this technique. The MC simulations reveal that for the light intensities and optical densities encountered in these experiments, the accuracy of the DBI-TI technique is comparable or even better than other established optical thermometry techniques. Thus, DBI-TI promises to be an easily implementable extension to the existing DBI-EI technique, thereby extending its ability to provide comprehensive line-of-sight integrated information on soot.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Adsorption of Natural Gas in Metal–Organic Frameworks: Selectivity, Cyclability, and Comparison to Methane Adsorption

Evaluation of metal–organic frameworks (MOFs) for adsorbed natural gas (ANG) technology employs pure methane as a surrogate for natural gas (NG). This approximation is problematic, as it ignores the impact of other heavier hydrocarbons present in NG, such as ethane and propane, which generally have more favorable adsorption interactions with MOFs compared to methane. Herein, using quantitative Raman spectroscopic analysis and Monte Carlo calculations, we demonstrate the adsorption selectivity of high-performing MOFs, such as MOF-5, MOF-177, and SNU-70, for a methane and ethane mixture (95:5) that mimics the composition of NG. The impact of selectivity on the storage and deliverable capacities of these adsorbents during successive cycles of adsorption and desorption, simulating the filling and emptying of an ANG tank, is also demonstrated. Finally, the study reveals a gradual reduction in the storage performance of MOFs, particularly with smaller pore volumes, due to ethane accumulation over long-term cycling, until a steady state is reached with substantially degraded storage performance.

03 NATURAL GAS↗
DOE OSTI2026

Framework Short-Range Order Observed in a Spinel-Type Li Superionic Conductor

Solid-state superionic conductors are characterized by rich structural disorders. Though structural complexities are central to their functionalities, they often give rise to short-range order that eludes detection by conventional diffraction-based techniques and is thus overlooked in establishing precise structure-property relationships. In this work, we synthesized single crystals of a recently discovered lithium (Li) superionic conductor Li16.2(1)In9.00(2)Sn1.10(1)O23.8 (LISO) for in-depth characterizations of structural subtleties. LISO exhibits an unusual spinel-like phase with significant Li overstoichiometry and a face-sharing Li network. Single-crystal neutron diffraction confirms significant Li disorder, as manifested in Li site splitting and partial occupancy. More importantly, synchrotron diffuse scattering combined with 3D-ΔPDF analysis and Monte Carlo simulations reveal short-range order in the nonalkali framework that might contribute to the phase stability and ionic conductivity. This work showcases an example in which subtle local energetics can be directly visualized in structurally disordered ionic conductors.

Chen, Yu↗
DOE OSTI2021

Photoinduced anisotropic lattice dynamic response and domain formation in thermoelectric SnSe

Identifying and understanding the mechanisms behind strong phonon–phonon scattering in condensed matter systems is critical to maximizing the efficiency of thermoelectric devices. To date, the leading method to address this has been to meticulously survey the full phonon dispersion of the material in order to isolate modes with anomalously large linewidth and temperature-dependence. Here we combine quantitative MeV ultrafast electron diffraction (UED) analysis with Monte Carlo based dynamic diffraction simulation and first-principles calculations to directly unveil the soft, anharmonic lattice distortions of model thermoelectric material SnSe. A small single-crystal sample is photoexcited with ultrafast optical pulses and the soft, anharmonic lattice distortions are isolated using MeV-UED as those associated with long relaxation time and large displacements. We reveal that these modes have interlayer shear strain character, induced mainly by c -axis atomic displacements, resulting in domain formation in the transient state. These findings provide an innovative approach to identify mechanisms for ultralow and anisotropic thermal conductivity and a promising route to optimizing thermoelectric devices.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

First principles reaction discovery: from the Schrodinger equation to experimental prediction for methane pyrolysis

Our recent success in exploiting graphical processing units (GPUs) to accelerate quantum chemistry computations led to the development of the ab initio nanoreactor, a computational framework for automatic reaction discovery and kinetic model construction. In this work, we apply the ab initio nanoreactor to methane pyrolysis, from automatic reaction discovery to path refinement and kinetic modeling. Elementary reactions occurring during methane pyrolysis are revealed using GPU-accelerated ab initio molecular dynamics simulations. Subsequently, these reaction paths are refined at a higher level of theory with optimized reactant, product, and transition state geometries. Reaction rate coefficients are calculated by transition state theory based on the optimized reaction paths. The discovered reactions lead to a kinetic model with 53 species and 134 reactions, which is validated against experimental data and simulations using literature kinetic models. We highlight the advantage of leveraging local brute force and Monte Carlo sensitivity analysis approaches for efficient identification of important reactions. Both sensitivity approaches can further improve the accuracy of the methane pyrolysis kinetic model. The results in this work demonstrate the power of the ab initio nanoreactor framework for computationally affordable systematic reaction discovery and accurate kinetic modeling.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2026

Dark matter clumps as sources of gravitational-wave glitches in LIGO-Virgo-KAGRA data

We consider the hypothetical possibility that nonstationary glitch features in the noise of ground-based gravitational-wave detectors could be produced by small dark matter clumps that pass through the earth in the vicinity of gravitational-wave detectors. Here, we first derive the gravitational-wave strain that would be generated by the passage of such a dark matter clump. We find that the strain is primarily sourced by the Newtonian gravitational acceleration of the mirrors toward the clump and by the Shapiro time delay of the photons in the laser beams as they pass through the gravitational potential created by the dark matter clump. We also find that the Newtonian acceleration effect dominates the gravitational-wave strain for both ground and space-based interferometers. We then compare our dark matter clump, gravitational-wave strain model to 84 Koi-Fish glitches detected during the second observing run of the LIGO/Virgo/KAGRA Collaboration through a Markov chain Monte Carlo Bayesian analysis. We find that all glitches but nine can be confidently rejected as having originated from dark matter clumps. For the remaining glitches, the dark matter hypothesis cannot be excluded, and the maximum a posteriori parameters yield minimum densities of about 10 −7 g⁡/cm 3 , within the model. These results allow us to place the first direct upper limits with gravitational-wave detectors on the local over-density of dark matter in the form of clumps in the local neighborhood of Earth, namely 𝜌 DM ⁢clumps ≲ 10 −15 g⁡/cm −3 .

Astronomy and AstroPhysics↗
DOE OSTI2024

CMOS-Based Single-Cycle in-Memory XOR/XNOR

Big data applications are on the rise, and so is the number of data centers. The ever-increasing massive data pool needs to be periodically backed up in a secure environment. Moreover, a massive amount of securely backed-up data is required for training binary convolutional neural networks for image classification. XOR and XNOR operations are essential for large-scale data copy verification, encryption, and classification algorithms. The disproportionate speed of existing compute and memory units makes the von Neumann architecture inefficient to perform these Boolean operations. Compute-in-memory (CiM) has proved to be an optimum approach for such bulk computations. The existing CiM-based XOR/XNOR techniques either require multiple cycles for computing or add to the complexity of the fabrication process. Here, we propose a CMOS-based hardware topology for single-cycle in-memory XOR/XNOR operations. Our design provides at least 2× improvement in the latency compared with other existing CMOS-compatible solutions. We verify the proposed system through circuit/system-level simulations and evaluate its robustness using a 5000-point Monte Carlo variation analysis. This all-CMOS design paves the way for practical implementation of CiM XOR/XNOR at scaled technology nodes.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Bidding Curve Design for Hybrid Power Plants with Uncertain Solar Forecast

This paper presents a novel bidding curve design algorithm tailored for hybrid power plants (HPPs) to participate in the wholesale electricity market. Utilizing forecasts for photovoltaic (PV) generation and available battery power, our algorithm strategically computes the bidding curve to maximize HPP profit while adeptly managing the inherent uncertainty associated with PV power generation. In addition, the introduction of the penalty cost in HPP bidding curves provides the system operator a tool to effectively manage the system-level uncertainty that caused by HPPs. Numerical analysis through Monte Carlo simulations confirms that our bidding curve methodology outperforms the benchmark across various scenarios.

bidding curve↗
DOE OSTI2022

A Sparse and Low Rank Penalized Signal Decomposition Model with Constraints: Anomaly Detection in PV Systems

Recently, robust PCA has seen its wide application in various industries for its ability to perform the task of anomaly detection. The essence of robust PCA approach is to break down the signal into a low rank component and sparse component. In many applications, a simple breakdown of the signal without accounting for the signs of low rank components and sparse components would violate the physical constraints of the decomposed signal. In addition, often times, the signals in the real world collected for a long duration has smooth changes within a day and between days. As an example, the power signals collected in a photovoltaic (PV) system are cyclostationary, exhibiting these characteristics. Neglecting the smoothness of signals would result in miss detection of anomalous signals which are smooth within a day but non-smooth between days and vice versa. In this paper, we developed a signal decomposition approach for the purpose of anomaly detection based on the idea of low rank and sparse decomposition taking into consideration the signs of the decomposed low rank and sparse components and the within-day and between-day smooth changes in the original signals. The proposed unsupervised approach for fault detection eliminates the need for faulty samples required by other machine learning methods. It does not require the full I-V characteristics to work. Furthermore, there is no need for complex modelling of PV systems as in the case of power loss analysis. Using Monte Carlo simulations, we demonstrate the ability of our proposed approach for detecting anomalies of different duration and severity in PV systems.

14 SOLAR ENERGY↗
DOE OSTI2026

Efficient Probabilistic Visualization of Local Divergence of 2D Vector Fields with Independent Gaussian Uncertainty

This work focuses on visualizing uncertainty of local divergence of two-dimensional vector fields. Divergence is one of the fundamental attributes of fluid flows, as it can help domain scientists analyze potential positions of sources (positive divergence) and sinks (negative divergence) in the flow. However, uncertainty inherent in vector field data can lead to erroneous divergence computations, adversely impacting downstream analysis. While Monte Carlo (MC) sampling is a classical approach for estimating divergence uncertainty, it suffers from slow convergence and poor scalability with increasing data size and sample counts. Thus, we present a two-fold contribution that tackles the challenges of slow convergence and limited scalability of the MC approach. (1) We derive a closed-form approach for highly efficient and accurate uncertainty visualization of local divergence, assuming independently Gaussian-distributed vector uncertainties. (2) We further integrate our approach into Viskores, a platform-portable parallel library, to accelerate uncertainty visualization. In our results, we demonstrate significantly enhanced efficiency and accuracy of our serial analytical (speed-up up to 1946×) and parallel Viskores (speed-up up to 19698×) algorithms over the classical serial MC approach. We also demonstrate qualitative improvements of our probabilistic divergence visualizations over traditional mean-field visualization, which disregards uncertainty. We validate the accuracy and efficiency of our methods on wind forecast and ocean simulation datasets.

Ouermi, Timbwaoga [University of Utah]↗