SEARCH · Engineering Papers

Results for “Response surface model”

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,297 records · Page 72

All sources NASA NTRS DOE OSTI

NASA NTRS2016

Strata-1: An International Space Station Experiment into Fundamental Regolith Processes in Microgravity

The Strata-1 experiment will study the evolution of asteroidal regolith through long-duration exposure of simulant materials to the microgravity environment on the International Space Station (ISS). Many asteroids feature low bulk densities, which implies high values of porosity and a mechanical structure composed of loosely bound particles, (i.e. the "rubble pile" model), a prime example of a granular medium. Even the higher-density, mechanically coherent asteroids feature a significant surface layer of loose regolith. These bodies are subjected to a variety of forces and will evolve in response to very small perturbations such as micrometeoroid impacts, planetary flybys, and the YORP effect. Our understanding of this dynamical evolution and the inter-particle forces involved would benefit from long-term observations of granular materials exposed to small vibrations in microgravity. A detailed understanding of asteroid mechanical evolution is needed in order to predict the surface characteristics of as-of-yet unvisited bodies, to understand the larger context of samples collected by missions such as OSIRIS-REx and Hayabusa 1 and 2, and to mitigate risks for both manned and unmanned missions to asteroidal bodies. Understanding regolith dynamics will inform designs of how to land and set anchors, safely sample/move material on asteroidal surfaces, process large volumes of material for in situ resource utilization (ISRU) purposes, and, in general, predict behavior of large and small particles on disturbed asteroid surfaces.

Fries, M.↗

NASA NTRS1993

Similarity and scale in catchment storm response

Until recently, very little progress had been made in understanding the relationship between small-scale variability of topography, soil, and rainfalls and the storm response seen at the catchment scale. The work reviewed here represents the first attempt at a systematic theoretical framework for such understanding in the context of surface runoff generation by different processes. The parameterization of hydrological processes over a range of scales is examined, and the concept of the 'representative elementary area' (REA) is introduced. The REA is a fundamental scale for catchment modeling at which continuum assumptions can be applied for the spatially variable controls and parameters, and spatial patterns no longer have to be considered explicitly. The investigation of scale leads into the concept of hydrologic similarity in which the effects of the environmental controls on runoff generation and flood frequency response be investigated independently of catchment scale. The paper reviews the authors' initial results and hopefully will motivate others to also investigate the issues of hydrologic scale and similarity.

Wood, Eric F.↗

DOE OSTI2021

Role of underlying surface, rainstorm and antecedent wetness condition on flood responses in small and medium sized watersheds in the Yangtze River Delta region, China

Floods are becoming more frequent and are hard to control due to the shortage of water conservancy projects in small and medium sized basins, especially in developing regions. Understanding the hydrologic responses and estimated flood characteristics to storm events can help to predict flood disasters and form better effective mitigation and adaptation strategies. Here, using observations of several representative watersheds and Artificial Neural Networks (ANN) and Principal Component Regression (PCR) models, we conducted data-driven analyses to examine the flood responding characteristics. We quantified the relative contributions of the influencing factors to the variation of each flood characteristic. Statistical analysis of the observations shows that the drainage area plays a key role in determining the distribution of lag time and peak discharge. Rainstorm variability has direct influence on floods, and typhoon-induced rainstorms with high total rainfall and rainfall intensity generate higher lag time, flood peak, unit discharge and runoff depth, but lower runoff coefficient. The ANN and PCR models accurately predicted the variations of flood features using the driving factors including physical geographical characteristics, rainstorm features, and antecedent condition. Physical geographic characteristics are key influential factors of lag time, flood peak and runoff coefficient, while the rainstorm features control the magnitude of unit discharge and runoff depth. These results indicate that floods are mainly affected by rainstorm features and physical geographic characteristics in the Yangtze River delta, and it might become more damaging with the increasing rainfall extremes and sprawling impervious surfaces in a changing environment.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

Autonomous Inversion of In Situ Deformation Measurement Data for Injection-Induced Stress Change

