Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Steady state process”

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 829 records · Page 46
All sourcesNASA NTRSDOE OSTI
NASA NTRS2000

Dynamics in Layer Models of Solid Flame Propagation

Self-propagating high-temperature synthesis (SHS) is a process in which combustion waves, e.g., "solid flames", which are considered here, are employed to synthesize desired materials. Like many other systems, SHS is a pattern forming system. The problem of describing experimentally observed patterns and of predicting new, as yet unobserved, patterns continues to attract the attention of scientists and mathematicians due to the fundamental significance of the phenomena in combustion in particular, and in nonlinear science in general. Here, we analyze the dynamics of solid flame propagation in a 2D region by considering the region to be composed of parallel, identical layers aligned along the direction of propagation and having thermal contact. Each layer is then described by wave propagation in 1D, with the transverse Laplacian replaced by a term describing heat exchange between neighboring layers. This configuration is the simplest model of a 2D system because it accounts, in a simple way, for the principal feature of the problem, i.e., heat exchange between neighbors in the transverse direction. For simplicity, we describe the situation for two layers. Because the layers are identical, uniformly propagating waves in each layer must be identical, independent of the heat exchange rate alpha. When the Zeldovich number Z exceeds a critical value Z(sub c), which depends on alpha, uniformly propagating waves become unstable. The stability diagram for the two coupled layers reproduces that for the full 2D problem after appropriate identification of parameters in the two problems. Depending on parameter values, we determine three different steady-state dynamical behaviors (though additional behaviors are also expected to occur). The three behaviors are: (i) waves in each layer which pulsate in phase as they propagate, so that together they form a single pulsating propagating wave; (ii) the waves in each layer are no longer identical, and antiphase pulsations occur, with the two waves alternately advancing and receding as they propagate. This wave is the analog of the spinning wave on the surface of a circular cylinder; (iii) finally, there is a region of bistability between the in phase and antiphase waves. with each having its own domain of attraction, so that which of the two behaviors occur depends on the condition of initiation of the wave. The results of our computations indicate a qualitative similarity in the behavior of combustion waves in the layer model and in the full 2D model. Specifically, due to the similarity between the small alpha wave behavior in the layer model and the large diameter behavior in the model of waves on the surface of a cylinder, we are able to predict the behavior of the mean velocity for the waves on the cylinder, where computations of the full problem can be rather difficult. We also compute solutions for three or more layers. The results of our computations prompt us to predict that, while planar uniformly propagating waves are unstable, the wave will be quasiplanar, i.e., the resulting spinning waves have very low amplitude hot spots, and travel with the velocity close to that of the uniformly propagating wave. Such waves may be difficult to distinguish from uniformly propagating waves in experiments. We also find that for both the layer model and full 2D problem, steady-state time-dependent waves, e.g., pulsating and spinning wave, have a conserved quantity H which characterizes the excess energy in the wave, just as in the case of uniformly propagating waves. The quantity H, which is generated by dissipation, does not vary in time and is proportional to the diffusivity and caloricity of the system, and inversely proportional to the mean wave velocity.

Aldushin, A. P.↗
NASA NTRS1997

Eye Movement Latencies to Direction Change for Different Classes of Motion

