Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “variational Monte Carlo”

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

A probabilistic creep model incorporating test condition, initial damage, and material property uncertainty

Uncertainty is prevalent in the creep resistance of alloys, where at elevated temperature and low pressure, rupture can range across logarithmic decades. In this study, a probabilistic continuum-damage-mechanics (CDM)-based model is derived to capture the uncertainty of creep resistance. To meet this objective, creep data for alloy 304 Stainless Steel is gathered. A constitutive model, “Sinh”, is calibrated deterministically to determine the statistical variability of the material properties. Three sources of uncertainty are injected into the model: test condition (stress and temperature), initial damage, and material properties. Probabilistic simulations are carried out by (a) calibrating probability distribution functions (pdfs) for each source of uncertainty (b) randomly sampling the pdfs using Monte Carlo methods and (c) executing simulations to replicate the uncertain creep behavior. A sensitivity analysis is performed to evaluate the relative effect of each source of uncertainty. In full probabilistic simulations, the cumulative uncertainty of creep behavior is evaluated. The probabilistic model accurately predicts the creep deformation and rupture of the available experiments. The probabilistic model is validated for interpolation but lacks extrapolation ability. Several future works are proposed to further improve the model.

36 MATERIALS SCIENCE↗

Uncertainty Quantification of Bifacial Performance Modeling

Analysis on uncertainty in the annual energy of PV systems that can be attributed to parameters of particular importance to bifacial PV modules is presented. Monte Carlo simulations are used to evaluate the effect of uncertain module bifaciality factors, module transmission fractions, albedo values, and ground clearance. The analyses cover a wide spectrum of potential PV array archetypes through variation of installation parameters. The results of the Monte Carlo analysis reveal that the uncertainty is largely dependent on albedo uncertainty, but more simulations are needed to identify trends across system archetypes. The simulations are aimed at attributing an annual energy uncertainty factor for bifacial considerations that can be applied in post-processing of project probability of exceedance analysis.

energy modeling↗

Uncertainty Quantification of Bifacial Performance Modeling

Analysis on uncertainty in the annual energy of PV systems that can be attributed to parameters of particular importance to bifacial PV modules is presented. Monte Carlo simulations are used to evaluate the effect of uncertain module bifaciality factors, module transmission fractions, albedo values, and ground clearance. The analyses cover a wide spectrum of potential PV array archetypes through variation of installation parameters. The results of the Monte Carlo analysis reveal that the uncertainty is largely dependent on albedo uncertainty, but more simulations are needed to identify trends across system archetypes. The simulations are aimed at attributing an annual energy uncertainty factor for bifacial considerations that can be applied in post-processing of project probability of exceedance analysis.

bifacial↗

Uncertainty Quantification of Bifacial Performance Modeling: Preprint

Analysis on uncertainty in the annual energy of PV systems that can be attributed to parameters of particular importance to bifacial PV modules is presented. Monte Carlo simulations are used to evaluate the effect of uncertain module bifaciality factors, module transmission fractions, albedo values, and ground clearance. The analyses cover a wide spectrum of potential PV array archetypes through variation of installation parameters. The results of the Monte Carlo analysis reveal that the uncertainty is largely dependent on albedo uncertainty, but more simulations are needed to identify trends across system archetypes. The simulations are aimed at attributing an annual energy uncertainty factor for bifacial considerations that can be applied in post-processing of project probability of exceedance analysis.

energy modeling↗

Bioproducts from high-protein algal biomass: An economic and environmental sustainability review and risk analysis

