Search NASA⌕ Search

SEARCH · Search NASA

Results for “spherical geometry”

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 145 records · Page 8

Disruption of coronal magnetic field arcades

The ideal and resistive properties of isolated large-scale coronal magnetic arcades are studied using axisymmetric solutions of the time-dependent magnetohydrodynamic (MHD) equations in spherical geometry. We examine how flares and coronal mass ejections may be initiated by sudden disruptions of the magnetic field. The evolution of coronal arcades in response to applied shearing photospheric flows indicates that disruptive behavior can occur beyond a critical shear. The disruption can be traced to ideal MHD magnetic nonequilibrium. The magnetic field expands outward in a process that opens the field lines and produces a tangential discontinuity in the magnetic field. In the presence of plasma resistivity, the resulting current sheet is the site of rapid reconnection, leading to an impulsive release of magnetic energy, fast flows, and the ejection of a plasmoid. We relate these results to previous studies of force-free fields and to the properties of the open-field configuration. We show that the field lines in an arcade are forced open when the magnetic energy approaches (but is still below) the open-field energy, creating a partially open field in which most of the field lines extend away from the solar surface. Preliminary application of this model to helmet streamers indicates that it is relevant to the initiation of coronal mass ejections.

Mikic, Zoran↗

A Very Large Array search for emission from HI associated with nearby Lyman alpha absorbers

We present a sensitive Very Large Array (VLA) search for H I emission from the vicinity of the Lyman alpha clouds in the Virgo Cluster, which were recently discovered with Hubble Space Telescope (HST) in absorption toward the nearest quasar 3C273. We searched an area of 40 by 40 arcmin centered on 3C273, covering a velocity range from 840 to 1840 km/s. The bandpass was self-calibrated on 3C273 leading to a spectral dynamic range of better than 10(exp 5) to 1. No H I was detected. The rms noise in the final images corresponds to a 3 sigma column density sensitivity of 2.8 x 10(exp 19) sq cm on scales of a few kpc. Small H I clouds could have been detected down to a few times 10(exp 6) solar mass. Our failure to detect H I emission at the higher column densities sets a lower limit to the radius of the Lyman alpha clouds of 3.9 kpc, assuming a spherical geometry.

Van Gorkom, J. H.↗

A comparison of eastern North American seismic strain-rates to glacial rebound strain-rates

Glacial rebound strain-rates computed using a simple Laurentide glacial loading model are of the order of 10(exp -9) per year within the region of glaciation and extending several hundred kilometers beyond. The horizontal strain-rates receive approximately equal contributions from horizontal and vertical velocities, a consequence of the spherical geometry adopted for the Earth model. In the eastern United States and southeastern Canada the computed strain-rates are 1-3 orders of magnitude greater than an estimate of the average seismic strain-rate (Anderson, 1986) and approximately 1 order of magnitude greater than predicted erosional strain-rates. The predicted glacial rebound strain-rates are not, in general, oriented in such a way as to augment the observed state of deviatoric stress, possibly explaining why the seismic strain-rates are much smaller than the glacial rebound strain-rates. An exception to this may be seismically active regions in the St. Lawrence valley.

James, Thomas S.↗

Photochemical modeling of Titan's atmosphere

We have developed a new photochemical model of Titan's atmosphere which includes all the important compounds and reactions in spherical geometry from the surface to 1240 km. Compared to the previous model of Yung et al. (1984), the most significant recent change in the reactions used is the updated methane dissociation scheme (Mordaunt et al. 1993). Moreover, the transfer of the solar radiation in the atmosphere and the photolysis rates have been calculated by using a Monte Carlo code. Finally, the eddy diffusion coefficient profile is adjusted in order to fit the mean vertical distribution of HCN retrieved from millimeter groundbased observations of Tanguy et al. (1990); using new values for the boundary flux of atomic nitrogen (Strobel et al. 1992). We have run the model in both steady-state and diurnal modes, with 62 species involved in 249 reactions. There is little difference between diurnal and steady-state results. Overall our results are in a closer agreement with the abundances inferred from the Voyager infrared measurements at the equator than the Yung et al. results. We find that the catalytic scheme for H recombination invoked by Yung et al. only slightly improves the model results and we conclude that this scheme is not essential to fit observations.

Toublanc, D.↗

A Large-Particle Monte Carlo Code for Simulating Non-Linear High-Energy Processes Near Compact Objects

