Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Formal analysis”

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 1,225 records · Page 68

Cross Cutting Structural Design for Exploration Systems

The challenge of our new National Space Policy and NASA's Vision for Space Exploration (VSE) is keyed to the development of more effective space access and transportation systems. Optimizing in-space systems through innovative cross cutting structural designs that reduce mass, combine functional requirements and improve performance can significantly advance spacecraft designs to meet the ever growing demands of our new National Space Policy. Dependence on limited structural designs is no longer an option. We must create robust materials, forms, function and evolvable systems. We must advance national policy objectives in the design, development, test and operation of multi-billion dollar new generation crew capsules by enabling them to evolve in meeting the requirements of long duration missions to the moon and mars. This paper discusses several current issues and major design drivers for consideration in structural design of advanced spacecraft systems. Approaches to addressing these multifunctional requirements is presented as well as a discussion on utilizing Functional Analysis System Technique (FAST) in developing cross cutting structural designs for future spacecraft. It will be shown how easy it is to deploy such techniques in any conceptual architecture definition or ongoing preliminary design. As experts in merging mission, safety and life support requirements of the frail human existence into robust vehicle and habitat design, we will conquer the final frontier, harness new resources and develop life giving technologies for mankind through more innovative designs. The rocket equation tells us that a reduction in mass optimizes our propulsive results. Primary and secondary structural elements provide for the containment of gases, fluids and solids; translate and sustain loads/impacts; conduct/radiate thermal energy; shield from the harmful effects of radiation; provide for grounding/bonding of electrical power systems; compartmentalize operational functions; and provide physical interface with multiple systems. How can we redefine, combine, substitute, rearrange and otherwise modify our structural systems to reduce mass? New technologies will be needed to fill knowledge gaps and propagate new design methods. Such an integrated process is paramount in maintaining U.S. leadership and in executing our national policy goals. The cross cutting process can take many forms, but all forms will have a positive affect on the demanding design environment through initial radical thinking. The author will illustrate such cross cutting results achievable through a formal process called FAST. The FAST example will be used to show how a multifunctional structural system concept for long duration spacecraft might be generated.

Semmes, Edmund B.↗

Reply to “A Comment on ‘an Evolutionary System of Mineralogy: Proposal for A Classification of Planetary Materials Based on Natural Kind Clustering’”

I welcome the “Comment” from Hatert et al. (2021) related to the proposal for an “Evolutionary system of mineralogy” (Hazen 2019) and thank them for their historically informed, conceptually nuanced, and consistently constructive contribution. They offer corrections related to two facets of my paper that seemed unfairly to criticize aspects of the International Mineralogical Association’s Commission on New Minerals, Nomenclature and Classification (IMA-CNMNC) protocols for classifying minerals. First, they note an unfortunate inferred ambivalence with respect to the relationship between the IMA system and the new evolutionary system. If I was once ambivalent, that view has changed. Having spent the past two years in an ongoing effort to develop this new historical approach, I am struck every day by the power of the IMA-CNMNC system of species classification and nomenclature, which is fundamental and indispensable to the science of mineralogy. As Hatert et al. suggest, any new approach to organizing natural solids, including one focused on planetary evolution, must rest on the foundation provided by the IMA-CNMNC and its many volunteers who selflessly bring order to the mineral kingdom. In the best scenario, the evolutionary system may one day emerge as one of several useful approaches that complement and amplify but in no way replace this core IMA-CNMNC foundation, as clearly stated in the abstract of Hazen (2019). Second, Hatert et al. (2021) offer corrections regarding the IMA-CNMNC approach to classification, in particular a mischaracterization of the formal process to incorporate amorphous phases, poorly crystalline materials, and loosely defined “mineraloids.” I am grateful for the clarifications, as well as the implication that IMA protocols may facilitate the embrace of additional such phases in the future. Finally, I welcome the chance to explore further the emerging concept of “natural kinds” as applied to the mineral kingdom. Here, our thoughts differ. I suggest that minerals, considered in their information-rich, idiosyncratic, paragenetic contexts (in contrast to IMA-CNMNC species), have the potential to represent quintessential examples of “natural kinds.” Furthermore, when viewed in their evolutionary context, minerals offer an intriguing opportunity to expand the concept of “historical natural kinds” beyond its present limited and, at times, controversial use in biology, into the realm of the co-evolving geosphere and biosphere.

