Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “D CODES”

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 397 records · Page 22

Testing the pressure-confined Ly alpha cloud model

The Ly-alpha absorption line forest, seen in quasar spectra, is generally interpreted as being due to cosmologically distributed 'clouds' of primordial gas. Analyses of the observations reveal that the number distribution can be described by power laws: dN/dz is proportional to (1 + z)(sup gamma) and dN/dN(sub HI) is proportional to N(sub HI)(sup -beyda), where N(sub HI) is the HI column density. The typical values for power law indices range between 2 is approximately less than gamma is approximately less than 2.6 and 1.7 is approximately less than gamma is approximately less than 1.9. One model postulates that the Ly-alpha clouds are optically thin entities, photoionized by the background UV flux, J(sub nu) is proportional to (1 + z)(sup j), and confined by an adiabatically evolving intercloud medium (ICM): P(z) is proportional to (1 + z)(sup 5). Analytic studies of this model suggest that the ensuing Ly-alpha line statistics can account for the observations (in particular, the dN/dz and the dN/dN(sub HI) distributions) if the cloud mass spectrum is a power law dN/dN is proportional to M(sup -delta), delta is approximately 1.9, and j is approximately 4. One of the simplifying assumptions incorporated into these studies is the existence of a large mass range for the clouds at all epochs, the validity of which is questionable. The pressure-confined model is investigated using a 1-D spherically symmetric hydrodynamical code to simulate cloud evolution over the epoch 1.8 less than z less than 6. This enables us to relax many of the assumptions incorporated in the analytic studies.

Babul, Arif↗

The first ultraviolet spectropolarimetric study of NGC 1068

We present the first ultraviolet spectropolarimetric observations of the Type 2 Seyfert galaxy NGC 1068 obtained by the Wisconsin Ultraviolet Photo Polarimeter Experiment (WUPPE). In this communication we discuss the continuum polarization in the ultraviolet. Shortward of 2500 A the polarization becomes constant at 16 percent at a position angle of 112 deg. This is in general agreement with the 1983 prediction of Miller and Antonucci (1983) and their thick torus model, in which a hidden Type 1 Seyfert nucleus is seen in light scattered by electrons outside the obscured nuclei.

Code, Arthur D.↗

Ultraviolet spectropolarimetry of the Be star PP Carinae with the Wisconsin Ultraviolet Photo-Polarimeter Experiment

We present the first ultraviolet spectropolarimetric observations of the Be star PP Car, obtained with the Wisconsin Ultraviolet Photo-Polarimeter Experiment (WUPPE) aboard the Astro 1 mission. Usable polarization data were obtained from 1400 to 2330 A, along with a good spectrum from 1400 to 3200 A. These data show a lower polarization shortward of the Balmer jump than had been predicted by standard models, and a broad UV polarization dip around 1900 A is seen. These results are in agreement with those found from the WUPPE observations of two other Be stars, Xi Tau and Pi Aqr, which were published earlier. All these observations are an important probe of the Be circumstellar envelopes and demonstrate the need for the inclusion of metal-line effects in circumstellar disk models of Be star UV polarization.

Bjorkman, K. S.↗

Advanced electromagnetic methods for aerospace vehicles

The Advanced Helicopter Electromagnetics (AHE) Industrial Associates Program has fruitfully completed its fourth year. Under the support of the AHE members and the joint effort of the research team, new and significant progress has been achieved in the year. Following the recommendations by the Advisory Task Force, the research effort is placed on more practical helicopter electromagnetic problems, such as HF antennas, composite materials, and antenna efficiencies. In this annual report, the main topics to be addressed include composite materials and antenna technology. The research work on each topic has been driven by the AHE consortium members' interests and needs. The remarkable achievements and progresses in each subject is reported respectively in individual sections of the report. The work in the area of composite materials includes: modeling of low conductivity composite materials by using Green's function approach; guidelines for composite material modeling by using the Green's function approach in the NEC code; development of 3-D volume mesh generator for modeling thick and volumetric dielectrics by using FD-TD method; modeling antenna elements mounted on a composite Comanche tail stabilizer; and antenna pattern control and efficiency estimate for a horn antenna loaded with composite dielectric materials.

Balanis, Constantine A.↗

Low Reynolds number kappa-epsilon and empirical transition models for oscillatory pipe flow and heat transfer

