Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical computation”

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 469 records · Page 26

Effect of convection on the microstructure of a lamellar eutectic growing with a stepped interface

A two dimensional model developed to study the influence of convective flow on the concentration field ahead of a growing lamellar eutectic is described when one phase projects out into the melt creating a stepped interface. The two dimensional convective flow field, which is periodic in the horizontal direction, was computed numerically using the software FLUENT. The velocity field generated due to the flow of melt over the steps was then incorporated into a finite difference scheme employed to solve the concentration field. The average interfacial composition was calculated and converted to lamellar spacing using the Jackson and Hunt minimum supercooling criterion. It was found that a stepped interface is more sensitive to convection than a planar interface.

Seth, Jayshree↗

Welding arc plasma physics

The problems of weld quality control and weld process dependability continue to be relevant issues in modern metal welding technology. These become especially important for NASA missions which may require the assembly or repair of larger orbiting platforms using automatic welding techniques. To extend present welding technologies for such applications, NASA/MSFC's Materials and Processes Lab is developing physical models of the arc welding process with the goal of providing both a basis for improved design of weld control systems, and a better understanding of how arc welding variables influence final weld properties. The physics of the plasma arc discharge is reasonably well established in terms of transport processes occurring in the arc column itself, although recourse to sophisticated numerical treatments is normally required to obtain quantitative results. Unfortunately the rigor of these numerical computations often obscures the physics of the underlying model due to its inherent complexity. In contrast, this work has focused on a relatively simple physical model of the arc discharge to describe the gross features observed in welding arcs. Emphasis was placed of deriving analytic expressions for the voltage along the arc axis as a function of known or measurable arc parameters. The model retains the essential physics for a straight polarity, diffusion dominated free burning arc in argon, with major simplifications of collisionless sheaths and simple energy balances at the electrodes.

Cain, Bruce L.↗

A finite element-boundary integral method for scattering and radiation by two- and three-dimensional structures

A review of a hybrid finite element-boundary integral formulation for scattering and radiation by two- and three-dimensional composite structures is presented. In contrast to other hybrid techniques involving the finite element method, the proposed one is in principle exact and can be implemented using a low O(N) storage. This is of particular importance for large scale applications and is a characteristic of the boundary chosen to terminate the finite element mesh, usually as close to the structure as possible. A certain class of these boundaries lead to convolutional boundary integrals which can be evaluated via the fast Fourier transform (FFT) without a need to generate a matrix; thus, retaining the O(N) storage requirement. The paper begins with a general description of the method. A number of two- and three-dimensional applications are then given, including numerical computations which demonstrate the method's accuracy, efficiency, and capability.

Jin, Jian-Ming↗

Pressure measurements in a low-density nozzle plume for code verification

Measurements of Pitot pressure were made in the exit plane and plume of a low-density, nitrogen nozzle flow. Two numerical computer codes were used to analyze the flow, including one based on continuum theory using the explicit MacCormack method, and the other on kinetic theory using the method of direct-simulation Monte Carlo (DSMC). The continuum analysis was carried to the nozzle exit plane and the results were compared to the measurements. The DSMC analysis was extended into the plume of the nozzle flow and the results were compared with measurements at the exit plane and axial stations 12, 24 and 36 mm into the near-field plume. Two experimental apparatus were used that differed in design and gave slightly different profiles of pressure measurements. The DSMC method compared well with the measurements from each apparatus at all axial stations and provided a more accurate prediction of the flow than the continuum method, verifying the validity of DSMC for such calculations.

Penko, Paul F.↗

Recent manufacturing advances for spiral bevel gears

The U.S. Army Aviation Systems Command (AVSCOM), through the Propulsion Directorate at NASA LRC, has recently sponsored projects to advance the manufacturing process for spiral bevel gears. This type of gear is a critical component in rotary-wing propulsion systems. Two successfully completed contracted projects are described. The first project addresses the automated inspection of spiral bevel gears through the use of coordinate measuring machines. The second project entails the computer-numerical-control (CNC) conversion of a spiral bevel gear grinding machine that is used for all aerospace spiral bevel gears. The results of these projects are described with regard to the savings effected in manufacturing time.

Handschuh, Robert F.↗

Arcjet plasma plume effect on a microwave reflector antenna

The effect of plasma produced by an arc-jet thruster on the performance of on-board reflector antennas was investigated using a configuration in which the arc-jet plume is modeled as an inhomogeneous scatterer with a smoothly varying refractive index. On the basis of laws of geometrical optics, used to calculate the wave propagation, the trajectory, phase, amplitude, and polarization of the high-frequency ray field are computed numerically. The experiments performed to validate the calculated results experimentally are described.

Birkner, B. W.↗

Dynamics of surface tension in microgravity environment

