Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “conditionally specified models”

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 523 records · Page 29
All sourcesNASA NTRSDOE OSTI
NASA NTRS1998

Liquid Motion in a Rotating Tank Experiment (LME)

The Liquid Motion Experiment (LME), which flew on STS 84 in May 1997, was an investigation of liquid motions in spinning, nutating tanks. LME was designed to quantify the effects of such liquid motions on the stability of spinning spacecraft, which are known to be adversely affected by the energy dissipated by the liquid motions. The LME hardware was essentially a spin table which could be forced to nutate at specified frequencies at a constant cone angle, independently of the spin rate. Cylindrical and spherical test tanks, partially filled with liquids of different viscosities, were located at the periphery of the spin table to simulate a spacecraft with off-axis propellant tanks; one set of tanks contained generic propellant management devices (PMDs). The primary quantitative data from the flight tests were the liquid-induced torques exerted on the tanks about radial and tangential axes through the center of the tank. Visual recordings of the liquid oscillations also provided qualitative information. The flight program incorporated two types of tests: sine sweep tests, in which the spin rate was held constant and the nutation frequency varied over a wide range; and sine dwell test, in which both the spin rate and the nutation frequency were held constant. The sine sweep tests were meant to investigate all the prominent liquid resonant oscillations and the damping of the resonances, and the sine dwell tests were meant to quantify the viscous energy dissipation rate of the liquid oscillations for steady state conditions. The LME flight data were compared to analytical results obtained from two companion IR&D programs at Southwest Research Institute. The comparisons indicated that the models predicted the observed liquid resonances, damping, and energy dissipation rates for many test conditions but not for all. It was concluded that improved models and CFD simulations are needed to resolve the differences. This work is ongoing under a current IR&D program.

Deffenbaugh, D. M.↗
NASA NTRSDate not supplied

Probabalistic Risk Analysis and Thermal Margin Process for an Inflatable Aeroshell

Uncertainties always exist in atmospheric entry aeroheating environments and the thermal response of thermal protection system (TPS) material. These uncertainties are mitigated in the design by ap-plying margin and factors of safety to the TPS. Entry vehicle TPS is often conservatively over-sized for the heat loads that are experienced along the entry trajectory by stacking worst-case scenarios together. Additionally, the current TPS design and margin process used by NASA offers very little insight into the risk of over-temperature during flight and the reliability of the heat shield performance [1,3]. A probabilistic margin process can be used to calculate the amount of TPS margin necessary to survive a given entry heat load at a specified level of risk [2,3,4]. The vehicle’s initial entry state (entry velocity, flight path angle, and entry mass) determines the expected atmospheric entry environmental conditions and resulting heat load that the entry vehicle will experience. If there is flexibility in the entry state, then this process can be used to select an appropriate combination of entry state parameters and TPS size to target a desired reentry reliability. This probabilistic margin process allows engineers to make informed aeroshell design, entry-trajectory design, and TPS performance risk trades while preventing excessive TPS margin from being applied. The probabilistic TPS margin process has been performed to determine TPS thickness and entry heating constraints given an acceptable risk level for the Low Earth Orbit Flight Experiment of an Inflatable Decelerator (LOFTID) flight project. The process is used in a manner to size the entry heat load for a given flexible TPS (FTPS) thickness so that it meets project reliability standards while allowing the FTPS and the underlying inflatable structure (IS) to be pushed to adequately high temperatures. Since the LOFTID project is an experimental flight demonstration, it is de-sired to drive the FTPS and IS to temperatures that cover a large range of their thermal response models’ applicability. This will allow the thermal response models to be better improved and validated post-flight using LOFTID’s extensive instrumentation embedded within the aeroshell. The presentation demonstrates how uncertainty analysis is carried out using an end-to-end Monte Carlo process where three separate Monte Carlo simulations are run in sequence. The first Monte Carlo simulation operates on the entry trajectory model to generate trajectory parameter dispersions that are fed into the second Monte Carlo simulation. The second Monte Carlo simulation operates on the aerothermodynamics model to generate aeroheating parameter dispersions that are fed into the third Monte Carlo simulation. The third Monte Carlo simulation operates on the FTPS material thermal response model to generate the final FTPS/IS thermal response dispersions. The end-to-end Monte Carlo simulation propagates the uncertainties of each model into the next to quantify the resulting uncertainty of the FTPS/IS thermal response. The fractional contributions of the uncertain parameters in the trajectory, aerothermal, and thermal response models to the variance in the FTPS/IS thermal response is determined as a byproduct of the Monte Carlo analysis. The structural uncertainty of the FTPS thermal response model is evaluated by flight relevant ground testing and model error analysis using test measurements. This probabilistic TPS margin process had never been applied to an entry vehicle and it is one of the LOFTID project’s goals to demonstrate its merits.

Steven A. Tobin↗
NASA NTRS2002

Program Instrumentation and Trace Analysis

