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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
NASA NTRS1980

A three-dimensional body-fitted coordinate system for flow field calculations on asymmetric nosetips

A three dimensional body-fitted coordinate system developed for use in the calculation of inviscid flows over ablated, asymmetric reentry vehicle nosetips is described. Because of the potential geometric asymmetries, no standard coordinate system (e.g., spherical, axisymmetric reference surface-normal) is capable of being closely aligned with the nosetip surface. To generate a 3-D, body-fitted coordinate system an analytic mapping procedure is applied that is conformal within each meridional plane of the nosetip; these transformations are then coupled circumferentially to yield a three dimensional coordinate system. The mappings used are defined in terms of hinge points, which are points selected to approximate the body contours in each meridional plane. The selection of appropriate hinge points was automated to facilitate the use of the resulting nosetip flow field code.

Hall, D. W.↗
NASA NTRS1983

An interactive system for compositing digital radar and satellite data

This paper describes an approach for compositing digital radar data and GOES satellite data for meteorological analysis. The processing is performed on a user-oriented image processing system, and is designed to be used in the research mode. It has a capability to construct PPIs and three-dimensional CAPPIs using conventional as well as Doppler data, and to composite other types of data. In the remapping of radar data to satellite coordinates, two steps are necessary. First, PPI or CAPPI images are remapped onto a latitude-longitude projection. Then, the radar data are projected into satellite coordinates. The exact spherical trigonometric equations, and the approximations derived for simplifying the computations are given. The use of these approximations appears justified for most meteorological applications. The largest errors in the remapping procedure result from the satellite viewing angle parallax, which varies according to the cloud top height. The horizontal positional error due to this is of the order of the error in the assumed cloud height in mid-latitudes. Examples of PPI and CAPPI data composited with satellite data are given for Hurricane Frederic on 13 September 1979 and for a squall line on 2 May 1979 in Oklahoma.

Heymsfield, G. M.↗
NASA NTRS1992

Tests of a simple data merging algorithm for the GONG project

The GONG (Global Oscillation Network Group) project proposes to reduce the impact of diurnal variations on helioseismic measurements by making long-term observations of solar images from six sites placed around the globe. The sun will be observed nearly constantly for three years, resulting in the acquisition of l+ terabyte of image data. To use the solar network to maximum advantage, the images from the sites must be combined into a single time series to determine mode frequencies, amplitudes, and line widths. Initial versions of combined, i.e., merged, time series were made using a simple weighted average of data from different sites taken simultaneously. In order to accurately assess the impact of the data merge on the helioseismic measurements, a set of artificial solar disk images was made using a standard solar model and containing a well known set of oscillation modes and frequencies. This undegraded data set and data products computed from it were used to judge the relative merits of various data merging schemes. The artificial solar disk images were subjected to various instrumental and atmospheric degradations, dependent on site and time, in order to create a set of images simulating those likely to be taken at the site. The degraded artificial solar disk images for the six observing sites were combined in various ways to form merged time series of images and mode coefficients. Various forms of a weighted average were used, including an equally-weighted average, an average with weights dependent upon air mass and averages with weights dependent on various quality assurance parameters. Both the undegraded solar disk image time series and several time series made up of various combinations of the degraded solar disk images from the six sites were subjected to standard helioseismic measurement processing. This processing consisted of coordinate remapping, detrending, spherical harmonic transformation, computation of power series for the oscillation mode coefficients, and mode frequency identification. Visual and statistical evaluation of the merged data sets themselves and differences between the merged and undegraded data set shows good agreement between the two data sets. Some slight differences in image scale and registration appear between the undegraded data set and the various merged data sets. In the set of power series made from the mode coefficients of the merged data sets, some power leakage is observed into the background and into slightly lower l-value modes, especially at higher l-value mode frequencies. The results of the comparison of time series and mode oscillation frequencies of the undegraded data with those of the data merged using weighted averages indicate that, at least for p-mode solar oscillations, a weighted average of either the detrended remapped images or the mode coefficients gives good determination of the mode frequencies and adequate-to-good determination of the amplitudes and widths of the mode frequency lines. This conclusion is most advantageous to the analysis of the massive amounts of data to be received by the GONG network.