High-energy radiation processes in compact cosmic objects are often expected to have a strongly non-linear behavior. Such behavior is shown, for example, by electron-positron pair cascades and the time evolution of relativistic proton distributions in dense radiation fields. Three independent techniques have been developed to simulate these non-linear problems: the kinetic equation approach; the phase-space density (PSD) Monte Carlo method; and the large-particle (LP) Monte Carlo method. In this paper, we present the latest version of the LP method and compare it with the other methods. The efficiency of the method in treating geometrically complex problems is illustrated by showing results of simulations of 1D, 2D and 3D systems. The method is shown to be powerful enough to treat non-spherical geometries, including such effects as bulk motion of the background plasma, reflection of radiation from cold matter, and anisotropic distributions of radiating particles. It can therefore be applied to simulate high-energy processes in such astrophysical systems as accretion discs with coronae, relativistic jets, pulsar magnetospheres and gamma-ray bursts.

Stern, Boris E.↗

Experiments on Thermal Convection in Rotating Spherical Shells With Radial Gravity: The Geophysical Fluid Flow Cell

Experiments designed to study the fluid dynamics of buoyancy driven circulations in rotating spherical shells were conducted on the United States Microgravity Laboratory 2 spacelab mission. These experiments address several aspects of prototypical global convection relevant to large scale motions on the Sun, Earth, and on the giant planets. The key feature is the consistent modeling of radially directed gravity in spherical geometry by using dielectric polarization forces. Imagery of the planforms of thermally driven flows for rapidly-rotating regimes shows an initial separation and eventual merger of equatorial and polar convection as the heating (i.e. the Rayleigh number) is increased. At low rotation rates, multiple-states of motion for the same external parameters were observed.

Hart, John E.↗

Global Dynamic Modeling of Space-Geodetic Data

The proposal had outlined a year for program conversion, a year for testing and debugging, and two years for numerical experiments. We kept to that schedule. In first (partial) year, author designed a finite element for isostatic thin-shell deformation on a sphere, derived all of its algebraic and stiffness properties, and embedded it in a new finite element code which derives its basic solution strategy (and some critical subroutines) from earlier flat-Earth codes. Also designed and programmed a new fault element to represent faults along plate boundaries. Wrote a preliminary version of a spherical graphics program for the display of output. Tested this new code for accuracy on individual model plates. Made estimates of the computer-time/cost efficiency of the code for whole-earth grids, which were reasonable. Finally, converted an interactive graphical grid-designer program from Cartesian to spherical geometry to permit the beginning of serious modeling. For reasons of cost efficiency, models are isostatic, and do not consider the local effects of unsupported loads or bending stresses. The requirements are: (1) ability to represent rigid rotation on a sphere; (2) ability to represent a spatially uniform strain-rate tensor in the limit of small elements; and (3) continuity of velocity across all element boundaries. Author designed a 3-node triangle shell element which has two different sets of basis functions to represent (vector) velocity and all other (scalar) variables. Such elements can be shown to converge to the formulas for plane triangles in the limit of small size, but can also applied to cover any area smaller than a hemisphere. The difficult volume integrals involved in computing the stiffness of such elements are performed numerically using 7 Gauss integration points on the surface of the sphere, beneath each of which a vertical integral is performed using about 100 points.

Bird, Peter↗

Experimental Observations of PMMA Spheres Burning at Reduced Gravity

Polymer combustion is a very complicated process which involves the coupling of gas-phase combustion, melting, pyrolysis, and possible charring of the condensed phase. Although only a few studies on the combustion of a spherically shaped polymer have been conducted, there is renewed interest in the subject. At the Third International Microgravity Combustion Workshop, we presented a preliminary experimental plan and proposed an apparatus to study the combustion of a PolyMethylMethAcrylate (PMMA) sphere at reduced gravity. In this paper, we describe the experimental hardware in detail, summarize our observations since the last workshop, and describe future studies. The main objective is to determine the burning rate of PMMA spheres at reduced gravity under different ambient oxygen concentrations and total pressures. The dependence of the burning rate on the initial sphere diameter will be examined. It is anticipated that the simple spherical geometry in conjunction with reduced gravity will facilitate an assessment of the effect of condensed phase behavior on polymer burning processes and can be used as a means to test polymer combustion models.

Yang, Ji-Ann C.↗

Signatures of Nonlinear Waves in Coronal Plumes and Holes