Several attempts have been made recently to apply techniques such as model checking and theorem proving to the analysis of programs. This shall be seen as a current trend to analyze real software systems instead of just their designs. This includes our own effort to develop a model checker for Java, the Java PathFinder 1, one of the very first of its kind in 1998. However, model checking cannot handle very large programs without some kind of abstraction of the program. This paper describes a complementary scalable technique to handle such large programs. Our interest is turned on the observation part of the equation: How much information can be extracted about a program from observing a single execution trace? It is our intention to develop a technology that can be applied automatically and to large full-size applications, with minimal modification to the code. We present a tool, Java PathExplorer (JPaX), for exploring execution traces of Java programs. The tool prioritizes scalability for completeness, and is directed towards detecting errors in programs, not to prove correctness. One core element in JPaX is an instrumentation package that allows to instrument Java byte code files to log various events when executed. The instrumentation is driven by a user provided script that specifies what information to log. Examples of instructions that such a script can contain are: 'report name and arguments of all called methods defined in class C, together with a timestamp'; 'report all updates to all variables'; and 'report all acquisitions and releases of locks'. In more complex instructions one can specify that certain expressions should be evaluated and even that certain code should be executed under various conditions. The instrumentation package can hence be seen as implementing Aspect Oriented Programming for Java in the sense that one can add functionality to a Java program without explicitly changing the code of the original program, but one rather writes an aspect and compiles it into the original program using the instrumentation. Another core element of JPaX is an observation package that supports the analysis of the generated event stream. Two kinds of analysis are currently supported. In temporal analysis the execution trace is evaluated against formulae written in temporal logic. We have implemented a temporal logic evaluator on finite traces using the Maude rewriting system from SRI International, USA. Temporal logic is defined in Maude by giving its syntax as a signature and its semantics as rewrite equations. The resulting semantics is extremely efficient and can handle event streams of hundreds of millions events in few minutes. Furthermore, the implementation is very succinct. The second form of even stream analysis supported is error pattern analysis where an execution trace is analyzed using various error detection algorithms that can identify error-prone programming practices that may potentially lead to errors in some different executions. Two such algorithms focusing on concurrency errors have been implemented in JPaX, one for deadlocks and the other for data races. It is important to note, that a deadlock or data race potential does not need to occur in order for its potential to be detected with these algorithms. This is what makes them very scalable in practice. The data race algorithm implemented is the Eraser algorithm from Compaq, however adopted to Java. The tool is currently being applied to a code base for controlling a spacecraft by the developers of that software in order to evaluate its applicability.

Havelund, Klaus↗
NASA NTRS2012

Creation of an Upper Stage Trajectory Capability Boundary to Enable Booster System Trade Space Exploration

The problem of trajectory optimization is important in all space missions. The solution of this problem enables one to specify the optimum thrust steering program which should be followed to achieve a specified mission objective, simultaneously satisfying the constraints.1 It is well known that whether or not the ascent trajectory is optimal can have a significant impact on propellant usage for a given payload, or on payload weight for the same gross vehicle weight.2 Consequently, ascent guidance commands are usually optimized in some fashion. Multi-stage vehicles add complexity to this analysis process as changes in vehicle properties in one stage propagate to the other stages through gear ratios and changes in the optimal trajectory. These effects can cause an increase in analysis time as more variables are added and convergence of the optimizer to system closure requires more analysis iterations. In this paper, an approach to simplifying this multi-stage problem through the creation of an upper stage capability boundary is presented. This work was completed as part of a larger study focused on trade space exploration for the advanced booster system that will eventually form a part of NASA s new Space Launch System.3 The approach developed leverages Design of Experiments and Surrogate Modeling4 techniques to create a predictive model of the SLS upper stage performance. The design of the SLS core stages is considered fixed for the purposes of this study, which results in trajectory parameters such as staging conditions being the only variables relevant to the upper stage. Through the creation of a surrogate model, which takes staging conditions as inputs and predicts the payload mass delivered by the SLS upper stage to a reference orbit as the response, it is possible to identify a "surface" of staging conditions which all satisfy the SLS requirement of placing 130 metric tons into low-Earth orbit (LEO).3 This identified surface represents the 130 metric ton capability boundary for the upper stage, such that if the combined first stage and boosters can achieve any one staging point on that surface, then the design is identified as feasible. With the surrogate model created, design and analysis of advanced booster concepts is streamlined, as optimization of the upper stage trajectory is no longer required in every design loop.

Walsh, Ptrick↗
NASA NTRS2015

Experimental and Finite Element Modeling of Near-Threshold Fatigue Crack Growth for the K-Decreasing Test Method

The experimental methods to determine near-threshold fatigue crack growth rate data are prescribed in ASTM standard E647. To produce near-threshold data at a constant stress ratio (R), the applied stress-intensity factor (K) is decreased as the crack grows based on a specified K-gradient. Consequently, as the fatigue crack growth rate threshold is approached and the crack tip opening displacement decreases, remote crack wake contact may occur due to the plastically deformed crack wake surfaces and shield the growing crack tip resulting in a reduced crack tip driving force and non-representative crack growth rate data. If such data are used to life a component, the evaluation could yield highly non-conservative predictions. Although this anomalous behavior has been shown to be affected by K-gradient, starting K level, residual stresses, environmental assisted cracking, specimen geometry, and material type, the specifications within the standard to avoid this effect are limited to a maximum fatigue crack growth rate and a suggestion for the K-gradient value. This paper provides parallel experimental and computational simulations for the K-decreasing method for two materials (an aluminum alloy, AA 2024-T3 and a titanium alloy, Ti 6-2-2-2-2) to aid in establishing clear understanding of appropriate testing requirements. These simulations investigate the effect of K-gradient, the maximum value of stress-intensity factor applied, and material type. A material independent term is developed to guide in the selection of appropriate test conditions for most engineering alloys. With the use of such a term, near-threshold fatigue crack growth rate tests can be performed at accelerated rates, near-threshold data can be acquired in days instead of weeks without having to establish testing criteria through trial and error, and these data can be acquired for most engineering materials, even those that are produced in relatively small product forms.

