SEARCH · Engineering Papers

Results for “data conditioning”

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 1,189 records · Page 66

All sources NASA NTRS DOE OSTI

DOE OSTI2022

Photovoltaic System Health-State Architecture for Data-Driven Failure Detection

The timely detection of photovoltaic (PV) system failures is important for maintaining optimal performance and lifetime reliability. A main challenge remains the lack of a unified health-state architecture for the uninterrupted monitoring and predictive performance of PV systems. To this end, existing failure detection models are strongly dependent on the availability and quality of site-specific historic data. The scope of this work is to address these fundamental challenges by presenting a health-state architecture for advanced PV system monitoring. The proposed architecture comprises of a machine learning model for PV performance modeling and accurate failure diagnosis. The predictive model is optimally trained on low amounts of on-site data using minimal features and coupled to functional routines for data quality verification, whereas the classifier is trained under an enhanced supervised learning regime. The results demonstrated high accuracies for the implemented predictive model, exhibiting normalized root mean square errors lower than 3.40% even when trained with low data shares. The classification results provided evidence that fault conditions can be detected with a sensitivity of 83.91% for synthetic power-loss events (power reduction of 5%) and of 97.99% for field-emulated failures in the test-bench PV system. Finally, this work provides insights on how to construct an accurate PV system with predictive and classification models for the timely detection of faults and uninterrupted monitoring of PV systems, regardless of historic data availability and quality. Such guidelines and insights on the development of accurate health-state architectures for PV plants can have positive implications in operation and maintenance and monitoring strategies, thus improving the system’s performance.

photovoltaics↗

DOE OSTI2020

NSFnets (Navier-Stokes flow nets): Physics-informed neural networks for the incompressible Navier-Stokes equations

In the last 50 years there has been a tremendous progress in solving numerically the Navier-Stokes equations using finite differences, finite elements, spectral, and even meshless methods. Yet, in many real cases, we still cannot incorporate seamlessly (multi-fidelity) data into existing algorithms, and for industrial-complexity applications the mesh generation is time consuming and still an art. Moreover, solving ill-posed problems (e.g., lacking boundary conditions) or inverse problems is often prohibitively expensive and requires different formulations and new computer codes. Here, we employ physics-informed neural networks (PINNs), encoding the governing equations directly into the deep neural network via automatic differentiation, to overcome some of the aforementioned limitations for simulating incompressible laminar and turbulent flows. We develop the Navier-Stokes flow nets (NSFnets) by considering two different mathematical formulations of the Navier-Stokes equations: the velocity-pressure (VP) formulation and the vorticity-velocity (VV) formulation. Since this is a new approach, we first select some standard benchmark problems to assess the accuracy, convergence rate, computational cost and flexibility of NSFnets; analytical solutions and direct numerical simulation (DNS) databases provide proper initial and boundary conditions for the NSFnet simulations. The spatial and temporal coordinates are the inputs of the NSFnets, while the instantaneous velocity and pressure fields are the outputs for the VP-NSFnet, and the instantaneous velocity and vorticity fields are the outputs for the VV-NSFnet. This is unsupervised learning and, hence, no labeled data are required beyond boundary and initial conditions and the fluid properties. The residuals of the VP or VV governing equations, together with the initial and boundary conditions, are embedded into the loss function of the NSFnets. No data is provided for the pressure to the VP-NSFnet, which is a hidden state and is obtained via the incompressibility constraint without extra computational cost. Unlike the traditional numerical methods, NSFnets inherit the properties of neural networks (NNs), hence the total error is composed of the approximation, the optimization, and the generalization errors. Here, we empirically attempt to quantify these errors by varying the sampling (“residual”) points, the iterative solvers, and the size of the NN architecture. For the laminar flow solutions, we show that both the VP and the VV formulations are comparable in accuracy but their best performance corresponds to different NN architectures. The initial convergence rate is fast but the error eventually saturates to a plateau due to the dominance of the optimization error. For the turbulent channel flow, we show that NSFnets can sustain turbulence at , but due to expensive training we only consider part of the channel domain and enforce velocity boundary conditions on the subdomain boundaries provided by the DNS data base. We also perform a systematic study on the weights used in the loss function for balancing the data and physics components, and investigate a new way of computing the weights dynamically to accelerate training and enhance accuracy. In the last part, we demonstrate how NSFnets should be used in practice, namely for ill-posed problems with incomplete or noisy boundary conditions as well as for inverse problems. We obtain reasonably accurate solutions for such cases as well without the need to change the NSFnets and at the same computational cost as in the forward well-posed problems. As a result, we also present a simple example of transfer learning that will aid in accelerating the training of NSFnets for different parameter settings.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2023

