Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “parameter space exploration”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Gravitational waves and galaxies cross-correlations: a forecast on GW biases for future detectors

ABSTRACT Gravitational waves (GWs) have rapidly become important cosmological probes since their first detection in 2015. As the number of detected events continues to rise, upcoming instruments like Einstein Telescope (ET) and Cosmic Explorer (CE) will observe millions of compact binary mergers. These detections, coupled with galaxy surveys by instruments such as the Dark Spectroscopic Energy Instrument (DESI), Euclid, and the Vera Rubin Observatory, will provide unique information on the large-scale structure of the universe by cross-correlating GWs with the distribution of galaxies hosting them. In this paper, we focus on how cross-correlations constrain the clustering bias of GWs emitted by the coalescence of binary black holes (BBHs). This parameter links BBHs to the underlying dark matter distribution, hence informing us how they populate galaxies. Using a multitracer approach, we forecast the precision of these measurements under different survey combinations. Our results indicate that current GW detectors will have limited precision, with measurement errors as high as $\displaystyle \sim 50~{{\ \rm per\ cent}}$. However, third-generation detectors like ET, when cross-correlated with Legacy Survey of Space and Time (LSST) data, can improve clustering bias measurements to within 2.5 per cent. Furthermore, we demonstrate that these cross-correlations can enable a per cent-level measurement of the magnification lensing effect on GWs. Despite this, there is a degeneracy between magnification and evolution biases, which hinders the precision of both. This degeneracy is most effectively addressed by assuming knowledge of one bias or targeting an optimal redshift range of $\displaystyle 1 \lt z \lt 2.5$. Our analysis opens new avenues for studying the distribution of BBHs and testing the nature of gravity through large-scale structure.

Zazzera, Stefano (ORCID:0000000158979221)↗
DOE OSTI2024

Chapter 4: Physically informed deep learning networks for simulating microstructure evolution of 3D polycrystals

As discussed in the previous chapter, high energy diffraction microscopy (HEDM) is used to study the micromechanical evolution of a material during in situ loading. HEDM experiments have been used to verify crystal plasticity (CP) simulations [119, 91, 90, 120], for experimental planning, material design, and to further analyze experimental results. However, Fast Fourier transform-based CP (CP-FFT) or finite element-based CP (CP-FE) methods are often too slow to be used in real-time during an experiment. CP-FFT is faster than CP-FE simulations due to the absence of meshing, but can still take hours to simulate the response of a single volume depending on the size and number of strain steps [127]. Reducing computation time would create a larger exploration space in planning and design, and enable faster analysis of experimental results and real-time feedback during an experiment. This research expands upon previous works to develop a workflow for predicting the full-field evolution of a 3D polycrystal. The workflow is simplified from previous works to predict only orientation and elastic strain tensors (from which stress tensors are calculated). The network is physically informed through loss functions and network architecture for a more robust model. The orientation predictions are informed about the cubic crystal symmetry of the material by incorporating disorientation and misorientation information into the network architecture and loss. The Von Mises stress is used to enforce the correct stress-strain trends in the strain tensor predictions. Additional total strain steps from the elastic and elastoplastic region are included to better capture the stress-strain evolution at smaller total strain steps. Material and hardening parameters are additional inputs into the networks to further inform the network and to study the network’s ability to predict different materials other than those used for training.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Fully Distributed Multi-Material Magnetic Sensing Structures for Multiparameter DAS Applications