Stirling engine heat exchangers are shell-and-tube type with oscillatory flow (zero-mean velocity) for the inner fluid. This heat transfer process involves laminar-transition turbulent flow motions under oscillatory flow conditions. A low Reynolds number kappa-epsilon model, (Lam-Bremhorst form), was utilized in the present study to simulate fluid flow and heat transfer in a circular tube. An empirical transition model was used to activate the low Reynolds number k-e model at the appropriate time within the cycle for a given axial location within the tube. The computational results were compared with experimental flow and heat transfer data for: (1) velocity profiles, (2) kinetic energy of turbulence, (3) skin friction factor, (4) temperature profiles, and (5) wall heat flux. The experimental data were obtained for flow in a tube (38 mm diameter and 60 diameter long), with the maximum Reynolds number based on velocity being Re(sub max) = 11840, a dimensionless frequency (Valensi number) of Va = 80.2, at three axial locations X/D = 16, 30 and 44. The agreement between the computations and the experiment is excellent in the laminar portion of the cycle and good in the turbulent portion. Moreover, the location of transition was predicted accurately. The Low Reynolds Number kappa-epsilon model, together with an empirical transition model, is proposed herein to generate the wall heat flux values at different operating parameters than the experimental conditions. Those computational data can be used for testing the much simpler and less accurate one dimensional models utilized in 1-D Stirling Engine design codes.

Bauer, Christopher↗

Numerical simulation of cosmogenic nuclide production in lunar rocks

The production rates of cosmogenic nuclides depend on the primary cosmic-ray particles, the irradiated-body's bulk composition, size, and shape, and the sample's composition and shielding depth. Although much work has been done on some of these dependencies, more detailed studies still need to be done on others. This work describes the influence of irradiation geometry on nuclide production in lunar rocks. In most cases, computer simulations of cosmogenic nuclide production were restricted to spherical objects irradiated with a 4 pi isotropic flux (meteoroids) or in lunar core samples irradiated by a 2 pi flux incident on semi-infinite layers or cylinders of huge sizes. Many lunar samples are rocks found on top of the lunar surface. For these rocks, neither of the above-mentioned models correspond to the real conditions. We present results of our simulations of cosmogenic nuclide production in models simulating the irradiation of rocks sitting on top of the lunar surface. The Galactic Cosmic Rays (GCR) production profiles in lunar rocks were calculated using the Los Alamos 3-D Monte Carlo LAHET Code System (LCS). The irradiated object was modeled as the union of a sphere with the radius of the Moon and a small hemisphere with radii varying from 10 to 100 g/sq cm simulating the lunar rock. These calculations for the production of cosmogenic nuclides in lunar rocks by GCR particle show that there are important differences between the results obtained by commonly used geometric irradiation models and the lunar-rock models presented. The steeper GCR production profiles for a rock could help to explain the poor agreement for Be-10 in rock 68815, where slab models give GCR profiles flatter than the observed profiles.

Reedy, R. C.↗

Polarization from scattering in blobs

We describe a Monte Carlo method for calculating polarization due to scattering in various geometries. We apply the method in this paper to the case of scattering in a spherical cloud illuminated by parallel rays. The polarization is largest in small optical depth blobs, where single scattering dominates. The polarization decreases as optical depth increases until a limit is reached where scattering occurs on the surface, and increasing optical depth has little effect on the results. For a given optical depth, decreasing the albedo leads to an increase of the degree of polarization of the scattered light because single scattering increases relative to multiple scattering. As the optical depth of the blob increases, the scattered flux becomes increasingly backward throwing. In the high albedo blobs, the maximum polarization becomes skewed towards the illuminated hemisphere at high optical depths. We discuss the applications to polarization in supergiants and R CrB stars, and scattering in a clumpy interstellar medium.

Code, Arthur D.↗

Astro-1 and ground-based observations of Markarian 335: Evidence for an accretion disk

