Search NASA⌕ Search

SEARCH · Search NASA

Results for “Simulation methods and programs”

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 541 records · Page 30

New Developments in Retropropulsion Testing for Mars Entry, Descent and Landing

NASA’s plans for landing human-scale payloads on Mars in the next decade require that retrorockets be used to decelerate the atmospheric entry vehicle continuously from supersonic conditions through soft touchdown. Conventional Mars entry vehicle architectures that include a single parachute for supersonic-to-subsonic descent are not scalable to the sizes needed to land humans on Mars (~20 metric tons). The major aerosciences risks are the uncertainties in predicting aerodynamic stability and performance during powered free-flight and landing. These risks are influenced partially by current limitations in relevant data and testing methods. Consequently, trajectory simulations currently depend on unvalidated powered descent and landing aerodynamics models. NASA engineers have identified gaps in testing methods that, if addressed, would improve the ability to validate these models. There are gaps in capabilities to test multi-engine hot-gas retropropulsion systems in US wind tunnels. This is partially due to the successful use of parachutes as decelerators for human spaceflight at Earth and for the entire Mars lander program to date. Retropulsion test data historically and to this day are limited to using high pressure air jets, at comparatively low temperatures,as engine plume simulants on subscale wind tunnel models. Additionally, the ability to directly measure aerodynamic interference force and moments is limited by existing flow-through balance capabilities. This paper briefly covers historical and recent test data, and identifies new ground test techniques as a means to provide more relevant test data for powered flight and landing aerodynamic model validation.These techniques include using heated inert gases as a substitute for combustion products, additively manufactured 6-component flow-through force and moment balances, and off-body quantitative diagnostic measurements.

Retropropulsion↗

Implicit Coupling Approach for Simulation of Charring Carbon Ablators

This study demonstrates that coupling of a material thermal response code and a flow solver with nonequilibrium gas/surface interaction for simulation of charring carbon ablators can be performed using an implicit approach. The material thermal response code used in this study is the three-dimensional version of Fully Implicit Ablation and Thermal response program, which predicts charring material thermal response and shape change on hypersonic space vehicles. The flow code solves the reacting Navier-Stokes equations using Data Parallel Line Relaxation method. Coupling between the material response and flow codes is performed by solving the surface mass balance in flow solver and the surface energy balance in material response code. Thus, the material surface recession is predicted in flow code, and the surface temperature and pyrolysis gas injection rate are computed in material response code. It is demonstrated that the time-lagged explicit approach is sufficient for simulations at low surface heating conditions, in which the surface ablation rate is not a strong function of the surface temperature. At elevated surface heating conditions, the implicit approach has to be taken, because the carbon ablation rate becomes a stiff function of the surface temperature, and thus the explicit approach appears to be inappropriate resulting in severe numerical oscillations of predicted surface temperature. Implicit coupling for simulation of arc-jet models is performed, and the predictions are compared with measured data. Implicit coupling for trajectory based simulation of Stardust fore-body heat shield is also conducted. The predicted stagnation point total recession is compared with that predicted using the chemical equilibrium surface assumption

Ablation↗

Simulation and analyses of the aeroassist flight experiment attitude update method

A method which will be used to update the alignment of the Aeroassist Flight Experiment's Inertial Measuring Unit is simulated and analyzed. This method, the Star Line Maneuver, uses measurements from the Space Shuttle Orbiter star trackers along with an extended Kalman filter to estimate a correction to the attitude quaternion maintained by an Inertial Measuring Unit in the Orbiter's payload bay. This quaternion is corrupted by on-orbit bending of the Orbiter payload bay with respect to the Orbiter navigation base, which is incorporated into the payload quaternion when it is initialized via a direct transfer of the Orbiter attitude state. The method of updating this quaternion is examined through verification of baseline cases and Monte Carlo analysis using a simplified simulation, The simulation uses nominal state dynamics and measurement models from the Kalman filter as its real world models, and is programmed on Microvax minicomputer using Matlab, and interactive matrix analysis tool. Results are presented which confirm and augment previous performance studies, thereby enhancing confidence in the Star Line Maneuver design methodology.

Carpenter, J. R.↗

Recent developments of the NESSUS probabilistic structural analysis computer program

