Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “dynamic model learning”

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.

391 records · Page 22

High-Performance Computing Optimization for Aladyn – Adaptive Neural Network Molecular Dynamics Mini-Application

This report provides a description and performance evaluation of the optimization techniques for high performance computing (HPC) implementation of the open source Computational Materials mini-application Aladyn (https://github.com/nasa/aladyn). Aladyn is a basic molecular dynamics code written in FORTRAN 2003, which is designed to demonstrate the use of adaptive neural networks (ANNs) in atomistic simulations. The role of ANNs is to efficiently reproduce the very complex energy landscape resulting from the atomic interactions in materials with the accuracy of the more expensive quantum mechanics-based calculations. The ANN is trained on a large set of atomic structures calculated using the density functional theory (DFT) method. While achieving orders of magnitude faster computational performance than DFT, the ANN-based approach was still very computationally demanding compared to the conventional approach of using empirically fitted energy functions. After its initial development, Aladyn was evaluated and optimized by experts at the NASA Advanced Supercomputing (NAS) division to exploit modern supercomputer architectures. The code has been optimized for execution on multicore central processing units (CPUs), including Intel® Skylake microarchitecture, and on graphic accelerators, such as Nvidia® V100 graphic processing units (GPUs), using Open Multi-Processing (OpenMP) and Open Accelerators (OpenACC) programming interfaces. The optimization achieved a speedup of 4.7 times the baseline version on CPU performance and an additional 2.4 times on CPU+GPU performance. Atomistic computer simulations are a fundamental tool in materials research to model material properties form physics-based first principles. Atomic interaction, governed by Quantum Mechanics (QM) require sophisticated and highly computationally demanding mathematical models to calculate [1]. Classical methods use approximate functional forms, empirically fitted through a set of variable parameters to emulate atomic energies as direct functions of atomic coordinates [2]. While empirical potentials are computationally much simpler, allowing simulations of large-scale systems of up to a trillion (1012) atoms [3], they are substantially less accurate compared to quantum calculations and applicable only to very specific atomic configurations or predefined crystallographic phases. A recently suggested approach is to use heuristic machine learning methods [4], such as those based on Adaptive Neural Networks (ANNs) to predict atomic energies, after being trained on a sufficiently large database of QM-calculated structures [5,6]. This approach reduces significantly the computational complexity, allowing for simulations of orders of magnitude larger systems compared to QM-based methods without compromising accuracy. Still, compared to classical methods using empirical energy functions, ANN methods remain two- to three orders of magnitude more computationally demanding. Hence, the computational cost of simulations, together with the need for extensive training of ANNs, still makes the practical implementation of ANN-based methods quite challenging. The purpose of the Aladyn mini-application software [7], available as open source at https://github.com/nasa/aladyn, is to be a testbed for exploring possible optimization strategies to develop highly scalable parallel algorithms for ANN-based atomistic simulations. Aladyn is aimed at utilizing the architecture of the high-end modern highperformance computing (HPC) hardware based on multicore central processing units (CPUs) equipped with graphic processing unit (GPU) accelerators. Specifically, the goal is to optimize the performance on a single HPC compute node, before implementing scaling to multi-node parallelization using message passing interface (MPI). At the same time, the open source code of Aladyn can serve as a training model for students and professors in academia.

Yamakov, Vesselin I.↗

Multi-Zone Liquid Thrust Chamber Performance Code with Domain Decomposition for Parallel Processing

Computational Fluid Dynamics (CFD) has considerably evolved in the last decade. There are many computer programs that can perform computations on viscous internal or external flows with chemical reactions. CFD has become a commonly used tool in the design and analysis of gas turbines, ramjet combustors, turbo-machinery, inlet ducts, rocket engines, jet interaction, missile, and ramjet nozzles. One of the problems of interest to NASA has always been the performance prediction for rocket and air-breathing engines. Due to the complexity of flow in these engines it is necessary to resolve the flowfield into a fine mesh to capture quantities like turbulence and heat transfer. However, calculation on a high-resolution grid is associated with a prohibitively increasing computational time that can downgrade the value of the CFD for practical engineering calculations. The Liquid Thrust Chamber Performance (LTCP) code was developed for NASA/MSFC (Marshall Space Flight Center) to perform liquid rocket engine performance calculations. This code is a 2D/axisymmetric full Navier-Stokes (NS) solver with fully coupled finite rate chemistry and Eulerian treatment of liquid fuel and/or oxidizer droplets. One of the advantages of this code has been the resemblance of its input file to the JANNAF (Joint Army Navy NASA Air Force Interagency Propulsion Committee) standard TDK code, and its automatic grid generation for JANNAF defined combustion chamber wall geometry. These options minimize the learning effort for TDK users, and make the code a good candidate for performing engineering calculations. Although the LTCP code was developed for liquid rocket engines, it is a general-purpose code and has been used for solving many engineering problems. However, the single zone formulation of the LTCP has limited the code to be applicable to problems with complex geometry. Furthermore, the computational time becomes prohibitively large for high-resolution problems with chemistry, two-equation turbulence model, and two-phase flow. To overcome these limitations, the LTCP code is rewritten to include the multi-zone capability with domain decomposition that makes it suitable for parallel processing, i.e., enabling the code to run every zone or sub-domain on a separate processor. This can reduce the run time by a factor of 6 to 8, depending on the problem.

Homayun K. Navaz↗

Automating Trend Analysis for Spacecraft Constellations

Spacecraft trend analysis is a vital mission operations function performed by satellite controllers and engineers, who perform detailed analyses of engineering telemetry data to diagnose subsystem faults and to detect trends that may potentially lead to degraded subsystem performance or failure in the future. It is this latter function that is of greatest importance, for careful trending can often predict or detect events that may lead to a spacecraft's entry into safe-hold. Early prediction and detection of such events could result in the avoidance of, or rapid return to service from, spacecraft safing, which not only results in reduced recovery costs but also in a higher overall level of service for the satellite system. Contemporary spacecraft trending activities are manually intensive and are primarily performed diagnostically after a fault occurs, rather than proactively to predict its occurrence. They also tend to rely on information systems and software that are oudated when compared to current technologies. When coupled with the fact that flight operations teams often have limited resources, proactive trending opportunities are limited, and detailed trend analysis is often reserved for critical responses to safe holds or other on-orbit events such as maneuvers. While the contemporary trend analysis approach has sufficed for current single-spacecraft operations, it will be unfeasible for NASA's planned and proposed space science constellations. Missions such as the Dynamics, Reconnection and Configuration Observatory (DRACO), for example, are planning to launch as many as 100 'nanospacecraft' to form a homogenous constellation. A simple extrapolation of resources and manpower based on single-spacecraft operations suggests that trending for such a large spacecraft fleet will be unmanageable, unwieldy, and cost-prohibitive. It is therefore imperative that an approach to automating the spacecraft trend analysis function be studied, developed, and applied to missions such as DRACO with the intent that mission operations costs be significantly reduced. The goal of the Constellation Spacecraft Trend Analysis Toolkit (CSTAT) project is to serve as the pathfinder for a fully automated trending system to support spacecraft constellations. The development approach to be taken is evolutionary. In the first year of the project, the intent is to significantly advance the state of the art in current trending systems through improved functionality and increased automation. In the second year, the intent is to add an expert system shell, likely through the adaptation of an existing commercial-off-the-shelf (COTS) or government-off-the-shelf (GOTS) tool to implement some level of the trending intelligence that humans currently provide in manual operations. In the third year, the intent is to infuse the resulting technology into a near-term constellation or formation-flying mission to test it and gain experience in automated trending. The lessons learned from the real missions operations experience will then be used to improve the system, and to ultimately incorporate it into a fully autonomous, closed-loop mission operations system that is truly capable of supporting large constellations. In this paper, the process of automating trend analysis for spacecraft constellations will be addressed. First, the results of a survey on automation in spacecraft mission operations in general, and in trending systems in particular will be presented to provide an overview of the current state of the art. Next, a rule-based model for implementing intelligent spacecraft subsystem trending will be then presented, followed by a survey of existing COTS/GOTS tools that could be adapted for implementing such a model. The baseline design and architecture of the CSTAT system will be presented. Finally, some results obtained from initial software tests and demonstrations will be presented.

Davis, George↗

TRMM and Its Connection to the Global Water Cycle

The importance of quantitative knowledge of tropical rainfall, its associated latent heating and variability is summarized in the context of the global hydrologic cycle. Much of the tropics is covered by oceans. What land exists, is covered largely by rainforests that are only thinly populated. The only way to adequately measure the global tropical rainfall for climate and general circulation models is from space. The TRMM orbit is inclined 35 degrees leading to good sampling in the tropics and a rapid precession to study the diurnal cycle of precipitation. The precipitation instrument complement consists of the first rain radar to be flown in space (PR), a multi-channel passive microwave sensor (TMI) and a five-channel VIS/IR (VIRS) sensor. The precipitation radar operates at a frequency of 13.6 GHz. The swath width is 220 km, with a horizontal resolution of 4 km and the vertical resolution of 250 m. The minimum detectable signal from the precipitation radar has been measured at 17 dBZ. The TMI instrument is designed similar to the SSM/I with two important changes. The 22.235 GHz water vapor absorption channel of the SSM/I was moved to 21.3 GHz in order to avoid saturation in the tropics and 10.7 GHz V&H polarized channels were added to expand the dynamic range of rainfall estimates. The resolution of the TMI varies from 4.6 km at 85 GHz to 36 km at 10.7 GHz. The visible and infrared sensor (VIRS) measures radiation at 0.63, 1.6, 3.75, 10.8 and 12.0 microns. The spatial resolution of all five VIRS channels is 2 km at nadir. In addition to the three primary rainfall instruments, TRMM will also carry a Lightning Imaging Sensor (LIS) and a Clouds and the Earth's Radiant Energy System (CERES) instrument. This presentation will focus primarily on the advances in our understanding of tropical rain systems needed to interpret the TRMM data. Global averages, as well as case studies from TRMM radar (PR), the TRMM Microwave Imager (TMI) and Visible and Infrared Sensor (VIRS) will be presented. Comparisons and contrasts among the different sensors will be drawn. Results will also be compared to previous rainfall climatologies generated from the SSM/I instrument. In particular this paper will focus on the synergy between the TRMM radar and passive microwave radiometer and what we have learned from is synergy.

Kummerow, Chiristian↗

TRMM and its Connection to the Global Water Cycle

The importance of quantitative knowledge of tropical rainfall, its associated latent heating and variability is summarized in the context of the global hydrologic cycle. Much of the tropics is covered by oceans. What land exists, is covered largely by rainforests that are only thinly populated. The only way to adequately measure the global tropical rainfall for climate and general circulation models is from space. The Tropical Rainfall Measuring Mission (TRMM) orbit is inclined 35 degrees leading to good sampling in the tropics and a rapid precession to study the diurnal cycle of precipitation. The precipitation instrument complement consists of the first rain radar to be flown in space (PR), a multi-channel passive microwave sensor (TMI) and a five-channel VIS/IR (VIRS) sensor. The precipitation radar operates at a frequency of 13.6 GHz. The swath width is 220 km, with a horizontal resolution of 4 km and the vertical resolution of 250 m. The minimum detectable signal from the precipitation radar has been measured at - 17 dBZ. The TMI instrument is designed similar to the SSM/I with two important changes. The 22.235 GHz water vapor absorption channel of the SSM/I was moved to 21.3 GHz in order to avoid saturation in the tropics and 10.7 GHz V&H polarized channels were added to expand the dynamic range of rainfall estimates. The resolution of the TMI varies from 4.6 km at 85 GHz to 36 km at 10.7 GHz. The visible and infrared sensor (VIRS) measures radiation at 0.63, 1.6, 3.75, 10.8 and 12.0 microns. The spatial resolution of all five VIRS channels is 2 km at nadir. In addition to the three primary rainfall instruments, TRMM will also carry a Lightning Imaging Sensor (LIS) and a Clouds and the Earth's Radiant Energy System (CERES) instrument. This presentation will focus primarily on the advances in our understanding of tropical rain systems needed to interpret the TRMM data. Global averages, as well as case studies from TRMM radar (PR), the TRMM Microwave Imager (TMI) and Visible and Infrared Sensor (VIRS) will be presented. Comparisons and contrasts among the different sensors will be drawn. Results will also be compared to previous rainfall climatologies generated from the SSM/I instrument. In particular this paper will focus on the synergy between the TRMM radar and passive microwave radiometer and what we have learned from its synergy.

Kummerow, Christian↗

Utilization of the Space Vision System as an Augmented Reality System For Mission Operations

Augmented reality is a technique whereby computer generated images are superimposed on live images for visual enhancement. Augmented reality can also be characterized as dynamic overlays when computer generated images are registered with moving objects in a live image. This technique has been successfully implemented, with low to medium levels of registration precision, in an NRA funded project entitled, "Improving Human Task Performance with Luminance Images and Dynamic Overlays". Future research is already being planned to also utilize a laboratory-based system where more extensive subject testing can be performed. However successful this might be, the problem will still be whether such a technology can be used with flight hardware. To answer this question, the Canadian Space Vision System (SVS) will be tested as an augmented reality system capable of improving human performance where the operation requires indirect viewing. This system has already been certified for flight and is currently flown on each shuttle mission for station assembly. Successful development and utilization of this system in a ground-based experiment will expand its utilization for on-orbit mission operations. Current research and development regarding the use of augmented reality technology is being simulated using ground-based equipment. This is an appropriate approach for development of symbology (graphics and annotation) optimal for human performance and for development of optimal image registration techniques. It is anticipated that this technology will become more pervasive as it matures. Because we know what and where almost everything is on ISS, this reduces the registration problem and improves the computer model of that reality, making augmented reality an attractive tool, provided we know how to use it. This is the basis for current research in this area. However, there is a missing element to this process. It is the link from this research to the current ISS video system and to flight hardware capable of utilizing this technology. This is the basis for this proposed Space Human Factors Engineering project, the determination of the display symbology within the performance limits of the Space Vision System that will objectively improve human performance. This utilization of existing flight hardware will greatly reduce the costs of implementation for flight. Besides being used onboard shuttle and space station and as a ground-based system for mission operational support, it also has great potential for science and medical training and diagnostics, remote learning, team learning, video/media conferencing, and educational outreach.

Maida, James C.↗

Structural Dynamics and Data Analysis

This project consists of two parts, the first will be the post-flight analysis of data from a Delta IV launch vehicle, and the second will be a Finite Element Analysis of a CubeSat. Shock and vibration data was collected on WGS-5 (Wideband Global SATCOM- 5) which was launched on a Delta IV launch vehicle. Using CAM (CAlculation with Matrices) software, the data is to be plotted into Time History, Shock Response Spectrum, and SPL (Sound Pressure Level) curves. In this format the data is to be reviewed and compared to flight instrumentation data from previous flights of the same launch vehicle. This is done to ensure the current mission environments, such as shock, random vibration, and acoustics, are not out of family with existing flight experience. In family means the peaks on the SRS curve for WGS-5 are similar to the peaks from the previous flights and there are no major outliers. The curves from the data will then be compiled into a useful format so that is can be peer reviewed then presented before an engineering review board if required. Also, the reviewed data will be uploaded to the Engineering Review Board Information System (ERBIS) to archive. The second part of this project is conducting Finite Element Analysis of a CubeSat. In 2010, Merritt Island High School partnered with NASA to design, build and launch a CubeSat. The team is now called StangSat in honor of their mascot, the mustang. Over the past few years, the StangSat team has built a satellite and has now been manifested for flight on a SpaceX Falcon 9 launch in 2014. To prepare for the final launch, a test flight was conducted in Mojave, California. StangSat was launched on a Prospector 18D, a high altitude rocket made by Garvey Spacecraft Corporation, along with their sister satellite CP9 built by California Polytechnic University. However, StangSat was damaged during an off nominal landing and this project will give beneficial insights into what loads the CubeSat experienced during the crash. During the year, the MIHS students generated a SolidWorks (CAD software) geometry model of StangSat. This model will be imported into FEMAP (Finite Element Analysis (FEA) Software) and a finite element model wiiJ be created to predict the loads encountered during the crash of this rocket. This analysis will require learning how to import CAD models into the FEM, mesh and add constraints and concentrated masses to represent components inside the CubeSat frame, such as circuit boards, batteries and accelerometers. During the analysis the loads will be varied, in effort to duplicate the damage to the CubeSat. Results will then be peer reviewed and documented.

Luthman, Briana L.↗

Technology Alignment and Portfolio Prioritization (TAPP): Advanced Methods in Strategic Analysis, Technology Forecasting and Long Term Planning for Human Exploration and Operations, Advanced Exploration Systems and Advanced Concepts

The Advanced Concepts Office (ACO) at NASA, Marshall Space Flight Center is expanding its current technology assessment methodologies. ACO is developing a framework called TAPP that uses a variety of methods, such as association mining and rule learning from data mining, structure development using a Technological Innovation System (TIS), and social network modeling to measure structural relationships. The role of ACO is to 1) produce a broad spectrum of ideas and alternatives for a variety of NASA's missions, 2) determine mission architecture feasibility and appropriateness to NASA's strategic plans, and 3) define a project in enough detail to establish an initial baseline capable of meeting mission objectives ACO's role supports the decision­-making process associated with the maturation of concepts for traveling through, living in, and understanding space. ACO performs concept studies and technology assessments to determine the degree of alignment between mission objectives and new technologies. The first step in technology assessment is to identify the current technology maturity in terms of a technology readiness level (TRL). The second step is to determine the difficulty associated with advancing a technology from one state to the next state. NASA has used TRLs since 1970 and ACO formalized them in 1995. The DoD, ESA, Oil & Gas, and DoE have adopted TRLs as a means to assess technology maturity. However, "with the emergence of more complex systems and system of systems, it has been increasingly recognized that TRL assessments have limitations, especially when considering [the] integration of complex systems." When performing the second step in a technology assessment, NASA requires that an Advancement Degree of Difficulty (AD2) method be utilized. NASA has used and developed or used a variety of methods to perform this step: Expert Opinion or Delphi Approach, Value Engineering or Value Stream, Analytical Hierarchy Process (AHP), Technique for the Order of Prioritization by Similarity to Ideal Solution (TOPSIS), and other multi­‐criteria decision-making methods. These methods can be labor-intensive, often contain cognitive or parochial bias, and do not consider the competing prioritization between mission architectures. Strategic Decision-Making (SDM) processes cannot be properly understood unless the context of the technology is understood. This makes assessing technological change particularly challenging due to the relationships "between incumbent technology and the incumbent (innovation) system in relation to the emerging technology and the emerging innovation system." The central idea in technology dynamics is to consider all activities that contribute to the development, diffusion, and use of innovations as system functions. Bergek defines system functions within a TIS to address what is actually happening and has a direct influence on the ultimate performance of the system and technology development. ACO uses similar metrics and is expanding these metrics to account for the structure and context of the technology. At NASA technology and strategy is strongly interrelated. NASA's Strategic Space Technology Investment Plan (SSTIP) prioritizes those technologies essential to the pursuit of NASA's missions and national interests. The SSTIP is strongly coupled with NASA's Technology Roadmaps to provide investment guidance during the next four years, within a twenty-year horizon. This paper discusses the methods ACO is currently developing to better perform technology assessments while taking into consideration Strategic Alignment, Technology Forecasting, and Long Term Planning.