Physics model validation of propane and methane for Hydrogen Plus Other Alternative Fuels Risk Assessment Models (HyRAM+)

HyRAM+ is a toolkit that includes fast-running models for the unconstrained (i.e., no wall interactions) dispersion and flames for non-premixed fuels. The models were developed for use with hydrogen, but the toolkit was expanded to include propane and methane in a recent release. Here, in this work we validate the dispersion and flame models for these additional fuels, based on reported literature data. The validation efforts spanned a range of release conditions, from subsonic to underexpanded jets and flames for a range of mass flow rates. In general, the dispersion model works well for both propane and methane although the width of the jet/plume is predicted to be wider than observed in some cases. The flame model tends to over-predict the induced buoyancy for low-momentum flames, while the radiative heat flux agrees with the experimental data reasonably well, for both fuels. The models could be improved but give acceptable predictions for propane and methane behavior for the purposes of risk assessment.

33 ADVANCED PROPULSION SYSTEMS↗

DOE OSTI2025

Towards provision of regularly updated climate data from the Coupled Model Intercomparison Project

The Coupled Model Intercomparison Project (CMIP) is a flagship of the World Climate Research Programme (WCRP). CMIP has become a recognised ‘brand’ in climate circles evolving over the last thirty years from a targeted research activity by a small number of climate modelling centres intercomparing their Earth System Model (ESM) simulations to a broad international coordinated research effort (Durack et al, 2025). CMIP is organized as a research activity leveraging funded and in-kind contributions from experts within modelling centres and the broader scientific community supported more recently by a fully-funded International Project Office. Within CMIP, Model Intercomparison Projects (MIPs) are community-designed to understand past, present and future climate. CMIP data provides a valuable resource for climate research and is routinely used to assess model representation of climate processes and test scientific hypotheses in the context of model uncertainty and (forced and internal) variability as evident from its prolific use in scientific publications1 . The impact relies on enabling infrastructure (most prominently via the Earth System Grid Federation (ESGF)), which allows sharing of simulation output, provision of the boundary conditions used in each simulation, and definition of the data standards that are essential to facilitating wide use of the data. The impact is supplemented by the wide-ranging scrutiny to which model simulations are subjected. Beyond its use in research, CMIP data is a key resource for communities producing derived climate information from downscaling and impact studies, such as the Coordinated Regional Downscaling Experiment (CORDEX; Gutowski et al., 2016) and the Intersectoral Impacts MIP (ISIMIP; Frieler et al., 2024). Government, academic and commercial entities also increasingly rely on CMIP and its downstream data for climate risk assessments and climate services (for example, Copernicus Climate Change Service and World Bank portal). This means that, although CMIP is a research activity, it increasingly serves a secondary and very relevant role as a provider of climate data – a long-recognised dichotomy (Stevens, 2024). Research and applications have distinct needs, with the former requiring flexibility and generality and the latter consistency. Here we explain how the design of the research activity has been adapted to reduce the burdens imposed by applications and how the research infrastructure might evolve to further enable scientific inquiry. We propose one possible approach to consistently providing model information and projections for applications in the future.

Environmental sciences↗

DOE OSTI2026

Conditional distribution estimation of building characteristics with diffusion models for urban energy modeling

Understanding current energy consumption behavior in communities is critical for informing future energy use decisions and enabling efficient energy management. Urban energy models, which are used to simulate these energy use patterns, require large datasets with detailed building characteristics for accurate outcomes. However, such detailed characteristics at the individual building level are often unknown and costly to acquire, or unavailable. Through this work, we propose using a generative modeling approach to generate realistic building attributes to fill in the data gaps and finally provide complete characteristics as inputs to energy models. Our model learns complex, building-level patterns from training on a large-scale residential building stock model containing 2.2 million buildings. We employ a tabular diffusion-based framework that is designed to handle heterogeneous (discrete and continuous) features in tabular building data, such as occupancy, floor area, heating, cooling, and other equipment details. We develop a capability for conditional diffusion, enabling the imputation of missing building characteristics conditioned on known attributes. We conduct a comprehensive validation of our conditional diffusion model, firstly by comparing the generated conditional distributions against the underlying data distribution, and secondly, by performing a case study for a Baltimore residential region, showing the practical utility of our approach. Our work is one of the first to demonstrate the potential of generative modeling to accelerate building energy modeling workflows.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2021

