SEARCH · Search NASA

Results for “CARTESIAN COORDINATE”

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 109 records · Page 6

All sources NASA NTRS DOE OSTI

NASA NTRS1981

Two-boundary grid generation for the solution of the three dimensional compressible Navier-Stokes equations

A grid generation technique called the two boundary technique is developed and applied for the solution of the three dimensional Navier-Stokes equations. The Navier-Stokes equations are transformed from a cartesian coordinate system to a computational coordinate system, and the grid generation technique provides the Jacobian matrix describing the transformation. The two boundary technique is based on algebraically defining two distinct boundaries of a flow domain and the distribution of the grid is achieved by applying functions to the uniform computational grid which redistribute the computational independent variables and consequently concentrate or disperse the grid points in the physical domain. The Navier-Stokes equations are solved using a MacCormack time-split technique. Grids and supersonic laminar flow solutions are obtained for a family of three dimensional corners and two spike-nosed bodies.

Smith, R. E.↗

NASA NTRS1981

Three methods for performing Hankel transforms

Three methods for performing Hankel transforms with optical or digital processors are described. The first method is applicable when the input data is available in Cartesian (x-y) format and uses the close connection between generalized Hankel transform and the two dimensional Fourier transform in Cartesian coordinates. The second method is useful when the input data is in polar (r - theta) format and uses change of variables to perform the nth order Hankel transform as a correlation integral. The third method utilizes the von Neumann addition theorem for Bessel functions to extract the Hankel coefficients from a correlation between the radial part of the input and a Bessel function. Initial experimental results obtained for optical implementation of the first two methods are presented.

Athale, R. A.↗

NASA NTRS1981

Eruptive solar magnetic fields

The quasi-steady evolution of solar magnetic fields in response to gradual photospheric changes is considered, with particular attention given to the threshold of a sudden eruption in the solar atmosphere. The formal model of an evolving, force-free field dependent on two Cartesian coordinates is extended to a field which is not force free but in static equilibrium with plasma pressure and gravity. The basic physics is illustrated through the evolution of a loop-shaped electric current sheet enclosing a potential bipolar field with footpoints rooted in the photosphere. A free-boundary problem is posed and then solved for the equilibrium configuration of the current sheet in a hydrostatically supported isothermal atmosphere. As the footpoints move apart to spread a constant photospheric magnetic flux over a larger region, the equilibria available extend the field to increasing heights.

Low, B. C.↗

NASA NTRS1982

Bubble migration in a rotating, liquid-filled sphere

Results and analysis of ground-based experiments performed to aid in designing experiments on the behavior of bubbles in a rotating liquid body on board the Shuttle in free fall are presented. Spherical shells filled with silicone oil containing a small gas bubble were spun and filmed by high speed motion picture photography. The rotation of the shell and the trajectory of the bubble motion were recorded and the film was exposed to a motion analyzer connected to a keypunch. The analyzer measured Cartesian coordinates as well as angle, frame number, and rotation rate. Optical correction equations were employed to determine the apparent bubble trajectory relative to an inertial frame of reference. An analytical model for the bubble motion was defined, yielding predictions of velocity and position at different times. Rotation of the fluid container is concluded to aid in centering the bubbles.

Annamalai, P.↗

NASA NTRS1983

Application of second-order-accurate Total Variation Diminishing (TVD) schemes to the Euler equations in general geometries

A one-parameter family of explicit and implicit second-order-accurate, entropy satisfying, total variation diminishing (TVD) schemes was developed by Harten. These TVD schemes were the property of not generating spurious oscillations for one-dimensional nonlinear scalar hyperbolic conservation laws and constant coefficient hyperbolic systems. Application of these methods to one- and two-dimensional fluid flows containing shocks (in Cartesian coordinates) yields highly accurate nonoscillatory numerical solutions. The goal of this work is to expand these methods to the multidimensional Euler equations in generalized coordinate systems. Some numerical results of shock waves impinging on cylindrical bodies are compared with MacCormack's method.

Yee, H. C.↗

NASA NTRS1983

High resolution applications of the Osher upwind scheme for the Euler equations

The 'Osher' scheme was introduced by Osher (1981). It represents an upwind finite-difference method for hyperbolic systems of conservation laws, including the Euler equations. In studies conducted by Osher (1981) and Osher and Solomon (1982), the method was applied to the nonisentropic form of the Euler equations in one dimension and the isentropically restricted form in two spatial dimensions, both in Cartesian coordinates. Chakravarthy and Osher (1982) have shown an approach for extending the Osher scheme to the Euler equations written for general geometries, taking into account the use of mappings to arbitrary curvilinear coordinate systems. The present investigation is concerned with the high resolution extension of the Osher scheme to second-order accuracy. Results are presented for several example problems, giving attention to quasi-one-dimensional Laval nozzle flow, a one-dimensional shock tube problem, and supersonic flow over a cylinder.

Chakravarthy, S. R.↗

NASA NTRS1983

A method for solving the transonic full-potential equation for general configurations

