SEARCH · Search NASA
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.
Porosity in composite media using ultrasonic waves is characterized. The derivation of local approximations to the Kramers-Kronig relations are presented and it is shown that they may also be applicable to systems that could conceivably exhibit considerable dispersion such as composite laminates containing porosity.
A new dispersion relation for an unbounded degenerate electron plasma imbedded in a superstrong magnetic field is derived on the basis of the Canuto-Chiu equation of state. The properties of the various types of waves are discussed for two density regimes of astrophysical interest. Propagation parallel and perpendicular to the magnetic field are examined in detail.
A three dimensional fluid model which includes the dispersive effect of electron inertia is used to study the nonlinear macroscopic plasma dynamics of small scale discrete auroral arcs within the auroral acceleration zone and ionosphere. The motion of the Alfven wave source relative to the magnetospheric and ionospheric plasma forms an oblique Alfven wave which is reflected from the topside ionosphere by the negative density gradient. The superposition of the incident and reflected wave can be described by a steady state analytical solution of the model equations with the appropriate boundary conditions. This two dimensional discrete auroral arc equilibrium provides a simple explanation of auroral acceleration associated with the parallel electric field. Three dimensional fully nonlinear numerical simulations indicate that the equilibrium arc configuration evolves three dimensionally through collisionless tearing and reconnection of the current layer. The interaction of the perturbed flow and the transverse magnetic field produces complex transverse structure that may be the origin of the folds and curls observed to be associated with small scale discrete arcs.
Solar wind electrons are observed often to consist of two components: a core and a halo. The anisotropics and relative average speeds of these components correspond to a heat flux that has the potential to excite several different electromagetic instabilities; wave-particle scattering by the resulting enhanced fluctuations can limit this heat flux. This manuscript describes theoretical studies using the linear Vlasco dispersion equation for drifting bi-Maxwellian component distributions in a homogeneous plasma to examine the threshold of the whistler heat flux instability. Expressions for this threshold are obtained from two different parametric baselines: a local model that yields scalings as functions of local dimensionless plasma paramaters, and a global model based on average electron properties observed during the in-eliptic phase of the Ulysses mission. The latter model yields an expression for the heat flux at threshold of the whistler instability as a function of helisopheric radius that scales in the same way as the average heat flux observed form Ulysses and that provides an approximate upper bound for that same quantity. This theoretical scaling is combined with the observational results to yield a semi-empirical closure relation for the average electron heat flux in the solar wind between 1 and 5 AU.
We consider the stability of a circularly polarized Alfven wave (the pump wave) which propagates parallel to the ambient magnetic field. Only parallel-propagating perturbations are considered, and we ignore dispersive effects due to the ion cyclotron frequency. The dissipationless MHD equations are used throughout; thus possibibly important effects arising from Landau and transit time damping are omitted. We derive a series of analytical approximations to the dispersion relation using A = (Delta B/B(sub O))(exp 2) as a small expansion parameter; Delta B is the pump amplitude, and B(sub O) is the ambient magnetic field strength. We find that the plasma beta (the square of the ratio of the sound speed to the Alfven speed) plays a crucial role in determining the behavior of the parametric instabilities of the pump. If 0 less than beta less than 1 we find the familiar result that the pump decays into a forward propagating sound wave and a backward propagating Alfven wave with maximum growth rate gamma(sub max) varies A(sup 1/2), but beta cannot be too close to 0 or to 1. If beta approx. 1, we find gamma(sub max) varies A(sup 3/4), if beta greater than 1, we find gamma(sub max) varies A(sup 3/2), while if beta approx. 0, we obtain gamma(sub max) varies A(sup 1/3); moreover, if beta approx. 0 there is a nearly purely growing instability. In constrast to the familiar decay instability, for which the backward propagating Alfven wave has lower frequency and wavenumber than the pump, we find that if beta greater than or approx. equal to 1 the instability is really a beat instability which is dominated by a transverse wave which is forward propagating and has frequency and wavenumber which are nearly twice the pump values. Only the decay instability for 0 less than beta less than 1 can be regarded as producing two recognizable normal modes, namely, a sound wave and an Alfven wave. We discuss how the different characteristics of the instabilities may affect the evolution of Alfven waves in the solar wind. However, for a solar wind in which beta approx. 1 the growth times of the instabilities are probably too long for these instabilities to have an appreciable effect inside 1 AU.
To the first order of approximation, the complex amplitude of a wave packet in an anisotropic and dispersive medium is convected with the group of velocity. However, a gravity wave is a vector wave. Its wave packet must be formed by superposition of various wave numbers with corresponding frequencies, as is the case for scalar waves, and additionally by superposing many eigenmodes which also depend on the wave number. To represent the vector wave packet self-consistently, it is found that a gradient term must be included in the expansion. For a Guassian wave packet, this gradient term is shown to have important implications on the velocity vector as represented by its hodograph. Numerical results show that the hodograph is influenced by the location of the relative position of interest from the center of a Gaussian pulse. Higher order expansion shows that an initial Gaussian wave packet will retain its Gaussian shape as it propagates, but the pulse will spread in all directions with its major axis undergoing a rotation. Numerical results indicate that these higher order dispersive effects may be marginally observable in the atmosphere.
We reexamine both numerically and analytically the collapse of the singular isothermal sphere in the context of low-mass star formation. We consider the case where the onset of collapse is initiated by some arbitrary process which is accompanied by a central output of either heat or kinetic energy. We find two classes of numerical solutions describing this manner of collapse. The first approaches in time the expansion wave solution of Shu, while the second class is characterized by an ever-decreasing central accretion rate and the presence of an outwardly propagating weak shock. The collapse solution which represents the dividing case between these two classes is determined analytically by a similarity analysis. This solution shares with the expansion wave solution the properties that the gas remains stationary with an r(exp -2) density profile at large radius and that, at small radius, the gas free-falls onto a nascent core at a constant rate which depends only on the isothermal sound speed. This accretion rate is a factor of approx. 0.1 that predicted by the expansion wave solution. This reduction is due in part to the presence of a weak shock which propagates outward at 1.26 times the sound speed. Gas in the postshock region first moves out subsonically but is then decelerated and begins to collapse. The existence of two classes of numerical collapse solutions is explained in terms of the instability to radial perturbations of the analytic solution. Collapse occurring in the manner described by some of our solutions would eventually unbind a finite-sized core. However, this does not constitute a violation of the instability properties of the singular isothermal sphere which is unstable both to collapse and to expansion. To emphasize this, we consider a purely expanding solution for isothermal spheres. This solution is found to be self-similar and results in a uniform density core in the central regions of the gas. Our solutions may be relevant to the 'luminosity' problem of protostellar cores since the predicted central accretion rates are significantly reduced relative to that of the expansion wave solution. Furthermore, our calculations indicate that star-forming cloud cores are not very tightly bound and that modest disturbances can easily result in both termination of infall and dispersal of unaccreted material.
We reexamine both numerically and analytically the collapse of the singular isothermal sphere in the context of low-mass star formation. We consider the case where the onset of collapse is initiated by some arbitrary process which is accompanied by a central output of either heat or kinetic energy. We find two classes of numerical solutions describing this manner of collapse. The first approaches in time the expansion wave solution of Shu, while the second class is characterized by an ever-decreasing central accretion rate and the presence of an outwardly propagating weak shock. The collapse solution which represents the dividing case between these two classes is determined analytically by a similarity analysis. This solution shares with the expansion wave solution the properties that the gas remains stationary with an r(sup -2) density profile at large radius and that, at small radius, the gas free-falls onto a nascent core at a constant rate which depends only on the isothermal sound speed. This accretion rate is a factor of approx. 0.1 that predicted by the expansion wave solution. This reduction is due in part to the presence of a weak shock which propagates outward at 1.26 times the sound speed. Gas in the postshock region first moves out subsonically but is then decelerated and begins to collapse. The existence of two classes of numerical collapse solutions is explained in terms of the instability to radial perturbations of the analytic solution. Collapse occurring in the manner described by some of our solutions would eventually unbind a finite-sized core. However, this does not constitute a violation of the instability properties of the singular isothermal sphere which is unstable both to collapse and to expansion. To emphasize this, we consider a purely expanding solution for isothermal spheres. This solution is found to be self-similar and results in a uniform density core in the central regions of the gas. Our solutions may be relevant to the 'luminosity' problem of protostellar cores since the predicted central accretion rates are significantly reduced relative to that of the expansion wave solution. Furthermore, our calculations indicate that star-forming cloud cores are not very tightly bound and that modest disturbances can easily result in both termination of infall and dispersal of unaccreted material.
Dispersive Alfven wave activity is identified in four dayside auroral oval events measured by the Freja satellite. The events are characterized by ion injection, bursty electron precipitation below about I keV, transverse ion heating and broadband extremely low frequency (ELF) emissions below the lower hybrid cutoff frequency (a few kHz). The broadband emissions are observed to become more electrostatic towards higher frequencies. Large-scale density depletions/cavities, as determined by the Langmuir probe measurements, and strong electrostatic emissions are often observed simultaneously. A correlation study has been carried out between the E- and B-field fluctuations below 64 Hz (the dc instrument's upper threshold) and the characteristics of the precipitating electrons. This study revealed that the energization of electrons is indeed related to the broadband ELF emissions and that the electrostatic component plays a predominant role during very active magnetospheric conditions. Furthermore, the effect of the ELF electromagnetic emissions on the larger scale field-aligned current systems has been investigated, and it is found that such an effect cannot be detected. Instead, the Alfvenic activity creates a local region of field-aligned currents. It is suggested that dispersive Alfven waves set up these local field-aligned current regions and in turn trigger more electrostatic emissions during certain conditions. In these regions ions are transversely heated, and large-scale density depletions/cavities may be created during especially active periods.
This final report contains reports of research related to the tasks "Scalable High Performance Computing: Direct and Lark-Eddy Turbulent FLow Simulations Using Massively Parallel Computers" and "Devleop High-Performance Time-Domain Computational Electromagnetics Capability for RCS Prediction, Wave Propagation in Dispersive Media, and Dual-Use Applications. The discussion of Scalable High Performance Computing reports on three objectives: validate, access scalability, and apply two parallel flow solvers for three-dimensional Navier-Stokes flows; develop and validate a high-order parallel solver for Direct Numerical Simulations (DNS) and Large Eddy Simulation (LES) problems; and Investigate and develop a high-order Reynolds averaged Navier-Stokes turbulence model. The discussion of High-Performance Time-Domain Computational Electromagnetics reports on five objectives: enhancement of an electromagnetics code (CHARGE) to be able to effectively model antenna problems; utilize lessons learned in high-order/spectral solution of swirling 3D jets to apply to solving electromagnetics project; transition a high-order fluids code, FDL3DI, to be able to solve Maxwell's Equations using compact-differencing; develop and demonstrate improved radiation absorbing boundary conditions for high-order CEM; and extend high-order CEM solver to address variable material properties. The report also contains a review of work done by the systems engineer.
Wave measurements in planetary foreshocks and cometary environments show the sporadic occurrence of magnetic spectra with harmonic structure related to ion-cyclotron frequencies. Dilute populations of anisotropic and/or drifting charged particles can excite obliquely propagating modes with spacecraft frequencies close to the observed harmonics. Previous analyses of this generation mechanism are extended to drifting and nondrifting loss-cone-type distributions of heavy ions in a dense hydrogen magnetoplasma, characterizing the complex (real frequency and growth rate) dispersion, polarization, and compressibility of the unstable cyclotron harmonic waves. Solution of the full kinetic dispersion equation shows that it is possible to attain harmonic excitation, both in the drifting and nondrifting regimes. However, the bandwidth inherent to frequency Doppler shifts of obliquely propagating waves might preclude the observation of spectral structure in the spacecraft frame. The Giotto observations in the upstream region of comet Halley provide a reference to discuss the results.
Climate simulations by 16 atmospheric general circulation models (AGCMs) are compared on an aqua-planet, a water-covered Earth with prescribed sea surface temperature varying only in latitude. The idealised configuration is designed to expose differences in the circulation simulated by different models. Basic features of the aqua-planet climate are characterised by comparison with Earth. The models display a wide range of behaviour. The balanced component of the tropospheric mean flow, and mid-latitude eddy covariances subject to budget constraints, vary relatively little among the models. In contrast, differences in damping in the dynamical core strongly influence transient eddy amplitudes. Historical uncertainty in modelled lower stratospheric temperatures persists in APE.Aspects of the circulation generated more directly by interactions between the resolved fluid dynamics and parameterized moist processes vary greatly. The tropical Hadley circulation forms either a single or double inter-tropical convergence zone (ITCZ) at the equator, with large variations in mean precipitation. The equatorial wave spectrum shows a wide range of precipitation intensity and propagation characteristics. Kelvin mode-like eastward propagation with remarkably constant phase speed dominates in most models. Westward propagation, less dispersive than the equatorial Rossby modes, dominates in a few models or occurs within an eastward propagating envelope in others. The mean structure of the ITCZ is related to precipitation variability, consistent with previous studies.The aqua-planet global energy balance is unknown but the models produce a surprisingly large range of top of atmosphere global net flux, dominated by differences in shortwave reflection by clouds. A number of newly developed models, not optimised for Earth climate, contribute to this. Possible reasons for differences in the optimised models are discussed.The aqua-planet configuration is intended as one component of an experimental hierarchy used to evaluate AGCMs. This comparison does suggest that the range of model behaviour could be better understood and reduced in conjunction with Earth climate simulations. Controlled experimentation is required to explore individual model behavior and investigate convergence of the aqua-planet climate with increasing resolution.
Spacecraft observations in the earth's magnetotail at distances of 30 to 40 R(E) have revealed the presence of broadband electrostatic waves. These waves are generally most intense in the regions just outside of the plasma sheet and are correlated with the observations of relatively cold and energetic ion beams traveling in either the earthward or the tailward direction. These waves are observed to propagate obliquely to the geomagnetic field with wave normal angles around 70 deg. Because the broadband electrostatic noise is the most intense of the waves observed in the magnetotail, it is important to understand the generation mechanism of these waves. The purpose of this study is to provide for the first time a correct solution to the dispersion equation for ion beams observed in the magnetotail. By numerically solving this equation, it is shown that obliquely propagating waves have growth rates that can be an order of magnitude larger than those of parallel propagating waves, in agreement with observations. In addition, the effect of beam temperature on the ion beam instability is studied, and it is shown that this instability can be a viable generation mechanism only when the ion beam has a relatively small thermal spread.
The dispersion equation is obtained analytically of waves along an infinitely long cylindrical ?grid-like' antenna in an anisotropic compressible homogeneous plasma. Numerical solutions of this equation show that an external static magnetic field has great influence upon the anomalously dispersive region of the dispersion relation around the plasma frequency.
A biaxial graded-index fiber with a homogeneous cladding is studied. Two methods, wave equation and matrix differential equation, of formulating the problem and their respective solutions are discussed. For the wave equation formulation of the problem it is shown that for the case of a diagonal permittivity tensor the longitudinal electric and magnetic fields satisfy a pair of coupled second-order differential equations. Also, a generalized dispersion relation is derived in terms of the solutions for the longitudinal electric and magnetic fields. For the case of a step-index fiber, either isotropic or uniaxial, these differential equations can be solved exactly in terms of Bessel functions. For the cases of an istropic graded-index and a uniaxial graded-index fiber, a solution using the Wentzel, Krammers and Brillouin (WKB) approximation technique is shown. Results for some particular permittivity profiles are presented. Also the WKB solutions is compared with the vector solution found by Kurtz and Streifer. For the matrix formulation it is shown that the tangential components of the electric and magnetic fields satisfy a system of four first-order differential equations which can be conveniently written in matrix form. For the special case of meridional modes, the system of equations splits into two systems of two equations. A general iterative technique, asymptotic partitioning of systems of equations, for solving systems of differential equations is presented. As a simple example, Bessel's differential equation is written in matrix form and is solved using this asymptotic technique. Low order solutions for particular examples of a biaxial and uniaxial graded-index fiber are presented. Finally numerical results obtained using the asymptotic technique are presented for particular examples of isotropic and uniaxial step-index fibers and isotropic, uniaxial and biaxial graded-index fibers.
We treat RR flux backgrounds of type II string theory in the framework of closed superstring field theory based on the NSR formalism, focusing on two examples: (1) the pp-wave background supported by 5-form flux, and (2) AdS 3 × S 3 × M 4 supported by mixed 3-form fluxes. In both cases, we analyze the classical string field solution perturbatively, and compute the correction to the dispersion relation of string states to quadratic order in the RR flux. In the first example, our result is in a delicate way consistent with that obtained from lightcone quantization of the Green-Schwarz string. In the second example, we will obtain numerically the mass corrections to pulsating type IIB strings in AdS 3 × S 3 × M 4 . Our results, valid at finite AdS radius, agree with previously known answers in the semiclassical limit and in the BMN limit respectively.
The physical mechanisms which potentially lead to instabilities in line driven winds, the drift instability and the line shape instability, are discussed. A general three dimensional treatment of the stability problem of line driven winds which leads to the general dispersion equation is proposed. From this dispersion equation automatically a third physical mechanism driving instability in stellar winds is deduced; the thermal drift instability which is related to changes in absorption of radiation caused by temperature perturbations. This mechanism results in growing inwardly propagating sound waves.
The burning of liquid propellants is a fundamental combustion problem that is applicable to various types of propulsion and energetic systems. The deflagration process is often rather complex, with vaporization and pyrolysis occurring at the liquid/gas interface and distributed combustion occurring either in the gas phase or in a spray. Nonetheless, there are realistic limiting cases in which combustion may be approximated by an overall reaction at the liquid/gas interface. In one such limit, the gas flame occurs under near-breakaway conditions, exerting little thermal or hydrodynamic influence on the burning propellant. In another such limit, distributed combustion occurs in an intrusive regime, the reaction zone lying closer to the liquid/gas interface than the length scale of any disturbance of interest. Finally, the liquid propellant may simply undergo exothermic decomposition at the surface without any significant distributed combustion, such as appears to occur in some types of HydroxylAmmonium Nitrate (HAN)-based liquid propellants at low pressures. Such limiting models have recently been formulated,thereby significantly generalizing earlier classical models that were originally introduced to study the hydrodynamic stability of a reactive liquid/gas interface. In all of these investigations, gravity appears explicitly and plays a significant role, along with surface tension, viscosity, and, in the more recent models, certain reaction-rate parameters associated with the pressure and temperature sensitivities of the reaction itself. In particular, these parameters determine the stability of the deflagration with respect to not only classical hydrodynamic disturbances, but also with respect to reactive/diffusive influences as well. Indeed, the inverse Froude number, representing the ratio of buoyant to inertial forces, appears explicitly in all of these models, and consequently, in the dispersion relation that determines the neutral stability boundaries beyond which steady, planar burning is unstable to nonsteady, and/or nonplanar (cellular) modes of burning. These instabilities thus lead to a number of interesting phenomena, such as the sloshing type of waves that have been observed in mixtures of HAN and TriEthanolAmmonium Nitrate (TEAN) with water. Although the Froude number was treated as an O(1) quantity in these studies, the limit of small inverse Froude number corresponding to the microgravity regime is increasingly of interest and can be treated explicitly, leading to various limiting forms of the models, the neutral stability boundaries, and, ultimately, the evolution equations that govern the nonlinear dynamics of the propagating reaction front. In the present work, we formally exploit this limiting parameter regime to compare some of the features of hydrodynamic instability of liquid-propellant combustion at reduced gravity with the same phenomenon at normal gravity.