Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “transient dynamic 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 613 records · Page 34

Development, Verification and Validation of Parallel, Scalable Volume of Fluid CFD Program for Propulsion Applications

There are many instances involving liquid/gas interfaces and their dynamics in the design of liquid engine powered rockets such as the Space Launch System (SLS). Some examples of these applications are: Propellant tank draining and slosh, subcritical condition injector analysis for gas generators, preburners and thrust chambers, water deluge mitigation for launch induced environments and even solid rocket motor liquid slag dynamics. Commercially available CFD programs simulating gas/liquid interfaces using the Volume of Fluid approach are currently limited in their parallel scalability. In 2010 for instance, an internal NASA/MSFC review of three commercial tools revealed that parallel scalability was seriously compromised at 8 cpus and no additional speedup was possible after 32 cpus. Other non-interface CFD applications at the time were demonstrating useful parallel scalability up to 4,096 processors or more. Based on this review, NASA/MSFC initiated an effort to implement a Volume of Fluid implementation within the unstructured mesh, pressure-based algorithm CFD program, Loci-STREAM. After verification was achieved by comparing results to the commercial CFD program CFD-Ace+, and validation by direct comparison with data, Loci-STREAM-VoF is now the production CFD tool for propellant slosh force and slosh damping rate simulations at NASA/MSFC. On these applications, good parallel scalability has been demonstrated for problems sizes of tens of millions of cells and thousands of cpu cores. Ongoing efforts are focused on the application of Loci-STREAM-VoF to predict the transient flow patterns of water on the SLS Mobile Launch Platform in order to support the phasing of water for launch environment mitigation so that vehicle determinantal effects are not realized.

West, Jeff↗

Acoustic Database for Turbofan Engine Core-Noise Sources: Volume - I

In this program, a database of dynamic temperature and dynamic pressure measurements were acquired inside the core of a TECH977 turbofan engine to support investigations of indirect combustion noise. Dynamic temperature and pressure measurements were recorded for engine gas dynamics up to temperatures of 3100 degrees Fahrenheit and transient responses as high as 1000 hertz. These measurements were made at the entrance of the high pressure turbine (HPT) and at the entrance and exit of the low pressure turbine (LPT). Measurements were made at two circumferential clocking positions. In the combustor and inter-turbine duct (ITD), measurements were made at two axial locations to enable the exploration of time delays. The dynamic temperature measurements were made using dual thin-wire thermocouple probes. The dynamic pressure measurements were made using semi-infinite probes. Prior to the engine test, a series of bench, oven, and combustor rig tests were conducted to characterize the performance of the dual wire temperature probes and to define and characterize the data acquisition systems. A measurement solution for acquiring dynamic temperature and pressure data on the engine was defined. A suite of hardware modifications were designed to incorporate the dynamic temperature and pressure instrumentation into the TECH977 engine. In particular, a probe actuation system was developed to protect the delicate temperature probes during engine startup and transients in order to maximize sensor life. A set of temperature probes was procured and the TECH977 engine was assembled with the suite of new and modified hardware. The engine was tested at four steady state operating speeds, with repeats. Dynamic pressure and temperature data were acquired at each condition for at least one minute. At the two highest power settings, temperature data could not be obtained at the forward probe locations since the mean temperatures exceeded the capability of the probes. The temperature data were processed using software that accounts for the effects of convective and conductive heat transfer. The software was developed under previous NASA sponsored programs. Compensated temperature spectra and compensated time histories corresponding to the dynamic temperature of the gas stream were generated. Auto-spectral and cross-spectral analyses of the data were performed to investigate spectral features, acoustic circumferential mode content, signal coherence, and time delays. The dynamic temperature data exhibit a wideband and fairly flat spectral content. The temperature spectra do not change substantially with operating speed. The pressure spectra in the combustor and ITD exhibit generally similar shapes and amplitudes, making it difficult to identify any features that suggest the presence of indirect combustion noise. Cross-spectral analysis reveal a strong correlation between pressure and temperature fluctuations in the ITD, but little correlation between temperature fluctuations at the entrance of the HPT and pressure fluctuations downstream of it. Temperature fluctuations at the entrance of the low pressure turbine were an order of magnitude smaller than those at the entrance to the high pressure turbine. Time delay analysis of the temperature fluctuations in the combustor was inconclusive, perhaps due to the substantial mixing that occurs between the upstream and downstream locations. Time delay analysis of the temperature fluctuations in the ITD indicate that they convect at the mean flow speed. Analysis of the data did not reveal any convincing indications of the presence of indirect combustion noise. However, this analysis has been preliminary and additional exploration of the data is recommended including the use of more sophisticated signal processing to explore subtle issues that have been revealed but which are not yet fully understood or explained.