A method is developed for solving the full-potential equation for two-dimensional and axisymmetric flow which retains the grid and boundary condition simplicity of the transonic small-disturbance codes. The method is based on a finite-volume formulation of the mass conservation equation in a Cartesian coordinate system, and is an extension of the method of Purvis and Burkhalter (1979). This finite-volume approach, combined with the simple boundary treatment, is shown to result in a highly robust method applicable to a wide range of geometries and flow conditions. The accuracy of the method is demonstrated for general geometries in two-dimensional and axisymmetric flows. The use of this method results in significant gains in convergence rate over the vertical-line over-relaxation scheme by incorporating an AF2-type algorithm (Ballhaus et al., 1978). It is suggested that the simplicity of this method shold allow a relatively easy extension to complex geometries in three-dimensional flows, and complex two-dimensional configurations such as multielement airfoils should be amenable to this method.

Wedan, B.↗

NASA NTRS1984

Tortuosity of lightning return stroke channels

Data obtained from photographs of lightning are presented on the tortuosity of return stroke channels. The data were obtained by making piecewise linear fits to the channels, and recording the cartesian coordinates of the ends of each linear segment. The mean change between ends of the segments was nearly zero in the horizontal direction and was about eight meters in the vertical direction. Histograms of these changes are presented. These data were used to create model lightning channels and to predict the electric fields radiated during return strokes. This was done using a computer generated random walk in which linear segments were placed end-to-end to form a piecewise linear representation of the channel. The computer selected random numbers for the ends of the segments assuming a normal distribution with the measured statistics. Once the channels were simulated, the electric fields radiated during a return stroke were predicted using a transmission line model on each segment. It was found that realistic channels are obtained with this procedure, but only if the model includes two scales of tortuosity: fine scale irregularities corresponding to the local channel tortuosity which are superimposed on large scale horizontal drifts. The two scales of tortuosity are also necessary to obtain agreement between the electric fields computed mathematically from the simulated channels and the electric fields radiated from real return strokes. Without large scale drifts, the computed electric fields do not have the undulations characteristics of the data.

Levine, D. M.↗

NASA NTRS1984

Three-dimensional transonic potential flow about complex 3-dimensional configurations

An analysis has been developed and a computer code written to predict three-dimensional subsonic or transonic potential flow fields about lifting or nonlifting configurations. Possible condfigurations include inlets, nacelles, nacelles with ground planes, S-ducts, turboprop nacelles, wings, and wing-pylon-nacelle combinations. The solution of the full partial differential equation for compressible potential flow written in terms of a velocity potential is obtained using finite differences, line relaxation, and multigrid. The analysis uses either a cylindrical or Cartesian coordinate system. The computational mesh is not body fitted. The analysis has been programmed in FORTRAN for both the CDC CYBER 203 and the CRAY-1 computers. Comparisons of computed results with experimental measurement are presented. Descriptions of the program input and output formats are included.

Reyhner, T. A.↗

NASA NTRS1984

Observable in the reversible reaction gamma + d yields n + p

Calculations are presented on the effect of the two-body charge and current density effects on the observables in the photodisintegration of a deuteron. The computations are performed in a cartesian coordinate system, initially for gamma rays plane-polarized in the x-direction. Unpolarized photons are treated with an incoherent superposition of transition amplitudes for gamma rays polarized in the x- and y- directions. Account is taken of recoil and interaction energies after an initial state is defined. Consideration is focused on the cross section, polarization, and asymmetry of the reaction at low and medium energy incident gamma rays.

Rustgi, M. L.↗

NASA NTRS1984

Finite-amplitude models of convection in the early mantle

Models of finite-amplitude, time-dependent mantle convection at Reynolds number ten million are calculated. The models are based on the Boussinesq approximation for convection at infinite Prandtl number and treat convection in a two-dimensional, Cartesian coordinate box with an aspect ratio of 1.4. The initial conditions consist of a static fluid with a purely conductive temperature profile, while the boundary conditions are taken to be free slip, with impenetrable insulating side walls. The initial conditions are found to result in an initial burst of convection that reverses the horizontally averaged temperature gradient, making the middle layers stable with respect to convection. With bottom heating and cooling of the top, two distinct layers of convective cells occur; with top cooling only, the upper layer alone convects.

Boss, A. P.↗

NASA NTRS1985

A direct-inverse method for transonic and separated flows about airfoils

A direct-inverse technique and computer program called TAMSEP that can be sued for the analysis of the flow about airfoils at subsonic and low transonic freestream velocities is presented. The method is based upon a direct-inverse nonconservative full potential inviscid method, a Thwaites laminar boundary layer technique, and the Barnwell turbulent momentum integral scheme; and it is formulated using Cartesian coordinates. Since the method utilizes inverse boundary conditions in regions of separated flow, it is suitable for predicing the flowfield about airfoils having trailing edge separated flow under high lift conditions. Comparisons with experimental data indicate that the method should be a useful tool for applied aerodynamic analyses.

Carlson, K. D.↗

NASA NTRS1985

Formulation and application of Russell's method