Smith, Stephen W.↗
NASA NTRS1996

A Sundial-Atlas Precursor to the TIMED Mission: A Quick-Response Global Investigation into Coupled Lower Thermospheric, Ionospheric, and Mesospheric Physics

The SUNDIAL-ATLAS effort was a global-scale investigation which responded to the science priorities of the ITM Panel, the Integrated SPD Strategy Implementation Plan as a whole, and the need for potential cost-saving design criteria for the TIMED mission. The investigation focused on coupling processes in the ionospheric-thermospheric system, taking advantage of the timelines of the ATLAS-1 mission (March 1992), and the global-scale ground-based measurement and modeling activities of the SUNDIAL program. The collaborative SUNDIAL-ATLAS activity was the first opportunity for global measurements of the chemistry, kinetics, and electrodynamics which couple the E-, Fl-, and F2-regions into a single interactive system. As such, the program represented an important first step in studying global issues; and accordingly, was an important proof of concept experiment relevant to the strategic mission plans for the ITM community and the upcoming intermediate class satellite program called TIMED. To meet its projected goals, TIMED must perform a number of critical measurements and execute a number of correlations that were to be tried and tested for the first time in the SUNDIAL-ATLAS investigation. This was designed to include global correlations of thermospheric and ionospheric composition during quiet and disturbed conditions and the co-registration of global-scale ground-based measurements with along-track satellite diagnostics. The SUNDIAL component of the current investigation addressed this need by acquiring, reducing, and analyzing a multi-sensor database that complemented and extended that which was generated in the ATLAS mission (Atmospheric Laboratory for Applications and Science). The SUNDIAL data defined the state and condition of the global-scale ionosphere in the altitude range from 100 km to the F2-peak. These data specified the peak heights and densities of the E-, Fl-, and F2-regions, along with the global distributions of intermediate, descending, and sequential layers which play a critical role in the dynamo region of the lower ionospheric-thermospheric domain. The data were collected by the SUNDIAL network of more than 50 ground-based stations utilizing ionosondes, radars, photometers, Fabry-Perot interferometers, and total electron content measurements. The data were acquired during a three-week period centered on the eight-day ATLAS-1 mission, which provided image and photometric sensing of the altitude distributions of the major and minor ions and neutrals in the ITM system. This report focuses on the scientific contributions of the SUNDIAL component of the overall investigation. Specific findings are described in seven papers (attached) published in the Journal of Geophysical Research.

Szuszczewicz, E. P.↗
NASA NTRS2015

4.4 Development of a 30-Year Soil Moisture Climatology for Situational Awareness and Public Health Applications

This paper provided a brief background on the work being done at NASA SPoRT and the CDC to create a soil moisture climatology over the CONUS at high spatial resolution, and to provide a valuable source of soil moisture information to the CDC for monitoring conditions that could favor the development of Valley Fever. The soil moisture climatology has multi-faceted applications for both the NOAA/NWS situational awareness in the areas of drought and flooding, and for the Public Health community. SPoRT plans to increase its interaction with the drought monitoring and Public Health communities by enhancing this testbed soil moisture anomaly product. This soil moisture climatology run will also serve as a foundation for upgrading the real-time (currently southeastern CONUS) SPoRT-LIS to a full CONUS domain based on LIS version 7 and incorporating real-time GVF data from the Suomi-NPP Visible Infrared Imaging Radiometer Suite (Vargas et al. 2013) into LIS-Noah. The upgraded SPoRT-LIS run will serve as a testbed proof-of-concept of a higher-resolution NLDAS-2 modeling member. The climatology run will be extended to near real-time using the NLDAS-2 meteorological forcing from 2011 to present. The fixed 1981-2010 climatology shall provide the soil moisture "normals" for the production of real-time soil moisture anomalies. SPoRT also envisions a web-mapping type of service in which an end-user could put in a request for either an historical or real-time soil moisture anomaly graph for a specified county (as exemplified by Figure 2) and/or for local and regional maps of soil moisture proxy percentiles. Finally, SPoRT seeks to assimilate satellite soil moisture data from the current Soil Moisture Ocean Salinity (SMOS; Blankenship et al. 2014) and the recently-launched NASA Soil Moisture Active Passive (SMAP; Entekhabi et al. 2010) missions, using the EnKF capability within LIS. The 9-km combined active radar and passive microwave retrieval product from SMAP (Das et al. 2011) has the potential to provide valuable information about the near-surface soil moisture state for improving land surface modeling output.

land surface modeling↗
NASA NTRS2009

Protocol for Communication Networking for Formation Flying