Funaro, Gregory V.↗

NASA's Earth Science Flight Program Meets the Challenges of Today and Tomorrow

NASA's Earth science flight program is a dynamic undertaking that consists of a large fleet of operating satellites, an array of satellite and instrument projects in various stages of development, a robust airborne science program, and a massive data archiving and distribution system. Each element of the flight program is complex and present unique challenges. NASA builds upon its successes and learns from its setbacks to manage this evolving portfolio to meet NASA's Earth science objectives. NASA fleet of 16 operating missions provide a wide range of scientific measurements made from dedicated Earth science satellites and from instruments mounted to the International Space Station. For operational missions, the program must address issues such as an aging satellites operating well beyond their prime mission, constellation flying, and collision avoidance with other spacecraft and orbital debris. Projects in development are divided into two broad categories: systematic missions and pathfinders. The Earth Systematic Missions (ESM) include a broad range of multi-disciplinary Earth-observing research satellite missions aimed at understanding the Earth system and its response to natural and human-induced forces and changes. Understanding these forces will help determine how to predict future changes, and how to mitigate or adapt to these changes. The Earth System Science Pathfinder (ESSP) program provides frequent, regular, competitively selected Earth science research opportunities that accommodate new and emerging scientific priorities and measurement capabilities. This results in a series of relatively low-cost, small-sized investigations and missions. Principal investigators whose scientific objectives support a variety of studies lead these missions, including studies of the atmosphere, oceans, land surface, polar ice regions, or solid Earth. This portfolio of missions and investigations provides opportunity for investment in innovative Earth science that enhances NASA's capability for better understanding the current state of the Earth system. ESM and ESSP projects often involve partnerships with other US agencies and/or international organizations. This adds to the complexity of mission development, but allows for a greater scientific return on NASA's investments. The Earth Science Airborne Science Program provides manned and unmanned aircraft systems that further science and advance the use of satellite data. NASA uses these assets worldwide in campaigns to investigate extreme weather events, observe Earth system processes, obtain data for Earth science modeling activities, and calibrate instruments flying aboard Earth science spacecraft. The Airborne Science Program has six dedicated aircraft and access to many other platforms. The Earth Science Multi-Mission Operations program acquires, preserves, and distributes observational data from operating spacecraft to support Earth Science research focus areas. The Earth Observing System Data and Information System (EOSDIS), which has been in operations since 1994, primarily accomplishes this. EOSDIS acquires, processes, archives, and distributes Earth Science data and information products. The archiving of NASA Earth Science information happens at eight Distributed Active Archive Centers (DAACs) and four disciplinary data centers located across the United States. The DAACs specialize by topic area, and make their data available to researchers around the world. The DAACs currently house over 9 petabytes of data, growing at a rate of 6.4 terabytes per day. NASA's current Earth Science portfolio is responsive to the National Research Council (NRC) 2007 Earth Science Decadal Survey and well as the 2010 NASA Response to President Obama's Climate Plan. As the program evolves into the future it will leverage the lessons learned from the current missions in operations and development, and plan for adjustments to future objectives in response to the anticipated 2017 NRC Decadal Survey.