The NESSUS probabilistic structural analysis computer program combines state-of-the-art probabilistic algorithms with general purpose structural analysis methods to compute the probabilistic response and the reliability of engineering structures. Uncertainty in loading, material properties, geometry, boundary conditions and initial conditions can be simulated. The structural analysis methods include nonlinear finite element and boundary element methods. Several probabilistic algorithms are available such as the advanced mean value method and the adaptive importance sampling method. The scope of the code has recently been expanded to include probabilistic life and fatigue prediction of structures in terms of component and system reliability and risk analysis of structures considering cost of failure. The code is currently being extended to structural reliability considering progressive crack propagation. Several examples are presented to demonstrate the new capabilities.

Millwater, H.↗

Problems Involved in an Emergency Method of Guiding a Gliding Vehicle from High Altitudes to a High Key Position

An investigation has been conducted to determine the problems involved in an emergency method of guiding a gliding vehicle from high altitudes to a high key position (initial position) above a landing field. A jet airplane in a simulated flameout condition, conventional ground-tracking radar, and a scaled wire for guidance programming on the radar plotting board were used in the tests. Starting test altitudes varied from 30,000 feet to 46,500 feet, and starting positions ranged 8.4 to 67 nautical miles from the high key. Specified altitudes of the high key were 12,000, 10,000 or 4,000 feet. Lift-drag ratios of the aircraft of either 17, 16, or 6 were held constant during any given flight; however, for a few flights the lift-drag ratio was varied from 11 to 6. Indicated airspeeds were held constant at either 160 or 250 knots. Results from these tests indicate that a gliding vehicle having a lift-drag ratio of 16 and an indicated approach speed of 160 knots can be guided to within 800 feet vertically and 2,400 feet laterally of a high key position. When the lift-drag ratio of the vehicle is reduced to 6 and the indicated approach speed is raised to 250 knots, the radar controller was able to guide the vehicle to within 2,400 feet vertically and au feet laterally of the high key. It was also found that radar stations which give only azimuth-distance information could control the glide path of a gliding vehicle as well as stations that receive azimuth-distance-altitude information, provided that altitude information is supplied by the pilot.

Jewel, Joseph W., Jr.↗

Atomistic Method Applied to Computational Modeling of Surface Alloys

The formation of surface alloys is a growing research field that, in terms of the surface structure of multicomponent systems, defines the frontier both for experimental and theoretical techniques. Because of the impact that the formation of surface alloys has on surface properties, researchers need reliable methods to predict new surface alloys and to help interpret unknown structures. The structure of surface alloys and when, and even if, they form are largely unpredictable from the known properties of the participating elements. No unified theory or model to date can infer surface alloy structures from the constituents properties or their bulk alloy characteristics. In spite of these severe limitations, a growing catalogue of such systems has been developed during the last decade, and only recently are global theories being advanced to fully understand the phenomenon. None of the methods used in other areas of surface science can properly model even the already known cases. Aware of these limitations, the Computational Materials Group at the NASA Glenn Research Center at Lewis Field has developed a useful, computationally economical, and physically sound methodology to enable the systematic study of surface alloy formation in metals. This tool has been tested successfully on several known systems for which hard experimental evidence exists and has been used to predict ternary surface alloy formation (results to be published: Garces, J.E.; Bozzolo, G.; and Mosca, H.: Atomistic Modeling of Pd/Cu(100) Surface Alloy Formation. Surf. Sci., 2000 (in press); Mosca, H.; Garces J.E.; and Bozzolo, G.: Surface Ternary Alloys of (Cu,Au)/Ni(110). (Accepted for publication in Surf. Sci., 2000.); and Garces, J.E.; Bozzolo, G.; Mosca, H.; and Abel, P.: A New Approach for Atomistic Modeling of Pd/Cu(110) Surface Alloy Formation. (Submitted to Appl. Surf. Sci.)). Ternary alloy formation is a field yet to be fully explored experimentally. The computational tool, which is based on the BFS (Bozzolo, Ferrante, and Smith) method for the calculation of the energetics, consists of a small number of simple PCbased computer codes that deal with the different aspects of surface alloy formation. Two analysis modes are available within this package. The first mode provides an atom-by-atom description of real and virtual stages 1. during the process of surface alloying, based on the construction of catalogues of configurations where each configuration describes one possible atomic distribution. BFS analysis of this catalogue provides information on accessible states, possible ordering patterns, and details of island formation or film growth. More importantly, it provides insight into the evolution of the system. Software developed by the Computational Materials Group allows for the study of an arbitrary number of elements forming surface alloys, including an arbitrary number of surface atomic layers. The second mode involves large-scale temperature-dependent computer 2. simulations that use the BFS method for the energetics and provide information on the dynamic processes during surface alloying. These simulations require the implementation of Monte-Carlo-based codes with high efficiency within current workstation environments. This methodology capitalizes on the advantages of the BFS method: there are no restrictions on the number or type of elements or on the type of crystallographic structure considered. This removes any restrictions in the definition of the configuration catalogues used in the analytical calculations, thus allowing for the study of arbitrary ordering patterns, ultimately leading to the actual surface alloy structure. Moreover, the Monte Carlo numerical technique used for the large-scale simulations allows for a detailed visualization of the simulated process, the main advantage of this type of analysis being the ability to understand the underlying features that drive these processes. Because of the simplicity of the BFS method for e energetics used in these calculations, a detailed atom-by-atom analysis can be performed at any point in the simulation, providing necessary insight on the details of the process. The main objective of this research program is to develop a tool to guide experimenters in understanding and interpreting often unexpected results in alloy formation experiments. By reducing the computational effort without losing physical accuracy, we expect that powerful simulation tools will be developed in the immediate future, which will allow material scientists to easily visualize and analyze processes at a level not achievable experimentally.