Aeroacoustics↗

Basic studies of baroclinic flows

A fully nonlinear 3-dimensional numerical model (GEOSIM), previously developed and validated for several cases of geophysical fluid flow, has been used to investigate the dynamical behavior of laboratory experiments of fluid flows similar to those of the Earth's atmosphere. The phenomena investigated are amplitude vacillation, and the response of the fluid system to uneven heating and cooling. The previous year's work included hysteresis in the transition between axisymmetric and wave flow. Investigation is also continuing of the flows in the Geophysical Fluid Flow Cell (GFFC), a low-gravity Spacelab experiment. Much of the effort in the past year has been spent in validation of the model under a wide range of external parameters including nonlinear flow regimes. With the implementation of a 3-dimensional upwind differencing scheme, higher spectral resolution, and a shorter time step, the model has been found capable of predicting the majority of flow regimes observed in one complete series of baroclinic annulus experiments of Pfeffer and co-workers. Detailed analysis of amplitude vacillation has revealed that the phase splitting described in the laboratory experiments occurs in some but not all cases. Through the use of animation of the models output, a vivid 3-dimensional view of the phase splitting was shown to the audience of the Southeastern Geophysical Fluid Dynamics Conference in March of this year. A study on interannual variability was made using GEOSIM with periodic variations in the thermal forcing. Thus far, the model has not predicted a chaotic behavior as observed in the experiments, although there is a sensitivity in the wavenumber selection to the initial conditions. Work on this subject, and on annulus experiments with non-axisymmetric thermal heating, will continue. The comparison of GEOSIM's predictions will result from the Spacelab 3 GFFC experiments continued over the past year, on a 'back-burner' basis. At this point, the study (in the form of a draft of a journal article) is nearly completed. The results from GEOSIM compared very well with the experiments, and the use of the model allows the demonstration of flow mechanics that were not possible with the experimental data. For example, animation of the model output shows that the forking of the spiral bands is a transient phenomenon, due to the differential east-west propagation of convection bands from different latitudes.

Miller, Timothy L.↗

Features of steady magnetospheric convection

The large-scale patterns of ionospheric convection and particle precipitation are described during two intervals of steady magnetospheric convection (SMC) on November 24, 1981. The unique data set used in the analysis includes recordings from the worldwide network of magnetometers and all-sky cameras, global auroral images from the Dynamics Explorer (DE) 1 spacecraft, and particle precipitation data from low-altitude National Oceanic and Atmospheric Administration (NOAA) 6 and NOAA 7 spacecraft. The data show that intense magnetospheric convection continued during more than 10 hours under the steady southward interplanetary magnetic field without any distinct substorm signatures. All data sets available confirmed the stable character of the large-scale magnetospheric configuration during this period. In particular, the magnetic flux threading the polar cap was stable (within 10%) during 3.5 hours of continued DE 1 observations. The dayside cusp was located at an unusually low latitude (70 deg CGL). The nightside auroral pattern consisted of two distinct regions. The diffuse aurora in the equatorward half of the expanded (10 deg wide) auroral oval was well-separated from the bright, active auroral forms found in the vicinity of the poleward boundary of the oval. The twin-vortex convection pattern had no signature of the Harang discontinuity; its nightside 'convection throat' was spatially coincident with the poleward active auroras. This region of the auroral oval was identified as the primary site of the short-lived transient activations during the SMC intervals. The energetic particle observations show that the auroral precipitation up to its high-latitude limit is on closed field lines and that particle acceleration up to greater than 30-keV energy starts close to this limit. The isotropic boundaries of the greater than 30-keV protons and electrons were found close to each other, separating regions of discrete and diffuse precipitation. This suggests that these precipitation types originate on the very taillike and very dipolelike field lines, respectively.

Yahnin, A.↗

Damage Based Analysis (DBA): Theory, Derivation and Practical Application - Using Both an Acceleration and Pseudo-Velocity Approach

