SEARCH · Engineering Papers

Results for “Global 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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 163 records · Page 9

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Improving North American Wildfire Prediction by Integrating a Machine-Learning Fire Model in a Land Surface Model

Wildfires have shown increasing trends in both frequency and severity across the Contiguous United States (CONUS). However, process-based fire models have difficulties in accurately simulating the burned area over the CONUS due to a simplification of the physical process and cannot capture the interplay among fire, ignition, climate, and human activities. The deficiency of burned area simulation deteriorates the description of fire impact on energy balance, water budget, and carbon fluxes in the Earth System Models (ESMs). Alternatively, machine learning (ML) based fire models, which capture statistical relationships between the burned area and environmental factors, have shown promising burned area predictions and corresponding fire impact simulation. We develop a hybrid framework (ML4Fire-XGB) that integrates a pretrained eXtreme Gradient Boosting (XGBoost) wildfire model with the Energy Exascale Earth System Model (E3SM) land model (ELM) version 2.1. A Fortran-C-Python deep learning bridge is adapted to support online communication between ELM and the ML fire model. Specifically, the burned area predicted by the ML-based wildfire model is directly passed to ELM to adjust the carbon pool and vegetation dynamics after disturbance, which are then used as predictors in the ML-based fire model in the next time step. Evaluated against the historical burned area from Global Fire Emissions Database 5 from 2001-2020, the ML4Fire-XGB model outperforms process-based fire models in terms of spatial distribution and seasonal variations. Sensitivity analysis confirms that the ML4Fire-XGB well captures the responses of the burned area to rising temperatures. The ML4Fire-XGB model has proved to be a new tool for studying vegetation-fire interactions, and more importantly, enables seamless exploration of climate-fire feedback, working as an active component in E3SM.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Simulated wildfire burned area over the CONUS during 2001-2020

Wildfires have shown increasing trends in both frequency and severity across the Contiguous United States (CONUS). However, process-based fire models have difficulties in accurately simulating the burned area over the CONUS due to a simplification of the physical process and cannot capture the interplay among fire, ignition, climate, and human activities. The deficiency of burned area simulation deteriorates the description of fire impact on energy balance, water budget, and carbon fluxes in the Earth System Models (ESMs). Alternatively, machine learning (ML) based fire models, which capture statistical relationships between the burned area and environmental factors, have shown promising burned area predictions and corresponding fire impact simulation. We develop a hybrid framework (ML4Fire-XGB) that integrates a pretrained eXtreme Gradient Boosting (XGBoost) wildfire model with the Energy Exascale Earth System Model (E3SM) land model (ELM). A Fortran-C-Python deep learning bridge is adapted to support online communication between ELM and the ML fire model. Specifically, the burned area predicted by the ML-based wildfire model is directly passed to ELM to adjust the carbon pool and vegetation dynamics after disturbance, which are then used as predictors in the ML-based fire model in the next time step. Evaluated against the historical burned area from Global Fire Emissions Database 5 from 2001-2020, the ML4Fire-XGB model outperforms process-based fire models in terms of spatial distribution and seasonal variations. Sensitivity analysis confirms that the ML4Fire-XGB well captures the responses of the burned area to rising temperatures. The ML4Fire-XGB model has proved to be a new tool for studying vegetation-fire interactions, and more importantly, enables seamless exploration of climate-fire feedback, working as an active component in E3SM.

Liu, Ye↗

DOE OSTI2024

Ignition Delay Times and Chemical Kinetic Model Validation for Hydrogen and Ammonia Blending With Natural Gas at Gas Turbine Relevant Conditions