Flight Program↗

FluxSat: Long-term Earth Science Data Record (ESDR) for Terrestrial Gross Primary Production (GPP) based on satellite data calibrated with eddy covariance data

Gross primary production (GPP), the amount of carbon dioxide (CO 2 ) assimilated by plants through photosynthesis, is one of the most variable and uncertain components of the global carbon cycle. Global GPP has been estimated with a number of process-based models, data-driven, and hybrid approaches. Dynamic global vegetation models (DGVMs), driven by observed environmental changes, are used for global carbon budget assessments and long-term (climate) prediction. Benchmarking these and other models globally with data-driven GPP estimates is critical for understanding the land sink and ensuring accurate forecasts of the carbon cycle. In addition, global data-driven GPP estimates are crucial for studies of interannual variability, including trends that are linked to mechanisms with large uncertainties, such as the indirect CO 2 fertilization effect related to greening. In response to a community need for a GPP data set that well captures spatio-temporal variability, we developed FluxSat, a data-driven approach that optimizes the use of satellite reflectance data from the NASA MODerate-resolution Imaging Spectroradiometer (MODIS) on the Terra and Aqua satellites, calibrated using ground-based eddy covariance (EC) data. We are enhancing (spatially, higher resolution) and extending FluxSat (in time, with additional sensors) to create a high quality long term GPP Earth System Data Record (ESDR) for use in model benchmarking, carbon cycle modeling, and studies of trends and interannual variability. Our team’s objectives are to: 1. Update and document the current MODIS FluxSat GPP (daily, 0.05o and 0.5o resolutions) products with latest available MODIS and EC data sets; 2. Extend FluxSat GPP record forward in time with the Visible Infrared Imaging Radiometer Suite (VIIRS) on operational weather satellites going forward; 3. Extend FluxSat GPP record backward in time using the Advanced Very High Resolution Radiometer (AVHRR) on weather satellites dating back to 1981; 4. Provide higher spatial resolution MODIS and VIIRS GPP (0.0083o). 5. Thoroughly evaluate all FluxSat products with independent data; and 6. Create a homogenized long-term GPP record spanning 40+ years. We will discuss plans for this long-term data set that is supported through the NASA Making Earth System Data Records for Use in Research Environments (MEaSUREs) program.

