Search NASA⌕ Search

SEARCH · Search NASA

Results for “Sensitivity 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 1,495 records · Page 83
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

SU-E-T-334: Dosimetric Impacts of Organ at Risk (OAR) Contour Errors in Prostate Volumetric Modulated Arc Therapy (VMAT) Planning: A Sensitivity Study

Purpose: To investigate the dosimetric impacts of OAR contour errors in prostate VMAT planning with a sensitivity analysis method. Methods: Ten randomly selected prostate VMAT patients were used to simulate OAR contour error in rectum and bladder. For each OAR 12 even spaced contour control points were selected in every three slices. For each simulation, only one of the control points was simulated to move inwards or outwards by up to 6 voxels and coordinates of adjacent voxels to that moving control point were moved accordingly with a 3D-spline smooth function. An in-house software was used to predict OARs dosimetric endpoints based on geometric relationship between PTV and OAR. In this study, the V75, V70, V65, V60 of rectum and V80, V75, V70, V65 of bladder with and without perturbed contours were calculated and compared. Results: The percentage of OAR dose volume difference around the reference OAR contours were plotted as iso-error lines overlaid on CT images. The significant difference was shown in OAR contours adjacent to the PTV where high dose gradient exists. When one of the 12 points adjacent to the PTV moved up to 6 voxels, maximum errors of 1.99%, 2.46%, 2.89%, and 3.19% were found in V75, V70, V65, and V60 of the rectum, and 1.59%, 1.67%, 1.87%, and 1.90% were found in V80, V75, V70, and V65 of the bladder, respectively. Conclusion: We can quantify the dosimetric impact of OAR contouring error by evaluating the percentage of OAR dose volume changes in prostate VMAT planning. For the adaptive radiation therapy, the iso-error lines can provide the planner a guideline that which portion of contour needed to be checked carefully because dosimetric sensitivity of contour errors are different for different part of contours.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2019

Mechanistic Modeling of Microtopographic Impacts on CO 2 and CH 4 Fluxes in an Alaskan Tundra Ecosystem Using the CLM-Microbe Model

Spatial heterogeneities in soil hydrology have been confirmed as a key control on CO 2 and CH 4 fluxes in the Arctic tundra ecosystem. In this study, we applied a mechanistic ecosystem model, CLM-Microbe, to examine the microtopographic impacts on CO 2 and CH 4 fluxes across seven landscape types in Utqia&#289vik, Alaska: trough, low-centered polygon (LCP) center, LCP transition, LCP rim, high-centered polygon (HCP) center, HCP transition, and HCP rim. We first validated the CLM-Microbe model against static-chamber measured CO 2 and CH 4 fluxes in 2013 for three landscape types: trough, LCP center, and LCP rim. Model application showed that low-elevation and thus wetter landscape types (i.e., trough, transitions, and LCP center) had larger CH 4 emissions rates with greater seasonal variations than high-elevation and drier landscape types (rims and HCP center). Sensitivity analysis indicated that substrate availability for methanogenesis (acetate, CO 2 + H 2 ) is the most important factor determining CH 4 emission, and vegetation physiological properties largely affect the net ecosystem carbon exchange and ecosystem respiration in Arctic tundra ecosystems. Modeled CH 4 emissions for different microtopographic features were upscaled to the eddy covariance (EC) domain with an area-weighted approach before validation against EC-measured CH 4 fluxes. The model underestimated the EC-measured CH 4 flux by 20% and 25% at daily and hourly time steps, suggesting the importance of the time step in reporting CH 4 flux. The strong microtopographic impacts on CO 2 and CH 4 fluxes call for a model-data integration framework for better understanding and predicting carbon flux in the highly heterogeneous Arctic landscape.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2018

Final Technical Report: Mathematical Foundations for Uncertainty Quantification in Materials Design