Ignition delay times from undiluted mixtures of natural gas (NG)/H 2 /Air and NG/NH 3 /Air were measured using a high-pressure shock tube at the University of Central Florida. The combustion temperatures were experimentally tested between 1000 and 1500 K near a constant pressure of 25 bar. As mentioned, mixtures were kept undiluted to replicate the same chemistry pathways seen in gas turbine combustion chambers. Recorded combustion pressures exceeded 200 bar due to the large energy release, hence why these were performed at the high-pressure shock tube facility. The data are compared to the predictions of the NUIGMech 1.1 mechanism for chemical kinetic model validation and refinement. An exceptional agreement was shown for stoichiometric conditions in all cases but strayed at lean and rich equivalence ratios, especially in the lower temperature regime of H 2 addition and all temperature ranges of the baseline NG mixture. Hydrogen addition also decreased ignition delay times by nearly 90%, while NH 3 fuel addition made no noticeable difference in ignition time. NG/NH 3 exhibited similar chemistry to pure NG under the same conditions, which is shown in a sensitivity analysis. Here, the reaction CH 3 + O 2 = CH 3 O + O is identified and suggested as a possible modification target to improve model performance. Increasing the robustness of chemical kinetic models via experimental validation will directly aid in designing next-generation combustion chambers for use in gas turbines, which in turn will greatly lower global emissions and reduce greenhouse effects.

33 ADVANCED PROPULSION SYSTEMS↗

DOE OSTI2025

Neural Active Manifolds: Nonlinear Dimensionality Reduction for Uncertainty Quantification

We present a new approach for nonlinear dimensionality reduction, specifically designed for computationally expensive mathematical models. We leverage autoencoders to discover a one-dimensional neural active manifold (NeurAM) capturing the model output variability, through the aid of a simultaneously learnt surrogate model with inputs on this manifold. Our method only relies on model evaluations and does not require the knowledge of gradients. The proposed dimensionality reduction framework can then be applied to assist outer loop many-query tasks in scientific computing, like sensitivity analysis and multifidelity uncertainty propagation. In particular, we prove, both theoretically under idealized conditions, and numerically in challenging test cases, how NeurAM can be used to obtain multifidelity sampling estimators with reduced variance by sampling the models on the discovered low-dimensional and shared manifold among models. Several numerical examples illustrate the main features of the proposed dimensionality reduction strategy and highlight its advantages with respect to existing approaches in the literature.

Autoencoders↗

DOE OSTI2022

Light-Duty Vehicle Trip Classification Using One-Class Novelty Detection and Exhaustive Feature Extraction

Travel mode classification within travel survey data sets, especially light-duty vehicle (LDV) trips, is foundational, though nontrivial, to emerging mobility systems, travel behavior analysis, and fuel consumption estimation. Current travel mode detection approaches require well-sampled and balanced data sets with ground truth travel mode labels. The detection approaches are rarely applied and validated on large-scale, real-world data sets, which may not satisfy the dataset requirements. This work proposes an LDV trip detection model as a supplement to current travel mode detection methods, for the case when the training set is highly (and/or completely) unbalanced, to the extent that classical machine-learning approaches become difficult or impossible to deploy. The proposed model uses a novelty detection technique - one-class support vector machines (OCSVMs) - and a novel exhaustive feature extraction (EFE) technique on continuous time series data (i.e., Global Positioning System [GPS] speed profiles) for single-mode trip trajectories. Training and validation of the model are conducted on a large-scale, real-world data set. The proposed method accurately identifies LDV trips from a broad set of multimodal trips by leveraging a wealth of preexisting in-vehicle GPS travel data. Additional sensitivity analysis sheds light on the optimal training size and feature selection, which will benefit applications limited by highly imbalanced data. The paper also discusses performance comparison with regular machine-learning approaches, the model's robustness, and the potential to extend the proposed model to multi-modal trip prediction.

33 ADVANCED PROPULSION SYSTEMS↗

DOE OSTI2020

Sensitivity Study for Forecasting Variables of WRF-Solar Using a Tangent Linear Approach

