Search NASASearch

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 91 records · Page 5

Numerical simulation of extended corona

A three-dimensional, time-dependent MHD model is presented for the study of coronal dynamics. The model, written in spherical coordinates, extends from the solar surface and was developed with two major issues in mind, namely for interpretation of various steady state and evolutionary dynamical structures in the corona. In order to achieve these objectives, two different numerical techniques are used to seek solutions for these two different, but related, problems; steady state structures and evolutionary structures. These two numerical techniques are: (1) relaxation technique for steady state structures; and (2) FICE (Full-Implicit-Continuous-Eulerian) technique for evolutionary structures. To illustrate this model, numerical results are given for examples of both the steady state and evolutionary structure of the corona which show the additional physical features which cannot be shown by a two-dimensional model.

Wu, S. T.

A new method for visualizing data on a sphere

A method for visualizing data on a globe or unit sphere is described. Information that is distributed over a sphere - global oceanographic or geographic measurements, all-sky astronomy observations, or any quantities that are best represented in spherical coordinates - can benefit from this technique. Retaining a better sense of the geometry and information content of the data, 3D graphics can provide an unobstructed view of the entire sphere, without undue deformation of its surface area. A 'parameterized ray trace' produces look-up tables (LUTs) that can be used for all visualizations. The ray-trace result shows one or more spheres with the data as a texture map and three reflecting rectangles that 'mirror' the far sides of the sphere(s) into view. The LUTs need only be created once. No special purpose hardware is required beyond a PC or workstation that supports color. Examples from astronomical and geophysical datasets, which are commonly displayed with an area deforming (2D) projection, are presented.

Hon, David

Application of a monotonic upstream-biased transport scheme to three-dimensional constituent transport calculations

The application of van Leer's scheme, a monotonic, upstream-biased differencing scheme, to three-dimensional constituent transport calculations is shown. The major disadvantage of the scheme is shown to be a self-limiting diffusion. A major advantage of the scheme is shown to be its ability to maintain constituent correlations. The scheme is adapted for a spherical coordinate system with a hybrid sigma-pressure coordinate in the vertical. Special consideration is given to cross-polar flow. The vertical wind calculation is shown to be extremely sensitive to the method of calculating the divergence. This sensitivity implies that a vertical wind formulation consistent with the transport scheme is essential for accurate transport calculations. The computational savings of the time-splitting method used to solve this equation are shown. Finally, the capabilities of this scheme are illustrated by an ozone transport and chemistry model simulation.

Allen, Dale J.

Wigner functions for nonclassical states of a collection of two-level atoms

The general theory of atomic angular momentum states is used to derive the Wigner distribution function for atomic angular momentum number states, coherent states, and squeezed states. These Wigner functions W(theta,phi) are represented as a pseudo-probability distribution in spherical coordinates theta and phi on the surface of a sphere of radius the square root of j(j +1) where j is the total angular momentum.

Agarwal, G. S.

Response-coefficient method for heat-conduction transients with time-dependent inputs

A theoretical overview of the response coefficient method for heat conduction transients with time-dependent input forcing functions is presented with a number of illustrative applications. The method may be the most convenient and economical if the same problem is to be solved many times with different input-time histories or if the solution time is relatively long. The method is applicable to a wide variety of problems, including irregular geometries, position-dependent boundary conditions, position-dependent physical properties, and nonperiodic irregular input histories. Nonuniform internal energy generation rates within the structure can also be handled by the method. The area of interest is long-time solutions, in which initial condition is unimportant, and not the early transient period. The method can be applied to one dimensional problems in cartesian, cylindrical, and spherical coordinates as well as to two dimensional problems in cartesian and cylindrical coordinates.

Ceylan, Tamer

A global 3-D MHD model of the solar wind with Alfven waves

