Search NASA⌕ Search

SEARCH · Search NASA

Results for “spherical coordinates”

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 19 records
All sourcesNASA NTRSDOE OSTI
NASA NTRS2015

Spherical Coordinate Systems for Streamlining Suited Mobility Analysis

Introduction: When describing human motion, biomechanists generally report joint angles in terms of Euler angle rotation sequences. However, there are known limitations in using this method to describe complex motions such as the shoulder joint during a baseball pitch. Euler angle notation uses a series of three rotations about an axis where each rotation is dependent upon the preceding rotation. As such, the Euler angles need to be regarded as a set to get accurate angle information. Unfortunately, it is often difficult to visualize and understand these complex motion representations. It has been shown that using a spherical coordinate system allows Anthropometry and Biomechanics Facility (ABF) personnel to increase their ability to transmit important human mobility data to engineers, in a format that is readily understandable and directly translatable to their design efforts. Objectives: The goal of this project was to use innovative analysis and visualization techniques to aid in the examination and comprehension of complex motions. Methods: This project consisted of a series of small sub‐projects, meant to validate and verify a new method before it was implemented in the ABF's data analysis practices. A mechanical test rig was built and tracked in 3D using an optical motion capture system. Its position and orientation were reported in both Euler and spherical reference systems. In the second phase of the project, the ABF estimated the error inherent in a spherical coordinate system, and evaluated how this error would vary within the reference frame. This stage also involved expanding a kinematic model of the shoulder to include the rest of the joints of the body. The third stage of the project involved creating visualization methods to assist in interpreting motion in a spherical frame. These visualization methods will be incorporated in a tool to evaluate a database of suited mobility data, which is currently in development. Results: Initial results demonstrated that a spherical coordinate system is helpful in describing and visualizing the motion of a space suit. The system is particularly useful in describing the motion of the shoulder, where multiple degrees of freedom can lead to very complex motion paths.

Benson, Elizabeth↗
NASA NTRS1987

A simplified derivation of leaf normal spherical coordinates

The paper presents formulas for the direction in spherical coordinates of the normal vector to a planar leaf in terms of the coordinates of incident and reflected rays. A simple idea from vector analysis is applied. These formulas are important in plant canopy modeling and experimental data collection in situations where specular reflection at optical wavelengths plays a part.

Card, Don H.↗
NASA NTRS2014

Spherical Coordinate Systems for Streamlining Suited Mobility Analysis

When describing human motion, biomechanists generally report joint angles in terms of Euler angle rotation sequences. However, there are known limitations in using this method to describe complex motions such as the shoulder joint during a baseball pitch. Euler angle notation uses a series of three rotations about an axis where each rotation is dependent upon the preceding rotation. As such, the Euler angles need to be regarded as a set to get accurate angle information. Unfortunately, it is often difficult to visualize and understand these complex motion representations. One of our key functions is to help design engineers understand how a human will perform with new designs and all too often traditional use of Euler rotations becomes as much of a hindrance as a help. It is believed that using a spherical coordinate system will allow ABF personnel to more quickly and easily transmit important mobility data to engineers, in a format that is readily understandable and directly translatable to their design efforts. Objectives: The goal of this project is to establish new analysis and visualization techniques to aid in the examination and comprehension of complex motions. Methods: This project consisted of a series of small sub‐projects, meant to validate and verify the method before it was implemented in the ABF's data analysis practices. The first stage was a proof of concept, where a mechanical test rig was built and instrumented with an inclinometer, so that its angle from horizontal was known. The test rig was tracked in 3D using an optical motion capture system, and its position and orientation were reported in both Euler and spherical reference systems. The rig was meant to simulate flexion/extension, transverse rotation and abduction/adduction of the human shoulder, but without the variability inherent in human motion. In the second phase of the project, the ABF estimated the error inherent in a spherical coordinate system, and evaluated how this error would vary within the reference frame. This stage also involved expanding a kinematic model of the shoulder, to include the torso, knees, ankle, elbows, wrists and neck. Part of this update included adding a representation of 'roll' about an axis, for upper arm and lower leg rotations. The third stage of the project involved creating visualization methods to assist in interpreting motion in a spherical frame. This visualization method will be incorporated in a tool to evaluate a database of suited mobility data, which is currently in development.