An application-layer protocol and a network architecture have been proposed for data communications among multiple autonomous spacecraft that are required to fly in a precise formation in order to perform scientific observations. The protocol could also be applied to other autonomous vehicles operating in formation, including robotic aircraft, robotic land vehicles, and robotic underwater vehicles. A group of spacecraft or other vehicles to which the protocol applies could be characterized as a precision-formation- flying (PFF) network, and each vehicle could be characterized as a node in the PFF network. In order to support precise formation flying, it would be necessary to establish a corresponding communication network, through which the vehicles could exchange position and orientation data and formation-control commands. The communication network must enable communication during early phases of a mission, when little positional knowledge is available. Particularly during early mission phases, the distances among vehicles may be so large that communication could be achieved only by relaying across multiple links. The large distances and need for omnidirectional coverage would limit communication links to operation at low bandwidth during these mission phases. Once the vehicles were in formation and distances were shorter, the communication network would be required to provide high-bandwidth, low-jitter service to support tight formation-control loops. The proposed protocol and architecture, intended to satisfy the aforementioned and other requirements, are based on a standard layered-reference-model concept. The proposed application protocol would be used in conjunction with conventional network, data-link, and physical-layer protocols. The proposed protocol includes the ubiquitous Institute of Electrical and Electronics Engineers (IEEE) 802.11 medium access control (MAC) protocol to be used in the datalink layer. In addition to its widespread and proven use in diverse local-area networks, this protocol offers both (1) a random- access mode needed for the early PFF deployment phase and (2) a time-bounded-services mode needed during PFF-maintenance operations. Switching between these two modes could be controlled by upper-layer entities using standard link-management mechanisms. Because the early deployment phase of a PFF mission can be expected to involve multihop relaying to achieve network connectivity (see figure), the proposed protocol includes the open shortest path first (OSPF) network protocol that is commonly used in the Internet. Each spacecraft in a PFF network would be in one of seven distinct states as the mission evolved from initial deployment, through coarse formation, and into precise formation. Reconfiguration of the formation to perform different scientific observations would also cause state changes among the network nodes. The application protocol provides for recognition and tracking of the seven states for each node and for protocol changes under specified conditions to adapt the network and satisfy communication requirements associated with the current PFF mission phase. Except during early deployment, when peer-to-peer random access discovery methods would be used, the application protocol provides for operation in a centralized manner.

Jennings, Esther↗
NASA NTRS2003

Formation of Supermassive Black Holes in Galactic Bulges: A Rotating Collapse Model Consistent with the M(sub BH-sigma) Relation

Motivated by the observed correlation between black hole masses M(sub BH) and the velocity dispersion sigma of host galaxies, we develop a theoretical model of black hole formation in galactic bulges (this paper generalizes an earlier ApJ Letter). The model assumes an initial state specified by a uniform rotation rate OMEGA and a density distribution of the form rho = a(sup 2)(sub eff)per2piGR(sup 2)(so that a(sub eff)is an effective transport speed). The black hole mass is determined when the centrifugal radius of the collapse flow exceeds the capture radius of the central black hole (for Schwarzschild geometry). This model reproduces the observed correlation between the estimated black hole masses and the velocity dispersions of galactic bulges, i.e., M(sub BH) approximately equal to 10(sup 8) solar mass(sigma per 200 kilometers per second)(sup 4) where sigma = the square root of 2a(sub eff). To obtain this normalization, the rotation rate OMEGA approximately equal to 2 x 10(exp -15) rad per second. The model also defines a bulge mass scale M(sub B). If we identify the scale M(sub B) with the bulge mass, the model determines the ratio mu(sub B) of black hole mass to the host mass: mu(sub B) approximately equal to 0.0024(sigma per 200 kilometer per second), again in reasonable agreement with observed values. In this scenario, supermassive black holes form quickly (in approximately 10(exp 5) yr) and are born rapidly rotating (with a per M approximately 0.9). This paper also shown how these results depend on the assumed initial conditions; the most important quantity is the initial distribution of specific angular momentum in the precollapse state.

Adams, Fred C.↗
NASA NTRS1996

Conversion-Integration of MSFC Nonlinear Signal Diagnostic Analysis Algorithms for Realtime Execution of MSFC's MPP Prototype System

NASA's advanced propulsion system Small Scale Magnetic Disturbances/Advanced Technology Development (SSME/ATD) has been undergoing extensive flight certification and developmental testing, which involves large numbers of health monitoring measurements. To enhance engine safety and reliability, detailed analysis and evaluation of the measurement signals are mandatory to assess its dynamic characteristics and operational condition. Efficient and reliable signal detection techniques will reduce the risk of catastrophic system failures and expedite the evaluation of both flight and ground test data, and thereby reduce launch turn-around time. During the development of SSME, ASRI participated in the research and development of several advanced non- linear signal diagnostic methods for health monitoring and failure prediction in turbomachinery components. However, due to the intensive computational requirement associated with such advanced analysis tasks, current SSME dynamic data analysis and diagnostic evaluation is performed off-line following flight or ground test with a typical diagnostic turnaround time of one to two days. The objective of MSFC's MPP Prototype System is to eliminate such 'diagnostic lag time' by achieving signal processing and analysis in real-time. Such an on-line diagnostic system can provide sufficient lead time to initiate corrective action and also to enable efficient scheduling of inspection, maintenance and repair activities. The major objective of this project was to convert and implement a number of advanced nonlinear diagnostic DSP algorithms in a format consistent with that required for integration into the Vanderbilt Multigraph Architecture (MGA) Model Based Programming environment. This effort will allow the real-time execution of these algorithms using the MSFC MPP Prototype System. ASRI has completed the software conversion and integration of a sequence of nonlinear signal analysis techniques specified in the SOW for real-time execution on MSFC's MPP Prototype. This report documents and summarizes the results of the contract tasks; provides the complete computer source code; including all FORTRAN/C Utilities; and all other utilities/supporting software libraries that are required for operation.