Simultaneous UV and optical observations of the Seyfert galaxy Markarian 335 (z = 0.026) during the Astro-1 mission yield a spectrum spanning the wavelength range of 912-8410 A. In the sub-Ly alpha region a prominent blended emission feature of O VI lambda lambda 1032, 1038, and Ly beta is nearly as strong as C IV wavelength 1549. The continuum flux extends beyond the redshifted Luman limit without a noticeable discontinuity, but a siginificant change in slope exists near the redhsifted Lyman edge. We suggest that such a change may be the signature of a Lyman edge in an accretion disk seen at a small inclination angle. Using a disk model including such an edge, we fit the spectrum with a central black hole mass of 5 x 10(exp 7) solar mass, an accretion rate of 0.07 solar mass/yr, and an optical depth at the Lyman edge of 0.4. To account for the strong O VI emission as well as the soft X-ray excess, we consider the effects of Comptonization on the disk spectrum, which can produce a high-energy tail for the UV bump and also smooth the Lyman edge feature.

Zheng, W.↗

Estimation of 3-D Cloud Effects on TOMS Satellite Retrieval of Surface UV Irradiance

To improve surface UV irradiance retrieval from the Total Ozone Mapping Spectrometer (TOMS) we simulate errors of the TOMS cloud correction algorithm for summertime broken cloud conditions. Cloud scenes (50 km by 50 km) are modeled by a normal random (Gaussian) field with a fixed lower boundary and conservative scattering. The model relates stochastic field characteristics with the cloud amount, mean cloud diameter and aspect ratio. Clouds are embedded into Rayleigh atmosphere with standard ozone profile. Radiative transfer calculations of the radiance at the top of the atmosphere and irradiance at the surface were performed using 3-D Monte Carlo (MC) code. The results are averaged over the satellite field of view on the surface (50 km by 50 km) and compared with TOMS predicted surface irradiance for the same scene reflectance. The TOMS algorithm assumes horizontally homogeneous Cl-type cloud between 3 km and 5.5 km. The effective optical depth is determined by fitting observed (MC) radiance at 380 nm. Having the same radiance at the satellite the homogeneous and broken cloud models predict different average irradiances at the surface. This is due to the differences in Bidirectional Reflection Distribution Function (BRDF) for homogeneous and broken cloud scenes with the same hemispherical albedo. For typical TOMS observational geometry at mid-latitudes the simulated single pixels errors may be as large as +/- 20%. Qualitatively these errors are due to the dominance of the non-horizontal cloud surfaces, which are not accounted for in the homogeneous cloud model. However, due to high variability of the real cloud shapes and types it is unclear how these single pixel errors would affect TOMS time-integrated UV exposure over extended periods (weeks to months) for different regions.

Krotkov, Nickolay A.↗

Parametric Weight Comparison of Advanced Metallic, Ceramic Tile, and Ceramic Blanket Thermal Protection Systems

A parametric weight assessment of advanced metallic panel, ceramic blanket, and ceramic tile thermal protection systems (TPS) was conducted using an implicit, one-dimensional (I-D) finite element sizing code. This sizing code contained models to account for coatings fasteners, adhesives, and strain isolation pads. Atmospheric entry heating profiles for two vehicles, the Access to Space (ATS) vehicle and a proposed Reusable Launch Vehicle (RLV), were used to ensure that the trends were not unique to a certain trajectory. Ten TPS concepts were compared for a range of applied heat loads and substructural heat capacities to identify general trends. This study found the blanket TPS concepts have the lightest weights over the majority of their applicable ranges, and current technology ceramic tiles and metallic TPS concepts have similar weights. A proposed, state-of-the-art metallic system which uses a higher temperature alloy and efficient multilayer insulation was predicted to be significantly lighter than the ceramic tile stems and approaches blanket TPS weights for higher integrated heat loads.

Myers, David E.↗

The Effect of Heat Treatment on Residual Stress and Machining Distortions in Advanced Nickel Base Disk Alloys

This paper describes an extension of NASA's AST and IDPAT Programs which sought to predict the effect of stabilization heat treatments on residual stress and subsequent machining distortions in the advanced disk alloy, ME-209. Simple "pancake" forgings of ME-209 were produced and given four heat treats: 2075F(SUBSOLVUS)/OIL QUENCH/NO AGE; 2075F/OIL QUENCH/1400F@8HR;2075F/OIL QUENCH/1550F@3HR/l400F@8HR; and 2160F(SUPERSOLVUS)/OIL QUENCH/1550F@3HR/ 1400F@8HR. The forgings were then measured to obtain surface profiles in the heat treated condition. A simple machining plan consisting of face cuts from the top surface followed by measurements of the surface profile opposite the cut were made. This data provided warpage maps which were compared with analytical results. The analysis followed the IDPAT methodology and utilized a 2-D axisymmetric, viscoplastic FEA code. The analytical results accurately tracked the experimental data for each of the four heat treatments. The 1550F stabilization heat treatment was found to significantly reduce residual stresses and subsequent machining distortions for fine grain (subsolvus) ME209, while coarse grain (supersolvus) ME209 would require additional time or higher stabilization temperatures to attain the same degree of stress relief.

