Search NASA⌕ Search

SEARCH · Search NASA

Results for “surface parameters”

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

Spin and lattice dynamics of the two-dimensional van der Waals ferromagnet CrI3

Chromium tri-iodide (CrI 3 ) is a prototypical ferromagnetic van der Waals insulator with its genuinely two-dimensional (2D) long-range magnetic order below 45 K demonstrated recently. The underlying magnetic anisotropy has not been completely understood while both the Dzyaloshinskii-Moriya (DM) interaction and the Kitaev—Γ type interaction have been considered as the relevant magnetic Hamiltonian. In addition, the relation between the crystal structure and the magnetic order needs to be further elucidated concerning their possible coupling in few-layer samples and in the topmost surface layers of bulk samples. Here, we investigate these issues via temperature- and magnetic field-dependent terahertz spectroscopy on bulk CrI 3 single crystals, focusing on the dynamics of ferromagnetic resonances (FMRs) and optical phonons in the terahertz (THz) region from 4 to 120 cm -1 (from 0.5 to 15 meV). We narrow down the possible ranges of the interaction parameters such as the off-diagonal symmetric terms and the single-ion anisotropy. The accurate values of these parameters significantly constrain the magnitude of possible Kitaev—Γ exchange interaction and the topological magnon gap. Moreover, the structure-magnetism relationship was critically analyzed based on the temperature- and field-dependences of two E u in-plane optical phonon modes, which shows that the commonly believed structural phase diagram of CrI 3 , derived from surface-preferential data, has to be seriously modified.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Macroscale properties and atomic-scale mechanisms of ash removal in low-temperature hydrothermal carbonization

Biogenic ash is a significant impediment to the utilization of agricultural residues in biofuel production. Such challenge can be addressed by various treatments, as demonstrated in this study on the experimental and computational mechanisms involved in the hydrothermal treatment (HT) of wheat straw. A combination of classical (all-atom) molecular dynamics simulations of cellulose carrying silica and calcium species, along with first principles quantum chemical calculations, indicates the dissociation of inorganics from the cellulose with increased HT temperature. Here, this observation is confirmed by experimental evidence of effective ash removal by HT, showing at least 50% removal of sulfur, chlorine, potassium, and calcium, and 12.5% of silica, leading to a reduced total ash content (from 6.7% to 4.2%).Changes in structural features upon HT, such as surface cellular structure and porosity, were revealed, accompanied by an increased specific surface area (from 1.17 to 6.34 m 2 /g). Our simulations suggest that silica binds tightly to the hydrophobic face of cellulose at room temperature, but HT significantly reduces the binding free energy of association with both hydrophobic and hydrophilic surfaces. Most significantly, ash removal leads to an increased calorific value, rising from approximately 16 MJ/kg to about 19 MJ/kg, along with improved thermal behavior. The improved integration combustion index parameter S indicates that the combustion properties improve with ash removal efficiency. The proposed atomic-level mechanism for the observed removal of inorganics during mild HT underscores the potential of such treatment in producing energy-dense wheat straw, a widely available agricultural residue.

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

Relative permeabilities for two-phase flow through wellbore cement fractures

Multiple fluids are likely to exist in fractures and flow paths associated with leaky wellbores, including liquids (e.g., crude oil) and gases (e.g., gas exsolved from liquid). These fluids occupy and move through different portions of the pore spaces within the fractures depending on many factors, including fluid properties, fracture size, and the amount of the different fluids. Upward leakage of any phase, through the fracture, can contaminate water-bearing formations, create hazardous surface conditions, and compromise the functionality of the wellbore. Early signs of wellbore leaks may be expressed by anomalous pressure behavior at surface monitoring points on cavern storage wells. These pressure anomalies are difficult to interpret, necessitating knowledge of the factors that affect the multiphase flow in fractures and porous media. These parameters are critical to modeling multiphase flow in fractures. This insight can guide further diagnosis and maximize leak remediation. Here, our study focuses on the relationship of the liquid–gas relative permeabilities for representative variable-aperture wellbore cement fracture. To obtain the relative permeability of each phase, two-phase flow tests were conducted where both fluids were flowing simultaneously through a fractured wellbore cement specimen under a range of factors, namely (1) aperture size, (2) capillary numbers, and (3) viscosity ratio. The flow experiments were conducted under a range of confining stresses and flow velocities, using nitrogen gas and silicone oils (of different viscosities) in a specially designed pressure vessel. The sum of gas and oil relative permeabilities were found to be less than one under all conditions, which indicates that the presence of one phase affects the permeability of the other phase, and vice versa. Since the gas phase flow conditions include a significant inertial flow component in addition to viscous flow, the inertial flow coefficients at different saturation states are presented. The factors affecting the relationship between the relative permeabilities are discussed in detail. A new mathematical model for estimating the relative permeability of wellbore cement fracture is presented and experimentally validated.

