Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “1D”

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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
NASA NTRS1995

The middle atmospheric response to short and long term solar UV variations: analysis of observations and 2D model results

We have investigated the middle atmospheric response to the 27-day and 11-yr solar UV flux variations at low to middle latitudes using a two-dimensional photochemical model. The model reproduced most features of the observed 27-day sensitivity and phase lag of the profile ozone response in the upper stratosphere and lower mesosphere, with a maximum sensitivity of +0.51% per 1% change in 205 nm flux. The model also reproduced the observed transition to a negative phase lag above 2 mb, reflecting the increasing importance with height of the solar modulated HO(x) chemistry on the ozone response above 45 km. The model revealed the general anti-correlation of ozone and solar UV at 65-75 km, and simulated strong UV responses of water vapor and HO(x) species in the mesosphere. Consistent with previous 1D model studies, the observed upper mesospheric positive ozone response averaged over +/- 40 was simulated only when the model water vapor concentrations above 75 km were significantly reduced relative to current observations. In agreement with observations, the model computed a low to middle latitude total ozone phase lag of +3 days and a sensitivity of +0.077% per 1% change in 205 nm flux for the 27-day solar variation, and a total ozone sensitivity of +0.27% for the 11-yr solar cycle. This factor of 3 sensitivity difference is indicative of the photochemical time constant for ozone in the lower stratosphere which is comparable to the 27-day solar rotation period but is much shorter than the 11-yr solar cycle.

Fleming, Eric L.↗
NASA NTRS1995

The reaction O((3)P) + HOBr: Temperature dependence of the rate constant and importance of the reaction as an HOBr stratospheric loss process

The absolute rate constant for the reaction O((3)P) + HOBr has been measured between T = 233K and 423K using the discharge-flow kinetic technique coupled to mass spectrometric detection. The value of the rate coefficient at room temperature is (2.5 +/- 0.6) x 10(exp -11)cu cm/molecule/s and the derived Arrhenius expression is (1.4 +/- 0.5) x 10(exp -10) exp((-430 +/- 260)/T)cu cm/molecule/s. From these rate data the atmospheric lifetime of HOBr with respect to reaction with O((3)P) is about 0.6h at z = 25 km which is comparable to the photolysis lifetime based on recent measurements of the UV cross section for HOBr. Implications for HOBr loss in the stratosphere have been tested using a 1D photochemical box model. With the inclusion of the rate parameters and products for the O + HOBr reaction, calculated concentration profiles of BrO increase by up to 33% around z = 35 km. This result indicates that the inclusion of the O + HOBr reaction in global atmospheric chemistry models may have an impact on bromine partitioning in the middle atmosphere.

Nesbitt, F. L.↗
NASA NTRS1995

Efficient implementation of weighted ENO schemes

In this paper, we further analyze, test, modify and improve the high order WENO (weighted essentially non-oscillatory) finite difference schemes of Liu, Osher and Chan. It was shown by Liu et al. that WENO schemes constructed from the r-th order (in L1 norm) ENO schemes are (r+1)-th order accurate. We propose a new way of measuring the smoothness of a numerical solution, emulating the idea of minimizing the total variation of the approximation, which results in a 5-th order WENO scheme for the case r = 3, instead of the 4-th order with the original smoothness measurement by Liu et al. This 5-th order WENO scheme is as fast as the 4-th order WENO scheme of Liu et al., and both schemes are about twice as fast as the 4-th order ENO schemes on vector supercomputers and as fast on serial and parallel computers. For Euler systems of gas dynamics, we suggest computing the weights from pressure and entropy instead of the characteristic values to simplify the costly characteristic procedure. The resulting WENO schemes are about twice as fast as the WENO schemes using the characteristic decompositions to compute weights, and work well for problems which do not contain strong shocks or strong reflected waves. We also prove that, for conservation laws with smooth solutions, all WENO schemes are convergent. Many numerical tests, including the 1D steady state nozzle flow problem and 2D shock entropy wave interaction problem, are presented to demonstrate the remarkable capability of the WENO schemes, especially the WENO scheme using the new smoothness measurement, in resolving complicated shock and flow structures. We have also applied Yang's artificial compression method to the WENO schemes to sharpen contact discontinuities.

