Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Simulation and Modeling Tools”

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 343 records · Page 19

Modeling and Control Design for a Turboelectric Single Aisle Aircraft Propulsion System

A nonlinear dynamic model with full flight envelope controller is developed for the propulsion system of a partially turboelectric single-aisle aircraft. The propulsion system model consists of two turbofan engines with a large percentage of power extraction, feeding an electric tail fan for boundary layer ingestion. The dynamic model is compared against an existing steady state design model. An electrical system model using a simple power flow approach is integrated into existing modeling tools used for dynamic simulation of the turbomachinery of the vehicle. In addition to the simple power flow model of the electrical system, a more detailed model is used for comparison at a key vehicle transient flight condition. The controller is a gain scheduled proportional-integral type that is examined throughout the flight envelope for performance metrics such as rise time and operability margins. Potential improvements in efficiency for the vehicle are explored by adjusting the power split between the energy used for thrust by the turbofans and that extracted to supply power to the tail fan. Finally, an operability study of the vehicle is conducted using a 900 nautical mile mission profile for a nominal vehicle configuration, a deteriorated propulsion system at the end of its operating life, and an optimized power schedule with improved efficiency.

Turboelectric↗

Model-free estimation of completeness, uncertainties, and outliers in atomistic machine learning using information theory

Abstract An accurate description of information is relevant for a range of problems in atomistic machine learning (ML), such as crafting training sets, performing uncertainty quantification (UQ), or extracting physical insights from large datasets. However, atomistic ML often relies on unsupervised learning or model predictions to analyze information contents from simulation or training data. Here, we introduce a theoretical framework that provides a rigorous, model-free tool to quantify information contents in atomistic simulations. We demonstrate that the information entropy of a distribution of atom-centered environments explains known heuristics in ML potential developments, from training set sizes to dataset optimality. Using this tool, we propose a model-free UQ method that reliably predicts epistemic uncertainty and detects out-of-distribution samples, including rare events in systems such as nucleation. This method provides a general tool for data-driven atomistic modeling and combines efforts in ML, simulations, and physical explainability.

36 MATERIALS SCIENCE↗

Core Modeling & Decision Support Capabilities: FASTSim, RouteE, T3CO & OpenPATH

These are the slides on NREL's Core Modeling and Decision Support Capabilities project for presentation at the 2023 Annual Merit Review for the U.S. Department of Energy, Vehicle Technologies Office. The project is funded by the Energy-Efficient Mobility Systems program to support four core tools: the Future Automotive Systems Technology Simulator (FASTSim), the Route Energy Prediction Model (RouteE), the Transportation Technology Total Cost of Ownership (T3CO) tool, and the Open Platform for Agile Trip Heuristics (OpenPATH).

ADVANCED PROPULSION SYSTEMS↗

Adaptation of G-TAG Software for Validating Touch-and-Go Comet Surface Sampling Design Methodology

The G-TAG software tool was developed under the R&TD on Integrated Autonomous Guidance, Navigation, and Control for Comet Sample Return, and represents a novel, multi-body dynamics simulation software tool for studying TAG sampling. The G-TAG multi-body simulation tool provides a simulation environment in which a Touch-and-Go (TAG) sampling event can be extensively tested. TAG sampling requires the spacecraft to descend to the surface, contact the surface with a sampling collection device, and then to ascend to a safe altitude. The TAG event lasts only a few seconds but is mission-critical with potentially high risk. Consequently, there is a need for the TAG event to be well characterized and studied by simulation and analysis in order for the proposal teams to converge on a reliable spacecraft design. This adaptation of the G-TAG tool was developed to support the Comet Odyssey proposal effort, and is specifically focused to address comet sample return missions. In this application, the spacecraft descends to and samples from the surface of a comet. Performance of the spacecraft during TAG is assessed based on survivability and sample collection performance. For the adaptation of the G-TAG simulation tool to comet scenarios, models are developed that accurately describe the properties of the spacecraft, approach trajectories, and descent velocities, as well as the models of the external forces and torques acting on the spacecraft. The adapted models of the spacecraft, descent profiles, and external sampling forces/torques were more sophisticated and customized for comets than those available in the basic G-TAG simulation tool. Scenarios implemented include the study of variations in requirements, spacecraft design (size, locations, etc. of the spacecraft components), and the environment (surface properties, slope, disturbances, etc.). The simulations, along with their visual representations using G-View, contributed to the Comet Odyssey New Frontiers proposal effort by indicating problems and/or benefits of different approaches and designs.