Philosophy of mineralogy↗

Satellite Data Assimilation

Satellite altimetry provides the only routine observation of a dynamic variable of the global ocean. For instance, measurements from the TOPEX/POSEIDON altimeter have provided dramatic basin-wide images of the space-time evolution of sea level associated with the 1997-98 El Nino event. The nature of large-scale ocean circulation is being studied that underlies such sea level changes measured by satellite altimetry. To first approximation, properties of large-scale (order 1000 km and larger) sea level variability is latitudinally dependent. An analysis using a general circulation model shows that sea level changes are largely due to wind-driven baroclinic (depth dependent) circulation in the tropics (within 20 degrees of the equator), but are primarily due to the expansion and contraction of near surface water in temperate latitudes (between 20 deg. and 40 deg.) forced by seasonal heating and cooling. In contrast, wind-driven barotropic (depth independent) circulation dominates sea level variability in high latitudes (40 deg. and higher) characterized by periods that are as short as a few days. The presence of such inhomogeneity and the significance of high-frequency, large-scale sea level changes had not been fully recognized prior to this study, and are summarized in Fukumori et al. Numerical models provide theoretical relationships among properties that can be inverted using observations so as to estimate the entire state of the ocean, including properties that are otherwise difficult to measure remotely. The process is data assimilation. An approximate Kalman filter and smoother have been devised to assimilate three years of TOPEX/POSEIDON sea level data into a global ocean general circulation model. The figure below demonstrates the skill of the assimilation, and shows the altimeter assimilated estimate being in closer agreement than the simulation is with independent in situ measurements of subsur ace temperature and velocity, consistent with formal uncertainty estimates. Results from this study have been summarized and submitted for publication. The study demonstrates the feasibility of global ocean data assimilation and illustrates applications in monitoring and understanding of processes controlling the evolution of the ocean. Additional information is contained in the original.

Fukumori, Ichiro↗

A Factorial Data Rate and Dwell Time Experiment in the National Transonic Facility

This report is an introductory tutorial on the application of formal experiment design methods to wind tunnel testing, for the benefit of aeronautical engineers with little formal experiment design training. It also describes the results of a Study to determine whether increases in the sample rate and dwell time of the National Transonic Facility data system Would result in significant changes in force and moment data. Increases in sample rate from 10 samples per second to 50 samples per second were examined, as were changes in dwell time from one second per data point to two seconds. These changes were examined for a representative aircraft model in a range of tunnel operating conditions defined by angles of attack from 0 deg to 3.8 degrees, total pressure from 15.0 psi to 24.1 psi, and Mach numbers from 0.52 to 0.82. No statistically significant effect was associated with the change in sample rate. The change in dwell time from one second to two seconds affected axial force measurements, and to a lesser degree normal force measurements. This dwell effect comprises a "rectification error" caused by incomplete cancellation of the positive and negative elements of certain low frequency dynamic components that are not rejected by the one-Hz low-pass filters of the data system. These low frequency effects may be due to tunnel circuit phenomena and other sources. The magnitude of the dwell effect depends on dynamic pressure, with angle of attack and Mach number influencing the strength of this dependence. An analysis is presented which suggests that the magnitude of the rectification error depends on the ratio of measurement dwell time to the period of the low-frequency dynamics, as well as the amplitude of the dynamics The essential conclusion of this analysis is that extending the dwell time (or, equivalently, replicating short-dwell data points) reduces the rectification error.

DeLoach, R.↗

Nuclear-level effective theory of 𝜇→𝑒 conversion: Formalism and applications

