Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “coordinate system”

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 1,243 records · Page 69

Computational aeroacoustics on massively parallel computers

A computational fluid dynamics (CFD) technique is employed to solve aeroacoustics problems on massively parallel computers. The algorithm is based on a 4th order accurate central finite difference and 4th order accurate 4-stage Runge-Kutta time integration method that solves the 3-D full Navier-Stokes/Euler equations on a curvilinear coordinate system. The code has been developed to predict noise radiation from ducted fans. A specific application of this code is made here to the oscillating circular piston problem (Workshop Category 4, Problem 2) by solving the full Euler equations with nonreflecting boundary conditions in the axisymmetric mode. The features of the code that apply to this particular problem are described in this paper.

Oezyoeruek, Yusuf↗

Growth-phase thinning of the near-Earth current sheet during the CDAW 6 substorm

The thinning of the near-Earth current sheet during the growth phase of the Coordinated Data Analysis Workshop (CDAW) 6 magnetospheric substorm is studied. The expansion onset of the substorm occurred at 1054 UT, March 22, 1979. During the growth phase, two spacecraft, International Sun Earth Explorer (ISEE) 1 and ISEE 2, were within the current sheet approximately 13 R(sub E) from the Earth and obtained simultaneous high-resolution magnetic data at two points in the current sheet. Plasma data were also provided by the ISEE spacecraft and solar wind data by IMP 8. To facilitate the analysis, the GSM magnetic field data are transformed to a 'neutral sheet coordinate system' in which the new x axis is parallel to the average magnetic field above and below the neutral sheet and the new y axis lies in the GSM equatorial plane. A model based on the assumption that the current sheet is a time-invariant structure fails to predict neutral sheet crossing times. Consequently, the Harris sheet model, which allows one to remove the restriction of time invariancy, is used instead. It is found that during the growth phase, a model parameter corresponding to the thickness of the current sheet decreased exponentially from about 5 R(sub E) to 1 R(sub E) with a time constant of about 14 min. In addition, the ISEE 1 and ISEE 2 neutral sheet crossings after expansion onset indicate that the neutral sheet was moving upward at 7 km/s relative to the spacecraft. Since both crossings occurred in approximately 80 s, the current sheet thickness is estimated to be about 500 km. These results demonstrate that the near-Earth current sheet undergoes dramatic thinning during the substorm growth phase and expansion onset.

Sanny, Jeff↗

Interplanetary magnetic field dependency of stable Sun-aligned polar cap arcs

This is the first analysis, using a statistically significant data set, of the morphological dependence of the presence, orientation, and motion of stable sun-aligned polar cap arcs upon the vector interplanetary magnetic field (IMF). For the one winter season analyzed we had 1392 all-sky 630.0-nm images of 2-min resolution containing a total of 150 polar cap arcs, all with corresponding values of the IMF as measured by International Monitoring Platform (IMP) 8 or International Sun Earth Explorer (ISEE) 2. After demonstrating an unbiased data set with smooth normal distributions of events versus the dimensions of time, space, and IMF component, we examine IMF dependencies of the properties of the optical arcs. A well-defined dependence for B(sub z) is found for the presence/absence of stable Sun-aligned polar cap arcs. Consistent with previous statistical studies, the probability of observing polar cap aurora steadily increases for larger positive values of B(sub z), and linearly decreases when B(sub z) becomes more negative. The probability of observing Sun-aligned arcs within the polar cap is determined to vary sharply as a function of the arc location; arcs were observed 40% of the time on the dawnside and only 10% on the duskside. This implies an overall probability of at least 40% for the whole polar cap. 20% of the arcs were observed during 'southward IMF conditions,' but in fact under closer inspection were found to have been formed under northward IMF conditions; these 'residual' positive B(sub z) arcs ha d a delayed residence time in the polar cap of about what would be expected after a north to south transition of B(sub z). A firm dependence on B(sub y) is also found for both the orientation and the dawn-dusk direction of motion of the arcs. All the arcs are Sun-aligned to a first approximation, but present deviations from this orientation, depending primarily upon the location of the arc in corrected geomagnetic (CG) coordinates. The arcs populating the 06-12 and the 12-18 quadrants of the CG coordinate system point toward the cusp. The B(sub y) dependency of the arc alignment is consistent with a cusp displacement in local time according to the sign of B(sub y). We found that the arc direction of motion depended both on B(sub y) and the arc location within the polar cap. For a given value of B(sub y) two well-defined regions (or cells) exist. Within each cell the arcs move in the same direction toward the boundary between the cells. The arcs located in the duskside move dawnward; those in the dawnside move duskward. The relative size of these dusk and dawn regions (or cells) are controlled by the magnitude of B(sub y). This persistent dusk-dawn motion fo the polar cap arcs is interpreted in terms of newly open flux tubes entering the polar cap and exerting a displacement of the convective cells and the polar cap arcs that are embedded within them.

