Search NASA⌕ Search

SEARCH · Search NASA

Results for “Matlab data”

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 577 records · Page 32

Simulation and Analyses of Multi-Body Separation in Launch Vehicle Staging Environment

The development of methodologies, techniques, and tools for analysis and simulation of multi-body separation is critically needed for successful design and operation of next generation launch vehicles. As a part of this activity, ConSep simulation tool is being developed. ConSep is a generic MATLAB-based front-and-back-end to the commercially available ADAMS. solver, an industry standard package for solving multi-body dynamic problems. This paper discusses the 3-body separation capability in ConSep and its application to the separation of the Shuttle Solid Rocket Boosters (SRBs) from the External Tank (ET) and the Orbiter. The results are compared with STS-1 flight data.

Pamadi, Bandu N.↗

Development of a Piezoelectric Rotary Hammer Drill

The Piezoelectric Rotary Hammer Drill is designed to core through rock using a combination of rotation and high frequency hammering powered by a single piezoelectric actuator. It is designed as a low axial preload, low mass, and low power device for sample acquisition on future missions to extraterrestrial bodies. The purpose of this internship is to develop and test a prototype of the Piezoelectric Rotary Hammer Drill in order to verify the use of a horn with helical or angled cuts as a hammering and torque inducing mechanism. Through an iterative design process using models in ANSYS Finite Element software and a Mason's Equivalent Circuit model in MATLAB, a horn design was chosen for fabrication based on the predicted horn tip motion, electromechanical coupling, and neutral plane location. The design was then machined and a test bed assembled. The completed prototype has proven that a single piezoelectric actuator can be used to produce both rotation and hammering in a drill string through the use of a torque inducing horn. Final data results include bit rotation produced versus input power, and best drilling rate achieved with the prototype.

sampling↗

Stereo-Optic High Definition Imaging: A New Technology to Understand Bird and Bat Avoidance of Wind Turbines

There is a critical need to develop monitoring tools to track aerofauna (birds and bats) in three dimensions around wind turbines. New monitoring systems will reduce permitting uncertainty by increasing the understanding of how birds and bats are interacting with wind turbines, which will improve the accuracy of impact predictions. Biodiversity Research Institute (BRI), The University of Maine Orono School of Computing and Information Science (UMaine SCIS), HiDef Aerial Surveying Limited (HiDef), and SunEdison, Inc. (formerly First Wind) responded to this need by using stereo-optic cameras with near-infrared (nIR) technology to investigate new methods for documenting aerofauna behavior around wind turbines. The stereo-optic camera system used two synchronized high-definition video cameras with fisheye lenses and processing software that detected moving objects, which could be identified in post-processing. The stereo- optic imaging system offered the ability to extract 3-D position information from pairs of images captured from different viewpoints. Fisheye lenses allowed for a greater field of view, but required more complex image rectification to contend with fisheye distortion. The ability to obtain 3-D positions provided crucial data on the trajectory (speed and direction) of a target, which, when the technology is fully developed, will provide data on how animals are responding to and interacting with wind turbines. This project was focused on testing the performance of the camera system, improving video review processing time, advancing the 3-D tracking technology, and moving the system from Technology Readiness Level 4 to 5. To achieve these objectives, we determined the size and distance at which aerofauna (particularly eagles) could be detected and identified, created efficient data management systems, improved the video post-processing viewer, and attempted refinement of 3-D modeling with respect to fisheye lenses. The 29-megapixel camera system successfully captured 16,173 five-minute video segments in the field. During nighttime field trials using nIR, we found that bat-sized objects could not be detected more than 60 m from the camera system. This led to a decision to focus research efforts exclusively on daytime monitoring and to redirect resources towards improving the video post- processing viewer. We redesigned the bird event post-processing viewer, which substantially decreased the review time necessary to detect and identify flying objects. During daytime field trials, we determine that eagles could be detected up to 500 m away using the fisheye wide-angle lenses, and eagle-sized targets could be identified to species within 350 m of the camera system. We used distance sampling survey methods to describe the probability of detecting and identifying eagles and other aerofauna as a function of distance from the system. The previously developed 3-D algorithm for object isolation and tracking was tested, but the image rectification (flattening) required to obtain accurate distance measurements with fish-eye lenses was determined to be insufficient for distant eagles. We used MATLAB and OpenCV to improve fisheye lens rectification towards the center of the image, but accurate measurements towards the image corners could not be achieved. We believe that changing the fisheye lens to rectilinear lens would greatly improve position estimation, but doing so would result in a decrease in viewing angle and depth of field. Finally, we generated simplified shape profiles of birds to look for similarities between unknown animals and known species. With further development, this method could provide a mechanism for filtering large numbers of shapes to reduce data storage and processing. These advancements further refined the camera system and brought this new technology closer to market. Once commercialized, the stereo-optic camera system technology could be used to: a) research how different species interact with wind turbines in order to refine collision risk models and inform mitigation solutions; and b) monitor aerofauna interactions with terrestrial and offshore wind farms replacing costly human observers and allowing for long-term monitoring in the offshore environment. The camera system will provide developers and regulators with data on the risk that wind turbines present to aerofauna, which will reduce uncertainty in the environmental permitting process.