Over the next decade new 𝜇→𝑒 conversion searches at Fermilab (Mu2e) and J-PARC (COMET, DeeMe) are expected to advance limits on charged lepton flavor violation (CLFV) by more than four orders of magnitude. Here, by considering the consequence of 𝑃 and 𝐶⁢𝑃 on elastic 𝜇→𝑒 conversion and the structure of possible charge and current densities, we show that rates are governed by six nuclear responses and a single scale, 𝑞/𝑚 𝑁 , where 𝑞 ≈ 𝑚 𝜇 is the momentum transferred from the leptons to the nucleus. To relate this result to microscopic formulations of CLFV, we construct in nonrelativistic effective theory (NRET) the CLFV nucleon-level interaction, pointing out the relevance of the dimensionless scales 𝑦=($^{𝑞⁢𝑏}_2$) 2 ⁢ > |$\overrightarrow{𝑣}_N$⁢|>|$\overrightarrow{𝑣}_μ$⁢|>|$\overrightarrow{𝑣}_T$| , where 𝑏 is the nuclear size, $\overrightarrow{𝑣}_N$ and $\overrightarrow{𝑣}_μ$ are the nucleon and muon intrinsic velocities, and $\overrightarrow{𝑣}_T$ is the target recoil velocity. We discuss previous work, noting the lack of a systematic treatment of the various small parameters. Because the parameter 𝑦 is not small, a proper calculation of 𝜇→𝑒 conversion requires a full multipole expansion of the nuclear response functions, an apparently daunting task with Coulomb-distorted electron partial waves. We demonstrate that the multipole expansion can be carried out to high precision by introducing a simplifying local momentum 𝑞 eff for the electron. Previous work has been limited to simple charge or spin interactions, thereby treating the nucleus effectively as a point particle. We show that such formulations are not compatible with the general form of the 𝜇→𝑒 conversion rate, failing to generate three of the six allowed nuclear response functions. The inclusion of the nucleon velocity $\overrightarrow{𝑣}_N$ yields an NRET with 16 operators and a rate of the general form. Consequently, in the current discovery era for CLFV, it provides the most sensible starting point for experimental analysis, defining what can and cannot be determined about CLFV from the highly exclusive process of 𝜇→𝑒 conversion. Finally, we expand the NRET operator basis to account for the effects of $\overrightarrow{𝑣}_μ$, associated with the muon's lower component, generating corrections to the CLFV coefficients of the point-nucleus response functions. Using advanced shell-model methods, we compute 𝜇→𝑒 conversion rates for a series of experimental targets, deriving bounds on the coefficients of the CLFV operators. These calculations are the first to include a general basis of CLFV operators, full evaluation of the associated nuclear response functions, and an accurate treatment of electron and muon Coulomb effects. We discuss target selection as an experimental “knob” that can be turned to probe the microscopic origins of CLFV. We describe two types of coherence that enhance certain CLFV operators and selection rules that blind elastic 𝜇→𝑒 conversion to others. We discuss the matching of the NRET onto higher level effective field theories, such as those constructed at the light quark level, noting opportunities to build on existing work in direct detection of dark matter. We discuss the relation of 𝜇→𝑒 conversion to 𝜇→𝑒+𝛾 and 𝜇→3⁢𝑒, showing how MEG II and Mu3e results will complement those of Mu2e and COMET. Finally we describe a accompanying script—in Mathematica and Python versions—that can be used to compute 𝜇→𝑒 conversion rates in various nuclear targets for the full set of NRET operators.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Derivation of A Representative Elementary Volume (REV) for Upscaled Two-Phase Flow in Porous Media

