Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical relativity”

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 523 records · Page 29

Modulation of energetic particle fluxes by a mixed mode of transverse and compressional waves

Modulation characteristics of particle fluxes in the presence of a mixed mode of compressional and transverse magnetic waves at hydromagnetic frequencies are investigated through kinetic perturbation of the distribution function. The magnetospheric medium where the particles are modulated contains both the magnetic and pressure gradients. The modulation features are found to be strongly dependent on the energy and pitch angle of the particles. Drifting particles can resonate with waves whose phase velocities are close to their drift velocities. When this occurs, the modulation amplitudes become significantly large and large phase shifts will occur. It is pointed out that resonance is important for particles with mid pitch angles (40-70 deg). The phase shift between the particle modulations and the magnetic field oscillations are strongly controlled by the combined effects of transverse and compressional wave components and/or the occurrence of drift resonance. Numerical calculations are performed using the dispersion relation of drift mirror Alfven waves as an example of waves with both compressional and transverse components.

Lin, C. S.↗

Optimization of fringe-type laser anemometers for turbine engine component testing

The fringe type laser anemometer is analyzed using the Cramer-Rao bound for the variance of the estimate of the Doppler frequency as a figure of merit. Mie scattering theory is used to calculate the Doppler signal wherein both the amplitude and phase of the scattered light are taken into account. The noise from wall scatter is calculated using the wall bidirectional reflectivity and the irradiance of the incident beams. A procedure is described to determine the optimum aperture mask for the probe volume located a given distance from a wall. The expected performance of counter type processors is also discussed in relation to the Cramer-Rao bound. Numerical examples are presented for a coaxial backscatter anemometer.

Seasholtz, R. G.↗

Optimization of fringe-type laser anemometers for turbine engine component testing

The fringe type laser anemometer is analyzed using the Cramer-Rao bound for the variance of the estimate of the Doppler frequency as a figure of merit. Mie scattering theory is used to calculate the Doppler signal wherein both the amplitude and phase of the scattered light are taken into account. The noise from wall scatter is calculated using the wall bidirectional reflectivity and the irradiance of the incident beams. A procedure is described to determine the optimum aperture mask for the probe volume located a given distance from a wall. The expected performance of counter type processors is also discussed in relation to the Cramer-Rao bound. Numerical examples are presented for a coaxial backscatter anemometer. Previously announced in STAR as N84-25019

Seasholtz, R. G.↗

Some features of auroral electric fields as seen in 2D numerical simulations

Results of 2D plasma simulations are presented and related to auroral observations. The formation of V-shaped potentials is studied with a 2 1/2 dimensional electrostatic particle-in-cell code for a magnetized plasma. It is shown that amplitudes for perpendicular electric fields are larger than for parallel electric fields, and for Te less than 100 eV, the amplitudes are comparable to the electric fields associated with the electrostatic shocks observed from the S3-3 satellite. The excitation of electrostatic ion-cyclotron EIC waves which occurs in the region below the parallel potential drop is discussed. In auroral plasmas EIC waves are observed above the V-shaped double layers in association with ion beams and field-aligned currents. The results also show that oppositely directed electric fields in the center and at the edges of the simulation region produce oppositely directed currents. Precipitating auroral ions in association with electron inverted-V events are seen by the DMSP-F6 satellite.

Thiemann, H.↗

NASA-Chinese Aeronautical Establishment (CAE) Symposium

Several topics relative to combustion research are discussed. A numerical study of combustion processes in afterburners; the modeling of turbulent, reactive flow; gas turbine research; modeling of dilution jet flow fields; and chemical shock tubes as tools for studying high-temperature chemical kinetics are among the topics covered.

Source record↗

Calculation of external-internal flow fields for mixed-compression inlets

Supersonic inlet flows with mixed external-internal compressions were computed using a combined implicit-explicit (Beam-Warming-Steger/MacCormack) method for solving the three-dimensional unsteady, compressible Navier-Stokes equations in conservation form. Numerical calculations were made of various flows related to such inlet operations as the shock-wave intersections, subsonic spillage around the cowl lip, and inlet started versus unstarted conditions. Some of the computed results were compared with wind tunnel data.

Chyu, W. J.↗

Calculation of external-internal flow fields for mixed-compression inlets