We developed path-wise information theory-based and goal-oriented sensitivity analysis and parameter identification methods for complex high-dimensional dynamics and in particular of non-equilibrium extended molecular systems. The combination of these novel methodologies provided the first methods in the literature which are capable to handle UQ questions for stochastic complex systems with some or all of the following features: (a) multi-scale stochastic models such as (bio)chemical reaction networks, with a very large number of parameters, (b) spatially distributed systems such as Kinetic Monte Carlo or Langevin Dynamics, (c) non-equilibrium processes typically associated with coupled physico-chemical mechanisms, driven boundary conditions, hybrid micro-macro systems, etc. A particular computational challenge arises in simulations of multi-scale reaction networks and molecular systems. Mathematical techniques were applied to in silico prediction of novel materials with emphasis on the effect of microstructure on model uncertainty quantification (UQ). We outline acceleration methods to make calculations of real chemistry feasible followed by two complementary tasks on structure optimization and microstructure-induced UQ.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

Multi-region Modeling of Daily Global Solar Radiation with Artificial Intelligence Ensemble

Solar radiation data are crucial for the design and evaluation of solar energy systems, climatic studies, water resources management, estimating crop productivity, etc. As so, for locations where direct measurements are not available, reliable models may be developed to estimate solar radiation from more readily available data. In this study, two artificial intelligence (AI)-based models including artificial neural network and adaptive neuro-fuzzy inference systems, three temperature-based empirical models including Meza–Varas, Hargreaves–Samani, and Chen, and a conventional multi-linear regression (MLR) model were employed for multi-region daily global solar radiation estimation for Iraq. To ensure appropriate selection of input variables, sensitivity analysis was conducted to determine the dominant parameters. Finally, two ensemble approaches, neural average ensemble and simple average ensemble, were applied to improve the performance of the single models. For this purpose, daily meteorological data of maximum temperature (T{sub max}), minimum temperature (T{sub min}), mean temperature (T{sub mean}), relative humidity (R{sub H}), and wind speed (U{sub 2}) were obtained from January 2006 to December 2016 from four major cities in Iraq representing, north, west, south, and east regions. The results revealed that temperatures (T{sub max}, T{sub mean}, T{sub min}) and relative humidity are the dominant parameters. While temperature-based empirical models and MLR model could be employed to achieve reliable results, AI-based models are superior in performance to other models. Also promising improvement in daily global solar radiation modeling could be achieved by model ensemble. The results of this study affirmed that the provided ensemble approaches can increase the performance of single models up to 19.19%, 7.59%, and 16.81% in training, validation, and testing steps, respectively.

58 GEOSCIENCES↗
DOE OSTI2019

Optimization of cool roof and night ventilation in office buildings: A case study in Xiamen, China

Increasing roof albedo (using a “cool” roof) and night ventilation are passive cooling technologies that can reduce the cooling loads in buildings, but existing studies have not comprehensively explored the potential benefits of integrating these two technologies. This study combines an experiment in the summer and transition seasons with an annual simulation so as to evaluate the thermal performance, energy savings and thermal comfort improvement that could be obtained by coupling a cool roof with night ventilation. A holistic approach integrating sensitivity analysis and multi-objective optimization is developed to explore key design parameters (roof albedo, night ventilation air change rate, roof insulation level and internal thermal mass level) and optimal design options for the combined application of the cool roof and night ventilation. The proposed approach is validated and demonstrated through studies on a six-storey office building in Xiamen, a cooling-dominated city in southeast China. Simulations show that combining a cool roof with night ventilation can significantly decrease the annual cooling energy consumption by 27% compared to using a black roof without night ventilation and by 13% compared to using a cool roof without night ventilation. Roof albedo is the most influential parameter for both building energy performance and indoor thermal comfort. Optimal use of the cool roof and night ventilation can reduce the annual cooling energy use by 28% during occupied hours when air-conditioners are on and reduce the uncomfortable time slightly during occupied hours when air-conditioners are off.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2015

Interpretation of thermoreflectance measurements with a two-temperature model including non-surface heat deposition