Jiang, Guang-Shan↗
NASA NTRS1995

High efficiency solar cells for laser power beaming applications

Understanding solar cell response to pulsed laser outputs is important for the evaluation of power beaming applications. The time response of high efficiency GaAs and silicon solar cells to a 25 nS monochromatic pulse input is described. The PC-1D computer code is used to analyze the cell current during and after the pulse for various conditions.

Jain, Raj K.↗
NASA NTRS1995

A dynamic model of the three-fluid solar wind: Equilibrium structures and fluxes

We report results from a dynamic approach to equilibrium for a three fluid solar-wind model in spherical 1D. We advance the mass, momentum and energy transport equations for protons, alpha particles and electrons, with full interspecies coulomb coupling and the simplifications recommended by Burgi (1992). They include no ion heat conduction and an ambipolar electric field in place of the negligible electron momentum coupling. At the inner boundary (the solar surface) we specify the inflow velocity, mass flux, and temperature for each species, and at the outer boundary at large radius we allow the fluids to flow freely from the computational domain. Starting from an initial wind model (e.g.. single fluid) we evolve to a three-fluid steady state. We will describe comparisons with Burgi's static equilibrium results. and examine the stability of the discontinuities which appear in his derived parameter variations.

Ruden, S. P.↗
NASA NTRS1995

A Large-Particle Monte Carlo Code for Simulating Non-Linear High-Energy Processes Near Compact Objects

High-energy radiation processes in compact cosmic objects are often expected to have a strongly non-linear behavior. Such behavior is shown, for example, by electron-positron pair cascades and the time evolution of relativistic proton distributions in dense radiation fields. Three independent techniques have been developed to simulate these non-linear problems: the kinetic equation approach; the phase-space density (PSD) Monte Carlo method; and the large-particle (LP) Monte Carlo method. In this paper, we present the latest version of the LP method and compare it with the other methods. The efficiency of the method in treating geometrically complex problems is illustrated by showing results of simulations of 1D, 2D and 3D systems. The method is shown to be powerful enough to treat non-spherical geometries, including such effects as bulk motion of the background plasma, reflection of radiation from cold matter, and anisotropic distributions of radiating particles. It can therefore be applied to simulate high-energy processes in such astrophysical systems as accretion discs with coronae, relativistic jets, pulsar magnetospheres and gamma-ray bursts.

Stern, Boris E.↗
NASA NTRS1997

Temperature Dependence of the O + HO2 Rate Coefficient

A pulsed laser photolysis technique has been employed to investigate the kinetics of the radical-radical reaction O((sup 3)P) + HO2 OH + O2 over the temperature range 266-391 K in 80 Torr of N2 diluent gas. O((sup 3)P) was produced by 248.5-nm KrF laser photolysis of O3 followed by rapid quenching of O(1D) to O((sup 3)P) while HO2 was produced by simultaneous photolysis of H2O2 to create OH radicals which, in turn, reacted with H2O2 to yield HO2. The O((sup 3)P) temporal profile was monitored by using time-resolved resonance fluorescence spectroscopy. The HO2 concentration was calculated based on experimentally measured parameters. The following Arrhenius expression describes our experimental results: k(sub 1)(T) equals (2.91 +/- 0.70) x 10(exp -11) exp[(228 +/- 75)/T] where the errors are 2 sigma and represent precision only. The absolute uncertainty in k, at any temperature within the range 266-391 K is estimated to be +/- 22 percent. Our results are in excellent agreement with a discharge flow study of the temperature dependence of k(sub 1) in 1 Torr of He diluent reported by Keyser, and significantly reduce the uncertainty in the rate of this important stratospheric reaction at subambient temperatures.

Nicovich, J. M.↗
NASA NTRS1994

Stabilization Mechanisms and Burning Rates of Cylindrical Burner Flames