Supersonic inlet flows with mixed external-internal compressions were computed using a combined implicit-explicit (Beam-Warming-Steger/MacCormack) method for solving the three-dimensional unsteady, compressible Navier-Stokes equations in conservation form. Numerical calculations were made of various flows related to such inlet operations as the shock-wave intersections, subsonic spillage around the cowl lip, and inlet started versus unstarted conditions. Some of the computed results were compared with wind tunnel data.

Chyu, W. J.↗

Magnetohydrodynamic waves and particle acceleration in a coronal loop

A model is presented for the acceleration of electrons in a flaring solar coronal loop. It is suggested that the primary energy release event in the flare establishes a spectrum of MHD waves within the loop which accelerate electrons by wave-particle resonant interactions. The nature of the MHD oscillations are described and numerical and analytic solution of the dispersion relation are examined. It is found that the behavior of a particle depends on whether or not it is in resonance with the wave spectrum. Resonant particles follow chaotic trajectories and can diffuse in velocity space between wave effective potentials, thereby gaining energy.

De La Beaujardiere, Jean-Francois↗

Development of a defect stream function, law of the wall/wake method for compressible turbulent boundary layers

The method presented is designed to improve the accuracy and computational efficiency of existing numerical methods for the solution of flows with compressible turbulent boundary layers. A compressible defect stream function formulation of the governing equations assuming an arbitrary turbulence model is derived. This formulation is advantageous because it has a constrained zero-order approximation with respect to the wall shear stress and the tangential momentum equation has a first integral. Previous problems with this type of formulation near the wall are eliminated by using empirically based analytic expressions to define the flow near the wall. The van Driest law of the wall for velocity and the modified Crocco temperature-velocity relationship are used. The associated compressible law of the wake is determined and it extends the valid range of the analytical expressions beyond the logarithmic region of the boundary layer. The need for an inner-region eddy viscosity model is completely avoided. The near-wall analytic expressions are patched to numerically computed outer region solutions at a point determined during the computation. A new boundary condition on the normal derivative of the tangential velocity at the surface is presented; this condition replaces the no-slip condition and enables numerical integration to the surface with a relatively coarse grid using only an outer region turbulence model. The method was evaluated for incompressible and compressible equilibrium flows and was implemented into an existing Navier-Stokes code using the assumption of local equilibrium flow with respect to the patching. The method has proven to be accurate and efficient.

Wahls, Richard A.↗

A comparative study between experimental results and numerical predictions of multi-wall structural response to hypervelocity impact

Over the last three decades, multiwall structures have been analyzed extensively, primarily through experiment, as a means of increasing the protection afforded to spacecraft structure. However, as structural configurations become more varied, the number of tests required to characterize their response increases dramatically. As an alternative, numerical modeling of high-speed impact phenomena is often being used to predict the response of a variety of structural systems under impact loading conditions. This paper presents the results of a preliminary numerical/experimental investigation of the hypervelocity impact response of multiwall structures. The results of experimental high-speed impact tests are compared against the predictions of the HULL hydrodynamic computer code. It is shown that the hypervelocity impact response characteristics of a specific system cannot be accurately predicted from a limited number of HULL code impact simulations. However, if a wide range of impact loadings conditions are considered, then the ballistic limit curve of the system based on the entire series of numerical simulations can be used as a relatively accurate indication of actual system response.

Schonberg, William P.↗

A combined algorithm for minimum time slewing of flexible spacecraft

The use of Pontryagin's Maximum Principle for the large-angle slewing of large flexible structures usually results in the so-called two-point boundary-value problem (TPBVP), in which many requirements (e.g., minimum time, small flexible amplitude, and limited control powers, etc.) must be satisfied simultaneously. The successful solution of this problem depends largely on the use of an efficient numerical computational algorithm. There are many candidate algorithms available for his problem (e.g., quasilinearization, gradient, and shooting, etc.). In this paper, a proposed algorithm, which combines the quasilinearization method with a time shortening technique and a shooting method, is applied to the minimum-time, three-dimensional, and large-angle maneuver of flexible spacecraft, particularly the orbiting Spacecraft Control Laboratory Experiment (SCOLE) configuration. Theoretically, the nonlinear TPBVP can be solved only through the shooting method to find the 'exact' switching times for the bang-bang controls. However, computationally, a suitable guess for the missing initial costates is crucial because the convergence range of the unknown initial costates is usually narrow, especially for systems with high dimensions and when a multi-bang-bang control strategy is needed. On the other hand, the problems of near minimum time attitude maneuver of general rigid spacecraft and fast slewing of flexible spacecraft have been examined by the authors through a numerical approach based on the quasilinearization algorithm with a time shortening technique. Computational results have demonstrated its broad convergence range and insensitivity to initial costate choices. Consequently, a combined approach is naturally suggested here to solve the minimum time slewing problem. That is, in the computational process, the quasilinearization method is used first to obtain a near minimum time solution. Then, the acquired converged initial costates from the quasilinearization approach are transformed (tailored) to and used as the initial costate guess for starting the shooting method. Finally, the shooting method takes over the remaining calculations until the minimum-time solution converges. The nonlinear equations of motion of the SCOLE are formulated by using Lagrange's equations, with the mast modeled as a continuous beam subject to three-dimensional deformations. The numerical results will be presented and some related computational issues will also be discussed.