We develop a solution to the two-temperature diffusion equation in axisymmetric cylindrical coordinates to model heat transport in thermoreflectance experiments. Our solution builds upon prior solutions that account for two-channel diffusion in each layer of an N-layered geometry, but adds the ability to deposit heat at any location within each layer. We use this solution to account for non-surface heating in the transducer layer of thermoreflectance experiments that challenge the timescales of electron-phonon coupling. A sensitivity analysis is performed to identify important parameters in the solution and to establish a guideline for when to use the two-temperature model to interpret thermoreflectance data. We then fit broadband frequency domain thermoreflectance (BB-FDTR) measurements of SiO{sub 2} and platinum at a temperature of 300 K with our two-temperature solution to parameterize the gold/chromium transducer layer. We then refit BB-FDTR measurements of silicon and find that accounting for non-equilibrium between electrons and phonons in the gold layer does lessen the previously observed heating frequency dependence reported in Regner et al. [Nat. Commun. 4, 1640 (2013)] but does not completely eliminate it. We perform BB-FDTR experiments on silicon with an aluminum transducer and find limited heating frequency dependence, in agreement with time domain thermoreflectance results. We hypothesize that the discrepancy between thermoreflectance measurements with different transducers results in part from spectrally dependent phonon transmission at the transducer/silicon interface.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Nested polynomial trends for the improvement of Gaussian process-based predictors

The role of simulation keeps increasing for the sensitivity analysis and the uncertainty quantification of complex systems. Such numerical procedures are generally based on the processing of a huge amount of code evaluations. When the computational cost associated with one particular evaluation of the code is high, such direct approaches based on the computer code only, are not affordable. Surrogate models have therefore to be introduced to interpolate the information given by a fixed set of code evaluations to the whole input space. When confronted to deterministic mappings, the Gaussian process regression (GPR), or kriging, presents a good compromise between complexity, efficiency and error control. Such a method considers the quantity of interest of the system as a particular realization of a Gaussian stochastic process, whose mean and covariance functions have to be identified from the available code evaluations. In this context, this work proposes an innovative parametrization of this mean function, which is based on the composition of two polynomials. This approach is particularly relevant for the approximation of strongly non linear quantities of interest from very little information. After presenting the theoretical basis of this method, this work compares its efficiency to alternative approaches on a series of examples.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2013

Adaptive polynomial chaos techniques for uncertainty quantification of a gas cooled fast reactor transient

Uncertainty quantification plays an increasingly important role in the nuclear community, especially with the rise of Best Estimate Plus Uncertainty methodologies. Sensitivity analysis, surrogate models, Monte Carlo sampling and several other techniques can be used to propagate input uncertainties. In recent years however polynomial chaos expansion has become a popular alternative providing high accuracy at affordable computational cost. This paper presents such polynomial chaos (PC) methods using adaptive sparse grids and adaptive basis set construction, together with an application to a Gas Cooled Fast Reactor transient. Comparison is made between a new sparse grid algorithm and the traditionally used technique proposed by Gerstner. An adaptive basis construction method is also introduced and is proved to be advantageous both from an accuracy and a computational point of view. As a demonstration the uncertainty quantification of a 50% loss of flow transient in the GFR2400 Gas Cooled Fast Reactor design was performed using the CATHARE code system. The results are compared to direct Monte Carlo sampling and show the superior convergence and high accuracy of the polynomial chaos expansion. Since PC techniques are easy to implement, they can offer an attractive alternative to traditional techniques for the uncertainty quantification of large scale problems. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2016

Inverse Uncertainty Quantification of Reactor Simulation with Polynomial Chaos Surrogate Model

Reliable uncertainty quantification (UQ) and sensitivity analysis (SA) is significant part of modelling and simulation. This is particularly important in the nuclear community due to the fact that the licensing of nuclear installations has shifted from conservative to Best Estimate plus Uncertainty (BEPU) methodologies. UQ can be broadly classified as forward UQ and inverse UQ. Forward UQ is the process of quantifying the uncertainties in simulation outputs by propagating the uncertainties in input parameters. It requires knowledge in the input uncertainty, which is often specified by ad-hoc expert judgment which lacks mathematical rigor. The inverse UQ is the process to determine the uncertainty in input parameters given observations of the output. The purpose of this paper is to replace such ad-hoc expert judgment of the statistical properties of input model parameters in nuclear reactor simulation. This summary presents work on inverse UQ of a nonlinear Point Reactor Kinetics Equations - lumped thermal-hydraulics parameter coupled problem. Bayesian inference and Markov Chain Monte Carlo (MCMC) are used for the inverse UQ. Because MCMC requires a large number of simulations which is impractical for models that are computationally expensive, surrogate models is constructed by generalized Polynomial Chaos (PC) expansion technique. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS1992