Geologic carbon storage (GCS) is likely to play a key part of the global effort to dramatically reduce CO2 emissions and perhaps even reduce atmospheric CO2 concentrations through carbon negative operations. A critical part of effort to commercialize and widely deploy this technology is developing the capability to rapidly assimilate real-time monitoring data into a form that will enable site operators to make decisions to manage the safe and efficient operations. Two of the risks associate with GCS are the risk of inducing fractures in the sealing formations that can create leakage pathways and the risk of inducing earthquakes of sufficient magnitude to cause public concern, property damage, or safety risks. To properly manage these risks the site operator needs to know the initial state of stress, the change in stress induced by injection, and the relationship between operational parameters such as injection rate and pressure and the change in stress. Current methods of estimating the change in stress require choosing the type of constitutive model and the model parameters based on core, log, and geophysical data during the characterization phase, with little feedback from operational observations to validate or refine these choices. These characterization methods interrogate the geologic formations using length scales, loading rates or magnitudes that are quite different from those encountered by the actual storage system. It is shown that errors in the assumed constitutive response, even when informed by laboratory tests on core samples, are likely to be common, large, and underestimate the magnitude of stress change caused by injection. Recent advances in borehole-based strain instruments and borehole and surface-based tilt and displacement instruments have now enabled monitoring of the deformation of the storage system throughout its operational lifespan. This data can enable validation and refinement of the knowledge of the geomechanical properties and state of the system, but brings with it a challenge to transform the raw data into actionable knowledge. We demonstrate a method that uses automatic differentiation and a finite-element based geomechanical model perform a gradient-based deterministic inversion of geomechanical monitoring data. This approach allows autonomous integration of the instrument data without the need for time consuming manual interpretation and selection of updated model parameters. Furthermore, only isotropic linear elasticity is considered in this paper, the approach presented is very flexible as to what type of geomechanical constitutive response can be used. The approach is easily adaptable to nonlinear physics-based constitutive models to account for common rock behaviors such as creep and plasticity. The approach also enables training of machine learning-based constitutive models by allowing back propagation of errors through the finite element calculations. This enables strongly enforcing known physics, such as conservation of momentum and continuity, while allowing data-driven models to learn the truly unknown physics such as the constitutive or petrophysical responses.

Burghardt, Jeffrey A.↗

DOE OSTI2024

River–aquifer interactions enhancing evapotranspiration in a semiarid riparian zone: A modelling study

The hydrologic flows across the river–aquifer interface play an important role in groundwater dynamics and biogeochemical reactions within the subsurface; however, little is known about the effects of river–aquifer interactions on land surface processes. In this study, we developed a fully coupled three-dimensional (3D) land surface and subsurface model at a high resolution (~1 km) that accounts for high-frequency hydrologic exchange flow conditions to investigate how river–aquifer interactions modulate surface water budgets in the Upper Columbia-Priest Rapids watershed, a typical semiarid watershed located in the northwestern United States where river stage fluctuates in response to reservoir releases changing. Our results show that the spatiotemporal dynamics of river–aquifer interactions are highly heterogeneous, driven mainly by river-stage fluctuations. Adding 6.64 × 10 6 m 3 year –1 of water over the watershed from the river to groundwater owing to the lateral flow, river–aquifer interactions led to an increase in soil evaporation and transpiration supplied by higher soil moisture content, particularly in deeper subsurface. In a hypothetic future scenarios where a 5-m rise in river stage was assumed, the hydrologic flow exchange rates were intensified, resulting in higher surface water over the entire watershed. Overall, lateral flow induced by river–aquifer exchanges leads to an increase in evapotranspiration of ~75% in the historical period and of ~83% in the hypothetical future scenario. Finally, our study demonstrates the potential of coupled model as an effective tool for understanding river–aquifer–land surface interactions, and indicates that river–aquifer interactions fundamentally alter the water balance of the riparian zone for the semiarid watershed and will likely become more frequent and intense in the future under the effects of climate change.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Revealing Defect-Seeded and Interfacial Generation Mechanisms of Photoinduced Coherent Phonons with 4D Ultrafast Electron Microscopy

Ultrafast photoexcitation of coherent phonons is driven by an impulsive, collective displacement of constituent atoms from their average equilibrium lattice positions [1]. Models describing the generation of coherent acoustic modes typically invoke the creation of an anisotropic strain profile arising from the relatively instantaneous absorption of an ultrafast laser pulse [2]. If the skin depth is shallow relative to the specimen thickness, a steep tensile strain gradient perpendicular to the surface ($\frac{∂ε}{∂z}$) results. Initial relaxation occurs via rapid contraction of the surface layers followed by subsequent coherent oscillations of the lattice and launch of a train of coherent elastic strain waves [i.e., coherent acoustic phonons (CAPs)]. Macroscopically, responses are generally well-described by constitutive relations as gleaned from data gathered using ultrasonic methods or ultrafast spectroscopies. Furthermore, at the atomic to nanoscale level, individual lattice discontinuities and their impact on CAP behaviors can be modeled using multiscale methods [3,4]. Further, average unit-cell level responses on ultrafast timescales can be probed using femtosecond electron and X-ray scattering [5,6].