Gayda, John↗

Progress on the Plasmoid Thruster Experiment (PTX)

A plasmoid is a compact plasma structure with an integral magnetic field, that may be categorized according to the relative strength of the poloidal and toroidal magnetic field (B(sub p) and B(sub t), respectively). An object with B(sub p)/B(sub t), much much more than 1 is called a Field Reverse Configuration (FRC); if B(sub p) approximately equal to B(sub t), it is called a Spheromak. The thruster operates by repetitively producing plasmoids that are accelerated and ejected at high velocity. As this process is inductive, there are no electrodes. Also, the magnetic structure of the plasmoid should suppress thermal and mass losses to the wall, and improve detachment of the plasma exhaust from the thruster. This concept should be capable of producing an Isp in the range of 5,000 - 10,000 s with thrust densities of order 10(exp 5) N per square meters. The current experiment is designed to produce jet powers in the range of 5-10 kW, although the concept should be scalable into the MW range. In PTX, the plasmoid is formed inside of a single turn conical theta-pinch coil (17.58 cone angle). The coil is driven by a 640 nF, 35 kV capacitor bank, which rings at a frequency of 500 kHz. Previous experiments on PTX were conducted with a static-fill of propellant gas (6% H2 in He), and demonstrated reliable ionization over a pressure range of 40 - 200 mTorr. We are now adding a fast gas-puff valve to load the propellant, and a ringing pre-ionization circuit (f = 5 Mhz) to better control the plasmoid formation. An alternate coil (8.58 cone angle) will also be used, so as to investigate the effect of coil shape on performance. In addition, a variety of propellants will be used, including hydrogen, nitrogen, and argon. The plasmoid mass and velocity will be measured with a variety of diagnostics, including external B-dot probes and flux loops, a high-speed framing camera, and a HeNe laser interferometer. Internal B-dot probes and a quadruple Langmuir probe will provide additional measurements of the plasmoid shape and structure, as well as density, and will be described in a companion paper. The experimental results will be compared to calculations made using a 0-D coil-gun model and also a 2-D time dependent MHD code.

Martin, Adam↗

The Production of Jets From Magnetic Accretion Disks: Simulation of the Blandford-Payne Mechanism

We have performed magnetohydrodynamic (MRD) simulations of the production of jets from magnetized accretion disks with a factor of 5 greater extent in space and time, and with more models, than any study published so far. We find that jets are produced by such disks in a broad range of parameter space, and by at least two different mechanisms. We also are able to follow the propagation of the jet well beyond the accretion disk into the region of hydrodynamic collimation. The code used is our MHD simulation code FLOW (K. Lind, D. Payne, D. Meier, and R. Blandford, 1989), converted to run on Caltech's massively parallel Intel Touchstone Delta supercomputer. Some of these models may be directly applicable to observed radio sources.

Blandford-Payne↗

Recent Simulation Results on Ring Current Dynamics Using the Comprehensive Ring Current Model

Plasma sheet conditions and electromagnetic field configurations are both crucial in determining ring current evolution and connection to the ionosphere. In this presentation, we investigate how different conditions of plasma sheet distribution affect ring current properties. Results include comparative studies in 1) varying the radial distance of the plasma sheet boundary; 2) varying local time distribution of the source population; 3) varying the source spectra. Our results show that a source located farther away leads to a stronger ring current than a source that is closer to the Earth. Local time distribution of the source plays an important role in determining both the radial and azimuthal (local time) location of the ring current peak pressure. We found that post-midnight source locations generally lead to a stronger ring current. This finding is in agreement with Lavraud et al.. However, our results do not exhibit any simple dependence of the local time distribution of the peak ring current (within the lower energy range) on the local time distribution of the source, as suggested by Lavraud et al. [2008]. In addition, we will show how different specifications of the magnetic field in the simulation domain affect ring current dynamics in reference to the 20 November 2007 storm, which include initial results on coupling the CRCM with a three-dimensional (3-D) plasma force balance code to achieve self-consistency in the magnetic field.