Relative permeability plays an important role in the upscaling of multiphase flow in porous media from the pore scale to the Darcy scale. The entire concept of relative permeability is contingent on the existence of a representative elementary volume (REV). As we move to smaller samples to measure relative permeability, such as with digital core analysis, the concept of a classical REV has become increasingly unlikely when using the conventional approach to defining a representative volume. The “‘conventional”’ understanding of an REV is that a large enough volume must be considered such that spatial variability averages out. In digital rock methods, such as pore-scale simulations based on micro-computed tomography (CT) images, the domain size is typically 2 to 4 mm. This is approximately the length scale of a single-phase flow REV using the classic REV approach. However, the single-phase perspective does not consider the complex dynamics and fluctuations often observed in multiphase flow systems, even at centimeter-scale experiments and/or simulations. A fundamental question is, therefore, whether the domain size commonly used in digital rock simulations can provide a consistent energy budget such that the concept of relative permeability exists. Based on first principles, relative permeability accounts for the rate of energy dissipated in a stationary process. If the dynamics are fluctuating, the energy dissipated can vary but will average out over a long enough timescale. The key to determining the validity of the relative permeability is the timescale of the measurement, not the spatial scale. The conventional REV theory assumes that spatial, temporal, and ensemble averages are equivalent in an ergodic system, but it does not provide a way to test this assumption. We provide a formal way to identify the timescale where the relative permeability accurately captures energy dissipation as a way to validate relative permeability measurements and quantitatively assess their accuracy. This result will be tested for a practical SCAL test, determining how long a flow experiment needs to be run to accurately characterize the rate of energy dissipation by the flow. The outcome will be a best practice guide for the determination of relative permeability from core-scale experiments and/or digital core simulations that ensure the energy budget is fully accounted for in the relative permeability coefficient.

Mcclure, James [Virginia Tech, Blacksburg]↗

An XML Representation for Crew Procedures

NASA ensures safe operation of complex systems through the use of formally-documented procedures, which encode the operational knowledge of the system as derived from system experts. Crew members use procedure documentation on the ground for training purposes and on-board space shuttle and space station to guide their activities. Investigators at JSC are developing a new representation for procedures that is content-based (as opposed to display-based). Instead of specifying how a procedure should look on the printed page, the content-based representation will identify the components of a procedure and (more importantly) how the components are related (e.g., how the activities within a procedure are sequenced; what resources need to be available for each activity). This approach will allow different sets of rules to be created for displaying procedures on a computer screen, on a hand-held personal digital assistant (PDA), verbally, or on a printed page, and will also allow intelligent reasoning processes to automatically interpret and use procedure definitions. During his NASA fellowship, Dr. Simpson examined how various industries represent procedures (also called business processes or workflows), in areas such as manufacturing, accounting, shipping, or customer service. A useful method for designing and evaluating workflow representation languages is by determining their ability to encode various workflow patterns, which depict abstract relationships between the components of a procedure removed from the context of a specific procedure or industry. Investigators have used this type of analysis to evaluate how well-suited existing workflow representation languages are for various industries based on the workflow patterns that commonly arise across industry-specific procedures. Based on this type of analysis, it is already clear that existing workflow representations capture discrete flow of control (i.e., when one activity should start and stop based on when other activities start and stop), but do not capture the flow of data, materials, resources or priorities. Existing workflow representation languages are also limited to representing sequences of discrete activities, and cannot encode procedures involving continuous flow of information or materials between activities.

Simpson, Richard C.↗

Probing New Physics with Tau Leptons using the CMS Detector

The Standard Model (SM) of particle physics is a thoroughly well-tested theory of all electroweak scale phenomena to date—the most recent test being the discovery of the long-predicted Higgs boson. However, there is much evidence to suggest that it is incomplete. On the formal side, the mass of the Higgs particle itself is perplexingly small, while the largest theoretical energy scale of the SM is 17 orders of magnitude larger, a puzzle known as the hierarchy problem. On the experimental side, among other things, the SM lacks a candidate for dark matter, which is now known to play an indispensable role in astrophysical large scale structure formation. Many theories that address the fundamental puzzles of the SM predict the existence of multiple Higgs bosons, including a very light boson that couples preferentially to third-generation fermions. Observation of this new light boson, either directly produced or in the decay of the recently discovered SM Higgs, would provide an unmistakable sign of new dynamics in nature. The research supported by this award used data collected by the Compact Muon Solenoid (CMS) detector at the Large Hadron Collider (LHC) to explore boosted di-tau signatures of a new light boson as a tool in the quest to understand the shortcomings of the SM. With a proton-proton collision energy of 13 trillion electron-volts (TeV) and peak instantaneous proton luminosity of over 10^34 collisions per square centimeter per second (cm –2 s –1 ), the LHC is the only facility in the world that could enable this line of research. The research sought to address the hierarchy problem and other shortcomings of the SM through a mixture of LHC data analysis, maintenance of the CMS detector to ensure high quality data, and contributions to the Phase 2 upgrade detector for the upcoming High Luminosity phase of the LHC (HL-LHC).

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Approximating the linear response of physical chaos