A study is conducted of the structure and response of curved (but unstretched), cylindrically-symmetric 1D premixed flames from a cylindrical porous burner. The study has employed (1) activation-energy asymptotics with one-step reaction constant and constant properties; (2) a numerical computation which encompassed detailed chemistry and transport behavior, and (3) drop-tower microgravity tests. Attention was given to the relative importance of heat loss vs. flow divergence as the dominant mechanism for flame stabilization; the results show that, with increasing flow discharge rate, the dominant flame stabilization mechanism changes from heat loss to flow divergence.

Eng, J. A.↗
NASA NTRS1998

The Wavelet Element Method: Realization and Additional Features in 2D and 3D - Part 2

The Wavelet Element Method (WEM) provides a construction of multiresolution systems and biorthogonal wavelets on fairly general domains. These are split into subdomains that are mapped to a single reference hypercube. Tensor products of scaling functions and wavelets defined on the unit interval are used on the reference domain. By introducing appropriate matching conditions across the interelement boundaries, a globally continuous biorthogonal wavelet basis on the general domain is obtained. This construction does not uniquely define the basis functions but rather leaves some freedom for fulfilling additional features. In this paper we detail the general construction principle of the WEM to the 1D, 2D and 3D cases. We address additional features such as symmetry, vanishing moments and minimal support of the wavelet functions in each particular dimension. The construction is illustrated by using biorthogonal spline wavelets on the interval.

Canuto, Claudio↗
NASA NTRS1996

An Independent Evaluation of the FMEA/CIL Hazard Analysis Alternative Study

The present instruments of safety and reliability risk control for a majority of the National Aeronautics and Space Administration (NASA) programs/projects consist of Failure Mode and Effects Analysis (FMEA), Hazard Analysis (HA), Critical Items List (CIL), and Hazard Report (HR). This extensive analytical approach was introduced in the early 1970's and was implemented for the Space Shuttle Program by NHB 5300.4 (1D-2. Since the Challenger accident in 1986, the process has been expanded considerably and resulted in introduction of similar and/or duplicated activities in the safety/reliability risk analysis. A study initiated in 1995, to search for an alternative to the current FMEA/CIL Hazard Analysis methodology generated a proposed method on April 30, 1996. The objective of this Summer Faculty Study was to participate in and conduct an independent evaluation of the proposed alternative to simplify the present safety and reliability risk control procedure.

Ray, Paul S.↗
NASA NTRS1999

Simulation of the December 1998 Stratospheric Major Warming

Prior to 1991, major warmings (defined by increasing zonal mean temperatures and zonal mean easterly winds from 60degN to the pole at 10 hPa) typically occurred approximately once every two Arctic winters; a major warming in mid-Dec. 1998 was the first since Feb. 1991. The Dec. 1998 warming was also the second earliest on record. The earliest, and the only other major warming on record before the end of Dec. was in early Dec 1987; prior to that, the earliest was in late Dec./early Jan. 1984-85. The 1984-85 and 1987 warmings resulted in the warmest and weakest lower stratospheric polar vortices in the 20 years before 1998-99. Fig. 1 compares temperatures and vortex strength in 1998-99 with those in the previous 20 years, using the US National Center for Environmental Prediction (NCEP) record; 1987-88 and 1984-85 are also highlighted. The Dec. 1998 warming had a more pronounced effect on mid-stratospheric temperatures than the Dec. 1987 warming (Fig. 1a), although smaller than that of warmings later in winter (e.g., 1984-85). 10-hPa temperatures fell well below average again in late Jan. 1999 and remained unusually low until an early final warming began in late Feb. 840 K PV gradients (Fig. 1c) set a record minimum in Jan. 1999, but were near average in Feb before the final warming. The effect of the Dec. 1998 warming on lower stratospheric temperatures was comparable to that of other major warmings; there was a brief period of record-high minimum 46-hPa temperatures in early Jan 1999 (Fig. 1b), and temperatures then fell to near average for a short period in mid-Feb. Lower stratospheric PV gradients were the weakest on record during the 1998-99 winter (Fig. 1d). The evolution of the vortex and minimum temperatures during 1998-99 was remarkably similar to that during 1987-88, the only previous year when a major warming was observed before the end of Dec.

