Search NASA⌕ Search

SEARCH · Search NASA

Results for “dispersive equations”

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 163 records · Page 9

Analysis of Kelvin Helmholtz Instabilities of Plasma Jets

Ulysses data indicates density fluctuations which axe theorized to be the result of shear between a solar jet and its ambient. The MHD Kelvin-Helmholtz ("KH") instability causes such fluctuations as observed by Ulysses. A new dispersion relationship which accounts for this KH instability is derived via the linearization of the MHD equations. This generalizes an earlier result by Hardee. This dispersion relationship has the form of eight non-linear equations with nine unknowns.

Friedman, Mark J.↗

Transformation of apparent ocean wave spectra observed from an aircraft sensor platform

The problem considered was transformation of a unidirectional apparent ocean wave spectrum observed from an aircraft sensor platform into the true spectrum that would be observed from a stationary platform. Spectral transformation equations were developed in terms of the linear wave dispersion relationship and the wave group speed. An iterative solution to the equations was outlined and used to transform reference theoretical apparent spectra for several assumed values of average water depth. Results show that changing the average water depth leads to a redistribution of energy density among the various frequency bands of the transformed spectrum. This redistribution is most severe when much of the energy density is expected, a priori, to reside at relatively low true frequencies.

Poole, L. R.↗

Direct Measurements of Eddy Transport and Thermal Dispersion in a High Porosity Matrix

Thermal losses from the hot end to the cold end of a Stirling cycle regenerator due to thermal dispersion through the regenerator matrix may significantly degrade the performance of the machine. Because of poor access to the void spaces within the porous medium, no direct measurements of thermal dispersion have been made and dispersion models have been derived indirectly. This is done by measuring the overall thermal performance of the regenerator and subtracting off the energy transfer caused by molecular conduction and advected enthalpy flows computed from volume-averaged fluid velocity and temperature. In the current program, a large-scale porous matrix consisting of stacked wire screens with a porosity of 90% is installed in a flow rig which is operated in a Reynolds number range that represents Stirling engine regenerator flow. Experiments are conducted to measure turbulent transport of momentum at the exit phase using hot-wires. The relationship of such turbulent transport terms to the thermal dispersion term in the volumetric-averaged energy equation for the regenerator matrix is developed and the measurements are used to determine cross-stream thermal dispersion. A dispersion model based upon the measurements is proposed and compared with models documented in the literature.

Niu, Yi↗

Dispersive waves in a seeded MHD generator.

The equations giving the response of a slightly ionized plasma with monatomic components to sinusoidal perturbations have been formulated. Included in the model equations were the electron Hall effect, electron thermal diffusion, radiation, and electron-atom rate processes. Plasma conditions were limited to those where viscous effects, the induced magnetic field, ion slip, and atom-atom inelastic processes can be neglected. Presented are results of numerical calculations for MHD generators with a working fluid of potassium seeded argon.

Harstad, K. G.↗

Vorticity Transport and Wave Emission in the Protoplanetary Nebula

Higher order numerical algorithms (4th order in time, 3rd order in space) are applied to the Euler/Energy equations and are used to examine vorticity transport and wave motion in a non-self gravitating, initially isentropic Keplerian disk. In this talk we will examine the response of the nebula to an isolated vortex with a circulation about equal to the rotation rate of Jupiter. The vortex is located on the 4 AU circle and the nebula is simulated from 1 to 24 AU. We show that the vortex emits pressure-supported density and Rossby-type wave packets before it decays within a few orbits. The acoustic density waves evolve into weak (non entropy preserving) shock waves that propagate over the entire disk. The Rossby waves remain in the vicinity of the initial vortex disturbance, but are rapidly damped. Temporal frequencies and spatial wavenumbers are derived using the simulation data and compared with analytical dispersion relations from the linearized Euler/Energy equations.

Davis, S. S.↗

Rheological properties of nondilute suspensions of deformable particles