This dissertation demonstrates the first of its kind distributed magnetic field sensor based on a fiber optic distributed acoustic sensing (DAS) scheme. Ferromagnetic nickel and Metglas® were dispersed internally within a fiber optic preform and then drawn on an in-house fiber optic draw tower to lengths in the kilometers. Due to the close proximity of the ferromagnetic metals and fiber optic core, the magnetostrictive strain response of the ferromagnetic materials when exposed to a magnetic field would perturbate within the fiber cladding and transfer that strain, internally, to the fiber optic core. Strain resulting from the magnetostrictive effect allows the DAS based sensor to accurately translate strain into readable magnetic field data. Due to the high sensitivity seen in this sensor design, multiparameter sources, acoustic and magnetic fields, were tested and validated and a three dimensional magnetic-field vector sensor was proposed. Numerical analysis of the novel sensor design was first implemented using COMSOL Multiphysics, where inputs such as magnetostrictive element shape, size, distance, and number were first investigated. Upon optimizing system constraints, the sensor design was further modified such that single mode operation was consistent across multiple fiber draws while retaining high strain transfer from the ferromagnetic elements to the fiber optic core. Ferromagnetic material selection was evaluated as a function of the saturation magnetostriction constants and a total of 4 modules were used to fully characterize the complex physics involved in this sensor design. All fabrication and testing were performed in-house using a full scale 3-story fiber draw tower and custom environmental testing stations to imitate naturally occurring events such as magnetic or acoustic point sources. A unique stacking method was used to embed ferromagnetic nickel and Metglas® into a fiber optic preform which when combined with a custom fiber draw process resulted in consistent multi-material fibers drawn to lengths of 1-km. In-house testing facilities included different types of electromagnetic generators, in addition to a soil test bed, and an outdoor test bed which allowed 100 meters of fiber to be tested simultaneously. All tested sensors demonstrated high strain transfer capabilities on the order of 0.01-10 μϵ depending on the materials used, ferromagnetic rod number, and core to metal spacing. Due to the sensitivity of the system the difference between AC and DC was distinct, and directional magnetostriction was studied. Transverse and longitudinal magnetic wave propagation was controlled through a solenoid and rectangular Helmholtz coil, both built in-house. A three-dimensional magnetic field vector sensor was proposed due to the success of the magnetic field sensor, and a design was proposed and initially tested to validate direction as a function of field strength and distance. To summarize, this dissertation explores the first fully distributed magnetic field sensor using DAS based techniques and one of the first multi-material fiber draw processes which can produce consistent single mode fiber up to 1-km. Due to extensive FEA modeling, multiple iterations of the magnetic sensor were fully characterized and an equation describing the relationship between sensor design and strain transfer has been created and validated experimentally. Multi-parameter tests including acoustic and magnetic fields were implemented and an algorithm was developed to separate the mixed signals. Finally, a test was performed to demonstrate the feasibility of sensing magnetic fields directionally. Cumulative results demonstrate a high-quality sensor alternative to current designs which may surpass other magnetic sensors due to innate multi-parameter capabilities, in addition to the inexpensive production cost and extremely long operating lengths.

02 PETROLEUM↗
DOE OSTI2024

Hijacking a rapid and scalable metagenomic method reveals subgenome dynamics and evolution in polyploid plants

Premise: The genomes of polyploid plants archive the evolutionary events leading to their present forms. However, plant polyploid genomes present numerous hurdles to the genome comparison algorithms for classification of polyploid types and exploring genome dynamics. Methods: Here, the problem of intra- and inter-genome comparison for examining polyploid genomes is reframed as a metagenomic problem, enabling the use of the rapid and scalable MinHashing approach. To determine how types of polyploidy are described by this metagenomic approach, plant genomes were examined from across the polyploid spectrum for both k-mer composition and frequency with a range of k-mer sizes. In this approach, no subgenome-specific k-mers are identified; rather, whole-chromosome k-mer subspaces were utilized. Results: Given chromosome-scale genome assemblies with sufficient subgenome-specific repetitive element content, literature-verified subgenomic and genomic evolutionary relationships were revealed, including distinguishing auto- from allopolyploidy and putative progenitor genome assignment. The sequences responsible were the rapidly evolving landscape of transposable elements. An investigation into the MinHashing parameters revealed that the downsampled k-mer space (genomic signatures) produced excellent approximations of sequence similarity. Furthermore, the clustering approach used for comparison of the genomic signatures is scrutinized to ensure applicability of the metagenomics-based method. Discussion: The easily implementable and highly computationally efficient MinHashing-based sequence comparison strategy enables comparative subgenomics and genomics for large and complex polyploid plant genomes. Such comparisons provide evidence for polyploidy-type subgenomic assignments. In cases where subgenome-specific repeat signal may not be adequate given a chromosomes' global k-mer profile, alternative methods that are more specific but more computationally complex outperform this approach.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Effect of the space weather conditions on the Earth magnetosphere soft X-ray emissivity