gross Primary Production↗

Visualization of High Latitude Ion Upflow in Support of the Image Mission

The study of the magnetosphere is a 400 year old science that began with the publication by Gilbert, in 1600, of his hypotheses that the Earth was a giant magnet. Since then we have learned many things about the magnetosphere, particularly in the last 40 years of the space age, but we still have many unanswered questions. In spite of the many thousands of observations of this system we still lack a global understanding of how it works. This is due to its large size and tenuous nature that mean that any measurement made of the fields or particles involved only give one a knowledge of the local conditions at a given time. To gain a global perspective through such observations would require the simultaneous operation of thousands of satellites spread throughout the magnetospheric system in addition to observations made on the ground. Such a program would be impractical at least from financial considerations. What is needed for the advancement of magnetospheric physics is to develop the same capabilities that astrophysicists, solar physicists and meteorologists have been using for years --- the ability to stand back from the object under study and see it in its entirety. The challenge for doing this for the magnetosphere is that the particle densities are very low and the material is, for the most part, not luminous. In the last 25 years several ideas have been proposed that would allow at least the imaging of certain portions of the magnetosphere. These include imaging of the plasmasphere through the resonant scattering of solar 304 A from He+ ions, imaging of various hot plasma populations (i.e. the ring current, plasmasheet, upflowing ionospheric ions, etc.) from the neutral atoms that result when ions of these populations charge exchange with the hydrogen geocorona, and imaging the aurora at various wavelengths in the far ultraviolet. In addition, a novel technique for probing various boundaries in the magnetosphere by bouncing low frequency radio waves off of them has been extensively studied. Such a technique is analogous to the way the under water world can be probed with sonar. About five years ago NASA convened a science working group to study the possibility of flying a magnetospheric imaging mission. This resulted in a number of proposals for such a mission, one of which was selected to be the first MIDEX mission, to be launched in early 2000. The mission is called IMAGE (Imager for Magnetopause to Aurora Global Exploration) and its P.I. is J. Burch at SwRI. The IMAGE spacecraft will carry imagers to view the plasmasphere, aurora, ring current, inner plasmasheet, and upflowing ionospheric ions as well as a radio sounder to probe the location, shape and dynamics of the magnetopause, plasmapause, etc. Between its selection last April and the non advocacy mission review, which takes place next spring, the IMAGE teams needs to further refine the design of the mission and its instruments. The theory and modeling (T&M) subgroup of this team has the task of demonstrating what kind of images the instruments on IMAGE will see as well as showing that useful scientific information can be extracted from such images. As a central element to the efforts of the T&M subgroup we have decided to simulate and create synthetic images for the magnetic cloud event of October, 1995. In this event a large cloud, with high plasma densities and strong magnetic fields, ejected from the sun collided with the earth's magnetosphere triggering a three day period of intense magnetic storms and substorms. This event was observed from a number of different spacecraft and on the ground so we have a good data set to work with. In our work we will place the IMAGE spacecraft in the magnetosphere on its proposed orbit, with its proposed instruments, to see what it would see had it been there. Existing models of the plasmasphere, ring current and magnetopause will be run for this event to give the structures for the imaging instruments. There are several models which are lacking and which need to be developed. These include a model for the cusp, the inner plasmasheet and the upflowing ions. My task this summer was to develop the upflowing ion model and use it to create synthetic images.