Bozzolo, Guillermo H.↗

Computer simulation of shock wave boundary layer interaction

Development of a boundary layer correction for an analysis program previously used by Bauer et al. (1972) to study two-dimensional flow past transonic airfoils at off-design conditions where shocks which interact with the boundary layer appear. Using a method which combines conformal mapping with a finite-difference scheme due to Murman and Cole (1971), an analysis is made of the interaction between shock waves on a supercritical wing section and a turbulent boundary layer, assuming that the shocks are weak and that separation is insignificant. It is shown that relatively simple shock wave/boundary layer interactions of the kind envisioned in this study can be treated by determining the displacement thickness iteratively together with the flow and altering the profile accordingly. It is concluded that the proposed boundary layer correction appears to furnish an essential improvement in cases where both pronounced aft loading and a large enclosed supersonic region of flow are present.

Bauer, F.↗

Data Selection Improvement For MicroBooNE

Data selection is an extremely important part of data analysis for any experiment. Finding a physics result is often the result of sifting through a massive amount of data, keeping data that we believe to be signal, and throwing out data we do not. This process is called data selection. Creating a selection algorithm is an intensive process that must balance keeping enough data to have statistics and maximizing the signal purity of that data. We also need to choose the right reconstruction method, a tool to take raw data from the detector and convert it into physics results. In this study, we used three different reconstruction tools, Pandora, WireCell, and LANTERN, for the MicroBooNE experiment in conjunction to improve the selection algorithm for analysis. For the case of this study, we look into the charged current N proton 0 pions (CCNp0$\pi$) interaction channel. This is the dominant channel for the Short Baseline Neutrino (SBN) program and is expected to be a large contributor to the Deep Underground Neutrino Experiment (DUNE). We first investigated each of the three tools to find out more about their strengths and weaknesses as reconstructions, and compared them to the truth information directly from the MicroBooNE simulation pipeline. We then put together a direct comparison of the three methods to find which method or combination of methods would return the best result for us. While the study is ongoing, we have learned a lot about data selection for the experiment and the differences between the reconstruction tools.

Dillon, Brayden [Fermilab]↗

HRP's Healthcare Spin-Offs Through Computational Modeling and Simulation Practice Methodologies