The aim of the study is to model and characterize the soft X-ray emissivity on the Earth magnetosphere for different space weather conditions (SWC), providing information to interpret the soft X-ray measurements of the Solar wind Magnetosphere Ionosphere Link Explorer space mission. The MHD code pluto in spherical coordinates is used to perform parametric studies with respect to the solar wind (SW) dynamic pressure (considering density and velocity effects independently) as well as the IMF intensity and orientation, predicting the soft X-ray emissivity for different SWC. The integrated soft X-ray emissivity inside the magnetosheath is calculated as a proxy of the soft X-ray emission dependencies with the SWC independently of the satellite orbit and camera line of sight. The analysis indicates fluctuations of the interplanetary magnetic field (IMF) orientation and magnitude may significantly affect the measured soft X-ray emission although changes in the SW dynamic pressure should be the main source of variability. The southward IMF orientation leads to the configuration with the largest soft X-ray emissivity and northward to the lowest. Strongly distorted magnetospheres explored in configurations showing SW and IMF parameters comparable to the impact of interplanetary coronal mass ejections may show a decrease of the soft X-ray emissivity as the IMF magnitude increases, explained by the strong magnetosphere compression and constriction of the magnetosheath region where the soft X-ray emissivity maximum is located. The simulations also indicate large excursions of the soft X-ray emissivity maximum inside the magnetosheath as the IMF magnitude and SW dynamic pressure fluctuate particularly for radial and ecliptic IMF orientations.

Earth↗
DOE OSTI2025

Dark baryon black holes

We explore a novel mechanism for dark matter production through the formation of light black holes from the collapse of dark baryons in confining SU(N) gauge theories in the large- N limit. While glueballs and mesons cannot form black holes under physically reasonable conditions, we prove that for appropriate ranges of the confinement scale, quark masses, number of colors N , and dark sector temperature, dark baryons can produce Planck-scale black hole relics in the early Universe. Assuming the relics are stable, the abundance of both the dark baryon black hole population directly arising at confinement and that frozen in from dark glueball and meson pair annihilation are exponentially suppressed in N , leading to an upper limit N ≲ 100 and of a few hundred Planck units in mass for models where the black hole relics are the entirety of the dark matter. We present a detailed numerical study of the parameter space where this scenario is realized. Published by the American Physical Society 2025

Profumo, Stefano (ORCID:0000000291597556)↗
DOE OSTI2024

Dark dimension and decaying dark matter gravitons

We explore the cosmology of the dark-dimension scenario taking into account perturbations in the linear regime. In the context of the dark-dimension scenario, a natural candidate for dark matter in our Universe is the excitations of a tower of massive spin-2 Kaluza-Klein (KK) gravitons. These dark gravitons are produced in the early Universe and decay to lighter KK gravitons during the course of cosmological evolution. The decay causes the average dark matter mass to decrease as the Universe evolves. In addition, the kinetic energy liberated in each decay leads to a kick velocity for the dark matter particles, leading to a suppression of structure formation. Using current cosmic microwave background (), baryon acoustic oscillation, and cosmic shear (KiDS-1000) data, we put a bound on the dark matter kick velocity today v today ≤ 2.2 × 10 − 4 c at 95% CL. This leads to rather specific regions of parameter space for the dark-dimension scenario. The combination of the experimental bounds from cosmology, astrophysics, and table-top experiments leads to the range l 5 ∼ 1 – 10 μ m for the size of the dark dimension. The dark-dimension scenario is found to be remarkably consistent with current observations and provides signatures that are within reach of near-future experiments. Published by the American Physical Society 2024

Astronomy & Astrophysics↗
DOE OSTI2021

Cosmology from clustering, cosmic shear, CMB lensing, and cross correlations: combining Rubin observatory and Simons Observatory

ABSTRACT In the near future, the overlap of the Rubin Observatory Legacy Survey of Space and Time (LSST) and the Simons Observatory (SO) will present an ideal opportunity for joint cosmological data set analyses. In this paper, we simulate the joint likelihood analysis of these two experiments using six two-point functions derived from galaxy position, galaxy shear, and CMB lensing convergence fields. Our analysis focuses on realistic noise and systematics models and we find that the dark energy Figure of merit (FoM) increases by 53 per cent (92 per cent) from LSST-only to LSST+SO in Year 1 (Year 6). We also investigate the benefits of using the same galaxy sample for both clustering and lensing analyses, and find the choice improves the overall signal to noise by ${\sim}30{-}40{{\ \rm per\ cent}}$, which significantly improves the photo-z calibration and mildly improves the cosmological constraints. Finally, we explore the effects of catastrophic photo-z outliers finding that they cause significant parameter biases when ignored. We develop a new mitigation approach termed ‘island model’, which corrects a large fraction of the biases with only a few parameters while preserving the constraining power.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Axion strings are superconducting