In recent Ultraviolet Coronagraph Spectrometer/Solar and Heliospheric Observatory (UVCS/SOHO) White Light Channel (WLC) observations we found quasi-periodic variations in the polarized brightness (pB) in the polar coronal holes at heliocentric distances of 1.9-2.45 solar radii. The motivation for the observation is the 2.5D Magnetohydrodynamics (MHD) model of solar wind acceleration by nonlinear waves, that predicts compressive fluctuations in coronal holes. To help identify the waves observed with the UVCS/WLC we model the propagation and dissipation of slow magnetosonic waves in polar plumes using 1D MHD code in spherical geometry, We find that the slow waves nonlinearly steepen in the gravitationally stratified plumes. The nonlinear steepening of the waves leads to enhanced dissipation due to compressive viscosity at the wave-fronts.

Ofman, Leon↗

Theoretical and Observational Determination of Global and Regional Radiation Budget, Forcing and Feedbacks at the Top-of-Atmosphere and Surface

Report consists of: 1. List of accomplishments 2. List of publications 3. Abstracts of published or submitted papers and 4. Subject invention disclosure. The accomplishments of the grant listed are: 1. Improved the third-order turbulence closure in cloud resolving models to remove the liquid water oscillation. 2. Used the University of California-Los Angeles (UCLA) large-eddy simulation (LES) model to provide data for radiation transfer testing. 3. Revised shortwave k-distribution models based on HITRAN 2000. 4. Developed a gamma-weighted two-stream radiative transfer model for radiation budget estimate applications. 5. Estimated the effect of spherical geometry to the earth radiation budget. 6. Estimated top-of-atmosphere irradiance over snow and sea ice surfaces. 7. Estimated the aerosol direct radiative effect at the top of the atmosphere. 8. Estimated the top-of-atmosphere reflectance of the clear-sky molecular atmosphere over ocean. 9. Developed and validated new set of Angular Distribution Models for the CERES TRMM satellite instrument (tropical) 10. Developed and validated new set of Angular Distribution Models for the CERES Terra satellite instrument (global) 11. Quantified the top-of-atmosphere direct radiative effect of aerosols over global oceans from merged CERES and MODIS observations 12 Clarified the definition of TOA flux reference level for radiation budget studies 13. Developed new algorithm for unfaltering CERES measured radiances 14. Used multiangle POLDER measurements to produce narrowband angular distribution models and examine the effect of scene identification errors on TOA albedo estimates 15. Developed and validated a novel algorithm called the Multidirectional Reflectance Matching (MRM) model for inferring TOA albedos from ice clouds using multi-angle satellite measurements. 16. Developed and validated a novel algorithm called the Multidirectional Polarized Reflectance Matching (MPRM) model for inferring particle shapes from ice clouds using multi-angle polarized satellite measurements. 17. Developed 4 advanced light scattering models including the three-dimensional (3D) uniaxial perfectly matched layer (UPML) finite-difference time-domain (FDTD) model. 18. Develop sunglint in situ measurement and study reflectance distribution in the sunglint area. 19. Lead a balloon-borne radiometer TOA albedo validation effort. 20. Developed a CERES surface UVB, UVA, and UV index product.

Loeb, Norman G.↗

INTEGRAL Observations of the Be/X-ray binary EX0 2030+375 During Outburst

We present a type-I outburst of the high-mass X-ray binary EX0 2030+375, detected during INTEGRAL'S Performance and Verification phase in December 2002 (on-source time about 10(exp 6) seconds). In addition, six more outbursts have been observed during INTEGRAL'S Galactic Plane Scans. X-ray pulsations have been detected with a pulse period of 41.691798 plus or minus 0.000016 s. The X-ray luminosity in the 5-300 keV energy range was 9.7 x 10 (exp 36) erg per second, for a distance of 7.1 kpc. Two unusual features were found in the light curve, with an initial peak before the main outburst and another possible spike after the maximum. RXTE observations confirm only the existence of the initial spike. Although the initial peak appears to be a recurrent feature, the physical mechanisms producing it and the possible second spike are unknown. Moreover, a four-day delay between periastron passage and the peak of the outburst is observed. We present for the first time a 5-300 keV broad-band spectrum of this source. It can be modelled by the sum of a disk black body (kT(sub BB) approximately 8 keV) and either with a power law model with Gamma=2.04 plus or minus 0.11 keV or a Comptonized component (spherical geometry, kT(sub e).=30 keV, tau = 2.64, kT(sub w)=1.5 keV).