Integrating aerodynamics and structures in the minimum weight design of a supersonic transport wing

An approach is presented for determining the minimum weight design of aircraft wing models which takes into consideration aerodynamics-structure coupling when calculating both zeroth order information needed for analysis and first order information needed for optimization. When performing sensitivity analysis, coupling is accounted for by using a generalized sensitivity formulation. The results presented show that the aeroelastic effects are calculated properly and noticeably reduce constraint approximation errors. However, for the particular example selected, the error introduced by ignoring aeroelastic effects are not sufficient to significantly affect the convergence of the optimization process. Trade studies are reported that consider different structural materials, internal spar layouts, and panel buckling lengths. For the formulation, model and materials used in this study, an advanced aluminum material produced the lightest design while satisfying the problem constraints. Also, shorter panel buckling lengths resulted in lower weights by permitting smaller panel thicknesses and generally, by unloading the wing skins and loading the spar caps. Finally, straight spars required slightly lower wing weights than angled spars.

Barthelemy, Jean-Francois M.↗
DOE OSTI2015

A feasibility study: Selection of a personalized radiotherapy fractionation schedule using spatiotemporal optimization

Purpose: To investigate the impact of using spatiotemporal optimization, i.e., intensity-modulated spatial optimization followed by fractionation schedule optimization, to select the patient-specific fractionation schedule that maximizes the tumor biologically equivalent dose (BED) under dose constraints for multiple organs-at-risk (OARs). Methods: Spatiotemporal optimization was applied to a variety of lung tumors in a phantom geometry using a range of tumor sizes and locations. The optimal fractionation schedule for a patient using the linear-quadratic cell survival model depends on the tumor and OAR sensitivity to fraction size (α/β), the effective tumor doubling time (T{sub d}), and the size and location of tumor target relative to one or more OARs (dose distribution). The authors used a spatiotemporal optimization method to identify the optimal number of fractions N that maximizes the 3D tumor BED distribution for 16 lung phantom cases. The selection of the optimal fractionation schedule used equivalent (30-fraction) OAR constraints for the heart (D{sub mean} ≤ 45 Gy), lungs (D{sub mean} ≤ 20 Gy), cord (D{sub max} ≤ 45 Gy), esophagus (D{sub max} ≤ 63 Gy), and unspecified tissues (D{sub 05} ≤ 60 Gy). To assess plan quality, the authors compared the minimum, mean, maximum, and D{sub 95} of tumor BED, as well as the equivalent uniform dose (EUD) for optimized plans to conventional intensity-modulated radiation therapy plans prescribing 60 Gy in 30 fractions. A sensitivity analysis was performed to assess the effects of T{sub d} (3–100 days), tumor lag-time (T{sub k} = 0–10 days), and the size of tumors on optimal fractionation schedule. Results: Using an α/β ratio of 10 Gy, the average values of tumor max, min, mean BED, and D{sub 95} were up to 19%, 21%, 20%, and 19% larger than those from conventional prescription, depending on T{sub d} and T{sub k} used. Tumor EUD was up to 17% larger than the conventional prescription. For fast proliferating tumors with T{sub d} less than 10 days, there was no significant increase in tumor BED but the treatment course could be shortened without a loss in tumor BED. The improvement in the tumor mean BED was more pronounced with smaller tumors (p-value = 0.08). Conclusions: Spatiotemporal optimization of patient plans has the potential to significantly improve local tumor control (larger BED/EUD) of patients with a favorable geometry, such as smaller tumors with larger distances between the tumor target and nearby OAR. In patients with a less favorable geometry and for fast growing tumors, plans optimized using spatiotemporal optimization and conventional (spatial-only) optimization are equivalent (negligible differences in tumor BED/EUD). However, spatiotemporal optimization yields shorter treatment courses than conventional spatial-only optimization. Personalized, spatiotemporal optimization of treatment schedules can increase patient convenience and help with the efficient allocation of clinical resources. Spatiotemporal optimization can also help identify a subset of patients that might benefit from nonconventional (large dose per fraction) treatments that are ineligible for the current practice of stereotactic body radiation therapy.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