Integrating solar generation in recent years has highlighted the need for improved accuracy in predicting solar power. Confidence in solar power forecasting can be achieved by designing an ensemble that provides reliable probabilistic information for solar radiation with reduced uncertainty and error. Ideally, ensemble members are created through the optimized perturbation of the initial conditions in numerical weather prediction (NWP) models. Tangent linear models are capable of efficiently investigating the sensitivity of solar radiation to model input parameters because they do not require individual perturbation of each variable. This sensitivity study using tangent linear models provide us the capability to identify the right variables to perturb in an ensemble prediction system. In this study, we developed tangent linear models for WRF-Solar modules that directly impact the computation of solar radiation and the simulation of cloud formation and dissipation including the Fast All-sky Model for Solar Applications (FARMS), the Noah land surface model (LSM), the Thompson microphysics, the Mello-Yamada-Nakanishi-Niino (MYNN) boundary layer parameterization, and the Deng scheme for a shallow-convection parameterization. A sensitivity analysis was conducted under various scenarios based on satellite observations and model simulations from the National Solar Radiation Data Base (NSRDB) and WRF-Solar, respectively. Critical forecasting variables that are highly sensitive to the forecasting of global horizontal irradiance (GHI), direct normal irradiance (DNI), cloud mixing ratio, cloud tendency, cloud fraction, and sensible and latent heat fluxes were determined using the relevant WRF-Solar module. This study will be used as a guidance on future research leading to high-quality probabilistic solar forecasting. In this presentation, we discuss the validation of tangent linear approach for WRF-Solar modules and illustrate how the sensitivity results are valuable in the improvement of probabilistic solar prediction.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2021

Global Microphysical Sensitivity of Superparameterized Precipitation Extremes

Abstract A recent study found statistically significant differences in extreme precipitation distributions over the contiguous United States (CONUS) when changing the microphysics scheme in a superparameterized global climate model. Here, we repeat the analysis globally and similarly find that differences are widespread when varying the number of predicted moments in the microphysics parameterization, but not when comparing variants of the double‐moment scheme. However, contrary to the previous study in which differences largely disappeared over CONUS when 5‐day simulations were conducted, we found that the signal in these shorter integrations remains within the tropics, implying a direct local effect of microphysics on precipitation extremes in these regions. The effect on precipitation is traced back to changes in vertical velocity profiles changes that are then amplified in the climatological simulations compared to the 5‐day ones. Finally, the superparameterized extremes, regardless of the microphysics scheme, are shown to be larger than those from the Global Precipitation Climatology Project One‐Degree Daily data set and generally smaller than those from the Tropical Rainfall Measuring Mission 3B42 data set.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

How Well do Earth System Models Capture Apparent Relationships Between Phytoplankton Biomass and Environmental Variables?

Abstract As phytoplankton form the base of the marine food web, understanding the controls on their abundance is fundamental to understanding marine ecology and its sensitivity to global climate change. While many Earth System Models (ESMs) predict phytoplankton biomass, it is unclear whether they properly capture the mechanistic relationships that control this quantity in the real ocean. We used Random Forest analysis to analyze the output of 13 ESMs as well as two observational data sets. The target variable was phytoplankton carbon and the predictors included environmental parameters known to influence phytoplankton, including nutrients, light, mixed layer depth, salinity, temperature, and upwelling. We examined the following: (a) What fractions of variability in ESMs and observations can be linked to the large‐scale environmental variables simulated by ESMs? (b) What are the dominant predictors and relationships affecting phytoplankton biomass? (c) How well do ESMs simulate phytoplankton carbon and do they simulate the relationships we see in observations? About 88%–96% of the variability in observational data sets and greater than 98% in the ESMs was accounted for by environmental variables known to influence phytoplankton biomass. The dominant predictors in the observational data sets were shortwave radiation and dissolved iron, with temperature and ammonium also relatively important. All the ESMs show that shortwave radiation is the most important variable and most of them predict the right sign of sensitivity to most variables. However, the models predict that biomass reaches maximum levels at unrealistically low levels of iron and unrealistically high levels of light.

Environmental Sciences & Ecology↗

DOE OSTI2025

Files and scripts to support manuscript Needham et al. Canopy Gradients of Respiration