Valladares, C. E.↗

Statistical study of ULF wave occurrence in the dayside magnetosphere

Ultralow-frequency (ULF) waves are observed almost everywhere in the dayside magnetosphere. The mechanism by which these waves are generated and transformed in the dayside magnetosphere is still not understood. Here we report a statistical study of these waves based on magnetic field data from the International Sun-Earth Explorer 1 (ISEE 1) spacecraft. Data from the first traversal of the spacecraft through the entire dayside magnetosphere have been examined to determine the spatial distribution of wave occurrence. Successive 20-min segments of data were transformed to a field-aligned coordinate system. The parallel component was detrended and all three components of the field spectrally analyzed. Wave occurrence was defined by the presence of significant peaks in the power spectra. Wave events were categorized by three wave frequency bands: Pc 3 with T approximately 10-45 s; Pc 4 with T approximately 45-150 s; the short-period part of the Pc 5 wave band with T approximately 150-324 s. Properties of the spectral peaks were then entered into a data base. The data base was next sorted to determine the spatial occurrence pattern for the waves. Our results show that Pc 3 waves most frequently occur just outside synchronous orbit and are approximately centered on local noon. Pc 4 waves have a similar distribution with its peak further out. Pc 5 waves have high occurrence rate at the two flanks of the magnetosphere. Peaks in spectra obtained near the magnetopause are less clearly defined than those deeper in the magnetosphere.

Cao, M.↗

Implications of the BATSE data for a helocentric origin of gamma-ray bursts

Several possible models have been suggested to explain the observed distribution of gamma-ray bursts: heliocentric distributions such as the Oort cloud, large galactic halos, and cosmological models. We report here on an investigation into the implications of the Burst and Transient Source Experiment (BATSE) gamma-ray burst distribution (Meegan et al. 1992a) data on the possible helocentric origin of gamma-ray bursts. We find no statistically significant anisotropy in the angular distribution of the bursts in a Sun-referenced coordinate system; there is no dipole moment in the direction of the Sun, and no quardrupole moment associated with the ecliptic plane. We have employed direct analytic calculations and Monte Carlo simulations of sources in the Oort cloud to constrain possible helicentric burst distributions. These can produce distributions consistent with the observed angular isotropy, the meal value of V/V(sub max), and the observed C/C(sub min) distribution of BATSE, and provide limits to burst energy of a few times approximately 10(exp 27) ergs. However, the agreement of the heliocentric C/C(sub min) distributions with the BATSE data is attributable to the relatively limited sampling of strong, nearby bursts. These bursts are known from observation to be homogeneously distributed, yet the density of sources in the Oort cloud is not constant in this region. Integral number-intensity distributions from the Oort cloud for larger numbers of bursts cannot reproduce the known homogeneity of the strong bursts without modification to the computed cometary number density and are therefore unlikely explanations of the gamma-ray burst distribution.

Horack, J. M.↗

Evidence for Langmuir oscillations and a low density cavity in the Venus magnetotail

We report the discovery of Langmuir oscillations in a very low plasma density region in the Venus magnetotail. These waves are observed more often at 30 kHz, but also at 5.4 kHz indicating densities as low as 0.3 cu cm in the central tail lobe. The Langmuir probe on board the Pioneer Venus Orbiter cannot resolve such a low plasma density. We use the magnetic field strength and the assumption of total pressure balance to infer the electron temperature as a test of the Langmuir wave interpretation. By investigating the spatial distribution of this wave activity we find that the plasma cavity is ordered in a coordinate system defined by the interplanetary magnetic field and is found at either side of the central tail current sheet.

Ho, C. M.↗

An interpretation of a mysterious 3.0- to 4.6-kHz emission band observed on Voyager 2 near Neptune

