SEARCH · Engineering Papers

Results for “numeric profiling”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,405 records · Page 78

All sources NASA NTRS DOE OSTI

NASA NTRS2013

Nuclear Thermal Rocket Simulation in NPSS

Four nuclear thermal rocket (NTR) models have been created in the Numerical Propulsion System Simulation (NPSS) framework. The models are divided into two categories. One set is based upon the ZrC-graphite composite fuel element and tie tube-style reactor developed during the Nuclear Engine for Rocket Vehicle Application (NERVA) project in the late 1960s and early 1970s. The other reactor set is based upon a W-UO2 ceramic- metallic (CERMET) fuel element. Within each category, a small and a large thrust engine are modeled. The small engine models utilize RL-10 turbomachinery performance maps and have a thrust of approximately 33.4 kN (7,500 lbf ). The large engine models utilize scaled RL-60 turbomachinery performance maps and have a thrust of approximately 111.2 kN (25,000 lbf ). Power deposition profiles for each reactor were obtained from a detailed Monte Carlo N-Particle (MCNP5) model of the reactor cores. Performance factors such as thermodynamic state points, thrust, specific impulse, reactor power level, and maximum fuel temperature are analyzed for each engine design.

Belair, Michael L.↗

NASA NTRS2012

Influence of Alternative Engine Concepts on LCTR2 Sizing and Mission Profile

The Large Civil Tiltrotor (LCTR) was developed as part of the NASA Heavy Lift Rotorcraft Systems Investigation in order to establish a consistent basis for evaluating the benefits of advanced technology for large tiltrotors. The concept has since evolved into the second-generation LCTR2, designed to carry 90 passengers for 1,000 nm at 300 knots, with vertical takeoff and landing. This paper examines the impact of advanced propulsion system concepts on LCTR2 sizing. Two concepts were studied: an advanced, single-speed engine with a conventional power turbine layout (Advanced Conventional Engine, or ACE), and a variable-speed power turbine engine (VSPT). The ACE is the lighter engine, but requires a multi-speed (shifting) gearbox, whereas the VSPT uses a lighter, fixed-ratio gearbox. The NASA Design and Analysis of Rotorcraft (NDARC) design code was used to study the trades between rotor and engine efficiency and weight. Rotor performance was determined by Comprehensive Analytical Model of Rotorcraft Aerodynamics and Dynamics (CAMRAD II), and engine performance was estimated with the Numerical Propulsion System Simulation (NPSS). Design trades for the ACE vs. VSPT are presented in terms of vehicle weight empty for variations in mission altitude and range; the effect of different One Engine Inoperative (OEI) criteria are also examined. Because of its strong effect on gearbox weight and on both rotor and engine efficiency, rotor speed was chosen as the reference design variable for comparing design trades. The two propulsion concepts had nearly identical vehicle weights and mission fuel consumption, and their relative advantages varied little with cruise altitude, mission range, or OEI criteria; high cruise altitude and low cruise tip speed were beneficial for both concepts.

Acree, C. W., Jr.↗

NASA NTRS2013

Nuclear Thermal Rocket Simulation in NPSS

Four nuclear thermal rocket (NTR) models have been created in the Numerical Propulsion System Simulation (NPSS) framework. The models are divided into two categories. One set is based upon the ZrC-graphite composite fuel element and tie tube-style reactor developed during the Nuclear Engine for Rocket Vehicle Application (NERVA) project in the late 1960s and early 1970s. The other reactor set is based upon a W-UO2 ceramic-metallic (CERMET) fuel element. Within each category, a small and a large thrust engine are modeled. The small engine models utilize RL-10 turbomachinery performance maps and have a thrust of approximately 33.4 kN (7,500 lbf ). The large engine models utilize scaled RL-60 turbomachinery performance maps and have a thrust of approximately 111.2 kN (25,000 lbf ). Power deposition profiles for each reactor were obtained from a detailed Monte Carlo N-Particle (MCNP5) model of the reactor cores. Performance factors such as thermodynamic state points, thrust, specific impulse, reactor power level, and maximum fuel temperature are analyzed for each engine design.

Belair, Michael L.↗

NASA NTRS1994

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.↗

NASA NTRS1988

