Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “surrogates”

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 433 records · Page 24

Stress Intensity Factors for Layered Pressure Vessel Inner Layer Through Cracks

A difficulty encountered when performing Fitness-for- Service assessments for layered pressure vessels (LPVs) is the lack of stress intensity factor solution in literature that produce accurate results for inner layer longitudinal through cracks. Using surrogate solutions such as a through crack in a plate or cylinder produce results that can be overly conservative especially for longer cracks. This is largely due to the ability of a layered pressure vessel to redistribute hoop load to other layers, the restricted radial movement of the cracked layer, and the friction forces applied in the cracked region. To understand this problem, a parametric finite element model (FEM) generator was developed that is capable of producing layered pressure vessel models with inner layer through cracks. The results from the FEMs were used to create a dataset of inner layer through crack stress intensity factors (KI) for layered pressure vessels corresponding to variations of internal pressure, radius, layer thicknesses, friction factor, and crack length. The elastic modulus of the material also has an effect on KI but, for this dataset, the elastic modulus was fixed at the typical value for steel - 29,500 ksi (203 GPa). Finally, a non-dimensional model was developed and calibrated using the dataset. This allows KI to be calculated without the need of a FEM using a closed-form equation. The results of the closed-form solution were then compared to FEM results showing accuracy was generally within 10%.

Hobbs, Joel R.↗

Investigation of Reduced-Order Modeling for Aircraft Stability and Control Prediction

High-fidelity computational fluid dynamics tools offer the potential to approximate increments for ground-to-flight scaling effects, as well as to augment the dynamic damping derivative data for motion-based flight simulators. Unfortunately, the computational expense is currently prohibitive for populating a complete simulator database. This work investigates an existing surrogate-based, indicial response reduced-order model methodology as a means to efficiently augment a flight simulator database with high-fidelity nonlinear aerodynamic damping derivatives. Creation of the reduced-order model is based on the superposition integrals of the step response with the derivative of its corresponding input signal. Step responses are calculated using a computational grid motion approach that separates the effects of angle of attack and sideslip angle from angular rates, and rates from angle of attack and sideslip. It is demonstrated that the transients produced during the start of a forced-oscillation motion are captured by the reduced-order model to the level of fidelity of a comparable computational solution. Aerodynamic coefficients computed within minutes by the reduced-order model for an aircraft undergoing an 18-second half Lazy-8 maneuver and a 25-second Immelmann turn maneuver are compared with those from full computational flight solutions that required days to complete. Finally, a cost-benefit assessment is included that demonstrates a compelling advantage for this approach. d for maneuvering, flexible vehicles.

Frink, Neal T.↗

Simulating Fleet Noise for Notional UAM Vehicles and Operations in New York

This paper presents the results of systems-level simulations using Metrosim that were conducted for notional Urban Air Mobility (UAM)-style vehicles analyzed for two different scenarios for New York (NY). UAM is an aviation industry term for passenger or cargo-carrying air transportation services, which are often automated, operating in an urban/city environment. UAM-style vehicles are expected to use vertical takeoff and landing with fixed wing cruise flight. Metrosim is a metroplex-wide route and airport planning tool that can also be used in standalone mode as a simulation tool. The scenarios described and reported in this paper were used to evaluate a fleet noise prediction capability for this tool. The work was a collaborative effort between the National Aeronautics and Space Administration (NASA), Intelligent Automation, Inc (IAI), and the Port Authority of New York and New Jersey (PANYNJ). One scenario was designed to represent an expanded air-taxi operation from existing helipads around Manhattan to the major New York airports. The other case represented a farther term vision case with commuters using personal air vehicles to hub locations just outside New York, with an air-taxi service running frequent connector trips to a few key locations inside Manhattan. For both scenarios, the trajectories created for the entire fleet were passed to the Aircraft Environmental Design Tool (AEDT) to generate Day-Night Level (DNL) noise contours for inspection. Without data for actual UAM vehicles available, surrogate AEDT empirical Noise-Power-Distance (NPD) tables used a similar sized current day helicopter as the Baseline, and a version of that same data linearly scaled as a first guess at possible UAM noise data. Details are provided for each of the two scenario configurations, and the output noise contours are presented for the Baseline and reduced noise DNL cases.

