Search NASA⌕ Search

SEARCH · Search NASA

Results for “variable density”

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 721 records · Page 40
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

In Quest of Calibration Density and Consistency in Hydrologic Modeling: Distributed Parameter Calibration against Streamflow Characteristics

Conventional basin-by-basin approaches to calibrate hydrologic models are limited to gauged basins and typically result in spatially discontinuous parameter fields. Moreover, the consequent low calibration density in space falls seriously behind the need from present-day applications like high resolution river hydrodynamic modeling. In this study we calibrated three key parameters of the Variable Infiltration Capacity (VIC) model at every 1/8° grid-cell using machine learning-based maps of four streamflow characteristics for the conterminous United States (CONUS), with a total of 52,663 grid-cells. This new calibration approach, as an alternative to parameter regionalization, applied to ungauged regions too. A key difference made here is that we tried to regionalize physical variables (streamflow characteristics) instead of model parameters whose behavior may often be less well understood. The resulting parameter fields no longer presented any spatial discontinuities and the patterns corresponded well with climate characteristics, such as aridity and runoff ratio. The calibrated parameters were evaluated against observed streamflow from 704/648 (calibration/validation period) small-to-medium-sized catchments used to derive the streamflow characteristics, 3941/3809 (calibration/validation period) small-to-medium-sized catchments not used to derive the streamflow characteristics as well as five large basins. Comparisons indicated marked improvements in bias and Nash-Sutcliffe efficiency. Model performance was still poor in arid and semiarid regions, which is mostly due to both model structural and forcing deficiencies. Although the performance gain was limited by the relative small number of parameters to calibrate, the study and results here served as a proof-of-concept for a new promising approach for fine-scale hydrologic model calibrations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2017

Variable range hopping electric and thermoelectric transport in anisotropic black phosphorus

Black phosphorus (BP) is a layered semiconductor with a high mobility of up to ~1000 cm 2 V -1 s -1 and a narrow bandgap of ~0.3 eV, and shows potential applications in thermoelectrics. In stark contrast to most other layered materials, electrical and thermoelectric properties in the basal plane of BP are highly anisotropic. In order to elucidate the mechanism for such anisotropy, we fabricated BP nanoribbons (~100 nm thick) along the armchair and zigzag directions, and measured the transport properties. It is found that both the electrical conductivity and Seebeck co efficient increase with temperature, a behavior contradictory to that of traditional semiconductors. The three-dimensional variable range hopping model is adopted to analyze this abnormal temperature dependency of electrical conductivity and Seebeck coefficient. Furthermore, the hopping transport of the BP nanoribbons, attributed to high density of trap states in the samples, provides a fundamental understanding of the anisotropic BP for potential thermoelectric applications.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Effects of macroparticle weighting in axisymmetric particle-in-cell Monte Carlo collision simulations

In this paper, the effects of macroparticle (MP) weighting on the plasma discharge, particularly near the centerline, are investigated using a two-dimensional axisymmetric particle-in-cell Monte Carlo collision (PIC/MCC) model. A variable MP weight according to the radial position of the MPs is employed to maintain sufficient number of MPs near the centerline of the plasma source. Finally, the plasma density obtained from the PIC/MCC simulations for low-pressure (25–100 mTorr) capacitively coupled plasmas is found to be artificially large when the MP weight near the centerline is not well resolved, demonstrating the need for particle convergence studies for axisymmetric PIC/MCC simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

A Hybrid Modulation Method for Single-Stage Soft-Switching Inverter Based on Series Resonant Converter

Pseudo-dc-link inverters are promising due to the single-stage configuration. The inverters can be isolated by high-frequency transformers, which also provide the freedom of step-up/down through winding turns ratio. Conventionally, a pseudo-dc-link inverter is based upon a pulsewidth modulation converter like flyback converter, but the accompanying switching losses can degrade efficiency and overall system performance. Alternatively, resonant converters (e.g., series resonant converter) seem to fit the need as they typically display soft-switching characteristics. However, the relevant research results are scarce because of two main reasons: high resonant current and limited soft-switching range. To overcome the difficulties, a hybrid modulation method is proposed in this article for a pseudo-dc-link inverter based on series resonant converter. In high-line region, the inverter operates with variable-frequency modulation, which leads to zero-voltage switching and relatively low-resonant current. In low-line region, the inverter operates with a short-pulse density modulation so that switching losses are greatly curtailed without excessive resonant current. The merits in both regions together contribute toward the high efficiency of the inverter. This article presents the operational principle of the proposed modulation as well as a comparison with other modulation methods. In addition, a 400 V-230 Vac, 2-kW prototype was built to verify the proposed modulation and the peak efficiency reaches 98.0%.