Subcritical Copper-Reflected-α-phase Plutonium (SCRαP) Integral Experiment Design

A Subcritical Copper-Reflected α-phase Plutonium (SCRαP) integral benchmark experiment is proposed where multiplication will be approximated using correlated neutron data from a detector system consisting of {sup 3}He tubes inside polyethylene. Measurements will be performed on various subcritical experimental configurations consisting of a weapons-grade plutonium sphere surrounded by different Cu thicknesses. In addition to the proposed base experimental configurations with Cu, additional configurations are planned with the plutonium ball nested in various thicknesses of interleaved polyethylene spherical shells mixed in with the Cu shells. The polyethylene is intended to provide fast neutron moderation and reflection, resulting in additional measurements with differing multiplication, spectra, and nuclear data sensitivity. The proposed experiments will performed at the National Criticality Experiments Research Center (NCERC). NCERC is home to a 4.5-kg α-phase stainless steel clad plutonium sphere that is referred to as the BeRP (Beryllium-Reflected Plutonium) ball due to its historical use in a beryllium-reflected critical experiment. In 2014, similar subcritical measurements were performed with the BeRP ball surrounded by nickel and tungsten. The nickel-reflected measurements were documented as a benchmark evaluation which was accepted by the International Criticality Safety Benchmark Evaluation Project (ICSBEP). Similar to past measurements, the proposed work will help both identify deficiencies and quantify uncertainties in nuclear data and validate computational methods related to neutron multiplication inference for subcritical benchmark evaluations. Monte Carlo simulations were performed for the preliminary design of a subcritical Copper-Reflected - phase Plutonium (SCRP) integral benchmark experiment. Specifically, MCNP{sup R} version 6 was used to determine the k{sub eff} and associated nuclear data sensitivities for different preliminary configurations (reflector thickness, material type). The final experimental design process will include continuous-energy nuclear data sensitivity analysis. This will help assess the potential impact of this integral measurement as it relates to international efforts to continuously improve current libraries (ENDF, JEFF, JENDL etc) that still tend to over/under-estimate the measured results of integral experiments, sometimes significantly. Deficiencies in underlying nuclear data quantities such as Nubar have been shown to have an effect on inferred values from subcritical measurements as well. The proposed experiments and subsequent computational validation will help identify such deficiencies. It will also help validate new nuclear data evaluations, including one for Cu that was recently performed in the resolved region up to 100 MeV for {sup 65}Cu and {sup 63}Cu which is expected to improve Cu-related benchmark simulation performance. The final experimental design phase will also include simulation results that will incorporate the detectors that will be used for the subcritical measurements. For the benchmark portion of the measurements, the main detectors that are being considered are portable neutron multiplicity counters made up of {sup 3}He tubes embedded in polyethylene and lined with cadmium. These detectors employ a data acquisition system that produces list-mode data (correlated information where each event is listed sequentially) that can be analyzed using the Feynman variance-to-mean method. Other detector systems including SNAP (single {sup 3}He tube embedded in polyethylene), and a high-purity germanium (HPGe) detector are being considered for supplementary measurements only and are not intended to be used for the final benchmark measurements. Decoupling the supplementary detectors from the system eliminates any contribution to the total system multiplication, and greatly simplifies the final computational benchmark model. Therefore, an optimization of the portable neutron multiplicity counter placement in the benchmark measurements will be performed to ensure that there is minimal contribution to the total system multiplication. In addition, new temperature measurement techniques and hardware will be explored. (authors)