Glaab, Patricia↗

Parametric Cost Estimates of Four 20 Ton Payload Mars EDL Vehicle Concepts

Cost is one of the biggest obstacles to sending humans to Mars. Many studies have examined various vehicles and/or technologies which could be used to send humans to Mars. However, nearly all of these studies have either ignored cost or have used the mass of the vehicles as a surrogate for cost. This work presents parametric cost estimates of four Mars Entry, Descent, and Landing vehicle concepts capable of landing 20 tons of payload on the surface of Mars. The concepts are the Co-Optimization Blunt-body Re-entry Analysis-Mid lift to drag ratio Rigid Vehicle (Cobra–MRV), Adaptive Deployable Entry and Placement Technology (ADEPT), Hypersonic Inflatable Aerodynamic Decelerator (HIAD), and Capsule concepts. Two variants of the HIAD concept are studied, one with a single HIAD heatshield and one with two HIADs. The study presents the development and production cost of each vehicle and the total cost to build five flight units and fifteen flight units. Five flight units is representative of a single mission, or a "Flags and Footprints" style campaign, whereas fifteen flight units is representative of multiple missions going to multiple sites on Mars. All costs presented in this work are normalized to the average of the first unit production cost of the lander concepts. There is a large amount of uncertainty both in the design of the vehicles and the cost of the vehicles and it is impossible to claim with certainty which vehicle will cost the least. However, it is likely that the total cost to produce five flight units will be in the following ranking from least costly to most costly, Single HIAD, Capsule, Dual HIAD, Cobra, ADEPT. The ranking for fifteen units is similar but with the positions of the Dual HIAD and Cobra reversed. This study only considers the cost of the entry vehicle and does not take into account the effect of the entry vehicle’s mass on the cost of the in-space propulsion stage or launch vehicle.

Paul D Friz↗

Toward Transition Modeling in a Hypersonic Boundary Layer at Flight Conditions

An accurate physics-based transition prediction method integrated with computational fluid dynamics (CFD) solvers is pursued for hypersonic boundary layer flows over slender hypersonic vehicles at flight conditions. The geometry and flow conditions are selected to match relevant trajectory locations from the ascent phase of the HIFiRE-1 flight experiment, namely, a 7-degree half-angle cone with 2.5 mm nose radius, freestream Mach numbers in the range of 3.8 – 5.5 and freestream unit Reynolds numbers in the range of 3.3 × 10(exp 6) – 21.4 × 10(exp 6) m(exp -1). Earlier research had shown that the onset of transition during the HIFiRE-1 flight experiment correlated with an amplification factor of N ≈ 13.5 for the planar Mack modes. However, to incorporate the N-factor correlations into a CFD code, we investigate surrogate models for disturbance amplification that avoid the direct computation of stability characteristics. A commonly used approach for low-speed flows is based on an a priori database of stability characteristics for locally similar profiles. However, the results presented in this paper demonstrate that the application of this approach to hypersonic boundary layers over blunt spherical nose-tip cones leads to large, unacceptable errors in the predictions of amplification factors, mainly due to its failure in accounting for the effects of the entropy layer on the boundary-layer profiles along the length of the model. We propose and demonstrate an alternate approach that employs the stability computations for a canonical set of blunt cone configurations to train a physics-informed convolutional neural network model that is shown to provide substantially improved transition predictions for hypersonic flow configurations with entropy-layer effects. Furthermore, the excellent performance of the neural network model is also confirmed for cone configurations with nose radius and half-angle values that do not correspond to those used to build the database. Finally, the convolutional neural network model is shown to outperform the linear stability calculations for underresolved basic states.