Williams, W. E.↗
NASA NTRS1996

An MHD Study of the Interaction Between the Solar Wind and the Interstellar Medium

The overall objective of this research program is to obtain a better understanding of the interaction between the solar wind and the interstellar medium through the use of numerical solutions of the time-dependent magnetohydrodynamic (MHD) equations. The simulated results have been compared with observations where possible and with the results from previous analytic and numerical studies. The primary accomplishment of this project has been the development of codes for 2-D models in both spherical and cylindrical coordinates and the application of the codes to the solar wind/interstellar medium interaction. Computations have been carried out for both a relatively simple gas-dynamic interaction and a flow-aligned interstellar magnetic field. The results have been shown to compare favorably with models that use more approximations and to modify and extend the previous results as would be expected. The simulations have also been used along with a data analysis study to provide a quantitative estimate of the distance to the termination and bow shocks. Some of the specific topics that have been studied are: (1) gas dynamic models of the solar wind/interstellar medium interaction, (2) termination shock response to large-scale solar wind fluctuations, and (3) distances to the termination shock and heliopause. The main results from each of these studies are summarized. The results were published in three papers which are included as attachments.

Steinolfson, R. S.↗
NASA NTRS1998

An Alternative Lunar Ephemeris Model for On-Board Flight Software Use

In calculating the position vector of the Moon in on-board flight software, one often begins by using a series expansion to calculate the ecliptic latitude and longitude of the Moon, referred to the mean ecliptic and equinox of date. One then performs a reduction for precession, followed by a rotation of the position vector from the ecliptic plane to the equator, and a transformation from spherical to Cartesian coordinates before finally arriving at the desired result: equatorial J2000 Cartesian components of the lunar position vector. An alternative method is developed here in which the equatorial J2000 Cartesian components of the lunar position vector are calculated directly by a series expansion, saving valuable on-board computer resources.

Simpson, David G.↗
NASA NTRS1999

An Alternative Lunar Ephemeris Model for On-Board Flight Software Use

In calculating the position vector of the Moon in on-board flight software, one often begins by using a series expansion to calculate the ecliptic latitude and longitude of the Moon, referred to the mean ecliptic and equinox of date. One then performs a reduction for precession, followed by a rotation of the position vector from the ecliptic plane to the equator, and a transformation from spherical to Cartesian coordinates before finally arriving at the desired result: equatorial J2000 Cartesian components of the lunar position vector. An alternative method is developed here in which the equatorial J2000 Cartesian components of the lunar position vector are calculated directly by a series expansion, saving valuable onboard computer resources.

Simpson, David G.↗
NASA NTRS1980

Orthogonality of spherical harmonic coefficients

Orthogonality relations are obtained for the spherical harmonic coefficients of functions defined on the surface of a sphere. Following a brief discussion of the orthogonality of Fourier series coefficients, consideration is given to the values averaged over all orientations of the coordinate system of the spherical harmonic coefficients of a function defined on the surface of a sphere that can be expressed in terms of Legendre polynomials for the special case where the function is the sum of two delta functions located at two different points on the sphere, and for the case of an essentially arbitrary function. It is noted that the orthogonality relations derived have found applications in statistical studies of the geomagnetic field.

Mcleod, M. G.↗
NASA NTRS1976

A tensor formulation of the equation of transfer for spherically symmetric flows

A tensor formulation of the equation of radiative transfer is derived in a seven-dimensional Riemannian space such that the resulting equation constitutes a divergence in any coordinate system. After being transformed to a spherically symmetric comoving coordinate system, the transfer equation contains partial derivatives in angle and frequency, as well as optical depth due to the effects of aberration and the Doppler shift. However, by virtue of the divergence form of this equation, the divergence theorem may be applied to yield a numerical differencing scheme which is expected to be stable and to conserve luminosity. It is shown that the equation of transfer derived by this method in a Lagrangian coordinate system may be reduced to that given by Castor (1972), although it is, of course, desirable to leave the equation in divergence form.