42 ENGINEERING↗
DOE OSTI2021

Kilonovae Across the Nuclear Physics Landscape: The Impact of Nuclear Physics Uncertainties on r-process-powered Emission

Merging neutron stars produce “kilonovae”—electromagnetic transients powered by the decay of unstable nuclei synthesized via rapid neutron capture (the r-process) in material that is gravitationally unbound during inspiral and coalescence. Kilonova emission, if accurately interpreted, can be used to characterize the masses and compositions of merger-driven outflows, helping to resolve a long-standing debate about the origins of r-process material in the Universe. We explore how the uncertain properties of nuclei involved in the r-process complicate the inference of outflow properties from kilonova observations. Using r-process simulations, we show how nuclear physics uncertainties impact predictions of radioactive heating and element synthesis. For a set of models that span a large range in both predicted heating and final abundances, we carry out detailed numerical calculations of decay product thermalization and radiation transport in a kilonova ejecta with a fixed mass and density profile. The light curves associated with our models exhibit great diversity in their luminosities, with peak brightness varying by more than an order of magnitude. We also find variability in the shape of the kilonova light curves and their color, which in some cases runs counter to the expectation that increasing levels of lanthanide and/or actinide enrichment will be correlated with longer, dimmer, redder emission.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2017

Seamless Staircase Electrical Contact to Semiconducting Graphene Nanoribbons

Electrical contact to low-dimensional (low-D) materials is a key to their electronic applications. Traditional metal contacts to low-D semiconductors typically create gap states that can pin the Fermi level (E F ). Yet, low-D metals possessing a limited density of states at E F can enable gate-tunable work functions and contact barriers. Furthermore, a seamless contact with native bonds at the interface, without localized interfacial states, can serve as an optimal electrode. To realize such a seamless contact, one needs to develop atomically precise heterojunctions from the atom up. In this work, we demonstrate an all-carbon staircase contact to ultranarrow armchair graphene nanoribbons (aGNRs). The coherent heterostructures of width-variable aGNRs, consisting of 7, 14, 21, and up to 56 carbon atoms across the width, are synthesized by a surface-assisted self-assembly process with a single molecular precursor. The aGNRs exhibit characteristic vibrational modes in Raman spectroscopy. A combined scanning tunneling microscopy and density functional theory study reveals the native covalent-bond nature and quasi-metallic contact characteristics of the interfaces. Our electronic measurements of such seamless GNR staircase constitute a promising first step toward making low resistance contacts.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2024

Laser pulse induced second- and third-harmonic generation of gold nanorods with real-time time-dependent density functional tight binding (RT-TDDFTB) method

In this paper, we investigate second- and third-harmonic generation processes in Au nanorod systems using the real-time time-dependent density functional tight binding method. Our study focuses on the computation of nonlinear signals based on the time dependent dipole response induced by linearly polarized laser pulses interacting with nanoparticles. We systematically explore the influence of various laser parameters, including pump intensity, duration, frequency, and polarization directions, on harmonic generation. We demonstrate all the results using Au nanorod dimer systems arranged in end-to-end configurations, and disrupting the spatial symmetry of regular single nanorod systems is crucial for second-harmonic generation processes. Furthermore, we study the impact of nanorod lengths, which lead to variable plasmon energies, on harmonic generation, and estimates of polarizabilities and hyper-polarizabilities are provided.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2018

Optical Tomograms of Multiple-Photon-Added Gaussian States via the Intermediate State Representation Theory

