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 487 records · Page 27

Courant Number and Mach Number Insensitive CE/SE Euler Solvers

It has been known that the space-time CE/SE method can be used to obtain ID, 2D, and 3D steady and unsteady flow solutions with Mach numbers ranging from 0.0028 to 10. However, it is also known that a CE/SE solution may become overly dissipative when the Mach number is very small. As an initial attempt to remedy this weakness, new 1D Courant number and Mach number insensitive CE/SE Euler solvers are developed using several key concepts underlying the recent successful development of Courant number insensitive CE/SE schemes. Numerical results indicate that the new solvers are capable of resolving crisply a contact discontinuity embedded in a flow with the maximum Mach number = 0.01.

Chang, Sin-Chung↗

Energetic Metastable Oxygen and Nitrogen Atoms in the Terrestrial Atmosphere

This report summarizes our research performed under NASA Grant NAG5-11857. The three-year grant have been supported by the Geospace Sciences SR&T program. We have investigated the energetic metastable oxygen and nitrogen atoms in the terrestrial stratosphere, mesosphere and thermosphere. Hot atoms in the atmosphere are produced by solar radiation, the solar wind and various ionic reactions. Nascent hot atoms arise in ground and excited electronic states, and their translational energies are larger by two - three orders of magnitude than the thermal energies of the ambient gas. The relaxation kinetics of hot atoms determines the rate of atmospheric heating, the intensities of aeronomic reactions, and the rate of atom escape from the planet. Modeling of the non-Maxwellian energy distributions of metastable oxygen and nitrogen atoms have been focused on the determination of their impact on the energetics and chemistry of the terrestrial atmosphere between 25 and 250 km . At this altitudes, we have calculated the energy distribution functions of metastable O and N atoms and computed non-equilibrium rates of important aeronomic reactions, such as destruction of the water molecules by O(1D) atoms and production of highly excited nitric oxide molecules. In the upper atmosphere, the metastable O(lD) and N(2D) play important role in formation of the upward atomic fluxes. We have computed the upward fluxes of the metastable and ground state oxygen atoms in the upper atmosphere above 250 km. The accurate distributions of the metastable atoms have been evaluated for the day and night-time conditions.

Kharchenko, Vasili↗

Single Axis Flywheel IPACS @1300W, 0.8 N-m

NASA Glenn Research Center is developing flywheels for space systems. A single axis laboratory version of an integrated power and attitude control (IPACs) system has been experimentally demonstrated. This is a significant step on the road to a flight qualified three axes IPACS system. The presentation outlines the flywheel development process at NASA GRC, the experimental hardware and approach, the IPACS control algorithm that was formulated and the results of the test program and then proposes a direction for future work. GRC has made progress on flywheel module design in terms of specific energy density and capability through a design and test program resulting in three flywheel module designs. Two of the flywheels are used in the 1D-IPACS experiment with loads and power sources to simulate a satellite power system. The system response is measured in three power modes: charge, discharge, and charge reduction while simultaneously producing a net output torque which could be used for attitude control. Finally, recommendations are made for steps that should be taken to evolve from this laboratory demonstration to a flight like system.

Jansen, Ralph↗

Measurements of Thermal Conductivity of Superfluid Helium Near its Transition Temperature T(sub lambda) in a 2D Confinement

We report our recent experiments on thermal conductivity measurements of superfluid He-4 near its phase transition in a two-dimensional (2D) confinement under saturated vapor pressure. A 2D confinement is created by 2-mm- and 1-mm-thick glass capillary plates, consisting of densely populated parallel microchannels with cross-sections of 5 x 50 and 1 x 10 microns, correspondingly. A heat current (2 < Q < 400 nW/sq cm) was applied along the channels long direction. High-resolution measurements were provided by DC SQUID-based high-resolution paramagnetic salt thermometers (HRTs) with a nanokelvin resolution. We might find that thermal conductivity of confined helium is finite at the bulk superfluid transition temperature. Our 2D results will be compared with those in a bulk and 1D confinement.

Jerebets, Sergei↗