Benson, Elizabeth↗
NASA NTRS1993

Density-conserving shape factors for particle simulations in cylindrical and spherical coordinates

It is established that a particle-to-grid weighting which is based on both conventional particle-in-cell (PIC) and cloud-in-cell (CIC) shape factors will yield nonuniform grid densities, even in the case of uniform particle distributions, in both cylindrical and spherical coordinates. In view of this, alternative density-conserving weighting schemes are discussed which encompass modified PIC and CIC weighting. These alternative shape factors are mathematically no more complicated than the conventional nondensity-conserving expressions.

Ruyten, Wilhelmus M.↗
NASA NTRS1996

Marching Cubes in Cylindrical and Spherical Coordinates

Isosurface extraction is a common analysis and visualization technique for three-dimensional scalar data. Marching Cubes is the most commonly-used algorithm for finding polygonal representations of isosurfaces in such data. We extend Marching Cubes to produce geometry for data sets that lie in spherical and cylindrical coordinate systems as well as show the steps for derivation of transformations for other coordinate systems.

techniques marching cubes isosurfaces coordinate s↗
NASA NTRS1992

Protostellar hydrodynamics: Constructing and testing a spatially and temporally second-order-accurate method. I - Spherical coordinates

Two related numerical schemes for calculating the 3D collapse of protostellar clouds are defined, developed, and checked on a wide variety of test problems in spherical symmetry and multiple dimensions. One scheme is first-order accurate in time (code S), and the other second-order accurate in time (code ST). Through convergence testing, the codes are shown to be second-order accurate in spatial differences. Compared with the previous 3D code, the combination of reduced numerical dissipation through second-order accuracy and of removing the systematic bias toward central concentrations implies that the tendency for fragmentation into binary or multiple protostars should increase. A reinvestigation of fragmentation as a mechanism for forming binary stars is expected to yield an even more favorable evaluation.

Boss, Alan P.↗
NASA NTRS1980

Errors incurred in a plane-wave-type expansion of a Gaussian beam

The multipole expansion obtained by Morita et al. (1968) of the Gaussian laser beam used to levitate an aerosol particle in order that its complete phase matrix may be measured is compared with that of Tsai and Pogorzelski (1975) in order to demonstrate the effect of the incorrect expansion used by Morita. Errors incurred by the use of an equation in which one side satisfies the scalar wave equation while the other side does not and can be reduced to a plane wave amplitude are calculated as functions of the inverse of the wave number times the beam waist, the wave number times the radial spherical coordinate and the angular spherical coordinate. Errors on the order of a few percent, considered undetectable are obtained in the squared-field amplitudes due to the expansion, however, they are found to become significant (several tens of percent) when the angle is zero. It is concluded that the expansion of Morita should only be used in the regions where the spherical angle is less than 0.01 and its product with the wave number and the radial spherical coordinate is less than unity.

Kattawar, G. W.↗
NASA NTRS1995

Radiation Pressure, Poynting-Robertson Drag, and Solar Wind Drag in the Restricted Three-Body Problem

