Search NASA⌕ Search

SEARCH · Search NASA

Results for “automatic differentiation”

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.

323 records · Page 18

Canonical forms for nonlinear systems

Necessary and sufficient conditions for transforming a nonlinear system to a controllable linear system have been established, and this theory has been applied to the automatic flight control of aircraft. These transformations show that the nonlinearities in a system are often not intrinsic, but are the result of unfortunate choices of coordinates in both state and control variables. Given a nonlinear system (that may not be transformable to a linear system), we construct a canonical form in which much of the nonlinearity is removed from the system. If a system is not transformable to a linear one, then the obstructions to the transformation are obvious in canonical form. If the system can be transformed (it is called a linear equivalent), then the canonical form is a usual one for a controllable linear system. Thus our theory of canonical forms generalizes the earlier transformation (to linear systems) results. Our canonical form is not unique, except up to solutions of certain partial differential equations we discuss. In fact, the important aspect of this paper is the constructive procedure we introduce to reach the canonical form. As is the case in many areas of mathematics, it is often easier to work with the canonical form than in arbitrary coordinate variables.

Su, R.↗

New integration techniques for chemical kinetic rate equations. I - Efficiency comparison

A comparison of the efficiency of several recently developed numerical techniques for solving chemical kinetic rate equations is presented. The solution procedures examined include two general-purpose codes, EPISODE and LSODE, developed as multipurpose differential equation solvers, and three specialzed codes, CHEMEQ, CREK1D, and GCKP84, developed specifically for chemical kinetics. The efficiency comparison is made by applying these codes to two practical combustion kinetics problems. Both problems describe adiabatic, constant-pressure, gas-phase chemical reactions and include all three combustion regimes: induction, heat release, and equilibration. The comparison shows that LSODE is the fastest routine currently available for solving chemical kinetic rate equations. An important finding is that an iterative solution of the algebraic enthalpy conservation equation for temperature can be significantly faster than evaluation of the temperature by integration of its time derivative. Significant increases in computational speed are realized by updating the reaction rate constants only when the temperature change exceeds an amount Delta-T that is problem dependent. An approximate expression for the automatic evaluation of Delta-T is presented and is shown to result in increased computational speed.

Radhakrishnan, K.↗

Aerodynamic Shape Optimization Using Hybridized Differential Evolution

An aerodynamic shape optimization method that uses an evolutionary algorithm known at Differential Evolution (DE) in conjunction with various hybridization strategies is described. DE is a simple and robust evolutionary strategy that has been proven effective in determining the global optimum for several difficult optimization problems. Various hybridization strategies for DE are explored, including the use of neural networks as well as traditional local search methods. A Navier-Stokes solver is used to evaluate the various intermediate designs and provide inputs to the hybrid DE optimizer. The method is implemented on distributed parallel computers so that new designs can be obtained within reasonable turnaround times. Results are presented for the inverse design of a turbine airfoil from a modern jet engine. (The final paper will include at least one other aerodynamic design application). The capability of the method to search large design spaces and obtain the optimal airfoils in an automatic fashion is demonstrated.

Madavan, Nateri K.↗

Comparison of numerical techniques for integration of stiff ordinary differential equations arising in combustion chemistry