Benchmark study for charge deposition by high energy electrons in thick slabs

The charge deposition profiles created when highenergy (1, 10, and 100 MeV) electrons impinge ona thick slab of elemental aluminum, copper, andtungsten are presented in this paper. The chargedeposition profiles were computed using existing representative Monte Carlo codes: TIGER3.0 (1D module of ITS3.0) and MCNP version 4B. The results showed that TIGER3.0 and MCNP4B agree very well (within ~20% of each other) in the majority of the problem geometry. The TIGER results were considered to be accurate based on previous studies. Thus, it was demonstrated that MCNP, with its powerful geometry capability and flexible source and tally options, could be used in calculations of electron charging in high energy electron-rich space radiation environments.

charging high energy electrons ITS MCNP↗

Runtime and Pressurization Analyses of Propellant Tanks

Multi-element unstructured CFD has been utilized at NASA SSC to carry out analyses of propellant tank systems in different modes of operation. The three regimes of interest at SSC include (a) tank chill down (b) tank pressurization and (c) runtime propellant draw-down and purge. While tank chill down is an important event that is best addressed with long time-scale heat transfer calculations, CFD can play a critical role in the tank pressurization and runtime modes of operation. In these situations, problems with contamination of the propellant by inclusion of the pressurant gas from the ullage causes a deterioration of the quality of the propellant delivered to the test article. CFD can be used to help quantify the mixing and propellant degradation. During tank pressurization under some circumstances, rapid mixing of relatively warm pressurant gas with cryogenic propellant can lead to rapid densification of the gas and loss of pressure in the tank. This phenomenon can cause serious problems during testing because of the resulting decrease in propellant flow rate. With proper physical models implemented, CFD can model the coupling between the propellant and pressurant including heat transfer and phase change effects and accurately capture the complex physics in the evolving flowfields. This holds the promise of allowing the specification of operational conditions and procedures that could minimize the undesirable mixing and heat transfer inherent in propellant tank operation. It should be noted that traditional CFD modeling is inadequate for such simulations because the fluids in the tank are in a range of different sub-critical and supercritical states and elaborate phase change and mixing rules have to be developed to accurately model the interaction between the ullage gas and the propellant. We show a typical run-time simulation of a spherical propellant tank, containing RP-1 in this case, being pressurized with room-temperature nitrogen at 540 R. Nitrogen, shown in blue on the right-hand side of the figures, enters the tank from the diffuser at the top of the figures and impinges on the RP-1, shown in red, while the propellant is being continuously drained at the rate of 1050 lbs/sec through a pipe at the bottom of the tank. The sequence of frames in Figure 1 shows the resultant velocity fields and mixing between nitrogen and RP-1 in a cross-section of the tank at different times. A vortex is seen to form in the incoming nitrogen stream that tends to entrain propellant, mixing it with the pressurant gas. The RP-1 mass fraction contours in Figure 1 are also indicative of the level of mixing and contamination of the propellant. The simulation is used to track the propagation of the pure propellant front as it is drawn toward the exit with the evolution of the mixing processes in the tank. The CFD simulation modeled a total of 10 seconds of run time. As is seen from Figure 1d, after 5.65 seconds the propellant front is nearing the drain pipe, especially near the center of the tank. Behind this pure propellant front is a mixed fluid of compromised quality that would require the test to end when it reaches the exit pipe. Such unsteady simulations provide an estimate of the time that a high-quality propellant supply to the test article can be guaranteed at the modeled mass flow rate. In the final paper, we will discuss simulations of the LOX and propellant tanks at NASA SSC being pressurized by an inert ullage. Detailed comparisons will be made between the CFD simulations and lower order models as well as with test data. Conditions leading to cryo collapse in the tank will also be identified.

Field, Robert E.↗

Aeromechanics Analysis of a Compound Helicopter