Spaceflight missions expose astronauts to novel operational and environmental conditions that pose health risks that are currently not well understood, and perhaps unanticipated. Furthermore, given the limited number of humans that have flown in long duration missions and beyond low Earth-orbit, the amount of research and clinical data necessary to predict and mitigate these health and performance risks are limited. Consequently, NASA's Human Research Program (HRP) conducts research and develops advanced methods and tools to predict, assess, and mitigate potential hazards to the health of astronauts. In this light, NASA has explored the possibility of leveraging computational modeling since the 1970s as a means to elucidate the physiologic risks of spaceflight and develop countermeasures. Since that time, substantial progress has been realized in this arena through a number of HRP funded activates such as the Digital Astronaut Project (DAP) and the Integrated Medical Model (IMM). Much of this success can be attributed to HRP's endeavor to establish rigorous verification, validation, and credibility (VV&C) processes that ensure computational models and simulations (M&S) are sufficiently credible to address issues within their intended scope. This presentation summarizes HRP's activities in credibility of modeling and simulation, in particular through its outreach to the community of modeling and simulation practitioners. METHODS: The HRP requires all M&S that can have moderate to high impact on crew health or mission success must be vetted in accordance to NASA Standard for Models and Simulations, NASA-STD-7009 (7009) [5]. As this standard mostly focuses on engineering systems, the IMM and DAP have invested substantial efforts to adapt the processes established in this standard for their application to biological M&S, which is more prevalent in human health and performance (HHP) and space biomedical research and operations [6,7]. These methods have also generated substantial interest by the broader medical community though institutions like the National Institutes of Health (NIH) and the Food and Drug Administration (FDA) to develop similar standards and guidelines applicable to the larger medical operations and research community. DISCUSSION: Similar to NASA, many leading government agencies, health institutions and medical product developers around the world are recognizing the potential of computational M&S to support clinical research and decision making. In this light, substantial investments are being made in computational medicine and notable discoveries are being realized [8]. However, there is a lack of broadly applicable practice guidance for the development and implementation of M&S in clinical care and research in a manner that instills confidence among medical practitioners and biological researchers [9,10]. In this presentation, we will give an overview on how HRP is working with the NIH's Interagency Modeling and Analysis Group (IMAG), the FDA and the American Society of Mechanical Engineers (ASME) to leverage NASA's biomedical VV&C processes to establish a new regulatory standard for Verification and Validation in Computational Modeling of Medical Devices, and Guidelines for Credible Practice of Computational Modeling and Simulation in Healthcare.

Mulugeta, Lealem↗

Probabilistic structural analysis algorithm development for computational efficiency

The PSAM (Probabilistic Structural Analysis Methods) program is developing a probabilistic structural risk assessment capability for the SSME components. An advanced probabilistic structural analysis software system, NESSUS (Numerical Evaluation of Stochastic Structures Under Stress), is being developed as part of the PSAM effort to accurately simulate stochastic structures operating under severe random loading conditions. One of the challenges in developing the NESSUS system is the development of the probabilistic algorithms that provide both efficiency and accuracy. The main probability algorithms developed and implemented in the NESSUS system are efficient, but approximate in nature. In the last six years, the algorithms have improved very significantly.

Wu, Y.-T.↗

Particle/Continuum Hybrid Simulation in a Parallel Computing Environment

The objective of this study was to modify an existing parallel particle code based on the direct simulation Monte Carlo (DSMC) method to include a Navier-Stokes (NS) calculation so that a hybrid solution could be developed. In carrying out this work, it was determined that the following five issues had to be addressed before extensive program development of a three dimensional capability was pursued: (1) find a set of one-sided kinetic fluxes that are fully compatible with the DSMC method, (2) develop a finite volume scheme to make use of these one-sided kinetic fluxes, (3) make use of the one-sided kinetic fluxes together with DSMC type boundary conditions at a material surface so that velocity slip and temperature slip arise naturally for near-continuum conditions, (4) find a suitable sampling scheme so that the values of the one-sided fluxes predicted by the NS solution at an interface between the two domains can be converted into the correct distribution of particles to be introduced into the DSMC domain, (5) carry out a suitable number of tests to confirm that the developed concepts are valid, individually and in concert for a hybrid scheme.

Baganoff, Donald↗

Sensor Lead Wires Positioned on SiC-based Monolithic Ceramic and Fiber- reinforced Ceramic Matrix Composite Subcomponents with Flat and Curved Surfaces