Manney, G. L.↗
NASA NTRS1999

A Theory of Oscillating Edge Flames

It has been known for some years that when a near-limit flame spreads over a liquid pool of fuel, the edge of the flame can oscillate relative to a frame moving with the mean speed. Each period of oscillation is characterized by long intervals of modest motion during which the edge gases radiate like those of a diffusion flame, punctuated by bursts of rapid advance during which the edge gases radiate like those in a deflagration. Substantial resources have been brought to bear on this issue within the microgravity program, both experimental and numerical. It is also known that when a near-asphyxiated candle-flame burns at zero gravity, the edge of the (hemispherical) flame can oscillate violently prior to extinction. Thus a web-surfer, turning to the NASA web-site at http://microgravity.msfc.nasa.gov, and following the trail combustion science/experiments/experimental results/candle flame, will find photographs and a description of candle burning experiments carried out on board both the Space-shuttle and the Russian space station Mir. A brief report can also be found in the proceedings of the Fourth Workshop. And recently, in a third microgravity program, the leading edge of the flame supported by injection of ethane through the porous surface of a plate over which air is blown has been found to oscillate when conditions are close to blow-off. A number of important points can be made with respect to these observations: It is the edge itself which oscillates, advancing and retreating, not the diffusion flame that trails behind the edge; oscillations only occur under near limit conditions; in each case the Lewis number of the fuel is significantly larger than 1; and because of the edge curvature, the heat losses from the reacting edge structure are larger than those from the trailing diffusion flame. We propose a general theory for these oscillations, invoking Occam's 'Law of Parsimony' in an expanded form, to wit: The same mechanism is responsible for the oscillations in all three experiments; and no new mechanism is invoked (Occam's original 'Razor'). Such a strategy eliminates Marangoni effects as the source, for these are absent in the second and third experiments. And it eliminates arguments that point to numerically predicted gas eddies as the source, a new mechanism, unelucidated. Indeed, we hypothesize that the essential driving mechanism for the instability is a combination of large Lewis number and heat losses from the reacting structure near the flame edge. Instabilities driven by these mechanisms are commonplace in 1D configurations. Chemical reactor theory, for example, leads to system responses which mimic the response of the candle flame - steady flame, oscillations, extinction. In a combustion context, oscillating instabilities were first reported for diffusion flames in a theoretical study by Kirkby and Schmitz, and here also the instabilities are associated with near-extinction conditions, large Lewis numbers, and heat losses. And deflagrations will oscillate if the Lewis number is large enough, oscillations that are exacerbated when heat losses are present, whether global or to a surface.

Buckmaster, J.↗
NASA NTRS1999

Microwave Radiometric Signatures of Different Surface Types in Deserts

In arid environments, specific microwave signatures have been observed with the Special Sensor Microwave/Imager (SSM/I). For a given diurnal change in surface skin temperature, the corresponding change in the microwave brightness temperature is smaller than expected. With the help of a 1D, time-dependent heat conduction model, this behavior is explained by microwave radiation coming from different depths in the soil, depending on the soil type and on the microwave radiation frequency. Using the eight-times daily estimates of the surface skin temperature by the International Satellite Cloud Climatology Project (ISCCP) and a simple Fresnel model, collocated month-long time series of the SSM/I brightness temperatures and the surface skin temperatures give a consistent estimate of the effective microwave emissivity and penetration depth parameters. Results are presented and analyzed for the Sahara and the Arabian Peninsula, for July and November 1992. The case of the Australian desert is also briefly mentioned. Assuming a reasonable thermal diffusivity for the soil in desert areas, the microwave radiation is estimated to come from soil layers down to depths of at least five wavelengths in some locations. Regions where the microwave radiation comes from deeper soil layers also have large microwave emissivity polarization differences and large visible reflectances, suggesting that these areas correspond to sand dune fields.

Prigent, Catherine↗
NASA NTRS1999

Low-Storage, Explicit Runge-Kutta Schemes for the Compressible Navier-Stokes Equations