58 GEOSCIENCES↗
DOE OSTI2024

Surface anchoring requirements for vanadia clusters on titanium oxide surfaces and their impact on activity for oxidative dehydrogenation of ethanol

Titanium oxide site requirements for the anchoring of catalytically active forms of dispersed vanadia species are investigated. Using a series of oxygen-modified titanium nitride materials, we systematically modify the available vanadia-anchoring moieties on the titanium-bearing support surface. This strategy in turn results in a very narrow size distribution of well-dispersed vanadia clusters anchored on the support. This synthesis technique ensures that only highly dispersed vanadia is present on the catalyst surface even at nominal monolayer coverages. A deeper insight into the catalytic behavior for ethanol partial oxidation of the active vanadia species is obtained by comparing intrinsic kinetic and thermodynamic parameters of kinetically relevant reaction steps (heat of ethanol adsorption, hydrogen abstraction activation energy, and the energetics of catalyst reoxidation). In conclusion, it is found that these parameters are dependent on the relative amount of oxygen in the titanium-bearing support and the resulting distribution of vanadia species, which validates the premise of a potential participation of lattice oxygen atoms of the titania support into the catalytic activity of active vanadia species.

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

Correlating Nb-SRF Surface Processing with Evolution of Surface Electronic States

The few nanometers of the surface exposed to RF field plays a major role in defining the RF performance of superconducting cavities. Over the past two decades, several pioneering surface treatment and processing methods have emerged, enabling remarkable improvements in cavity performance by simultaneously achieving high Q with increasing Eacc. These processing methods include: thermal treatment under ultra-high vacuum (UHV) conditions across lo¬¬¬¬¬¬¬w-, mid-, and high-temperature ranges and high temperature treatments under controlled N2 atmosphere. These processes also produce distinct surface oxide configurations with different valence states, thicknesses, and uniformity, as well as different oxygen concentration profiles in bulk Nb. In this work, we are trying to understand how do surface-processing methods and the resulting oxide/oxygen profiles affect the electronic structure of surface and the mechanism of superconductivity? With the help of Fermilab’s in-house X-ray photoemission facility and, in collaboration with the synchrotron-based angle-resolved photoemission (ARPES) facility at Argonne National Laboratory, we are investigating how the valence band structure and density of states (DoS) near the Fermi level modify with different surface treatments. Our observations show that different surface-processing methods lead to distinct evolutions of the valence-band states near the Fermi level during the superconducting transition. This behavior suggests variations in Nb-O orbital hybridizations and points towards the possibility of different underlying mechanisms of superconductivity governed by the surface chemistry and oxide configuration. We also correlate these distinct superconducting mechanisms with RF cavity performance, specifically focusing on measured surface resistance, the nature of the Q-slope, and quench fields observed in SRF measurements. These results will enable us to identify the potential limiting factors and relevant controllable parameters that can be further optimized to improve the performance of SRF cavities.

Tripathi, Malvika [Fermilab]↗
DOE OSTI2024

Deep-learning-derived planetary boundary layer height from conventional meteorological measurements