The objective of this presentation is to give a brief overview of the theory behind the (DBA) method, an overview of the derivation and a practical application of the theory using the Python computer language. The Theory and Derivation will use both Acceleration and Pseudo Velocity methods to derive a series of equations for processing by Python. We will take the results and compare both Acceleration and Pseudo Velocity methods and discuss implementation of the Python functions. Also, we will discuss the efficiency of the methods and the amount of computer time required for the solution. In conclusion, (DBA) offers a powerful method to evaluate the amount of energy imparted into a system in the form of both Amplitude and Duration during qualification testing and flight environments. Many forms of steady state and transient vibratory motion can be characterized using this technique. (DBA) provides a more robust alternative to traditional methods such Power Spectral Density (PSD) using a Maximax approach.

Pseudo Velocity↗

Data and Modeling of Dendrites Subject to A Step Change in Pressure (TDSE)

There is considerable interest in dendritic solidification because of the influence dendrites have in the determination of microstructure, and thereby in the physical properties of cast metals and alloys. Current theories and models of dendritic growth generally couple diffusion effects in the melt with the physics of the interface. Data and subsequent analysis prior of the tip growth speed and radii of thermal succinonitrile dendrites in the near-convection free, on-orbit, free-fall environment demonstrate that these theories yield predictions that are reasonably in agreement with the results of experiment. However, data and analysis for assessing the interfacial physics component of theory are not sufficiently detailed or definitive. To study fundamental aspects of dendritic interface stability, we are measuring and modeling the kinetics and morphology of dendrites as they evolve from one well-defined steady state at a pre-set supercooling, through a transient stage, to a different well-defined steady state. More specifically, we subject succinonitrile dendrites, growing under steady-state conditions, to a rapid change in pressure. This leads to a rapid change in thermal driving force from the corresponding change in both the equilibrium melting temperature due to the Clapeyron effect, and a change in the far-field temperature due to adiabatic temperature changes in the bulk liquid and solid. Subsequently, we observe transformations from a well-characterized initial state into a new steady-state. Initial data reveal that the dendrite tip velocity changes almost as fast as the pressure charges, while the tip radius changes occur more slowly, taking from 10 60 seconds depending on the size of the step change and the final supercooling. Computer modeling of this process shows both agreements and disagreements with the experimental data. In making these observations and measurements, we are gaining new understandings of interfacial dynamics and state-selection physics.

Koss, Matthew B.↗

Thermosphere Global Time Response to Geomagnetic Storms Caused by Coronal Mass Ejections

We investigate, for the first time with a spatial superposed epoch analysis study, the thermosphere global time response to 159 geomagnetic storms caused by coronal mass ejections (CMEs) observed in the solar wind at Earth's orbit during the period of September 2001 to September 2011. The thermosphere neutral mass density is obtained from the CHAMP (CHAllenge Mini-Satellite Payload) and GRACE (Gravity Recovery Climate Experiment) spacecraft. All density measurements are intercalibrated against densities computed by the Jacchia-Bowman 2008 empirical model under the regime of very low geomagnetic activity. We explore both the effects of the pre-CME shock impact on the thermosphere and of the storm main phase onset by taking their times of occurrence as zero epoch times (CME impact and interplanetary magnetic field Bz southward turning) for each storm. We find that the shock impact produces quick and transient responses at the two high-latitude regions with minimal propagation toward lower latitudes. In both cases, thermosphere is heated in very high latitude regions within several minutes. The Bz southward turning of the storm onset has a fast heating manifestation at the two high-latitude regions, and it takes approximately 3 h for that heating to propagate down to equatorial latitudes and to globalize in the thermosphere. This heating propagation is presumably accomplished, at least in part, with traveling atmospheric disturbances and complex meridional wind structures. Current models use longer lag times in computing thermosphere density dynamics during storms. Our results suggest that the thermosphere response time scales are shorter and should be accordingly adjusted in thermospheric empirical models.

mass ejections↗

Interannual Variability in Soil Trace Gas (CO2, N2O, NO) Fluxes and Analysis of Controllers