Adams, Evan↗

Highly Available Nuclear Power for Mission-Critical Applications

Some energy consumers require power on an anytime, all-year-round basis with a high level of certainty, including defense installations, isolated communities, and some industrial processes. For these customers, interruptions in electricity or heat can mean substantial financial loss or even loss of life. In the absence of grid-scale energy storage, a high level of power availability can be accomplished only through the robustness and redundancy of power generators. The NuScale small modular reactor design is well suited to provide highly available power because of several features related to both the nuclear steam supply system and the overall plant design. In analogy to Redundant Array of Independent Disks (RAID) systems used to provide highly reliable data storage, a NuScale plant can assure sustained power generation by virtue of its Redundant Array of Integral Reactors (RAIR).This paper describes the NuScale RAIR plant features and summarizes the results of a rigorous analysis of RAIR availability as a function of power or, conversely, the RAIR plant output power as a function of power availability. The analysis utilized the Matrix Laboratory code (MATLAB) and included probability distributions for the frequency and duration of module outages due to planned and unplanned events. Further, the study evaluated the impact of implementing turbine bypass rather than module shutdown and using one or more modules to supply house loads in the case of loss of off-site power (LOOP). Availability results are presented for a 12-module RAIR plant with and without turbine bypass enabled during a LOOP and for different possible connections to the off-site power distribution grid and dedicated service loads. Results indicate that a very high level of availability can be achieved at relatively high power output levels, regardless of turbine bypass and dedicated load connection, compared to the operating fleet.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

SU-F-R-50: Radiation-Induced Changes in CT Number Histogram During Chemoradiation Therapy for Pancreatic Cancer

Purpose: In an effort of early assessment of treatment response, we investigate radiation induced changes in CT number histogram of GTV during the delivery of chemoradiation therapy (CRT) for pancreatic cancer. Methods: Diagnostic-quality CT data acquired daily during routine CT-guided CRT using a CT-on-rails for 20 pancreatic head cancer patients were analyzed. All patients were treated with a radiation dose of 50.4 in 28 fractions. On each daily CT set, the contours of the pancreatic head and the spinal cord were delineated. The Hounsfiled Units (HU) histogram in these contourswere extracted and processed using MATLAB. Eight parameters of the histogram including the mean HU over all the voxels, peak position, volume, standard deviation (SD), skewness, kurtosis, energy, and entropy were calculated for each fraction. The significances were inspected using paired two-tailed t-test and the correlations were analyzed using Spearman rank correlation tests. Results: In general, HU histogram in pancreatic head (but not in spinal cord) changed during the CRT delivery. Changes from the first to the last fraction in mean HU in pancreatic head ranged from −13.4 to 3.7 HU with an average of −4.4 HU, which was significant (P<0.001). Among other quantities, the volume decreased, the skewness increased (less skewed), and the kurtosis decreased (less sharp) during the CRT delivery. The changes of mean HU, volume, skewness, and kurtosis became significant after two weeks of treatment. Patient pathological response status is associated with the changes of SD (ΔSD), i.e., ΔSD= 1.85 (average of 7 patients) for good reponse, −0.08 (average of 6 patients) for moderate and poor response. Conclusion: Significant changes in HU histogram and the histogram-based metrics (e.g., meam HU, skewness, and kurtosis) in tumor were observed during the course of chemoradiation therapy for pancreas cancer. These changes may be potentially used for early assessment of treatment response.