Abstract. The planetary boundary layer (PBL) height (PBLH) is an important parameter for various meteorological and climate studies. This study presents a multi-structure deep neural network (DNN) model, which can estimate PBLH by integrating the morning temperature profiles and surface meteorological observations. The DNN model is developed by leveraging a rich dataset of PBLH derived from long-standing radiosonde records augmented with high-resolution micro-pulse lidar and Doppler lidar observations. We access the performance of the DNN with an ensemble of 10 members, each featuring distinct hidden-layer structures, which collectively yield a robust 27-year PBLH dataset over the southern Great Plains from 1994 to 2020. The influence of various meteorological factors on PBLH is rigorously analyzed through the importance test. Moreover, the DNN model's accuracy is evaluated against radiosonde observations and juxtaposed with conventional remote sensing methodologies, including Doppler lidar, ceilometer, Raman lidar, and micro-pulse lidar. The DNN model exhibits reliable performance across diverse conditions and demonstrates lower biases relative to remote sensing methods. In addition, the DNN model, originally trained over a plain region, demonstrates remarkable adaptability when applied to the heterogeneous terrains and climates encountered during the GoAmazon (Green Ocean Amazon; tropical rainforest) and CACTI (Cloud, Aerosol, and Complex Terrain Interactions; middle-latitude mountain) campaigns. These findings demonstrate the effectiveness of deep learning models in estimating PBLH, enhancing our understanding of boundary layer processes with implications for improving the representation of PBL in weather forecasting and climate modeling.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

LeaPP: Learning Pathways to Polymorphs through Machine Learning Analysis of Atomic Trajectories

Understanding the mechanisms underlying crystal nucleation and growth is crucial for many technological applications. Due to the short length and time scales involved, crystal nucleation is often studied using molecular simulations. Most existing approaches to extract the nucleation mechanism from simulations focus on the analysis of static snapshots of the configurations, potentially overlooking subtle local fluctuations and the history of the particles involved in the formation of solid nuclei. Here, in this work, we propose a novel methodology called LeaPP that categorizes nucleation trajectories based on the temporal information of their constituent particles. We leverage the time evolution of the local environment of the crystallizing particles to encapsulate the relationship between the structure and dynamics and distinguish between different evolving particle paths. Identification of the distinct particle paths further enables characterizing the nucleation trajectories into different pathways. Collectively, LeaPP provides a more nuanced understanding of nucleation through an unsupervised approach with lesser dependence on traditional order parameters. Furthermore, the pathways identified by LeaPP are predictive of the resulting polymorph. We demonstrate LeaPP on three different systems─Lennard-Jones-like particles, Ni 3 Al, and water on surfaces. The general methodology underlying LeaPP─considering the time evolution of the building blocks─applies to a wide range of self-assembly problems.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Crystallographic Evidence of Size-Dependent Bond Flexibility in Metal–Organic Framework Nanocrystals

Size-dependent electronic, magnetic, and optical behavior suggests that metal–organic frameworks become softer materials as their particle sizes decrease, but direct evidence is lacking. Here, we report variable-temperature powder X-ray diffraction data of Fe(1,2,3-triazolate)2 particles that offer crystallographic insight into size-dependent bond flexibility. Rietveld refinement reveals size-dependent positive thermal expansion upon downsizing the crystalline domains from 178 to 9 nm, with a 6-fold increase from 16 MK –1 to 96 MK –1 . Here, this behavior occurs in tandem with size-dependent elongation of metal–ligand bonds and increasing thermal displacement parameters, consistent with pronounced metal-linker bond lability. We propose that these effects, as well as size-dependent annealing of crystallite sizes, originate from the high charge density and surface stress of smaller particles. Taken together, these results provide structural evidence that size reduction serves as a synthetic route to controlling the dynamic response of materials to external stimuli.

Crystals↗
DOE OSTI2025

A Machine Learning Framework for Predicting Microphysical Properties of Ice Crystals From Cloud Particle Imagery