Time-dependent evolutions of the profile of free surface (bubble shapes) for a cylindrical container partially filled with a Newtonian fluid of constant density, rotating about its axis of symmetry, in low and microgravity environments, have been studied. Numerical computations of the dynamics of bubble shapes have been carried out (1) linear time-dependent functions of spin-up and spin-down in low and microgravity environments, (2) linear time-dependent functions of increasing and decreasing gravity environment at high and low rotating cylinder speeds, (3) time-dependent step functions of spin-up and spin-down in a low-gravity environment, and (4) sinusoidal function oscillation of the gravity environment at high and low rotating cylinder speeds.

Hung, R. J.↗

Three dimensional numerical investigation of gravitational and solutal effects in a cylindrical cell

Three-dimensional numerical computations of natural convection in a cylindrical ampoule (L/R = 4) are presented for a Ga-doped Ge melt. The ampoule is maintained at isothermal end conditions, and the g vector is oriented at an angle gamma with the cylinder axis. Detailed velocity, temperature, and concentration field distributions are presented for different gravity levels and a range of inclination angles (gamma between 0 and 180 deg). For terrestrial conditions, complex, multicellular flow is found to occur for gamma between 0 and 180 deg. The strong convection results in significant isotherm distortions and enhanced heat transfer. As the gravity level is reduced to 0.001 g(0), the overall convection strength decays but the 3D multicellular flow persists, causing appreciable mass transfer. For a further reduction in the gravity level to 0.00001 g(0), the system thermo/solutal characteristics are close to a purely diffusive behavior. Flow and thermal and solutal distributions are presented in the principal, orthogonal, and cross sectional planes, and the average system heat and mass transfer are calculated.

Ramachandran, N.↗

Pressure measurements in a low-density nozzle plume for code verification

Measurements of Pitot pressure were made in the exit plane and plume of a low-density, nitrogen nozzle flow. Two numerical computer codes were used to analyze the flow, including one based on continuum theory using the explicit MacCormack method, and the other on kinetic theory using the method of direct-simulation Monte Carlo (DSMC). The continuum analysis was carried to the nozzle exit plane and the results were compared to the measurements. The DSMC analysis was extended into the plume of the nozzle flow and the results were compared with measurements at the exit plane and axial stations 12, 24 and 36 mm into the near-field plume. Two experimental apparatus were used that differed in design and gave slightly different profiles of pressure measurements. The DSMC method compared well with the measurements from each apparatus at all axial stations and provided a more accurate prediction of the flow than the continuum method, verifying the validity of DSMC for such calculations.

Penko, Paul F.↗

Energy deposition by X-ray photoelectrons into interstellar molecular clouds

The deposition of energy by 100-10,000 eV electrons into molecular hydrogen and helium gas of cosmic abundances and fractional ionization ranging from 10 to the -12th to 0.1 is considered. In order to account for energy deposition by both discrete and continuous degradation processes, the loss spectrum, a form of the Spencer-Fano degradation spectrum particularly suited to describing electron energy deposition in a partially ionized gas is considered. The calculated excitation and ionization yields drop as the fractional ionization increases. The dependences of these yields on fractional ionization are fitted to a form suitable for swift numerical computation.

Voit, G. M.↗

Gust response analysis for cascades operating in nonuniform mean flows

The linearized unsteady aerodynamic response of a cascade of airfoils subjected to entropic, vortical, and acoustic gusts is analyzed. Field equations for the first-order unsteady perturbation flow are obtained by linearizing the full time-dependent mass, momentum, and energy conservation equations about a nonlinear, isentropic, and irrotational mean or steady flow. A splitting technique is then used to decompose the unsteady velocity field into irrotational and rotational parts leading to field equations for the unsteady entropy, rotational velocity, and irrotational velocity fluctuations that are coupled only sequentially. The entropic and rotational velocity fluctuations can be described in terms of the mean-flow drift and stream functions which can be computed numerically. The irrotational unsteady velocity is described by an inhomogeneous linearized potential equation which contains a source term that depends on the rotational velocity field. This equation is solved via a finite difference technique. Results are presented to indicate the status of the numerical solution procedure and to demonstrate the impact of blade geometry and mean blade loading on the aerodynamic response of cascades to vortical gust excitations. The analysis described leads to very efficient predictions of cascade unsteady aerodynamics phenomena making it useful for turbomachinery aeroelastic and aeroacoustic design applications.

Hall, Kenneth C.↗

Detonation of cryogenic gaseous hydrogen-oxygen mixtures

The accidental mixing and detonation of oxygen-hydrogen mixtures is a serious aerospace industry hazard. The detonation characteristics of cryogenic mixtures of gaseous hydrogen and oxygen have accordingly been studied at various initial pressures and equivalence ratios at 100 K, and the results compared with numerical computations of idealized, steady-state, one-dimensional Zeldovich-von Neumann-Doering detonation structure with detailed chemical reaction kinetics. The predictions thus obtained for critical tube diameter, on the basis of calculated reaction-zone thickness, are found to agree reasonably well with experimental data.

Plaster, M.↗