Jong, Jen-Yi↗
NASA NTRS2009

Conical-Domain Model for Estimating GPS Ionospheric Delays

The conical-domain model is a computational model, now undergoing development, for estimating ionospheric delays of Global Positioning System (GPS) signals. Relative to the standard ionospheric delay model described below, the conical-domain model offers improved accuracy. In the absence of selective availability, the ionosphere is the largest source of error for single-frequency users of GPS. Because ionospheric signal delays contribute to errors in GPS position and time measurements, satellite-based augmentation systems (SBASs) have been designed to estimate these delays and broadcast corrections. Several national and international SBASs are currently in various stages of development to enhance the integrity and accuracy of GPS measurements for airline navigation. In the Wide Area Augmentation System (WAAS) of the United States, slant ionospheric delay errors and confidence bounds are derived from estimates of vertical ionospheric delay modeled on a grid at regularly spaced intervals of latitude and longitude. The estimate of vertical delay at each ionospheric grid point (IGP) is calculated from a planar fit of neighboring slant delay measurements, projected to vertical using a standard, thin-shell model of the ionosphere. Interpolation on the WAAS grid enables estimation of the vertical delay at the ionospheric pierce point (IPP) corresponding to any arbitrary measurement of a user. (The IPP of a given user s measurement is the point where the GPS signal ray path intersects a reference ionospheric height.) The product of the interpolated value and the user s thin-shell obliquity factor provides an estimate of the user s ionospheric slant delay. Two types of error that restrict the accuracy of the thin-shell model are absent in the conical domain model: (1) error due to the implicit assumption that the electron density is independent of the azimuthal angle at the IPP and (2) error arising from the slant-to-vertical conversion. At low latitudes or at mid-latitudes under disturbed conditions, the accuracy of SBAS systems based upon the thin-shell model suffers due to the presence of complex ionospheric structure, high delay values, and large electron density gradients. Interpolation on the vertical delay grid serves as an additional source of delay error. The conical-domain model permits direct computation of the user s slant delay estimate without the intervening use of a vertical delay grid. The key is to restrict each fit of GPS measurements to a spatial domain encompassing signals from only one satellite. The conical domain model is so named because each fit involves a group of GPS receivers that all receive signals from the same GPS satellite (see figure); the receiver and satellite positions define a cone, the satellite position being the vertex. A user within a given cone evaluates the delay to the satellite directly, using (1) the IPP coordinates of the line of sight to the satellite and (2) broadcast fit parameters associated with the cone. The conical-domain model partly resembles the thin-shell model in that both models reduce an inherently four-dimensional problem to two dimensions. However, unlike the thin-shell model, the conical domain model does not involve any potentially erroneous simplifying assumptions about the structure of the ionosphere. In the conical domain model, the initially four-dimensional problem becomes truly two-dimensional in the sense that once a satellite location has been specified, any signal path emanating from a satellite can be identified by only two coordinates; for example, the IPP coordinates. As a consequence, a user s slant-delay estimate converges to the correct value in the limit that the receivers converge to the user s location (or, equivalently, in the limit that the measurement IPPs converge to the user s IPP).

Sparks, Lawrence↗
NASA NTRS2014

Quantification and Mitigation of Long-Term Impacts of Urbanization and Climate Change in the Tropical Coastal City of San Juan, Puerto Rico

Urbanization, along with other cases of land cover and land use changes, has significant climate impacts in tropical regions with the added complexity of occurring within the context of global warming. The individual and combined effects of these two factors on the surface energy balance of a tropical city are investigated by use of an integrated atmospheric modeling approach, taking the San Juan Metropolitan Area (SJMA), Puerto Rico as the test case. To achieve this goal, an ensemble of climate and weather simulations is performed, with the climate scenarios combining urban development and sprawl with regional climate change over the past 50 years, and the short-term simulations designed to test the sensitivity to different urban vegetation configurations as mitigating alternatives. As indicator of change, we use the thermal response number (TRN), which is a measure of the sensible heating to the thermal storage of a surface or region, and the Bowen ratio, which is defined as the ratio of sensible to latent heat fluxes. The TRN of the area occupied by the SJMA has decreased as a consequence of replacing the low land coastal plain vegetation with man made materials, indicating that it takes less energy to raise the surface temperature of the urban area, whereas the TRN of forested regions has remained virtually unchanged. The global warming signal also has effects on the thermal response of the SJMA, where dryer current conditions generate lower TRN values. Differences due to global warming are more evident in the Bowen ratio pattern, mostly associated with the drier present conditions observed and its effects on sensible and latent heat fluxes. In terms of testing different mitigation strategies, the short-term simulations show that the urban area is more efficient in partitioning surface energy balance terms when green roofs are specified, as opposed to including vegetation inside the urban core.