Pedro Paredes↗

Analytical and Experimental Demonstration of an Alternate Mixing Performance Metric for High-Speed Fuel Mixing Studies

To experimentally assess the fuel/air mixing performance of high-speed fuel injectors, one-dimensional metrics that quantify the degree of mixing completeness downstream of the fuel injection location are required. The most accurate assessment of mixing performance is achieved with the mixing efficiency parameter. In order to experimentally determine the mixing efficiency parameter, the spatial distributions of both mass flux and fuel mass fraction must be measured. In-stream gas sampling techniques are commonly used to measure the fuel mass fraction distribution; however, the mass flux distribution is not easily determined because it requires the measurement of three independent aerothermodynamic variables in addition to the gas composition. Therefore, to experimentally determine the mixing efficiency parameter, the spatial distributions of four independent properties must be measured, with each property generally requiring its own unique probe. Because of this difficulty, it is commonly assumed that alternate metrics, which rely solely on the fuel distribution, are good indicators of mixing performance. However, since these alternate metrics do not provide a mass flux-weighted measure of mixing completeness, they can lead to incorrect conclusions being drawn about the mixing performance of the studied fuel injector configuration. Recognizing this shortcoming, this work proposes two new alternate mixing performance metrics that are easier to obtain than the mixing efficiency parameter. The analytical development of the new metrics, as well as their application to relevant CFD and experimental data of high-speed fuel injector configurations, is presented in this work. For two different experiments, the new metrics are shown to provide an excellent representation of the true mass flux-weighted mixing performance, unlike the traditionally-used alternatives. The results presented herein suggest that the new metrics can serve as accurate surrogates for mixing efficiency in future high-speed fuel/air mixing studies.

Cody R Ground↗

Tuning Neural Network Models for Improved Prediction of Boundary Layer Transition

Boundary layer transition can strongly impact flight vehicle performance as it influences surface skin friction and aerodynamic heating, making accurate transition prediction a key to designing next generation aircraft. Artificial neural networks (ANNs) have shown promise toward predicting laminar-turbulent transition based on linear stability correlations. The computational efficiency of ANNs and the substantially reduced user involvement in relation to direct computations based on the linear stability theory (LST) makes them an attractive methodology for integrating the LST based correlations in computational fluid dynamics codes. Tollmien-Schlichting (TS) waves correspond to the dominant transition mechanism in 2D or weakly 3D subsonic boundary layers, such as those encountered in general aviation applications. Improvements to neural network model accuracy in predicting the amplification rates of TS instability waves have been investigated by leveraging recent machine learning developments in conjunction with surrogate optimization techniques and via suitable augmentation of the data used to train the networks. The optimized models trained on the modified dataset reduced the average transition location errors on different airfoils at several flow conditions by 51% of the original manually-tuned network’s errors on the same flow cases. The actual transition locations were derived from the Langley Stability and Transition Analysis Code (LASTRAC).

Machine Learning↗

Improved Benchmarking of Cohesive Elements in Abaqus Standard for Predicting Disbond and Delamination in Composite Structures