Plug-and-play polymer microfluidic chips for hydrated, room temperature, fixed-target serial crystallography

The practice of serial X-ray crystallography (SX) depends on efficient, continuous delivery of hydrated protein crystals while minimizing background scattering. Of the two major types of sample delivery devices, fixed-target devices offer several advantages over widely adopted jet injectors, including: lower sample consumption, clog-free delivery, and the ability to control on-chip crystal density to improve hit rates. In this study we present our development of versatile, inexpensive, and robust polymer microfluidic chips for routine and reliable room temperature serial measurements at both synchrotrons and X-ray free electron lasers (XFELs). Our design includes highly X-ray-transparent enclosing thin film layers tuned to minimize scatter background, adaptable sample flow layers tuned to match crystal size, and a large sample area compatible with both raster scanning and rotation based serial data collection. The optically transparent chips can be used both for in situ protein crystallization (to eliminate crystal handling) or crystal slurry loading, with prepared samples stable for weeks in a humidified environment and for several hours in ambient conditions. Serial oscillation crystallography, using a multi-crystal rotational data collection approach, at a microfocus synchrotron beamline (SSRL, beamline 12-1) was used to benchmark the performance of the chips. Furthermore, high-resolution structures (1.3–2.7 Å) were collected from five different proteins – hen egg white lysozyme, thaumatin, bovine liver catalase, concanavalin-A (type VI), and SARS-CoV-2 nonstructural protein NSP5. Overall, our modular fabrication approach enables precise control over the cross-section of materials in the X-ray beam path and facilitates chip adaption to different sample and beamline requirements for user-friendly, straightforward diffraction measurements at room temperature.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2026

Low-energy 17 O(𝑛,𝛾)⁢ 18 O reaction within the microscopic potential model and its role for the weak 𝑟 process

The neutron radiative capture reaction 17 O ⁡(𝑛,𝛾) ⁢18 O plays a pivotal role in both nuclear structure studies and astrophysical nucleosynthesis, particularly in the formation of elements during hydrostatic and explosive stellar environments. We calculated the 17 O ⁡(𝑛,𝛾) ⁢18 O cross section within the Skyrme Hartree-Fock potential model and analyzed electric dipole 𝐸⁢1 transitions to both positive- and negative-parity states below the α-decay threshold in 18 O. Our cross sections are significantly different from the data available in commonly used libraries. We further investigate the impact of the new calculated cross section on weak 𝑟-process nucleosynthesis using large-scale reaction network calculations across a wide range of electron fractions and entropies. Our results show that the 17 O ⁡(𝑛,𝛾) ⁢18 O reaction rate significantly influences the production of first 𝑟-process peak elements, such as strontium, under specific astrophysical conditions. This study highlights the importance of accurate nuclear data for light isotopes in modeling heavy-element synthesis and provides updated reaction rates for future nucleosynthesis simulations.

6 ≤ A ≤ 19↗

DOE OSTI2020

Effective Missing Value Imputation Methods for Building Monitoring Data

To understand behaviors of natural and man-made events, such as energy consumption of buildings, which accounts for 40% of energy uses in the US, we deploy automated monitoring devices to record periodic observations. However, such experimental and observation data often contains problems and irregularities that have to be cleaned up before analyses. Due to various conditions affecting sensor operations, the communication channels, recording steps, or the recording media, the recorded data might have missing values, errors, or anomalous values. An effective way to clean up these problems is to replace these missing values, errors and anomalous values with expected values, a process generally known as imputation. In this work, we survey commonly used missing value imputation techniques and compare their performance on a set of building monitoring data. To compare the different types of sensor measurements with widely varying characteristics, we use normalized root mean squared error (NRMSE) as the key metric for the effectiveness of the imputation methods. We additionally consider periodicity and run time when considering comparing methods. Through extensive testing, we find that for small gap sizes, up to 8 consecutive missing values, linear interpolation performs the best; for larger gaps stretching up to 48 consecutive missing values, K-nearest neighbors provides the most accurate imputations; for even larger gaps, more computational intensive methods, such as matrix factorization, achieve the smallest NRMSE. Additionally, we observe that these computationally intensive algorithms not only provide accurate imputations for large gaps, but are also more robust across all types of sensors.