The derivation of storage explicit Runge-Kutta (ERK) schemes has been performed in the context of integrating the compressible Navier-Stokes equations via direct numerical simulation. Optimization of ERK methods is done across the broad range of properties, such as stability and accuracy efficiency, linear and nonlinear stability, error control reliability, step change stability, and dissipation/dispersion accuracy, subject to varying degrees of memory economization. Following van der Houwen and Wray, 16 ERK pairs are presented using from two to five registers of memory per equation, per grid point and having accuracies from third- to fifth-order. Methods have been assessed using the differential equation testing code DETEST, and with the 1D wave equation. Two of the methods have been applied to the DNS of a compressible jet as well as methane-air and hydrogen-air flames. Derived 3(2) and 4(3) pairs are competitive with existing full-storage methods. Although a substantial efficiency penalty accompanies use of two- and three-register, fifth-order methods, the best contemporary full-storage methods can be pearl), matched while still saving two to three registers of memory.

Kennedy, Chistopher A.↗
NASA NTRS1998

Temperature Dependence of Diffusivities in Liquid Elements (LMD)

This research was to advance the understanding of diffusion mechanisms in liquid metals and alloys through accurate diffusivity measurements over a wide range of temperatures, including the proximity of the materials melting points. Specifically, it was driven towards developing a methodology (and subsequent flight hardware) to enable several diffusion coefficient measurements (i.e., at several different temperatures) to be performed using a single sample. The Liquid Metal Diffusion (LMD) was funded as a Flight Definition Project in February 1993 in response to NRA 91-OSSA-20 (Microgravity Science and Applications Division). The Science Concept Review for LAID was held during April 1994. In January 1995 we were informed that we had failed this review and the project was change to ground-based activities only. A new proposal was submitted for the next NRA addressing the panels concerns. As part of NASA's Risk Mitigation program, a scaled-down version of the hardware was funded in July of 1995 for a flight opportunity utilizing experiment on the Microgravity Isolation Mount. This experiment was to determine the self-diffusivity of indium at 185 C. The LMD was transferred to the Mir Space Station in STS-81 and returned on STS-84 (January - May 1997). Three, out of five, self-diffusion data sets were returned. A description of this experiment/hardware is included below. This summary is only intended to give the reader an overview of the results obtained for the tasks outlined in the original proposal. Research that was not published is explained in more detail. At the end of this report is a list of refereed publications and invited talks that were given as a result of this work. The reader is directed to these for further details. Attachment: Real-time diffusivity measurements in liquids at several temperatures with one sample, On the insensitivity of liquid diffusivity measurements to deviations from 1D transport, and Numerical simulations of the convective contamination of diffusivity measurements in liquids.

Banish, R. Michael↗
NASA NTRS2000

Multiresolution With Super-Compact Wavelets

The solution data computed from large scale simulations are sometimes too big for main memory, for local disks, and possibly even for a remote storage disk, creating tremendous processing time as well as technical difficulties in analyzing the data. The excessive storage demands a corresponding huge penalty in I/O time, rendering time and transmission time between different computer systems. In this paper, a multiresolution scheme is proposed to compress field simulation or experimental data without much loss of important information in the representation. Originally, the wavelet based multiresolution scheme was introduced in image processing, for the purposes of data compression and feature extraction. Unlike photographic image data which has rather simple settings, computational field simulation data needs more careful treatment in applying the multiresolution technique. While the image data sits on a regular spaced grid, the simulation data usually resides on a structured curvilinear grid or unstructured grid. In addition to the irregularity in grid spacing, the other difficulty is that the solutions consist of vectors instead of scalar values. The data characteristics demand more restrictive conditions. In general, the photographic images have very little inherent smoothness with discontinuities almost everywhere. On the other hand, the numerical solutions have smoothness almost everywhere and discontinuities in local areas (shock, vortices, and shear layers). The wavelet bases should be amenable to the solution of the problem at hand and applicable to constraints such as numerical accuracy and boundary conditions. In choosing a suitable wavelet basis for simulation data among a variety of wavelet families, the supercompact wavelets designed by Beam and Warming provide one of the most effective multiresolution schemes. Supercompact multi-wavelets retain the compactness of Haar wavelets, are piecewise polynomial and orthogonal, and can have arbitrary order of approximation. The advantages of the multiresolution algorithm are that no special treatment is required at the boundaries of the interval, and that the application to functions which are only piecewise continuous (internal boundaries) can be efficiently implemented. In this presentation, Beam's supercompact wavelets are generalized to higher dimensions using multidimensional scaling and wavelet functions rather than alternating the directions as in the 1D version. As a demonstration of actual 3D data compression, supercompact wavelet transforms are applied to a 3D data set for wing tip vortex flow solutions (2.5 million grid points). It is shown that high data compression ratio can be achieved (around 50:1 ratio) in both vector and scalar data set.