There is strong interest in the development of silicon carbide-based monolithic ceramic and composite materials and components for demanding, high-temperature applications. Thorough characterization of material properties, including high-temperature testing under simulated or actual operating conditions, is a high priority for programs involved in developing these silicon carbide- (SiC) based materials and components. Members of the Sensors and Electronics Technology Branch at the NASA Lewis Research Center are developing minimally intrusive methods of measuring the properties (such as the surface temperature, strain, and heat flux characteristics) of components and subelements that are being tested or operated in hostile, high-temperature environments. Their primary goal is to instrument the test article or operating component with durable sensors that have a minimal effect on test conditions such as the gas flow across the surface of the item and the material response (including the through-thickness conduction of heat). Therefore, the main thrust of their work has been the development of thin-film sensors (e.g., thermocouples or strain gauges) for use on various advanced material test articles, including SiC/SiC composite components. There was a need for a better method of securing sensor lead wires on SiC-based components and subelements that would be tested at temperatures to 1000 C (or higher), to enhance the durability of the overall minimally intrusive sensor system. To address this need, Lewis researchers devised an alternative approach for positioning the sensor lead wires (which are connected to the thin-film sensors) on SiC or SiC/SiC components. A reaction-forming method of joining was used to strongly bond hoop-shaped monolithic SiC and SiC/SiC composite attachments of various sizes to both flat and curved surfaces of SiC/SiC composite subelements (see the photos). This approach is based on an affordable, robust ceramic joining technology, named ARCJoinT, which was developed at Lewis for the joining of SiC-based ceramics and fiber-reinforced composites.

Kiser, James D.↗

Human Performance in Simulated Reduced Gravity Environments

NASA is currently designing a new space suit capable of working in deep space and on Mars. Designing a suit is very difficult and often requires trade‐offs between performance, cost, mass, and system complexity. Our current understanding of human performance in reduced gravity in a planetary environment (the moon or Mars) is limited to lunar observations, studies from the Apollo program, and recent suit tests conducted at JSC using reduced gravity simulators. This study will look at our most recent reduced gravity simulations performed on the new Active Response Gravity Offload System (ARGOS) compared to the C‐9 reduced gravity plane. Methods: Subjects ambulated in reduced gravity analogs to obtain a baseline for human performance. Subjects were tested in lunar gravity (1.6 m/sq s) and Earth gravity (9.8 m/sq s) in shirt‐sleeves. Subjects ambulated over ground at prescribed speeds on the ARGOS, but ambulated at a self‐selected speed on the C‐9 due to time limitations. Subjects on the ARGOS were given over 3 minutes to acclimate to the different conditions before data was collected. Nine healthy subjects were tested in the ARGOS (6 males, 3 females, 79.5 +/- 15.7 kg), while six subjects were tested on the C‐9 (6 males, 78.8 +/- 11.2 kg). Data was collected with an optical motion capture system (Vicon, Oxford, UK) and was analyzed using customized analysis scripts in BodyBuilder (Vicon, Oxford, UK) and MATLAB (MathWorks, Natick, MA, USA). Results: In all offloaded conditions, variation between subjects increased compared to 1‐g. Kinematics in the ARGOS at lunar gravity resembled earth gravity ambulation more closely than the C‐9 ambulation. Toe‐off occurred 10% earlier in both reduced gravity environments compared to earth gravity, shortening the stance phase. Likewise, ankle, knee, and hip angles remained consistently flexed and had reduced peaks compared to earth gravity. Ground reaction forces in lunar gravity (normalized to Earth body weight) were 0.4 +/- 0.2 on the ARGOS, but only 0.2 +/- 0.1 on the C‐9. Discussion: Gait analysis showed differences in joint kinematics and temporal‐spatial parameters between the reduced gravity simulators and with respect to earth gravity. Although most of the subjects chose a somewhat unique ambulation style as a result of learning to ambulate in a new environment, all but two were consistent with keeping an Earth‐like gait. Learning how reduced gravity affects ambulation will help NASA to determine optimal suit designs, influence mission planning, help train crew, and may shed light on the underlying methods the body uses to optimize gait for energetic efficiency. Conclusion: Kinematic and kinetic analysis demonstrated noteworthy differences between an offloaded environment and 1‐g, as would be expected. The analysis showed a trend to change the ambulation style in an offloaded environment to a rolling‐loping walk (resembling crosscountry skiing) with increased swing time. This ambulation modification, particularly in the ARGOS, indicated that the relative kinetic energy of the subject was increased, on average, per the static body weight compared to the 1‐g condition. How much of this was influenced by the active offloading of the ARGOS system is unknown.