Traditional approaches for aircraft certification require the assumption of an initial flaw condition, either represented as barely visible impact damage (BVID) or through inclusion of a Teflon insert to serve as surrogate damage. Based on the initial composite damage state, the structure must be shown to demonstrate structural durability and damage tolerance (DaDT) according to the following criteria: a. Damage displays no detrimental growth under cyclic loading b. The structure is able to sustain design limit load (DLL) Currently, the only available manner for validating structural performance is through test. Since damage can occur over a wide variety of areas within a structure, this approach has proven to be increasingly expensive and time consuming for composite airframes and acreage structure within the design-test-certification building block. A further complicating factor is the requirement to accurately capture the most critical damage morphologies as a starting condition. To understand the severity of the damage, it is either required to experimentally determine the most critical areas at tremendous expense or rely on legacy data of similar structural testing, which limits design space expansion. A preferred solution is to use advanced analysis to provide improved understanding of load margins for critical locations based on a wide variety of potential starting damage conditions. The standard industry approach for DaDT certification adheres to the use of the traditional virtual crack closure technique (VCCT) method. VCCT is generally a preferred method because it conforms to the current certification principles of damage from a known flaw, and when used correctly, can be effective at predicting delamination propagation under static and cyclic loading. The VCCT method requires the inclusion of an initial flaw in the finite element (FE) model requiring a-priori knowledge of the flaw location. This in turn requires a plethora of analysis cases to be examined to cover a reasonable span of potential damage states. Additionally, the VCCT approach requires node-to-node connectivity rendering it incompatible with the best practices and approaches for using continuum damage mechanics (CDM) based progressive damage and failure analysis (PDFA) tools within a typical FE solver. Alternatives to VCCT have emerged in the form of cohesive elements which utilize the cohesive zone model (CZM). Unlike VCCT which models linear elastic fracture mechanics, cohesive elements couples continuum and fracture based responses through the use of bilinear traction separation laws. These laws are defined based on a penalty stiffness, a cohesive strength, and a strain energy release rate. The approach can be mesh regularized with native cohesive elements within many FE solvers such as Abaqus and LS-DYNA. In Phase I of the NASA Advanced Composites Consortium (ACC) post-buckled stiffened panel with BVID, Strength and Life [1], the performance of cohesive elements were benchmarked in comparison to VCCT and LEFM solutions and showed good agreement using Abaqus explicit [2]. To realize savings on current and future programs, it is still necessary to close technical gaps related to the use of cohesive elements with Abaqus Standard. Within a program environment, standard finite element analysis is the preferred analytical capability for quasi-static loading as it eliminates uncertainty due to oscillatory behavior commonly seen with explicit analysis. This oscillatory behavior creates difficulties in writing margins of safety based on the analysis. The use of negative tangent stiffness material models complicates convergence which typically requires the use of numerical controls such as viscous damping to overcome. To date, there has not been a comprehensive study on how to establish best practices for cohesive element convergence for predictive capability within the Abaqus implicit solver. In pursuit of these goals, under the NASA ACC program, several numerical benchmark problems were proposed including pure mode I (double cantilevered beam – DCB), pure mode II (end notch flexure – ENF), and symmetric/unsymmetric evolving mixed mode (single leg bend – SLB). This paper focuses on the use of cohesive elements to model the delamination through the use of CZM. Specifically, finite element models for the DCB, ENF, symmetric SLB, and unsymmetric SLB, are developed and various solution controls for convergence are studied to develop a best practice. Once the best practice has been developed, the predictive capability of the objective CZM model is used to analyze the hat pull-off strength of a standard hat stiffened configuration under various loading conditions.

Abaqus↗

Fast and Precise Trajectory Simulation for Entry, Descent, and Landing Using A Multi-Model Monte Carlo Approach

Predicting landing radius and other quantities of interest (QoI) for entry, descent, and landing (EDL) applications requires a viable uncertainty propagation method for quantifying the impact of uncertainties in wind pattern variations, atmospheric uncertainties, etc. While standard MC simulation is the defacto standard for providing robust and unbiased predictions,it is often infeasible for expensive, high-fidelity EDL models. Low-fidelity models are commonly constructed to replace the high-fidelity model in MC simulation for computational speedup,but at the expense of accuracy and unbiasedness. Emerging multi-model MC methods are bridging this gap by combining predictions from two or more models of varying fidelity and computational cost for efficient and unbiased uncertainty propagation. This work explores the use of multi-model MC for increasing the speed and precision of trajectory simulation for EDL. It is shown that combining a high-fidelity EDL model with low-fidelity models (e.g,data-driven, reduced physics) yields substantial computational speedup versus standard MCwith only the high-fidelity model. Moreover, the unbiasedness of multi-model MC predictions is highlighted by showing increased accuracy versus an approach that leverages a low-fidelity surrogate model alone.