Abstract Parametric derivatives of statistics are highly desired quantities in prediction, design optimization and uncertainty quantification. In the presence of chaos, the rigorous computation of these quantities is certainly possible, but mathematically complicated and computationally expensive. Based on Ruelle’s formalism, this paper shows that the sophisticated linear response algorithm can be dramatically simplified in higher-dimensional systems featuring statistical homogeneity in the physical space. We argue that the contribution of the SRB (Sinai–Ruelle–Bowen) measure gradient, which is an integral yet the most cumbersome part of the full algorithm, is negligible if the objective function is appropriately aligned with unstable manifolds. This abstract condition could potentially be satisfied by a vast family of real-world chaotic systems, regardless of the physical meaning and mathematical form of the objective function and perturbed parameter. We demonstrate several numerical examples that support these conclusions and that present the use and performance of a simplified linear response algorithm. In the numerical experiments, we consider physical models described by differential equations, including Lorenz 96 and Kuramoto–Sivashinsky.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

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↗

Shock Hugoniot calculations using on-the-fly machine learned force fields with ab initio accuracy

We present a framework for computing the shock Hugoniot using on-the-fly machine learned force field (MLFF) molecular dynamics simulations. In particular, we employ an MLFF model based on the kernel method and Bayesian linear regression to compute the free energy, atomic forces, and pressure, in conjunction with a linear regression model between the internal and free energies to compute the internal energy, with all training data generated from Kohn–Sham density functional theory (DFT). We verify the accuracy of the formalism by comparing the Hugoniot for carbon with recent Kohn–Sham DFT results in the literature. In so doing, we demonstrate that Kohn–Sham calculations for the Hugoniot can be accelerated by up to two orders of magnitude, while retaining ab initio accuracy. We apply this framework to calculate the Hugoniots of 14 materials in the FPEOS database, comprising 9 single elements and 5 compounds, between temperatures of 10 kK and 2 MK. We find good agreement with first principles results in the literature while providing tighter error bars. In addition, we confirm that the inter-element interaction in compounds decreases with temperature.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Symmetry-resolved entanglement entropy, spectra & boundary conformal field theory

We perform a comprehensive analysis of the symmetry-resolved (SR) entanglement entropy (EE) for one single interval in the ground state of a 1 + 1D conformal field theory (CFT), that is invariant under an arbitrary finite or compact Lie group, G. We utilize the boundary CFT approach to study the total EE, which enables us to find the universal leading order behavior of the SREE and its first correction, which explicitly depends on the irreducible representation under consideration and breaks the equipartition of entanglement. We present two distinct schemes to carry out these computations. The first relies on the evaluation of the charged moments of the reduced density matrix. This involves studying the action of the defect-line, that generates the symmetry, on the boundary states of the theory. This perspective also paves the way for discussing the infeasibility of studying symmetry resolution when an anomalous symmetry is present. The second scheme draws a parallel between the SREE and the partition function of an orbifold CFT. This approach allows for the direct computation of the SREE without the need to use charged moments. From this standpoint, the infeasibility of defining the symmetry-resolved EE for an anomalous symmetry arises from the obstruction to gauging. Finally, we derive the symmetry-resolved entanglement spectra for a CFT invariant under a finite symmetry group. We revisit a similar problem for CFT with compact Lie group, explicitly deriving an improved formula for U(1) resolved entanglement spectra. Using the Tauberian formalism, we can estimate the aforementioned EE spectra rigorously by proving an optimal lower and upper bound on the same. In the abelian case, we perform numerical checks on the bound and find perfect agreement.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Goal-Function Tree Modeling for Systems Engineering and Fault Management