Arranz, A. Camero↗

Metal Hydride and Alkali Halide Opacities in Extrasolar Giant Planets and Cool Stellar Atmospheres

The lack of accurate and complete molecular line and continuum opacity data has been a serious limitation to developing atmospheric models of cool stars and Extrasolar Giant Planets (EGPs). We report our recent calculations of molecular opacities resulting from the presence of metal hydrides and alkali halides. The resulting data have been included in the PHOENIX stellar atmosphere code (Hauschildt & Baron 1999). The new models, calculated using spherical geometry for all gravities considered, also incorporate our latest database of nearly 670 million molecular lines, and updated equations of state.

Weck, Philippe F.↗

1997 May 12 Coronal Mass Ejection Event. I. A Simplified Model of the Preeruptive Magnetic Structure

A simple model of the coronal magnetic field prior to the coronal mass ejection (CME) eruption on 1997 May 12 is developed. First, the magnetic field is constructed by superimposing a large-scale background field and a localized bipolar field to model the active region (AR) in the current-free approximation. Second, this potential configuration is quasi-statically sheared by photospheric vortex motions applied to two flux concentrations of the AR. Third, the resulting force-free field is then evolved by canceling the photospheric magnetic flux with the help of an appropriate tangential electric field applied to the central part of the AR. To understand the structure of the modeled configuration, we use the field line mapping technique by generalizing it to spherical geometry. We demonstrate that the initial potential configuration contains a hyperbolic flux tube (HFT) which is a union of two intersecting quasi-separatrix layers. This HFT provides a partition of the closed magnetic flux between the AR and the global solar magnetic field. Such a partition is approximate since the entire flux distribution is perfectly continuous. The vortex motions applied to the AR inter-lock the field lines in the coronal volume to form additionally two new HFTs pinched into thin current layers. Reconnection in these current layers helps to redistribute the magnetic flux and current within the AR in the flux-cancellation phase. In this phase, a magnetic flux rope is formed together with a bald patch separatrix surface wrapping around the rope. Other important implications of the identified structural features of the modeled configuration are also discussed.

Titov, V. S.↗

Topological Evolution of a Fast Magnetic Breakout CME in 3-Dimensions

W present the extension of the magnetic breakout model for CME initiation to a fully 3-dimensional, spherical geometry. Given the increased complexity of the dynamic magnetic field interactions in 3-dimensions, we first present a summary of the well known axisymmetric breakout scenario in terms of the topological evolution associated with the various phases of the eruptive process. In this context, we discuss the completely analogous topological evolution during the magnetic breakout CME initiation process in the simplest 3-dimensional multipolar system. We show that an extended bipolar active region embedded in an oppositely directed background dipole field has all the necessary topological features required for magnetic breakout, i.e. a fan separatrix surface between the two distinct flux systems, a pair of spine fieldlines, and a true 3-dimensional coronal null point at their intersection. We then present the results of a numerical MHD simulation of this 3-dimensional system where boundary shearing flows introduce free magnetic energy, eventually leading to a fast magnetic breakout CME. The eruptive flare reconnection facilitates the rapid conversion of this stored free magnetic energy into kinetic energy and the associated acceleration causes the erupting field and plasma structure to reach an asymptotic eruption velocity of greater than or approx. equal to 1100 km/s over an approx.15 minute time period. The simulation results are discussed using the topological insight developed to interpret the various phases of the eruption and the complex, dynamic, and interacting magnetic field structures.

Lynch, B. J.↗

Venusian Applications of 3D Convection Modeling

This study models mantle convection on Venus using the 'cubed sphere' code OEDIPUS, which models one-sixth of the planet in spherical geometry. We are attempting to balance internal heating, bottom mantle viscosity, and temperature difference across Venus' mantle, in order to create a realistic model that matches with current planetary observations. We also have begun to run both lower and upper mantle simulations to determine whether layered (as opposed to whole-mantle) convection might produce more efficient heat transfer, as well as to model coronae formation in the upper mantle. Upper mantle simulations are completed using OEDIPUS' Cartesian counterpart, JOCASTA. This summer's central question has been how to define a mantle plume. Traditionally, we have defined a hot plume the region with temperature at or above 40% of the difference between the maximum and horizontally averaged temperature, and a cold plume as the region with 40% of the difference between the minimum and average temperature. For less viscous cases (1020 Pa?s), the plumes generated by that definition lacked vigor, displaying buoyancies 1/100th of those found in previous, higher viscosity simulations (1021 Pa?s). As the mantle plumes with large buoyancy flux are most likely to produce topographic uplift and volcanism, the low viscosity cases' plumes may not produce observable deformation. In an effort to eliminate the smallest plumes, we experimented with different lower bound parameters and temperature percentages.