The microphysical properties of ice crystals are important because they significantly alter the radiative properties and spatiotemporal distributions of clouds, which in turn strongly affect Earth's climate. However, it is challenging to measure key properties of ice crystals, such as mass or morphological features. Here, we present a proof-of-concept framework for predicting three-dimensional (3D) microphysical properties of ice crystals from in situ two-dimensional (2D) imagery. First, we computationally generated synthetic ice crystals using 3D modeling software along with geometric parameters estimated from the 2021 Ice Cryo-Encapsulation Balloon (ICEBall) field campaign. Then, we used synthetic crystals to train machine learning (ML) models to predict effective density ($ρ_e$), effective surface area ($A_e$), and number of bullets ($N_b$) from synthetic rosette imagery. On unseen synthetic images, our ML models accurately predicted ice crystal properties. ResNet-18 performed best, achieving $R^2$ values of 0.99 and 0.98 for $ρ_e$ and $A_e$, respectively, and MAE of 0.10 for mathematical equation in single view tasks. Stereo view ResNet-18 further reduced RMSE by 40% for $ρ_e$ and $A_e$ and reduced MAE by 0.08 for $N_b$. This work provides a novel ML-driven framework for estimating ice microphysical properties from in situ imagery, which will allow for downstream constraints on microphysical parameterizations, such as the mass-size relationship.