Bainum, P. M.↗

Wave-current interaction near the Gulf Stream during the surface wave dynamics experiment

This paper presents a case study on the wave-current interaction near the local curvature of a Gulf Stream meander. The wave data were obtained from in situ measurements by a pitch-roll discus buoy during the Surface Wave Dynamics Experiment (SWADE) conducted off Wallops Island, Virginia, from October 1990 to March 1991. Owing to the advection of the Gulf Stream by the semidiurnal tide, the discus buoy was alternately located outside and inside the Gulf Stream. The directional wave measurements from the buoy show the changes in wave direction, wave energy, and directional spreading when waves encountered the current in the Gulf Stream meanders. A wave refraction model, using the ray-tracing method with an estimated Gulf Stream velocity field and meandering condition, was used to simulate wave refraction patterns and to estimate wave parameters at relative locations corresponding to buoy measurements. The numerical simulation shows that a focusing zone of wave rays was formed near the boundary and behind the crest of a simulated Gulf Stream meander. The focusing of wave rays causes changes in wave direction, increases in wave energy, and decreases in wave directional spreading, which are in good agreement with the results from the buoy measurements.

Wang, David W.↗

The effects of Kelvin-Helmholtz instability on resonance absorption layers in coronal loops

One of the long-standing uncertainties in the wave-resonance theory of coronal heating is the stability of the resonance layer. The wave motions in the resonance layer produce highly localized shear flows which vary sinusoidally in time with the resonance period. This configuration is potentially susceptible to the Kelvin-Helmholtz instability (KHI), which can enhance small-scale structure and turbulent broadening of shear layers on relatively rapid ideal timescales. We have investigated numerically the response of a characteristic velocity profile, derived from resonance absorption models, to finite fluid perturbations comparable to photospheric fluctuations. We find that the KHI primarily should affect long (approximately greater than 6 x 10(exp 4) km) loops where higher velocity flows (M approximately greater than 0.2) exist in resonance layers of order 100 km wide. There, the Kelvin-Helmholtz growth time is comparable to or less than the resonance quarter-period, and the potentially stabilizing magnetic effects are not felt until the instability is well past the linear growth stage. Not only is the resonance layer broadened by the KHI, but also the convective energy transport out of the resonance layer is increased, thus adding to the efficiency of the wave-resonance heating process. In shorter loops, e.g., those in bright points and compact flares, the stabilization due to the magnetic field and the high resonance frequency inhibit the growth of the Kelvin-Helmholtz instability beyond a minimal level.

Karpen, Judith T.↗

Computing supernova collapse to neutron stars and black holes

We present a new numerical code for spherical hydrodynamics in general relativity. The code can handle gravitational collapse to a neutron star or to a black hole without the appearance of singularities. Moreover, the variables and equations in the code are very similar to those appearing in traditional Lagrangian supernova codes. Any such existing code can thus be easily adapted to treat collapse where the final fate is uncertain and may be either a neutron star or a black hole. The code is based on the formulation of Hernandez & Misner, in which retarded time is used as coordinate. This prevents the computational grid from penetrating inside any black hole that may form. We present the equations and a complete finite difference scheme for the adiabatic evolution of a fluid that obeys a gamma-law equation of state. We summarize the results of several testbed calculations performed to check our code. We also give the transformation of the analytic Oppenheimer-Snyder solution for homogeneous dust collapse to our coordinate system.

Baumgarte, Thomas W.↗