A design and aeromechanics investigation was conducted for a 100,000-lb compound helicopter with a single main rotor, which is to cruise at 250 knots at 4000 ft/95 deg F condition. Performance, stability, and control analyses were conducted with the comprehensive rotorcraft analysis CAMRAD II. Wind tunnel test measurements of the performance of the H-34 and UH-1D rotors at high advance ratio were compared with calculations to assess the accuracy of the analysis for the design of a high speed helicopter. In general, good correlation was obtained with the increase of drag coefficients in the reverse flow region. An assessment of various design parameters (disk loading, blade loading, wing loading) on the performance of the compound helicopter was made. Performance optimization was conducted to find the optimum twist, collective, tip speed, and taper using the comprehensive analysis. Blade twist was an important parameter on the aircraft performance and most of the benefit of slowing the rotor occurred at the initial 20 to 30% reduction of rotor tip speed. No stability issues were observed with the current design and the control derivatives did not change much with speed, but did exhibit significant coupling.

Yeo, Hyeonsoo↗

Optimum Design of a Compound Helicopter

A design and aeromechanics investigation was conducted for a 100,000-lb compound helicopter with a single main rotor, which is to cruise at 250 knots at 4000 ft/95 deg F condition. Performance, stability, and control analyses were conducted with the comprehensive rotorcraft analysis CAMRAD II. Wind tunnel test measurements of the performance of the H-34 and UH-1D rotors at high advance ratio were compared with calculations to assess the accuracy of the analysis for the design of a high speed helicopter. In general, good correlation was obtained when an increase of drag coefficients in the reverse flow region was implemented. An assessment of various design parameters (disk loading, blade loading, wing loading) on the performance of the compound helicopter was conducted. Lower wing loading (larger wing area) and higher blade loading (smaller blade chord) increased aircraft lift-to-drag ratio. However, disk loading has a small influence on aircraft lift-to-drag ratio. A rotor parametric study showed that most of the benefit of slowing the rotor occurred at the initial 20 to 30% reduction of the advancing blade tip Mach number. No stability issues were observed with the current design. Control derivatives did not change significantly with speed, but the did exhibit significant coupling.

Yeo, Hyeonsoo↗

Electron Collisional Excitation Rates for OI USING the B-Spline R-Matrix Approach

The B-spline R-matrix approach has been used to calculate electron collisional excitation strengths and rates for transitions between the 3P, 1D, and 1S states of ground configuration and from these states to the states of the excited 2s(sup 2)2p(sup 3)ns (n = 3-5), 2s(sup 2)2p(sup 3)np (n = 3-4), 2s(sup 2)2p(sup 3)nd (n = 3-4), 2s(sup 2)2p(sup 3)4f, and 2s2p(sup 5) configurations. The nonorthogonal orbitals are used for an accurate description of both the target wave functions and the R-matrix basis functions. The thermally averaged collision strengths are obtained from the collision strengths by integrating over a Maxwellian velocity distribution of electron energies, and these are tabulated over a temperature range from 1000 to 60,000 K. The parametric functions of scaled energy have also been obtained to represent collision strengths over a wide energy range or thermally averaged collision strengths at any desired temperature.

Zatsarinny, O.↗

On the Stratospheric Chemistry of Hydrogen Cyanide

HCN profiles measured by solar occultation spectrometry during 10 balloon flights of the JPL MkIV instrument are presented. The HCN profiles reveal a compact correlation with stratospheric tracers. Calculations with a 2D-model using established rate coefficients for the reactions of HCN with OH and O(1D) severely underestimate the measured HCN in the middle and upper stratosphere. The use of newly available rate coefficients for these reactions gives reasonable agreement of measured and modeled HCN. An HCN yield of approx.30% from the reaction of CH3CN with OH is consistent with the measurements.

photochemical properties↗

The Attenuation of a Detonation Wave by an Aircraft Engine Axial Turbine Stage

A Constant Volume Combustion Cycle Engine concept consisting of a Pulse Detonation Combustor (PDC) followed by a conventional axial turbine was simulated numerically to determine the attenuation and reflection of a notional PDC pulse by the turbine. The multi-stage, time-accurate, turbomachinery solver TURBO was used to perform the calculation. The solution domain consisted of one notional detonation tube coupled to 5 vane passages and 8 rotor passages representing 1/8th of the annulus. The detonation tube was implemented as an initial value problem with the thermodynamic state of the tube contents, when the detonation wave is about to exit, provided by a 1D code. Pressure time history data from the numerical simulation was compared to experimental data from a similar configuration to verify that the simulation is giving reasonable results. Analysis of the pressure data showed a spectrally averaged attenuation of about 15 dB across the turbine stage. An evaluation of turbine performance is also presented.