James E. Warner↗

Development and Validation of Generic Maneuvering Flight Noise Abatement Guidance for Helicopters

An extensive flight test campaign has been conducted to look into developing actionable advice for pilots of today’s vehicles to reduce their acoustic footprints. Ten distinct vehicles were tested at three different test ranges, with nine of the vehicles’ data being documented here. Twelve pairs of turning conditions were tested to determine their effect on blade-vortex interaction noise. Each turning flight condition was evaluated using the peak A-weighted, band-limited (50 Hz - 2500 Hz), sound pressure level measured throughout the maneuver. This metric was a surrogate for blade-vortex interaction noise, and the difference between the peak values of each turning pair was investigated. That peak value difference was subsequently corrected by the offset from the intended vehicle altitude at turn initiation from the actual altitude at initiation. The corrected amplitudes were investigated and grouped into six validated actionable guidance principles that can be given to pilots to immediately reduce their acoustic footprint during operations.

Flight Test↗

Parametric Studies of Human Mars Entry, Descent, and Landing Systems

This paper will present a parametric analysis for entry, descent, and landing (EDL) concepts, enabling rapid systems assessment and tradespace exploration. The entry system uses a hypersonic inflatable aerodynamic decelerator (HIAD) technology. The baseline system includes elements for Mars aerocapture (AC) and EDL segments of the mission. The mission concept of operations (ConOps) begins at Mars arrival in a polar inclination. After performing an aerocapture (AC) maneuver into a pre-defined Mars parking orbit, the AC HIAD is jettisoned. The rest of thevehicle stays in the parking orbit for up to one year. The EDL sequence starts with a deorbit burn at the apoapsis of the parking orbit. After hypersonic entry with the EDL HIAD, the entry system uses a supersonic retropropulsion maneuver to slow the vehicle for the descent and landing segments of the mission. The vehicle will maintain a constant velocity of 2.5 m/s for 5 seconds prior to landing.The system includes a Mars Ascent Vehicle (MAV), Mars Descent Module (MDM), and two HIADs. The MDM includes a primary structure, tanks, engines, and radiators. The primary structure is an aluminum-lithium (Al-Li)cruciform design similar to the structural design of the Apollo Lunar Module. The cruciform planform layout results in four outer bays, with adequate volume in the corners between outer bays to package four landing gears. The central bay is reserved for packaging the MAV and the recessed MAV engines. Two of the outer bays accommodate main propellant tanks, with one LOX and one CH4 tank in each bay. The two remaining outer bays each house four rocket engine systems and associated support structure. Each HIAD comprises an inflatable structure, flexible thermal protection system, gas, and gas generators. The HIAD design used in this study is a stacked-toroid concept with pairing loop straps and radial/chevron straps. The baseline system lands a 20-t payload on the Mars surface. It is assumed the vehicle arrives at 6.2 km/s relative velocity at 90° inclination and is captured to a 1-Sol parking orbit.The baseline design includes many assumptions such as margins, arrival state, ConOps options, parking orbit, physical dimensions, propellant options, and technology concepts. The impact of these parameters are quantified through systems-level sensitivity analyses, which capture the global impact—not at a component level—but at the systems level. The systems-level sensitivities expose major design drivers and importance of each assumption for a design.Through tradespace exploration, a wide range of systems parameters are examined and compared for several feasible design options. Studies have been completed for the following input parameters: payload mass, propellant options, AC/EDL ballistic coefficient, lander thrust to weight ratio (T/W) (surrogate for the maximum EDL g’s),engine specific impulse (Isp), parking orbit, and inclination.The final paper will present and discuss the parametric approach used in the study. It will also include the results of recent systems analyses, sensitivity analyses, and tradespace exploration

Jamshid Samareh↗