Cowley, Matthew↗

Single and Multi-Node Modeling of Direct, Submerged, and Self-Pressurization of A Cryogenic Propellant Tank Using Nodal Tools

The pressurization of cryogenic propellant tanks will always be an important process so long as cryogenic liquids are being considered as fuel sources or used for other in-space applications. Pressure control of the tank ullage is necessary for managing propellant flowrates to an engine or a receiver tank, and modeling of the process is used to predict the pressurant requirements and the amount of propellant boiloff. Direct ullage pressurization is the more traditional approach to tank pressurization, as the physics are straight-forward, and ample test (flight) data have been collected and analyzed over the past several decades. Submerged injection pressurization is an alternate method for tank pressurization and has been shown to reduce pressurant requirements, subcool the propellant, and reduce the risk of ullage collapse. Additionally, the pressurant gas entering the ullage is usually much colder when using the submerged pressurization approach, resulting in reduced propellant boiloff. These benefits are at the expense of vaporizing a small percentage of the propellent. Both tank pressurization methods are viable options for current and future space missions, and it is important to have the capability of analyzing the tank ullage conditions for both approaches. Our previous work has demonstrated the development of a Generalized Fluid System Simulation Program (GFSSP) model, which contains a thermodynamic equilibrium heat and mass transfer subroutine capable of effectively analyzing both direct and submerged pressurization systems [1-2]. This subroutine has most recently been enhanced to include the non-equilibrium effect of pressurant dissolution into the propellant. To date the ullage has always been represented as a single node, and although the simulated single-node temperatures have good comparison with the volume-averaged ullage temperatures computed from test data, the physics of the thermal stratification in the ullage were never captured, and adjustment factors in the model were required. The purpose of this paper is to introduce the development of a multi-node ullage model using GFSSP and to discuss the improvements of the simulated ullage temperature distribution and its resulting effects on ullage heat transfer processes. Test data from the Cryogenic Propellant Storage and Transfer Engineering Developmental Unit (CPST EDU) was used for model validation. For additional comparison, a Thermal Desktop (TD) model was also developed to analyze the CPST EDU direct ullage pressurization tests using both a single node and multi-node approach. The model includes the direct pressurant line, vent line, fill/drain line, and a TD FloCAD Compartment. The TD FloCAD Compartment is employed to represent the liquid and ullage as single volumes inside the tank, to include a liquid/vapor interface, and to generate network level objects such as lumps (analogous to nodes in GFSSP), paths, and ties between the fluid and thermal elements. An established heat load on the model tank was leveraged from a pre-existing higher-fidelity model correlated to CPST EDU test data.

pressurization↗

Comparing Control Performance Between Simulation and Experiment using the Microreactor Automated Control System Testbed

In the advanced reactor domain, a flexible and scalable software/hardware infrastructure is crucial for integrating and validating various control technologies. This study used the Microreactor Automated Control System (MACS) hardware platform as a testbed. MACS was originally designed to mirror Idaho National Laboratory (INL)'s Microreactor Applications Research Validation and Evaluation (MARVEL), a 85-kW thermal fission microreactor. It features control drums for simulated reactivity control; lights that function as a surrogate reactor core, with the brightness being proportional to the reactor power; and light sensors that emulate neutron detectors. To transform MACS into a physical twin of MARVEL for evaluating control methods, the Control and Optimization Modular Modeling Application for Nuclear Deployment (COMMAND) software was employed. This software integrated the hardware with two models of the MARVEL core, based on Reactor Excursion and Leak Analysis Program (RELAP5-3D) and Monte Carlo N-Particle (MCNP) models. The study aimed to demonstrate the gap between control theory and actual practice—a gap that often necessitates empirical adjustments such as control gain retuning, filters, time discretization, and integrator anti-windup measures. Controllers were developed based on increasingly complex simulations without hardware, starting from the base MARVEL model and then introducing actuator saturation constraints and sensor noise. The final control strategy was then tested using MACS, and a comparative performance analysis was conducted.

46 - INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AN↗

Measurement of Insulation Compaction in the Cryogenic Fuel Tanks at Kennedy Space Center by Fast/Thermal Neutron Techniques