We explore the cosmological consequences of the superconductivity of QCD axion strings. Axion strings can support a sizeable chiral electric current and charge density, which alters their early universe dynamics. We examine the possibility that shrinking axion string loops can become effectively stable remnants called vortons, supported by the repulsive electromagnetic force of the string current. We find that vortons in our scenario are generically unstable, and so do not pose a cosmological difficulty. Furthermore, if a primordial magnetic field (PMF) exists in the early universe, a large current is induced on axion strings, creating a significant drag force from interactions with the surrounding plasma. As a result, the strings are slowed down, which leads to an orders of magnitude enhancement in the number of strings per Hubble volume. Finally, we study potential implications for the QCD axion relic abundance. The QCD axion window is shifted by orders of magnitude in some parts of our parameter space.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Using a Coarse-Grained Modeling Framework to Identify Oligomeric Motifs with Tunable Secondary Structure

Coarse-grained modeling can be used to explore general theories that are independent of specific chemical detail. In this paper, we present cg_openmm, a Python-based simulation framework for modeling coarse-grained hetero-oligomers and screening them for structural and thermodynamic characteristics of cooperative secondary structures. cg_openmm facilitates the building of coarse-grained topology and random starting configurations, setup of GPU-accelerated replica exchange molecular dynamics simulations with the OpenMM software package, and features a suite of postprocessing thermodynamic and structural analysis tools. In particular, native contact analysis, heat capacity calculations, and free energy of folding calculations are used to identify and characterize cooperative folding transitions and stable secondary structures. In this work, we demonstrate the capabilities of cg_openmm on a simple 1–1 Lennard-Jones coarse-grained model, in which each residue contains 1 backbone and 1 side-chain bead. By scanning both nonbonded and bonded force-field parameter spaces at the coarse-grained level, we identify and characterize sets of parameters which result in the formation of stable helices through cooperative folding transitions. Furthermore, we show that the geometries and stabilities of these helices can be tuned by manipulating the force-field parameters.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Taking aim at the wino-Higgsino plane with the LHC

In this work we explore multiple search strategies for Higgsinos and mixed Higgsino-wino states in the minimal supersymmetric Standard Model (MSSM) and project the results onto the ( μ , M 2 ) plane. Assuming associated production of Higgsino-like pairs with a W / Z boson, we develop a search in a channel characterized by a hadronically tagged vector boson accompanied by missing energy. We use as our template an ATLAS search for dark matter produced in association with a hadronically decaying vector boson, but upgrade the search by implementing a joint-likelihood analysis, binning the missing transverse energy distribution, which greatly improves the search sensitivity. For Higgsino-like states (more than 96% admixture) we find sensitivity to masses up to 550 GeV. For well-mixed Higgsino-wino states (70–30% Higgsino) we still find sensitivities above 300 GeV. Using this newly proposed search, we draw a phenomenological map of the wino-Higgsino parameter space, recasting several complementary searches for disappearing tracks, soft leptons, trileptons, and hadronic diboson events in order to predict LHC coverage of the ( μ , M 2 ) mass plane at integrated luminosities of up to 3 ab − 1 . Altogether, the full run of the HL-LHC can exclude much of the “natural” ( μ , M 2 < 500 GeV ) wino-Higgsino parameter space. Published by the American Physical Society 2024

Astronomy & Astrophysics↗
DOE OSTI2025

Anisotropic electron damping and energy gap in Bi 2 ⁢Sr 2 ⁢CaCu 2 ⁢O 8+𝛿

The many-body electron-electron interaction in cuprates causes broadening of the electronic bands in 𝒌 space, leading to a deviation from the standard Fermi liquid. While a 𝒌-dependent anisotropic electronic scattering (𝒌-DAES) has been assessed by photoemission, its fingerprint in 𝑸 space has been scarcely considered. Here, we explore the 𝑸-dependent electron dynamics in optimally doped Bi 2⁢ Sr 2 ⁢CaCu 2 ⁢O 8+𝛿 through the evolution of low-energy charge excitations as measured by resonant inelastic x-ray scattering (RIXS). In the normal state, the RIXS spectra display a continuum of excitations down to 0 meV, while the superconducting state features a spectral weight suppression below 80 meV without any enhancement at higher energies. To interpret the energy and 𝑸 evolution of our data, we introduce a phenomenological expression of the charge susceptibility by including the 𝒌-DAES. We show that only the charge susceptibility with 𝒌-DAES captures the RIXS data, highlighting the importance of 𝒌-DAES when describing the 𝑸 dependence of charge excitations from 0 to a few eV scale. Furthermore, we also find that the inclusion of 𝒌-DAES is essential when quantitative parameters such as the electronic energy gap are extracted from RIXS data.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Targeted Adaptive Design