The draft NASA Fault Management (FM) Handbook (2012) states that Fault Management (FM) is a "part of systems engineering", and that it "demands a system-level perspective" (NASAHDBK- 1002, 7). What, exactly, is the relationship between systems engineering and FM? To NASA, systems engineering (SE) is "the art and science of developing an operable system capable of meeting requirements within often opposed constraints" (NASA/SP-2007-6105, 3). Systems engineering starts with the elucidation and development of requirements, which set the goals that the system is to achieve. To achieve these goals, the systems engineer typically defines functions, and the functions in turn are the basis for design trades to determine the best means to perform the functions. System Health Management (SHM), by contrast, defines "the capabilities of a system that preserve the system's ability to function as intended" (Johnson et al., 2011, 3). Fault Management, in turn, is the operational subset of SHM, which detects current or future failures, and takes operational measures to prevent or respond to these failures. Failure, in turn, is the "unacceptable performance of intended function." (Johnson 2011, 605) Thus the relationship of SE to FM is that SE defines the functions and the design to perform those functions to meet system goals and requirements, while FM detects the inability to perform those functions and takes action. SHM and FM are in essence "the dark side" of SE. For every function to be performed (SE), there is the possibility that it is not successfully performed (SHM); FM defines the means to operationally detect and respond to this lack of success. We can also describe this in terms of goals: for every goal to be achieved, there is the possibility that it is not achieved; FM defines the means to operationally detect and respond to this inability to achieve the goal. This brief description of relationships between SE, SHM, and FM provide hints to a modeling approach to provide formal connectivity between the nominal (SE), and off-nominal (SHM and FM) aspects of functions and designs. This paper describes a formal modeling approach to the initial phases of the development process that integrates the nominal and off-nominal perspectives in a model that unites SE goals and functions of with the failure to achieve goals and functions (SHM/FM). This methodology and corresponding model, known as a Goal-Function Tree (GFT), provides a means to represent, decompose, and elaborate system goals and functions in a rigorous manner that connects directly to design through use of state variables that translate natural language requirements and goals into logical-physical state language. The state variable-based approach also provides the means to directly connect FM to the design, by specifying the range in which state variables must be controlled to achieve goals, and conversely, the failures that exist if system behavior go out-of-range. This in turn allows for the systems engineers and SHM/FM engineers to determine which state variables to monitor, and what action(s) to take should the system fail to achieve that goal. In sum, the GFT representation provides a unified approach to early-phase SE and FM development. This representation and methodology has been successfully developed and implemented using Systems Modeling Language (SysML) on the NASA Space Launch System (SLS) Program. It enabled early design trade studies of failure detection coverage to ensure complete detection coverage of all crew-threatening failures. The representation maps directly both to FM algorithm designs, and to failure scenario definitions needed for design analysis and testing. The GFT representation provided the basis for mapping of abort triggers into scenarios, both needed for initial, and successful quantitative analyses of abort effectiveness (detection and response to crew-threatening events).

Patterson, Jonathan D.↗

The Geology of the Pluto System

