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

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↗

Latitude dependence of eddy variances

The eddy variance of a meteorological field must tend to zero at high latitudes due solely to the nature of spherical polar coordinates. The zonal averaging operator defines a length scale: the circumference of the latitude circle. When the circumference of the latitude circle is greater than the correlation length of the field, the eddy variance from transient eddies is the result of differences between statistically independent regions. When the circumference is less than the correlation length, the eddy variance is computed from points that are well correlated with each other, and so is reduced. The expansion of a field into zonal Fourier components is also influenced by the use of spherical coordinates. As is well known, a phenomenon of fixed wavelength will have different zonal wavenumbers at different latitudes. Simple analytical examples of these effects are presented along with an observational example from satellite ozone data. It is found that geometrical effects can be important even in middle latitudes.

Bowman, Kenneth P.↗

Protostellar hydrodynamics: Constructing and testing a spacially and temporally second-order accurate method. 2: Cartesian coordinates

In Boss & Myhill (1992) we described the derivation and testing of a spherical coordinate-based scheme for solving the hydrodynamic equations governing the gravitational collapse of nonisothermal, nonmagnetic, inviscid, radiative, three-dimensional protostellar clouds. Here we discuss a Cartesian coordinate-based scheme based on the same set of hydrodynamic equations. As with the spherical coorrdinate-based code, the Cartesian coordinate-based scheme employs explicit Eulerian methods which are both spatially and temporally second-order accurate. We begin by describing the hydrodynamic equations in Cartesian coordinates and the numerical methods used in this particular code. Following Finn & Hawley (1989), we pay special attention to the proper implementations of high-order accuracy, finite difference methods. We evaluate the ability of the Cartesian scheme to handle shock propagation problems, and through convergence testing, we show that the code is indeed second-order accurate. To compare the Cartesian scheme discussed here with the spherical coordinate-based scheme discussed in Boss & Myhill (1992), the two codes are used to calculate the standard isothermal collapse test case described by Bodenheimer & Boss (1981). We find that with the improved codes, the intermediate bar-configuration found previously disappears, and the cloud fragments directly into a binary protostellar system. Finally, we present the results from both codes of a new test for nonisothermal protostellar collapse.

Myhill, Elizabeth A.↗

Direct solution for thermal stresses in a nose cap under an arbitrary axisymmetric temperature distribution

The design of a nose cap for a hypersonic vehicle is an iterative process requiring a rapid, easy to use and accurate stress analysis. The objective of this paper is to develop such a stress analysis technique from a direct solution of the thermal stress equations for a spherical shell. The nose cap structure is treated as a thin spherical shell with an axisymmetric temperature distribution. The governing differential equations are solved by expressing the stress solution to the thermoelastic equations in terms of a series of derivatives of the Legendre polynomials. The process of finding the coefficients for the series solution in terms of the temperature distribution is generalized by expressing the temperature along the shell and through the thickness as a polynomial in the spherical angle coordinate. Under this generalization the orthogonality property of the Legendre polynomials leads to a sequence of integrals involving powers of the spherical shell coordinate times the derivative of the Legendre polynomials. The coefficients of the temperature polynomial appear outside of these integrals. Thus, the integrals are evaluated only once and their values tabulated for use with any arbitrary polynomial temperature distribution.

Davis, Randall C.↗

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 will be compared with observations where possible and with the results from previous analytic and numerical studies. The primary progress during the first two years has been to develop codes for 2-D models in both spherical and cylindrical coordinates and to apply them 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. Work has also been initiated on the development of a 3-D MHD code in spherical coordinates.

Steinolfson, R. S.↗

Gasdynamic models of the solar wind/interstellar medium interaction

The interaction between the solar wind and the interstellar medium is modeled self-consistently using numerical solutions of the time-dependent gasdynamic equations in spherical and cylindrical coordinates. For the results presented here it is assumed that the solar system moves through the surrounding medium with a supersonic velocity. After an initial (nonequilibrium) state has been specified, the numerical solution follows the evolution in time until the interaction relaxes to a dynamic equilibrium. As would be expected, the solutions show the formation of a bow shock upstream of the traveling solar system to deflect the interstellar plasma around the cavity created by the solar wind. A terminiation shock also forms to slow and compress the solar wind plasma. For the simulation in spherical coordinates, the downstream portion of the termination shock reaches equilibrium more than three times further from the Sun than the equilibrium distance to the termination shock on the upstream side.

Steinolfson, R. S.↗

A Lightning Channel Retrieval Algorithm for the North Alabama Lightning Mapping Array (LMA)

A new multi-station VHF time-of-arrival (TOA) antenna network is, at the time of this writing, coming on-line in Northern Alabama. The network, called the Lightning Mapping Array (LMA), employs GPS timing and detects VHF radiation from discrete segments (effectively point emitters) that comprise the channel of lightning strokes within cloud and ground flashes. The network will support on-going ground validation activities of the low Earth orbiting Lightning Imaging Sensor (LIS) satellite developed at NASA Marshall Space Flight Center (MSFC) in Huntsville, Alabama. It will also provide for many interesting and detailed studies of the distribution and evolution of thunderstorms and lightning in the Tennessee Valley, and will offer many interesting comparisons with other meteorological/geophysical wets associated with lightning and thunderstorms. In order to take full advantage of these benefits, it is essential that the LMA channel mapping accuracy (in both space and time) be fully characterized and optimized. In this study, a new revised channel mapping retrieval algorithm is introduced. The algorithm is an extension of earlier work provided in Koshak and Solakiewicz (1996) in the analysis of the NASA Kennedy Space Center (KSC) Lightning Detection and Ranging (LDAR) system. As in the 1996 study, direct algebraic solutions are obtained by inverting a simple linear system of equations, thereby making computer searches through a multi-dimensional parameter domain of a Chi-Squared function unnecessary. However, the new algorithm is developed completely in spherical Earth-centered coordinates (longitude, latitude, altitude), rather than in the (x, y, z) cartesian coordinates employed in the 1996 study. Hence, no mathematical transformations from (x, y, z) into spherical coordinates are required (such transformations involve more numerical error propagation, more computer program coding, and slightly more CPU computing time). The new algorithm also has a more realistic definition of source altitude that accounts for Earth oblateness (this can become important for sources that are hundreds of kilometers away from the network). In addition, the new algorithm is being applied to analyze computer simulated LMA datasets in order to obtain detailed location/time retrieval error maps for sources in and around the LMA network. These maps will provide a more comprehensive analysis of retrieval errors for LMA than the 1996 study did of LDAR retrieval errors. Finally, we note that the new algorithm can be applied to LDAR, and essentially any other multi-station TWA network that depends on direct line-of-site antenna excitation.

Koshak, William↗