VanZante, Dale↗

Time-resolved diffusion tomographic 2D and 3D imaging in highly scattering turbid media

A method for imaging objects in highly scattering turbid media. According to one embodiment of the invention, the method involves using a plurality of intersecting source/detectors sets and time-resolving equipment to generate a plurality of time-resolved intensity curves for the diffusive component of light emergent from the medium. For each of the curves, the intensities at a plurality of times are then inputted into the following inverse reconstruction algorithm to form an image of the medium: wherein W is a matrix relating output at source and detector positions r.sub.s and r.sub.d, at time t, to position r, .LAMBDA. is a regularization matrix, chosen for convenience to be diagonal, but selected in a way related to the ratio of the noise, to fluctuations in the absorption (or diffusion) X.sub.j that we are trying to determine: .LAMBDA..sub.ij =.lambda..sub.j .delta..sub.ij with .lambda..sub.j = /<.DELTA.Xj.DELTA.Xj> Y is the data collected at the detectors, and X.sup.k is the kth iterate toward the desired absorption information. An algorithm, which combines a two dimensional (2D) matrix inversion with a one-dimensional (1D) Fourier transform inversion is used to obtain images of three dimensional hidden objects in turbid scattering media.

Alfano, Robert R.↗

Time-resolved diffusion tomographic 2D and 3D imaging in highly scattering turbid media

A method for imaging objects in highly scattering turbid media. According to one embodiment of the invention, the method involves using a plurality of intersecting source/detectors sets and time-resolving equipment to generate a plurality of time-resolved intensity curves for the diffusive component of light emergent from the medium. For each of the curves, the intensities at a plurality of times are then inputted into the following inverse reconstruction algorithm to form an image of the medium: ##EQU1## wherein W is a matrix relating output at source and detector positions r.sub.s and r.sub.d, at time t, to position r, .LAMBDA. is a regularization matrix, chosen for convenience to be diagonal, but selected in a way related to the ratio of the noise, to fluctuations in the absorption (or diffusion) X.sub.j that we are trying to determine: .LAMBDA..sub.ij =.lambda..sub.j .delta..sub.ij with .lambda..sub.j = /<.DELTA.Xj.DELTA.Xj> Y is the data collected at the detectors, and X.sup.k is the kth iterate toward the desired absoption information. An algorithm, which combines a two dimensional (2D) matrix inversion with a one-dimensional (1D) Fourier transform inversion is used to obtain images of three dimensional hidden objects in turbid scattering media.

Alfano, Robert R.↗

Parallel Anisotropic Tetrahedral Adaptation

An adaptive method that robustly produces high aspect ratio tetrahedra to a general 3D metric specification without introducing hybrid semi-structured regions is presented. The elemental operators and higher-level logic is described with their respective domain-decomposed parallelizations. An anisotropic tetrahedral grid adaptation scheme is demonstrated for 1000-1 stretching for a simple cube geometry. This form of adaptation is applicable to more complex domain boundaries via a cut-cell approach as demonstrated by a parallel 3D supersonic simulation of a complex fighter aircraft. To avoid the assumptions and approximations required to form a metric to specify adaptation, an approach is introduced that directly evaluates interpolation error. The grid is adapted to reduce and equidistribute this interpolation error calculation without the use of an intervening anisotropic metric. Direct interpolation error adaptation is illustrated for 1D and 3D domains.

Park, Michael A.↗

Hybrid-dual-fourier tomographic algorithm for a fast three-dimensionial optical image reconstruction in turbid media