Zheng, Yihua↗

Simplified Analysis of Pulse Detonation Rocket Engine B1owdown Gasdynamics and Performance

Pulsed detonation rocket engines (PDREs) have generated considerable research interest in recent years as a chemical propulsion system potentially offering improved performance and reduced complexity compared to conventional rocket engines. The detonative mode of combustion employed by these devices offers a thermodynamic advantage over the constant-pressure deflagrative combustion mode used in conventional rocket engines and gas turbines. However, while this theoretical advantage has spurred a great deal of interest in building PDRE devices, the unsteady blowdown process intrinsic to the PDRE has made realistic estimates of the actual propulsive performance problematic. The recent review article by Kailasanath highlights some of the difficulties in comparing the available experimental measurements with numerical models. The goal of this paper is to improve understanding of PDRE blowdown gasdynamics and performance issues through use of a simplified model that captures the essential features of the unsteady blowdown process, and yet remains computationally inexpensive. The PDRE system studied here is highly idealized, consisting of a constant-area detonation tube with one end closed and the other end open to the environment. The tube is prefilled with a gaseous propellant mixture with no initial velocity or outflow to the environment. The detonation is initiated instantaneously at the closed end of the device. Chapman-Jouguet (C-J) post-detonation gas conditions are calculated using the CET89 version of the NASA thermochemical code. The I-D, unsteady method of characteristics is used to calculate the flowfield following the detonation front. See the compressible flow texts by Thompson and Zucrow and Hoffman for details of this method. Parametric studies of the effect of mixture stoichiometry, fill temperature, and blowdown pressure ratio on performance are reported. A comparison of the performance of an idealized straight-tube PDRE with a conventional steady-state rocket engine is provided. The effect of constant-gamma and equilibrium chemistry assumptions is also examined. Additionally, in order to form an assessment of the accuracy of the model, the flowfield time history is compared to experimental data from Stanford University.

Morris, Christopher I.↗

Numerical Simulations of the Partially-Ionized Gas in a 100-A LaB6 Hollow Cathode

Numerical simulations of a hollow cathode with a LaB6 emitter operating at 100 A have been performed for the first time using the 2-D Orificed Cathode (OrCa2D) code. Results for a variety of plasma properties are presented and compared with laboratory measurements. The large size of the device permits peak electron number densities in the cathode interior that are lower than those established in the NSTAR hollow cathode, which operates with a 7.3x lower discharge current and 3.2x lower mass flow rate. Also, despite the higher discharge current in the LaB6 cathode, the maximum electron current density is lower, by 4.2x, than that in the NSTAR cathode due to the larger orifice size. Simulations and direct measurements show that at 12 sccm of xenon flow the peak emitter temperature is in the range of 1594-1630 C. It is also found that the conditions for the excitement of current-driven streaming instabilities and ion-acoustic turbulence (IAT) are satisfied in this cathode, similarly to what was found in the past in its smaller counterparts like the NSTAR cathode. Based on numerical simulations, it has long been argued that these instabilities may be responsible for the anomalously large ion energies that have been measured in these discharges as well as for the enhancement of the plasma resistivity. Confirmation of the presence of IAT in this cathode is presented for the first time in a companion paper.

robotic exploration↗

Numerical investigations of background pressure effects and channel erosion in the SPT-140 Hall thruster for the Psyche mission

Four SPT-140 Hall thrusters are the main components of the electric propulsion system for the Psyche mission. We employ the 2-D (r-z) code Hall2De to conduct numerical simulations of the SPT-140 with the goal of estimating erosion of the acceleration channel rings for the Psyche mission profile. The simulations are informed by direct measurements of the plasma conditions in the acceleration channel that were obtained using laser-induced fluorescence. We validate our simulation results with additional plasma measurements, wear test erosion rates, and performance data as a function of background pressure. We then make use of a series of numerical simulations to compute the progressive erosion of the channel walls, using as input parameters a profile for the Psyche mission that includes operation at 4.5, 2.5 and 1.0 kW with duration of all segments increased by a 50% margin. We find that the channel walls are not completely eroded at the end of the mission (plus 50% margin). We also observe a progressive decrease in thrust as the channel walls erode that is consistent with the wear test and flight data, with a rate of ~4.5 mN/kh for the first 2,500 h and constant thrust after that. The reasons for these thrust trends are explained.

Lenguito, Giovanni↗