Optical tomogram of a quantum state can serve as an alternative to its density matrix since it contains all the information on this state. In this paper we use the explicit expression of the intermediate state |q, f, g〉(that is, an eigenvector of the coordinate and momentum operator fQ + gP) and the Radon transform between the Wigner operator Δ(p, q) and the pure-state density operator |q, f, g〉〈q, f, g| to calculate the specific tomograms for multiple-photon-added coherent states (MPACSs), thermal states (MPATSs) and displaced thermal states (MPADTSs). Their analytical tomograms are actually the finite-summation representations related to single-variable Hermite polynomials. Besides, the numerical results indicate that the larger photon-addition number can lead to the multi-peak tomographic distributions for these three states, and the measurable tomograms of MPACSs and MPADTSs exhibit π-periodic features but MPATSs can’t.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

A new spatially distributed added value index for regional climate models: the EURO-CORDEX and the CORDEX-CORE highest resolution ensembles

The added value of using regional climate models (RCMs) to downscale data from general circulation models (GCMs) has often been questioned and researched. Although several studies have used different methods to identify (and in some cases quantify) the added value, there is still a need to find a general metric that quantifies the added value of any variable. This paper builds on past studies to propose a new metric of added value in the simulation of present-day climate which measures the difference in the probability density functions (PDFs) at each grid-cell between a model and an observation source, and then compares the results of the RCM and GCM in order to spatially compute the added value index. The same method is also adapted to quantify the climate change downscaling signal in a way that is consistent with the present-day metric. These new metrics are tested on the daily precipitation output from the EURO-CORDEX and CORDEX-CORE projection ensembles and reveal an overall positive added value of RCMs, especially at the tail-end of the distribution. Higher added value is obtained in areas of complex topography and coast-lines, as well as in tropical regions. Areas with large added value in present-day climate are consistent with areas of significant climate change downscaling signal in the RCP 8.5 far future simulations, and when the analysis is repeated at a low-resolution. The use of different resolution observations shows that the added value tends to decrease when models are compared to low-resolution observation datasets.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Clean quartz matters for cosmogenic nuclide analyses: An exploration of the importance of sample purity using the CRONUS-N reference material

We report reference materials are key for assessing inter-laboratory variability and measurement quality, and for placing analytical uncertainty bounds on sample analyses. Here, we investigate four years of data resulting from repeated processing of the CRONUS-N reference material for cosmogenic 10 Be and 26 Al analyses. At University of Vermont, we prepared a CRONUS-N aliquot with most of our sample batches from 2013 to 2017; these reference material samples were then distributed to four different accelerator mass spectrometry facilities, yielding 73 10 Be analyses and 58 26 Al analyses. We determine CRONUS-N 10 Be concentrations of (2.26 ± 0.14) x 10 5 atoms g -1 (n = 73, mean, 1 SD) and 26 Al concentrations of (1.00 ± 0.08) x 10 6 atoms g -1 (n = 58, mean, 1 SD). We find a reproducibility of 6.3% for 10 Be and 7.7% for 26 Al (relative standard deviations). We also document highly variable 27 Al and Mg concentrations and a 10 Be dispersion twice as large as the mean AMS analytic uncertainty. Analyses of the CRONUS-N material with and without density separation demonstrate that non-quartz minerals are present in the material and have a large impact on measured concentrations of 27 Al, 10 Be, and impurities; these non-quartz minerals represent only a very small portion of the total mass (0.6–0.8%) but have a disproportionally large effect on the resulting data. Our results highlight the importance of completely removing all non-quartz mineral phases from samples prior to Be/Al extraction for the determination of in situ cosmogenic 10 Be and 26 Al concentrations.

58 GEOSCIENCES↗
DOE OSTI2020

Emergence of unprecedented climate change in projected future precipitation

The future time of emergence when precipitation changes due to anthropogenic influences begins to continuously exceed the previous maximum value is defined as the ‘tipping year’ Historical experiments and future experiments simulated by state-of-the-art climate models were utilized. A total of 510,000 time series from year 1856 to 2095 were generated by sampling the natural internal variability in precipitation. The time evolutions of internal variability in the whole time period were estimated from the combination of past and future experiments with preindustrial control experiments. A large ensemble size enabled an estimation of the probability density function of the tipping year at each grid point, providing precise information on the uncertainty of the projection. The tipping year of average precipitation emerges earlier in high latitudes than in lower latitudes. In some regions in lower latitudes and mid-latitudes, the tipping year of intense precipitation emerges faster than that of average precipitation. The tipping years of average and intense precipitation are earlier for higher anthropogenic forcing scenarios than for lower scenarios. The global average of the tipping year for intense precipitation might be attributed to the enhancement of the thermodynamic effect (moisture) rather than the dynamic effect (vertical motion).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Emerging investigator series: kinetics of diopside reactivity for carbon mineralization in mafic–ultramafic rocks