60 APPLIED LIFE SCIENCES↗

Use of a Microphone Phased Array to Determine Noise Sources in a Rocket Plume

A 70-element microphone phased array was used to identify noise sources in the plume of a solid rocket motor. An environment chamber was built and other precautions were taken to protect the sensitive condenser microphones from rain, thunderstorms and other environmental elements during prolonged stay in the outdoor test stand. A camera mounted at the center of the array was used to photograph the plume. In the first phase of the study the array was placed in an anechoic chamber for calibration, and validation of the indigenous Matlab(R) based beamform software. It was found that the "advanced" beamform methods, such as CLEAN-SC was partially successful in identifying speaker sources placed closer than the Rayleigh criteria. To participate in the field test all equipments were shipped to NASA Marshal Space Flight Center, where the elements of the array hardware were rebuilt around the test stand. The sensitive amplifiers and the data acquisition hardware were placed in a safe basement, and 100m long cables were used to connect the microphones, Kulites and the camera. The array chamber and the microphones were found to withstand the environmental elements as well as the shaking from the rocket plume generated noise. The beamform map was superimposed on a photo of the rocket plume to readily identify the source distribution. It was found that the plume made an exceptionally long, >30 diameter, noise source over a large frequency range. The shock pattern created spatial modulation of the noise source. Interestingly, the concrete pad of the horizontal test stand was found to be a good acoustic reflector: the beamform map showed two distinct source distributions- the plume and its reflection on the pad. The array was found to be most effective in the frequency range of 2kHz to 10kHz. As expected, the classical beamform method excessively smeared the noise sources at lower frequencies and produced excessive side-lobes at higher frequencies. The "advanced" beamform routine CLEAN-SC created a series of lumped sources which may be unphysical. We believe that the present effort is the first-ever attempt to directly measure noise source distribution in a rocket plume.

Panda, J.↗

ISSM: Ice Sheet System Model

In order to have the capability to use satellite data from its own missions to inform future sea-level rise projections, JPL needed a full-fledged ice-sheet/iceshelf flow model, capable of modeling the mass balance of Antarctica and Greenland into the near future. ISSM was developed with such a goal in mind, as a massively parallelized, multi-purpose finite-element framework dedicated to ice-sheet modeling. ISSM features unstructured meshes (Tria in 2D, and Penta in 3D) along with corresponding finite elements for both types of meshes. Each finite element can carry out diagnostic, prognostic, transient, thermal 3D, surface, and bed slope simulations. Anisotropic meshing enables adaptation of meshes to a certain metric, and the 2D Shelfy-Stream, 3D Blatter/Pattyn, and 3D Full-Stokes formulations capture the bulk of the ice-flow physics. These elements can be coupled together, based on the Arlequin method, so that on a large scale model such as Antarctica, each type of finite element is used in the most efficient manner. For each finite element referenced above, ISSM implements an adjoint. This adjoint can be used to carry out model inversions of unknown model parameters, typically ice rheology and basal drag at the ice/bedrock interface, using a metric such as the observed InSAR surface velocity. This data assimilation capability is crucial to allow spinning up of ice flow models using available satellite data. ISSM relies on the PETSc library for its vectors, matrices, and solvers. This allows ISSM to run efficiently on any parallel platform, whether shared or distrib- ISSM: Ice Sheet System Model NASA's Jet Propulsion Laboratory, Pasadena, California uted. It can run on the largest clusters, and is fully scalable. This allows ISSM to tackle models the size of continents. ISSM is embedded into MATLAB and Python, both open scientific platforms. This improves its outreach within the science community. It is entirely written in C/C++, which gives it flexibility in its design, and the power/speed that C/C++ allows. ISSM is svn (subversion) hosted, on a JPL repository, to facilitate its development and maintenance. ISSM can also model propagation of rifts using contact mechanics and mesh splitting, and can interface to the Dakota software. To carry out sensitivity analysis, mesh partitioning algorithms are available, based on the Scotch, Chaco, and Metis partitioners that ensure equal area mesh partitions can be done, which are then usable for sampling and local reliability methods.

