SEARCH · Engineering Papers
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
QSPR analyses can be used to identify useful descriptors leading to statistical models for membrane crossover. This data-driven approach can be used to evaluate ROMs for asymmetric non-aqueous redox flow batteries.
In this paper, we propose a variational approach to estimate eddy viscosity using forward sensitivity method (FSM) for closure modeling in nonlinear reduced order models. FSM is a data assimilation technique that blends model's predictions with noisy observations to correct initial state and/or model parameters. We apply this approach on a projection based reduced order model (ROM) of the one-dimensional viscous Burgers equation with a square wave defining a moving shock, and the two-dimensional vorticity transport equation formulating a decay of Kraichnan turbulence. We investigate the capability of the approach to approximate an optimal value for eddy viscosity with different measurement configurations. Specifically, we show that our approach can sufficiently assimilate information either through full field or sparse noisy measurements to estimate eddy viscosity closure to cure standard Galerkin reduced order model (GROM) predictions. Therefore, our approach provides a modular framework to correct forecasting error from a sparse observational network on a latent space. We highlight that the proposed GROM-FSM framework is promising for emerging digital twin applications, where real-time sensor measurements can be used to update and optimize surrogate model's parameters.
This study explores the added value (AV) of a regional Earth system model (ESM) compared to an atmosphere-only regional climate model (RCM) in simulating West African monsoon (WAM) rainfall. The primary goals are to foster discussions on the suitability of coupled RCMs for WAM projections and deepen our understanding of ocean–atmosphere coupling’s influence on the WAM system. The study employs results from dynamical downscaling of the ERA-Interim reanalysis and Max Plank Institute ESM, low resolution (MPI-ESM-LR), by two RCMs, atmosphere only (REMO) and REMO coupled with Max Planck Institute Ocean Model (MPIOM) (ROM), at ~25-km horizontal resolution. Results show that in regions distant from coupling domain boundaries such as West Africa (WA), constraint conditions from ERA-Interim are more beneficial than coupling effects. REMO, reliant on oceanic sea surface temperatures (SSTs) from observations and influenced by ERA-Interim, is biased under coupling conditions, although coupling offers potential advantages in representing heat and mass fluxes. Contrastingly, as intended, coupling improves SSTs and monsoon fluxes’ relationships under ESM-forced conditions. In this latter case, the coupling features a dipole-like spatial structure of AV, improving precipitation over the Guinea Coast but degrading precipitation over half of the Sahel. Our extensive examination of physical processes and mechanisms underpinning the WAM system supports the plausibility of AV. Additionally, we found that the monsoonal dynamics over the ocean respond to convective activity, with the Sahara–Sahel surface temperature gradient serving as the maintenance mechanism. While further efforts are needed to enhance the coupled RCM, we advocate for its use in the context of WAM rainfall forecasts and projections.
A challenge for geophysical research is the relatively short time for which continuous digital data have been available. There are many events of research interest where much of the data were collected by analog or temporary digital means. These legacy data require special efforts to find, organize, and preserve for future research. Here, we focus on one type of legacy data, nuclear explosions, and associated events (e.g., collapses, aftershocks) of high interest for nuclear monitoring research and development. For nuclear tests conducted in Nevada, the Lawrence Livermore Nevada Network recorded valuable data on a variety of formats. In this article, we briefly review this network, and describe one set of data that previously was available via CD-ROM by request, and since 2018 has become available through Incorporated Research Institutions for Seismology as assembled dataset 18-001.
This is a collaborative effort between Battelle Memorial Institute as manager and operator of Pacific Northwest National Laboratory (PNNL) and ESI North America Inc. (“ESI” or “Participant”) to apply computation and data analytics to the challenge of light weighting with a focus on the battery enclosures of electric vehicles (EVs). EVs use heavy batteries to increase range and power. A complex-shaped battery enclosure is required to meet a host of challenging performance requirements. The ability to virtually develop composite parts such as battery enclosure with tailored properties to meet required performance will be highly valuable to the automotive industry. However, efficient simulation of composite-manufacturing processes remains a challenging issue since simulation involves multiscale models in space and time, highly non-linear and anisotropic behavior, strongly coupled multi-physics, and complex geometries. This work will advance the state of the art by reducing the computational burden of composite optimization by using simulation data from a limited number of configurations off-line and then developing a reduced order model (ROM) using data analytics and machine learning (ML). Develop a data driven approach to link features of the material and manufacturing processes to the mechanical properties of thermoplastic composite parts.
This document pertains to the reporting requirements of DOE contract FE-0031767. The document covers the development of a model predictive control (MPC) library for implementing MPC for a general dynamic system. The library is implemented in a standardized manner in Matlab/Simulink, where core functions on model prediction, linearization and formulation and solution of a quadratic programming (QP) optimization problem is done in the core library - independent of the specific application. The user can provide the application-specific dynamic model in continuous and discrete time, to rapidly implement and test the MPC performance in a desktop simulation. The MPC optimization objective and constraints are also easily configured via an Excel file to allow iterative refinement as needed. Finally, the MPC library enables a rapid deployment to a target environment through auto C-code generation and containerization. The MPC library works seamlessly with the model based estimation (MBE) library to obtain the overall output feedback control solution. In this program, the reduced order model (ROM) of a coal-fired power plant (CFPP) is used to implement and test the MPC solution.
Integrated energy systems (IES) combine, in mutually beneficial ways, power from variable renewable energy sources and nuclear power plants (NPP) to improve economic viability under uncertain market and weather conditions. The open-source Framework for Optimization of Resources and Economics (FORCE) tool suite, developed at Idaho National Laboratory (INL), has enabled comprehensive modeling and simulation of IES. The capabilities within FORCE include grid portfolio optimization through the Holistic Energy Resource Optimization Network (HERON) and the transient process model analysis library HYBRID, among others. Continuous efforts and investments from the IES programs have been made to expand and improve the versatility of the FORCE toolset in fiscal year 2022. Code-coupling and cross-tool communication have been important methods for improving this versatility. This report focuses on an additional workflow in the HERON tool for capacity and dispatch stochastic optimization through integration with the external tool Design Integration and Synthesis Platform to Advance Tightly Coupled Hybrid Energy Systems (DISPATCHES). DISPATCHES was primarily developed by the National Energy Technology Laboratory, in collaboration with other national laboratories, which included INL, universities, and industry partners. It is coupled to a library of algebraic models for specific plant components, and to a framework for stochastic optimization different from that provided in the current Risk Analysis Virtual Environment (RAVEN)-running-RAVEN algorithm in HERON. HERON currently conducts stochastic optimization via an outer-inner loop: it optimizes over variable capacity on the outer loop, and at each step within the capacity parameter space, conducts an inner optimization over scenarios (of market signals, demand, and/or weather patterns) and hourly dispatch throughout a user-specified number of years. On the other hand, DISPATCHES conducts stochastic optimization via an “all-at-once” strategy in which capacity variables are optimized at the same level as dispatch variables, as all scenarios are considered at once. The latter method works especially well for projects of limited size and project length, as the necessary computational power and memory increases with the number of variables and scenarios. The new capability to use the DISPATCHES workflow in HERON enhances standalone simulations by leveraging FORCE tools—namely, the economic metrics from the Tool for Economic Analysis (TEAL) and reduced-order model (ROM) sampling from RAVEN. The initial demonstration of the DISPATCHES workflow simulates an existing nuclear-case flowsheet within the DISPATCHES repository—this models a NPP with a secondary revenue stream for hydrogen production. Electrical output from the plant is converted to hydrogen via a proton-exchange membrane (PEM) electrolyzer, hydrogen tanks are used for storage, and an additional turbine is added for hydrogen combustion. Continued work regarding this FORCE-DISPATCHES integration will include automatic generation of DISPATCHES models from HERON inputs, offering analysts the option of using either the RAVEN-runsRAVEN or DISPATCHES workflow to solve technoeconomic optimization problems.
This report provides a brief description on the use of the Seal_Flux computer program developed as part of the effort to quantify the risk of geologic storage of carbon dioxide (CO 2 ) under the U.S. Department of Energy’s (DOE) National Risk Assessment Partnership (NRAP). The Seal_Flux code simulates the flow of CO 2 through a low permeability rock horizon or seal formation overlying the storage reservoir into which CO 2 is injected. A two-phase, relative permeability approach with Darcy’s law is used for one-dimensional (1D) flow computations of CO 2 through the horizon in the vertical direction. The code also allows the simulation of time-dependent processes that can influence such flow. However, as part of its design, Seal_Flux is what can be termed a “reduced-order model” (ROM) and is not intended as a full-functioning flow code. The theory and simulation in the code is streamlined and directed towards the implementation of Monte Carlo risk analyses of CO 2 transport or as termed in this context as “leakage.” While presented in this report as a stand-alone tool, the Seal_Flux code is intended to function in the future as one of several models as part of an integrated, systems-level model of CO 2 storage performance. Finally, the code is written in Python 3.10 to provide an open framework for further development by others and to assist in linking the code with other modules in an integrated assessment model.
This work package is focused on creating a digital twin (DT) of the Facility to Alleviate Salt Technology Risks (FASTR) experiment. FASTR is a high-profile facility tasked with maturing molten salt applications such as liquid salt energy storage, concentrated solar power facilities, and molten salt cooling systems supporting other advanced reactor technologies, all of which may be incorporated as components within the Integrated Energy Systems (IES) program. As part of creating a digital twin (DT), verification and validation (V&V) has been performed for the physics-based model, and gaps have been identified. Additionally, a preliminary integration of the FASTR model has been completed with the existing physics-based model of the Thermal Energy Distribution System (TEDS)/Microreactor AGile Non-nuclear Experimental Testbed (MAGNET) facility at Idaho National Laboratory (INL). A Functional Mock-up Unit (FMU) of the physics-based system model and a reduced-order model (ROM) were both used to test for a preliminary hardware in the loop operation when the physical loop was offline. Future tasks could include improving the physics-based model to better match the available experimental data, improving the hardware in the loop integration, and improving the integration with the TEDS model.
This document demonstrates completion of the goals described in the technical narrative of the Department of Energy’s Industry Funding Opportunity Announcement (iFOA) project entitled “Modeling and Simulation Development Pathways to Accelerating KP-FHR Licensing,” which relates to the development and demonstration of capabilities for conducting engineering-scale simulations of Alloy 316H components under high-temperature conditions. This work encompassed two major aspects: integrating and testing constitutive models for the creep response of 316H at high temperatures, and developing tools for modeling creep crack growth (CCG) in 316H. Two classes of constitutive models were used in this effort: phenomenological models based on behavior observed at the engineering scale, and reduced-order models (ROMs) that represent the nonlinear response of mesoscale models that capture the sensitivity to material microstructure and processing. Likewise, the approaches employed for CCG modeling considered both simplified engineering approaches and detailed simulations of creep and damage ahead of the crack tip. These developments, which were performed by utilizing the Grizzly code as well as the open-source libraries it depends on, strengthen Grizzly’s ability to support licensing and safety analyses of high-temperature reactor components.
RTRC and UMD have developed novel machine learning based methods under the ARPA-E DIFFERENTIATE program for rapid acceleration of hypothesis generation in complex architecture design spaces involving both discrete choices of component inclusion and interconnection and continuous parametric decisions. The project named Learning Enabled Network Synthesis (LENS) further demonstrated the developed methods on challenging electrical power converter design problems by identifying the most suitable circuit topologies and simultaneously selecting the most appropriate components to achieve optimized design of power converter with improved performances. We demonstrated that LENS could enable exploration of very large design space of circuit topologies and components by addressing the limitations of conventional design process in non-linear, high switching speed, multi-dimensional power converter design and optimization. The key innovation developed in LENS is the seamless integration of statistical learning and logical reasoning techniques and building on the individual strengths of these techniques for rapid hypothesis discovery. The main component of LENS comprises of: 1) Graph Reasoning Engine (GRE) to enforce composition rules that rapidly reject all discrete architectures that are composed incorrectly and generates an adaptive database of feasible designs which can be used by ML modules, 2) Graph Generative Learning module which is a deep neural network based generative model for graph architectures which can enable design space exploration beyond the dataset generated by the GRE, 3) Graph Reduced Order Model (ROM) for graph domains for accelerating computation of output metrics, and 4) Active learning and Rule Discovery module for sample efficient learning and extracting logical rules from the learned ML models which will be integrated in the GRE to enhance the filtering effectiveness. LENS approach can be applied to any design domains where designs can be represented as multi-attribute graphs. The LENS team integrated the various technical innovations listed above into an optimization pipeline and exercised the optimization pipeline on the converter design problem. The LENS project demonstrated that the developed AI/ML technologies can be used to generate novel converter circuits >45x faster than experts on chosen use-cases. This can enable faster design space exploration and identification of new designs which are not considered by experts due to the increasing design space complexity. This has significant potential impact on the public and energy needs of the country. It is currently estimated that 30% of all electrical powers generated passes through power converters. The future estimate is that 80% of all power generated would be passing through converters. LENS fills a critical gap in this space since by accelerating the design process the designers would be able to generate more efficient converters which can lead to significant energy savings for the country.
Hydrogen utilization at Plug Power sites ranges between 70-93% with the remainder vented and lost to the atmosphere. The goal of this project was to improve hydrogen utilization at Plug Power fulfillment centers via the patented Heisenberg vortex tube (HVT). The HVT combines conventional counterflow vortex tubes with para-orthohydrogen conversion catalyst to allow the cooling of hydrogen flows with no moving parts or external energy inputs. In year 1 of this project the HVT was evaluated to improve hydrogen utilization via the following concepts: 1) increasing liquid hydrogen pump volumetric efficiency by 20% through vapor separation and subcooling of the liquid, 2) reducing liquid hydrogen storage tank boil-off losses by 20% through thermal vapor shielding (TVS), and 3) increase isentropic efficiency of supercritical hydrogen expansion at 40-50 K, from 31% to over 40%. The TVS concept was selected as the most favorable application for further development. While the conceptual analysis was underway, the cryocatalysis hydrogen experiment facility (CHEF) was retrofitted with a new cryocooler, higher pressure condenser tanks, and in-situ fiber-optic Raman probes for ortho-parahydrogen composition analysis at the inlets and outlets of the HVT. To our knowledge this is the first in-situ implementation of cryogenic Raman probes for ortho-parahydrogen analysis. Subsequent testing of catalyzed and non-catalyzed HVT identified a low flow rate, high-conversion efficacy regime well suited for TVS development. Computational fluid dynamics (CFD), Reduced order modeling (ROM), and quantum Monte Carlo atomistic simulations were applied to optimize the design of the HVT implemented for TVS. The simulations and analysis identified a ruthenium-based catalyst as the most optimal for high conversion with little pressure losses and was matched to experimental measurements. These results indicated the design could achieve the 20% boil-off reduction target. An HVT field-trial was designed and constructed by Plug Power and implemented at a customer site. However, budget limitations reduced the amount of catalyst that could be applied to the HVT. With the reduced amount of catalyst, CFD analysis estimated about 3% reduction in boil-off. A similar amount of boil-off reduction was observed during preliminary measurements during tank commissioning. However, the results have a large margin for error and were operated at relatively low tank temperatures when the HVT has the lowest potential for improvement. This project demonstrated the use of para-orthohydrogen catalysis for reducing liquid hydrogen tank boil-off losses. The field trial system is anticipated to provide long-term experimental measurements on tank operational performance that will allow Plug Power to instrument additional tanks in the future. Recommendations for future work include the development of higher activity, lower cost para-orthohydrogen conversion catalyst for use in the HVT. Tank process optimizations could maximize the efficacy of the HVT and further improve hydrogen utilization.
This report describes work performed during FY2024 under the auspices of the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program to inform and improve mechanistic models of fuel-cladding chemical interaction (FCCI) in metallic fuel. For fuel-side FCCI, atomistic simulations were performed to determine the diffusivity of iron (Fe) in the $\alpha$ and $\gamma$ phases of uranium (U). A model of liquid penetration of cladding due to melting of the fuel-side FCCI region was updated to account for the finite size of the FCCI region, and the model was validated through comparison with tests performed in the Fuel Behavior Test Apparatus (FBTA). For cladding wastage formation, a reduced-order model was improved by comparison with a multi-scale mechanistic model to better quantify the ROM parameters.
This study, performed on behalf of Washington River Protection Solutions (WRPS), examines the potential for making cost savings for the Immobilized Low Active Waste (ILAW) containers over the mission life of the Waste Treatment and Immobilization Plant (WTP). it builds on earlier study assessments 24590-LAW-RPT-M-01-001 Rev 0 ILAW Product Container Specifications Optimization Study and DOE/ORP-2013-02 Rev 0 Decision Analysis of Low-Active Waste Container Finishing Handling System Alternatives. The WTP design and installed equipment has evolved over an extended period ~20 years. Consequently, some of the existing requirements relating the ILAW containers may no longer be applicable under current mission plans where the glass filled containers are placed into the Integrated Disposal Facility (IDF) rather than stacked in an engineered purpose-built ventilated store as originally envisioned. The study team reviewed available WTP project reference documents to establish the functions the container must perform and assign requirements against each of these functions. With the assumptions that some of the existing WTP contract restrictions would not necessarily apply to the subsequent facility operations contract the container requirements were reviewed to identify those which were ‘hard constraints‘, imposed by regulatory requirements. Initially each of these hard constraints were reviewed, evaluated and documented to see if they precluded the use of an alternate less expensive material of construction, carbon steel, to the existing 304 L stainless steel in the current container specification. That review did not generate any fatal flaws for a change in material of construction to carbon steel allowing the study team to examine the impacts of material cost savings, potential changes to the design and requirements and impacts on interfaces of the containers with the existing WTP process equipment from container receipt, transfers, glass filling, cooling and export. In addition, a change in container target fill level was evaluated and qualitative, and quantitative where feasible, impacts to the container fabrication process assessed. Lastly, additional follow on development or testing activities to support the use of carbon steel were assessed and a recommended path forward with a ROM cost and timescale developed.
This document summarizes the findings of a strategic planning exercise commissioned by the Weapons Simulation and Computing, Computational Physics (WSC/CP) program at the Lawrence Livermore National Laboratory (LLNL) in FY24. During the year, the committee met with multiple stakeholder communities to gather input, opinions, suggestions and concerns which have been incorporated throughout this document. The key findings from this exercise are summarized: • The development of multiphysics modelling and simulation (mod/sim) codes and software technologies, their deployment on exascale compute platforms, and their broad adoption across the NNSA is a major success of the Advanced Simulation and Computing (ASC) program and the Exascale Computing Project (ECP). Sustained investment in these core technologies is essential. • Today’s state of the art involves running ensembles of O(100K) simulations to perform uncertainty quantification (UQ) and design studies using multiple statistical methods such as Bayesian optimization to understand sensitivities of our models and explore parameterized design spaces. Even with exascale computing, we are practically limited to O(10) parameters in these studies since the number of simulations required to sample the space scales exponentially with the number of design parameters. • The data from these simulation ensembles is increasingly being used to train machine learned (ML) surrogates (or reduced order models, ROMs) which can then be used for optimization or real time design exploration. However, the trained surrogates are still limited in the number of parameters they can represent due to the sampling limitations previously noted. • Augmenting our suite of integrated multiphysics simulation codes, both current and emerging, with the ability to compute gradients (solution derivatives) of arbitrary simulation outputs with respect to (some or all) simulation inputs would be a breakthrough technology, opening the door to a new era of efficient and automated inverse design based on verified and validated mod/sim capabilities. • This capability, which we refer to as differentiable multiphysics codes (DMCs), would revolutionize both UQ and optimization studies by breaking the curse of dimensionality that presently limits our “gradient-free” ensemble based computing approach. A similar breakthrough occurred in the AI/ML community once the ability to compute gradients of arbitrary loss functions using back-propagation became commonplace. Gradient information from the multiphysics codes can also be used to dramatically improve the efficiency and scale of training of ML/ROM surrogates for rapid assessments. • Achieving this in our suite of codes will be a grand challenge, similar to the amount of effort that was required to transition from CPU to GPU computing. It will require buy-in from the entire WSC/CP program and beyond, including all integrated codes, physics and engineering models, third-party library dependencies and performance portability abstractions. It will also require investment in research and development of numerical methods for computing adjoints of coupled physics across multiple adaptively refined moving meshes and of stochastic (Monte Carlo) and mesh free (SPH) methods. • New software and numerical techniques, largely pioneered by the AI/ML community, make this feasible. Chief among these is automatic differentiation (AD), the ability to employ AD at point-wise locations in a physics calculation (instead of traditional black-box approaches) and the ability to perform “back-propagation in time” (or reverse mode AD) for non-linear partial differential equations (PDEs). Fundamentally, the conclusion of this strategic planning exercise is that the time is right to undertake a large scale effort in WSC, centered on the existing integrated codes, to continue the natural evolution of mod/sim in the age of AI/ML. Instead of attempting to replace mod/sim with purely data driven AI/ML models, we believe the key to success is to integrate AI/ML by building on top of the decades of hard-won knowledge and the verified/validated multiphysics modelling capability that is the hallmark of the ASC program.
One of the practical challenges of the studies of rotating detonation engines (RDEs) is the direct estimation of losses from experimental measurements. This study attempts at resolving this limitation by combining a reduced order model (ROM) of the detonation wave characteristics with a Markov chain Monte Carlo parameter estimation framework. The model considers simple deflagration losses and the overall impact of deflagration on RDE performance. To evaluate this model, a Markov Chain Monte Carlo (MCMC) sampling technique was applied to estimate the loss parameters within the model for a set of conditions operated in hydrogen-air over a range of mass flow rates and equivalence ratios. The MCMC parameter estimation framework allowed for the determination of a posterior distribution of the loss parameters for each test condition, an examination of the correlation between the loss parameters and measured performance metrics of the RDE, and an uncertainty propagation of these parameters. The predicted model loss parameters were then compared to indirect experimental measurements of the deflagration combustion fractions to evaluating the validity of the approach, and shed light on the benefits and drawbacks of the model, measurement techniques, and the estimation framework.
In this work, a dataset including structural and mechanical properties of refractory multicomponent alloys was developed by fusing computations of phase diagram (CALPHAD) and density functional theory (DFT). The refractory multicomponent alloys, also named refractory complex concentrated alloys (CCAs) which contain 2–5 types of refractory elements were constructed based on Special Quasi-random Structure (SQS). The phase of alloys was predicted using CALPHAD and the mechanical property of alloys with stable and single body-centered cubic (BCC) at high temperature (over 1,500°C) was investigated using DFT-based simulation. As a result, a dataset with 393 refractory alloys and 12 features, including volume, melting temperature, density, energy, elastic constants, mechanical moduli, and hardness, were produced. To test the capability of the dataset on supporting machine learning (ML) study to investigate the property of CCAs, CALPHAD, and DFT calculations were compared with principal components analysis (PCA) technique and rule of mixture (ROM), respectively. It is demonstrated that the CALPHAD and DFT results are more in line with experimental observations for the alloy phase, structural and mechanical properties. Furthermore, the data were utilized to train a verity of ML models to predict the performance of certain CCAs with advanced mechanical properties, highlighting the usefulness of the dataset for ML technique on CCA property prediction.
Major concerns of carbon dioxide (CO2) sequestration in subsurface formations are knowledge of the well injectivity and gas storage capacity of the formation, the CO2 pressure and saturation plume extensions during and after injection, and the risks associated with CO2 leakage and fault reactivation. Saline reservoirs are considered as one of the target formations for CO2 sequestration through structural, residual, dissolution, and mineral trapping mechanisms. The boundary condition of the saline reservoir dictates the pressure and saturation plume extension of the injected supercritical CO2 that could expand over large distances. This can lead to sources of risk, e.g., leakage and/or fault reactivation due to presence of wells, thief zones, and geological discontinuities. Therefore, there is a critical need to develop a model that describes how risk-related performance metrics (i.e., the CO2 saturation plume size, the pressure differential plume area, and the pressure differential at specific locations) vary as a function of the size of injection, time following injection, injection operations, and geologic environment. In this study, a systematic reservoir modeling studies of anthropogenic CO2 sequestration in Citronelle dome, Alabama, was performed where all relevant scenarios and conditions to address the questions of the saturation and pressure plume size in the area of review (AoR) and post-injection site care (PISC) are considered. The objective for this study was firstly to systematically simulate CO2 sequestration, i.e., saturation dynamics, and pressure behavior over a range of operational and geological conditions and to derive conclusions about the factors influencing saturation and pressure plume size, post-injection behavior, and the risk associated with them, by developing third-generation reduced order models (ROMs) for reservoir behavior. Finally, to assess the uncertainty associated with our studies, Latin Hypercube Sampling (LHS) together with an experimental design technique, i.e., Plackett–Burman design, was used. Application of Pareto charts and respond surfaces enabled us to determine the most important parameters impacting saturation and pressure plume sizes and to quantify the auto- and cross-correlation among different parameters in both history-matched and upscaled models.