42 ENGINEERING↗
DOE OSTI2016

Emulation of dynamic simulators with application to hydrology

Many simulation-intensive tasks in the applied sciences, such as sensitivity analysis, parameter inference or real time control, are hampered by slow simulators. Emulators provide the opportunity of speeding up simulations at the cost of introducing some inaccuracy. An emulator is a fast approximation to a simulator that interpolates between design input–output pairs of the simulator. Increasing the number of design data sets is a computationally demanding way of improving the accuracy of emulation. We investigate the complementary approach of increasing emulation accuracy by including knowledge about the mechanisms of the simulator into the formulation of the emulator. To approximately reproduce the output of dynamic simulators, we consider emulators that are based on a system of linear, ordinary or partial stochastic differential equations with a noise term formulated as a Gaussian process of the parameters to be emulated. This stochastic model is then conditioned to the design data so that it mimics the behavior of the nonlinear simulator as a function of the parameters. The drift terms of the linear model are designed to provide a simplified description of the simulator as a function of its key parameters so that the required corrections by the conditioned Gaussian process noise are as small as possible. The goal of this paper is to compare the gain in accuracy of these emulators by enlarging the design data set and by varying the degree of simplification of the linear model. We apply this framework to a simulator for the shallow water equations in a channel and compare emulation accuracy for emulators based on different spatial discretization levels of the channel and for a standard non-mechanistic emulator. Our results indicate that we have a large gain in accuracy already when using the simplest mechanistic description by a single linear reservoir to formulate the drift term of the linear model. Adding some more reservoirs does not lead to a significant improvement in accuracy. However, the transition to a spatially continuous linear model leads again to a similarly large gain in accuracy as the transition from the non-mechanistic emulator to that based on one reservoir.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Radiation Transport in Random Media With Large Fluctuations

Strongly heterogeneous media arise in several applications that include radiation shields, nuclear fuel, BWR moderators, clouds, planetary and stellar atmospheres, turbulent gases and plasmas. Neutral and charged particle transport computations in such media have relied heavily on formulating transport equations with spatially random coefficients (physical data) and developing solution methods to deal with the additional stochastic dimensions. Attempts at developing approximate closures that yield only low order statistical information (e.g., mean and variance of the flux) have proved to be highly restrictive under real physics conditions or rely on techniques that require fluctuation amplitudes to be small for robustness. Recently, stochastic spectral methods such as polynomial chaos and stochastic collocation have been developed for aleatoric uncertainty quantification and sensitivity analysis, and successfully applied in radiation transport work. Advances in UQ techniques have tended to focus on efficiently handling large numbers of uncertain variables but the rigorous stochastic basis of the approach also promotes its use in situations where stochasticity is due to spatial heterogeneity and the associated uncertainty is large. Here we apply these techniques to radiation transport in media with spatially randomly varying cross sections without restriction on fluctuation amplitudes. Specifically, we represent the cross sections as a lognormal spatial random process with specified mean, variance and covariance function and use a Karhunen-Loeve (KL) decomposition to generate cross section realizations that are strictly positive. Woodcock Monte Carlo (WMC) is then used to simulate transport using random sampling of cross sections and deterministic sampling based on a stochastic collocation technique. Numerical results for the mean and variance in the scalar flux and leakage currents are obtained for weak and strong random variations. We focus in this work on demonstrating the approach and defer application to specific problems to a future investigation. We have demonstrated that a combination of stochastic spectral representation and Woodcock Monte Carlo simulation allows efficient computation of radiation transport in strongly random media. The use of a lognormal distribution to represent the fluctuations in medium properties facilitates enforcement of strict positivity of cross section realizations, a condition that has challenged standard approaches that assume Gaussian fluctuations. The fact that the lognormal random process is a memoryless nonlinear transformation of a Gaussian process further enables efficient and accurate reconstruction of the random cross section realizations using a Karhunen-Loeve representation of the cross section. Finally, the robustness of this model enables fluctuations of arbitrarily large amplitude to be studied, using both random sampling Monte Carlo and more efficient stochastic collocation techniques. The present work is being extended to allow reconstruction of probability density functions of output variables using polynomial chaos representation, and to accommodate general cross section covariance functions. Also, techniques are being explored to more judiciously select ceiling cross sections to enhance the efficiency of Woodcock Monte Carlo for transport in random media. (authors)

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2016