Larour, Eric↗

Development of a water source heat pump hardware-in-the-loop (HIL) testing facility for smart building applications

Over the last decade, the global fight against climate change through electrification has led to an increase in research on building heating, ventilation, and air conditioning (HVAC) systems that utilize intelligent control algorithms to provide demand-side grid service while maintaining the thermal comfort of building occupants. As the pivotalpoint between building electricity consumption and indoor thermal comfort, high-efficiency electrical heatpumps are at the center of these emerging studies, and various grid-interactive and occupant-comfort control algorithms have been developed for them. The impact of these algorithms on the heatpump operation andperformance under different weather, building load, and grid requests calls for investigation and verification via experimental tests with actual heat pumps integrated with real-time building and grid responses. This study presents a Water-Source Heat Pump Hardware-in-The-Loop (HIL) Test Facility developed with the capability to perform such tests. The hardware configuration for this testfacility introduces a hydronic system that emulates the conditions for the heat pump water-side, and a duct system that emulates conditions for the heat pump airside. Both data acquisition and emulator control are implemented through the National Instruments (NI) LabVIEW software running on an NI PXIplatform. The HIL mechanism based on the hardware-software integration that allows the testbed to communicate with a generic simulation environment is also discussed. Currently, the test facility setup includes a single heatpump and virtual building model in EnergyPlus coupled with an occupant behavioral model in MATLAB. Preliminary test results of the current setup demonstrate the building load emulator's ability to track the simulated gone temperature with a Root Mean Square Deviation (RSME) below 0.12°C (0.216°F). An uncertainty analysis based on sensor accuracies shows that the heat pump coefficient of performance (COP) can be measured with a relative uncertainty of 10.4% in cooling and 3.7% in heating. Apartfrom the current testing on a single heat pump, the test facility also provides the flexibility to include additional heat pumps to form a heat pump cluster, as well as coupling the heat pump with active thermal storage to provide enhanced demandflexibility.

Calfa, Caleb↗

Morphing: A Novel Approach to Astronaut Suit Sizing

The fitting of a spacesuit to an astronaut is an iterative process consisting of two parts. The first uses anthropometric data to provide an approximation of the suit components that will fit the astronaut. The second part is the subjective fitting, where small adjustments are made based on the astronaut s preference. By providing a better approximation of the correct suit components, the entire fit process time can be reduced significantly. The goals of this project are twofold: (1) To evaluate the effectiveness of the existing sizing algorithm for the Mark III Hybrid suit and (2) to determine what additional components are needed in order to provide adequate sizing for the existing astronaut population. A single subject was scanned using a 3D whole-body scanner (VITUS 3D) in the Mark III suit in eight different poses and four subjects in minimal clothing were also scanned in similar poses. The 3D external body scans of the suit and the subject are overlaid and visually aligned in a customized MATLAB program. The suit components were contracted or expanded linearly along the subjects limbs to match the subjects segmental lengths. Two independent measures were obtained from the morphing program on four subjects and compared with the existing sizing information. Two of the four subjects were in correspondence with the sizing algorithm and morphing results. The morphing outcome for a third subject, incompatible with the suit, suggested that an additional arm element at least 6 inches smaller than the existing smallest suit component would need to be acquired. The morphing result of the fourth subject, deemed incompatible with the suit using the sizing algorithm, indicated a different suit configuration which would be compatible. This configuration matched with the existing suit fit check data.