Comarazamy, Daniel↗
NASA NTRS1999

Characterization of the Heat Extraction Capability of a Compliant, Sliding, Thermal Interface for Use in a High Temperature, Vacuum Microgravity Furnace

Much of the material science gained in microgravity research requires processing a sample in a high performance furnace. One such furnace currently being designed is the Quench Module Insert (QMI). The Principle Investigators, for whom the furnace is proposed, require high temperature gradients in their cylindrical samples to achieve the science objectives. One of the components critical to achieving high sample axial temperature gradients in the Quench Module Insert is a high performance cold zone to extract the heat from the sample. This cold zone employs a compliant, sliding thermal interface based on a Vel-Therm felt. This felt provides a conductive path between the Sample Cartridge Assembly (SCA) exterior surface and the interior surface of the water cooled chill block while allowing movement of the sample relative to the chill block. The Vel-Therm felt is composed of long polymer-based fibers affixed to a thin flexible substrate layer. The fibers are oriented perpendicular to this substrate giving the felt the appearance of a velvet fabric. The Vel-Therm felt heat extraction capability was quantified in earlier tests performed in an inert gas environment. The current activity, described in this paper, is intended to characterize the extraction capability of Vel-Therm felt in a vacuum environment similar to the QMI environment. This testing is necessary to quantify the thermal performance of the Vel-Therm felt and the sensitivity of that performance to key variables. The data derived from these tests will be incorporated into the current thermal models to improve the quality of the models and reduce uncertainty of the analytical results. In addition, the data will be used to help select the appropriate Vel-Therm felt and set proper operating limits as well as assess the performance range of the furnace. The objective of this test is to measure the heat extraction rate of the Vel-Therm felt as specified by the effective heat transfer coefficient. Therefore, the test setup was designed to force the bulk of the heat transfer through the area where the Vel-Therm felt was applied. A heat source, consisting of a ceramic heating element encased in a Copper (Cu) housing is mounted on four isolated support rods. A 6-layer molybdenum radiation shield is used to insulate against heat loss from the heater and prevent heat exchange between the hot and cold sides of the test apparatus. The Vel-Therm felt is affixed to the surface of the cold sink, a water-cooled Cu chill block. An adjustable plate supports and isolates the cold sink from above and is used to control the amount of the deflection of the Vel- Therm when in contact with the Cu heating element housing. The primary means of establishing the power being conducted through the felt is to measure the energy being transferred to the water passing through the chill block. Analysis was performed to support the assumption that the source and sink surfaces were approximately isothermal under the specific test conditions. As a check on the amount of power passing through the felt, the power supplied to the heater was also measured. Thermocouples were strategically located throughout the test apparatus for measurement purposes. A bell jar was lowered over the assembly to impose vacuum conditions. Currently, variables tested have been fiber compression and fiber type and surface temperatures (both the hot and cold side temperatures are hypothesized to be important to the performance of the Vel-Therm.) Selected runs were repeated to ensure consistency and repeatability. Results obtained thus far reveal that Vel-Therm performance is significantly degraded by fibers being exposed to high compression. It also shows that performance is somewhat negatively impacted by previous compression, thereby, raising the question of repeatability. In addition, early results show a significant dependence on temperature. A computer aided mathematical analysis of the test setup is ongoing. The results will be correlated to actual results. The correlation will examine such details as parasitic loses, conduction down the power leads and many other concerns.

Bellomy-Ezell, Jenny↗
NASA NTRSDate not supplied

A Holistic DSMC Transport Database for Re-Entry and Ablation Modeling

Hybrid simulation frameworks combining Computational Fluid Dynamics (CFD) and Direct Simulation Monte Carlo (DSMC) are frequently employed to efficiently perform high-fidelity solutions of environments containing combined continuum/rarified flow. The use of DSMC, a stochastic, particle-based method, is necessary for high-Knudsen flow where continuum-based assumptions governing CFD break down. However, the DSMC methodology is generally very computationally inefficient to model the continuum regime. In a CFD/DSMC hybrid approach, obtaining an accurate, high-fidelity solution hinges on the consistent treatment of transport properties and the used thermo-chemical models employed within the two solvers. In principle, in regions where CFD and DSMC are both employed, the same gas mixture under the same conditions should have the same properties, regardless of simulation type. Observed differences should be due to non-equilibrium processes, rather than differences in physical models. While the transport models governing CFD and DSMC simulations are starkly different, they can effectively be linked via their use of reduced Chapman-Enskog collision integrals. In CFD, these integrals are typically stored as fitted polynomial expressions and used to directly compute gas transport properties via mixing rules or the full Chapman-Enskog formulation. In DSMC, they can be used to derive the collision parameters needed for the phenomenological collision cross-section models that govern particle interactions, via a Nelder-Mead optimization scheme. The goal of this work is to provide a unified DSMC transport database encompassing the vast majority of known gas species encountered during atmospheric entry, on Earth or any other Solar body. This goal is largely possible due to recently performed ab-initio quantum chemistry calculations. Combined with other high-fidelity literature sources, the planned database will consist of collision integral data for over 200 neutral and ionized species and over 17000 binary collisions. From these collision integrals, Nelder-Mead optimization is used to compute Variable Soft Sphere (VSS) collision model parameters for DSMC, fitted from 300 K to 20000 K. Initial comparisons of transport properties of relevant equilibrium gas mixtures show great agreement between CFD and DSMC-derived results. The completed database will be able to be readily applied to model binary collisions of any gas mixture containing the included species over the specified temperature range, making it a valuable tool for future planetary probe modeling efforts. An example is shown below. Equilibrium mixture transport properties for a 19-species Titan atmospheric model [4] are computed using both fitted VSS parameters and the original CFD collision integral values. Deviations in computed properties between the two approaches is less than 5% for the entire temperature range.