Mandic, Milan↗

Development of a Sizing and Modeling Platform for District Energy Systems with Geothermal Heat Pumps

Existing tools for community or urban scale energy system modeling and simulation are often limited in their capabilities and require expert-level modeling proficiency to develop system models. To fill this gap, this paper proposes an integrated sizing and modeling platform for district energy systems with geothermal heat pumps. The proposed platform takes in geometric and non-geometric user inputs related to the buildings, borefield, and district energy loop. Then, the platform sizes the geothermal heat exchanger, generates a corresponding district energy system model, and runs an annual simulation automatically. We validated the simulation performance of the borefield in our tool against EnergyPlus. A case study is provided in this paper to demonstrate the workflow and simulation result plausibility of the proposed platform.

decarbonization↗

The Circulation of the Plasmasphere Fluid during the Erosion Event on September 8, 2017

A strong solar wind pressure pulse triggered the magnetic storm on September 7-8, 2017. Near the end of September 7, the z-component of the IMF Bz dropped from 9 to -10 nT in 30 min and then had another rapid drop to -31 nT in 30 min. The sudden plunge of IMF Bz and the associated strong convection electric field stirred up the storm main phase with Dst falling from 0 to -122 nT from 2200 UT on September 7 to 0200 UT on the 8th. Severe plasmasphere erosion was observed on September 8 by multiple spacecraft, such as the Van Allen Probes and the Arase satellite. In this study, we examine the fate of the eroded plasmasphere particles during this event by model simulation as well as satellite data analysis. The simulation tool we use is the SWMF/BATS-R-US model coupled with the CIMI model. The cold plasmas in the plasmasphere is treated as a separate fluid in the MHD equations. As a result, the transport and circulation of the plasmasphere fluid in the global magnetosphere can be traced and the impacts of this cold fluid on the global magnetosphere can be evaluated. In this paper, we will show how the drainage plume is formed during the storm and how the plasmasphere fluid is transported to the flank and lobe regions and eventually to the plasma sheet and reenters into the inner magnetosphere. We will validate our simulation by plasmasphere signatures observed in both the inner and outer magnetosphere.

Mei-Ching Fok↗

CONFIG: Qualitative simulation tool for analyzing behavior of engineering devices

To design failure management expert systems, engineers mentally analyze the effects of failures and procedures as they propagate through device configurations. CONFIG is a generic device modeling tool for use in discrete event simulation, to support such analyses. CONFIG permits graphical modeling of device configurations and qualitative specification of local operating modes of device components. Computation requirements are reduced by focussing the level of component description on operating modes and failure modes, and specifying qualitative ranges of variables relative to mode transition boundaries. Simulation processing occurs only when modes change or variables cross qualitative boundaries. Device models are built graphically, using components from libraries. Components are connected at ports by graphical relations that define data flow. The core of a component model is its state transition diagram, which specifies modes of operation and transitions among them.

Malin, Jane T.↗

BDEM (Discrete-element-simulator for high-solids granular flows) [SWR-22-72]