Cho, B↗

DOE OSTI2025

Comparing Outdoor to Indoor Performance for Bifacial Modules Affected by Polarization-Type Potential-Induced Degradation

Bifacial photovoltaic (PV) modules have the advantage of using light reflected off of the ground to contribute to power production. Predicting the energy gain is challenging and requires complex models to do so accurately. Often, module degradation over time is neglected in models for the sake of simplicity or is underestimated. Comparing outdoor and indoor current–voltage (I–V) performance for bifacial modules is more challenging than for monofacial modules, as there are additional variables to consider such as rear albedo non-uniformity, cell mismatch, and their effects on temperature. This challenge is compounded when heterogeneous degradation modes occur, such as polarization-type potential-induced degradation (PID-p). To examine the effects of PID-p on I–V predictions using an empirical data-driven approach, 16 bifacial PERC modules are installed outdoors on racks with different albedo conditions. A subset is exposed to high-voltage biases of −1500 V or +1500 V. Outdoor data are traced at irradiance ranges of 150–250 W/m 2 , 500–600 W/m 2 , and 900–1000 W/m 2 . These curves are corrected using control module temperature, wire resistivity, and module resistance measured indoors. We examine several methods to transform indoor I–V curves to accurately, and more simply than existing methods, approximate outdoor performance for bifacial modules without and with varying levels of PID-p degradation. This way, bifacial performance modeling can be more accessible and informed by fielded, degraded modules. Distributions of percent errors between indoor and outdoor performance parameters and Mean Absolute Percent Errors (MAPEs) are used to assess method quality. Results including low-irradiance data (150–250 W/m 2 ) are discussed but are filtered for quantifying method quality as these data introduce substantial errors. The method with the most optimal tradeoff between low MAPE and analysis simplicity involves measuring the front side of a module indoors at an irradiance equal to plane-of-array irradiance plus the product of module bifaciality and albedo irradiance. This method gives MAPE values of 1–6.5% for non-degraded and 1.6–5.9% for PID-p degraded module performance.

14 SOLAR ENERGY↗

DOE OSTI2023

Kernel-based global sensitivity analysis obtained from a single data set

Results from global sensitivity analysis (GSA) often guide the understanding of complicated input–output systems. Kernel-based GSA methods have recently been proposed for their capability of treating a broad scope of complex systems. In this paper, we develop a new set of kernel GSA tools when only a single set of input–output data is available. Three key advances are made: (1) A new numerical estimator is proposed that demonstrates an empirical improvement over previous procedures. (2) A computational method for generating inner statistical functions from a single data set is presented. (3) A theoretical extension is made to define conditional sensitivity indices, which reveal the degree that the inputs carry shared information about the output when inherent input–input correlations are present. Utilizing these conditional sensitivity indices, a decomposition is derived for the output uncertainty based on what is called the optimal learning sequence of the input variables, which remains consistent when correlations exist between the input variables. Further, while these advances cover a range of GSA subjects, a common single data set numerical solution is provided by a technique known as the conditional mean embedding of distributions. The new methodology is implemented on benchmark systems to demonstrate the provided insights.

42 ENGINEERING↗

DOE OSTI2026

A large-scale benchmarking of deterministic and stochastic derivative-free optimization algorithms

This presentation summarizes our work in the PrOMMiS project on benchmarking of data-driven optimization algorithms and their applications in self-driving laboratories. This work supports the broader project goal of accelerating the identification of promising separation methods and operating conditions for critical minerals separation processes. We present a systematic benchmarking study of 42 data-driven optimization algorithms on a broad collection of 502 test problems. The results identify BAM, GLCCLUSTER, and MULTIMIN as the most effective optimization solvers, with BAM showing the highest overall performance and solving more than 80% of the benchmark problems. The study also shows that no single solver consistently outperforms the others across all problem types, indicating that our future laboratory applications may benefit from using a small set of strong solvers rather than relying on a single method. The presentation also illustrates an in-silico chemical reactor case study showing that data-driven optimization methods can guide autonomous experimentation in a self-driving laboratory and identify optimal operating conditions within a small number of experiments. Overall, the results provide a basis for selecting efficient optimization methods and demonstrate the practical use of data-driven optimization in self-driving laboratory workflows.