In the analysis of visual motion, local features such as orientation are analyzed early in the cortical processing stream (V1), while integration across orientation and space is thought to occur in higher cortical areas such as MT, MST, etc. If all areas provide inputs to eye movement control centers, we would expect that local properties would drive eye movements with relatively short latencies, while global properties would require longer latencies. When such latencies are observed, they can provide information about when (and where?) various stimulus properties are analyzed. To this end, a stimulus was employed in which local and global properties determining perceived direction-of-motion could be manipulated independently: an elliptical Gabor patch with a drifting carrier, with variable orientation of the carrier grating and the contrast window. We have previously demonstrated that the directional percepts evoked by this stimulus vary between the "grating direction" (the normal to the grating's orientation) and the "window direction" (ARVO 91, 94), and that similar effects can be observed in reflexive eye movements (ARVO 95). Subjects viewed such a stimulus while attempting to maintain steady fixation on the center of the pattern, and the small reflexive eye movements ("stare OKN") were recorded. In the middle of the trial, the orientation of either the grating or the window was rotated smoothly by 30 degrees. Responses to the shift of both grating orientation and window orientation are seen in the average OKN slow phase velocity. Grating rotations produce a rapid OKN rotation to the grating direction (100 ms latency, 300 ms time constant), followed by a slower rebound to the steady state perceived direction midway between the grating and window directions. Window rotations, on the other hand, evoke a slower response (200 ms latency, 500 ms time constant). The results demonstrate multiple cortical inputs to eye movement control: a taste early input driven by orientation, and a slower input from higher areas sensitive to global stimulus properties.

Mulligan, Jeffrey B.↗
NASA NTRS2017

Biosignatures in the Context of Low Energy Flux

Many of the features that are thought of as biosignatures - including the mediation of chemical and physical processes with speed, specificity, and selectivity - result directly or indirectly from life's unique capability to mediate and direct energy flux. As such, it is important to consider the impact that differences in energy flux may have on the quantity and quality of evidence for life. Earth differs from every other body in our solar system in the magnitude of biologically-usable energy flux into a liquid water environment. On a global basis, the capture of light energy into photosynthesis and the flux of chemical energy represented in the products of that photosynthesis (organic material + O2) are about six and four orders of magnitude larger, respectively, than the flux of energy represented in geochemical sources. Our conception of what an inhabited world "looks like" and our intuition about how to search for life are based in this high-energy context. Energy fluxes on worlds beyond Earth may be better approximated by the million-fold smaller flux provided to Earth's biosphere by geochemical sources. As a result, the nature, abundance, and quality of evidence for life that could be expected on an inhabited extraterrestrial world within our solar system may differ profoundly from that found on Earth. Understanding this potential difference in quantitative terms provides important context for the formulation of life detection strategies. The influence of energy flux on biosignatures can be evaluated through reference to the two basic purposes into which life partitions energy flux: (1) Life expends energy to sustain existing biomass in a metabolic steady state (metabolically functional but non-growing). The formal representation of this relationship in the traditional microbiology literature equates biomass directly with energy flux. The direct implication is that worlds having lower energy flux will have correspondingly lower potential to support biomass. Life detection strategies that directly target extant organisms should therefore be prepared to encounter average biomass densities that may be many orders of magnitude smaller than those found in most of Earth's surface environments (2) Life expends energy to synthesize new biomass. An end-member case in which new biomass is created at the energy-limited rate and the corresponding cells are immediately destroyed (so that the energy partitioned to cell maintenance is minimized) establishes an upper bound on the rate at which biological material can enter a bulk global pool. For a specified bulk concentration [i] of any particular biological compound, i, or for biologically produced matter overall, this synthesis rate, R (sub i), defines a characteristic time scale tau (sub i) equals [i] divided by R (sub i). tau (sub i) can be thought of as (a) the minimum time required for biosynthesis to yield a specific bulk concentration (e.g., a detection threshold) of i, and (b) the average residence time of i within a bulk pool when [i] is held in steady state through a balance between biosynthesis and attrition by physical, chemical, or biological consumption. tau (sub i) becomes an important quantity in considering the potential utility of enantiomeric excess (as a product of homochiral biosynthesis) as a biosignature. Spontaneous racemization of amino acids acts to "erase" the signature of homochiral synthesis over time scales that may range from hundreds to hundreds of thousands of years, depending on temperature. For environments in which low energy flux translates to low rates of biosynthesis, including the synthesis of homochiral amino acids, amino acid residence times in pools having detectable concentrations may compare to or significantly exceed the time scale for racemization. This and similar consequences of long residence times should be considered in the formulation of life detection strategies based on detection of biologically-produced species. Fluxes of biologically-useful energy on potentially habitable worlds within our solar system are, at present, not well constrained. Improving such constraint has the potential to inform priorities in the formulation and targeting of search-for-life strategies, based on the implications of energy flux for the abundance and quality of biosignatures overall, and in specific categories.

ENERGY FLUX↗
NASA NTRS2014

Landform Evolution Modeling of Sedimentary Processes on Icy Worlds: the Cases of Hyperion and Helene

Sedimentary activity (erosion, transportation and deposition) is a major landscape-shaping process on many icy worlds of the outer solar system for which an understanding of this process is central to their geologic characterizations. Several recent studies, lead by the authors, have identified the effects of sedimentary activity on icy satellite landforms and recognized the role of mass wasting and volatile loss and redistribution in the evolution of these features. We apply state of the art, physics based, landform evolution modeling to icy satellite sedimentary landforms in order to fully quantitatively characterize how their morphologies change with time, and what their present appearances imply for the initial abundance and distribution of loose material, volatiles, and refractories in original pristine landforms. We will present results of our landform evolution modeling to the development of the landscapes of the Saturnian moons of Hyperion and Helene. Our current conclusions for Hyperion are that its unique appearance can be explained in part by the loss to space of ballistic ejecta during impact events. In order to create the smooth surfaces and the reticulate, honeycomb pattern of narrow divides between old craters, appreciable subsequent modification of crater morphology must occur through mass-wasting processes accompanied by sublimation, probably facilitated by the loss of CO2 as a component of the relief-supporting matrix of the bedrock. This mass wasting effectively destroys small craters, at least in part accounting for the paucity of sub-kilometer craters on Hyperion. Helene’s unusual morphology consists of broad depressions (modified large craters) and a generally smooth surface patterned with streaks and grooves. The streaks appear to be oriented down-gradient, as are the grooves. This pattern suggests intensive mass-wasting as a dominant process. Our initial modeling of this surface suggests a Bingham-like rheological behavior for the loose down-slope moving material. Interestingly, as a Bingham flow, the models indicate that aperiodic "intermittent" behavior to be present, suggesting that periods of quiet steady landform evolution are punctuated short periods of active surface readjustment.

Saturnian moons↗
NASA NTRS2001

Fusion for Space Propulsion and Plasma Liner Driven MTF

The need for fusion propulsion for interplanetary flights is discussed. For a propulsion system, there are three important system attributes: (1) The absolute amount of energy available, (2) the propellant exhaust velocity, and (3) the jet power per unit mass of the propulsion system (specific power). For human exploration and development of the solar system, propellant exhaust velocity in excess of 100 km/s and specific power in excess of 10 kW/kg are required. Chemical combustion cannot meet the requirement in propellant exhaust velocity. Nuclear fission processes typically result in producing energy in the form of heat that needs to be manipulated at temperatures limited by materials to about 2,800 K. Using the energy to heat a low atomic weight propellant cannot overcome the problem. Alternatively the energy can be converted into electricity which is then used to accelerate particles to high exhaust velocity. The necessary power conversion and conditioning equipment, however, increases the mass of the propulsion system for the same jet power by more than two orders of magnitude over chemical system, thus greatly limits the thrust-to-weight ratio attainable. If fusion can be developed, fusion appears to have the best of all worlds in terms of propulsion - it can provide the absolute amount, the propellant exhaust velocity, and the high specific jet power. An intermediate step towards pure fusion propulsion is a bimodal system in which a fission reactor is used to provide some of the energy to drive a fusion propulsion unit. The technical issues related to fusion for space propulsion are discussed. There are similarities as well as differences at the system level between applying fusion to propulsion and to terrestrial electrical power generation. The differences potentially provide a wider window of opportunities for applying fusion to propulsion. For example, pulsed approaches to fusion may be attractive for the propulsion application. This is particularly so in the light of significant development of the enabling pulsed power component technologies that have occurred in the last two decades because of defense and other energy requirements. The extreme states of matter required to produce fusion reactions may be more readily realizable in the pulsed states with less system mass than in steady states. Significant saving in system mass may result in pulsed fusion systems using plasmas in the appropriate density regimes. Magnetized target fusion, which attempts to combine the favorable attributes of magnetic confinement and inertial compression-containment into one single integrated fusion scheme, appears to have benefits that are worth exploring for propulsion application.

Thio, Y.C. Francis↗
NASA NTRS2007

Nanorod-Based Fast-Response Pressure-Sensitive Paints

A proposed program of research and development would be devoted to exploitation of nanomaterials in pressuresensitive paints (PSPs), which are used on wind-tunnel models for mapping surface pressures associated with flow fields. Heretofore, some success has been achieved in measuring steady-state pressures by use of PSPs, but success in measuring temporally varying pressures has been elusive because of the inherent slowness of the optical responses of these materials. A PSP contains a dye that luminesces in a suitable wavelength range in response to photoexcitation in a shorter wavelength range. The luminescence is quenched by oxygen at a rate proportional to the partial pressure of oxygen and thus proportional to the pressure of air. As a result, the intensity of luminescence varies inversely with the pressure of air. The major problem in developing a PSP that could be easily applied to a wind-tunnel model and could be useful for measuring rapidly varying pressure is to provide very high gas diffusivity for rapid, easy transport of oxygen to and from active dye molecules. Most PSPs include polymer-base binders, which limit the penetration of oxygen to dye molecules, thereby reducing responses to pressure fluctuations. The proposed incorporation of nanomaterials (somewhat more specifically, nanorods) would result in paints having nanostructured surfaces that, relative to conventional PSP surfaces, would afford easier and more nearly complete access of oxygen molecules to dye molecules. One measure of greater access is effective surface area: For a typical PSP as proposed applied to a given solid surface, the nanometer-scale structural features would result in an exposed surface area more than 100 times that of a conventional PSP, and the mass of proposed PSP needed to cover the surface would be less than tenth of the mass of the conventional PSP. One aspect of the proposed development would be to synthesize nanorods of Si/SiO2, in both tangle-mat and regular- array forms, by use of chemical vapor deposition (CVD) and wet chemical processes, respectively. The rods would be coated with a PSP dye, and the resulting PSP signals would be compared with those obtained from PSP dye coats on conventional support materials. Another aspect of the proposed development would be to seek to exploit the quantum properties of nanorods of a suitable semiconductor (possibly GaN), which would be synthesized by CVD. These quantum properties of semiconductor nanorods include narrow-wavelength-band optical absorption and emission characteristics that vary with temperature. The temperature sensitivity might enable simultaneous measurement of fluctuating temperature and pressure and to provide a temperature correction for the PSP response.

Bencic, Timothy↗
NASA NTRS2007

Hybrid Automated Diagnosis of Discrete/Continuous Systems

A recently conceived method of automated diagnosis of a complex electromechanical system affords a complete set of capabilities for hybrid diagnosis in the case in which the state of the electromechanical system is characterized by both continuous and discrete values (as represented by analog and digital signals, respectively). The method is an integration of two complementary diagnostic systems: (1) beacon-based exception analysis for multi-missions (BEAM), which is primarily useful in the continuous domain and easily performs diagnoses in the presence of transients; and (2) Livingstone, which is primarily useful in the discrete domain and is typically restricted to quasi-steady conditions. BEAM has been described in several prior NASA Tech Briefs articles: "Software for Autonomous Diagnosis of Complex Systems" (NPO-20803), Vol. 26, No. 3 (March 2002), page 33; "Beacon-Based Exception Analysis for Multimissions" (NPO-20827), Vol. 26, No. 9 (September 2002), page 32; "Wavelet-Based Real-Time Diagnosis of Complex Systems" (NPO-20830), Vol. 27, No. 1 (January 2003), page 67; and "Integrated Formulation of Beacon-Based Exception Analysis for Multimissions" (NPO-21126), Vol. 27, No. 3 (March 2003), page 74. Briefly, BEAM is a complete data-analysis method, implemented in software, for real-time or off-line detection and characterization of faults. The basic premise of BEAM is to characterize a system from all available observations and train the characterization with respect to normal phases of operation. The observations are primarily continuous in nature. BEAM isolates anomalies by analyzing the deviations from nominal for each phase of operation. Livingstone is a model-based reasoner that uses a model of a system, controller commands, and sensor observations to track the system s state, and detect and diagnose faults. Livingstone models a system within the discrete domain. Therefore, continuous sensor readings, as well as time, must be discretized. To reason about continuous systems, Livingstone uses monitors that discretize the sensor readings using trending and thresholding techniques. In development of the a hybrid method, BEAM results were sent to Livingstone to serve as an independent source of evidence that is in addition to the evidence gathered by Livingstone standard monitors. The figure depicts the flow of data in an early version of a hybrid system dedicated to diagnosing a simulated electromechanical system. In effect, BEAM served as a "smart" monitor for Livingstone. BEAM read the simulation data, processed the data to form observations, and stored the observations in a file. A monitor stub synchronized the events recorded by BEAM with the output of the Livingstone standard monitors according to time tags. This information was fed to a real-time interface, which buffered and fed the information to Livingstone, and requested diagnoses at the appropriate times. In a test, the hybrid system was found to correctly identify a failed component in an electromechanical system for which neither BEAM nor Livingstone alone yielded the correct diagnosis.

Park, Han↗
NASA NTRS2004

Increasing Accuracy in Computed Inviscid Boundary Conditions

A technique has been devised to increase the accuracy of computational simulations of flows of inviscid fluids by increasing the accuracy with which surface boundary conditions are represented. This technique is expected to be especially beneficial for computational aeroacoustics, wherein it enables proper accounting, not only for acoustic waves, but also for vorticity and entropy waves, at surfaces. Heretofore, inviscid nonlinear surface boundary conditions have been limited to third-order accuracy in time for stationary surfaces and to first-order accuracy in time for moving surfaces. For steady-state calculations, it may be possible to achieve higher accuracy in space, but high accuracy in time is needed for efficient simulation of multiscale unsteady flow phenomena. The present technique is the first surface treatment that provides the needed high accuracy through proper accounting of higher-order time derivatives. The present technique is founded on a method known in art as the Hermitian modified solution approximation (MESA) scheme. This is because high time accuracy at a surface depends upon, among other things, correction of the spatial cross-derivatives of flow variables, and many of these cross-derivatives are included explicitly on the computational grid in the MESA scheme. (Alternatively, a related method other than the MESA scheme could be used, as long as the method involves consistent application of the effects of the cross-derivatives.) While the mathematical derivation of the present technique is too lengthy and complex to fit within the space available for this article, the technique itself can be characterized in relatively simple terms: The technique involves correction of surface-normal spatial pressure derivatives at a boundary surface to satisfy the governing equations and the boundary conditions and thereby achieve arbitrarily high orders of time accuracy in special cases. The boundary conditions can now include a potentially infinite number of time derivatives of surface-normal velocity (consistent with no flow through the boundary) up to arbitrarily high order. The corrections for the first-order spatial derivatives of pressure are calculated by use of the first-order time derivative velocity. The corrected first-order spatial derivatives are used to calculate the second- order time derivatives of velocity, which, in turn, are used to calculate the corrections for the second-order pressure derivatives. The process as described is repeated, progressing through increasing orders of derivatives, until the desired accuracy is attained.

Dyson, Roger↗
NASA NTRS2010

Adaptive Numerical Algorithms in Space Weather Modeling

Space weather describes the various processes in the Sun-Earth system that present danger to human health and technology. The goal of space weather forecasting is to provide an opportunity to mitigate these negative effects. Physics-based space weather modeling is characterized by disparate temporal and spatial scales as well as by different physics in different domains. A multi-physics system can be modeled by a software framework comprising of several components. Each component corresponds to a physics domain, and each component is represented by one or more numerical models. The publicly available Space Weather Modeling Framework (SWMF) can execute and couple together several components distributed over a parallel machine in a flexible and efficient manner. The framework also allows resolving disparate spatial and temporal scales with independent spatial and temporal discretizations in the various models. Several of the computationally most expensive domains of the framework are modeled by the Block-Adaptive Tree Solar wind Roe Upwind Scheme (BATS-R-US) code that can solve various forms of the magnetohydrodynamics (MHD) equations, including Hall, semi-relativistic, multi-species and multi-fluid MHD, anisotropic pressure, radiative transport and heat conduction. Modeling disparate scales within BATS-R-US is achieved by a block-adaptive mesh both in Cartesian and generalized coordinates. Most recently we have created a new core for BATS-R-US: the Block-Adaptive Tree Library (BATL) that provides a general toolkit for creating, load balancing and message passing in a 1, 2 or 3 dimensional block-adaptive grid. We describe the algorithms of BATL and demonstrate its efficiency and scaling properties for various problems. BATS-R-US uses several time-integration schemes to address multiple time-scales: explicit time stepping with fixed or local time steps, partially steady-state evolution, point-implicit, semi-implicit, explicit/implicit, and fully implicit numerical schemes. Depending on the application, we find that different time stepping methods are optimal. Several of the time integration schemes exploit the block-based granularity of the grid structure. The framework and the adaptive algorithms enable physics based space weather modeling and even forecasting.

Toth, Gabor↗
NASA NTRS1998

Cardiovascular function following reduced aerobic activity

PURPOSE: The aim of this study was to test the hypothesis that a sustained reduction of physical activity (deconditioning) would alter the cardiovascular regulatory function. METHODS: Nineteen young, healthy volunteers participated in physical deconditioning for a period of 8 wk. Before (pre) and following (post) physical deconditioning, the responses of heart rate (HR), mean arterial pressure (MAP, measured by Finapres), central venous pressure (CVP), stroke volume (SV, Doppler), and forearm blood flow (FBF, plethysmography) were determined during lower body negative pressure (LBNP). The carotid baroreflex (CBR) function was assessed using a train of pulsatile neck pressure (NP) and suction, and the aortic baroreflex control of HR was assessed during steady-state phenylephrine (PE) infusion superimposed by LBNP and NP to counteract the PE increased CVP and carotid sinus pressure, respectively. RESULTS: Active physical deconditioning significantly decreased maximal oxygen uptake (-7%) and LBNP tolerance (-13%) without a change in baseline hemodynamics. Plasma volume (-3% at P = 0.135), determined by Evans Blue dilution, and blood volume (-4% at P = 0.107) were not significantly altered. During LBNP -20 to -50 torr, there was a significantly greater drop of SV per unit decrease in CVP in the post- (14.7 +/- 1.6%/mm Hg) than predeconditioning (11.2 +/- 0.7%/mm Hg) test accompanied by a greater tachycardia. Deconditioning increased the aortic baroreflex sensitivity (pre vs post: -0.61 +/- 0.12 vs -0.84 +/- 0.14 bpm.mm-1 Hg, P = 0.009) and the slope of forearm vascular resistance (calculated from [MAP-CVP]/FBF) to CVP (-2.75 +/- 0.26 vs -4.94 +/- 0.97 PRU/mm Hg, P = 0.086). However, neither the CBR-HR (-0.28 +/- 0.03 VS -0.39 +/- 0.10 bpm.mm-1 Hg) nor the CBR-MAP (-0.37 +/- 0.16 vs -0.25 +/- 0.07 mm Hg/mm Hg) gains were statistically different between pre- and postdeconditioning. CONCLUSIONS: We concluded that the functional modification of the cardiac pressure-volume relationship resulted in the reduced LBNP tolerance, despite the accentuated aortic and cardiopulmonary baroreflex function following deconditioning.

Non-NASA Center↗
NASA NTRS2012

Simulations of Magnetic Reconnection - Kinetic Mechanisms Underlying the Fluid Description of Ions

Because of its ability to transfer the energy stored in magnetic field together with the breaking of the flux freezing constraint, magnetic reconnection is considered as one of the most important phenomena in plasma physics. When it happens in a collision less environment such as the terrestrial magnetosphere, it should a priori be modelled with in the framework of kinetic physics. The evidence of kinetic features has incidentally for a long time, been shown by researchers with the help of both numerical simulations and satellite observations. However, most of our understanding of the process comes from the more intuitive fluid interpretation with simple closure hypothesis which do not include kinetic effects. To what extent are these two separate descriptions of the same phenomenon related? What is the role of kinetic effects in the averaged/fluid dynamics of reconnection? This thesis addresses these questions for the proton population in the particular case of anti parallel merging with the help of 2D Hybrid simulations. We show that one can not assume, as is usually done, that the acceleration of the proton flow is only due to the Laplace force. Our results show, for symmetric and asymmetric connection, the importance of the pressure force, opposed to the electric one on the separatrices, in the decoupling region. In the symmetric case, we emphasize the kinetic origin of this force by analyzing the proton distribution functions and explain their structure by studying the underlying particle dynamics. Protons, as individual particles, are shown to bounce in the electric potential well created by the Hall effect. The spatial divergence of this well results in a mixing in phase space responsible for the observed structure of the pressure tensor. A detailed energy budget analysis confirms the role of the pressure force for the acceleration; but, contrary to what is sometimes assumed, it also reveals that the major part of the incoming Poynting flux is transferred to the thermal energy flux rather than to the convective kinetic energy flux, although the latter is generally supposed dominant. In the symmetric case, we propose the pressure tensor to be an additional proxy of the ion decoupling region in satellite data and verify this suggestion by studying a reconnection event encountered by the Cluster spacecrafts. Finally, the last part of this thesis is devoted to the study of the kinetic structure of asymmetric tangential current sheets where connection can develop. This theoretical part consists in finding a steady state solution to the Vlasov-Maxwell system for the protons in such a configuration. We present the theory and its first confrontation to numerical tests.

Aunai, icolas↗
NASA NTRS2005

Formation of a Tropopause Cirrus Layer Observed over Florida during CRYSTAL-FACE

On July 13, 2002 a widespread, subvisible tropopause cirrus layer occurred over the Florida region. This cloud was observed in great detail with the NASA Cirrus Regional Study of Tropical Anvils and Cirrus Layers-Florida Area Cirrus Experiment (CRYSTAL-FACE) instrumentation, including in situ measurements with the WB-57 aircraft. In this paper, we use the 13 July cloud as a case study to evaluate the physical processes controlling the formation and evolution of tropopause cirrus layers. Microphysics measurements indicate that ice crystal diameters in the cloud layer ranged from about 7 to 50 microns, and the peak number mode was about 10-25 microns. In situ water vapor and temperature measurements in the cloud indicated supersaturation with respect to ice throughout, with ice saturation ratios as large as 1.8. Even when the ice surface area density was as high as about 500 sq microns/cu cm, ice supersaturations of 20-30% were observed. Trajectory analysis shows that the air sampled near the tropopause on this day generally came from the north and cooled considerably during the previous few days. Examination of infrared satellite imagery along air parcel back trajectories from the WB-57 flight track indicates that the tropopause cloud layer formation was, in general, not simply left over ice from recently generated anvil cirrus. Simulations of cloud formation using time-height curtains of temperature along the trajectory paths show that the cloud could have formed in situ near the tropopause as the air was advected into the south Florida region and cooled to unusually low temperatures. If we assume a high threshold for ice nucleation via homogeneous freezing of aqueous sulfate aerosols, the model reproduces the observed cloud structure, ice crystal size distributions, and ice supersaturation statistics. Inclusion of observed gravity wave temperature perturbations in the simulations is essential to reproduce the observed cloud properties. Without waves, crystal number densities are too low, crystal sizes are too large, and the crystals fall out too fast, leaving very little cloud persisting at the end of the simulations. In the cloud simulations, coincidence of high supersaturations and high surface areas can be produced by either recent nucleation or sedimentation of crystals into supersaturated layers. The agreement between model results and observed supersaturations is improved somewhat if we assume that the steady state relative humidity within cirrus at T<200 K is enhanced by about 30%. The WB-57 measurements and the model results suggest that the cloud layer irreversibly dehydrated air near the tropopause.

Jensen, Eric↗
NASA NTRSDate not supplied

Effect of Electronic Shot Noise on Dynamic Measurements using Optical Techniques: Examples from Rayleigh Scattering and Unsteady PSP

Electronic shot noise is an unavoidable reality in all optical techniques that depend on measuring light intensity. For steady-state, time-averaged measurements the impact of shot noise can be easily reduced by increasing the exposure time, or by averaging over multiple-exposures. That is not the case for unsteady measurements, where time histories of light intensity variations need to be created either by high-speed photography (unsteady PSP) or via photoelectron counting (spectrally resolved Rayleigh scattering) over short time durations using a photo-multiplier tube (PMT) and photon counting electronics. Electronic shot noise introduces a fixed amount of random error, which can overwhelm the light intensity variation caused by turbulent fluctuations. Spectrum computed from such time series shows a fixed noise floor that is independent of the number of data points in the time series. For a fixed optical system, where the collected power of the luminescent light (uPSP), or the scattered light (Rayleigh) is fixed, one needs to resort to special techniques to obtain enough signal-to noise ratio (SNR). For the uPSP application it is shown that averaging of the adjacent pixels improves SNR; although this may lead to a sacrifice of spatial resolution. A second means is to increase the exposure time, which leads to a lowering of the frequency range. For the Rayleigh application, improvements of SNR can be achieved via two different cross-correlation based approaches. The first involves measuring the light intensity using two PMTs using short, contiguous gates; thereby, creating two time-series of data. The second one involves collecting one long time-series of data using one set of measurement device, followed by an odd-even splitting into two time series. When the two time-series are cross-correlated, and a power spectrum is calculated, a significant reduction in the shot noise floor can be achieved. Examples from measurements of density and velocity fluctuations spectra from two different Rayleigh setup are presented to demonstrate the process.

uPSP↗
NASA NTRS1998

Scaling Properties of Particle Density Fields Formed in Simulated Turbulent Flows

Direct numerical simulations (DNS) of particle concentrations in fully developed 3D turbulence were carried out in order to study the nonuniform structure of the particle density field. Three steady-state turbulent fluid fields with Taylor microscale Reynolds numbers (Re(sub lambda)) of 40, 80 and 140 were generated by solving the Navier-Stokes equations with pseudospectral methods. Large scale forcing was used to drive the turbulence and maintain temporal stationarity. The response of the particles to the fluid was parameterized by the particle Stokes number St, defined as the ratio of the particle's stopping time to the mean period of eddies on the Kolmogorov scale (eta). In this paper, we consider only passive particles optimally coupled to these eddies (St approx. = 1) because of their tendency to concentrate more than particles with lesser or greater St values. The trajectories of up to 70 million particles were tracked in the equilibrated turbulent flows until the particle concentration field reached a statistically stationary state. The nonuniform structure of the concentration fields was characterized by the multifractal singularity spectrum, f(alpha), derived from measures obtained after binning particles into cells ranging from 2(eta) to 15(eta) in size. We observed strong systematic variations of f(alpha) across this scale range in all three simulations and conclude that the particle concentration field is not statistically self similar across the scale range explored. However, spectra obtained at the 2(eta), 4(eta), and 8(eta) scales of each flow case were found to be qualitatively similar. This result suggests that the local structure of the particle concentration field may be flow-Independent. The singularity spectra found for 2n-sized cells were used to predict concentration distributions in good agreement with those obtained directly from the particle data. This Singularity spectrum has a shape similar to the analogous spectrum derived for the inertial-range energy dissipation fields of experimental turbulent flows at Re(sub lambda) = 110 and 1100. Based on this agreement, and the expectation that both dissipation and particle concentration are controlled by the same cascade process, we hypothesize that singularity spectra similar to the ones found in this work provide a good characterization of the spatially averaged statistical properties of preferentially concentrated particles in higher Re(sub lambda) turbulent flows.

Hogan, Robert C.↗
NASA NTRS1994

Magnetocentrifugally driven flows from young stars and disks. 1: A generalized model

We propose a generalized model for stellar spin-down, disk accretion, and truncation, and the origin of winds, jets, and bipolar outflows from young stellar objects. We consider the steady state dynamics of accretion of matter from a viscous and imperfectly conducting disk onto a young star with a strong magnetic field. For an aligned stellar magnetosphere, shielding currents in the surface layers of the disk prevent stellar field lines from penetrating the disk everywhere except for a range of radii about pi = R(sub x), where the Keplerian angular speed of rotation Omega(sub x) equals the angular speed of the star Omega(sub *). For the low disk accretion rates and high magnetic fields associated with typical T Tauri stars, R(sub x) exceeds the radius of the star R(sub *) by a factor of a few, and the inner disk is effectively truncated at a radius R(sub t) somewhat smaller than R(sub x). Where the closed field lines between R(sub t) and R(sub x) bow sufficiently inward, the accreting gas attaches itself to the field and is funneled dynamically down the effective potential (gravitational plus centrifugal) onto the star. Contrary to common belief, the accompanying magnetic torques associated with this accreting gas may transfer angular momentum mostly to the disk rather than to the star. Thus, the star can spin slowly as long as R(sub x) remains significantly greater than R(sub *). Exterior to R(sub x) field lines threading the disk bow outward, which makes the gas off the mid-plane rotate at super-Keplerian velocities. This combination drives a magnetocentrifugal wind with a mass-loss rate M(sub w) equal to a definite fraction f of the disk accretion rate M(sub D). For high disk accretion rates, R(sub x) is forced down to the stellar surface, the star is spun to breakup, and the wind is generated in a manner identical to that proposed by Shu, Lizano, Ruden, & Najita in a previous communication to this journal. In two companion papers (II and III), we develop a detailed but idealized theory of the magnetocentrifugal acceleration process.

Shu, Frank↗
NASA NTRS1995

Human parathyroid hormone-(1-38) restores cancellous bone to the immobilized, osteopenic proximal tibial metaphysis in rats

The purpose of this study was to determine if human parathyroid hormone-(1-38) (hPTH(1-38)) can restore cancellous bone mass to the established osteopenic, immobilized proximal tibial metaphyses of female rats. The right hindlimbs of 6-month-old female Sprague-Dawley rats were immobilized by bandaging the right hindlimbs to the abdomen. After 30 days of right hindlimb immobilization, the rats were subcutaneously injected with 200 micrograms hPTH(1-38)/kg/day for 15 days (short-term treatment) or 75 days (longer-term treatment). Static bone histomorphometry was performed on the primary spongiosa, and both static and dynamic histomorphometry were performed on the secondary spongiosa of the right proximal tibial metaphyses. Immobilization for 30 days without treatment decreased trabecular bone area, number, and thickness in both primary and secondary spongiosa, and induced an increase in eroded perimeter and a decrease in tissue referent-bone formation rate in the secondary spongiosa. These changes reached a new steady state thereafter. Treatment with 200 micrograms hPTH(1-38)/kg/day for 15 days, beginning 30 days after immobilization, significantly increased trabecular bone area, thickness, and number in both primary and secondary spongiosa despite continuous immobilization when compared with controls. The short-term PTH treatment (15 days) significantly increased labeling perimeter, mineral apposition rate, and tissue referent-bone formation rate in the secondary spongiosa and stimulated longitudinal bone growth as compared with the controls. Longer PTH treatment (75 days) further increased trabecular bone area, thickness, and number as compared with controls and groups given short-term PTH treatment (15 days). The bone formation indices in the secondary spongiosa of the longer-term treated rats were lower than those of the short-term treated group, but they were still higher than those of controls. Our findings indicate that PTH treatment stimulates cancellous bone formation, and restores and adds extra cancellous bone to the established, disuse-osteopenic proximal tibial metaphysis of female rats with continuously immobilized right hindlimbs. These results suggest that PTH may be useful in treating disuse-induced osteoporosis in humans.

NASA Program Space Physiology and Countermeasures↗
NASA NTRS2001

Runway Texture and Grid Pattern Effects on Rate-of-Descent Perception

To date, perceptual errors occur in determining descent rate from a computer-generated image in flight simulation. Pilots tend to touch down twice as hard in simulation than in flight, and more training time is needed in simulation before reaching steady-state performance. Barnes suggested that recognition of range may be the culprit, and he cited that problems such as collimated objects, binocular vision, and poor resolution lead to poor estimation of the velocity vector. Brown's study essentially ruled out that the lack of binocular vision is the problem. Dorfel added specificity to the problem by showing that pilots underestimated range in simulated scenes by 50% when 800 ft from the runway threshold. Palmer and Petitt showed that pilots are able to distinguish between a 1.7 ft/sec and 2.9 ft/sec sink rate when passively observing sink rates in a night scene. Platform motion also plays a role, as previous research has shown that the addition of substantial platform motion improves pilot estimates of vertical velocity and results in simulated touchdown rates more closely resembling flight. This experiment examined how some specific variations in the visual scene properties affect a pilot's perception of sink rate. It extended another experiment that focused on the visual and motion cues necessary for helicopter autorotations. In that experiment, pilots performed steep approaches to a runway. The visual content of the runway and its surroundings varied in two ways: texture and rectangular grid spacing. Four textures, included a no-texture case, were evaluated. Three grid spacings, including a no-grid case, were evaluated. The results showed that pilot better controlled their vertical descent rates when good texture cues were present. No significant differences were found for the grid manipulation. Using those visual scenes a simple psychophysics, experiment was performed. The purpose was to determine if the variations in the visual scenes allowed pilots to better perceive vertical velocity. To determine that answer, pilots passively viewed a particular visual scene in which the vehicle was descending at two different rates. Pilots had to select which of the two rates they thought was the fastest rate. The difference between the two rates changed using a staircase method, depending on whether or not the pilot was correct, until a minimum threshold between the two descent rates was reached. This process was repeated for all of the visual scenes to decide whether or not the visual scenes did allow pilots to perceive vertical velocity better among them. All of the data have yet to be analyzed; however, neither the effects of grid nor texture revealed any statistically significant trends. On further examination of the staircase method employed, a possibility exists that the lack of an evident trend may be due to the exit criterion used during the study. As such, the experiment will be repeated with an improved exit criterion in February. Results of this study will be presented in the submitted paper.

Schroeder, J. A.↗
NASA NTRS2006

The Solar-B Mission: First Light, Future Plans and Community Participation

The Solar-B spacecraft was launched from the Uchinoura Space Center into a circular, sun-synchronous, polar orbit by the Japanese Aerospace Exploration Agency in late September 2006. The spacecraft carries thee scientific instruments designed to follow the flow of magnetic energy from the photosphere to the corona to improve our understanding of both steady state and transient energy release. This goal will be achieved through coordinated observations of three highly advanced solar telescopes developed cooperatively by teams from Japan, the United States and the United Kingdom. The three telescopes are a 0.5m aperture, diffraction limited, solar optical telescope (SOT), an X-ray telescope (XRT) designed for full sun imaging with 1.0 arcsec pixels and an EUV imaging spectrometer (EIS) with an order of magnitude improvement in sensitivity over past instruments. The SOT focal plane contains three instruments, a spectropolarimeter for measuring vector magnetic fields, a broadband filter imager for recording images of the photosphere and chromosphere at the highest resolution the telescope is capable of, and a narrow band filter imager that will record Doppler grams and vector magnetograms. The XRT has broad temperature coverage and a spatial a resolution three times as high as Yohkoh. EIS covers a broad range of transition region and coronal temperatures in two spectral bands. Both XRT and EIS have 2 arcsec spatial resolution (1 arcsec pixels). Instrument first light occurred after five weeks on orbit to allow for out gassing and the opening of the telescopes doors. The initial observation sequences are designed to test the functionality of the different operating modes and for calibration. After this commissioning phase is complete a series of observations are planned to demonstrate the ability of the instruments to meet NASA's mission minimum success criteria. Data is downloaded every orbit to the Norwegian high latitude ground station at Svalbard. The data are transmitted to ISAS where they are reformatted into FITS files and archived as Level 0 data on the ISAS DARTS system. Once the initial observation period is complete, approximately six months after launch, the mission data will be open and freely available to researchers shortly after receipt at the DARTS data archive hosted in Japan and at NASA s Solar Data Analysis Center at the Goddard Space Flight Center. Scientific operations will be conducted from the ISAS facility in Sagamihara, Japan and the observatory will become available for performing joint operations with both ground and space based instruments and for conducting observing programs proposed by non-team members. This process will be described together with a status report from the initial operation of the observatory, showing examples of the first observations.

Davis, John M.↗