The creeping flow equations have been solved to account for effects of shear on drop shape when the influence of neighboring drops is considered. A cell model is used to account for finite concentrations of the dispersed phase. From these results a constitutive equation is derived and implications for several flows are discussed. The shape and orientation of drops are computed and compared with existing data. Finally, dynamic viscosity and rigidity are obtained for nondilute suspensions of viscoelastic particles in an oscillatory shear flow.-

Choi, S. J.↗

On the dynamics of agglomerating ensembles of clouds

It is shown that for dynamical time scales longer than 10 million years the ensemble of interstellar cold diffuse clouds may be described by equations of mass and momentum conservation where the pressure has the form p/v-squared/ where /v-squared/ is the mass averaged square velocity dispersion. Equilibrium cloud spectra are derived and an equation of state obtained. The equation of state undergoes a phase transition at a critical density close to the local interstellar gas density. At higher densities, more massive clouds are formed and the cloud fluid softens. It is suggested that this mechanism may naturally regulate the interstellar gas density if star formation is very efficient in the cooler high-density systems.

Cowie, L. L.↗

Numerical Boundary Conditions for Computational Aeroacoustics Benchmark Problems

Category 1, Problems 1 and 2, Category 2, Problem 2, and Category 3, Problem 2 are solved computationally using the Dispersion-Relation-Preserving (DRP) scheme. All these problems are governed by the linearized Euler equations. The resolution requirements of the DRP scheme for maintaining low numerical dispersion and dissipation as well as accurate wave speeds in solving the linearized Euler equations are now well understood. As long as 8 or more mesh points per wavelength is employed in the numerical computation, high quality results are assured. For the first three categories of benchmark problems, therefore, the real challenge is to develop high quality numerical boundary conditions. For Category 1, Problems 1 and 2, it is the curved wall boundary conditions. For Category 2, Problem 2, it is the internal radiation boundary conditions inside the duct. For Category 3, Problem 2, they are the inflow and outflow boundary conditions upstream and downstream of the blade row. These are the foci of the present investigation. Special nonhomogeneous radiation boundary conditions that generate the incoming disturbances and at the same time allow the outgoing reflected or scattered acoustic disturbances to leave the computation domain without significant reflection are developed. Numerical results based on these boundary conditions are provided.

Tam, Chritsopher K. W.↗

Low gravity transfer line chilldown

A code has been developed that solves for the transfer line chilldown time and flow and heat transfer characteristics in one-g environment. The code solves the transient, one dimensional, space averaged mass, momentum, and energy conservation equations for liquid-vapor two-phase flow in tubes. The physical configuration solved is that appropriate for bottom coolant injection in a vertically supported heated tube. Four distinct regions are considered consecutively: fully liquid, inverted annular, dispersed, and fully vapor flow. The conservation equations for both the liquid and the vapor are solved in each region separately. Also, in each region the mass and energy transport between each phase as well as the energy and momentum transport between the tube wall and the fluid are accounted for. A finite wall thickness is also considered. The model described above was solved numerically through a mixed finite difference scheme with forward time marching. The inverted annular regime was resolved using a semi-implicit finite differencing while the dispersed regime was solved explicitly. Also, a staggered mesh was used in which the velocity was resolved at mesh boundaries while all other field variables were resolved at the mesh centroids. Different mesh sizes were used depending on the region of solution. A coarse mesh was used in the dispersed flow region while a much finer mesh was used in both the inverted annular flow region and the tube walls.

B. N. Antar↗

Low gravity transfer line chilldown

A code has been developed that solves for the transfer line chilldown time and flow and heat transfer characteristics in one-g environment. The code solves the transient, one dimensional, space averaged mass, momentum and energy conservation equations for liquid-vapor two-phase flow in tubes. The physical configuration solved is that appropriate for bottom coolant injection in a vertically supported heated tube. Four distinct regions are considered consecutively: fully liquid. inverted annular. dispersed and fully vapor flow. The conservation equations for both the liquid and the vapor are solved in each region separately. Also, in each region the mass and energy transport between each phase as well as the energy and momentum transport between the tube wall and the fluid are accounted for. A finite wall thickness is also considered.

Antar, B. N.↗

Computational modeling of femtosecond optical solitons from Maxwell's equations