Hybrid accretion disks in active galactic nuclei. I - Structure and spectra

A unified treatment is presented of the two distinct states of vertically thin AGN accretion disks: a cool (about 10 to the 6th K) optically thick solution, and a hot (about 10 to the 9th K) optically thin solution. A generalized formalism and a new radiative cooling equation valid in both regimes are introduced. A new luminosity limit is found at which the hot and cool alpha solutions merge into a single solution of intermediate optical depth. Analytic solutions for the disk structure are given, and output spectra are computed numerically. This is used to demonstrate the prospect of fitting AGN broadband spectra containing both the UV bump as well as the hard X-ray and gamma-ray tail, using a single accretion disk model. Such models are found to make definite predictions about the observed spectrum, such as the relation between the hard X-ray spectral index, the UV-to-X-ray luminosity ratio, and a feature of about 1 MeV.

Wandel, Amri↗

An analysis of the Hubble Space Telescope fine guidance sensor fine lock mode

There are two guiding modes of the Hubble Space Telescope (HST) used for the acquisition of astronomical data by one of its six scientific instruments. The more precise one is called Fine Lock. Command and control problems in the onboard electronics has limited Fine Lock to brighter stars, V less than 13.0 mag, instead of fulfilling its goal of V = 14.5 mag. Consequently, the less precise guiding mode of Coarse Track (approximately 40 milli-arc seconds) has to be used fairly frequently. Indeed, almost half of the scientific observations to have been made with the HST will be compromised. The only realistic or extensive simulations of the Fine Lock guidance mode are reported. The theoretical analysis underlying the Monte Carlo experiments and the numerical computations clearly show both that the control electronics are severely under-engineered and how to adjust the various control parameters to successfully extend Fine Lock guiding performance back to V = 14.0 mag and sometimes beyond.

Taff, L. G.↗

Theoretical studies of lasers and converters

Discussed here is Doppler broadening and its effects upon the stimulated emission cross-section connecting an upper level with a lower level in an iodine laser. The stimulated emission cross section is given. The level transitions for the iodine laser are illustrated as are the relative intensities of these transitions. Also discussed is the Voigt profile, which considers the effect of both Doppler broadening and collision broadening upon the absorption line shape. The equations describing the Voigt profile were added to the continuous flow laser model laser simulation program. The results were compared with the standard absorption profile reported in an earlier study. There seems to be no advantage to using the Voigt profile as the laser power output is relatively insensitive to changes in the absorption cross section at the pressures being considered for a space laser. One disadvantage of using the Voigt profile is the excessive numerical computations required by the additional equations.

Heinbockel, John H.↗

Progress with multigrid schemes for hypersonic flow problems

Several multigrid schemes are considered for the numerical computation of viscous hypersonic flows. For each scheme, the basic solution algorithm uses upwind spatial discretization with explicit multistage time stepping. Two level versions of the various multigrid algorithms are applied to the two dimensional advection equation, and Fourier analysis is used to determine their damping properties. The capabilities of the multigrid methods are assessed by solving three different hypersonic flow problems. Some new multigrid schemes based on semicoarsening strategies are shown to be quite effective in relieving the stiffness caused by the high aspect ratio cells required to resolve high Reynolds number flows. These schemes exhibit good convergence rates for Reynolds numbers up to 200 x 10(exp 6) and Mach numbers up to 25.

Radespiel, R.↗

Planetary cratering mechanics

To obtain a quantitative understanding of the cratering process over a broad range of conditions, we have numerically computed the evolution of impact induced flow fields and calculated the time histories of the major measures of crater geometry (e.g., depth diameter, lip height ...) for variations in planetary gravity (0 to 10 exp 9 cm/sq seconds), material strength (0 to 140 kbar), thermodynamic properties, and impactor radius (0.05 to 5000 km). These results were fit into the framework of the scaling relations of Holsapple and Schmidt (1987). We describe the impact process in terms of four regimes: (1) penetration; (2) inertial; (3) terminal; and (4) relaxation.

Okeefe, John D.↗

Mutual impedance of parallel and perpendicular coplanar surface monopoles

One dimensional integral formulas are derived for mutual impedance of arbitrary size, coplanar, parallel, and perpendicular surface monopoles. The integrals in formulas are expressed as exponential integrals where possible. The mutual impedance expression for parallel monopoles is a summation of exponential integrals and one-dimensional integrals. For perpendicular monopoles, the mutual impedance is in closed form, containing exponential integrals only. The final expressions are in a form suitable for numerical computation. Since the expressions contain at most one-dimensional integrals, they can be utilized to reduce the matrix filling time in the moment method formulations, especially when inhomogeneous sectioning is preferred. Additionally, they can be used in rectangular surface patch modeling of conducting surfaces with edges which are at an angle to the surface patches, providing the angle is small. To this end, the expressions were utilized in the moment method analysis of linearly tapered slot antennas. Very good accuracy was obtained with a reduction in computer time.

Koksal, Adnan↗