36 MATERIALS SCIENCE↗

DOE OSTI2026

Algorithm-guided experimentation for autonomous AI systems in self-driving laboratories

36 MATERIALS SCIENCE↗

DOE OSTI2020

Physical drivers facilitating a toxigenic cyanobacterial bloom in a major Great Lakes tributary

The Maumee River is the primary source for nutrients fueling seasonal Microcystis -dominated blooms in western Lake Erie's open waters though such blooms in the river are infrequent. The river also serves as source water for multiple public water systems and a large food services facility in northwest Ohio. On 20 September 2017, an unprecedented bloom was reported in the Maumee River estuary within the Toledo metropolitan area, which triggered a recreational water advisory. Here we (1) explore physical drivers likely contributing to the bloom's occurrence, and (2) describe the toxin concentration and bacterioplankton taxonomic composition. A historical analysis using 10-years of seasonal river discharge, water level, and local wind data identified two instances when high-retention conditions occurred over ≥ 10 d in the Maumee River estuary: in 2016 and during the 2017 bloom. Observation by remote sensing imagery supported the advection of cyanobacterial cells into the estuary from the lake during 2017 and the lack of an estuary bloom in 2016 due to a weak cyanobacterial bloom in the lake. A rapid-response survey during the 2017 bloom determined levels of the cyanotoxins, specifically microcystins, in excess of recreational contact limits at sites within the lower 20 km of the river while amplicon sequencing found these sites were dominated by Microcystis . These results highlight the need to broaden our understanding of physical drivers of cyanobacterial blooms within the interface between riverine and lacustrine systems, particularly as such blooms are expected to become more prominent in response to a changing climate.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

Advanced wind turbine control development using field test analysis for generator overspeed mitigation

Abstract Turbulent and gusty wind conditions can cause generator overspeed peaks to exceed a threshold that then lead to wind turbine shutdowns, which then decrease the energy production of the wind turbines. We derive so‐called “gust measures” that predict when generator overspeed peaks may occur. These gust measures are then used to develop advanced controllers to mitigate generator overspeed peaks so that wind turbines can operate more robustly in difficult wind conditions without exceeding generator overspeed thresholds that would lead to turbine shutdown events. The advanced controllers are demonstrated in nonlinear aeroelastic simulations using the open‐source wind turbine simulation tool OpenFAST. To increase the realism of the simulations, they are run using field‐replicated wind conditions and a wind turbine model based on data from an experimental field campaign on a downscaled demonstrator of a novel extreme‐scale, two‐bladed, downwind rotor design.

17 WIND ENERGY↗

DOE OSTI2022

Application of the metabolic modeling pipeline in KBase to categorize reactions, predict essential genes, and predict pathways in an isolate genome

The DOE Systems Biology Knowledgebase (KBase) platform offers a range of powerful tools for the reconstruction, refinement, and analysis of genome-scale metabolic models built from microbial isolate genomes. In this chapter, we describe and demonstrate these tools in action with an analysis of isoprene production in the Bacillus subtilis DSM genome. Two different methods are applied to build initial metabolic models for the DSM genome, then the models are gapfilled in three different growth conditions. Next, flux balance analysis (FBA) and flux variability analysis (FVA) techniques are applied to both study the growth of these models in minimal media and classify reactions within each model based on essentiality and functionality. The models are applied with the FBA method to predict essential genes, which are then compared to an updated list of essential genes obtained for B. subtilis 168, a very similar strain to the DSM isolate. The models are also applied to simulate Biolog growth conditions, and these results are compared with Biolog data collected for B. subtilis 168. Finally, the DSM metabolic models are applied to explore the pathways and genes responsible for producing isoprene in this strain. These studies demonstrate the accuracy and utility of models generated from the KBase pipelines, as well as exploring the tools available for analyzing these models.

DOE knowledgebase↗

DOE OSTI2021

NO x adsorption with CO and C 2 H 4 on Pd/SSZ-13: Experiments and modeling