A whistler mode interpretation is provided for the narrowband signal (f approx. 3 - 4.6 kHz, Delta f approx. 200 - 800 Hz) detected by the plasma wave instrument on Voyager 2 during its encounter with Neptune. Our analysis indicates that this signal may have been generated in a limited spatial region and that it propagated to other regions of the Neptunian magnetosphere in the nonducted whistler mode with wave normal vectors lying close to the whistler mode resonance cone. The observed frequency variation of the emission along the Voyager 2 trajectory is consistent with this interpretation. The source location is estimated to be near the magnetic equator at L approx. 4 and dipole longitude of 111 deg W (260 deg W longitude in Neptune coordinate system). The source frequency and bandwidth are estimated to be 3.6 kHz and 300 Hz, respectively. The waves most likely would have been generated by energetic electrons with 2- to 20-keV parallel energy via a gyroresonance mechanism. Our interpretation of the narrowband emissions places the following limits on the Neptunian thermal plasma density and temperature: (1) N(sub e, min) greater than 0.16 el/cu cm for 1.2 R(sub N) less than R less than 5 R(sub N), (2) N(sub e, max) = 597.5/cu cm at R - 1.3 R(sub N), (3) T(sub e, max) less than 500-1000 K at R approx. 5 R(sub N). It is also possible that the weak UV aurora observed near Neptune could have been caused by the precipitation of energetic particles by the narrowband emission as a result of wave particle interactions.

Sonwalkar, Vikas S.↗

Use of Fourier transforms for asynoptic mapping: Applications to the Upper Atmosphere Research Satellite microwave limb sounder

Fourier analysis has been applied to data obtained from limb viewing instruments on the Upper Atmosphere Research Satellite. A coordinate system rotation facilitates the efficient computation of Fourier transforms in the temporal and longitudinal domains. Fields such as ozone (O3), chlorine monoxide (ClO), temperature, and water vapor have been transformed by this process. The transforms have been inverted to provide maps of these quantities at selected times, providing a method of accurate time interpolation. Maps obtained by this process show evidence of both horizontal and vertical transport of important trace species such as O3 and ClO. An examination of the polar regions indicates that large-scale planetary variations are likely to play a significant role in transporting midstratospheric O3 into the polar regions. There is also evidence that downward transport occurs, providing a means of moving O3 into the polar vortex at lower altitudes. The transforms themselves show the structure and propagation characteristics of wave variations.

Elson, Lee S.↗

Average motion, structure and orientation of the distant magnetotail determined from remote sensing of the edge of the plasma sheet boundary layer with E greater than 35 keV ions