It is shown that the numerical technique of Russell's momentum approach can be derived by using Hamilton's principle and Vance's numerical scheme. It results in a set of first order differnce equations for solving the angular velocities. The numerical examples show that the method is reliable. The algorithm is modified next to perform the analysis of N-body systems with closed loop topology. To increase the formulation flexibility, the equations of motion are represented by using Cartesian coordinates and Lagrange multipliers. The algorithm consists of two parts, Vance's scheme and an unconstrained minimization. The Vance's scheme is used to find the angular velocities, and the unconstrained minimization is applied to provide the correct angular displacements. The proposed scheme is further extended to find the design sensitivity of an N-body system with closed loop configuration, and to carry out the design optimization as well. The numerical example of a small-scaled mechanical system is presented to verify the proposed formulation.

Hou, J. W.↗

NASA NTRS1986

Analytical investigation of the dynamics of tethered constellations in Earth orbit, phase 2

The development of a two dimensional analytical model that describes the dynamics of an n-mass vertical tethered system is reported. Two different approaches are described: in the first one the control quantities are the independent variables while in the second one the Cartesian coordinates of each mass expressed in the orbiting reference frame are the independent variables. The latter model was used in the 3-mass version to simulate the dynamics of the tethered system in applications involving the displacement of the middle mass along the tether. In particular, issues related to reproducing predetermined acceleration profiles and g-tuning are reported.

Lorenzini, E.↗

NASA NTRS1986

Aeolian processes aboard a space station: Saltation and particle trajectory analysis

The Carousel wind tunnel (CWT) proposed to study aeolian processes aboard a space station consists of two concentric rotating drums. The space between the two drums comprises the wind tunnel test section. Differential rates of rotation of the two drums would provide a wind velocity with respect to either drum surface. Preliminary results of measured velocity profiles made in a CWT prototype indicate that the wall bounded boundary layer profiles are suitable to simulate flat plate turbulent boundary layer flow. The two dimensional flat plate Cartesian coordinate equations of motion of a particle moving through the air are explained. In order to assess the suitability of CWT in the analysis of the trajectories of windblown particles, a series of calculations were conducted comparing cases for gravity with those of zero gravity. Results from the calculations demonstrate that a wind tunnel of the carousel design could be fabricted to operate in a space station environment and that experiments could be conducted which would yield significant results contributing to the understanding of the physics of particle dynamics.

White, B. R.↗

NASA NTRS1986

A direct-inverse method for the prediction of transonic and separated flows about airfoils at high angles of attack

A direct-inverse technique and computer program called TAMSEP that can be used for the analysis of the flow about airfoils at subsonic and low transonic freestream velocities is presented. The method is based upon a direct-inverse nonconservative full potential inviscid method, a Thwaites laminar boundary layer technique, and the Barnwell turbulent momentum integral scheme; and it is formulated using Cartesian coordinates. Since the method utilizes inverse boundary conditions in regions of separated flow, it is suitable for predicting the flowfield about airfoils having trailing edge separated flow under high lift conditions. Comparisons with experimental data indicate that the method should be a useful tool for applied aerodynamic analyses.

Carlson, L. A.↗

NASA NTRS1985

Implicit finite-difference methods for the Euler equations

The present paper is concerned with two-dimensional Euler equations and with schemes which are in use of the time of this writing. Most of the development presented carries over directly to three dimensions. The characteristics of the two-dimensional Euler equations in Cartesian coordinates are considered along with generalized curvilinear coordinate transformations, metric relations, invariants of the transformation, flux Jacobian matrices and eigensystems, numerical algorithms, flux split algorithms, implicit and explicit nonlinear control (smoothing), upwind differencing in supersonic regions, unsteady and steady-state computation, the diagonal form of implicit algorithm, metric differencing and invariants, boundary conditions, geometry and mesh generation, and sample solutions.

Pulliam, T. H.↗

NASA NTRS1986

Static current-sheet models of quiescent prominences

A particular class of theoretical models idealize the prominence to be a discrete flat electric-current sheet suspended vertically in a potential magnetic field. The weight of the prominence is supported by the Lorentz force in the current sheet. These models can be extended to have curved electric-current sheets and to vary three-dimensionally. The equation for force balance is 1 over 4 pi (del times B) times Bdel p- p9 z=zero. Using Cartesian coordinates we take, for simplicity, a uniform gravity with constant acceleration g in the direction -z. If we are interested not in the detailed internal structure of the prominence, but in the global magnetic configuration around the prominence, we may take prominence plasma to be cold. Consideration is given to how such equilibrium states can be constructed. To simplify the mathematical problem, suppose there is no electric current in the atmosphere except for the discrete currents in the cold prominence sheet. Let us take the plane z =0 to be the base of the atmosphere and restrict our attention to the domain z greater than 0. The task we have is to solve for a magnetic field which is everywhere potential except on some free surface S, subject to suit able to boundary conditions. The surface S is determined by requiring that it possesses a discrete electric current density such that the Lorentz force on it is everywhere vertically upward to balance the weight of the material m(S). Since the magnetic field is potential in the external atmosphere, the latter is decoupled from the magnetic field and its plane parallel hydrostatic pressure and density can be prescribed.

Wu, F.↗