Ko, J. [Columbia Univ., New York, NY (United State↗
DOE OSTI2026

Modeling rf sheath formation in turbulent tokamak boundary plasma

During ICRF antenna operation, complex interactions between turbulent density profiles, nonlinear RF sheaths, and RF-induced convective transport are observed to alter plasma density in the tokamak edge [D’Ippolito et al., Nucl. Fusion 38, 1543 (1998)]. In this work, we explore the physics of such interactions via numerical modeling, using a nonlinear EM/plasma/sheath code (VSim) and profiles obtained from a fluid plasma turbulence code (Hermes) in a 3D slab domain containing biased side-wall limiters. RF-rectified sheath formation on antenna and limiter surfaces is observed as electromagnetic waves launched by the antenna are refracted through the turbulent density profile. On transport timescales, such sheath potentials have been shown to influence both the mean species density and its RMS fluctuation spectrum [Smithe et al., these proceedings]. On the faster RF timescales, we demonstrate that the converse is also true – regions of high plasma density near material surfaces give rise to the highest sheath potential amplitudes. When density is turbulent and spatially nonuniform, localized regions of high sheath potential (hotspots) may develop where high-density filaments intersect material surfaces. Such hotspots are of particular concern as sources of impurity sputtering, and we explore their behavior in response to changes both to the local plasma density and to antenna operating parameters and structure. Related results exploring the role of Faraday shields and/or enclosing structures in suppressing high sheath potentials for other devices (e.g. SPARC) will also be shown.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Correlating Surface Processing of Nb Superconducting RF Cavities with the Evolution of Surface Electronic States

Superconducting-radio frequency (RF) cavities provide an efficient way to accelerate particle beams with extremely high acceleration gradients while generating very small power dissipation. The few nanometers of the surface play a critical role in defining the RF performance of superconducting Nb based cavities. Over the past two decades, several pioneering surface treatment and processing methods have emerged, enabling remarkable improvements in superconducting cavity performance by simultaneously achieving high quality factors with increasing maximum acceleration gradients. These processing approaches include chemical polishing, distinct multi-step thermal treatments under ultra-high vacuum (UHV) conditions over low to high temperature regimes, as well as high-temperature treatments under controlled nitrogen atmospheres. Beyond their macroscopic impact on RF performance, these methods produce distinct surface oxide configurations characterized by different valence states, oxide thicknesses, chemical uniformity, and oxygen concentration profiles extending into the near-surface bulk of niobium. In this work, we are trying to understand how the surface-processing methods and the resulting oxide/oxygen profiles affect the electronic structure of surface and the mechanism of superconductivity. Using a combination of X-ray photoemission and X-ray absorption spectroscopies, we investigate how the valence-band structure and the electronic density of states (DoS) near the Fermi level evolve under different surface treatments. By employing tunable photon energies across multiple elemental absorption edges, we use resonant photoemission to disentangle and identify the elemental contributions to specific valence-band features. Our observations show that different surface-processing methods lead to distinct temperature evolutions of the DoS and valence-band states near the Fermi level. Our results suggest variations in different Nb-O orbital hybridizations in distinct processes and point towards the possibility of different underlying mechanisms of superconductivity governed by surface chemistry and oxide configuration. We also correlate these distinct superconducting mechanisms with RF cavity performance, specifically focusing on measured surface resistance, the nature of the Q-slope, and quench fields observed in superconducting RF measurements. These results will enable us to identify the potential limiting factors and relevant controllable parameters that can be further optimized to improve the performance of superconducting RF cavities.

Tripathi, Malvika [Fermilab] (ORCID:00000001989251↗
DOE OSTI2024

An Integrated Electrochemical Approach to the Precision Synthesis of Sustainable Catalyst Materials

There is a pressing need to replace critical materials such as platinum group elements (PGE) in applications related to energy. This work responds to the need for fundamental principles of materials design to replace these metals by addressing a gap in the understanding of how to precisely control the surface structure of nanoscale materials. Electrochemical nanomaterials synthesis, such as techniques developed in the PI’s research group, expands the toolbox of available synthetic handles to include both electrochemical and chemical parameters, providing access to control over synthetic conditions in ways that are not possible in purely chemical nanoparticle growth. However, electrochemical materials synthesis (electrodeposition) is inherently serial, limiting throughput and preventing widespread implementation. Initial efforts in the short period of this award resulted in an innovative, high throughput, parallel approach to synthetic discovery for the electrodeposition of shaped metal nanoparticles that overcomes this limitation. Looking ahead, this work establishes an important capability that will enable researchers to push boundaries in the shape control of nanoparticles composed of non-PGE such as copper, iron, nickel, or cobalt.

36 MATERIALS SCIENCE↗
DOE OSTI2024

A Deep Learning-Driven Sampling Technique to Explore the Phase Space of an RNA Stem-Loop

The folding and unfolding of RNA stem-loops are critical biological processes; however, their computational studies are often hampered by the ruggedness of their folding landscape, necessitating long simulation times at the atomistic scale. Here, we adapted DeepDriveMD (DDMD), an advanced deep learning-driven sampling technique originally developed for protein folding, to address the challenges of RNA stem-loop folding. Although tempering- and order parameter-based techniques are commonly used for similar rare-event problems, the computational costs or the need for a priori knowledge about the system often present a challenge in their effective use. DDMD overcomes these challenges by adaptively learning from an ensemble of running MD simulations using generic contact maps as the raw input. DeepDriveMD enables on-the-fly learning of a low-dimensional latent representation and guides the simulation toward the undersampled regions while optimizing the resources to explore the relevant parts of the phase space. We showed that DDMD estimates the free energy landscape of the RNA stem-loop reasonably well at room temperature. Our simulation framework runs at a constant temperature without external biasing potential, hence preserving the information on transition rates, with a computational cost much lower than that of the simulations performed with external biasing potentials. Here, we also introduced a reweighting strategy for obtaining unbiased free energy surfaces and presented a qualitative analysis of the latent space. This analysis showed that the latent space captures the relevant slow degrees of freedom for the RNA folding problem of interest. Finally, throughout the manuscript, we outlined how different parameters are selected and optimized to adapt DDMD for this system. We believe this compendium of decision-making processes will help new users adapt this technique for the rare-event sampling problems of their interest.

Gupta, Ayush↗
DOE OSTI2025

The relative importance of building design parameters in reducing energy use and sensible heat release from buildings in light of forecasted future weather data and building coverage ratio

Buildings typically have a 60-to-75-year lifespan before they require significant maintenance or modifications. However, most builders evaluate the performance of their new buildings using whole-building energy simulation tools based on the current typical meteorological year (TMY) file or actual meteorological year. The energy use consumption and sensible heat release pattern observed from buildings could potentially change based on shifting global climates. Therefore, the recommended energy-efficiency design parameters might also change during these periods. In this study, we evaluate the role of different building design parameters, such as material reflectivity, HVAC COP, and insulation values, on building energy usage and sensible heat release from buildings with different building coverage ratios (BCR), based on the current and future weather file TMY (fTMY) for the middle of the century (2040–2060). The role of sensible heat release from buildings is not accounted for accurately while estimating building energy usage in most whole-building energy simulations. The study conducts a series of whole-building energy simulation analyses using EnergyPlus to evaluate the role of different design parameters based on TMY and fTMY weather conditions. The analysis is conducted for two hot desert climatic cities: Phoenix (USA) and Abu Dhabi (UAE). The results show that, for the base case in a future climate, the sensible heat release is reduced by an average of 30% due to the reduced delta T between the surface and ambient air. Further, the results show an increase in total energy consumption by 5% annually. The results also show that, for buildings with traditional coatings, shorter buildings release more heat than taller buildings. On the other hand, for buildings with reflective paints, shorter buildings release less heat than taller buildings. The findings from this study can be used by policymakers, utility companies, and builders to better understand the relative role of different building design parameters while constructing new and retrofitting existing buildings.

Alhazmi, Mansour [King Fahd University of Petroleu↗
DOE OSTI2025

Amine Structure Governs Corrosion Rates of Copper Catalysts in Electrochemical Reactive Capture of CO 2

Reactive capture of CO 2 (RCC) offers an integrated approach that combines CO 2 capture with its direct electrochemical conversion, eliminating the need for CO 2 release from the capture agent. By avoiding the pH, pressure, and temperature swings required for the release step, RCC has the potential to reduce both energy consumption and capital costs compared to the conventional sequential process of CO 2 capture, release, concentration, and conversion. Amines, widely used in industrial CO 2 capture, face challenges in RCC systems due to their incompatibility with transition metal catalysts as well as their tendency to promote electrode corrosion and parasitic hydrogen evolution. Identifying suitable combinations of amines and catalysts is therefore critical to enabling integrated CO 2 capture and conversion. Here, this work systematically investigates the performance of four primary and four secondary amines for RCC on polycrystalline Cu catalysts. Among the eight tested amines, only dimethylamine showed no measurable Cu corrosion near the open circuit potential. In contrast, ammonia, methylamine, ethylamine, monoethanolamine, diethylamine, diethanolamine, and piperazine all induced Cu corrosion. Corrosion rates correlate with the pK a and steric hindrance of the amines, highlighting key parameters for catalyst–amine codesign. Grand canonical DFT calculations indicate a correlation between the adsorption strength of protonated amines, their pK a , and the extent of Cu corrosion, suggesting that both the surface binding of protonated amines and the lability of their protons play critical roles in corrosion acceleration near open circuit potentials. These finding suggest that amines with high pK a values and weak binding of their protonated forms to Cu surfaces are preferred, as they offer better corrosion resistance.

Choi, Jounghwan [Univ. of California, Los Angeles,↗
DOE OSTI2025

Microwave-Assisted Production of Polycarbonate Diols using Carbon Dioxide

Carbon dioxide (CO2) is a cheap and readily available resource that can be converted to value-added chemicals including fuels, polymers, and other products. One lucrative option is polycarbonate diols, which are the precursor to polyurethanes that find wide applications in automotive and aerospace industries. The traditional process for the synthesis of polycarbonate diols involves using hazardous phosgene and alcohols in a strongly basic medium that generates significant amounts of salt. Alternately, CO2 can be reacted with alkanediols over CeO2 catalysts to produce polycarbonate diols. For the current work, a microwave-assisted atmospheric flow system was investigated for converting CO2 into polycarbonate diols and was compared against a thermal system. Various parameters such as type and ratio of solvents, and the addition of dehydrating agents for the microwave system was tested. The different CeO2 catalysts used were characterized using Brunauer–Emmett–Teller (BET) surface area analysis, Raman spectroscopy, X-ray diffraction, and temperature programmed desorption (NH3-TPD and CO2-TPD).

CO2 utilization↗
DOE OSTI2024

Four Years of Atmospheric Boundary Layer Height Retrievals Using COSMIC-2 Satellite Data

This work aimed to study the atmospheric boundary layer height (ABLH) from COSMIC-2 refractivity data, endeavoring to refine existing ABLH detection algorithms and scrutinize the resulting spatial and seasonal distributions. Through validation analyses involving different ground-based methodologies (involving data from lidar, ceilometer, microwave radiometers, and radiosondes), the optimal ABLH determination relied on identifying the lowest refractivity gradient negative peak with a magnitude at least $τ$% times the minimum refractivity gradient magnitude, where $τ$ is a fitting parameter representing the minimum peak strength relative to the absolute minimum refractivity gradient. Different $τ$ values were derived accounting for the moment of the day (daytime, nighttime, or sunrise/sunset) and the underlying surface (land or sea). Results show discernible relations between ABLH and various features, notably, the land cover and latitude. On average, ABLH is higher over oceans (≈1.5 km), but extreme values (maximums > 2.5 km, and minimums < 1 km) are reached over intertropical lands. Variability is generally subtle over oceans, whereas seasonality and daily evolution are pronounced over continents, with higher ABLHs during daytime and local wintertime (summertime) in intertropical (middle) latitudes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Sampling Size Optimization for Bioburden Density Estimation in Planetary Protection

Planetary protection (PP) is a discipline that focuses on minimizing the biological contamination of spacecraft to ensure compliance with international policy. Precise estimation of bioburden - the total number of microbes in or on spacecraft hardware – and the bioburden density are of utmost importance for PP. Such estimation is the way concordance with requirements is demonstrated, and it is critical for quantifying the potential risk of inadvertently contaminating other planetary bodies. Although a suite of molecular techniques have been used to thoroughly characterize and profile the microbiome of various cleanroom environments and spacecraft, the gold standard remains the physical enumeration of microbes via culturing of samples directly taken from spacecraft and associated surfaces. However, due to technical, budgetary, and programmatic constraints, only a manageable portion (around 10%) of the entire spacecraft surface is directly sampled with cotton swabs or wipes. To generate the bioburden current best estimate (CBE) for components not directly verifiable, the accepted approach is to apply a NASA-defined bioburden estimate based on the components’ manufacturing or assembly environment. This approach utilizes a prespecified bioburden density estimation that applies a maximum value across the total surface area of the specified component. For hardware components that underwent similar assembly processes, an implied bioburden is adopted for all components, based on a direct verification of a representative component within the same lot. Once all components have a CBE, the bioburden estimates are generated. In previous publication [ 1], we have shown that statistical risks quantifying the accuracy of the estimates for sampled, prespecified, and implied components can be derived and ranked. For mean squared error (MSE) function, the risks are available analytically and hence a cost function can be obtained to optimize the risks with respect to the sampling area and sampling cost. Since the sampling area and sampling cost are two complimentary variables, their sum will have a well-defined minimum. This paper presents the multivariate optimization of the integrated risk of an empirical Bayes estimator to determine the optimal sampling schedule for a given number of components. It is assumed that given a number of components, N, the bioburden density for each component can either be sampled, implied, or prespecified. The multivariate optimization searches through different options to sample, imply or prespecify the bioburden density for a component, and account for the component’s surface area and cost of sampling. The idea of the optimization is based on the observation that the statistical risk of using an estimator is a monotonically decreasing function of the sampled area. The larger the sampled area, the lower the risk of using the estimator as the estimator becomes more and more accurate as the sampling area increases. On the other hand, the cost of sampling is monotonically increasing as the sampled surface grows. This makes the risk and total cost of sampling complimentary variables which can be counterbalanced to achieve an optimal overall value with respect to the sampled surface. In this paper, the integrated risk has been used to quantify the accuracy of the estimator. This risk has been selected because it depends on neither the true value of the parameter nor on the collected data. The cost of each sample was also available to obtain the total cost of sampling of N components. The paper will present the results based on computer-simulated data as well as the data collected during the InSight mission. The computer-simulated data have N components with randomly generated total areas and each component assigned to one of the three categories according to the method of estimating of bioburden density: sampled, implied, or prespecified. The cost of sampling is also available. The cost of sampling is estimated based on a cost model provided by the planetary protection group at JPL. For this paper, the overall cost was assumed to be a linear function of exposure. The optimization process finds the allocation of the components to the three categories that minimizes the tradeoff between integrated risk and total cost. For the InSight data, a set of components is selected representing all three categories, and optimization is performed to determine if the performed allocation was optimal or if a better allocation could have been obtained. To the best of our knowledge, this work is the first attempt not only perform an accurate estimation of bioburden density but also do it in an optimal way.

97 - MATHEMATICS AND COMPUTING↗