Modern advanced manufacturing and advanced materials design often require searches of relatively high-dimensional process control parameter spaces for settings that result in optimal structure, property, and performance parameters. The mapping from the former to the latter must be determined from noisy experiments or from expensive simulations. Here, we abstract this problem to a mathematical framework in which an unknown function from a control space to a design space must be ascertained by means of expensive noisy measurements, which locate control settings generating desired design features within specified tolerances, with quantified uncertainty. We describe targeted adaptive design (TAD), a new algorithm that performs this sampling task efficiently. TAD creates a Gaussian process surrogate model of the unknown mapping at each iterative stage, proposing a new batch of control settings to sample experimentally and optimizing the updated expected log-predictive probability density of the target design. TAD either stops upon locating a solution with uncertainties that fit inside the tolerance box or uses a measure of expected future information to determine that the search space has been exhausted with no solution. TAD thus embodies the exploration-exploitation tension in a manner that recalls, but is essentially different from, Bayesian optimization and optimal experimental design.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

$t\overline{t}t\overline{t}$ signatures through the lens of color-octet scalars

We reinterpret two recent LHC searches for events containing four top quarks ($t\overline{t}t\overline{t}$) in the context of supersymmetric models with Dirac gauginos and color-octet scalars (sgluons). We explore whether sgluon contributions to the four-top production cross section σ ($pp$ → $t\overline{t}t\overline{t}$) can accommodate an excess of four-top events recently reported by the ATLAS collaboration. We also study constraints on these models from an ATLAS search for new phenomena with high jet multiplicity and significant missing transverse energy ($E^{miss}_{T}$) sensitive to signals with four top quarks. We find that these two analyses provide complementary constraints, with the jets + $E^{miss}_{T}$ search exceeding the four-top cross section measurement in sensitivity for sgluons heavier than about 800 GeV. We ultimately find that either a scalar or a pseudoscalar sgluon can currently fit the ATLAS excess in a range of reasonable benchmark scenarios, though a pseudoscalar in minimal Dirac gaugino models is ruled out. We finally offer sensitivity projections for these analyses at the HL-LHC, mapping the 5σ discovery potential in sgluon parameter space and computing exclusion limits at 95% CL in scenarios where no excess is found.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

A Predictor-Corrector Strategy for Adaptivity in Dynamical Low-Rank Approximations

Here, in this paper, we present a predictor-corrector strategy for constructing rank-adaptive, dynamical low-rank approximations (DLRAs) of matrix-valued ODE systems. The strategy is a compromise between (i) low-rank step-truncation approaches that alternately evolve and compress solutions and (ii) strict DLRA approaches that augment the low-rank manifold using subspaces generated locally in time by the DLRA integrator. The strategy is based on an analysis of the error between a forward temporal update into the ambient full-rank space, which is typically computed in a step-truncation approach before recompressing, and the standard DLRA update, which is forced to live in a low-rank manifold. We use this error, without requiring its full-rank representation, to correct the DLRA solution. A key ingredient for maintaining a low-rank representation of the error is a randomized SVD, which introduces some degree of stochastic variability into the implementation. The strategy is formulated and implemented in the context of discontinuous Galerkin spatial discretizations of PDEs and applied to several versions of DLRA methods found in the literature as well as a new variant. Numerical experiments comparing the predictor-corrector strategy to other methods demonstrate robustness to overcome shortcomings of step truncation or strict DLRA approaches: The former may require more memory than is strictly needed, while the latter may miss transients solution features that cannot be recovered. The effect of randomization, tolerances, and other implementation parameters is also explored.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

New avenues for |∆ B | = 2 processes beyond neutron-antineutron oscillations