On the two-dimensional modelling of horizontal CVD reactors and its limitations

In order to clarify under what conditions two-dimensional theoretical models can realistically account for the actual growth rate distribution in horizontal CVD reactors, a systematic numerical study was performed, and the results were compared with published experimental data for the MOCVD of gallium arsenide. Three models of increasing complexity were studied. It is shown that two-dimensional models can produce realistic predictions only for reactors with large width-to-height (aspect) ratios, that are operated at subcritical Rayleigh number,i.e., under dominant forced flow conditions. But even then, depending on the carrier gas, thermal (Soret) diffusion must be included in the model. Furthermore, velocity corrections for finite aspect ratios must be made, and buoyancy effects can be significant in the entrance region. On the other hand, the results are relatively insensitive to axial diffusion, temperature dependence of transport properties (versus averaged properties), and an assumed Poiseuille velocity profile (versus a rigorous solution of the velocity field). Yet, boundary layer approximations are prone to give unrealistic results at the low Reynolds numbers typical for CVD operations.

Ouazzani, Jalil↗

NASA NTRS1992

Completing the evolution of supernova remnants and their bubbles

The filling fraction of hot gas in the ISM is reexamined with new calculations of the very long term evolution of SNRs and their fossil hot bubbles. Results are presented of a 1D numerical solution of the evolution of an SNR in a homogeneous medium with a nonthermal pressure corresponding to a 5-micro-G magnetic field and density of 0.2/cu cm. Comparison is made with a control simulation having no magnetic field pressure. It is found that the evolutions, once they have become radiative, differ in several significant ways, while both differ appreciably from qualitative pictures presented in the past. Over most of the evolution of either case, the hot bubble in the interior occupies only a small fraction of the shocked volume, the remainder in a thick shell of slightly compressed material. Column densities and radial distributions of O VI, N V, C IV, and Si IV as well as examples of absorption profiles for their strong UV lines are presented.

Slavin, Jonathan D.↗

NASA NTRS1992

Effects of temperature-dependent molecular absorption coefficients on the thermal infrared remote sensing of the earth surface

The effect of temperature-dependent molecular absorption coefficients on thermal infrared spectral signatures measured from satellite sensors is investigated by comparing results from the atmospheric transmission and radiance codes LOWTRAN and MODTRAN and the accurate multiple scattering radiative transfer model ATRAD for different atmospheric profiles. The sensors considered include the operational NOAA AVHRR and two research instruments planned for NASA's Earth Observing System (EOS): MODIS-N (Moderate Resolution Imaging Spectrometer-Nadir-Mode) and ASTER (Advanced Spaceborne Thermal Emission and Reflection Radiometer). The difference in band transmittance is as large as 6 percent for some thermal bands within atmospheric windows and more than 30 percent near the edges of these atmospheric windows. The effect of temperature-dependent molecular absorption coefficients on satellite measurements of sea-surface temperature can exceed 0.6 K. Quantitative comparison and factor analysis indicate that more accurate measurements of molecular absorption coefficients and better radiative transfer simulation methods are needed to achieve SST accuracy of 0.3 K, as required for global numerical models of climate, and to develop land-surface temperature algorithms at the 1-K accuracy level.

Wan, Zhengming↗

NASA NTRS2013

Hybrid-PIC Modeling of a High-Voltage, High-Specific-Impulse Hall Thruster

The primary life-limiting mechanism of Hall thrusters is the sputter erosion of the discharge channel walls by high-energy propellant ions. Because of the difficulty involved in characterizing this erosion experimentally, many past efforts have focused on numerical modeling to predict erosion rates and thruster lifespan, but those analyses were limited to Hall thrusters operating in the 200-400V discharge voltage range. Thrusters operating at higher discharge voltages (V(sub d) >= 500 V) present an erosion environment that may differ greatly from that of the lower-voltage thrusters modeled in the past. In this work, HPHall, a well-established hybrid-PIC code, is used to simulate NASA's High-Voltage Hall Accelerator (HiVHAc) at discharge voltages of 300, 400, and 500V as a first step towards modeling the discharge channel erosion. It is found that the model accurately predicts the thruster performance at all operating conditions to within 6%. The model predicts a normalized plasma potential profile that is consistent between all three operating points, with the acceleration zone appearing in the same approximate location. The expected trend of increasing electron temperature with increasing discharge voltage is observed. An analysis of the discharge current oscillations shows that the model predicts oscillations that are much greater in amplitude than those measured experimentally at all operating points, suggesting that the differences in oscillation amplitude are not strongly associated with discharge voltage.