Modeling and numerical simulation of microwave pulse propagation in air breakdown environment

It is shown that ionization occurs wherever the field intensity of the pulse exceeds the local breakdown threshold field of the background air. The produced plasma then attenuates the pulse and gives rise to a tail erosion phenomenon that plays the primary role in limiting the energy transfer of the pulse from source to destination. A theoretical model describing the propagation of an intense microwave pulse in an air breakdown environment is developed that includes the possible focusing effect introduced by either using phase array antennas or the other arrangements. The self-consistent description of the propagation process is provided by a set of two modal equations. These include a continuity equation (Poynting's equation) for the energy density of the pulse and a rate equation of the electron density. A forward wave approximation is used to simplify Poynting's equation, and a semiempirical formula is used for the ionization frequency, nu(sub i). This frequency provides the coupling between the two modal equations, and is used to express the electron rate equation explicitly. In terms of the relevant parameters of the atmosphere, these two equations are normalized for numerical analysis of pulse propagation in the atmosphere. The dependencies of the propagation characteristics of the pulse on intensity, frequency, width, and shape of the pulse are determined. The numerical simulations lead to a useful empirical relation p(exp 3)w = alpha = constant, where p and w are the incident power and width of the pulse and alpha depends on the percentage of pulse energy transferred from the source point to a destined position. The density distribution of the pulse's self-generated plasma is also evaluated. The results also show that for ionization caused by a single unfocused microwave pulse transmitted upwards from the ground, the maximum electron density produced at, for example, 50 km altitude is limited by the tail erosion effect to below 10(exp 6) cm(exp -3). Repetitive pulse and focused beam approaches are also examined. Both approaches can increase the maximum electron density by no more than an order of magnitude. A scheme using two obliquely propagating pulses intersecting at the destined height, e.g. 50 km, is considered. It is shown that the electron density generated at the lowest intersecting position can easily reach a value of 6.6 x 10(exp 8) cm(exp -3), which is considered to be high enough for artificial ionospheric mirror (AIM) application.

Kim, J.↗

Aerodynamic, unsteady, kinetic and heat loss effects on the dynamics and structure of weakly-burning flames

The first objective of the program is to introduce the meritorious counterflow methodology in microgravity in order to quantify the steady and unsteady characteristics of weakly-burning premixed and diffusion flames for a wide variety of conditions including elevated pressures. Subsequently, through detailed modeling and comparisons with the experimental data, to provide physical insight into the elementary mechanisms controlling the flame response. The configuration offers good control over the parameters of interest and can be modelled closely. The knowledge which will be gained from the counterflow flames will be subsequently used to analyze near-limit phenomena related to other configurations by conducting detailed numerical simulations including multidimensional ones. Among the problems to be analyzed are the downward and upward propagation of near-limit flames in tubes and phenomena observed in spherical and cylindrical geometries.

Egolfopoulos, Fokion N.↗

A Systematic Kernel Function Procedure for Determining Aerodynamic Forces on Oscillating or Steady Finite Wings at Subsonic Speeds

Details are given of a numerical solution of the integral equation which relates oscillatory or steady lift and downwash distributions in subsonic flow. The procedure has been programmed for the IBM 704 electronic data processing machine and yields the pressure distribution and some of its integrated properties for a given Mach number and frequency and for several modes of oscillation in from 3 to 4 minutes, results of several applications are presented.

Watkins, Charles E.↗

Interface Shape and Growth Rate Analysis of Se/GaAs Bulk Crystals Grown in the NASA Crystal Growth Furnace (CGF)

Selenium-doped gallium arsenide, Se/GaAs, bulk crystals have been grown on earth using NASA's crystal growth furnace (CGF) in preparation for microgravity experimentation on the USML-2 spacelab mission. Peltier cooling pulses of 50 ms duration, 2040 A magnitude, and 0.0033 Hz frequency were used to successfully demark the melt-solid interface at known times during the crystal growth process. Post-growth characterization included interface shape measurement, growth rate calculation, and growth rate transient determinations. It was found that the interface shapes were always slightly concave into the solid. The curvature of the seeding interfaces was typically 1.5 mm for the 15 mm diameter samples. This was in agreement with the predicted interface shapes and positions relative to the furnace determined using a numerical model of the sample/ampoule/cartridge assembly (SACA).

Bly, J. M.↗