Venus↗

A Suzaku X-ray Observation of One Orbit of the Supergiant Fast X-ray Transient IGR J16479-4514

We report on a 250 ks long X-ray observation of the supergiant fast X-ray transient (SFXT) IGR J16479-4514 performed with Suzaku in 2012 February. During this observation, about 80% of the short orbital period (P(sub orb) approximates 3.32 days) was covered as continuously as possible for the first time. The source light curve displays variability of more than two orders of magnitude, starting with a very low emission state (10(exp −13) erg / sq cm/s; 1-10 keV) lasting the first 46 ks, consistent with being due to the X-ray eclipse by the supergiant companion. The transition to the uneclipsed X-ray emission is energy dependent. Outside the eclipse, the source spends most of the time at a level of 6-7X10)(exp−12) erg/sq. cm/s) punctuated by two structured faint flares with a duration of about 10 and 15 ks, respectively, reaching a peak flux of 3-4X10(exp −11) erg/sq. cm./S, separated by about 0.2 in orbital phase. Remarkably, the first faint flare occurs at a similar orbital phase of the bright flares previously observed in the system. This indicates the presence of a phase-locked large scale structure in the supergiant wind, driving a higher accretion rate onto the compact object. The average X-ray spectrum is hard and highly absorbed, with a column density, NH, of 10*exp 23)/sq cm, clearly in excess of the interstellar absorption. There is no evidence for variability of the absorbing column density, except that during the eclipse, where a less absorbed X-ray spectrum is observed. A narrow Fe K-alpha emission line at 6.4 keV is viewed along the whole orbit, with an intensity which correlates with the continuum emission above 7 keV. The scattered component visible during the X-ray eclipse allowed us to directly probe the wind density at the orbital separation, resulting in rho(sub w)=7X10(exp −14) g/cubic cm. Assuming a spherical geometry for the supergiant wind, the derived wind density translates into a ratio M(sub w)/v(sub infinity) = 7X10(exp −17) Solar M/km which, assuming terminal velocities in a large range 500-3000 km/s, implies an accretion luminosity two orders of magnitude higher than that observed. As a consequence, a mechanism should be at work reducing the mass accretion rate. Different possibilities are discussed.

J164794514↗

3D Space Radiation Transport in a Shielded ICRU Tissue Sphere

A computationally efficient 3DHZETRN code capable of simulating High Charge (Z) and Energy (HZE) and light ions (including neutrons) under space-like boundary conditions with enhanced neutron and light ion propagation was recently developed for a simple homogeneous shield object. Monte Carlo benchmarks were used to verify the methodology in slab and spherical geometry, and the 3D corrections were shown to provide significant improvement over the straight-ahead approximation in some cases. In the present report, the new algorithms with well-defined convergence criteria are extended to inhomogeneous media within a shielded tissue slab and a shielded tissue sphere and tested against Monte Carlo simulation to verify the solution methods. The 3D corrections are again found to more accurately describe the neutron and light ion fluence spectra as compared to the straight-ahead approximation. These computationally efficient methods provide a basis for software capable of space shield analysis and optimization.

Wilson, John W.↗

A Study of Neutron Leakage in Finite Objects

A computationally efficient 3DHZETRN code capable of simulating High charge (Z) and Energy (HZE) and light ions (including neutrons) under space-like boundary conditions with enhanced neutron and light ion propagation was recently developed for simple shielded objects. Monte Carlo (MC) benchmarks were used to verify the 3DHZETRN methodology in slab and spherical geometry, and it was shown that 3DHZETRN agrees with MC codes to the degree that various MC codes agree among themselves. One limitation in the verification process is that all of the codes (3DHZETRN and three MC codes) utilize different nuclear models/databases. In the present report, the new algorithm, with well-defined convergence criteria, is used to quantify the neutron leakage from simple geometries to provide means of verifying 3D effects and to provide guidance for further code development.

Wilson, John W.↗