Haisch, B. M.↗
NASA NTRS1981

Spherical-earth Gravity and Magnetic Anomaly Modeling by Gauss-legendre Quadrature Integration

The usefulness of long-wavelength potential field anomalies in lithospheric interpretation is greatly increased with spherical Earth modeling techniques. Gauss-Legendre quadrature integration is used to calculate the anomalous potential of gravity and magnetic fields and their spatial derivatives on a spherical Earth for an arbitrary body represented by an equivalent point source distribution of gravity poles or magnetic dipoles. The distribution of equivalent point sources is determined directly from the coordinate limits of the source volume. Variable integration limits for an arbitrarily shaped body are derived from interpolation of points which approximate the body's surface envelope. The versatility of the method is enhanced by the ability to treat physical property variations within the source volume and to consider variable magnetic fields over the source and observation surface. A number of examples verify and illustrate the capabilities of the technique, including preliminary modeling of potential field signatures for Mississippi embayment crustal structure at satellite elevations.

Spherical coordinate gravity modeling↗
NASA NTRS1998

Mesogranulation as A Distinct Scale of Convection in the Sun

We present evidence for the existence of mesogranulation as a scale of convection distinct from granulation and supergranulation through analysis of full-disk Doppler velocity images of the Sun collected by the Michelson Doppler Imager (MDI) aboard the NASA/ESA Solar and Heliospheric Observatory (SOHO). Our analysis procedures isolate nearly steady flows in the solar photosphere and yield power spectra of convection for spherical harmonic degrees up to I = 1000. Each spectrum exhibits an obvious supergranulation peak at I approximately 130 and a broad secondary peak at I approximately 600 with a distinct break in the spectrum between these peaks at I approximately 300. We believe that this secondary peak is a signature of mesogranulation with typical cell diameters of about 7 Mm. Our standard analysis procedure is to first remove the p-mode oscillation signal by averaging individual Dopplergrams over 17-minute intervals. Next, by fitting to standard functional forms we remove Doppler signals due to the motion of the spacecraft, the convective blueshift, solar rotation including differential rotation, and the meridional circulation in order to produce Dopplergrams dominated by convective motions. By mapping these processed images onto heliographic coordinates and projecting onto spherical harmonics, we produce a power spectrum of solar convection for each 17-minute period. We construct synthetic images and pass them through the same analysis procedure in order to determine the actual solar convection spectrum that reproduces the analyzed results. We find that a small but increasing percentage of high-degree convective power is lost in the analysis as we approach the limit of resolution of the detector but'that the broad, mesogranulation peak at I approximately 600 must be included in the convection spectrum of the synthetic images.

Bachmann, Kurt T.↗
DOE OSTI2026

Entity—Hardware-agnostic Particle-in-cell Code for Plasma Astrophysics. I. Curvilinear Special Relativistic Module

Entity is a new-generation, fully open-source particle-in-cell (PIC) code developed to overcome key limitations in astrophysical plasma modeling, particularly the extreme separation of scales and the performance challenges associated with evolving, GPU-centric computing infrastructures. It achieves hardware-agnostic performance portability across various GPU and CPU architectures using the Kokkos library. Crucially, Entity maintains a high standard for usability, clarity, and customizability, offering a robust and easy-to-use framework for developing new algorithms and grid geometries, which allows extensive control without requiring edits to the core source code. This paper details the core general-coordinate special relativistic module. Entity is the first PIC code designed to solve the Vlasov–Maxwell system in general coordinates, enabling a coordinate-agnostic framework that provides the foundational structure for straightforward extension to arbitrary coordinate geometries. The core methodology achieves numerical stability by solving particle equations of motion in the global orthonormal Cartesian basis, despite using generalized coordinates like Cartesian, axisymmetric spherical, and quasi-spherical grids. Charge conservation is ensured via a specialized current deposition technique using conformal currents. The code exhibits robust scalability and performance portability on major GPU platforms (AMD MI250X, NVIDIA A100, and Intel Max Series), with the 3D particle pusher and the current deposition operating efficiently at about 2 ns per particle per time step. Functionality is validated through a comprehensive suite of standard Cartesian plasma tests and the accurate modeling of relativistic magnetospheres in curvilinear axisymmetric geometries.

