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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

AEOLUS: Advances in Experimental Design, Optimal Control, and Learning for Uncertain Complex Systems

The AEOLUS Center is dedicated to developing a unified optimization-under-uncertainty framework for (1) learning predictive models from data and (2) optimizing experiments, processes, and designs governed by these models, all driven by complex, uncertain energy systems. AEOLUS addressed the critical need for principled, rigorous, scalable, and structure-exploiting capabilities for exploring parameter and decision spaces of complex forward simulation models---the so-called outer loop. This report summarizes the work done under DE-SC0021077 on (1) nonlocal models for solidification problems, (2) a multifidelity method for a nonlocal diffusion model, and (3) multifidelity Monte Carlo methods.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Efficient and Selective Chemical Transformations in Highly Charged and Confined Nanodroplets

The acceleration of chemical reaction rates and the increased product selectivity in microdroplets compared to that in bulk solutions has become a topic of increasing interest that has been extensively characterized by electrospray ionization mass spectrometry (ESI-MS). However, the sources of this acceleration and the detailed relationships between droplet properties and resulting reaction rate acceleration are still under debate. Moreover, droplet properties are governed by multiple interrelated experimental parameters, i.e., electrospray voltage, solution flow rate, etc., which makes it difficult and time-consuming to explore this diverse parameter space using traditional manual experimental or computational approaches. In this work, we developed an automated experimental platform integrating reactions in controlled charged microdroplet environments with ESI-MS characterization and sequential hybrid Bayesian modeling, as well as an optimal experimental design framework, to achieve multidimensional parameter optimization for higher reaction turnover rates, based on a model reaction of tetraethylenepentamine (TEPA) with carbon dioxide. With the current platform, we have achieved automated scans with a range of electrospray voltages and solution flow rates, and determined and optimized parameter settings to achieve increased reaction turnovers. We have also linked this platform to the underlying properties of droplets via a hybrid model incorporating physics, high-level theoretical calculations, and machine learning (ML) approaches. The autonomous platform is broadly applicable to a range of chemical reactions relevant to DOE’s mission in chemical separations, catalysis, and materials synthesis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Advances in Experimental Design, Optimal Control, and Learning for Uncertain Complex Systems (Final Report for AEOLUS)

The AEOLUS Center is dedicated to developing a unified optimization-under-uncertainty framework for: (1) learning predictive models from data; and (2) optimizing experiments, processes, and designs governed by these models, all driven by complex, uncertain energy systems. AEOLUS addresses the critical need for principled, rigorous, scalable, and structure-exploiting capabilities for exploring parameter and decision spaces of complex forward simulation models. This report summarizes the key highlights of our research during the period of performance.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Preparing large area of thermochromic nanocomposite films for smart window application

The challenges posed by high building energy bill necessitate proactive implementation of energy-efficient strategies to minimize the uses of both electricity and heating gas and promote the use of free solar energy sources. "Smart" window films/glass, leveraging thermochromic vanadium dioxide (VO 2 ), offer an adaptive approach to harness solar energy, significantly reducing the thermal load of buildings. This is achieved by reflecting the infrared portion of the solar spectrum through a phase transition in the monoclinic M-phase (M) of VO 2 induced by heating. In this study, a scalable continuous flow synthesis process invented by Argonne was employed to explore a wide parameter space, targeting the controlled synthesis of monoclinic VO 2 (M) nanoparticles with well-controlled sizes and morphologies. Additionally, a doping strategy and surface modifications were utilized for high-throughput tuning of the metal-to-insulator transition temperature. Strategies to enhance the solar modulation properties of VO 2 nanoparticles in polymer films were investigated through (i) morphological transformation from spherical to nanorod structures, (ii) surface modification with low refractive index ligands, and (iii) incorporation of additional thermochromic materials for improved solar energy modulation and aesthetically appealing colors in smart films. Furthermore, the study outlines scaling-up synthesis methods for VO 2 nanoparticles in a continuous flow reactor using hydrazine monohydrate. This comprehensive investigation provides valuable insights into the upscaling synthesis and design of advanced VO 2 /polymer composite smart window films with enhanced functionality, solar energy modulation and visible light transmittance, driving the technology a step close for industrial application.