Margerum, Sarah↗

A High-Fidelity Simulation of a Generic Commercial Aircraft Engine and Controller

A new high-fidelity simulation of a generic 40,000 lb thrust class commercial turbofan engine with a representative controller, known as CMAPSS40k, has been developed. Based on dynamic flight test data of a highly instrumented engine and previous engine simulations developed at NASA Glenn Research Center, this non-proprietary simulation was created especially for use in the development of new engine control strategies. C-MAPSS40k is a highly detailed, component-level engine model written in MATLAB/Simulink (The MathWorks, Inc.). Because the model is built in Simulink, users have the ability to use any of the MATLAB tools for analysis and control system design. The engine components are modeled in C-code, which is then compiled to allow faster-than-real-time execution. The engine controller is based on common industry architecture and techniques to produce realistic closed-loop transient responses while ensuring that no safety or operability limits are violated. A significant feature not found in other non-proprietary models is the inclusion of transient stall margin debits. These debits provide an accurate accounting of the compressor surge margin, which is critical in the design of an engine controller. This paper discusses the development, characteristics, and capabilities of the C-MAPSS40k simulation

May, Ryan D.↗

Design and Testing of a Liquid Nitrous Oxide and Ethanol Fueled Rocket Engine

A small-scale, bi-propellant, liquid fueled rocket engine and supporting test infrastructure were designed and constructed at the Energetic Materials Research and Testing Center (EMRTC). This facility was used to evaluate liquid nitrous oxide and ethanol as potential rocket propellants. Thrust and pressure measurements along with high-speed digital imaging of the rocket exhaust plume were made. This experimental data was used for validation of a computational model developed of the rocket engine tested. The developed computational model was utilized to analyze rocket engine performance across a range of operating pressures, fuel-oxidizer mixture ratios, and outlet nozzle configurations. A comparative study of the modeling of a liquid rocket engine was performed using NASA CEA and Cantera, an opensource equilibrium code capable of being interfaced with MATLAB. One goal of this modeling was to demonstrate the ability of Cantera to accurately model the basic chemical equilibrium, thermodynamics, and transport properties for varied fuel and oxidizer operating conditions. Once validated for basic equilibrium, an expanded MATLAB code, referencing Cantera, was advanced beyond CEAs capabilities to predict rocket engine performance as a function of supplied propellant flow rate and rocket engine nozzle dimensions. Cantera was found to comparable favorably to CEA for making equilibrium calculations, supporting its use as an alternative to CEA. The developed rocket engine performs as predicted, demonstrating the developedMATLAB rocket engine model was successful in predicting real world rocket engine performance. Finally, nitrous oxide and ethanol were shown to perform well as rocket propellants, with specific impulses experimentally recorded in the range of 250 to 260 seconds.

33 ADVANCED PROPULSION SYSTEMS↗

Dynamic Characteristic Identification

This paper introduces a dynamic system parameter extraction using experimental data and a known characteristic equation for a single degree of freedom system. The efficacy of curve fitting to determine mass, damping, and stiffness for simulated underdamped, critically damped, and overdamped systems is explored. Two methods are explored to obtain these characteristics using a step response approach in the time domain or a transfer function and dynamic stiffness approach in the frequency domain for simulated systems. Commonly available software, MATLAB, is used for curve fitting. An examination of the effect of inaccurate seeds, or starting values, for mass, damping, and stiffness is included.

Jordan, Clay↗

Modelling of single walled carbon nanotube cylindrical structures with finite element method simulations