Hakobyan, Hayk [Flatiron Institute, New York, NY (↗
NASA NTRS1972

Conjugate gradient optimization programs for shuttle reentry

Two computer programs for shuttle reentry trajectory optimization are listed and described. Both programs use the conjugate gradient method as the optimization procedure. The Phase 1 Program is developed in cartesian coordinates for a rotating spherical earth, and crossrange, downrange, maximum deceleration, total heating, and terminal speed, altitude, and flight path angle are included in the performance index. The programs make extensive use of subroutines so that they may be easily adapted to other atmospheric trajectory optimization problems.

Powers, W. F.↗
NASA NTRS1987

Tail modeling in a stretched magnetosphere. I - Methods and transformations

A new method is developed for representing the magnetospheric field B as a distorted dipole field. Because Delta-B = 0 must be maintained, such a distortion may be viewed as a transformation of the vector potential A. The simplest form is a one-dimensional 'stretch transformation' along the x axis, concisely represented by the 'stretch function' f(x), which is also a convenient tool for representing features of the substorm cycle. One-dimensional stretch transformations are extended to spherical, cylindrical, and parabolic coordinates and then to arbitrary coordinates. It is shown that distortion transformations can be viewed as mappings of field lines from one pattern to another; the final result only requires knowledge of the field and not of the potentials. General transformations in Cartesian and arbitrary coordinates are derived, and applications to field modeling, field line motion, MHD modeling, and incompressible fluid dynamics are considered.

Stern, David P.↗
NASA NTRS2007

Generalized Squashing Factors for Covariant Description of Magnetic Connectivity in the Solar Corona

The study of magnetic connectivity in the solar corona reveals a need to generalize the field line mapping technique to arbitrary geometry of the boundaries and systems of coordinates. Indeed, the global description of the connectivity in the corona requires the use of the photospheric and solar wind boundaries. Both are closed surfaces and therefore do not admit a global regular system of coordinates. At least two overlapping regular systems of coordinates for each of the boundaries are necessary in this case to avoid spherical-pole-like singularities in the coordinates of the footpoints. This implies that the basic characteristic of magnetic connectivity-the squashing degree or factor Q of elemental flux tubes, according to Titov and coworkers-must be rewritten in covariant form. Such a covariant expression of Q is derived in this work. The derived expression is very flexible and highly efficient for describing the global magnetic connectivity in the solar corona. In addition, a general expression for a new characteristic Q1, which defines a squashing of the flux tubes in the directions perpendicular to the field lines, is determined. This new quantity makes it possible to filter out the quasi-separatrix layers whose large values of Q are caused by a projection effect at the field lines nearly touching the photosphere. Thus, the value Q1 provides a much more precise description of the volumetric properties of the magnetic field structure. The difference between Q and Q1 is illustrated by comparing their distributions for two configurations, one of which is the Titov-Demoulin model of a twisted magnetic field.

Titov, V. S.↗
NASA NTRS2018

Symmetric Equations on the Surface of a Sphere as Used by Model GISS:IB

Standard vector calculus formulas of Cartesian three space are projected onto the surface of a sphere. This produces symmetric equations with three nonindependent horizontal velocity components. Each orthogonal axis has a velocity component that rotates around its axis (eastward velocity rotates around the north–south axis) and a specific angular momentum component that is the product of the velocity component multiplied by the cosine of axis’ latitude. Angular momentum components align with the fixed axes and simplify several formulas, whereas the rotating velocity components are not orthogonal and vary with location. Three symmetric coordinates allow vector resolution and calculus operations continuously over the whole spherical surface, which is not possible with only two coordinates. The symmetric equations are applied to one-layer shallow water models on cubed-sphere and icosahedral grids, the latter being computationally simple and applicable to an ocean domain. Model results are presented for three different initial conditions and five different resolutions.

Icosahedral grid↗