In this paper, we examine the effects of radiation pressure, Poynting-Robertson (PR) drag, and solar wind drag on dust grains trapped in mean motion resonances with the Sun and Jupiter in the restricted (negligible dust mass) three-body Problem. We especially examine the evolution of dust grains in the 1:1 resonance. As a first step, the Sun and Jupiter are idealized to both be in circular orbit about a common center of mass (circular restricted three-body problem). From the equation of motion of the dust particle in the rotating reference frame, the drag-induced time rate of change of its Jacobi "constant," C, is then derived and expressed in spherical coordinates. This new mathematical expression in spherical coordinates shows that C, in the 1:1 resonance, both oscillates and secularly increases with increasing time. The new expression gives rise to an easy understanding of how an orbit evolves when the radiation force and solar wind drag are included. All dust grain orbits are unstable in time when PR and solar wind drag are included in the Sun-Jupiter-dust system. Tadpole orbits evolve into horseshoe orbits; and these orbits continuously expand in size to lead to close encounters with Jupiter. Permanent trapping is impossible. Orbital evolutions of a dust grain trapped in the 1:1 resonance in the planar circular, an inclined case, an eccentric case, and the actual Sun-Jupiter case are numerically simulated and compared with each other and show grossly similar time behavior. Resonances other than 1:1 are also explored with the new expression. Stable exterior resonance trapping may be possible under certain conditions. One necessary condition for such a trap is derived. Trapping in interior resonances is shown to be always unstable.

Liou, Jer-Chyi↗
NASA NTRS2001

Simple Analytic Expressions for the Magnetic Field of a Circular Current Loop

Analytic expressions for the magnetic induction (magnetic flux density, B) of a simple planar circular current loop have been published in Cartesian and cylindrical coordinates [1,2], and are also known implicitly in spherical coordinates [3]. In this paper, we present explicit analytic expressions for B and its spatial derivatives in Cartesian, cylindrical, and spherical coordinates for a filamentary current loop. These results were obtained with extensive use of Mathematica "TM" and are exact throughout all space outside of the conductor. The field expressions reduce to the well-known limiting cases and satisfy V · B = 0 and V x B = 0 outside the conductor. These results are general and applicable to any model using filamentary circular current loops. Solenoids of arbitrary size may be easily modeled by approximating the total magnetic induction as the sum of those for the individual loops. The inclusion of the spatial derivatives expands their utility to magnetohydrodynamics where the derivatives are required. The equations can be coded into any high-level programming language. It is necessary to numerically evaluate complete elliptic integrals of the first and second kind, but this capability is now available with most programming packages.

Simpson, James C.↗
NASA NTRS2001

A Highly Accurate Technique for the Treatment of Flow Equations at the Polar Axis in Cylindrical Coordinates using Series Expansions

Numerical methods for solving the flow equations in cylindrical or spherical coordinates should be able to capture the behavior of the exact solution near the regions where the particular form of the governing equations is singular. In this work we focus on the treatment of these numerical singularities for finite-differences methods by reinterpreting the regularity conditions developed in the context of pseudo-spectral methods. A generally applicable numerical method for treating the singularities present at the polar axis, when nonaxisymmetric flows are solved in cylindrical, coordinates using highly accurate finite differences schemes (e.g., Pade schemes) on non-staggered grids, is presented. Governing equations for the flow at the polar axis are derived using series expansions near r=0. The only information needed to calculate the coefficients in these equations are the values of the flow variables and their radial derivatives at the previous iteration (or time) level. These derivatives, which are multi-valued at the polar axis, are calculated without dropping the accuracy of the numerical method using a mapping of the flow domain from (0,R)*(0,2pi) to (-R,R)*(0,pi), where R is the radius of the computational domain. This allows the radial derivatives to be evaluated using high-order differencing schemes (e.g., compact schemes) at points located on the polar axis. The proposed technique is illustrated by results from simulations of laminar-forced jets and turbulent compressible jets using large eddy simulation (LES) methods. In term of the general robustness of the numerical method and smoothness of the solution close to the polar axis, the present results compare very favorably to similar calculations in which the equations are solved in Cartesian coordinates at the polar axis, or in which the singularity is removed by employing a staggered mesh in the radial direction without a mesh point at r=0, following the method proposed recently by Mohseni and Colonius (1). Extension of the method described here for incompressible flows or for any other set of equations that are solved on a non-staggered mesh in cylindrical or spherical coordinates with finite-differences schemes of various level of accuracy is immediate.