BDEM is a discrete element method based simulation tool developed specifically for modeling high-solids granular flows that include polydispersity, heat-transfer, moving boundaries and chemistry. Our solver provides facilities for simulating spherical/non-spherical particles with modified contact and friction models in complex dynamic geometries defined using level-sets or triangulated files. The solver is developed on top of NREL's open-source performance portable library, AMReX, providing parallel execution capabilities on current and upcoming high-performance-computing (HPC) architectures. Simulations at the scale of several millions to billion particles have been performed using this software on large scale computing resources. This software can be applied to non-reacting solids dominant flows in silos, hoppers and screw conveyors as well as in high temperature reacting systems such as screw kilns and auger reactors.

Sitaraman, Hariswaran↗

Computer system simulation and analysis

A simulation model has been developed as a tool for studying different computer system configurations. The model simulates a computer system in terms of the number and speed of hardware subsystems and the runs processed by the system in terms of equivalent adds for computation and data bits transferred for I/O. The model was written in the GPSS simulation language.

Williams, T. G.↗

A Survey of Modeling Activities by NASA’s Orbital Debris Program Office

NASA’s Orbital Debris Program Office (ODPO) develops and maintains a number of modeling tools to analyze and simulate the orbital debris environment. One of the most important products produced by the NASA ODPO is the Orbital Debris Engineering Model (ORDEM). This model can be used by satellite designers and operators to design missions for better protection against the debris environment. The ODPO is currently working on the next generation, designated ORDEM 4.0. ORDEM 4.0 will include many known features from previous models, such as the ability to input a spacecraft orbit and time and the ability to compute the flux as a function of debris size, impact speed, impact direction, and debris material densities, as well as uncertainty information on the flux. A new addition will be a parameterized debris shape model based on laboratory hypervelocity impact tests, including DebriSat. ORDEM is primarily based on dedicated debris measurements, such as by the Haystack Ultrawideband Satellite Imaging Radar (HUSIR), NASA’s Goldstone radar, and observations of geosynchronous orbits (GEO) using the Eugene Stansbery-Meter Class Autonomous Telescope (ES-MCAT). In addition to ORDEM, the ODPO also maintains other models, such as the LEO-to-GEO Environment Debris (LEGEND) model for studies of long-term evolution of Earth’s debris environment, with the ability to study various mitigation and remediation strategies. Another model, the Satellite Breakup Risk Assessment Model (SBRAM), is used to analyze how satellite breakups may affect critical space missions (such as the ISS) on short notice. In addition to these models, the ODPO maintains other secondary models used to model satellite explosions and collisions, analyze radar, optical, and in situ data, and to model such things as solar activity and orbit evolution. In this presentation, a survey of these models will be presented, showing how the different models are used together to create a comprehensive picture of Earth’s debris environment.

Mark Matney↗

Infrastructure improvements in the National Reactor Innovation Center Virtual Test Bed

The National Reactor Innovation Center’s (NRIC) mission is to support deployment of novel reactor concepts. This is achieved by providing physical and virtual spaces for building and testing various components, systems, and complete pilot plants. The Virtual Test Bed (VTB) represents the virtual counterpart to the physical test bed. The VTB is being developed in collaboration with the Department of Energy’s (DOE) Nuclear Energy Advanced Modeling and Simulation (NEAMS) program. The mission of the VTB is to accelerate the deployment and licensing of advanced reactors by leveraging state-of-the-art modeling and simulation (M&S) tools developed by the DOE NEAMS program. This is accomplished via three primary means: (1) openly hosting simulations that showcase analysis capabilities, (2) continuously testing the models hosted against code updates to avoid deprecation, and (3) filling key M&S gaps that are relevant for the physical NRIC test beds. The VTB repository consists of two sub-entities: • A documentation website detailing the models: https://mooseframework.inl.gov/virtual_test_bed • A GitHub repository that hosts the corresponding files: https://github.com/idaholab/virtual_test_bed. This paper presents the infrastructure added to the VTB in the last two fiscal years. Additional information about the VTB can be found in various publications [1, 2, 3, 4, 5, 6, 7, 8].

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗

Status of New Models Hosted on the Virtual Test Bed (VTB) in 2024

The National Reactor Innovation Center (NRIC) mission is to support deployment of novel reactor concepts. This is achieved by providing physical and virtual spaces for building and testing various components, systems, and complete pilot plants. The Virtual Test Bed (VTB) represents the virtual counterpart to the physical test bed. It is in development in collaboration with the Department of Energy’s (DOE) Nuclear Energy Advanced Modeling and Simulation (NEAMS) program. The mission of the VTB is to accelerate the deployment and licensing of advanced reactors by leveraging state-of-the-art modeling and simulation (M&S) tools developed by the DOE NEAMS program. This is accomplished by three primary means: (1) openly hosting simulations that showcase analysis capabilities, (2) continuously testing the models hosted against code updates to avoid deprecation, and (3) filling key M&S gaps that are relevant for the physical NRIC test beds. The VTB repository consists of two sub-entities: 1. A documentation website detailing the models (https://mooseframework.inl.gov/virtual_test_bed). 2. A GitHub repository that hosts the corresponding files (https://github.com/idaholab/virtual_test_bed). Previous documentation on the models hosted in the VTB can be found in [1,2,3,4]. These references also include additional background information on the various NEAMS codes showcased in the VTB (which is omitted here for brevity). This paper primarily provides a status update of the most recent additions to the repository.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗

Impact of Time-Dependent Reactor and Sensor Physics on Core Power Synthesis (Rev.1)

Online synthesis of power distribution is critical in the operation and control of nuclear power reactors to ensure that the core is operating within safety margins and to provide essential knowledge associated with the burnup of the fuel. In light-water reactors, power synthesis is achieved by using some a priori knowledge of the state of the reactor core and updating based on the signals coming from in-core sensors—namely, self-powered neutron detectors (SPNDs). This report examines the effects of fuel burnup and sensor degradation on the ability to accurately synthesize the power distribution in a pressurized water reactor (PWR), considering the typical low-enriched uranium (LEU, 3%-5% enrichment) fuel cycle as well as the higher enrichment LEU+ (5%-8% enrichment) fuel cycle. Several modeling tools were used to simulate power synthesis based on the responses of SPNDs, with emitters made out of Rh or V. A representative PWR LEU core was modeled using the Polaris/Purdue Advanced Reactor Core Simulator (PARCS) approach. The Monte Carlo N-Particle Transport 6 (MCNP6) code was used as well to calculate response functions between different segments of fuel to individual SPNDs; this is a crucial parameter for power synthesis. The Oak Ridge Isotope GENeration (ORIGEN) package in the Standardized Computer Analyses for Licensing Evaluation (SCALE) code was used to model the time-dependent isotopic transmutation in the SPND emitters. All these data were fed into a custom code that enacted the point-based iterative method to simulate power synthesis. Developmental work was also performed on high-fidelity SPND models in the GEometry ANd Tracking 4 (Geant4) code, which enables higher-accuracy modeling of the current responses from SPNDs. In this work, five sets of time-dependent power synthesis test cases were conducted. In these test cases, systematic changes in the input conditions enabled an analysis of the effect of (1) slightly inaccurate a priori power distribution assumptions with respect to fuel burnup, (2) highly inaccurate a priori power distribution assumptions with respect to fuel burnup (such that burnup is not included in the a priori assumed distribution), and (3) differences between Rh and V SPNDs in terms of downstream consequences of the transmutation in the emitters and the extended nature of the LEU+ fuel cycle in comparison with LEU. The authors discovered that one may permissibly have slightly inaccurate a priori assumptions of the fuel burnup (such that the level of burnup may be slightly underapproximated or overapproximated by the accumulated burnup in approximately 9.3 full power days), but to not account for burnup at all in the a priori assumptions leads to severe levels of error, approaching 25% at maximum (for LEU). The authors also discovered that V SPNDs are extraordinarily robust in both the LEU and LEU+ fuel cycles considered in this modeling work, whereas Rh SPNDs undergo significant transmutation that can result in large errors in the synthesized power distribution.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Quantum fluctuations in dense plasma simulations

Molecular dynamics (MD) simulations are a powerful tool for modeling warm and hot dense matter. Density functional theory (DFT) MD simulations are often preferred in dense plasmas in order to accurately model quantum electronic structure. However, DFT-MD simulations neglect interaction effects due to fluctuations in excited states. In this work, we present an MD approach that uses excited state method pseudoatoms to run dense plasma simulations with many different core-electron configurations at classical MD speeds. We also allow for transitions between different configurations in our simulations and find that these fluctuations are especially important for highly excited states. Our results suggest that finite configuration lifetimes that are comparable to the inverse ion plasma frequency need to be accounted for in order to accurately model ion distributions in dense plasma simulations. We also demonstrate that excited state fluctuations have a direct impact on ion plasma microfields, generate different plasma microfields for different excitation levels, and thereby induce absorption–emission line shape asymmetries even in steady-state plasmas.

36 MATERIALS SCIENCE↗

EMU Suit Performance Simulation

Introduction: Designing a planetary suit is very complex and often requires difficult trade‐offs between performance, cost, mass, and system complexity. To verify that new suit designs meet requirements, full prototypes must be built and tested with human subjects. However, numerous design iterations will occur before the hardware meets those requirements. Traditional draw‐prototype‐test paradigms for research and development are prohibitively expensive with today's shrinking Government budgets. Personnel at NASA are developing modern simulation techniques that focus on a human‐centric design paradigm. These new techniques make use of virtual prototype simulations and fully adjustable physical prototypes of suit hardware. This is extremely advantageous and enables comprehensive design down‐selections to be made early in the design process. Objectives: The primary objective was to test modern simulation techniques for evaluating the human performance component of two EMU suit concepts, pivoted and planar style hard upper torso (HUT). Methods: This project simulated variations in EVA suit shoulder joint design and subject anthropometry and then measured the differences in shoulder mobility caused by the modifications. These estimations were compared to human‐in‐the‐loop test data gathered during past suited testing using four subjects (two large males, two small females). Results: Results demonstrated that EVA suit modeling and simulation are feasible design tools for evaluating and optimizing suit design based on simulated performance. The suit simulation model was found to be advantageous in its ability to visually represent complex motions and volumetric reach zones in three dimensions, giving designers a faster and deeper comprehension of suit component performance vs. human performance. Suit models were able to discern differing movement capabilities between EMU HUT configurations, generic suit fit concerns, and specific suit fit concerns for crewmembers based on individual anthropometry

Cowley, Matthew S.↗

Supporting Hazard Analysis for Wildfire Response Using fmdtools and MIKA

The System Wide Safety (SWS) Safety Demonstrator (SD) Series drives development of an increasingly capable In-Time Aviation Safety Management System (IASMS) focusing on humanitarian applications, starting with wildfire response (SD-1). The goals of this report are to (1) provide an early hazard analysis and mitigation evaluation of wildfire response to support these efforts and (2) provide a demonstration of capabilities of the Fault Model Design Tools (fmdtools) and Manager for Intelligent Knowledge Access (MIKA) tools. fmdtools provides a modeling, simulation, and resiliency analysis framework in which a wildfire response model, the System Modeling and Analysis of Resiliency in Scalable Traffic Management for Emergency Response Operations (SMARt-STEReO), is built. MIKA is an intelligent knowledge manager with several capabilities, including assisting in hazard analysis by extracting and analyzing hazards from historical incident reports. The following topics are covered in the report: Understanding Wildfire Hazard Dynamics. We provide a description and simulated examples of how hazards occur in the SMARt-STEReO model of wildfire response and their effect on its outcome. This provides a common mental model and focuses the analysis presented in the remainder of the report. Wildfire Hazard Identification. MIKA identifies wildfire hazards from three relevant datasets: the ICS-209-PLUS, SAFECOM, and SAFENET. Hazards are manually organized into a taxonomy and MIKA analyzes each hazard’s effects, likelihood, severity, and risk. Evaluating Mitigation Strategies. The SMARt-STEReO wildfire response model built in fmdtools evaluates a subset of identified hazards. Specifically, we simulate the effect of communications faults and equipment faults on operator safety, the effect of changing winds and flammability, and a scenario with multiple ignition points and heavy smoke. Tool Limitations and Usage Considerations. We provide a discussion of appropriate tool use cases as well as limitations and considerations for usage. The tool findings are used to synthesize recommendations for wildfire response operations, which can be captured as part of an IASMS. Key recommendations are as follows: Hazards are identified from a broad spectrum of sources including aircraft subsystems, operational sources, and ground crew operations. Highest risk operational environment hazards identified are Evacuations. The highest risk manned aerial operations hazard categorized is Jumper Operations Mishap. Ground crew hazards that are highest risk are Burns, Cargo Operations Overhead, Dehydration, Entrapment, Falling Objects, Heart Attacks, Heat Exhaustion, Inadequate Training or Certification, Vehicle Breakdown, and Vehicle Collision. Modelled containment failures arise from a mismatch between the difficulty of the firefighting scenario and the capacity (e.g., speed, effectiveness, awareness) of the response. In firefighting scenarios where containment is possible (e.g., because the fire does not spread too quickly), these mismatches can occur because of a change in environmental conditions (e.g., wind, flammability, etc) or because of planning, equipment, or communications faults. Improvements to communications increase the capacity of the firefighting response by reducing the time needed to respond to the fire. While surveillance does not increase this capacity by itself, it increases operator safety by increasing state awareness, enabling firefighters to evade approaching fires. Increasing both has a synergistic effect. In general, these performance and resilience increases generalize over fault scenarios as well as unforeseen changes to circumstances (i.e., wind, aridity, etc.). However, these improvements need to be designed so as not to make the system prone to persistent large-scale communications outages, which can reduce performance.

Hazard analysis↗

Impact of Time Dependent Reactor and Sensor Physics on Core Power Synthesis

Online synthesis of the power distribution is critical in the operation and control of nuclear power reactors to ensure that the core is operating within safety margins, and to provide essential knowledge associated with the burnup of the fuel. In light water reactors (LWRs), power synthesis is achieved by using some a priori knowledge of the state of the reactor core and updating based on the signals coming from in-core sensors—namely, self-powered neutron detectors (SPNDs). This report aims to study the effects of fuel burnup and sensor degradation on the ability to accurately synthesize the power distribution in a LWR. Several modeling tools were used to simulate power synthesis based on the responses of SPNDs, with emitters made out of Rh or V. A representative pressurized water reactor low-enriched uranium (LEU) core was modeled using the Polaris/Purdue Advanced Reactor Core Simulator (PARCS) approach. The Monte Carlo N-Particle Transport 6 (MCNP6) code was used, as well, to calculate response functions between different segments of fuel to individual SPNDs; this is a crucial parameter for power synthesis. The Oak Ridge Isotope GENeration (ORIGEN) package in the Standardized Computer Analyses for Licensing Evaluation (SCALE) code was used to model the time-dependent isotopic transmutation in the SPND emitters. All these data were fed into a custom code that enacted the point-based iterative (PBI) method to simulate power synthesis. Developmental work was also performed on high-fidelity SPND models in the GEometry ANd Tracking 4 (Geant4) code, which enables higher-accuracy modeling of the current responses from SPNDs. In this work, five sets of time-dependent power synthesis test cases were conducted. In these test cases, systematic changes in the input conditions enabled an analysis of the effect of (1) slightly inaccurate a priori power distribution assumptions with respect to fuel burnup, (2) highly inaccurate a priori assumptions with respect to fuel burnup (such that burnup is not included in the a priori assumed distribution), and (3) differences between Rh and V SPNDs in terms of downstream consequences of the transmutation in the emitters. The authors discovered that one may permissibly have slightly inaccurate a priori assumptions of the fuel burnup (such that the level of burnup may be slightly under- or over-approximated by the accumulated burnup in approximately 9.3 full power days), but to not account for burnup at all in the a priori assumption leads to severe levels of error, approaching 25% at maximum. The authors also discovered that V SPNDs are extraordinarily robust in the low-enriched uranium fuel cycle considered in this modeling work, whereas Rh SPNDs undergo significant transmutation that can result in large errors in the synthesized power distribution.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Simulation Development and Analysis of Crew Vehicle Ascent Abort

NASA's Commercial Crew Program is an integral step in its journey to Mars as it would expedite development of space technologies and open up partnership with U.S. commercial companies. NASA reviews and independent assessment of Commercial Crew Program is fundamental to its success, and being able to model a commercial crew vehicle in a simulation rather than conduct a live test would be a safer, faster, and less expensive way to assess and certify the capabilities of the vehicle. To this end, my project was to determine the feasibility of using a simulation tool named SOMBAT version 2.0 to model a multiple parachute system for Commercial Crew Program simulation. The main tasks assigned to me were to debug and test the main parachute system model, (capable of simulating one to four main parachute bodies), and to utilize a graphical program to animate the simulation results. To begin tackling the first task, I learned how to use SOMBAT by familiarizing myself with its mechanics and by understanding the methods used to tweak its various parameters and outputs. I then used this new knowledge to set up, run, and analyze many different situations within SOMBAT in order to explore the limitations of the parachute model. Some examples of parameters that I varied include the initial velocity and orientation of the falling capsule, the number of main parachutes, and the location where the parachutes were attached to the capsule. Each parameter changed would give a different output, and in some cases, would expose a bug or limitation in the model. A major bug that I discovered was the inability of the model to handle any number of parachutes other than three. I spent quite some time trying to debug the code logically, but was unable to figure it out until my mentor taught me that digital simulation limitations can occur when some approximations are mistakenly assumed for certain in a physical system. This led me to the realization that unlike in all of the programming classes I have taken thus far that focus on pure logic, simulation code focuses on mimicking the physical world with some approximation and can have inaccuracies or numerical instabilities. Learning from my mistake, I adopted new methods to analyze these different simulations. One method the student used was to numerically plot various physical parameters using MATLAB to confirm the mechanical behavior of the system in addition to comparing the data to the output from a separate simulation tool called FAST. By having full control over what was being outputted from the simulation, I could choose which parameters to change and to plot as well as how to plot them, allowing for an in depth analysis of the data. Another method of analysis was to convert the output data into a graphical animation. Unlike the numerical plots, where all of the physical components were displayed separately, this graphical display allows for a combined look at the simulation output that makes it much easier for one to see the physical behavior of the model. The process for converting SOMBAT output for EDGE graphical display had to be developed. With some guidance from other EDGE users, I developed a process and created a script that would easily allow one to display simulations graphically. Another limitation with the SOMBAT model was the inability for the capsule to have the main parachutes instantly deployed with a large angle between the air speed vector and the chutes drag vector. To explore this problem, I had to learn about different coordinate frames used in Guidance, Navigation & Control (J2000, ECEF, ENU, etc.) to describe the motion of a vehicle and about Euler angles (e.g. Roll, Pitch, Yaw) to describe the orientation of the vehicle. With a thorough explanation from my mentor about the description of each coordinate frame, as well as how to use a directional cosine matrix to transform one frame to another, I investigated the problem by simulating different capsule orientations. In the end, I was able to show that this limitation could be avoided if the capsule is initially oriented antiparallel to its velocity vector.

Wong, Chi S.↗