Nuclear Data Uncertainty Quantification for Adjoint-Weighted Quantities via XGPT in Monte Carlo Codes

The Total Monte Carlo method (TMC) is an innovative approach to nuclear data uncertainty propagation. In TMC, a large number of independent ENDF files are randomly generated starting from nuclear models parameters and their uncertainties. These files are then processed with NJOY or similar codes to produce a set of formatted continuous-energy cross sections (i.e., ACE files), and independent Monte Carlo neutron transport simulations are run with these different ACE files as input. Finally the distributions of the considered response functions are obtained directly from the results of the Monte Carlo runs. Despite several improvements of the TMC method, this approach is often considered an inefficient way of propagating nuclear data uncertainties, as it requires several computationally expensive Monte Carlo runs. Covariance-based methods present important advantages for several applications, especially when different sources of uncertainties are to be investigated. For this reason, sensitivity analysis and uncertainty quantification capabilities are being implemented in several Monte Carlo neutron transport codes. On the other hand, the fact that no artificial multi-group covariance matrices are required and that higher moments of the responses distributions can be easily obtained, are definite advantages of the TMC method. In an effort to overcome the limitations of both approaches, a new method, called eXtended Generalized Perturbation Theory (XGPT) is under development. When performing nuclear data uncertainty propagation, such XGPT approach does not require generating discretized (i.e., multigroup) covariance matrices, and allows to estimate higher moments of the response probability distributions, via linear perturbation theory. One of the promising applications of the XGPT method involves its adoption in conjunction with the Proper Orthogonal Decomposition (POD) of nuclear data, to build a Reduced Order Model (ROM) for the approximation of thousands of independent TMC results in a single criticality source simulation. In the present work, this approach is employed for the first time to evaluate the impact of uncertainty in nuclear data on adjoint-weighted quantities. The method has been implemented in an extended version of the Serpent-2.1.24 Monte Carlo code and is demonstrated on two fast spectrum benchmarks, focusing on {sup 208}Pb nuclear data. A new method for propagating nuclear data uncertainties in criticality source Monte Carlo simulations has been developed. This eXtended Generalized Perturbation Theory (XGPT) approach allows to accurately reproduce reference Total Monte Carlo (TMC) uncertainty estimates in a single criticality source simulation. The present work aims to demonstrate that the proposed method can quantify nuclear data uncertainty based on adjoint-weighted quantities. The XGPT approach was applied to evaluate the impact of {sup 208}Pb nuclear data uncertainties on the effective prompt neutron lifetime (l{sub eff}) in two fast spectrum systems. The obtained XGPT results show a good agreement with the reference TMC results, obtained from independent Monte Carlo runs with {sup 208}Pb nuclear data from random TENDL-2013 evaluations. The XGPT method, implemented in a extended version of the Serpent Monte Carlo code and adopted along with the Proper Orthogonal Decomposition (POD) of the nuclear data, is able to provide an accurate reduced order model to approximate thousands TMC results within a single run, ensuring a large saving in computational requirements. The effect of {sup 208}Pb nuclear data uncertainties appears to be rather large. This is consistent with previous studies that adopted the same input random evaluations data from TENDL. (authors)

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Effects of radiative heat loss on extinction limits of counterflow premixed ammonia-air flames