In this study, the modulus of elasticity and shear modulus values of single-walled carbon nanotubes SWCNTs were modelled by using both finite element method and the Matlab code. Initially, cylindrical armchair and zigzag single walled 3D space frames were demonstrated as carbon nanostructures. Thereafter, macro programs were written by the Matlab code producing the space truss for zigzag and armchair models. 3D space frames were introduced to the ANSYS software and then tension, compression and additionally torsion tests were performed on zigzag and armchair carbon nanotubes with BEAM4 element in obtaining the exact values of elastic and shear modulus values. In this study, two different boundary conditions were tested and especially used in torsion loading. The equivalent shear modulus data was found by averaging the corresponding values obtained from ten different nodal points on the nanotube path. Finally, in this study it was determined that the elastic constant values showed proportional changes by increasing the carbon nanotube diameters up to a certain level but beyond this level these values remained stable.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Dynamic PET simulator via tomographic emission projection for kinetic modeling and parametric image studies

Purpose: To develop and evaluate a fast and simple tool called dPETSTEP (Dynamic PET Simulator of Tracers via Emission Projection), for dynamic PET simulations as an alternative to Monte Carlo (MC), useful for educational purposes and evaluation of the effects of the clinical environment, postprocessing choices, etc., on dynamic and parametric images. Methods: The tool was developed in MATLAB using both new and previously reported modules of PETSTEP (PET Simulator of Tracers via Emission Projection). Time activity curves are generated for each voxel of the input parametric image, whereby effects of imaging system blurring, counting noise, scatters, randoms, and attenuation are simulated for each frame. Each frame is then reconstructed into images according to the user specified method, settings, and corrections. Reconstructed images were compared to MC data, and simple Gaussian noised time activity curves (GAUSS). Results: dPETSTEP was 8000 times faster than MC. Dynamic images from dPETSTEP had a root mean square error that was within 4% on average of that of MC images, whereas the GAUSS images were within 11%. The average bias in dPETSTEP and MC images was the same, while GAUSS differed by 3% points. Noise profiles in dPETSTEP images conformed well to MC images, confirmed visually by scatter plot histograms, and statistically by tumor region of interest histogram comparisons that showed no significant differences (p < 0.01). Compared to GAUSS, dPETSTEP images and noise properties agreed better with MC. Conclusions: The authors have developed a fast and easy one-stop solution for simulations of dynamic PET and parametric images, and demonstrated that it generates both images and subsequent parametric images with very similar noise properties to those of MC images, in a fraction of the time. They believe dPETSTEP to be very useful for generating fast, simple, and realistic results, however since it uses simple scatter and random models it may not be suitable for studies investigating these phenomena. dPETSTEP can be downloaded free of cost from https://github.com/CRossSchmidtlein/dPETSTEP.

60 APPLIED LIFE SCIENCES↗

Benchmarking Gas Path Diagnostic Methods: A Public Approach

Recent technology reviews have identified the need for objective assessments of engine health management (EHM) technology. The need is two-fold: technology developers require relevant data and problems to design and validate new algorithms and techniques while engine system integrators and operators need practical tools to direct development and then evaluate the effectiveness of proposed solutions. This paper presents a publicly available gas path diagnostic benchmark problem that has been developed by the Propulsion and Power Systems Panel of The Technical Cooperation Program (TTCP) to help address these needs. The problem is coded in MATLAB (The MathWorks, Inc.) and coupled with a non-linear turbofan engine simulation to produce "snap-shot" measurements, with relevant noise levels, as if collected from a fleet of engines over their lifetime of use. Each engine within the fleet will experience unique operating and deterioration profiles, and may encounter randomly occurring relevant gas path faults including sensor, actuator and component faults. The challenge to the EHM community is to develop gas path diagnostic algorithms to reliably perform fault detection and isolation. An example solution to the benchmark problem is provided along with associated evaluation metrics. A plan is presented to disseminate this benchmark problem to the engine health management technical community and invite technology solutions.

Simon, Donald L.↗

Optical Communication Link Atmospheric Attenuation Model