This report has High-protein algal biomass as an important bio-commodity that has the potential to provide a new source of sustainable protein products. Herein is a critical review that identifies 1) the most relevant sustainability findings related to the processing of proteinaceous algal biomass to higher value protein products and 2) the potential pathways to improve life cycle assessment (LCA) and techno-economic analysis (TEA) metrics, including life-cycle carbon dioxide equivalent (CO 2 eq), life cycle energy, and minimum selling price (MSP) of these products. The critical review of the literature revealed a large variation in model input parameters relating to these metrics. Therefore, a Monte Carlo analysis was conducted to assess the risk associated with these input variations. To understand the uncertainties that propagate into high-protein algae to products’ systems, we reviewed more than 20 state-of-the-art unit operations for algal biomass processing., including cell disruption, protein solubilization, protein precipitation and purification, and protein concentration. We evaluated displacement of proteinaceous products by algal-bioproducts, including ruminant feed, aquaculture feed, protein tablets, and biopolymers and biopolyesters, with prices in the market ranging from 1.9 to 120 $\$$ kg -1 protein. This review realized that the MSP of ruminant and non-ruminant feed ranges from 0.65±0.56 to 2.9±1.1 $\$$ kg -1 protein, and bioplastics’ MSP ranges from 0.97 to 7.0 $\$$ kg -1 protein. Regarding LCA metrics, there is limited research on life cycle energy in proteinaceous biomass concentration and bioproduct systems, reported at 32.7 MJ kg protein -1 , for animal feed displacement. Animal feed emissions in the literature report negative fluxes, representing environmental benefits, as low as -3.7 kgCO 2 eq kg -1 protein and positive fluxes, i.e., global warming potential, as high as 12.8 kgCO 2 eq kg -1 protein. There is limited research on bioplastics life cycle emissions reported at 0.6 kgCO 2 eq kg -1 protein. In general, the studies to date of algae-derived protein bioproducts showed similar life cycle emissions to soybean meals, nylon, polymers, and polystyrenes. Our risk analysis realized that more than 50% of scenarios can result in negative-net life cycle CO 2 eq emissions. This review and risk analysis assess and demonstrate the scenarios that improve economic and environmental sustainability metrics in high-protein algal bioproducts systems.

09 BIOMASS FUELS↗

CoVVVR: Control Variates & Vegas Variance Reduction

This package is a wrapper over the vegas integration package. The control variate variance reduction method is applied to the function when integrated along with the techniques applied in vegas such as importance sampling. To understand control variates, lets first look at how Monte Carlo works.