The effect of radiative heat loss on counterflow premixed ammonia-air flames has been investigated as a function of equivalence ratio and system pressure. Non-adiabatic counterflow premixed flame simulations incorporating detailed chemistry, transport, and radiative heat transfer in the optically thin limit have been employed to elucidate important underlying physics and extinction characteristics. Similar to methane-air flame systems, non-adiabatic counterflow premixed computations show that the combined effects of positive flame stretch, sub-unity Lewis number, and radiative heat loss lead to the extension of the lean flammability limit for ammonia-air flames, revealing a C-shaped curve near the lean flammability limit that exhibits both radiation-induced and stretch-induced extinction states. Here, the computations compare favorably with the experimentally determined extinction stretch rate values for the stretch-induced extinction states over a wide range of equivalence ratios at different pressures up to 5 atm. A novel feature sensitivity analysis has also been developed to highlight important sensitive reactions for the stretch- and radiation-induced extinction states. Furthermore, the controlling chemistry at the dual-extinction states and at varying pressures are compared and discussed.

10 SYNTHETIC FUELS↗
DOE OSTI2016

Modeling of GaAs/Silicon PVMirror Tandem System: A Case Study

We propose a new, optically coupled tandem architecture with a silicon “PVMirror” as the bottom cell. A silicon PVMirror is a curved silicon module with an integrated dichroic mirror that transmits near-infrared light to the silicon cells and reflects and focuses visible light onto a GaAs receiver. Simulating such a system with record GaAs and silicon cells and an ideal dichroic mirror shows the tandem efficiency can reach over 36% with no diffuse light and over 30% with 20% diffuse light. A sensitivity analysis indicates that the PVMirror tandem is most sensitive to the GaAs cell performance and dichroic mirror cut-off wavelength but relatively insensitive to the silicon cell performance.

Yu, Zhengshan (Jason)↗
DOE OSTI2022

Fluid – Fluid Interfacial Area and Its Impact on Relative Permeability - A Pore Network Modeling Study

Abstract Relative permeability (kr) is commonly modeled as an empirical function of phase saturation. Although current empirical models can provide a good match of one or two measured relative permeabilities using saturation alone, they are unable to predict relative permeabilities well when there is hysteresis or when physical properties such as wettability change. Further, current models often result in relative permeability discontinuities that can cause convergence and accuracy problems in simulation. To overcome these problems, recent research has modeled relative permeability as a state function of both saturation (S) and phase connectivity (X). Pore network modeling (PNM) data, however, shows small differences in relative permeability for the same S-X value when approached from a different flow direction. This paper examines the impact of one additional Minkowski parameter (Mecke and Arns, 2005), the fluid-fluid interfacial area, on relative permeability to identify if that satisfactorily explains this discrepancy. We calculate the total fluid-fluid interfacial areas (IA) during two-phase (oil/water) flow in porous media using pore network modeling. The area is calculated from PNM simulations using the areas associated with corners and throats in pore elements of different shapes. The pore network is modeled after a Bentheimer sandstone, using square, triangular, and circular pore shapes. Simulations were conducted for numerous primary drainage and imbibition cycles at a constant contact angle of 0° for the wetting phase. Simultaneous measurements of capillary pressure, relative permeability, saturation, and phase connectivity are made for each displacement. Fluid-fluid interfacial area is calculated from the PNM capillary pressure, the fluid location in the pore elements, and the pore element dimensional data. The results show that differences in the relative permeability at the same (S,X) point is explained well by differences in the fluid-fluid interfacial area (IA). That is, for a larger change in IA at these intersection points, the permeability difference is greater. That difference in relative permeability approaches zero as the difference in IA approaches zero. This confirms that relative permeability can be modeled better as a unique function of S, X, and IA. The results also show that an increase in IA restricts flow decreasing the nonwetting (oil) phase permeability. This decrease is caused by an increase in the throat area fraction compared to the corner area as the total area IA increases. The wetting phase relative permeability, however, shows the inverse trend, in that its relative permeability is greater when IA becomes larger owing to a greater fraction of the total area associated with the corners. The area IA, however, impacts the nonwetting phase relative permeability more than the wetting phase relative permeability. Corner flow improves the wetting phase relative permeability because the wetting phase is continuous there. Finally, a sensitivity analysis shows that relative permeability a is more sensitive to change in S than they are for IA for the case studied implying that if only two parameters are used to model relative permeability it is better to choose S and X.

Engineering↗