M R Gosma↗
NASA NTRS1995

Solid surface combustion experiment flame spread in a quiescent, microgravity environment implications of spread rate and flame structure

A unique environment in which flame spreading, a phenomenon of fundamental, scientific interest, has importance to fire safety is that of spacecraft in which the gravitational acceleration is low compared with that of the Earth, i.e., microgravity. Experiments aboard eight Space Shuttle missions between October 1990 and February 1995 were conducted using the Solid Surface Combustion Experiment (SSCE) payload apparatus in an effort to determine the mechanisms of gas-phase flame spread over solid fuel surfaces in the absence of any buoyancy induced or externally imposed oxidizer flow. The overall SSCE effort began in December of 1984. The SSCE apparatus consists of a sealed container, approximately 0.039 cu m, that is filled with a specified O2/N2 mixture at a prescribed pressure. Five of the experiments used a thin cellulosic fuel, ashless filter paper, 3 cm wide x 10 cm long, 0.00825 cm half-thickness, ignited in five different ambient conditions. Three of the experiments, the most recent, used thick polymethylmethacrylate (PMMA) samples 0.635 cm wide x 2 cm long, 0.32 cm half-thickness. Three experiments, STS 41, 40 and 43, were designed to evaluate the effect of ambient pressure on flame spread over the thin cellulosic fuel while flights STS 50 and 47 were at the same pressure as two of the earlier flights but at a lower oxygen concentration in order to evaluate the effect of ambient oxygen level on the flame spread process at microgravity. For the PMMA flights, two experiments, STS 54 and 63, were at the same pressure but different oxygen concentrations while STS 64 was at the same oxygen concentration as STS 63 but at a higher pressure. Two orthogonal views of the experiments were recorded on 16 mm cine-cameras operating at 24 frames/s. In addition to filmed images of the side view of the flames and surface view of the burning samples, solid- and gas-phase temperatures were recorded using thermocouples. The experiment is battery powered and follows an automated sequence upon activation by the Shuttle Crew. In this study we separate the SSCE data into two groups according to the fuel type: (1) thin cellulose; and (2) thick PMMA. The experimental spread rates are compared with prediction from a number of models in an effort to uncover the important physics that characterize microgravity flame spread. Both steady and unsteady solutions are employed to explore the flame evolution, especially for thick fuels. Finally, the flame structure in downward spread is compared with the microgravity flame structure and modeling results to delineate the difference between the two configurations and the influence of normal gravity.

Bundy, Matthew↗
NASA NTRS2017

Summary of Data from the Sixth AIAA CFD Drag Prediction Workshop: CRM Cases 2 to 5

Results from the Sixth AIAA CFD Drag Prediction Workshop Common Research Model Cases 2 to 5 are presented. As with past workshops, numerical calculations are performed using industry-relevant geometry, methodology, and test cases. Cases 2 to 5 focused on force/moment and pressure predictions for the NASA Common Research Model wing-body and wing-body-nacelle-pylon configurations, including Case 2 - a grid refinement study and nacelle-pylon drag increment prediction study; Case 3 - an angle-of-attack buffet study; Case 4 - an optional wing-body grid adaption study; and Case 5 - an optional wing-body coupled aero-structural simulation. The Common Research Model geometry differed from previous workshops in that it was deformed to the appropriate static aeroelastic twist and deflection at each specified angle-of-attack. The grid refinement study used a common set of overset and unstructured grids, as well as user created Multiblock structured, unstructured, and Cartesian based grids. For the supplied common grids, six levels of refinement were created resulting in grids ranging from 7x10(exp 6) to 208x10(exp 6) cells. This study (Case 2) showed further reduced scatter from previous workshops, and very good prediction of the nacelle-pylon drag increment. Case 3 studied buffet onset at M=0.85 using the Medium grid (20 to 40x10(exp 6) nodes) from the above described sequence. The prescribed alpha sweep used finely spaced intervals through the zone where wing separation was expected to begin. Although the use of the prescribed aeroelastic twist and deflection at each angle-of-attack greatly improved the wing pressure distribution agreement with test data, many solutions still exhibited premature flow separation. The remaining solutions exhibited a significant spread of lift and pitching moment at each angle-of-attack, much of which can be attributed to excessive aft pressure loading and shock location variation. Four Case 4 grid adaption solutions were submitted. Starting with grids less than 2x10(exp 6) grid points, two solutions showed a rapid convergence to an acceptable solution. Four Case 5 coupled aerostructural solutions were submitted. Both showed good agreement with experimental data. Results from this workshop highlight the continuing need for CFD improvement, particularly for conditions with significant flow separation. These comparisons also suggest the need for improved experimental diagnostics to guide future CFD development.