14 SOLAR ENERGY↗
DOE OSTI2022

Parametric model-order-reduction development for unsteady convection

A time-averaged error indicator with POD- h Greedy is developed to drive parametric model order reduction (pMOR) for 2D unsteady natural convection in a high-aspect ratio slot parameterized with the Prandtl number, Rayleigh number, and slot angle with respect to the gravity. The error indicator is extended to accommodate the energy equation and Leray regularization. Despite being two-dimensional and laminar, the target flow regime presents several challenges: 1) there is a bifurcation in the angle parameter space; 2) the solution can be multivalued, even at steady state; and 3) the solution exhibits spatio-temporal chaos at several points in the parameter space. The authors explore several reduced-order models (ROMs) and demonstrate that Leray-regularized Galerkin ROMs provide a robust solution approach for this class of flows. They further demonstrate that error-indicated pMOR can efficiently predict several QOIs, such as mean flow, mean Nusselt number and mean turbulent kinetic energy, even in the presence of a bifurcation. Finally, they show that spatio-temporal chaos can lead to lack of reproducibility in both the full-order model and the reduced-order model and that the variance in the full-order model provides a lower bound on the pMOR error in these cases.

leray regularization↗
DOE OSTI2025

Numerical simulations of liquid jetting with solid inclusions

The dynamics of finite-sized particles in fluids, and their influence on the overall flow, are of great interest across several industrial, environmental, and medical fields. In the context of inkjet printing, the presence of solid inclusions can be either intentional, as in additive manufacturing, or unintentional, as in standard printing processes. These inclusions can strongly impact the jetting process, causing effects such as jet asymmetry, bubble entrapment, and the formation of satellite droplets. Understanding and controlling particle behavior is therefore essential, particularly to predict how and when particles are ejected over multiple jetting cycles. It is therefore critical to develop reliable models that allow for a deeper understanding of the complex interplay between particle and fluid during the whole printing process. To address this, we present a tailored implementation of the Color-Gradient multicomponent Lattice Boltzmann Method for fully resolved three-dimensional (3D) simulations of multicycle liquid jetting with particles. Our method supports realistic parameter settings aligned with industrial inkjet systems, and we provide both qualitative and quantitative validation against experimental data. Additionally, we introduce a simplified model based on the Stokes drag law, in which solid particles are represented as point particles and do not influence the fluid flow. Despite this limitation, the model offers a computationally efficient means to explore the vast parameter space typically encountered in industrial applications, allowing, e.g., identifying critical ejection regions and estimating the number of cycles required for particle release. These qualitative insights are valuable for guiding and complement fully two-way coupled simulations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

The Physics of Core-Collapse Supernovae: Explosion Mechanism and Explosive Nucleosynthesis

Recent developments in multi-dimensional simulations of core-collapse supernovae have considerably improved our understanding of this complex phenomenon. In addition to that, one-dimensional (1D) studies have been employed to study the explosion mechanism and its causal connection to the pre-collapse structure of the star, as well as to explore the vast parameter space of supernovae. Nonetheless, many uncertainties still affect the late stages of the evolution of massive stars, their collapse, and the subsequent shock propagation. In this review, we will briefly summarize the state-of-the-art of both 1D and 3D simulations and how they can be employed to study the evolution of massive stars, supernova explosions, and shock propagation, focusing on the uncertainties that affect each of these phases. Finally, we will illustrate the typical nucleosynthesis products that emerge from the explosion.

Boccioli, Luca (ORCID:000000024819310X)↗
DOE OSTI2021

Analysis of Previously Classified White Dwarf–Main-sequence Binaries Using Data from the APOGEE Survey