The Space Communications and Navigation (SCaN) Center for Networking, Integration, and Communications (SCENIC) user interface, which provides web accessible space mission simulation and communication system analysis capabilities using verified and validated analysis algorithms, can execute analyses including, but not limited to, line-of-sight, orbit propagation, and dynamic link budget calculations between sets of missions and/or assets. SCENIC's purpose is to provide NASA civil servants and contractors a user-friendly tool, integrated with model data, that can simulate and analyze a range of space mission architectures without the need for repeated and redundant modeling. Given the abundance and further future development of free space optical (FSO) communication channels within modern space infrastructure, the availability of a reliable optical link analysis capability is crucial for SCENIC users. The efforts outlined in this paper aim to provide a model for atmospheric attenuation of FSO communication links, both due to absorption/scattering and turbulence, to increase the accuracy of SCENIC's optical link assessment capabilities. A previous model existed for optical absorption/scattering within the SCaN Link Budget Tool, but it was not location specific for the Earth ground-based nodes, nor was the model optimized for run-time. The new model utilizes years of National Oceanic and Atmospheric Administration (NOAA) visibility data from ground station locations around the world. Visibility, along with the wavelength of the optical signal, are input parameters to calculate the optical specific attenuation, which is a parameter in the calculation of the slant-path attenuation. A final FSO atmospheric attenuation value is comprised of the absorption/scattering attenuation and the turbulence attenuation. A run-time efficient algorithm for the model was then developed and programmed in MATLAB ® . Due to the simple model and vectorization possible in MATLAB, the algorithm has an average run-time of less than one fourth of the run-time of the previous implementation.

Jack L Green↗

Sandia Data Archive (SDA) Specifications (Version 1.1)

The Sandia Data Archive (SDA) format is special implementation the HDF5 (Hierarchal Data Format, version 5) standard. Archive files store data in records identified with a unique text label. Primitive records store numeric, logical, and character arrays of arbitrary size and dimensionality. Compound records store composite MAT- LAB variables--cell arrays, structures, and objects--with arbitrary nesting. External records allow text/binary files to be stored alongside archived data or division of a large file into smaller archives for transmission. This report documents version 1.1 of the SDA standard, which adds support for structure and object arrays. The basic principles of SDA remain unchanged from version 1.0, with minor enhancements and bug fixes. The Sandia Mat lab AnalysiS Hierarchy (SMASH) toolbox uses SDA extensively; support utilities from the toolbox are highlighted here. Although the SDA format is designed around MATLAB, its use of HDF5 allows adoption in other computer languages.

97 MATHEMATICS AND COMPUTING↗

Reconstruction of the 1994 Pittsburgh Airplane Accident Using a Computer Simulation

On September 8, 1994, a Boeing 737-300 passenger airplane was on a downwind approach to the Pittsburgh International Airport at an altitude of 5000 feet above ground level (6000 feet MSL). While in a shallow left turn onto a downwind approach heading, the airplane crossed into the vortex trail of a Boeing 727 flying in the same approach pattern about 4 miles ahead. The B-737 airplane rolled and turned sharply to the left, exited the vortex wake and plunged into the ground. Weather was not a factor in the accident. The airplane was equipped with a 11+ channel digital Flight Data Recorder (FDR) and a multiple channel Cockpit Voice Recorder (CVR). Both recorders were recovered from the crash site and provided excellent data for the development of an accident scenario. Radar tracking of the two airplanes as well as the indicated air speed (IAS) perturbations clearly visible on the B-737 FDR recordings indicate that the upset was apparently initiated by the airplane's crossing into the wake of the B-727 flying ahead in the same traffic pattern. A 6 degree-of-freedom simulation program for the B-737 airplane using MATLAB and SIMULINK was constructed. The simulation was initialized at the stabilized flight conditions of the airplane about 13 seconds prior to its entry into the vortex trail of the B-727 airplane. By assuming a certain combination of control inputs, it was possible to produce a simulated motion that closely matched that recorded on the FDR.

Parks, Edwin K.↗