Scott, JacobL. [Univ. of Kansas, Lawrence, KS (Uni↗

GIGA-Lens: Fast Bayesian Inference for Strong Gravitational Lens Modeling

We present GIGA-Lens: a gradient-informed, GPU-accelerated Bayesian framework for modeling strong gravitational lensing systems, implemented in TensorFlow and JAX. The three components, optimization using multistart gradient descent, posterior covariance estimation with variational inference, and sampling via Hamiltonian Monte Carlo, all take advantage of gradient information through automatic differentiation and massive parallelization on graphics processing units (GPUs). We test our pipeline on a large set of simulated systems and demonstrate in detail its high level of performance. The average time to model a single system on four Nvidia A100 GPUs is 105 s. The robustness, speed, and scalability offered by this framework make it possible to model the large number of strong lenses found in current surveys and present a very promising prospect for the modeling of ${ \mathcal O }({10}^{5})$ lensing systems expected to be discovered in the era of the Vera C. Rubin Observatory, Euclid, and the Nancy Grace Roman Space Telescope.

79 ASTRONOMY AND ASTROPHYSICS↗

Sensitivity of modelled passive margin stratigraphy to variations in sea level, sediment supply and subsidence

Abstract We produced a 10 Myr synthetic stratigraphic section using a forward stratigraphic model that generates marine deltaic stratigraphy over geological timescales. We recursively fit the model using a Bayesian inversion algorithm to test: (1) if it could be accurately reconstructed; (2) if the parameters used to create it could be recovered; and (3) the sensitivity of the model output to given model parameters and the attendant physical processes. The original synthetic stratigraphic section was produced with cyclical sea‐level variations of 40 and 30 m with 2.4 and 10 Myr periods respectively. Sediment was also supplied cyclically, in 2.4 and 10 Myr cycles with amplitudes of 30 and 80 tons/100 kyr, respectively, varying from a mean of 232 tons/100 kyr. Parameter values were sampled to fit the model using a Markov chain Monte Carlo algorithm, resulting in a ±5 m (1σ) variation between the experimental output and the original. Sea level varied by ±7 m (1σ) within the posterior distribution of parameters. As a result, both the 10 Myr and 2.4 Myr sea‐level cycles could be extracted from the original output. The variation in sediment supply was approximately ±38 tons/100 kyr (1σ) and, as a result, only the larger long‐term supply variations could be accurately recovered in refitting the model. The variation in thermal, flexural and total subsidence across those parameter sets is less than ±10 m (1σ). The original section experienced 150 m of total subsidence at the depocentre. Our results demonstrate the distinct and interpretable imprint of sea level and subsidence on continental margin stratigraphy can be quantified. Moreover, we conclude that sea‐level change produces a defined effect on the geometries of stratigraphic architecture, and that techniques applied for the purpose of delineating sea‐level variation from continental margin strata have a well‐founded conceptual basis.

Schmelz, William J.↗

Uncertainties in greenhouse gas emission factors: A comprehensive analysis of switchgrass‐based biofuel production

Abstract This study investigates uncertainties in greenhouse gas (GHG) emission factors related to switchgrass‐based biofuel production in Michigan. Using three life cycle assessment (LCA) databases—US lifecycle inventory (USLCI) database, GREET, and Ecoinvent—each with multiple versions, we recalculated the global warming intensity (GWI) and GHG mitigation potential in a static calculation. Employing Monte Carlo simulations along with local and global sensitivity analyses, we assess uncertainties and pinpoint key parameters influencing GWI. The convergence of results across our previous study, static calculations, and Monte Carlo simulations enhances the credibility of estimated GWI values. Static calculations, validated by Monte Carlo simulations, offer reasonable central tendencies, providing a robust foundation for policy considerations. However, the wider range observed in Monte Carlo simulations underscores the importance of potential variations and uncertainties in real‐world applications. Sensitivity analyses identify biofuel yield, GHG emissions of electricity, and soil organic carbon (SOC) change as pivotal parameters influencing GWI. Decreasing uncertainties in GWI may be achieved by making greater efforts to acquire more precise data on these parameters. Our study emphasizes the significance of considering diverse GHG factors and databases in GWI assessments and stresses the need for accurate electricity fuel mixes, crucial information for refining GWI assessments and informing strategies for sustainable biofuel production.

Kim, Seungdo↗

A non-orthogonal variational quantum eigensolver

Variational algorithms for strongly correlated chemical and materials systems are one of the most promising applications of near-term quantum computers. We present an extension to the variational quantum eigensolver that approximates the ground state of a system by solving a generalized eigenvalue problem in a subspace spanned by a collection of parametrized quantum states. This allows for the systematic improvement of a logical wavefunction ansatz without a significant increase in circuit complexity. To minimize the circuit complexity of this approach, we propose a strategy for efficiently measuring the Hamiltonian and overlap matrix elements between states parametrized by circuits that commute with the total particle number operator. This strategy doubles the size of the state preparation circuits but not their depth, while adding a small number of additional two-qubit gates relative to standard variational quantum eigensolver. We also propose a classical Monte Carlo scheme to estimate the uncertainty in the ground state energy caused by a finite number of measurements of the matrix elements. We explain how this Monte Carlo procedure can be extended to adaptively schedule the required measurements, reducing the number of circuit executions necessary for a given accuracy. We apply these ideas to two model strongly correlated systems, a square configuration of H4 and the π-system of hexatriene (C6H8).

97 MATHEMATICS AND COMPUTING↗

Impacts of uncertain feedstock quality on the economic feasibility of fast pyrolysis biorefineries with blended feedstocks and decentralized preprocessing sites in the Southeastern United States

Abstract This study performs techno‐economic analysis and Monte Carlo simulations (MCS) to explore the effects that variations in biomass feedstock quality have on the economic feasibility of fast pyrolysis biorefineries using decentralized preprocessing sites (i.e., depots that produce pellets). Two biomass resources in the Southeastern United States, that is, pine residues and switchgrass, were examined as feedstocks. A scenario analysis was conducted for an array of different combinations, including different pellet ash control levels, feedstock blending ratios, different biorefinery capacities, and different biorefinery on‐stream capacities, followed by a comparison with the traditional centralized system. MCS results show that, with depot preprocessing, variations in the feedstock moisture and feedstock ash content can be significantly reduced compared with a traditional centralized system. For a biorefinery operating at 100% of its designed capacity, the minimum fuel selling price (MFSP) of the decentralized system is $3.97–$4.39 per gallon gasoline equivalent (GGE) based on the mean value across all scenarios, whereas the mean MFSP for the traditional centralized system was $3.79–$4.12/GGE. To understand the potential benefits of highly flowable pellets in decreasing biorefinery downtime due to feedstock handling and plugging problems, this study also compares the MFSP of the decentralized system at 90% of its designed capacity with a traditional system at 80%. The analysis illustrates that using low ash pellets mixed with switchgrass and pine residues generates a more competitive MFSP. Specifically, for a biorefinery designed for 2,000 oven dry metric ton per day, running a blended pellet made from 75% switchgrass and 25% pine residues with 2% ash level, and operating at 90% of designed capacity could make an MFSP between $4.49 and $4.71/GGE. In contrast, a traditional centralized biorefinery operating at 80% of designed capacity marks an MFSP between $4.72 and $5.28.

09 BIOMASS FUELS↗

Dynamical Large Deviations of Two-Dimensional Kinetically Constrained Models Using a Neural-Network State Ansatz

We use a neural-network ansatz originally designed for the variational optimization of quantum systems to study dynamical large deviations in classical ones. We use recurrent neural networks to describe the large deviations of the dynamical activity of model glasses, kinetically constrained models in two dimensions. We present the first finite size-scaling analysis of the large-deviation functions of the two-dimensional Fredrickson-Andersen model, and explore the spatial structure of the high-activity sector of the South-or-East model. These results provide a new route to the study of dynamical large-deviation functions, and highlight the broad applicability of the neural-network state ansatz across domains in physics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Shoulder fillet effects in strength distributions of microelectromechanical system components

The failure forces and fracture strengths of polysilicon microelectromechanical system (MEMS) components in the form of stepped tensile bars with shoulder fillets were measured using a sequential failure chain methodology. Approximately 150 specimens for each of four fillet geometries with different stress concentration factors were tested. The resulting failure force and strength distributions of the four geometries were related by a common sidewall flaw population existing within different effective stressed lengths. The failure forces, strengths, and flaw population were well described by a weakest-link based analytical framework. Finite element analysis was used to verify body-force based expressions for the stress concentration factors and to provide insight into the variation of specimen effective length with fillet geometry. Monte Carlo simulations of flaw size and location, based on the strength measurements, were also used to provide insight into fillet shape and size effects. The successful description of the shoulder fillet specimen strengths provides further empirical support for application of the strength and flaw framework in MEMS fabrication and design optimization.

36 MATERIALS SCIENCE↗

Geant4Reweight: a framework for evaluating and propagating hadronic interaction uncertainties in Geant4

Geant4Reweight is an open-source C++ framework that allows users to weight tracks produced by the Geant4 particle transport Monte Carlo simulation according to hadron interaction cross section variations and estimate uncertainties in Geant4 interaction models by comparing the simulation’s hadron interaction cross section predictions to data. The ability to weight hadron transport as simulated by Geant4 is crucial to the propagation of systematic uncertainties related to secondary hadronic interactions in current and upcoming neutrino oscillation experiments, including MicroBooNE, NOvA, and DUNE, aswell as hadron test beam experiments such as ProtoDUNE. Here, we provide motivation for weighting hadron tracks in Geant4 in the context of systematic uncertainty propagation, a description of Geant4’s transport simulation technique, and a description of our weighting technique and fitting framework in the momentum range 0–10 GeV/c, which is typical for the hadrons produced by neutrino interactions in these experiments.

97 MATHEMATICS AND COMPUTING↗

Fixed-phase diffusion Monte Carlo study of activation gap and skyrmion excitations of a ν = 1 system in the presence of charged impurities

The discrepancy between the theoretically calculated and experimentally measured activation gaps in quantum Hall effect has long been a puzzle. We revisit this issue in the context of the v = 1 quantum Hall state, while also incorporating the skyrmion physics. We find that the finite width and the Landau level mixing effects are not sufficient to explain the observed activation gap. We further show that the presence of charged impurities located adjacent to the quantum well can cause a significant reduction in the activation gap, while also causing a suppression of the skyrmion size.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Brazilian CBP - Technoeconomic analysis data

This data is related to the paper entitled "Techno-economic analysis of sugarcane bagasse and straw conversion into cellulosic ethanol via consolidated bioprocessing". That features the evaluation of sugarcane bagasse and straw conversion to ethanol at stand-alone facilities generating electricity from residues. The following scenarios were evaluated: Conventional, featuring hydrothermal pretreatment, fungal cellulase, and yeast fermentation (current commercial standard); Mid-term consolidated bioprocessing (CBP), relying on bagasse solubilization without pretreatment or cotreatment; and Mature CBP, incorporating cotreatment but no pretreatment and considering significant technological advance of the CBP. Available here are the spreadsheets used for Material and Energy balance calculation, Capital and Operational costs estimation and Cash flow analysis. Also available are the description and python code used for Monte Carlo analysis of the ethanol and capital investment variations. This data can be used as a source to implement other techno-economic analysis in the biorefinary context.

09 BIOMASS FUELS↗

Pilgrim : A thermal rate constant calculator and a chemical kinetics simulator

Pilgrim is a program written in Python and designed to use direct dynamics in the calculation of thermal rate constants of chemical reactions by the variational transition state theory (VTST), based on electronic structure calculations for the potential energy surface. Pilgrim can also simulate reaction mechanisms using kinetic Monte Carlo (KMC). For reaction processes with many elementary steps, the rate constant of each of these steps can be calculated by means of conventional transition state theory (TST) or by using VTST. In the current version, Pilgrim can evaluate thermal rates using the canonical version of reaction-path VTST, which requires the calculation of the minimum energy path (MEP) associated with each elementary step or transition structure. Multi-dimensional quantum effects can be incorporated through the small-curvature tunneling (SCT) approximation. These methodologies are available both for reactions involving a single structure of the reactants and the transition state and also for reactions involving flexible molecules with multiple conformations of the reactant and/or of the transition state. For systems with many conformers, the program can evaluate each of the elementary reaction rate constants by multipath canonical VTST or multi-structural VTST. Moreover, the reactant can be unimolecular or bimolecular. Torsional anharmonicity can be incorporated through either the MSTor or the Q2DTor programs. Dual-level calculations are also available in Pilgrim: automatic high-level single-point energies can be used to correct the energy of reactants, transition states, products, and MEP points using the interpolated single-point energies (ISPE) algorithm. When the rate constants of all the chemical processes of interest are known, by means of their calculation using Pilgrim or alternatively through analytical fits to the rate constants as functions of temperature, it is possible to simulate a multistep mechanism under specified laboratory conditions using KMC. Finally, this algorithm allows performing a kinetic simulation to monitor the evolution of each chemical species with time and obtain the product yields.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