We present analyses of near-infrared spectroscopic data from the Apache Point Observatory Galactic Evolution Experiment (APOGEE) survey for 45 previously confirmed or candidate white dwarf–main-sequence (WDMS) binaries identified by the optical Sloan Digital Sky Survey (SDSS) and LAMOST surveys. Among these 45 systems, we classify three as having red giant primaries in the LAMOST sample and 14 as young stellar object contaminants in the photometrically identified SDSS sample. From among the subsample of 28 systems that we confirm to have MS primaries, we derive and place limits on orbital periods and velocity amplitudes for 14. Seven systems have significant velocity variations that warrant a post-common-envelope (PCE) binary classification, four of which are newly classified, three of which are newly confirmed, and five for which we can derive full orbital parameters. If confirmed, one of these newly discovered systems (2M14544500+4626456) will have the second-longest orbital period reported for a typical compact PCE WDMS binary (P = 15.1 days). In addition to the seven above, we also recover and characterize with APOGEE data the well-known PCE WDMS systems EG UMa and HZ 9. We also investigate the overall metallicity distribution of the WDMS sample, which is a parameter space not often explored for these systems. Of note, we find one system (2M14244053+4929580) to be extremely metal-poor ([Fe/H] = − 1.42) relative to the rest of the near-solar sample. Additionally, the PCE systems in our sample are found to be, on average, higher in metallicity than their wide-binary counterparts, though we caution that with this small number of systems, the sample may not be representative of the overall distribution of WDMS systems.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Markov Chain Monte Carlo Predictions of Neutron-rich Lanthanide Properties as a Probe of r -process Dynamics

Lanthanide element signatures are key to understanding many astrophysical observables, from merger kilonova light curves to stellar and solar abundances. To learn about the lanthanide element synthesis that enriched our solar system, we apply the statistical method of Markov Chain Monte Carlo to examine the nuclear masses capable of forming the r-process rare-earth abundance peak. We describe the physical constraints we implement with this statistical approach and demonstrate the use of the parallel chains method to explore the multidimensional parameter space. We apply our procedure to three moderately neutron-rich astrophysical outflows with distinct types of r-process dynamics. We show that the mass solutions found are dependent on outflow conditions and are related to the r-process path. We describe in detail the mechanism behind peak formation in each case. We then compare our mass predictions for neutron-rich neodymium and samarium isotopes to the latest experimental data from the CPT at CARIBU. Here, we find our mass predictions given outflows that undergo an extended (n,γ)⇌(γ,n) equilibrium to be those most compatible with both observational solar abundances and neutron-rich mass measurements.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Search for Solar Flare Neutrinos with the KamLAND Detector

Abstract We report the result of a search for neutrinos in coincidence with solar flares from the GOES flare database. The search was performed on a 10.8 kton-year exposure of KamLAND collected from 2002 to 2019. This large exposure allows us to explore previously unconstrained parameter space for solar flare neutrinos. We found no statistical excess of neutrinos and established 90% confidence level upper limits of 8.4 × 10 7 cm −2 (3.0 × 10 9 cm −2 ) on the electron antineutrino (electron neutrino) fluence at 20 MeV normalized to the X12 flare, assuming that the neutrino fluence is proportional to the X-ray intensity.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Evolution of Main-sequence-like Surviving Companions in Type Ia Supernova Remnants

Recent theoretical and numerical studies of Type Ia supernovae (SNe Ia) explosions within the single-degenerate scenario suggest that the nondegenerate companions could survive during the supernova impact and could be detectable in nearby supernova remnants. However, observational efforts show less promising evidence of the existence of surviving companions from the standard single-degenerate channels. The spin-up/spin-down models are possible mechanisms to explain the nondetection of surviving companions. In these models, the spin-up phase could increase the critical mass for explosion, leading to a super-Chandrasekhar-mass explosion, and the spin-down phase could lead to extra mass loss and angular momentum redistribution. Since the spin-down timescale for the delayed explosion of a rotating white dwarf is unclear, in this paper we explore a vast parameter space of main-sequence-like surviving companions via two-dimensional hydrodynamic simulations of supernova impact and the subsequent stellar evolution of surviving companions. Tight universal relations to describe the mass-stripping effect, supernova kick, and depth of supernova heating are provided. Our results suggest that the not-yet-detected surviving companions from observations of nearby SN Ia remnants might favor low-mass companions, short binary separation, or stronger supernova explosion energies than the standard single-degenerate channels.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

BICEP/Keck. XVI. Characterizing Dust Polarization through Correlations with Neutral Hydrogen