Dose, LET, time and strain dependence of radiation-induced 53BP1 foci in 15 mouse strains ex vivo and associations to in vivo radiation susceptibility

We present a comparative analysis on the repair of radiation-induced DNA damage ex vivo in 15 strains of mice, including 5 inbred reference strains and 10 collaborative-cross strains, of both sexes. Non-immortalized primary skin fibroblasts derived from 76 mice were subjected to both low- and high-LET radiation (0.1, 1 and 4 Gy of X rays; 1.1 and 3 particles/100μm2 of 350 MeV/n 40Ar and 600 MeV/n 56Fe). Automated image quantification of 53BP1 radiation-induced foci (RIF) during the first 4-48 h post-irradiation was performed as a function of dose and LET. Similarly to what we had previously reported for immortalized human cell lines [1], we observed a saturation of RIF number with dose at 4h post-irradiation, with more RIF/Gy for lower LET (X rays and 40Ar) compared to 56Fe. However at later time points (24h and above), the trend was inverted with more RIF/Gy for higher LET. Our data suggest that multiple DSBs cluster into RIF: as the linear density of DSBs increases with LET, so does the probability of having more DSBs per RIF, which makes it more difficult for cells to fully resolve high-LET-induced RIF, explaining the hypersensitivity to high-LET radiation despite a low number of RIF. Taking into account the amount of clustering at a given dose and LET, but also the kinetics of DNA damage repair, we introduced a novel mathematical formalism to evaluate the number of remaining RIF over time. We showed that the newly introduced kinetic metrics can be used as surrogate biomarkers for in vivo radiation toxicity, with potential applications in radiotherapy and human space exploration. In particular, we observed an association between the repairable fraction of RIF measured in vitro and survival levels of immune cells collected from irradiated mice. Moreover, the speed of DNA damage repair correlated with spontaneous cancer incidence data collected from the Mouse Tumor Biology database, suggesting a relationship between the efficiency of DSB repair after irradiation and cancer risk. In addition to the efficacy of repair and persistent RIF levels, even the amount of spontaneous foci without irradiation was shown to be strain dependent, indicating that these phenotypes are at least partially driven by genetics, and supporting their potential as indicators of individual radiation sensitivity. [1] Neumaier, T., et al., PNAS, 2012 (8) 109:443

Radiation, DNA damage, repair kinetics↗

NASA and Blue Origin Collaborative Assessment of Precision Landing Algorithms and Computing

NASA’s Safe and Precise Landing Integrated Capabilities Evolution (SPLICE) project is developing sensor, algorithm, and compute technologies for precision landing and hazard avoidance. These technologies are being tested as an integrated Precision Landing and Hazard Avoidance (PL&HA) system on Blue Origin’s New Shephard suborbital vehicle. A key goal for the computing element of this technology development is to characterize the performance of the SPLICE software workloads on the project’s Descent and Landing Computer (DLC). The DLC is a multi-core processor designed as a surrogate for NASA’s High-Performance Space Computer (HPSC). Measurements of the SPLICE workload performance on the DLC provides NASA insight on how PL&HA capabilities will perform on the HPSC, and guidance on how the SPLICE algorithms can be implemented to best utilize the DLC platform. This insight can also be used to derive requirements to guide trade studies on candidate computing architectures, for use on platforms like Blue Moon. NASA and Blue Origin are collaborating under an agreement to pursue this mutual benefit. Performance metrics collected are based on measurement of common compute resources such as percentage used of memory bandwidth, I/O utilization, interrupt latency, and kernel vs. user space code residency. Where possible existing performance counters and metrics that are part of the operating system kernel are used. As the design has a significant FPGA component, performance counters are identified and instantiated in the fabric to measure DMA performance and interface metrics. Collection of metrics is performed on the DLC with a representative workload that simulates a full landing cycle of the Blue Origin New Shepard vehicle. Consideration is given to the other compute implementations and whether they can run SPLICE algorithms at the same rate and with the same latency as the DLC. One option being considered is the use of a RISC-V soft core instantiated in a radiation resilient FPGA fabric such as the Xilinx KU60. Select algorithms from the SPLICE code will be run for comparison with the DLC. This paper describes how the DLC is instrumented to collect performance measurements of the SPLICE workloads, preliminary results from these measurements, and their implications on SPLICE algorithm implementation. The results of experimentation to derive candidate requirements for architecture trades on a PL&HA computing system are also presented.