Tinoco, Edward N.↗
NASA NTRS1996

SSME Condition Monitoring Using Neural Networks and Plume Spectral Signatures

For a variety of reasons, condition monitoring of the Space Shuttle Main Engine (SSME) has become an important concern for both ground tests and in-flight operation. The complexities of the SSME suggest that active, real-time condition monitoring should be performed to avoid large-scale or catastrophic failure of the engine. In 1986, the SSME became the subject of a plume emission spectroscopy project at NASA's Marshall Space Flight Center (MSFC). Since then, plume emission spectroscopy has recorded many nominal tests and the qualitative spectral features of the SSME plume are now well established. Significant discoveries made with both wide-band and narrow-band plume emission spectroscopy systems led MSFC to develop the Optical Plume Anomaly Detection (OPAD) system. The OPAD system is designed to provide condition monitoring of the SSME during ground-level testing. The operational health of the engine is achieved through the acquisition of spectrally resolved plume emissions and the subsequent identification of abnormal emission levels in the plume indicative of engine erosion or component failure. Eventually, OPAD, or a derivative of the technology, could find its way on to an actual space vehicle and provide in-flight engine condition monitoring. This technology step, however, will require miniaturized hardware capable of processing plume spectral data in real-time. An objective of OPAD condition monitoring is to determine how much of an element is present in the SSME plume. The basic premise is that by knowing the element and its concentration, this could be related back to the health of components within the engine. For example, an abnormal amount of silver in the plume might signify increased wear or deterioration of a particular bearing in the engine. Once an anomaly is identified, the engine could be shut down before catastrophic failure occurs. Currently, element concentrations in the plume are determined iteratively with the help of a non-linear computer code called SPECTRA, developed at the USAF Arnold Engineering Development Center. Ostensibly, the code produces intensity versus wavelength plots (i.e., spectra) when inputs such as element concentrations, reaction temperature, and reaction pressure are provided. However, in order to provide a higher-level analysis, element concentration is not specified explicitly as an input. Instead, two quantum variables, number density and broadening parameter, are used. Past experience with OPAD data analysis has revealed that the region of primary interest in any SSME plume spectrum lies in the wavelength band of 3300 A to 4330 A. Experience has also revealed that some elements, such as iron, cobalt and nickel, cause multiple peaks over the chosen wavelength range whereas other elements (magnesium, for example) have a few, relatively isolated peaks in the chosen wavelength range. Iteration with SPECTRA as a part of OPAD data analysis is an incredibly labor intensive task and not one to be performed by hand. What is really needed is the "inverse" of the computer code but the mathematical model for the inverse mapping is tenuous at best. However, building generalized models based upon known input/output mappings while ignoring details of the governing physical model is possible using neural networks. Thus the objective of the research project described herein was to quickly and accurately predict combustion temperature and element concentrations (i.e., number density and broadening parameter) from a given spectrum using a neural network. In other words, a neural network had to be developed that would provide a generalized "inverse" of the computer code SPECTRA.

Hopkins, Randall↗
NASA NTRS1993

Liquid hydrogen and liquid oxygen feedline passive recirculation analysis

The primary goal of the National Launch System (NLS) program was to design an operationally efficient, highly reliable vehicle with minimal recurring launch costs. To achieve this goal, trade studies of key main propulsion subsystems were performed to specify vehicle design requirements. These requirements include the use of passive recirculation to thermally condition the liquid hydrogen (LH2) and liquid oxygen (LO2) propellant feed systems and Space Transportation Main Engine (STME) fuel pumps. Rockwell International (RI) proposed a joint independent research and development (JIRAD) program with Marshall Space Flight Center (MSFC) to study the LH2 feed system passive recirculation concept. The testing was started in July 1992 and completed in November 1992. Vertical and sloped feedline designs were used. An engine simulator was attached at the bottom of the feedline. This simulator had strip heaters that were set to equal the corresponding heat input from different engines. A computer program is currently being used to analyze the passive recirculation concept in the LH2 vertical feedline tests. Four tests, where the heater setting is the independent variable, were chosen. While the JIRAD with RI was underway, General Dynamics Space Systems (GDSS) proposed a JIRAD with MSFC to explore passive recirculation in the LO2 feed system. Liquid nitrogen (LN2) is being used instead of LO2 for safety and economic concerns. To date, three sets of calibration tests have been completed on the sloped LN2 test article. The environmental heat was calculated from the calibration tests in which the strip heaters were turned off. During the LH2 testing, the environmental heat was assumed to be constant. Therefore, the total heat was equal to the environmental heat flux plus the heater input. However, the first two sets of LN2 calibration tests have shown that the environmental heat flux varies with heater input. A Systems Improved Numerical Differencing Analyzer and Fluid Integrator (SINDA/FLUINT) model is currently being built to determine if this variation in environmental heat is due to a change in the wall temperature. During the third set of calibration tests, a faulty reference junction was found. Based on this anomaly with the reference junction, the heat flux calculations from the first two calibration sets are now considered questionable.

Holt, Kimberly Ann↗