We explore baryon-number-violating (|∆ B | = 2) processes beyond the well-known neutron-antineutron ($n - \bar{n}$) oscillations, focusing on the $Λ - \bar{Λ}$ system. The presence of a strange quark in the Λ baryon introduces a new set of six-quark operators roughly of the form (uds) 2 , which are different from the (udd) 2 operators responsible for oscillations. Using the Standard Model Effective Field Theory (SMEFT), we classify all dimension-9 operators that cause |∆ B | = 2 transitions and study their UV completions mediated by exotic scalar fields with trilinear interactions. We demonstrate that in these models, oscillations can occur at tree level, with $n - \bar{n}$ mixing potentially appearing at higher loop levels. We employ a chiral effective theory to constrain the effective mass mixing δm Λ , deriving bounds from current experimental limits on $n - \bar{n}$ oscillations and dinucleon decays such as pp → K + K + . These bounds indicate that $Λ - \bar{Λ}$ oscillations probe a complementary parameter space, sensitive to baryon-number violation at scales up to 10 2 − 10 3 TeV. We show that the existing indirect bounds make it challenging to provide a competitive bound on δm Λ at BESIII.

Baryon/Lepton Number Violation↗
DOE OSTI2023

Simulating organic aerosol in Delhi with WRF-Chem using the volatility-basis-set approach: exploring model uncertainty with a Gaussian process emulator

The nature and origin of organic aerosol in the atmosphere remain unclear. The gas–particle partitioning of semi-volatile organic compounds (SVOCs) that constitute primary organic aerosols (POAs) and the multigenerational chemical aging of SVOCs are particularly poorly understood. The volatility basis set (VBS) approach, implemented in air quality models such as WRF-Chem (Weather Research and Forecasting model with Chemistry), can be a useful tool to describe emissions of POA and its chemical evolution. However, the evaluation of model uncertainty and the optimal model parameterization may be expensive to probe using only WRF-Chem simulations. Gaussian process emulators, trained on simulations from relatively few WRF-Chem simulations, are capable of reproducing model results and estimating the sources of model uncertainty within a defined range of model parameters. In this study, a WRF-Chem VBS parameterization is proposed; we then generate a perturbed parameter ensemble of 111 model runs, perturbing 10 parameters of the WRF-Chem model relating to organic aerosol emissions and the VBS oxidation reactions. This allowed us to cover the model's uncertainty space and to compare outputs from each run to aerosol mass spectrometer observations of organic aerosol concentrations and O:C ratios measured in New Delhi, India. The simulations spanned the organic aerosol concentrations measured with the aerosol mass spectrometer (AMS). However, they also highlighted potential structural errors in the model that may be related to unsuitable diurnal cycles in the emissions and/or failure to adequately represent the dynamics of the planetary boundary layer. While the structural errors prevented us from clearly identifying an optimized VBS approach in WRF-Chem, we were able to apply the emulator in the following two periods: the full period (1–29 May) and a subperiod period of 14:00–16:00 h LT (local time) on 1–29 May. The combination of emulator analysis and model evaluation metrics allowed us to identify plausible parameter combinations for the analyzed periods. We demonstrate that the methodology presented in this study can be used to determine the model uncertainty and to identify the appropriate parameter combination for the VBS approach and hence to provide valuable information to improve our understanding of OA production.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

A Novel Active Optimization Approach for Rapid and Efficient Design Space Exploration Using Ensemble Machine Learning

In this work, a novel design optimization technique based on active learning, which involves dynamic exploration and exploitation of the design space of interest using an ensemble of machine learning algorithms, is presented. In this approach, a hybrid methodology incorporating an explorative weak learner (regularized basis function model) that fits high-level information about the response surface and an exploitative strong learner (based on committee machine) that fits finer details around promising regions identified by the weak learner is employed. For each design iteration, an aristocratic approach is used to select a set of nominees, where points that meet a threshold merit value as predicted by the weak learner are selected for evaluation. In addition to these points, the global optimum as predicted by the strong learner is also evaluated to enable rapid convergence to the actual global optimum once the most promising region has been identified by the optimizer. Additionally, this methodology is first tested by applying it to the optimization of a two-dimensional multi-modal surface and, subsequently, to a complex internal combustion (IC) engine combustion optimization case with nine control parameters related to fuel injection, initial thermodynamic conditions, and in-cylinder flow. It is found that the new approach significantly lowers the number of function evaluations that are needed to reach the optimum design configuration (by up to 80%) when compared to conventional optimization techniques, such as particle swarm and genetic algorithm-based optimization techniques.

24 POWER TRANSMISSION AND DISTRIBUTION↗