Hall Thrusters↗

NASA NTRS1988

Cosmic-ray-modified stellar winds. III - A numerical iterative approach

A numerical iterative method is used to determine the modification of a stellar wind flow with a termination shock by the galactic cosmic rays. A two-fluid model consisting of cosmic rays and thermal stellar wind gas is used in which the cosmic rays are coupled to the background flow via scattering with magnetohydrodynamic waves or irregularities. A polytropic model is used to describe the thermal stellar wind gas, and the cosmic-rays are modeled as a hot, low-density gas with negligible mass flux. The positive galactic cosmic-ray pressure gradient serves to brake the outflowing stellar wind gas, and the cosmic rays modify the location of the critical point of the wind, the location of the shock, the wind fluid velocity profile, and the thermal gas entropy constants on both sides of the shock. The transfer of energy to the cosmic rays results in an outward radial flux of cosmic-ray energy.

Ko, C. M.↗

NASA NTRS1994

Triggering of starbursts in galaxies by minor mergers

Using numerical simulation, we explore the triggering of starburst activity in disk galaxies which accrete low-mass dwarf companions. In response to the tidal perturbation of an infalling satellite, a disk galaxy develops a strong two-armed spiral pattern, which in turn drives large quantities of disk gas into its central regions. The global star formation rate stays constant during the early stages of an accretion, before rising rapidly by an order of magnitude when the central gas density becomes very large. The associated central starburst is quite compact. Models which include a bulge component in the disk galaxy show that the presence of a bulge can suppress the radial gas flow and limit the strength of the associated starburst, depending on the overall mass profile. The fact that such relatively common 'minor' mergers may trigger strong starburst activity suggests that many disk galaxies may have experienced starbursts at some point in their lifetime. Implications for galaxy evolution and formation are discussed.

Mihos, J. Christopher↗

NASA NTRS2010

Winds from Luminous Late-Type Stars: II. Broadband Frequency Distribution of Alfven Waves

We present the numerical simulations of winds from evolved giant stars using a fully non-linear, time dependent 2.5-dimensional magnetohydrodynamic (MHD) code. This study extends our previous fully non-linear MHD wind simulations to include a broadband frequency spectrum of Alfven waves that drive winds from red giant stars. We calculated four Alfven wind models that cover the whole range of Alfven wave frequency spectrum to characterize the role of freely propagated and reflected Alfven waves in the gravitationally stratified atmosphere of a late-type giant star. Our simulations demonstrate that, unlike linear Alfven wave-driven wind models, a stellar wind model based on plasma acceleration due to broadband non-linear Alfven waves, can consistently reproduce the wide range of observed radial velocity profiles of the winds, their terminal velocities and the observed mass loss rates. Comparison of the calculated mass loss rates with the empirically determined mass loss rate for alpha Tau suggests an anisotropic and time-dependent nature of stellar winds from evolved giants.

Airapetian, V.↗

NASA NTRS1981

Neutral winds and electric fields in the dust auroral oval. II - Theory and model

A two-dimensional numerical model of the thermosphere is applied to the auroral zone neutral wind, electric field, and plasma density data set, presented in an earlier paper. The model shows that the action of the Lorentz force can be responsible to a great extent for the large zonal velocities near the 150-km altitude. Model equations are described, an explanation of the use of the geophysical conditions is given, and model integrations are compared to the wind measurements. However, for the two-dimensional model to be effective, the atmosphere must not cross too many meridians of local time during the integration period, so that the background state should remain fairly uniform. It is concluded that the two-dimensional model cannot accurately explain the details of the wind profiles because of the three-dimensional character of the physical situation. Thus it is noted that the observed winds were part of a large-scale three-dimensional flow which is only weakly coupled to short-term variations in magnetospheric conditions.

Mikkelsen, I. S.↗

NASA NTRS1990

Theory of biaxial graded-index optical fiber

