Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Direct coupling analysis”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 505 records · Page 28

Mineralogical, Elemental and Tomographic Reconnaissance Investigation for CLPS (“METRIC”)

METRIC comprises a suite of two instruments in a Lunar lander that perform X-ray diffraction (XRD) for mineral structure, X-ray fluorescence (XRF) for elemental composition, and X-ray Micro Computed Tomography (XCT) for 3D internal micromorphology. The instruments are accompanied by optical-near IR cameras to provide local geologic context. The Honeybee Robotics PlanetVac pneumatic sampling and transfer system [1],positioned on a lander footpad, will deliver sieved regolith to the X-ray instruments for analysis. The instrument suite is intended for delivery to the lunar surface on a Commercial Lunar Payload Services (CLPS) lander. The METRICXRD/F instrument draws on heritage from the Mars Science Laboratory CheMin instrument [2] and improves upon the design in multiple ways [3]. Like CheMin, Rietveld refinement and full-pattern fitting of METRIC XRD data can identify minerals at a detection limit of ~1 wt.%, quantify their abundances when present at >3 wt.%, and determine mineral composition (e.g., Fo#in olivine) from lattice parameters for minerals present at >5 wt.%[4,5]. An optimized XRF geometry provides for improved detection and quantification of major, minor and trace elements. The METRIC XCT instrument is a miniaturized X-ray Computed Tomography scanner[6]. Lunar regolith is delivered to a 3 mm diameter, 10 mm long graphite tube inside the instrument. The tube is rotated through 360°in 0.9deg.incrementsand a divergent point source X-ray beam is directed through the material. A Charge Coupled Device (CCD)records attenuation images whose brightness and contrast area function of average atomic number and density. Quantitative data, including particle and void sizes, 3D particle shape parameters, modal volumes and pore geometry can be derived from the resulting 3D reconstructions(voxel resolution: 30 μm). Crystal morphologies derived from METRIC XCT data complement the bulk mineralogy determined by the METRIC XRD/F and provide a measure of grain size distribution for the different phases. Taken together, the METRIC instrument suite determines crystal structure, elemental composition and morphology, three principal characteristics of geological materials that are highly useful in determining the origin and subsequent processing of lunar regolith.

Moon↗

A Technology Pathway for Airbreathing, Combined-Cycle, Horizontal Space Launch Through SR-71 Based Trajectory Modeling

Access to space is in the early stages of commercialization. Private enterprises, mainly under direct or indirect subsidy by the government, have been making headway into the LEO launch systems infrastructure, of small-weight-class payloads of approximately 1000 lbs. These moderate gains have emboldened the launch industry and they are poised to move into the middle-weight class (roughly 5000 lbs). These commercially successful systems are based on relatively straightforward LOX-RP, two-stage, bi-propellant rocket technology developed by the government 40 years ago, accompanied by many technology improvements. In this paper we examine a known generic LOX-RP system with the focus on the booster stage (1st stage). The booster stage is then compared to modeled Rocket-Based and Turbine-Based Combined Cycle booster stages. The air-breathing propulsion stages are based on/or extrapolated from known performance parameters of ground tested RBCC (the Marquardt Ejector Ramjet) and TBCC (the SR-71/J-58 engine) data. Validated engine models using GECAT and SCCREAM are coupled with trajectory optimization and analysis in POST-II to explore viable launch scenarios using hypothetical aerospaceplane platform obeying the aerodynamic model of the SR-71. Finally, and assessment is made of the requisite research technology advances necessary for successful commercial and government adoption of combined-cycle engine systems for space access.

Kloesel, Kurt J.↗

Microscopic and Metallurgical Aspects of the Space Shuttle Columbia Accident Investigation and Reconstruction

The Space Shuttle Columbia was descending for a landing at the Kennedy Space Center (KSC) on February 1, 2003. Approximately 20 minutes prior to touchdown, the Columbia began disintegrating over the western United States; the majority of debris eventually impacted in eastern Texas and western Louisiana. A monumental effort eventually recovered approximately 84,000 pieces of debris, approximately 38% of the Orbiter's original dry weight. The debris was transported to KSC, where the items were catalogued and evaluated. Critical areas of interest, such as the left and right leading edge surfaces and the underside of the ship, were placed upon a grid to aid in the reconstruction. Items of interest included metallic structures, reinforced carbon-carbon composites, and ceramic heat insulation tiles. Many of the leading edge elements had re-solidified metallic deposits spattered on them. These deposits became known as slag and were one of the main focuses of the investigation. In order to help determine the sequence of events inside the left wing during the accident, the slag's composition, layering order, and directionality of deposition were studied. A myriad of analytical tests were performed in an attempt to ascertain the compositional and depositional characteristics of selected slag deposits, including the ordering of deposited layers within each individual slag deposit harvested. Initially, Scanning Electron Microscopy and Energy Dispersive X-Ray Spectroscopy (SEM/EDX) were performed to quickly characterize the overall composition of individual slag deposits: SEM utilizes a narrowlyfocused high-energy electron beam impinging upon a specimen. The incident beam excites and liberates lower energy secondary electrons, which are detected and analyzed, providing a visual representation of the sample's surface topography. EDX also relies on an incident electron beam, except an EDX unit measures X-ray energies generated by the impinging beam. Each element generates a unique X-ray signature; the EDX detector measures these discreet energies. EDX actually penetrates approximately 2 microns into the bulk of the sample. However, random examination of various portions of slag, coupled with the semiquantitative nature of the SEM/EDX analysis, did not yield convincingly pertinent data. Therefore, X-ray dot mapping was conducted, which provided more understandable data, both in terms of slag layering and composition. An X-ray dot map is generated by performing numerous EDX scans for individual elements, then compiling the scans in a visual representation. Eventually, specimens consisting of not only the slag, but of the adjacent RCC substrate as well were cross-sectioned. X-Ray dot mapping of the materialographicallymounted and -polished cross- sections provided a visual representation of both the layering sequence and compositional characteristics of the slag. Contemporaneously, Electron Spectroscopy for Chemical Analysis/X-Ray Photoelectron Spectroscopy (ESCA/XPS) and powdered X-Ray Diffraction (XRD) were performed to further characterize the deposits and to attempt to identify what, if any, compounds were present. The ESCA/XPS analysis allowed the analyst to "sputter" into the sample with an electron gun, aiding in the identification of the layering sequence. XPS uses photons, rather than electrons, which impinge upon the surface of the sample. XPS measures the electrons emitted from within the first 5 nm of the sample's surface. The XRD measures the scatter angles of incident X-rays; the angle and intensity of scatter depend upon the crystalline structure of the pulverized sample. XRD is considered a qualitative rather than quantitative technique. ESCA/XPS revealed that the final layer to deposit was predominantly carbonaceous. XRD was successful in identifying specific compounds, such as Al 2O3, Al and/or Al3 21SiO47, mullite (3(Al2)O3 -SiO2), and nickel-aluminides. Eventually, Electron MicroProbe Analysis (EMPA) was conducted on the marialographically-prepared cross- sections of selected slag deposits. Microprobe combines SEM and Wavelength Dispersive X-Ray Spectroscopy (WDS), and, like EDX, uses a narrowly-focused high-energy electron beam impinging upon a specimen to elicit, in the case of EPMA, characteristic X-rays with specific wavelengths. This quantitative, analytical tool proved the most useful in determining depositional layering and composition of the slag deposits. This information was utilized in verifying the location of the breach in the left leading edge of the wing of the Columbia.

McDaniels, Steven J.↗

Infrared Spectroscopy as a Chemical Fingerprinting Tool

Infrared (IR) spectroscopy is a powerful analytical tool in the chemical fingerprinting of materials. Any sample material that will interact with infrared light produces a spectrum and, although normally associated with organic materials, inorganic compounds may also be infrared active. The technique is rapid, reproducible and usually non-invasive to the sample. That it is non-invasive allows for additional characterization of the original material using other analytical techniques including thermal analysis and RAMAN spectroscopic techniques. With the appropriate accessories, the technique can be used to examine samples in liquid, solid or gas phase. Both aqueous and non-aqueous free-flowing solutions can be analyzed, as can viscous liquids such as heavy oils and greases. Solid samples of varying sizes and shapes may also be examined and with the addition of microscopic IR (microspectroscopy) capabilities, minute materials such as single fibers and threads may be analyzed. With the addition of appropriate software, microspectroscopy can be used for automated discrete point or compositional surface area mapping, with the latter providing a means to record changes in the chemical composition of a material surface over a defined area. Due to the ability to characterize gaseous samples, IR spectroscopy can also be coupled with thermal processes such as thermogravimetric (TG) analyses to provide both thermal and chemical data in a single run. In this configuration, solids (or liquids) heated in a TG analyzer undergo decomposition, with the evolving gases directed into the IR spectrometer. Thus, information is provided on the thermal properties of a material and the order in which its chemical constituents are broken down during incremental heating. Specific examples of these varied applications will be cited, with data interpretation and method limitations further discussed.

Huff, Timothy L.↗

Planetary Protection Bioburden Analysis Program

This program is a Microsoft Access program that performed statistical analysis of the colony counts from assays performed on the Mars Science Laboratory (MSL) spacecraft to determine the bioburden density, 3-sigma biodensity, and the total bioburdens required for the MSL prelaunch reports. It also contains numerous tools that report the data in various ways to simplify the reports required. The program performs all the calculations directly in the MS Access program. Prior to this development, the data was exported to large Excel files that had to be cut and pasted to provide the desired results. The program contains a main menu and a number of submenus. Analyses can be performed by using either all the assays, or only the accountable assays that will be used in the final analysis. There are three options on the first menu: either calculate using (1) the old MER (Mars Exploration Rover) statistics, (2) the MSL statistics for all the assays, or This software implements penetration limit equations for common micrometeoroid and orbital debris (MMOD) shield configurations, windows, and thermal protection systems. Allowable MMOD risk is formulated in terms of the probability of penetration (PNP) of the spacecraft pressure hull. For calculating the risk, spacecraft geometry models, mission profiles, debris environment models, and penetration limit equations for installed shielding configurations are required. Risk assessment software such as NASA's BUMPERII is used to calculate mission PNP; however, they are unsuitable for use in shield design and preliminary analysis studies. The software defines a single equation for the design and performance evaluation of common MMOD shielding configurations, windows, and thermal protection systems, along with a description of their validity range and guidelines for their application. Recommendations are based on preliminary reviews of fundamental assumptions, and accuracy in predicting experimental impact test results. The software is programmed in Visual Basic for Applications for installation as a simple add-in for Microsoft Excel. The user is directed to a graphical user interface (GUI) that requires user inputs and provides solutions directly in Microsoft Excel workbooks. This work was done by Shannon Ryan of the USRA Lunar and Planetary Institute for Johnson Space Center. Further information is contained in a TSP (see page 1). MSC- 24582-1 Micrometeoroid and Orbital Debris (MMOD) Shield Ballistic Limit Analysis Program Lyndon B. Johnson Space Center, Houston, Texas Commercially, because it is so generic, Enigma can be used for almost any project that requires engineering visualization, model building, or animation. Models in Enigma can be exported to many other formats for use in other applications as well. Educationally, Enigma is being used to allow university students to visualize robotic algorithms in a simulation mode before using them with actual hardware. This work was done by David Shores and Sharon P. Goza of Johnson Space Center; Cheyenne McKeegan, Rick Easley, Janet Way, and Shonn Everett of MEI Technologies; Mark Manning of PTI; and Mark Guerra, Ray Kraesig, and William Leu of Tietronix Software, Inc. For further information, contact the JSC Innovation Partnerships Office at (281) 483-3809. MSC-24211-1 Spitzer Telemetry Processing System NASA's Jet Propulsion Laboratory, Pasadena, California The Spitzer Telemetry Processing System (SirtfTlmProc) was designed to address objectives of JPL's Multi-mission Image Processing Lab (MIPL) in processing spacecraft telemetry and distributing the resulting data to the science community. To minimize costs and maximize operability, the software design focused on automated error recovery, performance, and information management. The system processes telemetry from the Spitzer spacecraft and delivers Level 0 products to the Spitzer Science Center. SirtfTlmProc is a unique system with automated error notification and recovery, with a real-time continuous service that can go quiescent after periods of inactivity. The software can process 2 GB of telemetry and deliver Level 0 science products to the end user in four hours. It provides analysis tools so the operator can manage the system and troubleshoot problems. It automates telemetry processing in order to reduce staffing costs. This work was done by Alice Stanboli, Elmain M. Martinez, and James M. McAuley of Caltech for NASA's Jet Propulsion Laboratory. For more information, contact iaoffice @jpl.nasa.gov. This software is available for commercial licensing. Please contact Dan Broderick at Daniel.F. Broderick@jpl.nasa.gov. Refer to NPO-47803. NASA Tech Briefs, September 2013 29 This rapid response computer program predicts Orbiter Wing Leading Edge (WLE) damage caused by ice or foam impact during a Space Shuttle launch (Program "IMPACT2"). The program was developed after the Columbia accident in order to assess quickly WLE damage due to ice, foam, or metal impact (if any) during a Shuttle launch. IMPACT2 simulates an impact event in a few minutes for foam impactors, and in seconds for ice and metal impactors. The damage criterion is derived from results obtained from one sophisticated commercial program, which requires hours to carry out simulations of the same impact events. The program was designed to run much faster than the commercial program with prediction of projectile threshold velocities within 10 to 15% of commercial-program values. The mathematical model involves coupling of Orbiter wing normal modes of vibration to nonlinear or linear springmass models. IMPACT2 solves nonlinear or linear impact problems using classical normal modes of vibration of a target, and nonlinear/ linear time-domain equations for the projectile. Impact loads and stresses developed in the target are computed as functions of time. This model is novel because of its speed of execution. A typical model of foam, or other projectile characterized by material nonlinearities, impacting an RCC panel is executed in minutes instead of hours needed by the commercial programs. Target damage due to impact can be assessed quickly, provided that target vibration modes and allowable stress are known. This work was done by Robert Clark, Jr., Paul Cotter, and Constantine Michalopoulos of The Boeing Company for Johnson Space Center. For further information, contact the JSC Innovation Partnerships Office at (281) 483-3809. MSC-24988-1 Wing Leading Edge RCC Rapid Response Damage Prediction Tool (IMPACT2) Lyndon B. Johnson Space Center, Houston, Texas (3) the MSL statistics for only the accountable assays. Other options on the main menu include a data editing form and utility programs that produce various reports requested by the microbiologists and the project, and tools to generate the groupings for the final analyses. The analyses can be carried out in three ways: Each assay can be treated separately, the assays can be collectively treated for the whole zone as a group, or the assays can be collected in groups designated by the JPL Planetary Protection Manager. The latter approach was used to generate the final report because assays on the same equipment or similar equipment can be assumed to have been exposed to the same environment and cleaning. Thus, the statistics are improved by having a larger population, thereby reducing the standard deviation by the square root of N. For each method mentioned above, three reports are available. The first is a detailed report including all the data. This version was very useful in verifying the calculations. The second is a brief report that is similar to the full detailed report, but does not print out the data. The third is a grand total and summary report in which each assay requires only one line. For the first and second reports, most of the calculations are performed in the report section itself. For the third, all the calculations are performed directly in the query bound to the report. All the numerical results were verified by comparing them with Excel templates, then exporting the data from the Planetary Protection Analysis program to Excel.

Beaudet, Robert A.↗

What the Long-Term Sunspot Record Tells Us About Space Climate

Direct observations of sunspots span the nearly 400 years since the time of Galileo. Dedicated observing programs at several observatories over the last 150 years have provided detailed information not only on the number of sunspots but on their sizes and positions as well. The data acquired by those original observers, and by those who have more recently brought those observations to light, provide important clues about the nature of the solar cycle and its contribution to space climate. The period of the cycle, the equator-ward drift of the active latitudes, the asymmetry between the rise to maximum and the fill to minimum, shifting asymmetries between northern and southern hemisphere activity, the tilt of active regions, and the increasing amplitude of the cycles since the Maunder Minimum are all well established. Other, less well established characteristics such as multi-cycle and short-term periodicities, often depend upon the method of data analysis. The strong correlation between sunspot statistics and other measures of solar activity, coupled with the length of the sunspot record, make these observations extremely valuable for characterizing and understanding space climate.

Hathaway, D. H.↗

Acoustic Testing of Flight Hardware Using Loudspeakers: How Much do We Know About This Method of Testing?

Loudspeakers have been used for acoustic qualification of spacecrafts, reflectors, solar panels, and other acoustically responsive structures for more than a decade. Even though a lot of hardware has been acoustic tested using this method, the nature of the acoustic field generated by controlling an ensemble of speakers with and without the hardware in the test volume has not been thoroughly investigated. Limited measurements from some of the recent speaker tests used to qualify flight hardware have indicated significant spatial variation of the acoustic field within the test volume. Also structural responses have been reported to differ when similar tests were performed using reverberant chambers. Unlike the reverberant chamber acoustic test, for which the acoustic field in most chambers is known to be diffuse except below several tens of Hz where acoustic standing waves and large spatial variations exist, the characteristics of the acoustic field within the speaker test volume has not been quantified. It has only been recently that a detailed acoustic field characterization of speaker testing has been made at Jet Propulsion Laboratory (JPL) with involvement of various organizations. To address the impact of non-uniform acoustic field on structures, a series of acoustic tests were performed using a flat panel and a 3-ft cylinder exposed to the field controlled by speakers and repeated in a reverberant chamber. The analysis of the data from this exercise reveals that there are significant differences both in the acoustic field and in the structural responses. In this paper the differences between the two methods are reviewed in some detail and the over- or under-testing of articles that could pose un-anticipated structural and flight qualification issues are discussed. A framework for discussing the validity of the speaker acoustic testing method with the current control system and a path forward for improving it will be provided.

Acoustics↗

Modeling of the Coupled Magnetospheric and Neutral Wind Dynamos

Over the past four years of funding, SRI, in collaboration with the University of Texas at Dallas, has been involved in assessing the influence of thermospheric neutral winds on the electric field and current systems at high latitudes. The initial direction of the project was to perform a set of numerical experiments concerning the contribution of the magnetospheric and neutral wind dynamo processes, under specific boundary conditions, to the polarization electric field and/or the field-aligned current distribution at high latitudes. To facilitate these numerical experiments we developed a numerical scheme that relied on using output from the NCAR Thermosphere-Ionosphere General Circulation Model (NCAR-TIGCM), expanding them in the form of spherical harmonics and solving the dynamo equations spectrally. Once initial calculations were completed, it was recognized that the neutral wind contribution could be significant but its actual contribution to the electric field or currents depended strongly on the generator properties of the magnetosphere. Solutions to this problem are not unique because of the unknown characteristics of the magnetospheric generator, therefore the focus was on two limiting cases. One limiting case was to consider the magnetosphere as a voltage generator delivering a fixed voltage to the high-latitude ionosphere and allowing for the neutral wind dynamo to contribute only to the current system. The second limiting case was to consider the magnetosphere as a current generator and allowing for the neutral wind dynamo to contribute only to the generation of polarization electric fields. This work was completed and presented at the l994 Fall AGU meeting. The direction of the project then shifted to applying the Poynting flux concept to the high-latitude ionosphere. This concept was more attractive as it evaluated the influence of neutral winds on the high-latitude electrodynamics without actually having to determine the generator characteristics of the magnetosphere. The influence of the neutral wind was then determined not by estimating how much electric potential or current density it provides, but by determining the contribution of the neutral wind to the net electromagnetic energy transferred between the ionosphere and magnetosphere. The estimate of the net electromagnetic energy transfer and the role of the neutral winds proves to be a more fundamental quantity in studies of magnetosphere- ionosphere coupling also showed that by using electric and magnetic field measurements from the HILAT satellite, the Poynting flux could be a measurable quantity from polar-orbiting, low- altitude spacecraft. Through collaboration with Dr. Heelis and others at UTD and their expertise of the electric field measurements on the DE-B satellite, an extensive analysis was planned to determine the Poynting flux from the DE-B measurements in combination with a modeling effort to help interpret the observations taking into account the coupled magnetosphere-ionosphere.

Thayer, Jeffrey P.↗

Control/structure interaction during Space Station Freedom-Orbiter berthing

The berthing maneuver is essential for the construction and assembly of Space Station Freedom (SSF) and has a direct effect on the SSF assembly build up and SSF/Orbiter operations. The effects of flexible body dynamics coupled with the available control system may impose new requirements on the maneuver. The problem is further complicated by the effect of the SSF control system on the Shuttle Remote Manipulator System (SRMS). These effects will play a major role in the development of operational requirements which need to be identified and validated in order to assure total safety and maneuver execution during SSF construction. This paper presents the results of ongoing studies to investigate the Control/Structure Interaction (CSI) during the berthing operations. The problem is formulated in terms of multi-flex body equations of motion for SSF and the SRMS and on-orbit flight control systems for the SRMS and the SSF, which includes the Control Moment Gyro (CMG) and Reaction Control System (RCS) Attitude Control Systems (ACS). The SSF control system designs are based on the Preliminary Design Review (PDR) version of the Honeywell design. The simulation tool used for the analysis is briefly described and the CSI results are presented for given berthing scenarios.

Hua, Tuyen↗

A high frequency analysis of electromagnetic plane wave scattering by perfectly-conducting semi-infinite parallel plate and rectangular waveguides with absorber coated inner walls

An approximate but sufficiently accurate high frequency solution which combines the uniform geometrical theory of diffraction (UTD) and the aperture integration (AI) method is developed for analyzing the problem of electromagnetic (EM) plane wave scattering by an open-ended, perfectly-conducting, semi-infinite hollow rectangular waveguide (or duct) with a thin, uniform layer of lossy or absorbing material on its inner wall, and with a planar termination inside. In addition, a high frequency solution for the EM scattering by a two dimensional (2-D), semi-infinite parallel plate waveguide with a absorber coating on the inner walls is also developed as a first step before analyzing the open-ended semi-infinite three dimensional (3-D) rectangular waveguide geometry. The total field scattered by the semi-infinite waveguide consists firstly of the fields scattered from the edges of the aperture at the open-end, and secondly of the fields which are coupled into the waveguide from the open-end and then reflected back from the interior termination to radiate out of the open-end. The first contribution to the scattered field can be found directly via the UTD ray method. The second contribution is found via the AI method which employs rays to describe the fields in the aperture that arrive there after reflecting from the interior termination. It is assumed that the direction of the incident plane wave and the direction of observation lie well inside the forward half space tht exists outside the half space containing the semi-infinite waveguide geometry. Also, the medium exterior to the waveguide is assumed to be free space.

Noh, H. M.↗

Mars Sample Return Using Commercial Capabilities: Propulsive Entry, Descent, and Landing of a Capsule Form Vehicle

This paper describes a critical portion of the work that has been done at NASA, Ames Research Center regarding the use of the commercially developed Dragon capsule as a delivery vehicle for the elements of a high priority Mars Sample Return mission. The objective of the investigation was to determine entry and landed mass capabilities that cover anticipated mission conditions. The "Red Dragon", Mars configuration, uses supersonic retro-propulsion, with no required parachute system, to perform Entry, Descent, and Landing (EDL) maneuvers. The propulsive system proposed for use is the same system that will perform an abort, if necessary, for a human rated version of the Dragon capsule. Standard trajectory analysis tools are applied to publically available information about Dragon and other legacy capsule forms in order to perform the investigation. Trajectory simulation parameters include entry velocity, flight path angle, lift to drag Ratio (L/D), landing site elevation, atmosphere density, and total entry mass, in addition engineering assumptions for the performance of the propulsion system are stated. Mass estimates for major elements of the overall proposed architecture are coupled to this EDL analysis to close the overall architecture. Three synodic launch opportunities, beginning with the 2022 opportunity, define the arrival conditions. Results state the relations between the analysis parameters as well as sensitivities to those parameters. The EDL performance envelope includes landing altitudes between 0 and -4 km referenced to the Mars Orbiter Laser Altimeter datum as well as minimum and maximum atmosphere density. Total entry masses between 7 and 10 mt are considered with architecture closure occurring between 9.0 and 10 mt. Propellant mass fractions for each major phase of the EDL - Entry, Terminal Descent, and Hazard Avoidance - have been derived. An assessment of the effect of the entry conditions on the Thermal Protection System (TPS) currently in use for Dragon missions shows no significant stressors. A useful payload mass of 2.0 mt is provided and includes mass and grow allowance for a Mars Ascent Vehicle (MAV), Earth Return Vehicle (ERV), and mission unique equipment. The useful payload supports an architecture that receives a sample from another surface asset and sends it directly back to Earth for recovery in a high Earth orbit. The work shows that emerging commercial capabilities as well as previously studied EDL methodologies can be used to efficiently support an important planetary science objective. The work also has applications for human exploration missions that will also use propulsive EDL techniques

EDL↗

Parametric Modeling Investigation of a Radially-Staged Low-Emission Aviation Combustor

Aviation gas-turbine combustion demands high efficiency, wide operability and minimal trace gas emissions. Performance critical design parameters include injector geometry, combustor layout, fuel-air mixing and engine cycle conditions. The present investigation explores these factors and their impact on a radially staged low-emission aviation combustor sized for a next-generation 24,000-lbf-thrust engine. By coupling multi-fidelity computational tools, a design exploration was performed using a parameterized annular combustor sector at projected 100% takeoff power conditions. Design objectives included nitrogen oxide emission indices and overall combustor pressure loss. From the design space, an optimal configuration was selected and simulated at 7.1, 30 and 85% part-power operation, corresponding to landing-takeoff cycle idle, approach and climb segments. All results were obtained by solution of the steady-state Reynolds-averaged Navier-Stokes equations. Species concentrations were solved directly using a reduced 19-step reaction mechanism for Jet-A. Turbulence closure was obtained using a nonlinear K-epsilon model. This research demonstrates revolutionary combustor design exploration enabled by multi-fidelity physics-based simulation.

design analysis↗

Satellite Attitude Control Using a Combination of Inertia Wheels and a Bar Magnet

The use of inertia wheels to control the attitude of a satellite has currently aroused much interest. The stability of such a system has been studied in this investigation. A single-degree-of-freedom analysis indicates that a response with suitable dynamic characteristics and precise control can be achieved by commanding the angular velocity of the inertia wheel with an error signal that is the sum of the attitude error, the attitude rate, and the integral of the attitude error. A digital computer was used to study the three-degree-of-freedom response to step displacements, and the results indicate that the cross-coupling effects of inertia coupling and precession coupling had no effect on system stability. A study was also made of the use of a bar magnet to supplement the inertia wheels by providing a means of removing any momentum introduced into the system by disturbances such as aerodynamic torques. A study of a case with large aerodynamic torques, with a typical orbit, indicated that the magnet was a suitable device for supplying the essential trimming force. Single-degree-of-freedom bench tests generally verified the dynamic response predicted by the analytical study. the test table to within plus or minus 9 arc-seconds of the reference direction, even though the hardware components that were used in these tests were not specifically designed for the control system.

Adams, James J.↗

Particle kinetic simulation of high altitude hypervelocity flight

In this grant period, the focus has been on enhancement and application of the direct simulation Monte Carlo (DSMC) particle method for computing hypersonic flows of re-entry vehicles. Enhancement efforts dealt with modeling gas-gas interactions for thermal non-equilibrium relaxation processes and gas-surface interactions for prediction of vehicle surface temperatures. Both are important for application to problems of engineering interest. The code was employed in a parametric study to improve future applications, and in simulations of aeropass maneuvers in support of the Magellan mission. Detailed comparisons between continuum models for internal energy relaxation and DSMC models reveals that several discrepancies exist. These include definitions of relaxation parameters and the methodologies for implementing them in DSMC codes. These issues were clarified and all differences were rectified in a paper (Appendix A) submitted to Physics of Fluids A, featuring several key figures in the DSMC community as co-authors and B. Haas as first author. This material will be presented at the Fluid Dynamics meeting of the American Physical Society on November 21, 1993. The aerodynamics of space vehicles in highly rarefied flows are very sensitive to the vehicle surface temperatures. Rather than require prescribed temperature estimates for spacecraft as is typically done in DSMC methods, a new technique was developed which couples the dynamic surface heat transfer characteristics into the DSMC flow simulation code to compute surface temperatures directly. This model, when applied to thin planar bodies such as solar panels, was described in AIAA Paper No. 93-2765 (Appendix B) and was presented at the Thermophysics Conference in July 1993. The paper has been submitted to the Journal of Thermophysics and Heat Transfer. Application of the DSMC method to problems of practical interest requires a trade off between solution accuracy and computational expense and limitations. A parametric study was performed and reported in AIAA Paper No. 93-2806 (Appendix C) which assessed the accuracy penalties associated with simulations of varying grid resolution and flow domain size. The paper was also presented at the Thermophysics Conference and will be submitted to the journal shortly. Finally, the DSMC code was employed to assess the pitch, yaw, and roll aerodynamics of the Magellan spacecraft during entry into the Venus atmosphere at off-design attitudes. This work was in support of the Magellan aerobraking maneuver of May 25-Aug. 3, 1993. Furthermore, analysis of the roll characteristics of the configuration with canted solar panels was performed in support of the proposed 'Windmill' experiment. Results were reported in AIAA Paper No. 93-3676 (Appendix D) presented at the Atmospheric Flight Mechanics Conference in August 1993, and were submitted to Journal of Spacecraft and Rockets.

Haas, Brian L.↗

WARP3D-Release 10.8: Dynamic Nonlinear Analysis of Solids using a Preconditioned Conjugate Gradient Software Architecture

This report describes theoretical background material and commands necessary to use the WARP3D finite element code. WARP3D is under continuing development as a research code for the solution of very large-scale, 3-D solid models subjected to static and dynamic loads. Specific features in the code oriented toward the investigation of ductile fracture in metals include a robust finite strain formulation, a general J-integral computation facility (with inertia, face loading), an element extinction facility to model crack growth, nonlinear material models including viscoplastic effects, and the Gurson-Tver-gaard dilatant plasticity model for void growth. The nonlinear, dynamic equilibrium equations are solved using an incremental-iterative, implicit formulation with full Newton iterations to eliminate residual nodal forces. The history integration of the nonlinear equations of motion is accomplished with Newmarks Beta method. A central feature of WARP3D involves the use of a linear-preconditioned conjugate gradient (LPCG) solver implemented in an element-by-element format to replace a conventional direct linear equation solver. This software architecture dramatically reduces both the memory requirements and CPU time for very large, nonlinear solid models since formation of the assembled (dynamic) stiffness matrix is avoided. Analyses thus exhibit the numerical stability for large time (load) steps provided by the implicit formulation coupled with the low memory requirements characteristic of an explicit code. In addition to the much lower memory requirements of the LPCG solver, the CPU time required for solution of the linear equations during each Newton iteration is generally one-half or less of the CPU time required for a traditional direct solver. All other computational aspects of the code (element stiffnesses, element strains, stress updating, element internal forces) are implemented in the element-by- element, blocked architecture. This greatly improves vectorization of the code on uni-processor hardware and enables straightforward parallel-vector processing of element blocks on multi-processor hardware.

Koppenhoefer, Kyle C.↗

Development of an Aerodynamic Analysis Method and Database for the SLS Service Module Panel Jettison Event Utilizing Inviscid CFD and MATLAB

This paper describes the development, testing, and utilization of an aerodynamic force and moment database for the Space Launch System (SLS) Service Module (SM) panel jettison event. The database is a combination of inviscid Computational Fluid Dynamic (CFD) data and MATLAB code written to query the data at input values of vehicle/SM panel parameters and return the aerodynamic force and moment coefficients of the panels as they are jettisoned from the vehicle. The database encompasses over 5000 CFD simulations with the panels either in the initial stages of separation where they are hinged to the vehicle, in close proximity to the vehicle, or far enough from the vehicle that body interference effects are neglected. A series of viscous CFD check cases were performed to assess the accuracy of the Euler solutions for this class of problem and good agreement was obtained. The ultimate goal of the panel jettison database was to create a tool that could be coupled with any 6-Degree-Of-Freedom (DOF) dynamics model to rapidly predict SM panel separation from the SLS vehicle in a quasi-unsteady manner. Results are presented for panel jettison simulations that utilize the database at various SLS flight conditions. These results compare favorably to an approach that directly couples a 6-DOF model with the Cart3D Euler flow solver and obtains solutions for the panels at exact locations. This paper demonstrates a method of using inviscid CFD simulations coupled with a 6-DOF model that provides adequate fidelity to capture the physics of this complex multiple moving-body panel separation event.

Applebaum, Michael P.↗

Origin of Ozone NO(x) in the Tropical Troposphere: A Photochemical Analysis of Aircraft Observations Over the South Atlantic Basin

The photochemistry of the troposphere over the South Atlantic basin is examined by modeling of aircraft observations up to 12-km altitude taken during the TRACE A expedition in September-October 1992. A close balance is found in the 0 to 12-km column between photochemical production and loss Of O3, with net production at high altitudes compensating for weak net loss at low altitudes. This balance implies that O3 concentrations in the 0-12 km column can be explained solely by in situ photochemistry; influx from the stratosphere is negligible. Simulation of H2O2, CH3OOH, and CH2O concentrations measured aboard the aircraft lends confidence in the computations of O3 production and loss rates, although there appears to be a major gap in current understanding of CH2O chemistry in the marine boundary layer. The primary sources of NO(x) over the South Atlantic Basin appear to be continental (biomass burning, lightning, soils). There is evidence that NO(x) throughout the 0 to 12-km column is recycled from its oxidation products rather than directly transported from its primary sources. There is also evidence for rapid conversion of HNO3 to NO(x) in the upper troposphere by a mechanism not included in current models. A general representation of the O3 budget in the tropical troposphere is proposed that couples the large scale Walker circulation and in situ photochemistry. Deep convection in the rising branches of the Walker circulation injects NO(x) from combustion, soils, and lightning to the upper troposphere, leading to O3 production; eventually, the air subsides and net O3 loss takes place in the lower troposphere, closing the O3 cycle. This scheme implies a great sensitivity of the oxidizing power of the atmosphere to NO(x) emissions in the tropics.

Jacob, D. J.↗

Aerodynamic Optimization of Rocket Control Surface Geometry Using Cartesian Methods and CAD Geometry

Aerodynamic design is an iterative process involving geometry manipulation and complex computational analysis subject to physical constraints and aerodynamic objectives. A design cycle consists of first establishing the performance of a baseline design, which is usually created with low-fidelity engineering tools, and then progressively optimizing the design to maximize its performance. Optimization techniques have evolved from relying exclusively on designer intuition and insight in traditional trial and error methods, to sophisticated local and global search methods. Recent attempts at automating the search through a large design space with formal optimization methods include both database driven and direct evaluation schemes. Databases are being used in conjunction with surrogate and neural network models as a basis on which to run optimization algorithms. Optimization algorithms are also being driven by the direct evaluation of objectives and constraints using high-fidelity simulations. Surrogate methods use data points obtained from simulations, and possibly gradients evaluated at the data points, to create mathematical approximations of a database. Neural network models work in a similar fashion, using a number of high-fidelity database calculations as training iterations to create a database model. Optimal designs are obtained by coupling an optimization algorithm to the database model. Evaluation of the current best design then gives either a new local optima and/or increases the fidelity of the approximation model for the next iteration. Surrogate methods have also been developed that iterate on the selection of data points to decrease the uncertainty of the approximation model prior to searching for an optimal design. The database approximation models for each of these cases, however, become computationally expensive with increase in dimensionality. Thus the method of using optimization algorithms to search a database model becomes problematic as the number of design variables is increased.

Nelson, Andrea↗