SEARCH · Engineering Papers
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.
A boundary-layer solution to a high magnetic Reynolds number R periodic dynamo model shows that: (1) flux expulsion forces the magnetic field into flux sheets; (2) the principal contribution to the alpha effect arises from regions of flow stagnation along a flux sheet; and (3) the alpha effect scales as R exp-1/2. Arguments for these effects persisting in turbulent dynamos are given.
Thermal instability is believed to determine the evolution and formation of cool structures in the solar atmosphere such as the transition region and prominences (or filaments). The linear modes that arise in a sheared, force-free, magnetic field due to thermal instability are studied numerically. Previous studies have considered separately modes that arise due to the effects of radiation, compression, anisotropic thermal conduction, and ohmic heating. Here the results of such studies are integrated, first by presenting simple arguments that illustrate the essential physics of ideal, sheared-field, condensation modes, and second by showing numerically how finite resistivity affects the condensational instability in parameter regimes applicable to the solar corona.
The Geodynamics Laser Ranging System (GLRS) is a spaceborne laser ranging instrument being developoed by NASA as a facility instrument for the Earth Observing System (EOS). GLRS is to be used to study regional and local crustal movements. It is designed to make highly precise range measurements to retroreflector targets located in geophysically interesting sites. Using a two-color ranging scheme, absolute range accuracies of several mm are expected. Simulations based on this accuracy and the EOS orbital parameters show that length of the intersite baseline between retroreflectors can be determined to several mm accuracy at distances from a few km to several hundred km with several passes of GLRS range data collected over a few-day interval. Short-arc techniques are used to minimize the effects of gravity field and other force model uncertainties. Relative heights can be determined to sub-cm accuracy over comparable distances.
Numerical solutions of the compressible MHD equations are used here to simulate the evolution of an initially force-free magnetic field in a static corona as a result of slow photospheric motion of the magnetic field footpoints. Simulations have been completed for values of plasma beta from 0.1 to 0.5, maximum shear velocities from 0.5 to 10.3 km/s, and with various amounts of resistive and viscous dissipation. In all cases the evolution proceeds in two qualitatively different stages. In the earlier stage, the field evolves gradually with the field lines, expanding outward at a velocity not unlike the shear velocity. Then, the field begins to expand much more rapidly until it reaches velocities exceeding a characteristic Alfven velocity. Inclusion of the thermodynamics, gravity, and compressibility is shown to have only a quantitative effect on the onset of the eruptive phase, illustrating that the primary interactions are between the dynamics and the magnetic field evolution.
We present a comparison between two numerical methods for the extrapolation of nonlinear force-free magnetic fields, the Iterative Method (IM) and the Progressive Extension Method (PEM). The advantages and disadvantages of these two methods are summarized and the accuracy and numerical instability are discussed. On the basis of this investigation, we claim that the two methods do resemble each other qualitatively.
The aim of the research activity was to increase our understanding of solar activity through data analysis, theoretical analysis, and computer modeling. Because the research subjects were diverse and many researchers were supported by this grant, a select few key areas of research are described in detail. Areas of research include: (1) energy storage and force-free magnetic field; (2) energy release and particle acceleration; (3) radiation by nonthermal electrons; (4) coronal loops; (5) flare classification; (6) longitude distributions of flares; (7) periodicities detected in the solar activity; (8) coronal heating and related problems; and (9) plasma processes.
The theory of axially symmetric self-focusing of whistler waves, based on the full system of Maxwell equations, is developed. The plasma is described by the magnetohydrodynamic equations including the ponderomotive force from RF field. The nonlinear Schrodinger equations (NSE) for arbitrary azimuthal modes of whistler waves are derived. It is shown that they differ from the NSE for a scalar field; this is connected with an intrinsic angular momentum due to the rotating polarization of whistlers. It is shown that the self-focusing, as described by the NSE, differs in its final stage from the results following the full set of Maxwell equations. The latter gives defocusing after sufficient narrowing of the initial wave beam, due to transformation of the trapped wave into a nontrapped branch which is not contained in the NSE description. The oscillatory character of the defocusing is demonstrated.
When a plasma is pushed across a magnetic field by some nonelectromagnetic force, ions and electrons get turned in opposite directions by the magnetic field. An exact analysis of that process is presented here for the internal region of the plasma. The energy provided by the initial push is used, in part, to create the electric field and in part to create some gyrations inside the plasma. When the rest energy density of the plasma exceeds twice the magnetic energy density (or when the Alfven speed is less than c), there will be enough energy to spare for the plasma to continue across the magnetic field at half its initial momentum. Two cases are considered: an impulsive start and a gentle push such as provided by gravity. The amplitude of the resulting internal gyrations becomes small in the second case. The frequencies of the gyrations are those of extraordinary modes of very long spatial wavelength.
We have constructed a new class of exact, nonlinear solutions to the Grad-Shafranov equation, representing force-free magnetic fields with translational symmetry. These exact solutions are pertinent to the study of magnetic structures in the solar corona that are subjected to photospheric shearing motions.
The interstellar medium (ISM) is inhomogeneous, with clouds of various temperatures and densities embedded in a tenuous intercloud medium. Shocks propagating through the ISM can ablate or destroy the clouds, at the same time significantly altering the properties of the intercloud medium. This paper presents a comprehensive numerical study of the simplest case of the interaction between a shock wave and a spherical cloud, in which the shock far from the cloud is steady and planar, and in which radiative losses, thermal conduction, magnetic fields, and gravitational forces are all neglected. As a result, the problem is completely specified by two numbers: the Mach number of the shock, M, and the ratio of the density of the cloud to that of the intercloud medium, Chi. For strong shocks we show that the dependence on M scales out, so the primary independent parameter is Chi. Variations from this simple case are also considered: the potential effect of radiative losses is assessed by calculations in which the ratio of specific heats in the cloud is 1.1 instead of 5/3; the effect of the initial shape of the cloud is studied by using a cylindrical cloud instead of a spherical one; and the role of the initial shock is determined by considering the case of a cloud embedded in a wind. Local adaptive mesh refinement techniques with a second-order, two-fluid, two-dimensional Godunov hydrodynamic scheme are used to address these problems, allowing heretofore unobtainable numerical resolution. Convergence studies to be described in a subsequent paper demonstrate that about 100 zones per cloud radius are needed for accurate results; previous calculations have generally used about a third of this number. The results of the calculations are analyzed in terms of global quantities which provide an overall description of te shocked cloud: the size and shape of the cloud, the mean density, the mean pressure, the mean velocity, the velocity dispersion, and the total circulation.
This paper outlines ISMAN suggestions for the joint NASA-RSA project 'Combustion and Structure formation in SHS Processes under Microgravity Conditions'. The basic ideas of this work naturally follow from our almost 30-year experience in the field of SHS. As a matter of fact, we have already obtained some results in the following two directions closely related to the microgravity problem. One is the studies on SHS processes in the field of centrifugal forces. These studies aimed at the intensification of gravity-sensitive SHS processes in multicomponent highly caloric systems forming melts at high overloads (up to 2000 g). In other words, these studies had the objectives that are inverse to those in the microgravity studies. The second group of results directly relates to the microgravity problem and the project under consideration. These experiments played the important role in establishing links between SHS and microgravity.
A three dimensional, unstructured-mesh methodology was developed to simulate unsteady flows past bodies in relative motion, where the trajectory was determined from the instantaneous aerodynamics. The method coupled the equations of fluid flow and those of rigid-body dynamics, and captured the time-dependent interference between stationary and moving boundaries. The unsteady, compressible Euler equations were solved on dynamic, unstructured meshes by an explicit, finite-volume, upwind method. The grid adaptation was performed within a window placed around the moving body. The Euler equations of dynamics were solved by a Runge-Kutta integration scheme. The flow solver and the adaptation scheme were validated by simulating the transonic, unsteady flow around a wing undergoing a forced, periodic pitching motion, then comparing the results with the experimental data. To validate the trajectory code, the six-degrees-of-freedom (DOF) motion of a store separating from a wing was computed using the experimentally determined force and moment fields, then comparing with an independently generated trajectory. Finally, the overall methodology was demonstrated by simulating the unsteady flowfield and the trajectory of a store dropped from a wing. The methodology, its computational cost notwithstanding, has proven to be accurate, automated, easy for dynamic gridding, and relatively efficient for the required man-hours.
A new methodology is developed to simulate unsteady flows about prescribed and aerodynamically determined moving boundary problems. The method couples the fluid dynamics and rigid-body dynamics equations to capture the time-dependent interference between stationary and moving boundaries. The unsteady, compressible, inviscid (Euler) equations are solved on dynamic, unstructured grids by an explicit, finite-volume, upwind method. For efficiency, the grid adaptation is performed within a window around the moving object. The Eulerian equations of the rigid-body dynamics are solved by a Runge-Kutta method in a non-inertial frame of reference. The two-dimensional flow solver is validated by computing the flow past a sinusoidally-pitching airfoil and comparing these results with the experimental data. The overall methodology is used for two two-dimensional examples: the flow past an airfoil which is performing a three-degrees-of-freedom motion in a transonic freestream, and the free-fall of a store after separation from a wing-section. Then the unstructured mesh methodology is extended to three-dimensions to simulate unsteady flow past bodies in relative motion, where the trajectory is determined from the instantaneous aerodynamics. The flow solver and the adaptation scheme in three dimensions are validated by simulating the transonic, unsteady flow around a wing undergoing a forced, periodic, pitching motion, and comparing the results with the experimental data. To validate the trajectory code, the six-degrees-of-freedom motion of a store separating from a wing was computed using the experimentally determined force and moment fields, then comparing with an independently generated trajectory. Finally, the overall methodology was demonstrated by simulating the unsteady flowfield and the trajectory of a store dropped from a wing. The methodology, its computational cost notwithstanding, has proven to be accurate, automated, easy for dynamic gridding, and relatively efficient for the required man-hours.
This research program involves the study of single droplets and linear arrays of droplets in weakly-buoyant and non-bouyant environments. The primary purpose of the single droplet work was to (1) provide a data base from which to compare droplet array results and (2) to correlate the effects of buoyancy on flame shape. Traditionally convective effects in droplet combustion are represented in terms of the Reynolds number, Re, for forced convection and the Grashof number, Gr, for natural convection. Typically, corrections to the burning rate constant for convective effects are written in terms of Re or Gr(exp 1). The Stefan velocity is not included in these correlations, even though from purely physical reasons, one would expect it to be important, especially at higher burning rates. The flame distortion due to convective effects is less documented quantitatively. Kumagai and Isoda do predict flame shape in natural and forced convective flow fields. Their focus, however, was to predict the actual flame dimensions. Law and co-workers used reduced pressure, high oxidizer ambients to obtain spherical flames. This implies that buoyant flows were reduced at the low pressures, as indicated by a very small Grashof number. Ross et al, however, using scaling arguments showed that reducing the pressure does not have a large effect on the magnitude of the buoyant velocity. Struk et al showed elongated flame shapes during simulated (porous sphere) droplet combustion. The elongation of the flames was due to residual gravity levels aboard the reduced gravity aircraft on which the experiments were conducted. These flame shapes, as well as some data from the literature were interpreted based on a dimensionless grouping called the sphericity parameter, Sp. Sp is the ratio of a characteristic computed buoyant velocity to the Stefan velocity at the flame front. One purpose of the droplet arrays work is to extend the database and theories that exist for single droplets into the regime where droplet interactions are important. The eventual goal being to use the results of this work as inputs to models on spray combustion where droplets seldom burn individually; instead the combustion history of a droplet is strongly influenced by the presence of the neighboring droplets. Recently, Annamali and Ryan have summarized he current status of droplet array, cloud and spray combustion. A number of simplified theories led numerical studies of droplet vaporization/combustion where multiple droplet effects are present are now available. These theories all neglect the effect of buoyancy. Experimentally, most studies to date suffer the effects of buoyancy. It is the dominant transport mechanism in the problem. Only the works of Law and co-worker and more recently by Mikami et al were performed in an environment where buoyancy effects were small. Law and co-workers were limited to high oxygen index, low pressure ambient environments since there studies were conducted in normal gravity.
The mission of this research is to be world-class creators and facilitators of innovative, intelligent, high performance, reliable information technologies that enable NASA missions to (1) increase software safety and quality through error avoidance, early detection and resolution of errors, by utilizing and applying empirically based software engineering best practices; (2) ensure customer software risks are identified and/or that requirements are met and/or exceeded; (3) research, develop, apply, verify, and publish software technologies for competitive advantage and the advancement of science; and (4) facilitate the transfer of science and engineering data, methods, and practices to NASA, educational institutions, state agencies, and commercial organizations. The goals are to become a national Center Of Excellence (COE) in software and system independent verification and validation, and to become an international leading force in the field of software engineering for improving the safety, quality, reliability, and cost performance of software systems. This project addresses the following problems: Ensure safety of NASA missions, ensure requirements are met, minimize programmatic and technological risks of software development and operations, improve software quality, reduce costs and time to delivery, and improve the science of software engineering
The purpose of this project was to develop rocket-borne probes to detect charged aerosol layers in the mesosphere. These include sporadic E layers, which have their origin in meteoric dust, and noctilucent clouds, which form in the arctic summer and are composed of ice crystals. The probe being developed consists of a charge collecting patch connected to a sensitive electrometer which measures the charge deposited on the patch by impacting aerosols. The ambient electrons and light ions in the mesosphere are prevented from being collected by a magnetic field. The magnetic force causes these lighter particles to turn so that they miss the collecting patch.
The University of Wisconsin and Louisiana State University have teamed to study the forcing of winter season cold frontal wind systems on sediment distribution patterns and geomorphology in the Louisiana coastal zone. Wind systems associated with cold fronts have been shown to modify coastal circulation and resuspend sediments along the microtidal Louisiana coast. The assessment includes quantifying the influence of cumulative winter season atmospheric forcing (through surface wind observations) from year to year in response to short term climate variability, such as El Nino events. A correlation between winter cyclone frequency and the strength of El Nino events has been suggested. The atmospheric forcing data are being correlated to geomorphic measurements along western Louisiana's prograding muddy coast. Remote sensing data is being used to map and track sediment distribution patterns for various wind conditions. Transferring a suspended sediment concentration (SSC) algorithm to EOS MODIS observations will enable estimates of SSC in case 2 waters over the global domain. Progress in Year 1 of this study has included data collection and analysis of wind observations for atmospheric forcing characterization, a field activity (TX-2001) to collect in situ water samples with co-incident remote sensing measurements from the NASA ER-2 based MODIS Airborne Simulator (MAS) and the EOS Terra based MODerate resolution Imaging Spectroradiometer (MODIS) instruments, aerial photography and of sediment burial pipe field measurements along the prograding muddy Chenier Plain coast of western Louisiana for documenting coastal change in that dynamic region, and routine collection of MODIS 250 in resolution data for monitoring coastal sediment patterns. The data sets are being used in a process to transfer an SSC estimation algorithm to the MODIS platform. Work is underway on assessing coastal transport for the winter 2000-01 season. Water level data for use in a Geomorphic Impact Index, which relates wind energy, water level conditions, and geomorphic change along the microtidal western Louisiana coastline is being assembled.
Our research is in the area of solid-state organic chemistry, which lies at the interface between physical organic chemistry and materials science. We use crystalline solids as models to probe fundamental issues about physical processes, molecular interactions and chemical reactions that are important for fabrication, stabilization and application of technological materials. Much of our most recent work has focused on the phenomena of ferroelastic and ferroelectric domain switching, in which application of an external force or electric field to a crystal causes the molecules inside the crystal to reorient, in tandem, to a new orientational state. To better understand and control the domain switching process, we have designed and synthesized over twenty closely related, ferroelastic organic crystals. Our approach has been to use crystalline inclusion compounds, in which one molecule (the guest) is trapped within the crystalline framework of a second molecule (the host). By keeping the host constant and varying the proportions and kinds of guests, it has been possible to tailor these materials so that domain switching is rapid and reversible (which is desirable for high technology applications). Inclusion compounds therefore serve as powerful systems for understanding the specific molecular mechanisms that control domain switching.