Flannigan, David J. [University of Minnesota, Minn↗

NASA NTRS1999

Seasonal-to-Interannual Variability in Antarctic Sea-Ice Dynamics, and Its Impact on Surface Fluxes and Water Mass Production

Strong seasonal and interannual signals in Antarctic bottom-water outflow remain unexplained yet are highly correlated with anomalies in net sea-ice growth in coastal polynyas. The mechanisms responsible for driving salination and replenishment and rejuvenation of the dense shelf "source" waters likely also generate pulses of bottom water outflow. The objective of this research is to investigate time-scales of variability in the dynamics of sea-ice in the Southern Ocean in order to determine the primary sites for production of dense shelf waters. We are using a merged satellite/buoy sea-ice motion data set for the period 1978-present day to compute the dynamics of opening and closing of coastal polynyas over the continental shelf. The Ocean Circulation and Climate Advanced Model (OCCAM) ocean general circulation model with coupled sea-ice dynamics is presently forced using National Center for Environmental Prediction (NCEP) data to simulate fluxes and the salination impact of the ocean shelf regions. This work is relevant in the context of measuring the influence of polar sea-ice dynamics upon polar ocean characteristics, and thereby upon global thermohaline ocean circulation. Interannual variability in simulated net freezing rate in the Southern Weddell Sea is shown for the period 1986-1993. There is a pronounced maximum of ice production in 1988 and minimum in 1991 in response to anomalies in equatorward meridional wind velocity. This follows a similar approximate 8-year interannual cycle in Sea Surface Temperature (SST) and satellite-derived ice-edge anomalies reported elsewhere as the "Antarctic Circumpolar Wave." The amplitude of interannual fluctuations in annual net ice production are about 40% of the mean value, implying significant interannual variance in brine rejection and upper ocean heat loss. Southward anomalies in wind stress induce negative anomalies in open water production, which are observed in passive microwave satellite images. Thus, cycles of enhanced poleward wind stress reduce ice growth by compacting the ice along the coastline and closing open water in leads and polynyas. Model simulations confirm that years of low ice production, such as 1991, coincide with years of lower than normal bottom water outflow. Future plans include the assimilation of satellite ice concentrations and ice drift dynamics to more accurately constrain boundary conditions in the model.

Drinkwater, Mark R.↗

NASA NTRS2011

Warm and Saline Events Embedded in the Meridional Circulation of the Northern North Atlantic

Ocean state estimates from 1958 to 2005 from the Simple Ocean Assimilation System (SODA) system are analyzed to understand circulation between subtropical and subpolar Atlantic and their connection with atmospheric forcing. This analysis shows three periods (1960s, around 1980, and 2000s) with enhanced warm, saline waters reaching high latitudes, alternating with freshwater events originating at high latitudes. It complements surface drifter and altimetry data showing the subtropical -subpolar exchange leading to a significant temperature and salinity increase in the northeast Atlantic after 2001. The warm water limb of the Atlantic meridional overturning cell represented by SODA expanded in density/salinity space during these warm events. Tracer simulations using SODA velocities also show decadal variation of the Gulf Stream waters reaching the subpolar gyre and Nordic seas. The negative phase of the North Atlantic Oscillation index, usually invoked in such variability, fails to predict the warming and salinization in the early 2000s, with salinities not seen since the 1960s. Wind stress curl variability provided a linkage to this subtropical/subpolar gyre exchange as illustrated using an idealized two ]layer circulation model. The ocean response to the modulation of the climatological wind stress curl pattern was found to be such that the northward penetration of subtropical tracers is enhanced when amplitude of the wind stress curl is weaker than normal. In this case both the subtropical and subpolar gyres weaken and the subpolar density surfaces relax; hence, the polar front moves westward, opening an enhanced northward access of the subtropical waters in the eastern boundary current.