This dataset includes the parameter files, relevant output files, and scripts to perform analysis with Jupyter notebooks that support the manuscript Needham et al 2025 “Canopy Gradients of Respiration Drive Plant Carbon Budgets and Leaf Area Index.” We add functionality to the Functionally Assembled Terrestrial Ecosystem Simulator (FATES) to allow flexible vertical gradients of leaf maintenance respiration (Rdark) and maximum carboxylation rate (Vcmax) through the canopy. We test the sensitivity of FATES to canopy gradients in Rdark, both in global simulations to assess broad scale impacts on leaf area index (LAI) and vegetation carbon, and in single site simulations where we assess impacts on plant functional type (PFT) competitive dynamics. Parameter files are netcdf files that can be converted to human readable .cdl files using NCO tools. Analysis scripts are Jupyter notebook files. These can be opened and run using the open source Jupyter notebook software. Model outputs are netcdf files.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Strangeness in the proton from $W+$ charm production and SIDIS data

We perform a global QCD analysis of unpolarized parton distribution functions (PDFs) in the proton, including new 𝑊+⁢ charm production data from 𝑝⁢𝑝 collisions at the LHC and semi-inclusive pion and kaon production data in lepton-nucleon deep-inelastic scattering, both of which have been suggested for constraining the strange quark PDF. Compared with a baseline global fit that does not include these datasets, the new analysis reduces the uncertainty on the strange quark distribution over the range 0.01 < 𝑥 < 0.3, and provides a consistent description of processes sensitive to strangeness in the proton. Including the new datasets, the ratio of strange to nonstrange sea quark distributions is $R_s = (s + \bar{s})/(\bar{u} +\bar{d})$ $=$ {$0.7⁢2^{+0.52}_{−0.34}, 0.4⁢6^{+0.30}_{−0.20}, 0.3⁢2^{+0.23}_{−0.15}$} for 𝑥 ={$0.01, 0.04, 0.1$} at 𝑄 2 $=$ 4 GeV 2 . The data place more stringent constraints on the strange asymmetry $(s - \bar{s})$, which is found to be consistent with zero in this range.