Wilson, Gordon R.↗

Spacecraft Modeling, Attitude Determination, and Control: Quaternion-Based Approach

Spacecraft attitude determination and control is an important part of a spacecraft to achieve its designed mission. As of today, many spacecrafts have been successfully launched, and most of them have performed well as they were designed. Many research papers have been published to address the attitude determination and control design problems. Several text books are available for students to learn the technology and for engineers to use as references. The most popular spacecraft models for attitude determination algorithms and control design methods are the Euler angle models and the quaternion models . The Euler angle models have been proved very efficient as the linearized models are controllable, and all standard linear control system design methods are directly applicable. The drawbacks related to the Euler angle methods are (a) the designs based on linearized models may not globally stabilize the original nonlinear spacecraft, i.e., the design may not work when the attitude of the spacecraft is far away from the point where the linearization is performed; (b) the models depend on the rotational sequences, this can be error prone if several teams work on the same project and they use different rotational sequences; (c) for any rotational sequence, there is a singular point where the model is not applicable; and (d) since most attitude determination methods use quaternion to represent the spacecraft attitude, there is a need to transform quaternion into Euler angles. On the other hand, for quaternion models, people have found controllers that can globally stabilize nonlinear spacecraft systems; the models do not depend on rotational sequences and they have no singular point; and the quaternion is provided by attitude determination system and ready to use. The main problem with the quaternion model based control system design is that the linearized quaternion model is not controllable. Therefore, most published design methods heavily rely on Lyapunov functions for the nonlinear spacecraft system. But there is no systematic way to obtain a desired Lyapunov functions. Moreover, the Lyapunov function based designs focus on the closed-loop system stability but pay little attention to the closed-loop system performance. In a series of papers, the author proposed some reduced quaternion models which lead to some controllable linearized spacecraft models. Therefore, all standard linear system theory can be directly applied to analyze and design the spacecraft control systems. We showed that, in some cases, the designed control system is not only optimal for the linearized system, but also globally stabilize the original nonlinear system . Clearly, the reduced quaternion models do not depend on rotational sequences. Due to the special structure of the linearized spacecraft model, some most important design methods, such as LQR design and robust pole assignment design are very simple, enjoy the analytical solutions for some problems, have direct connection to the performance measures, such as settling time, rising time , and percentage of overshoot . All these features are attractive for high quality control system designs. The idea mentioned above is then extended to more spacecraft control problems using specific actuators such as magnetic torque bars and control momentum gyroscopes. These types of actuators may not provide exactly desired torques. Most existing methods use different conversions to get approximate solutions, meaning that these actuators may generate a torque close to but not equal to the desired one. Using the reduced quaternion models that incorporate the actuators into the system model, the control inputs are not torques but the operational parameters. The main benefit of this idea is that the control actions are not approximate but accurate. As all actuators have their operational limit, design with input constraints are also considered in this book by using recently developed interior-point optimization techniques. This book grows up from my research on the spacecraft attitude determination and control design methods in more than a decade which is focused on using reduced quaternion models because of their merits stated above. The book provides all necessary background materials on orbital dynamics, rotations and quaternion, frequently used reference frames, transformations between reference frames, space environment and disturbance torques, ephemeris astronomical vector calculations and measurement instruments, spacecraft control actuators and their models, so that the readers will get a global picture and can apply all these information into the spacecraft system modeling, attitude determination, and spacecraft control system designs, which is the main purpose of this book. This book is different from existing books in that we focus on quaternion based spacecraft control system designs and we consider only attitude control system design related problems, from spacecraft modeling, to attitude determination and estimation, to control system design method selection, to control algorithm development, and to the simulation of the control system designs. Moreover, this book addresses different attitude control tasks in the spacecraft life cycle, including spacecraft maneuver, orbit raising, attitude control, and rendezvous. Finally, this book emphasizes the state space design methods rather than the classical frequency design methods.