A fully three-dimensional solar wind model that incorporates momentum and heat addition from Alfven waves is developed. The proposed model upgrades the previous one by considering self-consistently the total system consisting of Alfven waves propagating outward from the Sun and the mean polytropic solar wind flow. The simulation region extends from the coronal base (1 R(sub s) out to beyond 1 AU. The fully 3-D MHD equations written in spherical coordinates are solved in the frame of reference corotating with the Sun. At the inner boundary, the photospheric magnetic field observations are taken as boundary condition and wave energy influx is prescribed to be proportional to the magnetic field strength. The results of the model application for several time intervals are presented.

Usmanov, A. V.

The Development of Kolmogoroff-Like Power Spectra in the Expanding Solar Wind

Power spectra of the solar wind fluctuations consistently exhibit a -5/3 power-law slope consistent with the idea that the medium is undergoing a turbulent cascade as seen in ordinary fluids. This is surprising both because the radial streams and the magnetic field threading the plasma will induce anisotropies and because the expansion of the wind will tend to lead to the suppression of nonlinear cascades. These conditions violate the assumptions used by Kolmogoroff to derive the -5/3 law. We have studied this issue using a compressible IMHD code in spherical coordinates and have shown that a -5/3 spectrum results from a broad-band flat-spectrum input condition that is sheared and distorted by a current sheet as the wind expands. We determine spectra from time series taken at selected points in the domain as is done with observational spacecraft data. The spectra are very like those we have seen in nonexpanding runs and exhibit evolution and compressive characteristics very similar to those seen in observations. We will report on these results in addition to a new set of runs intended to constrain the necessary and sufficient conditions for the spectra to have this form. The simulation also allows us to examine the anisotropy for the spectra to attempt to determine why the result of an isotropic magnetofluid is obtained in a highly anisotropic situation.

Goldstein, M. L.

Short-Arc Analysis of Intersatellite Tracking Data in a Gravity Mapping Mission

A technique for the analysis of low-low intersatellite range-rate data in a gravity mapping mission is explored. The technique is based on standard tracking data analysis for orbit determination but uses a spherical coordinate representation of the 12 epoch state parameters describing the baseline between the two satellites. This representation of the state parameters is exploited to allow the intersatellite range-rate analysis to benefit from information provided by other tracking data types without large simultaneous multiple data type solutions. The technique appears especially valuable for estimating gravity from short arcs (e.g., less than 15 minutes) of data. Gravity recovery simulations which use short arcs are compared with those using arcs a day in length. For a high-inclination orbit, the short-arc analysis recovers low-order gravity coefficients remarkably well, although higher order terms, especially sectorial terms, are less accurate. Simulations suggest that either long or short arcs of GRACE data are likely to improve parts of the geopotential spectrum by orders of magnitude.

Rowlands, David D.

Large Eddy Simulation of a Turbulent Jet

Here we present the results of a Large Eddy Simulation of a non-buoyant jet issuing from a circular orifice in a wall, and developing in neutral surroundings. The effects of the subgrid scales on the large eddies have been modeled with the dynamic large eddy simulation model applied to the fully 3D domain in spherical coordinates. The simulation captures the unsteady motions of the large-scales within the jet as well as the laminar motions in the entrainment region surrounding the jet. The computed time-averaged statistics (mean velocity, concentration, and turbulence parameters) compare well with laboratory data without invoking an empirical entrainment coefficient as employed by line integral models. The use of the large eddy simulation technique allows examination of unsteady and inhomogeneous features such as the evolution of eddies and the details of the entrainment process.

Webb, A. T.

MHD Modeling of Differential Rotation in Coronal Holes

The photosphere and the magnetic flux therein undergo differential rotation. Coronal holes appear to rotate almost rigidly. Magnetic reconnection has been invoked to reconcile these phenomena. Mechanism relevant to the formation of the slow solar wind. We have used our MHD model in spherical coordinates to study the effect of differential rotation on coronal holes. We have imposed a magnetic flux distribution similar to and applied differential rotation for the equivalent of 5 solar rotations.

Lionello, Roberto

Solar Simulation

Simulation of the surface region and interior of the Sun's convection zone. The objective is to provide understanding of the tachocline at the base of the convection zone. of the differential rotation generated by global convective motions, and of supergranules that are observed on the solar surface. Two large-scale simulation codes are used, one for the global dynamics in spherical coordinates and another for the local surface events in Cartesian geometry. We work in close collaboration with the Helioseismic and Magnetic Imager instrument developers for the Solar Dynamics Observatory satellite. NASA Program: Living with a Star - Sun-Earth Connection.

Wray, Alan

MHD Simulation Heliospheric Magnetic Fields and Turbulence

This talk will present a summary of our results on simulations of heliospheric structure and dynamics. We use a three-dimensional MHD code in spherical coordinates to produce a solar wind containing a rotating, tilted heliospheric current sheet, fast-slow stream and microstream shear layers, waves, 2-D turbulence, and pressure balanced structures that are input to the inner (superAlfvenic) boundary. The evolution of various combinations of these has led to a deeper understanding of sector structure, magnetic holes, fluctuation anisotropies, and general turbulent evolution. We show how the sectors are likely to be connected, how spiral fields can arise, and how field line diffusion can be caused by waves with transverse structure and microstream shears.

Roberts, D. Aaron

Efforts to Simulate Solar Wind Turbulence

A three-dimensional integration of the MHD equations in spherical coordinates has been developed that attempts to simulate a variety of solar wind conditions. These include the interaction of Alfven wave packets and the development of a turbulent cascade, the role of the heliospheric current sheet, the role of quasi-two-dimensional fluctuations in determining how magnetic field lines meander throughout the heliosphere, and the role of interstellar pickup ions in perturbing the solar wind in the outer heliosphere.

Goldstein, Melvyn L.

Displaying and Analyzing Antenna Radiation Patterns

Radiant Energy Display and Analysis Software Package (REDAP) is a computer program for processing antenna-radiation- pattern data that have been preprocessed by a data-collection program and stored in a spherical-coordinate format. REDAP is designed specifically for application to data generated in testing of antennas in an anechoic chamber at Johnson Space Center; parts of REDAP may be reusable for processing antenna-test data collected elsewhere. REDAP provides a graphical user interface (GUI) and executes mathematical and plotting routines on a personal computer.

Siekierski, Jim

Interchange Reconnection and Coronal Hole Dynamics

We investigate the effect of magnetic reconnection between open and closed fields, often referred to as interchange reconnection, on the dynamics and topology of coronal hole boundaries. The most important and most prevalent three-dimensional topology of the interchange process is that of a small-scale bipolar magnetic field interacting with a large-scale background field. We determine the evolution of such a magnetic topology by numerical solution of the fully three-dimensional MHD equations in spherical coordinates. First, we calculate the evolution of a small-scale bipole that initially is completely inside an open field region and then is driven across a coronal hole boundary by photospheric motions. Next the reverse situation is calculated in which the bipole is initially inside the closed region and driven toward the coronal hole boundary. In both cases, we find that the stress imparted by the photospheric motions results in deformation of the separatrix surface between the closed field of the bipole and the background field, leading to rapid current sheet formation and to efficient reconnection. When the bipole is inside the open field region, the reconnection is of the interchange type in that it exchanges open and closed fields. We examine, in detail, the topology of the field as the bipole moves across the coronal hole boundary and find that the field remains well connected throughout this process. Our results, therefore, provide essential support for the quasi-steady models of the open field, because in these models the open and closed flux are assumed to remain topologically distinct as the photosphere evolves. Our results also support the uniqueness hypothesis for open field regions as postulated by Antiochos et al. On the other hand, the results argue against models in which open flux is assumed to diffusively penetrate deeply inside the closed field region under a helmet streamer. We discuss the implications of this work for coronal observations.

MHD (MAGNETOHYDRODYNAMICS)

Interchange Reconnection and Coronal Hole Dynamics

We investigate the effect of magnetic reconnection between open and closed field, (often referred to as "interchange" reconnection), on the dynamics and topology of coronal hole boundaries. The most important and most prevalent 3D topology of the interchange process is that of a small-scale bipolar magnetic field interacting with a large-scale background field. We determine the evolution of such a magnetic topology by numerical solution of the fully 3D MHD equations in spherical coordinates. First, we calculate the evolution of a small-scale bipole that initially is completely inside an open field region and then is driven across a coronal hole boundary by photospheric motions. Next the reverse situation is calculated in which the bipole is initially inside the closed region and driven toward the coronal hole boundary. In both cases we find that the stress imparted by the photospheric motions results in deformation of the separatrix surface between the closed field of the bipole and the background field, leading to rapid current sheet formation and to efficient reconnection. When the bipole is inside the open field region, the reconnection is of the interchange type in that it exchanges open and closed field. We examine, in detail, the topology of the field as the bipole moves across the coronal hole boundary, and find that the field remains well-connected throughout this process. Our results imply that open flux cannot penetrate deeply into the closed field region below a helmet streamer and, hence, support the quasi-steady models in which open and closed flux remain topologically distinct. Our results also support the uniqueness hypothesis for open field regions as postulated by Antiochos et al. We discuss the implications of this work for coronal observations. Subject Headings: Sun: corona Sun: magnetic fields Sun: reconnection Sun: coronal hole

Edmondson, J. K.

LIDAR Metrology for Prescription Characterization and Alignment of Large Mirrors

We describe the use of LIDAR, or "laser radar," (LR) as a fast, accurate, and non-contact tool for the measurement of the radius of curvature (RoC) of large mirrors. We report the results of a demonstration of this concept using a commercial laser radar system. We measured the RoC of a 1.4m x 1m spherical mirror with a nominal RoC of 4.6 m with a manufacturing tolerance of 4600mm +/- 6mm. The prescription of the mirror is related to its role as ground support equipment used in the test of part of the James Webb Space Telescope (JWST). The RoC of such a large mirror is not easily measured without contacting the surface. From a position near the center of curvature of the mirror, the LIDAR scanned the mirror surface, sampling it with 1 point per 3.5 sq cm. The measurement consisted of 3983 points and lasted only a few minutes. The laser radar uses the LIDAR signal to provide range, and encoder information from angular azimuth and elevation rotation stages provide the spherical coordinates of a given point. A best-fit to a sphere of the measured points was performed. The resulting RoC was within 20 ppm of the nominal RoC, also showing good agreement with the results of a laser tracker-based, contact metrology. This paper also discusses parameters such as test alignment, scan density, and optical surface type, as well as future possible application for full prescription characterization of aspherical mirrors, including radius, conic, off-axis distance, and aperture.

Eegholm, B.

A NICER view of PSR J0030+0541: Millisecond Pulsar Parameter Estimation

We report on Bayesian parameter estimation of the mass and equatorial radius of the millisecond pulsar PSRJ0030+0451, conditional on pulse-profile modeling of Neutron Star Interior Composition Explorer X-ray spectral timing event data. We perform relativistic ray-tracing of thermal emission from hot regions of the pulsar’s surface. We assume two distinct hot regions based on two clear pulsed components in the phase-folded pulse-profile data; we explore a number of forms (morphologies and topologies) for each hot region, inferring their parameters in addition to the stellar mass and radius. For the family of models considered, the evidence (prior predictive probability of the data) strongly favors a model that permits both hot regions to be located in the same rotational hemisphere. Models wherein both hot regions are assumed to be simply connected circular single-temperature spots, in particular those where the spots are assumed to be reflection-symmetric with respect to the stellar origin, are strongly disfavored. For the inferred configuration, one hot region subtends an angular extent of only a few degrees (in spherical coordinates with origin at the stellar center) and we are insensitive to other structural details; the second hot region is far more azimuthally extended in the form of a narrow arc, thus requiring a larger number of parameters to describe. The inferred mass M and equatorial radius Req are, respectively, -1.34+ M 0.16 0.15 and-12.71+ km 1.19 1.14 , while the compactness = -GM R c 0.156+ eq 2 0.010 0.008 is more tightly constrained; the credible interval bounds reported here are approximately the 16% and 84% quantiles in marginal posterior mass.

T E Riley