Interannual variability in flux rates of biogenic trace gases must be quantified in order to understand the differences between short-term trends and actual long-term change in biosphere-atmosphere interactions. We simulated interannual patterns (1983-1988) of global trace gas fluxes from soils using the NASA Ames model version of CASA (Carnegie-Ames-Stanford Approach) in a transient simulation mode. This ecosystem model has been recalibrated for simulations driven by satellite vegetation index data from the NOAA Advanced Very High Resolution Radiometer (AVHRR) over the mid-1980s. The predicted interannual pattern of soil heterotropic CO2 emissions indicates that relatively large increases in global carbon flux from soils occurred about three years following the strong El Nino Southern Oscillation (ENSO) event of 1983. Results for the years 1986 and 1987 showed an annual increment of +1 Pg (1015 g) C-CO2 emitted from soils, which tended to dampen the estimated global increase in net ecosystem production with about a two year lag period relative to plant carbon fixation. Zonal discrimination of model results implies that 80-90 percent of the yearly positive increments in soil CO2 emission during 1986-87 were attributable to soil organic matter decomposition in the low-latitudes (between 30 N and 30 S). Soils of the northern middle-latitude zone (between 30 N and 60 N) accounted for the residual of these annual increments. Total annual emissions of nitrogen trace gases (N2O and NO) from soils were estimated to vary from 2-4 percent over the time period modeled, a level of variability which is consistent with predicted interannual fluctuations in global soil CO2 fluxes. Interannual variability of precipitation in tropical and subtropical zones (30 N to 20 S appeared to drive the dynamic inverse relationship between higher annual emissions of NO versus emissions of N2O. Global mean emission rates from natural (heterotrophic) soil sources over the period modeled (1983-1988) were estimated at 57.1 Pg C-CO2yr-1, 9.8Tg (1012 g) N-NO yr-1, and 9.7 Tg N-N2O yr-1. Chemical fertilizer contributions to global soil N gas fluxes were estimated at between 1.3 to 7.3 Tg N-NO yr-1, and 1.2 to 4.0 Tg N-N2O yr-1.

Potter, C.↗

Report on New Mission Concept Study: Stereo X-Ray Corona Imager Mission

Studies of the three-dimensional structure and dynamics of the solar corona have been severely limited by the constraint of single viewpoint observations. The Stereo X-Ray Coronal Imager (SXCI) mission will send a single instrument, an X-ray telescope, into deep space expressly to record stereoscopic images of the solar corona. The SXCI spacecraft will be inserted into a approximately 1 AU heliocentric orbit leading Earth by approximately 25 deg at the end of nine months. The SXCI X-ray telescope forms one element of a stereo pair, the second element being an identical X-ray telescope in Earth orbit placed there as part of the NOAA GOES program. X-ray emission is a powerful diagnostic of the corona and its magnetic fields, and three dimensional information on the coronal magnetic structure would be obtained by combining the data from the two X-ray telescopes. This information can be used to address the major solar physics questions of (1) what causes explosive coronal events such as coronal mass ejections (CMEs), eruptive flares and prominence eruptions and (2) what causes the transient heating of coronal loops. Stereoscopic views of the optically thin corona will resolve some ambiguities inherent in single line-of-sight observations. Triangulation gives 3D solar coordinates of features which can be seen in the simultaneous images from both telescopes. As part of this study, tools were developed for determining the 3D geometry of coronal features using triangulation. Advanced technologies for visualization and analysis of stereo images were tested. Results of mission and spacecraft studies are also reported.

Liewer, Paulett C.↗

Spray Combustion Modeling with VOF and Finite-Rate Chemistry

A spray atomization and combustion model is developed based on the volume-of-fluid (VOF) transport equation with finite-rate chemistry model. The gas-liquid interface mass, momentum and energy conservation laws are modeled by continuum surface force mechanisms. A new solution method is developed such that the present VOF model can be applied for all-speed range flows. The objectives of the present study are: (1) to develop and verify the fractional volume-of-fluid (VOF) cell partitioning approach into a predictor-corrector algorithm to deal with multiphase (gas-liquid) free surface flow problems; (2) to implement the developed unified algorithm in a general purpose computational fluid dynamics (CFD) code, Finite Difference Navier-Stokes (FDNS), with droplet dynamics and finite-rate chemistry models; and (3) to demonstrate the effectiveness of the present approach by simulating benchmark problems of jet breakup/spray atomization and combustion. Modeling multiphase fluid flows poses a significant challenge because a required boundary must be applied to a transient, irregular surface that is discontinuous, and the flow regimes considered can range from incompressible to highspeed compressible flows. The flow-process modeling is further complicated by surface tension, interfacial heat and mass transfer, spray formation and turbulence, and their interactions. The major contribution of the present method is to combine the novel feature of the Volume of Fluid (VOF) method and the Eulerian/Lagrangian method into a unified algorithm for efficient noniterative, time-accurate calculations of multiphase free surface flows valid at all speeds. The proposed method reformulated the VOF equation to strongly couple two distinct phases (liquid and gas), and tracks droplets on a Lagrangian frame when spray model is required, using a unified predictor-corrector technique to account for the non-linear linkages through the convective contributions of VOF. The discontinuities within the sharp interface will be modeled as a volume force to avoid stiffness. Formations of droplets, tracking of droplet dynamics and modeling of the droplet breakup/evaporation, are handled through the same unified predictor-corrector procedure. Thus the new algorithm is non-iterative and is flexible for general geometries with arbitrarily complex topology in free surfaces. The FDNS finite-difference Navier-Stokes code is employed as the baseline of the current development. Benchmark test cases of shear coaxial LOX/H2 liquid jet with atomization/combustion and impinging jet test cases are investigated in the present work. Preliminary data comparisons show good qualitative agreement between data and the present analysis. It is indicative from these results that the present method has great potential to become a general engineering design analysis and diagnostics tool for problems involving spray combustion.

Chen, Yen-Sen↗

Control Design Strategies to Enhance Long-Term Aircraft Structural Integrity

Over the operational lifetime of both military and civil aircraft, structural components are exposed to hundreds of thousands of low-stress repetitive load cycles and less frequent but higher-stress transient loads originating from maneuvering flight and atmospheric gusts. Micro-material imperfections in the structure, such as cracks and debonded laminates, expand and grow in this environment, reducing the structural integrity and shortening the life of the airframe. Extreme costs associated with refurbishment of critical load-bearing structural components in a large fleet, or altogether reinventoring the fleet with newer models, indicate alternative solutions for life extension of the airframe structure are highly desirable. Increased levels of operational safety and reliability are also important factors influencing the desirability of such solutions. One area having significant potential for impacting crack growth/fatigue damage reduction and structural life extension is flight control. To modify the airframe response dynamics arising from command inputs and gust disturbances, feedback loops are routinely applied to vehicles. A dexterous flight control system architecture senses key vehicle motions and generates critical forces/moments at multiple points distributed throughout the airframe to elicit the desired motion characteristics. In principle, these same control loops can be utilized to influence the level of exposure to harmful loads during flight on structural components. Project objectives are to investigate and/or assess the leverage control has on reducing fatigue damage and enhancing long-term structural integrity, without degrading attitude control and trajectory guidance performance levels. In particular, efforts have focused on the effects inner loop control parameters and architectures have on fatigue damage rate. To complete this research, an actively controlled flexible aircraft model and a new state space modeling procedure for crack growth have been utilized. Analysis of the analytical state space model for crack growth revealed the critical mathematical factors, and hence the physical mechanism they represent, that influenced high rates of airframe crack growth. The crack model was then exercised with simple load inputs to uncover and expose key crack growth behavior. To characterize crack growth behavior, both "short-term" laboratory specimen test type inputs and "long-term" operational flight type inputs were considered. Harmonic loading with a single overload revealed typical exponential crack growth behavior until the overload application, after which time the crack growth was retarded for a period of time depending on the overload strength. An optimum overload strength was identified which leads to maximum retardation of crack growth. Harmonic loading with a repeated overload of varying strength and frequency again revealed an optimum overload trait for maximizing growth retardation. The optimum overload strength ratio lies near the range of 2 to 3 with dependency on frequency. Experimental data was found to correlate well with the analytical predictions.

Newman, Brett A.↗

Time-Frequency Analysis of Experimental and Analytical Hub Loads of a Rotor Undergoing a Rotor Speed Change

A two-part analytical study was conducted examining a small-scale, two-bladed rotor undergoing a change in rotor speed in a wind tunnel setting. First, a parametric study was conducted to understand how rotor transient response, determined using a time-marching solution, is affected by the choice of rotor wake model, blade elasticity, and inclusion of an elastic support structure model. The second part of the study compared results from an analytical model of the NASA Multirotor Test Bed (MTB) undergoing a rotor speed change to experimental data. Since the analytical and experimental time histories were non-stationary signals, a Stockwell transform was employed to analyze the data in lieu of traditional Fourier transform-based methods. The analysis included extraction of time-varying frequency content of the rotor thrust and hub motion, damping ratios, instantaneous rotor speed, and instantaneous phase difference between thrust and hub motion. Use of a free wake model was found to be necessary to predict higher-harmonic thrust, however, it underpredicted the amplitude of the unsteady loads. Modeling the elasticity of the stiff rotor blades installed on the MTB resulted in minimal impact on the vibratory loads compared to a rigid blade model. The current analytical model of the MTB, including an elastic support structure, properly predicted the measured trends in the frequency content of the thrust.

Rotor dynamics↗

Hydrogen leak detection in the Space Shuttle

This study focuses on a helium gas jet flowing into room air. Measurements of helium concentration and velocity in the jet-air mixture are reported. The objective is to learn about jet characteristics so that dynamically similar hydrogen leaks may be located in the Space Shuttle. The hazardous gas detection system (HGDS) in the mobile launch pad uses mass spectrometers to monitor the shuttle environment for leaks. The mass spectrometers are fed by long sample tubes which draw gas from the payload bay, mid body, aft engine compartment and external tank. The overall purpose of this study is to improve the HGDS especially in its potential for locating hydrogen leaks. A rapid-response leak detection experiment was designed, built, and tested, following on the work done in this program last summer. The apparatus included a Perkin Elmer MGA-1200 mass spectrometer and air velocity transducer, both monitored by a Macintosh IIFX computer using LabVIEW software. A jet of helium flowing into the lab air simulated a gas leak. Steady helium or hydrogen-nitrogen jets were logged for concentration and velocity, and the power spectral density of each was computed. Last year, large eddies and vortices were visually seen with Schlieren imaging, and they were detected in the time plots of the various instruments. The response time of the MGA-1200 was found in the range of 0.05 to 0.1 sec. Pulsed concentration waves were clearly detected at 25 cycles per sec by spectral analysis of MGA data. No peaks were detected in the power spectrum, so in the present study, 10 Hz bandwidth-averaged power levels were examined at regular frequency intervals. The practical consequences of last year's study are as follows: sampling frequency should be increased above the present rate of 1 sample per second so that transients could be observed and analyzed with frequency response methods. Many more experiments and conditions were observed in this second summer, including the effects of orifice diameter, jet velocity, sample tube design, radial effects, vertical flow, and low hydrogen concentrations (1 percent). A frequent observation was that the power spectrum, calculated from the Fourier transform of concentration fluctuations, gives a separate piece of information from concentration. Many of the tests suggest that power is high where mixing occurs at the helium-air interface. This fact is apparently independent of the concentration level, which could be high or low, but depends on the sample location relative to the jet (leak) origin, whereas high concentration may be due to a strong leak far away or a small leak close to the sample tube. If the power is low for any concentration level, this would signify helium is arriving at the sample tube by diffusion, not chaotic mixing caused by the jet interaction with air. The practical result is to propose a modification of the HGDL mass spectrometer data sampling and software so that sampling rates could be capable of observing at least 25 Hz fluctuations.

Barile, Ronald G↗

Heat Pipe Vapor Dynamics

The dynamic behavior of the vapor flow in heat pipes is investigated at startup and during operational transients. The vapor is modeled as two-dimensional, compressible viscous flow in an enclosure with inflow and outflow boundary conditions. For steady-state and operating transients, the SIMPLER method is used. In this method a control volume approach is employed on a staggered grid which makes the scheme very stable. It is shown that for relatively low input heat fluxes the compressibility of the vapor flow is low and the SIMPLER scheme is suitable for the study of transient vapor dynamics. When the input heat flux is high or the process under a startup operation starts at very low pressures and temperatures, the vapor is highly compressible and a shock wave is created in the evaporator. It is shown that for a wide range of input heat fluxes, the standard methods, including the SIMPLER scheme, are not suitable. A nonlinear filtering technique, along with the centered difference scheme, are then used for shock capturing as well as for the solution of the cell Reynolds-number problem. For high heat flux, the startup transient phase involves multiple shock reflections in the evaporator region. Each shock reflection causes a significant increase in the local pressure and a large pressure drop along the heat pipe. Furthermore, shock reflections cause flow reversal in the evaporation region and flow circulations in the adiabatic region. The maximum and maximum-averaged pressure drops in different sections of the heat pipe oscillate periodically with time because of multiple shock reflections. The pressure drop converges to a constant value at steady state. However, it is significantly higher than its steady-state value at the initiation of the startup transient. The time for the vapor core to reach steady-state condition depends on the input heat flux, the heat pipe geometry, the working fluid, and the condenser conditions. However, the vapor transient time, for an Na-filled heat pipe is on the order of seconds. Depending on the time constant for the overall system, the vapor transient time may be very short. Therefore, the vapor core may be assumed to be quasi-steady in the transient analysis of a heat pipe operation.

Issacci, Farrokh↗

A Modern Load Relief Guidance Scheme for Space Launch Vehicles

Launch vehicle load relief algorithms are concerned with realizing a reduction of transient bending moments near maximum dynamic pressure. Traditional approaches to load relief typically use inner-loop acceleration feedback to reduce the wind-induced angle of attack. When implemented in the inner loop, load relief bandwidth is necessarily limited by the achievable stability margins, and when acceleration feedback is employed, by the uncertainty associated with structural modes that couple with the body-mounted accelerometer. The structure of inner loop load relief increases the dimensionality of the flight control gain and filter optimization problem. Most importantly, classical load relief laws do not take advantage of high-rate and high-accuracy GPS-aided inertial velocity data that is readily available from modern strap down IMUs. In this paper, a novel load relief guidance scheme is described that uses direct angle-of-attack feedback in a clever mechanization. The steering commands are determined by examining the wind-perturbed dynamics of a launch vehicle with respect to a gravity turn ascent trajectory. An angle of attack estimate is derived from GPS-aided inertial data and pre-launch range wind measurements, and it is shown that a reduction worst-case rigid-body loads can be realized without requiring air data. The algorithm also includes a high-rate navigation data preprocessing scheme that operates directly on the IMU delta-theta and delta-velocity measurements in order to produce a filtered acceleration estimate at the vehicle center of mass. The outer-loop guidance scheme simplifies the design process for the classical inner-loop autopilot. Algorithm performance is demonstrated using Monte Carlo analysis of a representative liquid booster in a production high fidelity launch vehicle simulation.

NESC↗

Analysis and Design of Crew Sleep Station for ISS

This paper details the analysis and design of the Temporary Sleep Station (TeSS) environmental control system for International Space Station (ISS). The TeSS will provide crewmembers with a private and personal space, to accommodate sleeping, donning and doffing of clothing, personal communication and performance of recreational activities. The need for privacy to accommodate these activities requires adequate ventilation inside the TeSS. This study considers whether temperature, carbon dioxide, and humidity within the TeSS remain within crew comfort and safety levels for various expected operating scenarios. Evaluation of these scenarios required the use and integration of various simulation codes. An approach was adapted for this study, whereby results from a particular code were integrated with other codes when necessary. Computational Fluid Dynamics (CFD) methods were used to evaluate the flow field inside the TeSS, from which local gradients for temperature, velocity, and species concentration such as CO (sub 2) could be determined. A model of the TeSS, containing a human, as well as equipment such as a laptop computer, was developed in FLUENT, a finite-volume code. Other factors, such as detailed analysis of the heat transfer through the structure, radiation, and air circulation from the TeSS to the US Laboratory Aisle, where the TeSS is housed, were considered in the model. A complementary model was developed in G189A, a code which has been used by NASA/JSC for environmental control systems analyses since the Apollo program. Boundary conditions were exchanged between the FLUENT and G189A TeSS models. G189A provides human respiration rates to the FLUENT model, while the FLUENT model provides local convective heat transfer coefficients to G189A model. An additional benefit from using an approach with both a systems simulation and CFD model, is the capability to verify the results of each model by comparison to the results of the other model. The G189A and FLUENT models were used to evaluate various ventilation designs for the TeSS over a range of operating conditions with varying crew metabolic load, equipment operating modes, ventilation flow rates, and with the TeSS doors open and closed. Results from the study were instrumental in the optimization of a design for the TeSS ventilation hardware. A special case was considered where failure of the TeSS ventilation system occurred. In this case, a study was conducted in order to determine the time required for the CO (sub 2) concentration inside the TeSS to increase to ISS limit values under transient conditions. A lumped-capacitance code, SINDA-FLUINT was used in this case to provide accurate predictions of the human reaction to the TeSS cabin conditions including core and skin temperatures and body heat storage. A simple two-dimensional CFD model of a crewmember inside the TeSS was developed in FLUENT in order to determine the volume envelope of the respired air from the human, which maintained a minimum velocity profile. This volume was then used in the SINDA-FLUINT model to facilitate the calculations of CO (sub 2) concentrations, dry bulb temperatures and humidity levels inside the TeSS.

Keener, John F.↗

Studying X-Ray Binaries with High Energy Frequency Quasi-Periodic Oscillations

The goal of this investigation is to further our understanding of the dynamics of secreting neutron stars and black holes in the hope of using these systems as probes of the physics of strong gravitational fetus. The main focus of this work has been a multi-year program of simultaneous millisecond X-ray timing and spectral observations carried out with the Rossi X-Ray Timing Explorer (RXTE) to perform the X-ray timing and one of the satellites Asca, BeppoSAX, or Chandra to perform X-ray spectral measurements. With the advent of Chandra, we have extended our work to incLude extragalactic X-ray binaries. We conducted a comprehensive study of the X-ray and radio behavior of the Black Hole Candidate (BHC) X-ray transient XTE J1550-564 using RXTE, Chandra, and the Australian Telescope Compact Array (ATCA). We showed that strong radio emission is associated with major X-ray outbursts involving an X-ray state transition, while a compact radio jet is seen in the low/hard X-ray state found in the outburst decay. Interesting, the total energy required to produce the compact jet may be a substantial fraction of the total accretion energy of the system in that state. We also performed a detailed study of the spectral and timing properties of the decay. In joint RXTE/BeppoSAX observations of the neutron-star X-ray binary Cyg X-2, we discovered a correlation between the timing properties (the frequency of the horizontal branch oscillations) and the properties of a soft, thermal component of the X-ray spectrum. d e showed that more det<iiled measurements of this source may lead to constraints on models of the accretion flow and the emission of X- ray from accreting neutron stars. We have completed analysis of RXTE observations of the X-ray transient SAX J1750.8-2900 made after detection of X-ray bursts from the source with the BeppoSAX Wide-Field Camera. We discovered millisecond oscillations in both the persistent emission and in the X-ray bursts.

Kaaret, P.↗

Towards a Theory of Tropical/Midlatitude Mass Exchange from the Earth's Surface through the Stratosphere

The main findings of this research project have been the following: (1) there is a significant feedback from the stratosphere on tropospheric dynamics, and (2) a detailed analysis of the interaction between tropical and polar wave breaking in controlling stratospheric mixing. Two papers are were written and are included. The first paper is titled, "A New Perspective on the Dynamical Link Between the Stratosphere and Troposphere." Atmospheric processes of tropospheric origin can perturb the stratosphere, but direct feedback in the opposite direction is usually assumed to be negligible, despite the troposphere's sensitivity to changes in the release of wave activity into the stratosphere. Here, however, we present evidence that such a feedback exists and can be significant. We find that if the wintertime Arctic polar stratospheric vortex is distorted, either by waves propagating upward from the troposphere or by eastward-travelling stratospheric waves, then there is a concomitant redistribution of stratospheric potential vorticity that induces perturbations in key meteorological fields in the upper troposphere. The feedback is large despite the much greater mass of the troposphere: it can account for up to half of the geopotential height anomaly at the tropopause. Although the relative strength of the feedback is partly due to a cancellation between contributions to these anomalies from lower altitudes, our results imply that stratospheric dynamics and its feedback on the troposphere are more significant for climate modelling and data assimilation than was previously assumed. The second article is titled "Diagnosing the Polar Excitation of Subtropical Waves in the Stratosphere". The poleward migration of planetary scale tongues of subtropical air has often been associated with intense polar vortex disturbances in the stratosphere. This question of vortex influence is reexamined from a potential vorticity (PV) perspective. Anomalous geopotential height and wind fields associated solely with vortex PV anomalies are derived and their impact on the stratospheric subtropical circulation is evaluated. Combined PV inversion and Contour Advection (CA) calculations indicate that transient large scale disturbances of the polar vortex do have a far reaching impact that extends beyond the midlatitude surf zone all the way to the subtropics. This vortex influence is clearly non-local so that even simple wave 2 distortions that leave the vortex well confined within the midlatitudes are observed to excite subtropical waves. Treating subtropical PV as active tracers also showed that upon entrainment, these large scale tongues of low PV air also influenced the dynamics of their own poleward migration.

Hartley, Dana↗