Hakkinen, Sirpa↗

NASA NTRSDate not supplied

Mineralogical, Elemental, and Tomographic Reconnaissance Investigation for CLPS (METRIC): A Proposed Mission to the Lunar South Pole-Aitken Basin

The global Lunar Magma Ocean (LMO) hypothesis is based on analyses of ferroan anorthosites (FAN), magnesian plutonic rocks (Mg-suite), and KREEPy (Potassium - K; Rare Earth Elements - REE; Phosphorous - P) material present at the Apollo landing sites and the assumed global distribution of these lithologies. However, orbital spacecraft data over the last two decades have highlighted that the location of the Apollo missions within and near the Procellarum KREEP Terrane (PKT) is compositionally anomalous and not representative of the entire lunar surface and suggests that Apollo samples provide a biased view of the Moon. Further, meteorites, which originate from random areas on the lunar surface, including those not visited by Apollo, Luna, or Chang’E-5, provide complementary datasets that reveal a more complex lunar formation and evolution than the simple LMO onion skin model. While many recent investigations support and refine the LMO model, other recent studies have questioned whether this mechanism alone is responsible for the primordial differentiation of the Moon. Thus, relating the distinct and asymmetrically distributed geochemical terranes on the Moon, especially KREEP which is a predicted global LMO product, to lunar formation and differentiation remains a fundamental goal of lunar science. The South Pole-Aitken (SPA) impact carved a two-and-a-half thousand-kilometer basin into the Moon’s ancient crust, excavating material from as deep as the upper mantle. Orbital remote sensing of the basin indicates compositions with elevated thorium (Th) contents and unique mineralogies that provide a window into early lunar structure and LMO processes. We propose to land a powerful characterization laboratory at Birkeland crater on the far side of the Moon–where Th-rich SPA ejecta have been re-exposed near the surface–in order to address outstanding questions about the evolution of the Moon.

K M Cannon↗

DOE OSTI2024

Dynamics of Electric Polarization and Relaxation of Ions at Humid Calcite Surfaces

Mobile ions at mineral surfaces can respond to an applied electric field, adopting a new distribution that effectively represents polarization of the electrical double layer. When the field is released, the ions relax to their equilibrium distribution. In both cases, the dynamics are characteristic of the interface. However, current models of electrokinetic phenomena are not sufficiently robust to accurately predict collective ion dynamics at structurally and chemically specific mineral–water interfaces. Here, in this study, we use electrostatic force microscopy (EFM) to investigate the dynamics of ion relaxation at hydrated calcite (104) surfaces at controlled relative humidity (RH). Electrically biased probes are used to polarize the distributions of calcium and carbonate ions that are intrinsic to this interface across a range of RH values. Polarization kinetics are tracked by monitoring the tip–sample force gradient during charging, and EFM imaging is used to characterize the spatial relaxation dynamics after the applied field is released. Electrostatic finite element modeling of the sample/probe system across length-scales from nanometers to millimeters reproduces the observed stretched exponential charging response. Together, these results allow us to estimate the ion diffusivities at the interface across a wide range of RH values. These diffusivities increase by roughly 5 orders of magnitude as the RH is increased from 5 to 90%, highlighting the critical role of adsorbed water for surface ion solvation that enables ion mobility.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS1985

Numerical study of the response of a compressible, turbulent boundary layer to a short region of surface curvature

Calculations were performed to investigate the supersonic flow of a turbulent boundary layer over short regions of concave surface curvature. Upstream of each curved surface the freestream Mach number was 2.9. Three different constant radii of curvature models were investigated to cover a range of curvatures and turning angles. The numerical technique solved the full, Reynolds-averaged Navier-Stokes equations using two different turbulence models: the algebraic eddy viscosity model due to Baldwin and Lomax (1978), and the one-equation model to Rubesin (1976). The calculations were compared with the experimental data of Taylor and Smits (1984), and the agreement was very satisfactory, especially for the computations using the one-equation model.

Degani, D.↗

DOE OSTI2023

Surface enrichment dictates block copolymer orientation