Lee, Dohyung↗
NASA NTRS1999

Corals from Space

The goal of this research is to monitor the health and vigor of coral reef ecosystems, and their sensitivity to natural and anthropogenic climate changes. To achieve these lofty goals, this research is investigating the feasibility of using spaceborne high-resolution spectrometers (on the US Landsat, French Systeme Probatoire pour l'Observation de la Terre [SPOT] and/or the Indian Resources Satellite [IRS 1C & 1D] spacecraft) to first map the aerial extent of coral reef systems, and second separate the amount of particular corals. If this is successful, we could potentially provide a quantum leap in our understanding of coral reef systems, as well as provide much needed baseline data to measure future changes in global coral reef ecosystems. In collaboration with Tomas Tomascik, Yann Morel, and other colleagues, a series of experiments were planned to coordinate in situ coral observations, high-resolution spaceborne imagery (from Landsat, SPOT, and, possibly, IRS IC spacecraft), and NASA Space Shuttle photographs and digital images. Our eventual goal is to develop "coral health algorithms" that can be used to assess time series of imagery collected from satellite sensors (Landsat since 1972, SPOT since 1986) in concert with in situ observations. The bad news from last year was that from 1997 to mid- 1998, the extreme cloudiness over southeast Asia due to prolonged smoke from El Nino-related fires and the economic chaos in this region frustrated both our space and reef-based data collection activities. When this volatile situation stabilizes, we will restart these activities. The good news was that in collaboration with Al Strong at the National Oceanic and Atmospheric Administration (NOAA) we had an exciting year operationally using the NOAA's Advanced Very High Resolution Radiometer sensor derived sea surface temperature products to warn of coral "bleaching" at many locations throughout the tropics. Data from NOAA's satellites showed that during the El Nino of 1997 and the first half of 1998, more ocean area in the tropics experienced exceptionally high sea surface temperatures, or "hot spots," than have been observed in any fall year since the El Nino of 1982. From January to July, the coral bleaching events were concentrated in the Southern Hemisphere (during its warm season). Since July 1998, the reports of extensive coral bleaching have again spread into regions of the Northern Hemisphere following abnormally high sea surface temperatures, especially around the Philippines and throughout the Caribbean Basin, Bahamas, Bermuda and Florida Keys. These El Nino induced events clearly demonstrated that corals are the "canaries of the marine ecosystem," highly sensitive to short-term natural climate events (El Nino), and should be monitored as measures of longer-term environmental and climate change. Additional information is contained in the original.

Patzert, William C.↗
NASA NTRS2000

Isotopic NO as a Chemical Tracer in the Global Stratosphere

Stratospheric NO originates from nitrous oxide reacting with O(1D) and ion-molecule reactions. Most ionic reactions take place in the mesosphere and lower thermosphere. The resulting NO is transported into the stratosphere at high latitudes. Cosmic radiation and tropospheric lightning also produce nitric oxide. This NO originates from ion reactions involving N2. Ionic reactions preserve the N(15)/N(14) ratio present in atmospheric N2. Nitrous oxide has a mass-dependent sink that varies with altitude so that there is an altitude-dependent isotopic distinction in nitrous oxide. This difference will appear in NO formed from N2O. The expected NO isotopic distribution under different conditions will be a combination of NO derived from nitrous oxide with different masses and NO from ion reactions. The expected NO isotopic distribution will be presented talking into account the different processes, including particle events and downward transport in winter.

Aikin, A. C.↗