Anderson, Trey [College of William and Mary, Willi↗

DOE OSTI2026

Kinetics Modeling and Reactor Design Study of Glucose-to-Terpenes Cell-Free Conversion

Cell-free systems offer many advantages over traditional biological conversion by eliminating biological growth constraints. It also offers easy manipulation and finetuning of the reaction conditions for each individual enzyme. The conversion of cellulosic glucose to Limonene, a terpene, is a promising pathway for producing fuels and chemicals. Recent advances in developing cell-free systems focuses on bench scale optimization of terpene yield and to demonstrate its feasibility towards commercialization [1,2]. There is significant knowledge gap regarding reaction kinetics of these cell-free systems to further study how it will perform at larger scale. We present here, our studies on reaction kinetics and reactor design implications of cell-free glucose to Limonene conversion to facilitate the further development and commercialization of this process. We developed a novel kinetic model based on the metabolic-network structure of the cell-free system with multi-substrate reversible Michaelis-Menten rate law. To estimate kinetic parameters for this system of rate equations, we employed Bayesian optimization to perform global search with the assistance of gaussian processes to balance exploration and exploitation. The model parameters estimated showed good results compared with experimental data. The estimated parameters were used to perform sensitivity analysis. We found that Hexokinase is one of the most critical enzymes that affect the conversion of the glucose. We also observed that abundance of co-factors is also critical to the conversion of glucose to limonene. We investigated packed bed reactors with enzymes immobilized on the surface of particles to convert glucose stream into Limonene for larger scale production. The reactor design such as particle size, enzyme loading, and flow rate are found to be critical for improving yields. [1] Dudley, Q.M., Nash, C.J. and Jewett, M.C., 2019. Synthetic Biology, 4(1), p.ysz003. [2] Korman, T.P., Opgenorth, P.H. and Bowie, J.U., 2017. Nature communications, 8(1), p.15526.

09 BIOMASS FUELS↗

DOE OSTI2024

Significance of radiative corrections on measurements of the EMC effect

Deep inelastic scattering (DIS) from nuclear targets probes the parton distribution functions (PDFs) in nuclei. Comparisons of the PDFs from heavy nuclei and the deuteron show deviations that demonstrate a non-trivial nuclear dependence to these distributions, referred to as the EMC effect. A global analysis of the worlds data on the EMC effect reveals tensions between different extractions. Precise measurements at Jefferson Lab, studying the dependence on both the quark momentum fraction, x, and nuclear mass, show systematic discrepancies among experiments, making the extraction of the A dependence of the EMC effect sensitive to the selection of datasets. Further, by comparing various methods and assumptions used to calculate radiative corrections, we have identified differences that, while not large, significantly impact the EMC ratios and show that using a consistent radiative correction procedure resolves this discrepancy, leading to a more coherent global picture, and allowing for a more robust extraction of the EMC effect for infinite nuclear matter.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2023

A constraint on historic growth in global photosynthesis due to rising CO 2

Theory predicts that rising CO 2 increases global photosynthesis, a process known as CO 2 fertilization, and that this is responsible for much of the current terrestrial carbon sink. The estimated magnitude of the historic CO 2 fertilization, however, differs by an order of magnitude between long-term proxies, remote sensing-based estimates and terrestrial biosphere models. Here we constrain the likely historic effect of CO 2 on global photosynthesis by combining terrestrial biosphere models, ecological optimality theory, remote sensing approaches and an emergent constraint based on global carbon budget estimates. Our analysis suggests that CO 2 fertilization increased global annual terrestrial photosynthesis by 13.5 ± 3.5% or 15.9 ± 2.9 PgC (mean ± s.d.) between 1981 and 2020. Our results help resolve conflicting estimates of the historic sensitivity of global terrestrial photosynthesis to CO 2 and highlight the large impact anthropogenic emissions have had on ecosystems worldwide.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

FATES crown damage simulation outputs 2022

This dataset contains outputs from the Functionally Assembled Terrestrial Ecosystem Simulator (FATES) and accompanies the paper "Needham, J.F., Arellano, A., Davies, S.J., Fisher, R.A., Hammer, V., Knox, R., Mitre, D., Muller-Landau, H.C., Zuleta, D., Koven, C.D. Tree crown damage and its effects on forest carbon cycling in a tropical forest, 2022, Global Change Biology". Data are unprocessed netcdf file outputs from simulations that were run to test the effect of a new crown damage module in FATES. Specifically, this data package contains a sensitivity analysis to the carbon cushion parameter damage_Ccushion_ensemble_e1b5bd9_bf013ef_2021-09-02.h0.ensemble.sofar.nc, a sensitivity analysis to the root nitrogen stoichiometry parameter damage_Nstoich_ensemble_e1b5bd9_bf013ef_2021-09-02.h0.ensemble.sofar.nc, a sensitivity analysis to parameters controlling crown damage and recovery damage_recovery_ensemble_e1b5bd9_354f0b0_2021-09-02.h0.ensemble.sofar.nc, and a sensitivity analysis to the number of crown damage bins elm_fates_bci_*_damagebins.Eac53ccb80b-F8f994c29.2022-04-19.elm.h0.fullrun.nc. This data package also contains a high root nitrogen configuration of FATES, including both a control, and a crown damage simulation high_root_N_control_e1b5bd9_354f0b0_2021-09-02.clm2.h0.fullrun.nc and high_root_N_damage_e1b5bd9_354f0b0_2021-09-02.clm2.h0.fullrun.nc. There is an analogous low root nitrogen configuration of FATES, including a control, low_root_N_control_e1b5bd9_bf013ef_2021-09-02.clm2.h0.fullrun.nc a damage only simulation low_root_N_damageonly_e1b5bd9_bf013ef_2021-09-02.clm2.h0.fullrun.nc, a damage plus mortality simulation low_root_N_damage_mort_e1b5bd9_bf013ef_2021-09-02.clm2.h0.fullrun.nc, and a mortality only simulation low_root_N_mort_only_e1b5bd9_ef845c8_2021-09-02.clm2.h0.fullrun.nc. Finally, there is a two PFT simulation in which we test the effect of recovery on competitive dynamics, low_root_N_damage_two_pfts_stoichastic_e1b5bd9_bf013ef_2021-09-10.clm2.h0.fullrun.nc. These simulations test the effect of representing crown damage in FATES, compared with simulations that have an equivalent increase in mortality. Jupyter notebooks to analyse these files can be found at https://github.com/JessicaNeedham/Needham_etal_GCB_2022_FATES_crown_damage.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Joint effects of climate, tree size, and year on annual tree growth derived from tree–ring records of ten globally distributed forests

Tree rings provide an invaluable long-term record for understanding how climate and other drivers shape tree growth and forest productivity. However, conventional tree-ring analysis methods were not designed to simultaneously test effects of climate, tree size, and other drivers on individual growth. This has limited the potential to test ecologically relevant hypotheses on tree growth sensitivity to environmental drivers and their interactions with tree size. Here, we develop and apply a new method to simultaneously model nonlinear effects of primary climate drivers, reconstructed tree diameter at breast height (DBH), and calendar year in generalized least squares models that account for the temporal autocorrelation inherent to each individual tree's growth. We analyze data from 3811 trees representing 40 species at 10 globally distributed sites, showing that precipitation, temperature, DBH, and calendar year have additively, and often interactively, influenced annual growth over the past 120 years. Growth responses were predominantly positive to precipitation (usually over ≥3-month seasonal windows) and negative to temperature (usually maximum temperature, over ≤3-month seasonal windows), with concave-down responses in 63% of relationships. Climate sensitivity commonly varied with DBH (45% of cases tested), with larger trees usually more sensitive. Trends in ring width at small DBH were linked to the light environment under which trees established, but basal area or biomass increments consistently reached maxima at intermediate DBH. Accounting for climate and DBH, growth rate declined over time for 92% of species in secondary or disturbed stands, whereas growth trends were mixed in older forests. These trends were largely attributable to stand dynamics as cohorts and stands age, which remain challenging to disentangle from global change drivers. By providing a parsimonious approach for characterizing multiple interacting drivers of tree growth, our method reveals a more complete picture of the factors influencing growth than has previously been possible.

generalized least squares (GLS)↗

DOE OSTI2026

Dual-loop Solvent-based CCS for Net Negative CO 2 Emissions with Lower Cost

This final technical report details the successful design, construction, and operational validation of an innovative dual-loop CO 2 capture technology designed to achieve deep decarbonization (>99%) from Natural Gas Combined Cycle (NGCC) power plants that results in the electricity with carbon intensity of approximate 42 kg CO 2 -eq/MWh, less than the electricity produced by solar PV. The integrated process couples a primary aqueous solvent (in this project, a water lean solvent – WLS) absorption loop for bulk CO 2 removal with a secondary potassium hydroxide (KOH) polishing loop featuring electrochemical regeneration. This architecture leverages the higher exergy efficiency of the primary loop while utilizing the fastest kinetic of the secondary loop to capture dilute residual CO 2 , achieving an overall capture rate of 99.9% and co-producing pure hydrogen after moisture being condensed and dehydrated. Technical feasibility was established through a comprehensive 2,000-hour experimental campaign on a bench-scale fully-integrated unit (using 4” absorber and 4” stripper) with the feeding flue gas flowrate in the range of 8-20 cfm, confirming the attainment of Technology Readiness Level (TRL) 4. The project executed extensive parametric testing followed by 1,000 hours of continuous steady-state testing, demonstrating exceptional process stability with the electrochemical regenerator exhibiting less than a 10% reduction in electrical conductivity over the duration of the campaign. Operational characterization gathering on the bench unit revealed distinct energy profiles for the hybrid system: the primary loop required approximately 280 kJ mol -1 for bulk removal, while the polishing loop required approximately 1,600 kJ mol -1 specifically when reducing dilute CO 2 concentrations from ~740 ppm down to <50 ppm. (Please note those energy values/numbers can only be viewed as relative relationship and should not be extrapolated as absolute values required for CO 2 capture). Furthermore, dynamic testing validated the system’s flexibility for utility applications, demonstrating a rapid process response time of <30 minutes to changes in flue gas flowrate. Emission monitoring confirmed that the dual-loop architecture effectively mitigates solvent losses, utilizing a water wash to remove entrained aerosols to <1 ppm. The Techno-Economic Analysis (TEA) indicates a cost of capture of $\$$59.3/tonne and a Levelized Cost of Electricity (LCOE) of 71.7 $\$$/MWh at the overall capture efficiency of 99.8% of total carbon in the flue gas stream, with sensitivity analysis identifying an economic optimum when the primary loop captures 97% of the total CO 2 . A Life Cycle Assessment (LCA) confirms the technology’s potential for net-negative emissions, determining a Global Warming Potential (GWP) of 52 kg CO 2 -eq/MWh—significantly lower than the baseline—which further decreases to 42 kg CO 2 -eq/MWh when crediting the displacement of conventional hydrogen production.

03 NATURAL GAS↗

DOE OSTI2022

On the Convergence of Overlapping Schwarz Decomposition for Nonlinear Optimal Control

Here, we study the convergence properties of an overlapping Schwarz decomposition algorithm for solving nonlinear optimal control problems (OCPs). The algorithm decomposes the time domain into a set of overlapping subdomains, and solves all subproblems defined over subdomains in parallel. The convergence is attained by updating primal-dual information at the boundaries of overlapping subdomains. We show that the algorithm exhibits local linear convergence, and that the convergence rate improves exponentially with the overlap size. We also establish global convergence results for a general quadratic programming, which enables the application of the Schwarz scheme inside second-order optimization algorithms (e.g., sequential quadratic programming). The theoretical foundation of our convergence analysis is a sensitivity result of nonlinear OCPs, which we call "exponential decay of sensitivity" (EDS). Intuitively, EDS states that the impact of perturbations at domain boundaries (i.e., initial and terminal time) on the solution decays exponentially as one moves into the domain. Here, we expand a previous analysis available in the literature by showing that EDS holds for both primal and dual solutions of nonlinear OCPs, under uniform second-order sufficient condition, controllability condition, and boundedness condition. We conduct experiments with a quadrotor motion planning problem and a partial differential equations (PDE) control problem to validate our theory, and show that the approach is significantly more efficient than alternating direction method of multipliers and as efficient as the centralized interior-point solver.

42 ENGINEERING↗

DOE OSTI2020

Understanding cloud droplet spectral dispersion effect using empirical and semi-analytical parameterizations in NCAR CAM5.3

Five parameterizations of cloud droplet spectral shape are implemented in a global climate model to investigate the dispersion effect and aerosol indirect effect (AIE). We design a series of experiments by modifying the microphysical cloud scheme of NCAR CAM5.3 (Community Atmospheric Model Version 5.3). We employ four empirical (Martin94, RLiu03, PengL03, Liu08) and one semi-analytical (LiuLi15) expressions for cloud droplet spectral shape parameters. Analysis focuses on the instantaneous differences in the simulated cloud microphysical properties and the comparison between model output and satellite data. The results show that RLiu03, PengL03 and LiuLi15 produce wider droplet spectrum and faster autoconversion rate, but Liu08 has a narrower droplet spectrum and slower autoconversion rate than the default parameterization (Martin94) in CAM5.3. Global dispersion effects caused by the five parameterizations modify the aerosol indirect effect by -10% (counteract) to 13% (strengthen). The simulated AIEs and dispersion effects exhibit noticeably spatial inhomogeneity. In the sensitive regions of AIE (Southeast Asia, North Pacific and west coast of South America), we decompose the response of shortwave cloud forcing to the change in droplet number for analysis. The varying dispersion effects can be explained by different responses of cloud properties in different spectral parameterizations.

54 ENVIRONMENTAL SCIENCES↗