Orientation of block copolymer (BCP) morphology in thin films is critical to applications as nanostructured coatings. Despite being well-studied, the ability to control BCP orientation across all possible block constituents remains challenging. Here, in this study, we deploy coarse-grained molecular dynamics simulations to study diblock copolymer ordering in thin films, focusing on chain makeup, substrate surface energy, and surface tension disparity between the two constituent blocks. We explore the multi-dimensional parameter space of ordering using a machine-learning approach, where an autonomous loop using a Gaussian process (GP) control algorithm iteratively selects high-value simulations to compute. The GP kernel was engineered to capture known symmetries. The trained GP model serves as both a complete map of system response, and a robust means of extracting material knowledge. We demonstrate that the vertical orientation of BCP phases depends on several counter-balancing energetic contributions, including entropic and enthalpic material enrichment at interfaces, distortion of morphological objects through the film depth, and of course interfacial energies. BCP lamellae are found more resistant to these effects, and thus more robustly form vertical orientations across a broad range of conditions; while BCP cylinders are found to be highly sensitive to surface tension disparity.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

NASA NTRS1977

ALTKAL: An optimum linear filter for GEOS-3 altimeter data

ALTKAL is a computer program designed to smooth sea surface height data obtained from the GEOS 3 altimeter, and to produce minimum variance estimates of sea surface height and sea surface slopes, along with their standard derivations. The program operates by processing the data through a Kalman filter in both the forward and backward directions, and optimally combining the results. The sea surface height signal is considered to have a geoid signal, modeled by a third order Gauss-Markov process, corrupted by additive white noise. The governing parameters for the signal and noise processes are the signal correlation length and the signal-to-noise ratio. Mathematical derivations of the filtering and smoothing algorithms are presented. The smoother characteristics are illustrated by giving the frequency response, the data weighting sequence and the transfer function of a realistic steady-state smoother example. Based on nominal estimates for geoidal undulation amplitude and correlation length, standard deviations for the estimated sea surface height and slope are 12 cm and 3 arc seconds, respectively.

Fang, B. T.↗

NASA NTRS1979

Skynet satellite electron precharging experiments

Large surface areas (approximately 1 sq m) of the Skynet I qualification model satellite were exposed to the environments of a flux of monoenergetic electrons and pulsed high-intensity X-ray photons from an exploding-wire-radiator source, separately and simultaneously. Experiments were performed with both the cylindrical solar cell panels and the end thermal blanket exposed to these environments. The satellite was instrumented with fiber-optics isolated fast-response sensors capable of sensing and recording the time-varying electric fields and surface currents on the satellite. Spontaneous discharges of the two surfaces were characterized prior to measurements of the interaction of the system-generated electromagnetic pulse from the X-ray photon burst with the electron precharged surfaces. A two-dimensional self-consistent finite difference computer code was used to predict the currents on the Skynet both with and without precharging of the thermal-blanket. These calculations, which used an thermal-blanket voltage profile, predicted a factor of 2 increase in satellite exterior axial currents.

Vanlint, V. A. J.↗

NASA NTRS1989

A study of the effect of synoptic scale processes in GCM modelling

Research was conducted to help modeling groups at NASA to develop better weather forecasting and general circulation models (GCM) for activities relating to the meteorological uses of satellite data. The focus was on the physical processes that were being simulated by models: radiative effects and latent heat release associated with clouds; orographic influences; and heat transfer at the ocean and ice surfaces. An attempt was made to deduce the role of diabatic heating in North Atlantic cyclogenesis and in the global heat budget. Inferences were made in four studies: heat budget statistics from GCM assimilations; dynamics of north Atlantic cyclones; Cage-type energy budget calculations; and grid scale cloud formation. Mechanisms that were responsible for the variability and structure of the atmospheric on a hemispheric scale were studied by a hybrid of statistical analysis and theoretical modeling. Variability and structure are both related to synoptic scale processes through baroclinic and barotropic energy transformations.

Herman, Gerald F.↗

DOE OSTI2023

Evaluation of aerosol–cloud interactions in E3SM using a Lagrangian framework