computer performance↗

Improved Benchmarking of Cohesive Elements in Abaqus Standard for Predicting Disbond and Delamination in Composite Structures

Traditional approaches for aircraft certification require the assumption of an initial flaw condition, either represented as barely visible impact damage (BVID) or through inclusion of a Teflon insert to serve as surrogate damage. Based on the initial composite damage state, the structure must be shown to demonstrate structural durability and damage tolerance (DaDT) according to the following criteria: a. Damage displays no detrimental growth under cyclic loading. The structure is able to sustain design limit load (DLL). Currently, the only available manner for validating structural performance is through test. Since damage can occur over a wide variety of areas within a structure, this approach has proven to be increasingly expensive and time consuming for composite airframes and acreage structure within the design-test-certification building block. A further complicating factor is the requirement to accurately capture the most critical damage morphologies as a starting condition. To understand the severity of the damage, it is either required to experimentally determine the most critical areas at tremendous expense or rely on legacy data of similar structural testing, which limits design space expansion. A preferred solution is to use advanced analysis to provide improved understanding of load margins for critical locations based on a wide variety of potential starting damage conditions.

Eric M. Roberts↗

Improvements of On-orbit Characterization of Terra MODIS Short-wave Infrared Spectral Bands Out-of-Band Responses

The short-wave infrared (SWIR) bands (5-7, 26) of Terra MODIS, which are co-located with the mid-wav e infrared (MWIR) bands (20-25) on the short and mid-wave infrared (SMIR) Focal Plane Assembly (FPA),have a known issue related to 5.3 μm out-of-band (OOB) thermal leak and electronic crosstalk that was identified prelaunch. As a result, a crosstalk correction algorithm was designed and implemented in the MODIS Level 1B (L1B) calibration. Shortly after the Terra launch, extensive efforts were under taken to help characterize and mitigate the impact due to the OOB response and crosstalk on the SWIR on-orbit calibration and, consequently, the associated L1Bdata products. In addition, special night time day mode (NTDM) operations have been regularly scheduled to derive the crosstalk correction coefficients. Since MODIS does not have a spectral band centered at 5.3 μm, its band 28 (7.325 μm) was chosen as the surrogate sending band to simulate the OOB radiances at 5.3 μm. This was largely based on the measurements from the MODIS Airborne Simulator (MAS) spectrometer field campaigns in the early months after the Terra launch. In the case of Aqua MODIS, the magnitude of the SWIR crosstalk was much smaller and band 25 (4.52 μm) was found to be more effective as the sending band for the crosstalk correction. In recent years, the Terra MODIS photovoltaic (PV) long-wave infrared (LWIR) bands (27-30)electronic crosstalk has increased considerably, especially after the spacecraft safe-mode event occurred in February2016. This accentuated degradation in the PV LWIR performance has also impacted the performance of the SWIR crosstalk correction and thus its calibration and data quality. In this paper, we examine the use of band 25 as the sending band for the Terra MODIS SWIR crosstalk correction and compare its performance with that based on band 28 as the sending band. Results indicate an improvement of on-orbit gain stability for the SWIR calibration and reduced detector to detector and sub-frame to sub-frame striping in the calibrated L1B imagery, especially during the period when the PV LWIR electronic crosstalk has become more severe. This approach has been implemented in the forward production of Terra MODIS Collection 6 and Collection 6.1,starting from July 2019, and is planned to be used for the future reprocessing of MODIS L1B and to help improve the mission-long reflectance calibration and trending stability of the SWIR bands.