Attitude determination and control↗

Shockwave/Boundary-Layer Interaction Studies Performed in the NASA Langley 20-Inch Mach 6 Air Tunnel

This paper highlights results from a collaborative study performed by The University of Tennessee Space Institute (UTSI) and NASA Langley Research Center on the Shockwave/Boundary-Layer Interaction (SWBLI) generated by a cylindrical protuberance on a flat plate in a Mach 6 flow. The study was performed in the 20-Inch Mach 6 Air Tunnel at NASA Langley Research Center and consisted of two separate entries. In the first entry, simultaneous high-speed schlieren and high-speed pressure-sensitive paint (PSP) imaging – which was performed for the first time in the 20-Inch Mach 6 facility at NASA Langley – were performed as well as simultaneous high-speed schlieren and oil-flow imaging. In the second entry, the model configuration was modified to increase the size of the interaction region. High-speed schlieren and infrared thermography (IR) surface imaging were performed in this second entry. The goal of these tests was to characterize the SBLI in the presence of a laminar, transitional, and turbulent boundary layer using high-speed optical imaging techniques. AoA = sting angle-of-attack (°) dcylinder = cylinder diameter (mm) dtrip = cylindrical tripping element diameter (mm) Δshock = shock stand-off distance (mm) hcylinder = cylinder height (mm) htrip = cylindrical tripping element height (mm) HSS = high-speed schlieren M∞ = freestream Mach number PSP = pressure-sensitive paint Re∞ = freestream unit Reynolds number (m-1) SWBLI = shockwave/boundary-layer interaction θplate = model plate angle (°) Introduction his paper highlights two experimental entries performed in the 20-Inch Mach 6 Air Blowdown Tunnel at NASA Langley Research Center in collaboration with The University of Tennessee Space Institute (UTSI). The purpose of these entries was to characterize the dynamic shockwave/boundary-layer interaction (SWBLI) between a vertical cylinder on a flat plate and laminar, transitional (XSWBLI), and turbulent (SWTBLI) boundary layers with a freestream Mach number of 6 using non-intrusive optical diagnostics. Experiments performed by Murphree et al.1,2 were among the first to specifically characterize XSWBLI induced by a vertical cylinder on a flat plate geometry using several optical measurement techniques. Recent optical studies of XSWBLI phenomenon have been performed by UTSI at Mach 2 in their low-enthalpy blow wind tunnel3-8 and by Texas A&M University and UTSI at Mach numbers of 6 and 7 in their Adjustable Contour Expansion wind tunnel.9 The experiments described in this paper were intended to complement previous studies by expanding the freestream unit Reynolds number range, Re∞, over which the XSWBLI phenomena has been observed. Additionally these experiments, made possible under NASA’s new facility funding model under the Aeronautics Evaluation and Test Capabilities (AETC) project, promoted collaboration between university and NASA researchers. The initial entry in the 20-Inch Mach 6 Air Tunnel at NASA Langley occurred in December of 2016. Originally, testing was to occur in November of 2016 in the 31-Inch Mach 10 Air Tunnel at NASA Langley. This facility was chosen so that the XSWBLI phenomenon could be observed at much higher Mach numbers than had previously been attempted in ground test experiments. The model selected for this experiment, a 10° half-angle wedge with a sharp leading edge (described in detail in section II.B), had previously been used by Danehy et al. [10] for boundary layer transition studies using the nitric oxide planar laser-induced fluorescence (NO PLIF) flow visualization technique. In that work, it was determined that transition could be induced downstream of a single htrip = 1-mm tall, dtrip = 4-mm diameter cylindrical tripping element and that the streamwise location of the transition could be changed for a single Re∞ by changing the model angle-of-attack (AoA) (see Fig. A3 in Ref. [10] for more details). Based on the findings of that work, a decision was made to use the wedge model with the cylindrical tripping element to trip the boundary layer flow ahead of a cylindrical protuberance in order to achieve a XSWBLI. Unfortunately, the 31-Inch Mach 10 facility had been taken offline for repairs in October of 2016 and a decision was made to move the test to the 20-Inch Mach 6 facility. Since the behavior of the boundary layer with the chosen model configuration had not been studied before in that facility and the available test time was limited, the entry was considered to be exploratory and was used to collect spatially-resolved and time-resolved flow and surface visualization data that would be used to inform a second entry. Test techniques included simultaneous high-speed schlieren (HSS) captured at 160 kHz and high-speed pressure sensitive paint captured at 10 kHz as well as oil flow visualization, captured at 750 Hz. The second entry in the 20-Inch Mach 6 facility occurred in June and July of 2017. In this follow-on test, modifications to the wind tunnel model were made based on observations made during the first entry and included removing the cylindrical tripping element, increasing the size of the cylinder used to induce the SWBLI to increase the size of the interaction while simultaneously improving spatial resolution, and using a swept ramp array, similar to that described in Ref. [11], to trip the flow to turbulence. Simultaneous HSS (captured at 140 kHz, 100 kHz, and 40 kHz) and conventional IR thermography (captured at 30 Hz) imaging were performed simultaneously in this follow-on entry. This paper is intended to serve as a summary of the work performed during these two entries, to detail lessons learned from each entry, and to highlight some of the datasets acquired. Details on the experimental setup, model configuration, and techniques used are provided. Papers providing a more rigorous analysis of data acquired during the second entry, including statistical, spectral, and modal decomposition methods, can be found in Refs. [12,13]. An entry examining XSWBLI in the 31-Inch Mach 10 Blowdown Wind Tunnel facility is currently planned for mid-to-late calendar year 2019, pending the success of facility repairs. The work performed and described in this paper and the upcoming entry in the 31-Inch Mach 10 facility at NASA Langley have been made possible by NASA’s new facility funding model under the Aeronautics Evaluation and Test Capabilities (AETC) project. Wind Tunnel Facility All experiments discussed in this paper were performed in the 20-Inch Mach 6 Air Tunnel at NASA Langley Research Center. Specific details pertaining to this facility can be found in Refs. [14,15], with only a brief description of the facility provided here. For both entries, the nominal freestream unit Reynolds number was varied between 1.8×106 m-1 (0.5×106 ft-1) and 26.3×106 m-1 (8×106 ft-1). The nominal stagnation pressure was varied between 0.21 MPa and 3.33 MPa and the nominal stagnation temperature was varied between 480 K and 520 K to achieve the desired Re∞ condition. For all runs, the nominal freestream Mach number was 6. The nearly square test section is 520.7-mm (20.5-inches) wide by 508-mm (20-inches) high. Two 431.8-mm (17-inch) diameter windows made of Corning 7940, Grade 5F schlieren-quality glass serve as the side walls of the tunnel and provide optical access for the high-speed schlieren measurements. A rectangular window made of the same material as the side windows served as the top wall of the test section and provided optical access for the high-speed PSP and oil flow measurements. For the second entry, this top window was replaced with a Zinc Selenide (ZnSe) window with an anti-reflection coating capable of passing IR wavelengths between 8μm and 12μm with greater than 98% transmittance. The model was sting supported by a strut attached to a hydraulic system that allows for the model pitch angle to be adjusted between -5° to +55°. For the first entry, an initial pitch/pause sweep of the model AoA was performed to observe the resulting SWBLI. Ultimately, however, the sting pitch angle for this entry was fixed at +10.0° so that the angle of the top surface of the wedge relative to the streamwise axis of the tunnel (referred to herein as the plate angle, θplate), was θplate = 0°. For the second entry, θplate = 0° and θplate = -13.25° were initially tested with the swept ramp array (discussed in the following section) to determine which orientation produced conditions most favorable for XSWBLI to occur based on the heating signatures observed over the top surface of the model in the IR thermography images. Based on these initial tests, θplate = -13.25° was set for the remainder of the runs in the second entry. For both entries, any model changes were performed in a housing located beneath the closed test section. Prior to performing a run of the tunnel, the housing was sealed and the tunnel started. Once the appropriate freestream conditions were achieved, the model was injected into the test section using a hydraulic injection system. B. Model Geometry For all runs, a 10° half-angle (20° full-angle) wedge model with a sharp leading edge was used. The model is described in detail in Refs. [10,16]. The top surface of the sharp leading edge of the model extended 47.8 mm from its upstream-most edge to a junction with the upstream edge of a stainless steel top plate that then extended an (a) (c) (b) Fig. 1 (a) Schematic of top surface of wedge model with gas seeding insert, (b) perspective view of the model in the 20-Inch Mach 6 tunnel with centerline pressure orifices on sharp leading edge, and (c) a perspective view of the model with stainless steel (top) and SLA middle insert (bottom) during the first entry. Flow occurs from left to right.

Bathel, Brett F.↗