We characterize Galactic dust filaments by correlating BICEP/Keck and Planck data with polarization templates based on neutral hydrogen (H i) observations. Dust polarization is important for both our understanding of astrophysical processes in the interstellar medium (ISM) and the search for primordial gravitational waves in the cosmic microwave background (CMB). In the diffuse ISM, H i is strongly correlated with the dust and partly organized into filaments that are aligned with the local magnetic field. We analyze the deep BICEP/Keck data at 95, 150, and 220 GHz, over the low-column-density region of sky where BICEP/Keck has set the best limits on primordial gravitational waves. We separate the H i emission into distinct velocity components and detect dust polarization correlated with the local Galactic H i but not with the H i associated with Magellanic Stream i. We present a robust, multifrequency detection of polarized dust emission correlated with the filamentary H i morphology template down to 95 GHz. For assessing its utility for foreground cleaning, we report that the Hi morphology template correlates in B modes at a ~10%–65% level over the multipole range 20 < ℓ < 200 with the BICEP/Keck maps, which contain contributions from dust, CMB, and noise components. We measure the spectral index of the filamentary dust component spectral energy distribution to be β = 1.54 ± 0.13. We find no evidence for decorrelation in this region between the filaments and the rest of the dust field or from the inclusion of dust associated with the intermediate velocity H i. Finally, we explore the morphological parameter space in the H i-based filamentary model.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

From Feast to Famine: A Systematic Study of Accretion onto Oblique Pulsars with 3D GRMHD Simulations

Disk-fed accretion onto neutron stars can power a wide range of astrophysical sources ranging from X-ray binaries, to accretion-powered millisecond pulsars, ultraluminous X-ray sources, and gamma-ray bursts. A crucial parameter controlling the gas–magnetosphere interaction is the strength of the stellar dipole. In addition, coherent X-ray pulsations in many neutron star systems indicate that the star's dipole moment is oblique relative to its rotation axis. Therefore, it is critical to systematically explore the 2D parameter space of the star's magnetic field strength and obliquity, which is what this work does, for the first time, in the framework of 3D general-relativistic magnetohydrodynamics. If the accretion disk carries its own vertical magnetic field, this introduces an additional factor: the relative polarity of the disk and stellar magnetic fields. We find that depending on the strength of the stellar dipole and the star–disk relative polarity, the neutron star's jet power can either increase or decrease with increasing obliquity. For weak dipole strength (equivalently, high accretion rate), the parallel polarity results in a positive correlation between jet power and obliquity, whereas the antiparallel orientation displays the opposite trend. For stronger dipoles, the relative-polarity effect disappears, and jet power always decreases with increasing obliquity. The influence of the relative polarity gradually disappears as obliquity increases. Highly oblique pulsars tend to have an increased magnetospheric radius, a lower mass accretion rate, and enter the propeller regime at lower magnetic moments than aligned stars.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Physics makes the difference: Bayesian optimization and active learning via augmented Gaussian process

Abstract Both experimental and computational methods for the exploration of structure, functionality, and properties of materials often necessitate the search across broad parameter spaces to discover optimal experimental conditions and regions of interest in the image space or parameter space of computational models. The direct grid search of the parameter space tends to be extremely time-consuming, leading to the development of strategies balancing exploration of unknown parameter spaces and exploitation towards required performance metrics. However, classical Bayesian optimization (BO) strategies based on the Gaussian process (GP) do not readily allow for the incorporation of the known physical behaviors or past knowledge. Here we explore a hybrid optimization/exploration algorithm created by augmenting the standard GP with a structured probabilistic model of the expected system’s behavior. This approach balances the flexibility of the non-parametric GP approach with a rigid structure of physical knowledge encoded into the parametric model. The fully Bayesian treatment of the latter allows additional control over the optimization via the selection of priors for the model parameters. The method is demonstrated for a noisy version of a standard univariate test function used to evaluate optimization algorithms and further extended to physical lattice models. This methodology is expected to be universally suitable for injecting prior knowledge in the form of physical models and past data in the BO framework.

42 ENGINEERING↗
DOE OSTI2025

Computationally efficient subglacial drainage modelling using Gaussian process emulators: GlaDS-GP v1.0