Constantinescu, George S.↗
NASA NTRS1988

A numerical study of viscous vortex rings using a spectral method

Viscous, axisymmetric vortex rings are investigated numerically by solving the incompressible Navier-Stokes equations using a spectral method designed for this type of flow. The results presented are axisymmetric, but the method is developed to be naturally extended to three dimensions. The spectral method relies on divergence-free basis functions. The basis functions are formed in spherical coordinates using Vector Spherical Harmonics in the angular directions, and Jacobi polynomials together with a mapping in the radial direction. Simulations are performed of a single ring over a wide range of Reynolds numbers (Re approximately equal gamma/nu), 0.001 less than or equal to 1000, and of two interacting rings. At large times, regardless of the early history of the vortex ring, it is observed that the flow approaches a Stokes solution that depends only on the total hydrodynamic impulse, which is conserved for all time. At small times, from an infinitely thin ring, the propagation speeds of vortex rings of varying Re are computed and comparisons are made with the asymptotic theory by Saffman. The results are in agreement with the theory; furthermore, the error is found to be smaller than Saffman's own estimate by a factor square root ((nu x t)/R squared) (at least for Re=0). The error also decreases with increasing Re at fixed core-to-ring radius ratio, and appears to be independent of Re as Re approaches infinity). Following a single ring, with Re=500, the vorticity contours indicate shedding of vorticity into the wake and a settling of an initially circular core to a more elliptical shape, similar to Norbury's steady inviscid vortices. Finally, we consider the case of leapfrogging vortex rings with Re=1000. The results show severe straining of the inner vortex core in the first pass and merging of the two cores during the second pass.

Stanaway, S. K.↗
NASA NTRS1983

The NCOREL computer program for 3D nonlinear supersonic potential flow computations

An innovative computational technique (NCOREL) was established for the treatment of three dimensional supersonic flows. The method is nonlinear in that it solves the nonconservative finite difference analog of the full potential equation and can predict the formation of supercritical cross flow regions, embedded and bow shocks. The method implicitly computes a conical flow at the apex (R = 0) of a spherical coordinate system and uses a fully implicit marching technique to obtain three dimensional cross flow solutions. This implies that the radial Mach number must remain supersonic. The cross flow solutions are obtained by using type dependent transonic relaxation techniques with the type dependency linked to the character of the cross flow velocity (i.e., subsonic/supersonic). The spherical coordinate system and marching on spherical surfaces is ideally suited to the computation of wing flows at low supersonic Mach numbers due to the elimination of the subsonic axial Mach number problems that exist in other marching codes that utilize Cartesian transverse marching planes.

Siclari, M. J.↗
NASA NTRS1981

Magnetic and gravity anomalies in the Americas

The cleaning and magnetic tape storage of spherical Earth processing programs are reported. These programs include: NVERTSM which inverts total or vector magnetic anomaly data on a distribution of point dipoles in spherical coordinates; SMFLD which utilizes output from NVERTSM to compute total or vector magnetic anomaly fields for a distribution of point dipoles in spherical coordinates; NVERTG; and GFLD. Abstracts are presented for papers dealing with the mapping and modeling of magnetic and gravity anomalies, and with the verification of crustal components in satellite data.

Braile, L. W.↗
NASA NTRS1982

The reduction, verification and interpretation of MAGSAT magnetic data over Canada

Consideration is being given to representing the magnetic field in the area 40 deg N to 83 deg N by means of functions in spherical coordinates. A solution to Laplace's equation for the magnetic potential over a restricted area was found, and programming and testing are currently being carried out. Magnetic anomaly modelling is proceeding. The program SPHERE, which was adapted to function correctly on the Cyber computer, is now operational, for deriving gravity and magnetic models in a spherical coordinate system.

Coles, R. L.↗