A Lagrangian framework is used to evaluate aerosol–cloud interactions in the U.S. Department of Energy's Energy Exascale Earth System Model (E3SM) version 1 (E3SMv1) for measurements taken at Graciosa Island in the Azores where a U.S. Department of Energy Atmosphere Radiation Measurement (ARM) site is located. This framework uses direct measurements of cloud condensation nuclei (CCN) concentration (instead of relying on satellite retrievals of aerosol optical depth) and incorporates a suite of ground-based ARM measurements, satellite retrievals, and meteorological reanalysis products that when applied to over a 1500 trajectories provides key insights into the evolution of low-level clouds and aerosol radiative forcing that is not feasible from a traditional Eulerian analysis framework. Significantly lower concentrations (40 %) of surface CCN concentration are measured when precipitation rates in 48 h back trajectories average above 1.2 mm d -1 in the Integrated Multi-satellitE Retrievals for Global Precipitation Measurement (IMERG) product. The depletion of CCN concentration when precipitation rates are elevated is nearly twice as large in the ARM observations compared to E3SMv1 simulations. The model CCN concentration bias remains significant despite modifying the autoconversion and accretion rates in warm clouds. As the clouds in trajectories associated with larger surface-based CCN concentration advect away from Graciosa Island, they maintain higher values of droplet number concentrations (N d ) over multiple days in observations and E3SM simulations compared to trajectories that start with lower CCN concentrations. The response remains robust even after controlling for meteorological factors such as lower troposphere stability, the degree of cloud coupling with the surface, and island wake effects. E3SMv1 simulates a multi-day aerosol effect on clouds and a Twomey radiative effect that is within 30 % of the ARM and satellite observations. However, the mean cloud droplet concentration is more than 2–3 times larger than in the observations. While Twomey radiative effects are similar amongst autoconversion and accretion sensitivity experiments, the liquid water path and cloud fraction adjustments are positive when using a regression model as opposed to negative when using the present-day minus pre-industrial aerosol emissions approach. This result suggests that tuning the autoconversion and accretion alone is unlikely to produce the desired aerosol susceptibilities in E3SMv1.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS1994

Soft X-ray spectral fits of Geminga with model neutron star atmospheres

The spectrum of the soft X-ray pulsar Geminga consists of two components, a softer one which can be interpreted as thermal-like radiation from the surface of the neutron star, and a harder one interpreted as radiation from a polar cap heated by relativistic particles. We have fitted the soft spectrum using a detailed magnetized hydrogen atmosphere model. The fitting parameters are the hydrogen column density, the effective temperature T(sub eff), the gravitational redshift z, and the distance to radius ratio, for different values of the magnetic field B. The best fits for this model are obtained when B less than or approximately 1 x 10(exp 12) G and z lies on the upper boundary of the explored range (z = 0.45). The values of T(sub eff) approximately = (2-3) x 10(exp 5) K are a factor of 2-3 times lower than the value of T(sub eff) obtained for blackbody fits with the same z. The lower T(sub eff) increases the compatibility with some proposed schemes for fast neutrino cooling of neutron stars (NSs) by the direct Urca process or by exotic matter, but conventional cooling cannot be excluded. The hydrogen atmosphere fits also imply a smaller distance to Geminga than that inferred from a blackbody fit. An accurate evaluation of the distance would require a better knowledge of the ROSAT Position Sensitive Proportional Counter (PSPC) response to the low-energy region of the incident spectrum. Our modeling of the soft component with a cooler magnetized atmosphere also implies that the hard-component fit requires a characteristic temperature which is higher (by a factor of approximately 2-3) and a surface area which is smaller (by a factor of 10(exp 3), compared to previous blackbody fits.

Meyer, R. D.↗

DOE OSTI2026

Deep Gaussian process-based cost-aware batch Bayesian optimization for complex materials design campaigns

The accelerating pace and expanding scope of materials discovery demand optimization frameworks that efficiently navigate vast design spaces with complex response surfaces while judiciously allocating limited evaluation resources. We present a cost-aware, batch Bayesian optimization scheme powered by deep Gaussian process (DGP) surrogates and a heterotopic querying strategy. Our DGP surrogate, formed by stacking GP layers, models complex hierarchical relationships among high-dimensional compositional features and captures correlations across multiple target properties, propagating uncertainty through successive layers. We integrate evaluation cost into an upper-confidence-bound acquisition extension, which, together with heterotopic querying, proposes small batches of candidates in parallel, balancing exploration of under-characterized regions with exploitation of high-mean, low-variance predictions across correlated properties. Applied to refractory high-entropy alloys for high-temperature applications, our framework converges to optimal formulations in fewer iterations with cost-aware queries than conventional GP-based BO, highlighting the value of deep, uncertainty-aware, cost-sensitive strategies in materials campaigns.

36 MATERIALS SCIENCE↗