The ongoing use of fossil fuels to supply modern energy demands has necessitated research on combating carbon dioxide (CO 2 ) emissions and climate change. Carbon storage via mineral trapping in basalt and related rocks is a promising strategy. However, mineralization rates depend on the variable minerology that makes up these rock formations. Diopside (CaMgSi 2 O 6 ) is a common pyroxene mineral in ultramafic and mafic rocks including basalt, but relatively little work has been done to understand its carbon mineralization kinetics using hydrated supercritical CO 2 , which induces the formation of reactive nanoscale interfacial water films. Here, in situ XRD experiments at 50–110 °C and 90 bar indicate that diopside transforms into a myriad of Mg/Ca carbonates, including huntite [Mg 3 Ca(CO 3 ) 4 ] and very high magnesium calcite (VHMC, i.e., protodolomite). Through ex situ characterization, we were able to constrain reaction pathways for the dissolution–precipitation transformation process including metastable intermediate precipitates. Experiments performed at variable temperatures enabled Avrami-derived rate constants and an apparent activation energy of 97 ± 16 kJ mol –1 , implying the dissolution of diopside is the rate-limiting step. Density functional theory (DFT) calculations, used to gain molecular insight into the surface stability of the diopside during dissolution, suggest that exposed calcium cations are susceptible to dissolution when put in contact with water given their coordination environment. The collective results point to the high CO 2 mineralization potential of diopside in basalts, which could help guide parameterization of reactive transport models needed to design and permit commercial-scale subsurface carbon storage operations.

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

Action principles for relativistic extended magnetohydrodynamics: A unified theory of magnetofluid models

Two types of Eulerian action principles for relativistic extended magnetohydrodynamics (MHD) are formulated. With the first, the action is extremized under the constraints of density, entropy, and Lagrangian label conservation, which leads to a Clebsch representation for a generalized momentum and a generalized vector potential. The second action arises upon transformation to physical field variables, giving rise to a covariant bracket action principle, i.e., a variational principle in which constrained variations are generated by a degenerate Poisson bracket. Upon taking appropriate limits, the action principles lead to relativistic Hall MHD and well-known relativistic ideal MHD. For the first time, the Hamiltonian formulation of relativistic Hall MHD with electron thermal inertia (akin to Comisso et al. , Phys. Rev. Lett. 113 , 045001 (2014) for the electron–positron plasma) is introduced. This thermal inertia effect allows for violation of the frozen-in magnetic flux condition in marked contrast to nonrelativistic Hall MHD that does satisfy the frozen-in condition. We also find the violation of the frozen-in condition is accompanied by freezing-in of an alternative flux determined by a generalized vector potential. Finally, we derive a more general 3 + 1 Poisson bracket for nonrelativistic extended MHD, one that does not assume smallness of the electron ion mass ratio.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

Molecular-beam epitaxy of InGaAs/InAlAs/AlAs structures for heterobarrier varactors

The molecular-beam epitaxy of InGaAs/InAlAs/AlAs structures for heterobarrier varactors is studied and optimized. The choice of the substrate-holder temperature, growth rate and III/V ratio in the synthesis of individual heterostructure regions, the thickness of AlAs inserts and barrier-layer quality are critical parameters to achieve the optimal characteristics of heterobarrier varactors. The proposed triple-barrier structures of heterobarrier varactors with thin InGaAs strained layers immediately adjacent to an InAlAs/AlAs/InAlAs heterobarrier, mismatched with respect to the InP lattice constant at an AlAs insert thickness of 2.5 nm, provides a leakage current density at the level of the best values for heterobarrier varactor structures with 12 barriers and an insert thickness of 3 nm.

36 MATERIALS SCIENCE↗
DOE OSTI2017

On the Preconditioning of a Newton-Krylov Solver for a High-Order reconstructed Discontinuous Galerkin Discretization of All-Speed Compressible Flow with Phase Change for Application in Laser-Based Additive Manufacturing