MODIS, Terra, Aqua, MODIS, SWIR, calibration, out-↗

Comprehensive Measurement of Microbial Burden in Nutrient-Deprived Cleanrooms

Spacecraft surfaces that are destined to land on potential life-harboring celestial bodies are required to be rigorously cleaned and continuously monitored for spore bioburden as a proxy for spacecraft cleanliness. The NASA standard spore assay (NSA), used for spacecraft bioburden estimates, specifically measures spores that are cultivable, aerobic, resistant to heat shock, and grow at 30˚C in a nutrient-rich medium. Since the vast majority of microorganisms cannot be cultivated using the NSA assay, it is necessary to utilize state-of-the art molecular techniques to better understand the presence of all viable microorganisms, not just those measured with the NSA. In this study, the nutrient-deprived low biomass cleanrooms, where spacecraft are assembled, were used as a surrogate to spacecraft surfaces to measure the ratio of NSA spores in relation to the total viable microorganism population to compare with a 2006 space studies report that estimates that for every 1 spore there is approximately 50,000 viable organisms. Ninety-eight surface wipe samples were collected from the spacecraft assembly facility (SAF) cleanroom at the Jet Propulsion Laboratory (JPL) over a 6-month period. The samples were processed and analyzed using classical microbiology along with molecular assays. Traditional microbiology plating methods were used to determine the cultivable bacterial, fungal, and spore populations. Molecular assays were used to determine the total organisms (TO, dead and live) and the viable organisms (VO, live). The TO was measured using adenine triphosphate (ATP) and quantitative polymerase chain reaction (qPCR) assays. The VO was measured using internal ATP, propidium monoazide (PMA)-qPCR, and flow cytometry (after staining for viable microorganisms) assays. Based on the results, it was possible to establish a ratio between spore counts and VO for each viability assay. The ATP based spore to VO ratio ranged from 149 – 746 and the bacterial PMA-qPCR assay based ratio ranged from 314 – 1491 VO. The most conservative estimate came from FACS, which estimated the ratio to be 12,091 VO per 1 NSA spore. Since archaeal (<1%) and fungal (~2%) populations were negligible, the spore to VO ratios were based on bacterial population estimates. The most conservative ratio from this study can be used as a replacement for the SSB estimate on nutrient-deprived (oligotrophic) desiccated spacecraft surfaces, to estimate the VO from NSA measurements without utilizing state-of-the art molecular methods that are costly and require more biomass than is typically found of spacecraft surfaces.

Venkateswaran, Kasthuri↗

Understanding Resilience Optimization Architectures With an Optimization Problem Repository

Optimizing a system’s resilience can be challenging, especially when it involves considering both the inherent resilience of a robust design and the active resilience of a health management system to a set of computationally-expensive hazard simulations. While prior work has developed specialized architectures to effectively and efficiently solve combined design and resilience optimization problems, the comparison of these architectures has been limited to a single case study. To further study resilience optimization formulations, this work develops a problem repository which includes previously-developed resilience optimization problems and additional problems presented in this work: a notional system resilience model, a pandemic response model, and a cooling tank hazard prevention model. This work then uses models in the repository at large to understand the characteristics of resilience optimization problems and study the applicability of optimization architectures and decomposition strategies. Based on the comparisons in the repository, applying an optimization architecture effectively requires understanding the alignment and coupling relationships between the design and resilience models, as well as the efficiency characteristics of the algorithms. While alignment determines the necessity of a surrogate of resilience cost in the upper-level design problem, coupling determines the overall applicability of a sequential, alternating, or bilevel structure. Additionally, the application of decomposition strategies is dependent on there being limited interactions between variable sets, which often does not hold when a resilience policy is parameterized in terms of actions to take in hazardous model states rather than specific given scenarios.

Resilience↗