Subglacial drainage models represent water flow at the ice–bed interface through coupled distributed and channelized systems to determine water pressure, discharge, and drainage system geometry. While they are used to understand processes such as the relationship between surface melt and ice flow, the number of uncertain model parameters and the computational cost of running models makes it difficult to adequately explore the high-dimensional parameter space and evaluate uncertainty in model predictions. Here, we develop Gaussian process (GP) emulators that make fast predictions with associated uncertainty of subglacial drainage model outputs. Using a truncated principal component (PC) basis representation, we construct a GP emulator for diurnally averaged subglacial water pressure. We also explore emulation of scalar variables describing drainage efficiency and configuration. We train the emulators using ensembles of up to 512 simulations varying eight parameters of the Glacier Drainage System (GlaDS) model on a synthetic domain intended to represent an ice-sheet margin. The emulators make predictions ∼ 1000 times faster than GlaDS simulations, with errors <3 % for the water pressure field and ∼ 5 %–9 % for drainage efficiency and configuration. We apply the emulators to explore the eight-dimensional parameter space by computing variance-based parameter sensitivity indices, finding that three parameters (ice flow coefficient, bed bump aspect ratio, and the subglacial cavity system conductivity) explain 90 % of the variance in modelled water pressure in response to parameter changes. The GP emulator approach described here is well suited to integrating observational data with models to make calibrated, credible predictions of subglacial drainage.

58 GEOSCIENCES↗
DOE OSTI2021

Navigator function for the conformal bootstrap

Current numerical conformal bootstrap techniques carve out islands in theory space by repeatedly checking whether points are allowed or excluded. We propose a new method for searching theory space that replaces the binary information "allowed"/"excluded" with a continuous "navigator" function that is negative in the allowed region and positive in the excluded region. Such a navigator function allows one to efficiently explore high-dimensional parameter spaces and smoothly sail towards any islands they may contain. The specific functions we introduce have several attractive features: they are well-defined in large regions of parameter space, can be computed with standard methods, and evaluation of their gradient is immediate due to an SDP gradient formula that we provide. The latter property allows for the use of efficient quasi-Newton optimization methods, which we illustrate by navigating towards the 3d Ising island.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Toward Accelerating Discovery via Physics-Driven and Interactive Multifidelity Bayesian Optimization

Both computational and experimental material discovery bring forth the challenge of exploring multidimensional and often nondifferentiable parameter spaces, such as phase diagrams of Hamiltonians with multiple interactions, composition spaces of combinatorial libraries, processing spaces, and molecular embedding spaces. Often these systems are expensive or time consuming to evaluate a single instance, and hence classical approaches based on exhaustive grid or random search are too data intensive. This resulted in strong interest toward active learning methods such as Bayesian optimization (BO) where the adaptive exploration occurs based on human learning (discovery) objective. However, classical BO is based on a predefined optimization target, and policies balancing exploration and exploitation are purely data driven. In practical settings, the domain expert can pose prior knowledge of the system in the form of partially known physics laws and exploration policies often vary during the experiment. Here, we propose an interactive workflow building on multifidelity BO (MFBO), starting with classical (data-driven) MFBO, then expand to a proposed structured (physics-driven) structured MFBO (sMFBO), and finally extend it to allow human-in-the-loop interactive interactive MFBO (iMFBO) workflows for adaptive and domain expert aligned exploration. These approaches are demonstrated over highly nonsmooth multifidelity simulation data generated from an Ising model, considering spin–spin interaction as parameter space, lattice sizes as fidelity spaces, and the objective as maximizing heat capacity. Detailed analysis and comparison show the impact of physics knowledge injection and real-time human decisions for improved exploration with increased alignment to ground truth. Here, the associated notebooks allow to reproduce the reported analyses and apply them to other systems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Kitaev model on a quantum computer using VQE with Majorana fermions

We study the simulation of the Kitaev spin model on quantum computers. In particular we focus on the models defined on the honeycomb, and square-octagon lattices. Using a fermionic language to describe these models reveals a region of the parameter space that is exactly solvable. We explore an ansatz that is capable of expressing the ground state in the exactly solvable region of the parameter space and extend it outside this region with good accuracy. Doing the calculation using fermions, while requiring the introduction of a non-local map from the fermionic Hilbert space to that of qubits, offers the potentially interesting application of realizing non-abelian anyons on quantum computers, and can also lead to a reduction in the number of qubits required by half.

Jahin, Ammar↗