A reconstruction technique for reducing computation burden in the 3D image processes, wherein the reconstruction procedure comprises an inverse and a forward model. The inverse model uses a hybrid dual Fourier algorithm that combines a 2D Fourier inversion with a 1D matrix inversion to thereby provide high-speed inverse computations. The inverse algorithm uses a hybrid transfer to provide fast Fourier inversion for data of multiple sources and multiple detectors. The forward model is based on an analytical cumulant solution of a radiative transfer equation. The accurate analytical form of the solution to the radiative transfer equation provides an efficient formalism for fast computation of the forward model.

Alfano, Robert R.↗

The a(3) Scheme--A Fourth-Order Space-Time Flux-Conserving and Neutrally Stable CESE Solver

The CESE development is driven by a belief that a solver should (i) enforce conservation laws in both space and time, and (ii) be built from a non-dissipative (i.e., neutrally stable) core scheme so that the numerical dissipation can be controlled effectively. To initiate a systematic CESE development of high order schemes, in this paper we provide a thorough discussion on the structure, consistency, stability, phase error, and accuracy of a new 4th-order space-time flux-conserving and neutrally stable CESE solver of an 1D scalar advection equation. The space-time stencil of this two-level explicit scheme is formed by one point at the upper time level and three points at the lower time level. Because it is associated with three independent mesh variables (the numerical analogues of the dependent variable and its 1st-order and 2ndorder spatial derivatives, respectively) and three equations per mesh point, the new scheme is referred to as the a(3) scheme. Through the von Neumann analysis, it is shown that the a(3) scheme is stable if and only if the Courant number is less than 0.5. Moreover, it is established numerically that the a(3) scheme is 4th-order accurate.

Chang, Sin-Chung↗

Emission-Line Galaxies from the PEARS Hubble Ultra Deep Field: A 2-D Detection Method and First Results

The Hubble Space Telescope (HST) Advanced Camera for Surveys (ACS) grism PEARS (Probing Evolution And Reionization Spectroscopically) survey provides a large dataset of low-resolution spectra from thousands of galaxies in the GOODS North and South fields. One important subset of objects in these data are emission-line galaxies (ELGs), and we have investigated several different methods aimed at systematically selecting these galaxies. Here we present a new methodology and results of a search for these ELGs in the PEARS observations of the Hubble Ultra Deep Field (HUDF) using a 2D detection method that utilizes the observation that many emission lines originate from clumpy knots within galaxies. This 2D line-finding method proves to be useful in detecting emission lines from compact knots within galaxies that might not otherwise be detected using more traditional 1D line-finding techniques. We find in total 96 emission lines in the HUDF, originating from 81 distinct "knots" within 63 individual galaxies. We find in general that [0 1111 emitters are the most common, comprising 44% of the sample, and on average have high equivalent widths (70% of [0 1111 emitters having rest-frame EW> 100A). There are 12 galaxies with multiple emitting knots; several show evidence of variations in H-alpha flux in the knots, suggesting that the differing star formation properties across a single galaxy can in general be probed at redshifts approximately greater than 0.2 - 0.4. The most prevalent morphologies are large face-on spirals and clumpy interacting systems, many being unique detections owing to the 2D method described here, thus highlighting the strength of this technique.

Gardner, J. P.↗

The Formation of Coronal Loops by Thermal Instability in Three Dimensions

Plasma loops in solar active regions have been observed in EUV and soft X-rays for decades. Their formation mechanism and properties, however, are still not fully understood. Predictions by early models, based on 1D hydrostatic equilibria with uniform plasma heating, are not consistent with high-resolution measurements. In this Letter, we demonstrate, via 3D simulations, that a class of heating models can lead to the dynamic formation of plasma loops provided the plasma is heated sufficiently to match SXT soft X-ray measurements. We show that individual flux tubes in a 3D magnetic structure tend to stand out against their neighbors. The loops have large aspect ratios and nearly uniform cross sections in the corona, similar to those observed by EIT and TRACE. The coronal EUV emission from these thermally unstable solutions is roughly consistent with EIT measurements. The solution oscillates in time through a large-amplitude, nonlinear cycle, leading to repeated brightening and fading of the loops.

Mok, Yung↗