A biaxial graded-index fiber with a homogeneous cladding is studied. Two methods, wave equation and matrix differential equation, of formulating the problem and their respective solutions are discussed. For the wave equation formulation of the problem it is shown that for the case of a diagonal permittivity tensor the longitudinal electric and magnetic fields satisfy a pair of coupled second-order differential equations. Also, a generalized dispersion relation is derived in terms of the solutions for the longitudinal electric and magnetic fields. For the case of a step-index fiber, either isotropic or uniaxial, these differential equations can be solved exactly in terms of Bessel functions. For the cases of an istropic graded-index and a uniaxial graded-index fiber, a solution using the Wentzel, Krammers and Brillouin (WKB) approximation technique is shown. Results for some particular permittivity profiles are presented. Also the WKB solutions is compared with the vector solution found by Kurtz and Streifer. For the matrix formulation it is shown that the tangential components of the electric and magnetic fields satisfy a system of four first-order differential equations which can be conveniently written in matrix form. For the special case of meridional modes, the system of equations splits into two systems of two equations. A general iterative technique, asymptotic partitioning of systems of equations, for solving systems of differential equations is presented. As a simple example, Bessel's differential equation is written in matrix form and is solved using this asymptotic technique. Low order solutions for particular examples of a biaxial and uniaxial graded-index fiber are presented. Finally numerical results obtained using the asymptotic technique are presented for particular examples of isotropic and uniaxial step-index fibers and isotropic, uniaxial and biaxial graded-index fibers.

Kawalko, Stephen F.↗

NASA NTRS1987

A numerical study of three-dimensional vortex breakdown

A numerical simulation of bubble-type vortex breakdown using a unique discrete form of the full 3-D, unsteady incompressible Navier-Stokes equations was performed. The Navier-Stokes equations were written in a vorticity-velocity form and the physical problem was not restricted to axisymmetric flow. The problem was parametized on a Rossby- Reynolds-number basis. Utilization of this parameter duo was shown to dictate the form of the free-field boundary condition specification and allowed control of axial breakdown location within the computational domain. The structure of the breakdown bubble was studied through time evolution plots of planar projected velocity vectors as well as through plots of particle traces and vortex lines. These results compared favorably with previous experimental studies. In addition, profiles of all three velocity components are presented at various axial stations and a Fourier analysis was performed to identify the dominant circumferential modes. The dynamics of the breakdown process were studied through plots of axial variation of rate of change of integrated total energy and rate of change of integrated enstrophy, as well as through contour plots of velocity, vorticity and pressure.

Spall, Robert E.↗

NASA NTRS1994

Chaotic dynamics of corotating magnetospheric convection