NASA's New Horizons mission acquired a large set of images and other data making possible thorough geological analysis of landscapes in the Pluto System. Pluto and Charon exhibit strikingly different surface appearances, despite their similar densities and presumed bulk compositions. Systematic investigation, modelling and mapping revealed that much of Pluto's surface is attributed to surface-atmosphere interactions and the mobilization of volatile ices by insolation. Many mapped valley systems appear to be the consequence of glaciation involving nitrogen ice. Other geological activity requires or required long periods of internal heating, such as Pluto's extensive tectonic fabric. The convection and advection of volatile ices in Sputnik Planitia are thought to be powered by present-day radiogenic heat loss. The prominent mountains at the western margin of Sputnik Planitia, and the strange, multi-km-high mound features to the south, probably composed of H2O, are young geologically as inferred by light cratering and superposition relationships. These multi-km-high mound features might be cryo-volcanoes. Their origin, and what drove their formation so late in Solar System history, is under investigation. East of Sputnik Planitia are large fields of aligned ranks of sharp-crested ridges found only at high altitude that are apparently composed of massive deposits of CH4 ice, referred to as Bladed Terrain. New Horizons found evidence that Bladed Terrain may cover much of Pluto's low latitudes and may have originally formed there as a consequence of Pluto's very high obliquity. Currently Bladed Terrain is undergoing net erosion. This observation, along with evidence for formally more extensive nitrogen glaciation implies that Pluto undergoes significant climate evolution. The dynamic remolding of landscapes by volatile transport seen on Pluto is not unambiguously evident in the mapping of Charon. Charon does, however, display a large resurfaced plain and globally engirdling extensional tectonic network attesting to its early endogenic vigor.

Moore, J. M.↗

Normal stability of slow manifolds in nearly periodic Hamiltonian systems

Kruskal [J. Math. Phys. 3, 806 (1962)] showed that each nearly periodic dynamical system admits a formal U(1) symmetry, generated by the so-called roto-rate. We prove that such systems also admit nearly invariant manifolds of each order, near which rapid oscillations are suppressed. We study the nonlinear normal stability of these slow manifolds for nearly periodic Hamiltonian systems on barely symplectic manifolds—manifolds equipped with closed, non-degenerate 2-forms that may be degenerate to leading order. In particular, we establish a sufficient condition for long-term normal stability based on second derivatives of the well-known adiabatic invariant. We use these results to investigate the problem of embedding guiding center dynamics of a magnetized charged particle as a slow manifold in a nearly periodic system. Here, we prove that one previous embedding and two new embeddings enjoy long-term normal stability and thereby strengthen the theoretical justification for these models.

97 MATHEMATICS AND COMPUTING↗

Assessing Habitability: Lessons from the Phoenix Mission

The Phoenix mission's key objective was to search for a habitable zone. The Phoenix lander carried a robotic arm with digging scoop to collect soil and icy material for analysis with an instrument payload that included volatile mineral and organic analysis(3) and soil ionic chemistry analysis (4). Results from Phoenix along with theoretical modeling and other previous mission results were used to evaluate the habitability of the landing site by considering four factors that characterize the environments ability to support life as we know it: the presence of liquid water, the presence of an energy source to support metabolism, the presence of nutrients containing the fundamental building blocks of life, and the absence of environmental conditions that are toxic to or preclude life. Phoenix observational evidence for the presence of liquid water (past or present) includes clean segregated ice, chemical etching of soil grains, calcite minerals in the soil and variable concentrations of soluble salts5. The maximum surface temperature measured was 260K so unfrozen water can form only in adsorbed films or saline brines but warmer climates occur cyclically on geologically short time scales due to variations in orbital parameters. During high obliquity periods, temperatures allowing metabolism extend nearly a meter into the subsurface. Phoenix discovered ~1%w/w perchlorate salt in the soil, a chemical energy source utilized by a wide range of microbes. Nutrient sources including C, H, N, O, P and S compounds are supplied by known atmospheric sources or global dust. Environmental conditions are within growth tolerance for terrestrial microbes. Summer daytime temperatures are sufficient for metabolic activity, the pH is 7.8 and is well buffered and the projected water activity of a wet soil will allow growth. In summary, martian permafrost in the north polar region is a viable location for modern life. Stoker et al. presented a formalism for comparing the habitability of various regions visited to date on Mars that involved computing a habitability probability, defined as the product of probabilities for the presence of liquid water (P(sub lw)), energy (P(sub e)), nutrients (P(sub ch)), and a benign environment (P(sub b)). Using this formalism, they argued that the Phoenix site was the most habitable of any site visited to date by landed missions and warranted a follow up mission to search for modern evidence of life. This paper will review that conclusion in view of more recent information from the Mars Exploration Rovers and Mars Science Lander missions.