The efficiency and accuracy of several algorithms recently developed for the efficient numerical integration of stiff ordinary differential equations are compared. The methods examined include two general-purpose codes, EPISODE and LSODE, and three codes (CHEMEQ, CREK1D, and GCKP84) developed specifically to integrate chemical kinetic rate equations. The codes are applied to two test problems drawn from combustion kinetics. The comparisons show that LSODE is the fastest code currently available for the integration of combustion kinetic rate equations. An important finding is that an interactive solution of the algebraic energy conservation equation to compute the temperature does not result in significant errors. In addition, this method is more efficient than evaluating the temperature by integrating its time derivative. Significant reductions in computational work are realized by updating the rate constants (k = at(supra N) N exp(-E/RT) only when the temperature change exceeds an amount delta T that is problem dependent. An approximate expression for the automatic evaluation of delta T is derived and is shown to result in increased efficiency.

Radhakrishnan, K.↗

NASA's Next Generation Sunphotometer for Ground-based Remote Sensing Applications

Atmospheric aerosol concentrations and their optical properties, in terms of which cause differential warming/cooling effects in the atmosphere and on the surface, constitute one of the largest sources of uncertainty in current assessments and predictions of global climatic change. This is especially true over the regions of bright-reflecting surface, such as desert and urban areas. Under the name of AERONET since the 90's, Cimel's sunphotometers have been deployed worldwide as the standard instrument for aerosol monitoring network, developed to support NASA, CNES, and NASDA's Earth satellite systems. PREDE's skyradiometers, as deployed in SKYNET, serve the similar role. One of the key ingredients for achieving accurate aerosol retrievals from satellite observations is comprehensive understanding of surface spectral BRF's (Bidirectional Reflectance Factors), defined as a ratio of radiance measurements reflected from a targeted surface and from a spectral-angular featureless referencing plate. Although the weather-resistant, automatic, sun/sky-scanning spectroradiometers enable Frequent measurements of atmospheric aerosol optical properties and precipitable water at remote sites, they are too slow for surface BRF measurements (20-25 seconds per 360 degree scan, in addition to filter wheel rotation time). We have designed a next generation sun photometer whose sensor head has no moving parts. A dedicated detector for each channel enables 12 simultaneous measurements ranging from the UV (380 nm) to shortwave-IR (2.13 micron) regions. The scan platform will be capable of traveling 360 degrees in about 6 seconds. This is sufficient to finish a BRDF scan every 30 degrees in azimuth and 15 degrees in elevation in less than 4 minutes. More details about this instrument will be presented, together with its applications to the aerosol and trace gas studies. The current plan for this instrument is to deploy during the EAST-AIRE (East Asian Study of Tropospheric Aerosols: an International Regional Experiment) and BASE-ASIA (Biomass-burning Aerosols in South East-Asia: Smoke Impact Assessment).

Shu, Peter K.↗

Simulation Apps Bring Us Closer to Mars

My work with mathematical modeling and computer simulations began in earnest in 1987 when I signed on with the National Aeronautics and Space Administration, better known as NASA. At the time I was developing computer simulations to support the design of carbon dioxide (CO2) removal systems for the International Space Station (ISS) life support system. CO2 is a human metabolic waste, produced at a rate of a kilogram per day per crewmember. CO2 must be removed from the crew cabin, as it will quickly become detrimental to crew health. The systems used to remove CO2 are typically based on sorption processes, which include complex interactions of heat transfer, mass transfer, and gas flow through porous media. At this time, there were no commercial options for solving a set of coupled mathematical models such as partial differential equations (PDEs); you either force-fit the physics into a thermal analysis package or wrote your own code, complete with discretization, meshing, and solution algorithms. Unfortunately, coding CO2 removal processes from scratch did not allow an appropriate focus on the challenging yet crucial task of understanding and capturing the underlying physics via appropriate mathematical models. Configuration control was often simply ignored due to tight schedules, resulting in a code that would be quicker for a second party to rewrite than modify. In the early 2000s, I decided to move CO2 removal simulations to a platform with built-in meshing, solvers, and postprocessing capabilities, and that could solve user-defined multiphysics PDEs. The program chosen to meet these needs was COMSOL Multiphysics®. Along with freeing the engineer to focus on the underlying physics, a degree of configuration control was automatically achieved via a consistent user interface, thus allowing COMSOL users to share computer models. My team has developed simulations that are already providing a valuable debugging capability for the ISS CO2 removal system, and will provide guidance in upgrades to that system.

Knox, Jim↗

A Proxy Method to Bridge LCA Data Gaps Using Automated Material Classification and Probabilistic Under-Specification

Life cycle assessments (LCAs) are essential for understanding the environmental impacts of material production. However, gaps in life cycle inventory (LCI) data for material and chemical inputs present a key challenge for LCA practitioners, especially in the early design stages. Strategies for filling in these gaps require additional time and expertise, which can hinder the LCA’s completion. This study combined automatic material classification and probabilistic under-specification to create a time-efficient method to fill material LCI data gaps. To illustrate the proposed method, proxy environmental impact distributions were generated using publicly available material LCI data classified into the ChemOnt chemical taxonomy using the open-source chemical classification software ClassyFire. Input materials with data gaps were then classified into the same taxonomy, where proxy environmental impact values could be selected from the available distributions to quickly fill in any data gaps. Although these methods were applied to classify material production processes available in the Federal LCA Commons and Ecoinvent databases, they can be applied to any LCA database. This study shows that classifying materials by their chemical structure produces taxonomies with increased granularity relative to industrial classification, improving the ability of under-specified proxy data to be used for differentiating the environmental impacts of competing designs.

biological databases↗

Modeling, Simulation and Control of a Propellant Mixer

This paper addresses some modeling and control issues for a mixing chamber used in rocket engine testing at NASA Stennis Space Center. The mixer must combine high pressure liquid hydrogen (LH2) and gaseous hydrogen (GH2) to produce and output flow that meets certain thermodynamic properties before it is fed into a test article. More precisely, this paper considers that the quantities to be tracked and/or regulated are mixer internal pressure, exit mass flow, and exit temperature. The available control inputs are given by three value positions, namely those of the GH2, LH2 and exit valves. The mixer is modelled by a system of two nonlinear ordinary differential equations having density and internal energy as states. The model must be simple enough to lend itself to subsequent feedback controller design, yet its accuracy must be tested against real data. For this reason, the model includes function calls to thermodynamic property data. Some structural properties of the resulting model that pertain to controller design, such as controllability and uniqueness of the equilibrium point are shown to hold. Validation of the model against real data is also provided. As a first control approach, a small-signal (linear) model is developed based on the nonlinear model and simulated as well. Pulse disturbances are introduced to the valve positions and the quality of the linear model is ascertained by comparing its behavior against the nonlinear model simulations. Valve control strategies that simulate an operator-in-the-loop scenario are then explored demonstrating the need for automatic feedback control. Classical optimal single-output and multi-output Proportional/Integral controllers are designed based on the linear model and applied to the nonlinear model with excellent results to track simultaneous, constant setpoint changes in desired exit flow, exit temperature, and mixer pressure, as well as to reject unmeasurable but bounded additive step perturbations in the valve positions. A feedback linearization controller is designed and used to achieve tracking and regulation of the outputs over an extended range of the variables of interest.

Richter, Hanz↗

Thermomechanics coupling to Monte Carlo particle transport on unstructured mesh geometries using Cardinal

Geometry deformation due to thermal expansion influences neutron transport in many systems. Studying this phenomenon involves coupling models for neutronics, thermal hydraulics, and solid mechanics. To enable high fidelity modeling of these coupled physics, new capabilities were introduced in Cardinal, coupling OpenMC Monte Carlo particle transport models with MOOSE thermomechanical physics on unstructured moving-mesh geometries. In this work, we present a fully open-source capability leveraging on-the-fly mesh skinning to automatically regenerate OpenMC geometry, which allows multiphysics feedback from temperature, density, and geometry changes. The new capability is verified using an analytic benchmark slab problem, which couples S 2 neutron transport with thermal conduction, convective boundary conditions, Doppler-broadened cross sections, and nonlinear thermal expansion effects along the heated slab. Cardinal reproduces the analytic solutions for the neutron flux, heating, k eff , and temperature with demonstrated convergence in various error terms including mesh resolution and cross section temperature library spacing. For the nominal benchmark conditions and with a fine mesh, maximum relative errors for neutron flux, temperature, and heating are lower than 1%, while errors in integral quantities such as k eff and slab length are within 1 pcm and 48 µm, respectively. This work (i) presents a new numerical approach to thermomechanics coupling with OpenMC models, (ii) is the first (to our knowledge) to utilize a mechanical partial differential equation (PDE) solution to solve the (Griesheimer and Kooreman, 2022) analytic benchmark, and (iii) develops this verified capability within an open-source package.

97 - MATHEMATICS AND COMPUTING↗

Use of piloted simulation for studies of fighter departure/spin susceptibility

The NASA-Langley Research Center has incorporated into its stall/spin research program on military airplanes the use of piloted, fixed-base simulation to complement the existing matrix of unique research testing techniques. The piloted simulations of fighter stall/departure flight dynamics are conducted on the Langley Differential Maneuvering Simulator (DMS). The objectives of the simulation research are reviewed. The rationale underlying the simulation methods and procedures used in the evaluation of airplane characteristics is presented. The evaluation steps used to assess fighter stall/departure characteristics are discussed. Simulation results are presented to illustrate the flight dynamics phenomena dealt with. The considerable experience accumulated in the conduct of piloted stall/departure simulation indicates that simulation provides a realistic evaluation of an airplane's maneuverability at high angles of attack and an assessment of the departure and spin susceptibility of the airplane. This realism is obtained by providing the pilot a complete simulation of the airplane and control system which can be flown using a realistic cockpit and visual display in simulations of demanding air combat maneuvering tasks. The use of the piloted simulation methods and procedures described were found very effective in identifying stability and control problem areas and in developing automatic control concepts to alleviate many of these problems. A good level of correlation between simulated flight dynamics and flight test results were obtained over the many fighter configurations studied in the simulator.

Gilbert, W. P.↗

Vacuum Sealable Container (VSC) and Astronaut Lunar Drill (ALD) for Artemis

Introduction: NASA’s Artemis Program is under development to send first woman and next man to the Moon. Artemis will utilize a suite of new technology for Lunar exploration, including new space vehicles, new space suits, and new Astronaut Tools. Honeybee Robotics has been working with NASA JSC to develop a new Vacuum Sealable Container (VSC) and new Astronaut Lunar Drill (ALD) for the upcoming Artemis missions. Vacuum Sealable Container: Sample return continues to be the “Holy Grail” of space exploration, allowing for the analysis of materials using Earth-based laboratories instead of needing to miniaturize and ruggedize instrumentation for space. The Apollo missions to the Moon had several kinds of Sealable Containers which brought back Lunar samples for analysis [1]. These samples are still being analyzed, fifty years later. The VSC requirements are different from that for Apollo containers and as such, new development was required. One major difference between Artemis samples and those from Apollo is the desire to bring back volatiles which may be part of lunar regolith. The VSC is designed to withstand a high-pressure differential caused by sublimating volatiles. Because of the new, stricter sealing requirements, additional features have been added to the VSC. For example, the seal on the container is required to be more robust, thus required more force to actuate, and the seal must be locked in place with a secondary mechanism. Astronaut Lunar Drill: The ALD is designed to be a multi-functional platform for Lunar sample acquisition. The drill builds on lessons learned from the Apollo Lunar Surface Drill (ALSD), as well as Honeybee’s long history of mechanized sample acquisition devices for space [2]. The main functionality of the ALD is Deep Core Regolith Drilling. Additional functionality includes Surface Rock Coring (SRC), and GeoTech Tools (GTT). The ALD is a rotary-percussive drill designed with deep drilling in mind. The ALD is currently designed to have decoupled rotary and percussion subsystems to allow for maximum battery life and reduced fatigue on the crewmember. Honeybee drill technology will automatically engage the percussion when needed to drill at maximum efficiency. The mechanized drill stand helps improve drilling efficiency; the system utilizes advanced drilling algorithms which only require the crewmember to hold a single switch. Additionally, the stand aids in extraction of deep cores, something which was a problem on Apollo. The SRC functionality of the ALD utilizes Honeybee’s Eccentric Tube Core Breakoff technology to collect and retain rock core samples. This technology has also been infused into the Perseverance rover mission. The ALD is removable from the stand to allow crewmembers to collect samples from large boulders. Bringing back rock cores samples instead of full rocks allows for a wider variety of samples to be returned to Earth for study and puts them in a uniform form-factor for effective sealing and analysis. SRC bits will utilize the power of the drill’s percussion system to drill hard Lunar rocks and expedite sample acquisition. The mechanized stand on the ALD allows for additional attachments for taking geotechnical measurements with a Static Cone Penetrometer (SCP) and a Shear Vane (SV). With the stand, the ALD can take SCP measurements with the touch of a button, storing data for return to Earth. SV measurements utilize the ALD’s Rotary motor to spin the vanes in a controlled manner, getting clean data untampered by human error. References: [1] Bar Cohen and Zacny (2009), Drilling in Extreme Environments - Penetration and Sampling on Earth and Other Planets, Wiley. [2] Bar-Cohen and Zacny, Advances in Terrestrial and Extraterrestrial Drilling, CRC Press. [3] Myrick (2003), Core Break-off Mechanism. US Patent No. 6,550,549 Acknowledgements: This work has been supported by NASA via SBIR Phase 3.

Artemis↗

Constrained Multipoint Aerodynamic Shape Optimization Using an Adjoint Formulation and Parallel Computers

An aerodynamic shape optimization method that treats the design of complex aircraft configurations subject to high fidelity computational fluid dynamics (CFD), geometric constraints and multiple design points is described. The design process will be greatly accelerated through the use of both control theory and distributed memory computer architectures. Control theory is employed to derive the adjoint differential equations whose solution allows for the evaluation of design gradient information at a fraction of the computational cost required by previous design methods. The resulting problem is implemented on parallel distributed memory architectures using a domain decomposition approach, an optimized communication schedule, and the MPI (Message Passing Interface) standard for portability and efficiency. The final result achieves very rapid aerodynamic design based on a higher order CFD method. In order to facilitate the integration of these high fidelity CFD approaches into future multi-disciplinary optimization (NW) applications, new methods must be developed which are capable of simultaneously addressing complex geometries, multiple objective functions, and geometric design constraints. In our earlier studies, we coupled the adjoint based design formulations with unconstrained optimization algorithms and showed that the approach was effective for the aerodynamic design of airfoils, wings, wing-bodies, and complex aircraft configurations. In many of the results presented in these earlier works, geometric constraints were satisfied either by a projection into feasible space or by posing the design space parameterization such that it automatically satisfied constraints. Furthermore, with the exception of reference 9 where the second author initially explored the use of multipoint design in conjunction with adjoint formulations, our earlier works have focused on single point design efforts. Here we demonstrate that the same methodology may be extended to treat complete configuration designs subject to multiple design points and geometric constraints. Examples are presented for both transonic and supersonic configurations ranging from wing alone designs to complex configuration designs involving wing, fuselage, nacelles and pylons.

Reuther, James↗

Parallel architectures for iterative methods on adaptive, block structured grids

A parallel computer architecture well suited to the solution of partial differential equations in complicated geometries is proposed. Algorithms for partial differential equations contain a great deal of parallelism. But this parallelism can be difficult to exploit, particularly on complex problems. One approach to extraction of this parallelism is the use of special purpose architectures tuned to a given problem class. The architecture proposed here is tuned to boundary value problems on complex domains. An adaptive elliptic algorithm which maps effectively onto the proposed architecture is considered in detail. Two levels of parallelism are exploited by the proposed architecture. First, by making use of the freedom one has in grid generation, one can construct grids which are locally regular, permitting a one to one mapping of grids to systolic style processor arrays, at least over small regions. All local parallelism can be extracted by this approach. Second, though there may be a regular global structure to the grids constructed, there will be parallelism at this level. One approach to finding and exploiting this parallelism is to use an architecture having a number of processor clusters connected by a switching network. The use of such a network creates a highly flexible architecture which automatically configures to the problem being solved.

Gannon, D.↗

Discrete control of linear distributed systems with application to the deformable primary mirror of a large orbiting telescope

One of the more significant technological problems associated with the orbital operation of large astronomical telescope's is the fabrication and maintenance of the primary mirror surface to the tolerance required for diffraction-limited performance. An interesting approach to the solution of this problem involves continuously measuring and automatically correcting the optical surface of a thin deformable mirror by means of discrete actuators located on its rear surface: The realization of diffraction-limited performance from a telescope in space by this method rests on the ability of the designer to achieve extremely accurate control of a highly complex, interacting, multivariable system. This paper presents the results of a detailed study of the discrete control of linear distributed systems with specific application to the design of a practical controller for a plant representative of a telescope primary mirror for an orbiting astronomical observatory. The problem of controlling the distributed plant is treated by employing modal techniques to represent variations in the optical figure. Distortion of the mirror surface, which arises primarily from thermal gradients, is countered by actuators working against a backing structure to apply a corrective force distribution to the controlled surface. Each displacement actuator is in series with a spring attached to the mirror by means of a pad intentionally introduced to restrict the excitation of high-order modes. Control is then exerted over a finite number (equal.to the number of actuators) of the most significant modes. Through the application of the modal expansion technique the mirror equation of motion is transformed tb a set of uncoupled, linear, time-invariant, ordinary differential equations. The desired dynamic response and static accuracy may then be achieved by the application of classical single-variable design techniques. The formulation of a quadratic performance index which incorporates a measure of image quality permits determination of the trade-off between the-number of actuators and optical purity. A criterion for defining actuator placement and pad size is presented which minimizes the tendency of the controller to excite the unmonitored modes.

Jeremiah F Creedon↗

Venus tectonic styles and crustal differentiation

Two of the most important constraints are known from Pioneer Venus data: the lack of a system of spreading rises, indicating distributed deformation rather than plate tectonics; and the high gravity/topography ratio, indicating the absence of an asthenosphere. In addition, the high depth/diameter ratios of craters on Venus indicate that Venus probably has no more crust than Earth. The problems of the character of tectonics and crustal formation and recycling are closely coupled. Venus appears to lack a recycling mechanism as effective as subduction, but may also have a low rate of crustal differentiation because of a mantle convection pattern that is more distributed, less concentrated, than Earth's. Distributed convection, coupled with the nonlinear dependence of volcanism on heat flow, would lead to much less magmatism, despite only moderately less heat flow, compared to Earth. The plausible reason for this difference in convective style is the absence of water in the upper mantle of Venus. We have applied finite element modeling to problems of the interaction of mantle convection and crust on Venus. The main emphasis has been on the tectonic evolution of Ishtar Terra, as the consequence of convergent mantle flow. The early stage evolution is primarily mechanical, with crust being piled up on the down-stream side. Then the downflow migrates away from the center. In the later stages, after more than 100 m.y., thermal effects develop due to the insulating influence of the thickened crust. An important feature of this modeling is the entrainment of some crustal material in downflows. An important general theme in both convergent and divergent flows is that of mixing vs. stratification. Models of multicomponent solid-state flow obtain that lower-density crustal material can be entrained and recycled, provided that the ration of low-density to high-density material is small enough (as in subducted slabs on Earth). The same considerations should apply in upflows; a small percent of partial melt may be carried along with its matrix and never escape to the surface. Models that assume melt automatically rising to the crust and no entrainment or other mechanism of recycling lower-density material obtain oscillatory behavior, because it takes a long time for heat to build up enough to overcome a Mg-rich low-density residuum. However, these models develop much thicker crust than consistent with estimates from crater depth/diameter ratios.

Kaula, W. M.↗

Automated Fluidics Device for Extraction and Quantification of miRNA Biomarkers From Blood

Radiation Assessment DuRing Exposure And long-Duration Spaceflight (RADREADS) demonstrates space-compatible point-of-care technology for quantitative biological monitoring of blood miRNA biomarkers in response to long-term low dose radiation exposure. This individualized monitoring approach will inform targeted treatment strategies to maximize medical resource utilization by accounting for individual susceptibility to radiation-related illnesses. As human spaceflight progresses beyond Earth’s magnetic shielding, radiation exposure poses a significant risk to astronaut health and safety. Extended operation in this environment comes with an increased risk of radiation exposure, leading to higher risks of radiation sickness, cancer, central nervous system effects, and degenerative diseases. While conventional physical dosimetry techniques capture radiation dose, individualistic susceptibility to radiation damage is varied. Multiple characteristics, including age, body weight, sex, genetics, and immune status, have been found to influence radiosensitivity (Liu et al. 2011, and Bouffler 2016). This differential response necessitates individualized monitoring and targeted treatment strategies to maximize medical resource utilization; however, a practical diagnostic platform for quantifying long-term, low dose radiation-induced tissue damage does not currently exist. MicroRNAs (miRNAs) are a class of small, non-coding RNAs that regulate gene expression by mediating the degradation of messenger RNA. The levels of particular miRNAs are influenced by biological processes such as inflammation and serve as biomarkers for a variety of conditions including cancer (Singh et al. 2017). MicroRNAs are found in various bodily fluids and are amenable to collection via liquid biopsies, providing a minimally invasive and easily quantifiable readout for a variety of radiosensitive reporters. A preliminary signature of 15 spaceflight sensitive miRNA has been identified in rodent and human studies, including miR-21-5p, miR-24-3p, miR-92a-3p, miR-17-5p, miR-16a-3p, miR-34a-3p, and miR-223-3p. These targets generally increased expression with radiation dose and linear energy transfer, though variation between individuals is not yet described. Current gaps in the field include a lack of understanding of longitudinal biological responses to long-term, low dose radiation exposure and the absence of space-compatible point-of-care technology for quantitative biological monitoring. In this body of work, we aim to develop an automated bleed-to-read system to process whole blood for the detection of miRNA biomarkers in order to monitor individualistic responses to radiation exposure. This will be achieved via separating serum (or plasma) from whole blood, followed by extraction, amplification, and quantification of the miRNA using a RT-qPCR reaction. Previously, the WetLab-2 hardware enabled execution of a RT-qPCR reaction aboard ISS; however, it is a manual system that requires crew manipulation and bulky components (Parra et al. 2017). To address these issues, automated fluid handling hardware was developed for each stage of sample preparation. Extraction of total RNA is achieved by sequentially pumping reagents through an off-the-shelf nucleic acid binding column (miRNeasy Serum/Plasma Advanced Kit, Qiagen). This approach eliminates several manual pipetting and centrifuging steps and limits the use of toxic chemicals commonly found in other sample processing techniques. The resulting elution will then be automatically dispensed for RT-qPCR analysis using a compact rotary qPCR (Mic qPCR Cycler, Bio Molecular Systems) that will improve spaceflight compatibility by removing bubbles from the detection region, another challenge highlighted by WetLab-2 (Parra et al. 2017). Efforts are also being made to simplify the RT-qPCR reaction to a 1-step air-dryable mix to improve long-term reagent stability at room temperature and reduce system complexity. By automating the RT-qPCR processes via microfluidic manipulation, RADREADS will reduce crewmember hands-on time and enable the personalized detection of radiation-induced tissue damage during long duration missions. Minimally invasive, longitudinal monitoring of individual’s response to radiation exposure will inform how the physiological system responds to long-term low dose space radiation and enables development of targeted countermeasures by the medical team. Ultimately, this portable technology will require minimal technical expertise and can also be used to monitor miRNA biomarkers associated with other diseases.

Tristen Head↗

Template Matching Approach to Signal Prediction

A new approach to signal prediction and prognostic assessment of spacecraft health resolves an inherent difficulty in fusing sensor data with simulated data. This technique builds upon previous work that demonstrated the importance of physics-based transient models to accurate prediction of signal dynamics and system performance. While models can greatly improve predictive accuracy, they are difficult to apply in general because of variations in model type, accuracy, or intended purpose. However, virtually any flight project will have at least some modeling capability at its disposal, whether a full-blown simulation, partial physics models, dynamic look-up tables, a brassboard analogue system, or simple hand-driven calculation by a team of experts. Many models can be used to develop a predict, or an estimate of the next day s or next cycle s behavior, which is typically used for planning purposes. The fidelity of a predict varies from one project to another, depending on the complexity of the simulation (i.e. linearized or full differential equations) and the level of detail in anticipated system operation, but typically any predict cannot be adapted to changing conditions or adjusted spacecraft command execution. Applying a predict blindly, without adapting the predict to current conditions, produces mixed results at best, primarily due to mismatches between assumed execution of spacecraft activities and actual times of execution. This results in the predict becoming useless during periods of complicated behavior, exactly when the predict would be most valuable. Each spacecraft operation tends to show up as a transient in the data, and if the transients are misaligned, using the predict can actually harm forecasting performance. To address this problem, the approach here expresses the predict in terms of a baseline function superposed with one or more transient functions. These transients serve as signal templates, which can be relocated in time and space against the signal background. One then has the ability to reconstruct a signal regardless of the precise timing of the transients. During operation, one applies the actual start times of spacecraft activities as they occur, and produces a reconstructed, accurate predict in real-time. This general approach is valid under two important conditions. First, the transients themselves must be time-invariant. Second, the transients must be reasonably consistent with respect to different operating points. Both of these assumptions are generally valid, but in the case of a complicated system with numerous types of overlapping transients, this approach may not be effective. Fortunately, there tend to be few transients in spacecraft telemetry of sensor quantities or low-level health and status information because these signals rarely reflect multiple different types of operation. Furthermore, if the predict is at least reasonably close to actual operation, the shift either in time, or in the operating point when the transient occurs - is likely to be small. The proposed approach considers three ways to recognize a transient. The first and most reliable is to use a different signal that identifies operating mode - often the transient will be correlated to a change in operating mode, and this change can usually be detected positively from discrete signals in spacecraft telemetry. If there is no useful mode signal, the second option is to identify the transient template by hand, and detect the onset of the transient using a curve-fitting approach. Finally, if one elects not to choose by hand, there is an option for automatic selection. This technique has been applied to sensor data from several JPL missions and industrial applications, demonstrating an improvement in accurate prediction of future behavior and early detection of system problems. This technique is applicable to practically any time-varying, quantitative sensor measurement.

Mackey, Ryan↗