A transient monolith model containing microkinetic schemes for NO uptake and release over Pd/SSZ-13 without and with reductants CO or C2H4 is presented. The scheme involves three cationic Pd species (Z – [PdOH] + , Z – Pd 2+ Z – , Z – Pd + ) as the active sites during uptake of NO, CO and C 2 H 4 , and their desorption and conversion at higher temperature. Kinetic parameters are estimated through a combination of density functional theory (DFT) estimates and a fit of uptake, desorption and conversion data. The tuned model is validated at different uptake temperatures, ramp rates, and flowrates. A “degree of uptake control” parameter is defined that helps to identify the step(s) that are uptake controlling. The model helps to interpret the data features and is used to identify operating conditions to meet application-relevant performance metrics, including NO trapping efficiency and NO release temperature. Furthermore, the model demonstrates higher NO uptake on Pd/SSZ-13 compared to Pd/ZSM-5.

42 ENGINEERING↗

DOE OSTI2024

Exploring the potential and impact of single-crystal active materials on dry-processed electrodes for high-performance lithium-ion batteries

Roll-to-roll powder-to-film dry processing (DP) and single-crystal (SC) active materials (AMs) with many advantages are two hot topics of lithium-ion batteries (LIBs). However, DP of SC AMs for LIBs is rarely reported. Consequently, the impact of SC AMs on dry-processed LIBs is not well understood. Herein, for the first time, via a set of experimental and theoretical studies of the conventional polycrystalline-AM- and SC-AM-based DPed electrodes (DPEs), this work not only reports a high-performance dry SC-AM cathode for LIB manufacturing, but also establishes some fundamental understanding of SC-based dry-processed electrodes, including their morphology, structure, mechanical strength, electronic conductivity and LIB electrochemical behavior. Further, the results suggest that DP of SC AMs is promising, which can dramatically improve the electrochemical kinetics at electrode level and particle level. Specifically, for the rate capability and long-term cyclability in full cells, SC DPEs exhibit a discharge specific capacity of 152.1 mAh g -1 at 1C and a capacity retention rate of 79.9 % at C/3 over 500 cycles, which are superior to those of PC DPEs (135.6 mAh g -1 and 68.3 %) at the same conditions and are further confirmed by the simulation data from the theoretical modelling study. Therefore, this comprehensive work marks a significant milestone for DP strategy and SC AMs, enlightening future research and development of LIB manufacturing.

25 ENERGY STORAGE↗

DOE OSTI2022

Multi-stage heat release of multi-component fuels: Insights and implications for advanced engine operation

Multi-stage heat release (MSHR) is a unique phenomenon typically seen in lean/diluted and low- to intermediate-temperature combustion. Despite its relevance to advanced engine operation, the MSHR of multi-component fuels has barely been quantified. This study aims to characterize the MSHR of multi-component fuels in a rapid compression machine (RCM) at conditions representative of advanced combustion engines. New experimental data are first reported in the RCM at an equivalence ratio of 0.4, pressure of 60 bar and temperatures from 702 to 795 K for a research grade, multi-component gasoline surrogate, termed PACE-20. Here, experiments confirm the existence of MSHR for PACE-20 at all temperatures, and reveal the strong inhibiting effect of temperature on MSHR, where increasing temperature inhibits MSHR by suppressing first stage heat release and promoting second and third stage heat release. A detailed chemical kinetic model is also adopted to model the experiments, with good agreement observed. The response of MSHR characteristics to changes in different engine operating parameters (i.e., temperature, pressure, equivalence ratio and CO 2 dilution level) and fuel compositions (i.e., mole fraction of n-pentane, n-heptane, isooctane, cyclopentane, 1-hexene, ethanol and aromatics in PACE-20) is further evaluated via statistical analysis coupling quasi-random sampling, extensive computer experiments and high-dimensional model representation. The change in MSHR is characterized through 8 quantities of interest (QOIs), i.e., duration and extent of each heat release stage, and their standard deviations. The analysis highlights the dependence of the QOIs on the individual parameters as well as their interplays, with temperature exhibiting the strongest impact among all parameters. Furthermore, it is demonstrated that if MSHR is to be facilitated and combustion phasing is to be extended to enable engine operation at higher compression ratios, it is recommended to operate the engine at low intake temperatures, boosted intake pressures and high CO 2 dilutions (i.e., high EGR levels) with gasoline fuels containing more n-pentane and less aromatics.

42 ENGINEERING↗