An algorithm is developed that permits the direct time integration of full-vector nonlinear Maxwell's equations. This capability permits the modeling of both linear and nonlinear instantaneous and dispersive effects in the electric polarization in material media. The modeling of the optical carrier is retained. The fundamental innovation is to notice that it is possible to treat the linear and nonlinear convolution integrals, which describe the dispersion, as new dependent variables. A coupled system of nonlinear second-order ordinary differential equations can then be derived for the linear and nonlinear convolution integrals, by differentiating them in the time domain. These equations, together with Maxwell's equations, are solved to determine the electromagnetic fields in nonlinear dispersive media. Results are presented of calculations in one dimension of the propagation and collision of femtosecond electromagnetic solitons that retain the optical carrier, taking into account as the Kerr and Raman interactions.

Goorjian, Peter M.↗

Ultra-Broad-Band Optical Parametric Amplifier or Oscillator

A concept for an ultra-broad-band optical parametric amplifier or oscillator has emerged as a by-product of a theoretical study in fundamental quantum optics. The study was originally intended to address the question of whether the two-photon temporal correlation function of light [in particular, light produced by spontaneous parametric down conversion (SPDC)] can be considerably narrower than the inverse of the spectral width (bandwidth) of the light. The answer to the question was found to be negative. More specifically, on the basis of the universal integral relations between the quantum two-photon temporal correlation and the classical spectrum of light, it was found that the lower limit of two-photon correlation time is set approximately by the inverse of the bandwidth. The mathematical solution for the minimum two-photon correlation time also provides the minimum relative frequency dispersion of the down-converted light components; in turn, the minimum relative frequency dispersion translates to the maximum bandwidth, which is important for the design of an ultra-broad-band optical parametric oscillator or amplifier. In the study, results of an analysis of the general integral relations were applied in the case of an optically nonlinear, frequency-dispersive crystal in which SPDC produces collinear photons. Equations were found for the crystal orientation and pump wavelength, specific for each parametric-down-converting crystal, that eliminate the relative frequency dispersion of collinear degenerate (equal-frequency) signal and idler components up to the fourth order in the frequency-detuning parameter

Strekalov, Dmitry↗

Lagrangian description of warm plasmas

Efforts are described to extend the averaged Lagrangian method of describing small signal wave propagation and nonlinear wave interaction, developed by earlier workers for cold plasmas, to the more general conditions of warm collisionless plasmas, and to demonstrate particularly the effectiveness of the method in analyzing wave-wave interactions. The theory is developed for both the microscopic description and the hydrodynamic approximation to plasma behavior. First, a microscopic Lagrangian is formulated rigorously, and expanded in terms of perturbations about equilibrium. Two methods are then described for deriving a hydrodynamic Lagrangian. In the first of these, the Lagrangian is obtained by velocity integration of the exact microscopic Lagrangian. In the second, the expanded hydrodynamic Lagrangian is obtained directly from the expanded microscopic Lagrangian. As applications of the microscopic Lagrangian, the small-signal dispersion relations and the coupled mode equations are derived for all possible waves in a warm infinite, weakly inhomogeneous magnetoplasma, and their interactions are examined.

Kim, H.↗

Comparison of several finite-difference methods

The present numerical study establishes that, with the MacCormack and 'two-four' methods, the solution accuracy obtainable in problems involving wave propagation, shock-wave and contact discontinuities, and viscous effects, will be strongly dependent on Courant number. The application of flux correction to the MacCormack and two-four methods is noted to significantly attenuate dispersion errors; the ensuing solutions capture the discontinuities in the shock-tube problem with improved accuracy and resolution, and are free of dispersion errors for the viscous Burgers' equation.

Biringen, S.↗

Controlled Quantum Packets