The liquid hydrogen and oxygen cryogenic storage tanks at John F. Kennedy Space Center (KSC) use expanded perlite as thermal insulation. Th ere is evidence that some of the perlite has compacted over time, com promising the thermal performance and possibly also structural integr ity of the tanks. Therefore an Non-destructive Testing (NDT) method for measuring the perlite density or void fraction is urgently needed. Methods based on neutrons are good candidates because they can readil y penetrate through the 1.75 cm outer steel shell and through the ent ire 120 cm thickness of the perlite zone. Neutrons interact with the nuclei of materials to produce characteristic gamma rays which are the n detected. The gamma ray signal strength is proportional to the atom ic number density. Consequently, if the perlite is compacted then the count rates in the individual peaks in the gamma ray spectrum will i ncrease. Perlite is a feldspathic volcanic rock made up of the major elements Si, AI, Na, K and 0 along with some water. With commercially available portable neutron generators it is possible to produce simul taneously fluxes of neutrons in two energy ranges: fast (14 MeV) and thermal (25 meV). Fast neutrons produce gamma rays by inelastic scatt ering which is sensitive to Fe and O. Thermal neutrons produce gamma rays by radiative capture in prompt gamma neutron activation (PGNA) and this is sensitive to Si, AI, Na, Kand H. Thus the two energy ranges produce complementary information. The R&D program has three phases: numerical simulations of neutron and gamma ray transport with MCNP s oftware, evaluation of the system in the laboratory on test articles and finally mapping of the perlite density in the cryogenic tanks at KSC. The preliminary MCNP calculations have shown that the fast/therma l neutron NDT method is capable of distinguishing between expanded an d compacted perlite with excellent statistics.

Livingston, R. A.↗

Scenario Complexity for Unmanned Aircraft System Traffic

This work introduces an approach to estimate the complexity of a low-altitude air traffic scenario involving multiple UASs using mathematical programming. Given a set of multi-point UAS flight trajectories, vehicle dynamics, and a conflict resolution algorithm, an abstract model is developed such that it can be solved quickly using a mathematical programming optimization software without running high-fidelity simulations that can be computationally expensive and may not suit real-time apA quick and accurate assessment of complexity for a given traffic scenario can help plan and schedule flights to alleviate traffic bottleneck and mitigate operation risks, especially for unmanned aerial system traffic management where high traffic density or complexity is expected. This work introduces a traffic scenario complexity metric that was constructed based on the number of potential conflicts weighted by the conflict resolution cost associated. The cost associated with a conflict is calculated based on the corresponding conflict resolution maneuvers. To obtain the conflict resolution maneuvers, a MILP-based optimization was formulated with the vehicle model and conflict management parameters incorporated. To evaluate the complexity metrics, an approach of using measurements from high-fidelity simulations was proposed. The scenario complexity measurements for 920 random-generated scenarios were obtained through high-fidelity simulations and treated as the ground truth. Two statistics methods: Pearson and Alternative Conditional Expectations were applied for analysis. The results showed that the number of flights has low correlation with the scenario complexity according to the correlation coefficients calculated by both methods. The Alternative Conditional Expectations method shows that the proposed scenario complexity metric has better correlation with the ground truth than the number of potential conflicts.plications. In the abstract model, each vehicle is represented by a time-varied vector associated with position, speed, and heading information. The total extra distance that aircraft need to divert from their original routes to avoid collisions is computed and used to setup a quadratic programming formula. The metrics including the number of conflicts and extra distances travelled by all vehicles are then utilized to estimate the complexity of a given UAS flight scenario. Results and verification against high-fidelity simulations will be provided in the final draft.

traffic complexity↗

Methods for design and evaluation of parallel computating systems (The PISCES project)

The PISCES project started in 1984 under the sponsorship of the NASA Computational Structural Mechanics (CSM) program. A PISCES 1 programming environment and parallel FORTRAN were implemented in 1984 for the DEC VAX (using UNIX processes to simulate parallel processes). This system was used for experimentation with parallel programs for scientific applications and AI (dynamic scene analysis) applications. PISCES 1 was ported to a network of Apollo workstations by N. Fitzgerald.

Pratt, Terrence W.↗