We study gradients of the energetic ion intesity observed at the edge of the plasma sheet boundary layer (PSBL) by the energetic ion anisotropy spectrometer (EPAS) on International Sun Earth Explorer 3 (ISEE 3). In particular, we have determined the velocity of the boundary relative to the spacecraft in the direction perpendicular to the tail axis and the angle which the boundary normal makes to the spacecraft spin axis for 1160 PSBL encounters at X(sub GSM) greater than -240 R(sub E). By asuming that, on average, the edge of the PSBL is parallel to the cross-tail current sheet, we are then able to determine a number of properties of the structure, orientation and motion of the deep geomagnetic tail. We conclude the following: (1) Most crossings of the edge of the PSBL are caused by transverse motuion of the entire tail induced by solar wind direction variations, although some are caused by reconfiguration of the tail due to geomagnetic activity. (2) The typical velocity of the PSBL (and hence of the tail) in the direction perpendicular to the tail axis is 50-85 km/s. (3) The average twist of the tail is near zero, with the edge of the PSBL (and by inference the cross-tail current sheet) lying parallel to the ecliptic plane (however, large twists are found in individual events and the distribution of twists is broad, with one standard deviation of approximately 50 deg. (4) The width of the distribution decreases with downtail distance. (5) The variation of the distributions with cross-tail position reveals that this decreas in width is most likely due to the edge of the PSBL being concave, or significantly flared at the tail flanks, in the near-Earth region. (6) During days on which the Interplanetary Magnetic Field (IMF) has 'away' sector structure, the north lobe of the trail is twisted on average towards dawn by 7.0 +/-2.4 deg. (7) During days on which the IMF has 'toward' sector structure, the north lobe is tilted towards duskby 3.8 +/- 2.3. (8) A subset of events for which IMP 8 solar wind data are available show that, for southward IMF BH(sub z) the tail has a mean twist of -12.3 +/- 5.0 deg for IMF B(sub Y) greater than 0 and 5.5 +/- 3.8 deg for IMF B(sub Y) less than 0 (positive twist angles correspond to a tilt of the northern lobe towards dusk). (9) For northward IMF B(sub z) the tail has a twist of -23.9 +/- 5.0 deg for IMF B(sub Y) greater than 0 and 13.4 +/- 6.0 deg for IMF B(sub y) less than 0. Hence the tail appears more twisted on average for the IMF B(sub Z) northward case. (10) The distribution of tail twist is wider for lower levels of geomagnetic activity, indicating that the tail is able to twist more at lower levels of activity. (11) The data set reveals no evident effect of the earth's dipole wobble; tail orientation appears to be controlled by the solar wind and IMF, such that the GSE coordinate system may be appropriate for the study of field and plasma structures in the distant tail region.

Owen, C. J.↗

Interactive method for planning constrained, fuel-optimal orbital proximity operations

An interactive graphical method for planning fuel-efficient rendezvous trajectories in the multi-spacecraft environment of the space station is presented. The method allows the operator to compose a multi-burn transfer trajectory between arbitrary initial chaser and target trajectories. The available task time of the mission is limited and the maneuver is subject to various operational constraints, such as departure, arrival, plume impingement and spatial constraints. The maneuvers are described in terms of the relataive motion experienced in a Space-Station centered coordinate system. The optimization method is based on the primer vector and its extension to non-optimal trajectories. The visual feedback of trajectory shapes, operational constraints, and optimization functions, provided by user-transparaent and continuously active background computations, allows the operator to make fast, iterative design changes which rapidly converge to fuel-efficient solutions. The optimization functions are presented. A variety of simple design examples has been presented to demonstrate the usefulness of the method. In many cases the addition of a properly positioned intermediate waypoint resulted in fuel savings of up to 30%. Furthermore, due to the counter-intuitive character of the optimization functions, most fuel-optimal solutions could not have been found without the aid of the optimization tools. Operating the system was found to be very easy, and did not require any previous in-depth knowledge of orbital dynamics or trajectory. The planning tool is an example of operator assisted optimization of nonlinear cost-functions.

Abramovitz, Adrian↗

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

Automated Hardware-Identification System

"Compressed symbology" emerging technology involving one- and two-dimensional arrays of surface depressions to form optically readable dots. Patterns more durable and denser than common bar codes. Convey identification data in binary form and read by optoelectric sensors. Computers and compressed-symbology engraving machines they control constitute subsystems of "paperless" hardware-tracking and -identification systems coordinating flows of both identifying information and identified parts themselves, along with ancillary information like work orders. Modifications of software expected to accelerate marking operations, eliminate need for trial or practice marking, and reduce incidence of errors.

Schramm, Harry F., Jr.↗

Numerical Studies of Boundary-Layer Receptivity

Direct numerical simulations (DNS) of the acoustic receptivity process on a semi-infinite flat plate with a modified-super-elliptic (MSE) leading edge are performed. The incompressible Navier-Stokes equations are solved in stream-function/vorticity form in a general curvilinear coordinate system. The steady basic-state solution is found by solving the governing equations using an alternating direction implicit (ADI) procedure which takes advantage of the parallelism present in line-splitting techniques. Time-harmonic oscillations of the farfield velocity are applied as unsteady boundary conditions to the unsteady disturbance equations. An efficient time-harmonic scheme is used to produce the disturbance solutions. Buffer-zone techniques have been applied to eliminate wave reflection from the outflow boundary. The spatial evolution of Tollmien-Schlichting (T-S) waves is analyzed and compared with experiment and theory. The effects of nose-radius, frequency, Reynolds number, angle of attack, and amplitude of the acoustic wave are investigated. This work is being performed in conjunction with the experiments at the Arizona State University Unsteady Wind Tunnel under the direction of Professor William Saric. The simulations are of the same configuration and parameters used in the wind-tunnel experiments.

Reed, Helen L.↗

Sines and Cosines. Part 1 of 3

Applying the concept of similarities, the mathematical principles of circular motion and sine and cosine waves are presented utilizing both film footage and computer animation in this 'Project Mathematics' series video. Concepts presented include: the symmetry of sine waves; the cosine (complementary sine) and cosine waves; the use of sines and cosines on coordinate systems; the relationship they have to each other; the definitions and uses of periodic waves, square waves, sawtooth waves; the Gibbs phenomena; the use of sines and cosines as ratios; and the terminology related to sines and cosines (frequency, overtone, octave, intensity, and amplitude).

Apostol, Tom M.↗

A comparison of design variables for control theory based airfoil optimization

This paper describes the implementation of optimization techniques based on control theory for airfoil design. In our previous work in the area it was shown that control theory could be employed to devise effective optimization procedures for two-dimensional profiles by using either the potential flow or the Euler equations with either a conformal mapping or a general coordinate system. We have also explored three-dimensional extensions of these formulations recently. The goal of our present work is to demonstrate the versatility of the control theory approach by designing airfoils using both Hicks-Henne functions and B-spline control points as design variables. The research also demonstrates that the parameterization of the design space is an open question in aerodynamic design.

Reuther, James↗

Thermal and Structural Analysis of a Hollow Core Space Shuttle Main Engine (SSME) Turbine Blade

The influence of primary and secondary orientations on the elastic response of a hollow core, (001)-oriented nickel base single-crystal superalloy turbine blade, was investigated under combined thermal and mechanical conditions. Finite element techniques is employed through MARC finite element code to conduct the analyses on a hollow core SSME turbine blade made out of PWA 1480 single crystal material. Primary orientation of the single crystal superalloy was varied in increments of 2 deg, from 0 to 10 deg, from the (001) direction. Two secondary orientations (0 and 45 deg) were considered with respect to the global coordinate system, as the primary orientation angle was varied. The stresses developed within the single crystal blade were determined for different orientations of the blade. The influence of angular offsets such as the single crystal's primary and secondary orientations and the loading conditions on the elastic stress response of the PWA 1480 hollow blade are summarized. The influence of he primary orientation angle, when constrained between the bounds considered, was not found to be as significant as the influence of the secondary orientation angle.

Abdul-Aziz, Ali↗

Development of a two-equation turbulence model for hypersonic flows: Evaluation of a low Reynolds number correction to the Kappa - epsilon two equation compressible turbulence model - Volume 1

The objective of the current research is the development of an improved k-epsilon two-equation compressible turbulence model for turbulent boundary layer flows experiencing strong viscous-inviscid interactions. The development of an improved model is important in the design of hypersonic vehicles such as the National Aerospace Plane (NASP) and the High Speed Civil Transport (HSCT). Improvements have been made to the low Reynolds number functions in the eddy viscosity and dissipation of solenoidal dissipation of the k-epsilon turbulence mode. These corrections offer easily applicable modifications that may be utilized for more complex geometries. The low Reynolds number corrections are functions of the turbulent Reynolds number and are therefore independent of the coordinate system. The proposed model offers advantages over some current models which are based upon the physical distance from the wall, that modify the constants of the standard model, or that make more corrections than are necessary to the governing equations. The code has been developed to solve the Favre averaged, boundary layer equations for mass, momentum, energy, turbulence kinetic energy, and dissipation of solenoidal dissipation using Keller's box scheme and the Newton spatial marching method. The code has been validated by removing the turbulent terms and comparing the solution with the Blasius solution, and by comparing the turbulent solution with an existing k-epsilon model code using wall function boundary conditions. Excellent agreement is seen between the computed solution and the Blasius solution, and between the two codes. The model has been tested for both subsonic and supersonic flat-plate turbulent boundary layer flow by comparing the computed skin friction with the Van Driest II theory and the experimental data of Weighardt; by comparing the transformed velocity profile with the data of Weighardt, and the Law of the Wall and the Law of the Wake; and by comparing the computed results of an adverse pressure gradient with the experimental data of Fernando and Smits. Good agreement is obtained with the experimental correlations for all flow conditions.

Knight, Doyle D.↗

Development of Finite Elements for Two-Dimensional Structural Analysis Using the Integrated Force Method

The Integrated Force Method has been developed in recent years for the analysis of structural mechanics problems. This method treats all independent internal forces as unknown variables that can be calculated by simultaneously imposing equations of equilibrium and compatibility conditions. In this paper a finite element library for analyzing two-dimensional problems by the Integrated Force Method is presented. Triangular- and quadrilateral-shaped elements capable of modeling arbitrary domain configurations are presented. The element equilibrium and flexibility matrices are derived by discretizing the expressions for potential and complementary energies, respectively. The displacement and stress fields within the finite elements are independently approximated. The displacement field is interpolated as it is in the standard displacement method, and the stress field is approximated by using complete polynomials of the correct order. A procedure that uses the definitions of stress components in terms of an Airy stress function is developed to derive the stress interpolation polynomials. Such derived stress fields identically satisfy the equations of equilibrium. Moreover, the resulting element matrices are insensitive to the orientation of local coordinate systems. A method is devised to calculate the number of rigid body modes, and the present elements are shown to be free of spurious zero-energy modes. A number of example problems are solved by using the present library, and the results are compared with corresponding analytical solutions and with results from the standard displacement finite element method. The Integrated Force Method not only gives results that agree well with analytical and displacement method results but also outperforms the displacement method in stress calculations.

Kaljevic, Igor↗