We look at time evolution of a physical system from the point of view of dynamical control theory. Normally we solve motion equation with a given external potential and we obtain time evolution. Standard examples are the trajectories in classical mechanics or the wave functions in Quantum Mechanics. In the control theory, we have the configurational variables of a physical system, we choose a velocity field and with a suited strategy we force the physical system to have a well defined evolution. The evolution of the system is the 'premium' that the controller receives if he has adopted the right strategy. The strategy is given by well suited laboratory devices. The control mechanisms are in many cases non linear; it is necessary, namely, a feedback mechanism to retain in time the selected evolution. Our aim is to introduce a scheme to obtain Quantum wave packets by control theory. The program is to choose the characteristics of a packet, that is, the equation of evolution for its centre and a controlled dispersion, and to give a building scheme from some initial state (for example a solution of stationary Schroedinger equation). It seems natural in this view to use stochastic approach to Quantum Mechanics, that is, Stochastic Mechanics [S.M.]. It is a quantization scheme different from ordinary ones only formally. This approach introduces in quantum theory the whole mathematical apparatus of stochastic control theory. Stochastic Mechanics, in our view, is more intuitive when we want to study all the classical-like problems. We apply our scheme to build two classes of quantum packets both derived generalizing some properties of coherent states.

DeMartino, Salvatore↗

Failure Mode Analysis of V-Shaped Pyrotechnically Actuated Valves

Current V-shaped stainless steel pyrovalve initiators have rectified many of the deficiencies of the heritage Y-shaped aluminum design. However, a credible failure mode still exists for dual simultaneous initiator (NSI) firings in which low temperatures were detected at the booster cap and less consistent ignition was observed than when a single initiator was fired. In order to asses this issue, a numerical framework has been developed for predicting the flow through pyrotechnically actuated valves. This framework includes a fully coupled solution of the gas-phase equation with a non-equilibrium dispersed phase for solid particles as well as the capability to model conjugate gradient heat transfer to the booster cap. Through a hierarchy of increasingly complex simulations, a hypothesis for the failure mode of the nearly simultaneous dual NSI firings has been proven. The simulations indicate that the failure mode for simultaneous dual NSI firings may be caused by flow interactions between the flame channels. The shock waves from each initiator interact in the booster cavity resulting in a high pressure that prevents the gas and particulate velocity from rising in the booster cap region. This impedes the bulk of the particulate phase from impacting the booster cap and reduces the heat transfer to the booster cap since the particles do not impact it. Heat transfer calculations to the solid metal indicate that gas-phase convective heat transfer may not be adequate by itself and that energy transfer from the particulate phase may be crucial for the booster cap burn through.

Sachdev, Jai S.↗

I-BIEM calculations of the frequency dispersion and AC current distribution at disk and ring-disk electrodes

The Iterative Boundary Integral Equation Method (I-BIEM) has been applied to the problem of frequency dispersion at a disk electrode in a finite geometry. The I-BIEM permits the direct evaluation of the AC potential (a complex variable) using complex boundary conditions. The point spacing was made highly nonuniform, to give extremely high resolution in those regions where the variables change most rapidly, i.e., in the vicinity of the edge of the disk. Results are analyzed with respect to IR correction, equipotential surfaces, and reference electrode placement. The current distribution is also examined for a ring-disk configuration, with the ring and the disk at the same AC potential. It is shown that the apparent impedance of the disk is inductive at higher frequencies. The results are compared to analytic calculations from the literature, and usually agree to better than 0.001%.

Boris D Cahan↗

I-BIEM calculations of the frequency dispersion and ac current distribution at disk and ring-disk electrodes

The Iterative Boundary Integral Equation Method (I-BIEM) has been applied to the problem of frequency dispersion at a disk electrode in a finite geometry. The I-BIEM permits the direct evaluation of the AC potential (a complex variable) using complex boundary conditions. The point spacing was made highly nonuniform, to give extremely high resolution in those regions where the variables change most rapidly, i.e., in the vicinity of the edge of the disk. Results are analyzed with respect to IR correction, equipotential surfaces, and reference electrode placement. The current distribution is also examined for a ring-disk configuration, with the ring and the disk at the same AC potential. It is shown that the apparent impedance of the disk is inductive at higher frequencies. The results are compared to analytic calculations from the literature, and usually agree to better than 0.001 percent.

Cahan, Boris D.↗