This dissertation focuses on the development of a fully-implicit, high-order compressible ow solver with phase change. The work is motivated by laser-induced phase change applications, particularly by the need to develop large-scale multi-physics simulations of the selective laser melting (SLM) process in metal additive manufacturing (3D printing). Simulations of the SLM process require precise tracking of multi-material solid-liquid-gas interfaces, due to laser-induced melting/ solidi cation and evaporation/condensation of metal powder in an ambient gas. These rapid density variations and phase change processes tightly couple the governing equations, requiring a fully compressible framework to robustly capture the rapid density variations of the ambient gas and the melting/evaporation of the metal powder. For non-isothermal phase change, the velocity is gradually suppressed through the mushy region by a variable viscosity and Darcy source term model. The governing equations are discretized up to 4th-order accuracy with our reconstructed Discontinuous Galerkin spatial discretization scheme and up to 5th-order accuracy with L-stable fully implicit time discretization schemes (BDF2 and ESDIRK3-5). The resulting set of non-linear equations is solved using a robust Newton-Krylov method, with the Jacobian-free version of the GMRES solver for linear iterations. Due to the sti nes associated with the acoustic waves and thermal and viscous/material strength e ects, preconditioning the GMRES solver is essential. A robust and scalable approximate block factorization preconditioner was developed, which utilizes the velocity-pressure (vP) and velocity-temperature (vT) Schur complement systems. This multigrid block reduction preconditioning technique converges for high CFL/Fourier numbers and exhibits excellent parallel and algorithmic scalability on classic benchmark problems in uid dynamics (lid-driven cavity ow and natural convection heat transfer) as well as for laser-induced phase change problems in 2D and 3D.

42 ENGINEERING↗
DOE OSTI2016

UNIVERSE IN A BLACK HOLE IN EINSTEIN–CARTAN GRAVITY

The conservation law for the angular momentum in curved spacetime, consistent with relativistic quantum mechanics, requires that the antisymmetric part of the affine connection (torsion tensor) is a variable in the principle of least action. The coupling between the spin of elementary particles and torsion in the Einstein–Cartan theory of gravity generates gravitational repulsion at extremely high densities in fermionic matter, approximated as a spin fluid, and thus avoids the formation of singularities in black holes. The collapsing matter in a black hole should therefore bounce at a finite density and then expand into a new region of space on the other side of the event horizon, which may be regarded as a nonsingular, closed universe. We show that quantum particle production caused by an extremely high curvature near a bounce can create enormous amounts of matter, produce entropy, and generate a finite period of exponential expansion (inflation) of this universe. This scenario can thus explain inflation without a scalar field and reheating. We show that, depending on the particle production rate, such a universe may undergo several nonsingular bounces until it has enough matter to reach a size at which the cosmological constant starts cosmic acceleration. The last bounce can be regarded as the big bang of this universe.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Reference range of liver corrected T1 values in a population at low risk for fatty liver disease—a UK Biobank sub-study, with an appendix of interesting cases

Purpose: Corrected T1 (cT1) value is a novel MRI-based quantitative metric for assessing a composite of liver inflammation and fibrosis. It has been shown to distinguish between non-alcoholic fatty liver disease (NAFL) and non-alcoholic steatohepatitis. However, these studies were conducted in patients at high risk for liver disease. This study establishes the normal reference range of cT1 values for a large UK population, and assesses interactions of age and gender. Methods: MR data were acquired on a 1.5 T system as part of the UK Biobank Imaging Enhancement study. Measures for Proton Density Fat Fraction and cT1 were calculated from the MRI data using a multiparametric MRI software application. Data that did not meet quality criteria were excluded from further analysis. Inter and intra-reader variability was estimated in a set of data. A cohort at low risk for NAFL was identified by excluding individuals with BMI ≥ 25 kg/m{sup 2} and PDFF ≥ 5%. Of the 2816 participants with data of suitable quality, 1037 (37%) were classified as at low risk. Results: The cT1 values in the low-risk population ranged from 573 to 852 ms with a median of 666 ms and interquartile range from 643 to 694 ms. Iron correction of T1 was necessary in 36.5% of this reference population. Age and gender had minimal effect on cT1 values. Conclusion: The majority of cT1 values are tightly clustered in a population at low risk for NAFL, suggesting it has the potential to serve as a new quantitative imaging biomarker for studies of liver health and disease.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
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↗