The corotating plasma convection system of the Jovian magnetosphere is analyzed. The macroscopic (mhd) model introduced by Summers and Mu, (1992) that incorporates the effects of microdiffusion is extended by including previously neglected density effects. We reduce the governing partial differential equations to a third-order ordinary differential system by the Galerkin technique of mode truncation. We carry out such a severe truncation partly in the interests of tractability, and leave open the question of the efficacy of adding additional modes. Exhaustive numerical integrations are carried out to calculate the long-term solutions, and we discover that a rich array of plasma motions is possible, dependent on the value of the height-integrated ionospheric Pederson conductivity Sigma. If Sigma is less than a certain critical value Sigma(sub c), then plasma motion can be expected to be chaotic (or periodic), while if Sigma is greater than Sigma(sub c), then steady state convection is expected. In the former case, whether the plasma motion is chaotic or periodic (and, if periodic, the magnitude of the period) can be very sensitive to the value of Sigma. The value of Sigma(sub c), which is a function of a parameter q that occurs in the assumed form of the stationary radial profile (varies as L(exp -q) of the plasma mass per unit magnetic flux, lies well within the accepted range of values of Sigma for Jupiter, i.e. Sigma greater than or equal to 0.1 mho and less than or equal to 10 mho.

Summers, Danny↗

NASA NTRS2019

Single-mode Fiber and Few-Mode Fiber Photonic Lanterns Performance Evaluated for Use in a Scalable Real-Time Photon Counting Ground Receiver

Photonic lanterns provide an efficient way of coupling light from a single large-core fiber to multiple small-core fibers. This capability is of interest for space to ground communication applications. In these applications, the optical ground receivers require high-efficiency coupling from an atmospherically distorted focus spot to multiple fiber coupled single pixel super-conducting nanowire detectors. This paper will explore the use of photonic lanterns in a real-time ground receiver that is scalable and constructed with commercial parts. The number of small-core fibers that make a photonic lantern determines the number of spatial modes that they couple. For instance, lanterns made with n number of single-mode fibers can couple n number of spatial modes. Although the laser transmitted from a spacecraft originates as a Gaussian shape, the atmosphere distorts the beam profile by scattering energy into higher-order spatial modes. Therefore, if a ground receiver is sized for a target data rate with n number of detectors, the corresponding lantern made with single-mode fibers will couple n number of spatial modes. The energy of the transmitted beam scattered into spatial modes higher than n will be lost. This paper shows this loss may be reduced by making lanterns with few-mode fibers instead of single-mode fibers, increasing the number of spatial modes that can be coupled and therefore increasing the coupling efficiency to single pixel, single photon detectors. The free space to fiber coupling efficiency of these two types of photonic lanterns are compared over a range of the free-space coupling numerical apertures and mode field diameters. Results indicate the few mode fiber lantern has higher coupling efficiency for telescopes with longer focal lengths under higher turbulent conditions. Also presented is analysis of the jitter added to the system by the lanterns, showing the few-mode fiber photonic lantern adds more jitter than the single-mode fiber lantern, but less than a multimode fiber.

photon counting↗

NASA NTRSDate not supplied

Laminar to Turbulence Transition in Boundary Layers due to Tripping Devices

Direct numerical simulations (DNS)and linear stability analysis were performed to identify the mechanisms by which two-and three-dimensional trips promote early transition in subsonic boundary layers over a flat plate. Investigations were conducted for two boundary layer transition trip configurations: one is a two-dimensional rectangular-shaped rod and the other is a three-dimensional zigzag tape. Reynolds numbers based on the trip heights are 683 and 300 for the two-and three-dimensional trips, respectively. Simulations with the two-dimensional trip showed that the flow separates at the top of the trip and forms a long recirculation zone downstream. However, the flow remains steady in the simulation without any external disturbances. Linear stability analysis showed that the inflectional profile in the recirculation zone is strongly unstable to frequencies in the range of 100-300 Hz. Simulations with the external disturbances forced at the frequency of 200 Hz showed that disturbances grow strongly inside the recirculation zone as observed in the Klebanoff1 experiment.

Transition↗

NASA NTRS2008

Assessments of Hollow Cathode Wear in the Xenon Ion Propulsion System (XIPs(c)) by Numerical Analyses and Wear Tests

The standard approach presently followed by NASA to qualify electric propulsion for the required mission throughput has been based largely on life tests, which can be costly and time consuming. Revised electric propulsion lifequalification approaches are being formulated that combine analytical and/or computational methods with (shorter-duration) wear tests. As a model case, a wear test is being performed at JPL to assess the lifetime of the discharge hollow cathode in the Xenon Ion Propulsion System (XIPS(c)), a 25-cm ion engine developed by L-3 Communications Electron Technologies, Inc. for commercial applications. Wear and plasma data accumulated throughout this life-assessment program are being used to validate the existing 2-D hollow cathode code OrCa2D. We find that the OrCa2D steady-state solution predicts very well the time-averaged plasma data and the keeper voltage after 5500 hrs of operation in high-power mode. When the wave motion that occurs naturally in these devices is accounted for, based on an estimate of the maximum wave amplitude, the molybdenum-keeper erosion profile observed in the XIPS(c) discharge cathode is also reproduced within a factor of two of the observation. When the same model is applied to predict the erosion of a tantalum keeper we find that erosion is reduced by more than two orders of magnitude compared to the molybdenum keeper due the significantly lower sputtering yield of tantalum. A tantalum keeper would therefore allow keeper lifetimes that greatly exceed the present requirements for deep-space robotic missions considered by NASA. Moreover, such large reduction of the erosion renders the largest uncertainties in the models, which are associated with the wave amplitude estimates and the electron transport model, negligible.

plasma simulations↗