Stoker, Carol R.↗

Uncertainty quantification of optical models in fission fragment deexcitation

Here, we take the first step towards incorporating compound nuclear observables at astrophysically relevant energies into the experimental evidence used to constrain optical models, by propagating the uncertainty in two global optical potentials, one phenomenological and one microscopic, to correlated fission observables using the Monte Carlo Hauser-Feshbach formalism. We compare to a wide range of historic and recent experimental fission measurements, and discuss in detail regions of disagreement. We find that the parametric optical model uncertainty in neutron-fragment correlated observables involving neutron energy is significant. On the other hand, we observe that other experimental features, particularly neutron-fragment correlations near the 132 Sn shell closure and the high energy component of neutron spectra, are unlikely to be explained by the optical potential, and will require further experimental and theoretical effort to explain.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

The OpenSE Cookbook: A Practical, Recipe Based Collection of Patterns, Procedures, and Best Practices for Executable Systems Engineering for the Thirty Meter Telescope

The OpenSE Cookbook is an open-sourced collection of patterns, procedures, and best practices targeted for systems engineers who seek guidance on applying model-based and executable systems engineering (MBSE) using SysML. Its content has emerged from the system level modeling effort on the European Framework Program 6 (FP6) and the Thirty Meter Telescope (TMT). The TMT MBSE approach applied the Executable Systems Engineering Method (ESEM) and the open-source Engineering Environment (OpenMBEE) to specify, analyze, and verify requirements of TMT’s Alignment and Phasing System (APS) and the Narrow Field Infrared Adaptive Optics System (NFIRAOS). In these applications, implicit dependencies are made explicit in a formal model through the use of ESEM, OpenMBEE, and SysML modeling constructs. The value proposition for applying this MBSE approach was to establish precise requirements and fine-grained traceability to system designs, and to verify key requirements beginning early in development. The integration of ESEM and the OpenMBEE tooling infrastructure (providing linked-data and web-operability) is a significant added value for the MBSE approach. The APS is responsible for the overall pre-adaptive optics wavefront quality, using starlight to measure wavefront errors and align the TMT optics. In the formally integrated and executable SysML model, simulations are performed to analyze the impact of changed requirements and verify specified constraints for various operational scenarios. The APS team used several modeling patterns to capture information such as the requirements, the operational scenarios, involved subsystems and their interaction points, the estimated or required time durations, and the mass and power consumption. Adaptive optics systems are designed to sense real-time atmospheric turbulence and correct the telescope’s optical beam to remove its effect. The system model for the adaptive optics operational modes was developed to capture sequence behaviors and operational scenarios to run Monte-Carlo simulations for verifying acquisition time, observing efficiency, and operational behavior requirements. The model is particularly useful for investigating the effect of parallelization, identifying interface issues, and re-ordering sequence acquisition tasks. A former version of the Cookbook (which is now updated to MBSE challenges, goals, and lessons learned) included modeling guidelines and conventions for all system aspects, hierarchy levels, and views, which were developed during for the Active Phasing Experiment (APE), an opto-mechatronical system technology demonstrator for the Extremely Large Telescope (ELT). The Cookbook utilizes the above mentioned system models as real-world case-studies to demonstrate and document the applications of the recipes, providing also instructional examples and addressing the available tooling support. The Cookbook is accompanied by a number of SysML models and aodel libraries which facilitate model authoring and maintenance. The Cookbook covers the different aspects of Systems Engineering such as management of Requirements, Design (behavior and structure), Interfaces, Interdisciplinary Integration, Analysis, Trade Studies, and Technical Resources. This paper presents the background, motivation, architecture, and highlights some key content of the Cookbook. For example, interface management, error budget management, requirements verification, Monte Carlo driven analysis, and timing analysis of operational scenarios. The paper discusses how the capabilities of OpenMBEE contributed significantly to the adoption of executable systems engineering.

Brower, Eric↗