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

Comparison of Non-Parabolic Hydrodynamic Simulations for Semiconductor Devices

Parabolic drift-diffusion simulators are common engineering level design tools for semiconductor devices. Hydrodynamic simulators, based on the parabolic band approximation, are becoming more prevalent as device dimensions shrink and energy transport effects begin to dominate device characteristic. However, band structure effects present in state-of-the-art devices necessitate relaxing the parabolic band approximation. This paper presents simulations of ballistic diodes, a benchmark device, of Si and GaAs using two different non-parabolic hydrodynamic formulations. The first formulation uses the Kane dispersion relationship in the derivation of the conservation equations. The second model uses a power law dispersion relation {(hk)(exp 2)/2m = xW(exp Y)}. Current-voltage relations show that for the ballistic diodes considered. the non-parabolic formulations predict less current than the parabolic case. Explanations of this will be provided by examination of velocity and energy profiles. At low bias, the simulations based on the Kane formulation predict greater current flow than the power law formulation. As the bias is increased this trend changes and the power law predicts greater current than the Kane formulation. It will be shown that the non-parabolicity and energy range of the hydrodynamic model based on the Kane dispersion relation are limited due to the binomial approximation which was utilized in the derivation.

Smith, A. W.↗

Absorbing boundary conditions for second-order hyperbolic equations

A uniform approach to construct absorbing artificial boundary conditions for second-order linear hyperbolic equations is proposed. The nonlocal boundary condition is given by a pseudodifferential operator that annihilates travelling waves. It is obtained through the dispersion relation of the differential equation by requiring that the initial-boundary value problem admits the wave solutions travelling in one direction only. Local approximation of this global boundary condition yields an nth-order differential operator. It is shown that the best approximations must be in the canonical forms which can be factorized into first-order operators. These boundary conditions are perfectly absorbing for wave packets propagating at certain group velocities. A hierarchy of absorbing boundary conditions is derived for transonic small perturbation equations of unsteady flows. These examples illustrate that the absorbing boundary conditions are easy to derive, and the effectiveness is demonstrated by the numerical experiments.

Jiang, Hong↗

Absorbing boundary conditions for second-order hyperbolic equations

A uniform approach to construct absorbing artificial boundary conditions for second-order linear hyperbolic equations is proposed. The nonlocal boundary condition is given by a pseudodifferential operator that annihilates travelling waves. It is obtained through the dispersion relation of the differential equation by requiring that the initial-boundary value problem admits the wave solutions travelling in one direction only. Local approximation of this global boundary condition yields an nth-order differential operator. It is shown that the best approximations must be in the canonical forms which can be factorized into first-order operators. These boundary conditions are perfectly absorbing for wave packets propagating at certain group velocities. A hierarchy of absorbing boundary conditions is derived for transonic small perturbation equations of unsteady flows. These examples illustrate that the absorbing boundary conditions are easy to derive, and the effectiveness is demonstrated by the numerical experiments.

Jiang, Hong↗

Dispersion discontinuities of strong collisionless shocks

Linear fluid equations are used to estimate wave train properties of strong collisionless shocks. Fast shocks exhibit several dispersion changes with increasing Mach number. For perpendicular propagation into a finite-beta plasma, an ion cyclotron radius trailing wave train exists only for (M sub F)2 is smaller than 2. Oblique fast shocks have a leading ion inertia wave train if M sub A is smaller than root of M(+)/M(-) cos theta/2 and a trailing electron inertia train if M sub A is greater than root of M(+)/M(-) cos theta/2. If the downstream sound speed exceeds the flow speed, linear wave theory predicts a trailing ion acoustic structure which probably resides within the magnetic shock. For a turbulent shock model in which an effective electron-ion collision frequency exceeds the lower hybrid frequency, ions decouple from the magnetic field; the shock wave train now trails with electron inertia and electron gyroradius lengths. Comparisons of this turbulent model and observations on the earth's bow shock are made.

Coroniti, F. V.↗

Detecting Moving Targets by Use of Soliton Resonances

A proposed method of detecting moving targets in scenes that include cluttered or noisy backgrounds is based on a soliton-resonance mathematical model. The model is derived from asymptotic solutions of the cubic Schroedinger equation for a one-dimensional system excited by a position-and-time-dependent externally applied potential. The cubic Schroedinger equation has general significance for time-dependent dispersive waves. It has been used to approximate several phenomena in classical as well as quantum physics, including modulated beams in nonlinear optics, and superfluids (in particular, Bose-Einstein condensates). In the proposed method, one would take advantage of resonant interactions between (1) a soliton excited by the position-and-time-dependent potential associated with a moving target and (2) eigen-solitons, which represent dispersive waves and are solutions of the cubic Schroedinger equation for a time-independent potential.

Zak, Michael↗

Role of Compressive Viscosity and Thermal Conductivity on the Damping of Slow Waves in Coronal Loops with and Without Heating–Cooling Imbalance

In the present article, we derive a new dispersion relation for slow magnetoacoustic waves invoking the effect of thermal conductivity, compressive viscosity, radiation, and an unknown heating term along with the consideration of heating–cooling imbalance from linearized MHD equations. We solve the general dispersion relation to understand the role of compressive viscosity and thermal conductivity in the damping of slow waves in coronal loops with and without heating–cooling imbalance. We have analyzed the wave damping for the range of loop length L=50--500 Mm, temperature T=5--30 MK, and density ρ=10^(−11)--10^(−9) kg/cu.m. It was found that the inclusion of compressive viscosity along with thermal conductivity significantly enhances the damping of the fundamental mode oscillations in shorter (e.g. L=50 Mm) and super-hot (T>10 MK) loops. However, the role of viscosity in the damping is insignificant in longer (e.g. L=500 Mm) and hot loops (T≤10 MK) where, instead, thermal conductivity along with the presence of heating–cooling imbalance plays a dominant role. For shorter loops at a super-hot regime of temperature, the increment in the loop density substantially enhances the damping of the fundamental modes due to thermal conductivity when viscosity is absent, however, when the compressive viscosity is added the increase in density substantially weakens the damping. Thermal conductivity alone is found to play a dominant role in longer loops at lower temperatures (T≤10 MK), while compressive viscosity dominates the damping at super-hot temperatures (T>10 MK) in shorter loops. The predicted scaling law between damping time (τ) and wave period (P) is found to better match the observed SUMER (Solar Ultraviolet Measurements of Emitted Radiation) oscillations when the heating–cooling imbalance is taken into account in addition to thermal conductivity and compressive viscosity for the damping of the fundamental slow mode oscillations.

Abhinav Prasad↗

Fourier analysis of numerical algorithms for the Maxwell equations

The Fourier method is used to analyze the dispersive, dissipative, and isotropy errors of various spatial and time discretizations applied to the Maxwell equations on multi-dimensional grids. Both Cartesian grids and non-Cartesian grids based on hexagons and tetradecahedra are studied and compared. The numerical errors are quantitatively determined in terms of phase speed, wave number, propagation direction, gridspacings, and CFL number. The study shows that centered schemes are more efficient than upwind schemes. The non-Cartesian grids yield superior isotropy and higher accuracy than the Cartesian ones. For the centered schemes, the staggered grids produce less errors than the unstaggered ones. A new unstaggered scheme which has all the best properties is introduced. The study also demonstrates that a proper choice of time discretization can reduce the overall numerical errors due to the spatial discretization.

Liu, Yen↗

Invariant Theory for Dispersed Transverse Isotropy: An Efficient Means for Modeling Fiber Splay

Most soft tissues possess an oriented architecture of collagen fiber bundles, conferring both anisotropy and nonlinearity to their elastic behavior. Transverse isotropy has often been assumed for a subset of these tissues that have a single macroscopically-identifiable preferred fiber direction. Micro-structural studies, however, suggest that, in some tissues, collagen fibers are approximately normally distributed about a mean preferred fiber direction. Structural constitutive equations that account for this dispersion of fibers have been shown to capture the mechanical complexity of these tissues quite well. Such descriptions, however, are computationally cumbersome for two-dimensional (2D) fiber distributions, let alone for fully three-dimensional (3D) fiber populations. In this paper, we develop a new constitutive law for such tissues, based on a novel invariant theory for dispersed transverse isotropy. The invariant theory is based on a novel closed-form splay invariant that can easily handle 3D fiber populations, and that only requires a single parameter in the 2D case. The model is polyconvex and fits biaxial data for aortic valve tissue as accurately as the standard structural model. Modification of the fiber stress-strain law requires no re-formulation of the constitutive tangent matrix, making the model flexible for different types of soft tissues. Most importantly, the model is computationally expedient in a finite-element analysis.

Freed, alan D.↗

Long nonlinear waves in stratified shear flows

The propagation of finite-amplitude internal waves in a shear flow is considered for wavelengths that are long compared to the shear-layer thickness. Both singular and regular modes are investigated, and the equation governing the amplitude evolution is derived. The theory is generalized to allow for a radiation condition when the region outside the stratified shear layer is unbounded and weakly stratified. In this case, the evolution equation contains a damping term describing energy loss by radiation which can be used to estimate the persistence of solitary waves or nonlinear wave packets in realistic environments. A continuous three-layer model is studied in detail and closed-form expressions are obtained for the phase speed and the coefficients of the nonlinear and dispersive terms in the amplitude equation as a function of Richardson number.

Maslowe, S. A.↗

Plasma Dispersion Function for the Kappa Distribution

The plasma dispersion function is computed for a homogeneous isotropic plasma in which the particle velocities are distributed according to a Kappa distribution. An ordinary differential equation is derived for the plasma dispersion function and it is shown that the solution can be written in terms of Gauss' hypergeometric function. Using the extensive theory of the hypergeometric function, various mathematical properties of the plasma dispersion function are derived including symmetry relations, series expansions, integral representations, and closed form expressions for integer and half-integer values of K.

Podesta, John J.↗

Solitary-wave interaction

The interaction of solitary-wave solutions of a model equation for long waves in dispersive media is examined numerically. It is found that the waves do not emerge from the interaction unscathed. Instead, two new solitary waves, having slightly different amplitudes from the original waves, together with a small dispersive tail are generated as a result of the interaction.

Bona, J. L.↗

Theory of high-harmonic rectangular gyrotron for TE(mn) modes

In this paper, the linearized relativistic Vlasov equation is solved and a dispersion relation is calculated for the interaction between a relativistic electron beam and the electromagnetic fields of a rectangular waveguide supporting TE(mn) modes. The dispersion relation is simplified to the special case of a frame of reference moving with the electrons and the resulting coupling coefficient epsilon(mn)super l is calculated for various TE(mn) modes. The dependence of epsilon(mn)super l on the parity of the harmonic number l and various mn modes of the waveguide are discussed.

Ferendeci, A. M.↗

Wave modes of the Io plasma torus

The corotating magnetospheric convection model is analyzed with the objective of deriving the fundamental wave modes present in the Io torus. The time-dependent equations of motion are linearized for the case of the absence of plasma mass and momentum sources. By separating the magnetosphere-ionosphere coupling equation into two equations and by adopting a specific analytic form for the equilibrium plasma density, a cubic equation is derived for the dispersion relation defining the angular frequencies of the fundamental wave modes. A detailed analysis of the dispersion relation yields three fundamental wave modes. One of the modes is identified as the previously studied inertial interchange (unstable) mode, while the other two are newly derived decaying modes. The results obtained are discussed in the context of plasma transport in the Io torus.

Summers, D.↗

Nonlinear modes of the tensor Dirac equation and CPT violation

Recently, it has been shown that Dirac's bispinor equation can be expressed, in an equivalent tensor form, as a constrained Yang-Mills equation in the limit of an infinitely large coupling constant. It was also shown that the free tensor Dirac equation is a completely integrable Hamiltonian system with Lie algebra type Poisson brackets, from which Fermi quantization can be derived directly without using bispinors. The Yang-Mills equation for a finite coupling constant is investigated. It is shown that the nonlinear Yang-Mills equation has exact plane wave solutions in one-to-one correspondence with the plane wave solutions of Dirac's bispinor equation. The theory of nonlinear dispersive waves is applied to establish the existence of wave packets. The CPT violation of these nonlinear wave packets, which could lead to new observable effects consistent with current experimental bounds, is investigated.

Reifler, Frank J.↗

The possible role of high-frequency waves in heating solar coronal loops

We investigate the role of high-frequency waves in the heating of solar active region coronal loops. We assume a uniform background magnetic field, and we introduce a density stratification in a direction perpendicular to this field. We focus on ion compressive viscosity as the damping mechanism of the waves. We incorporate viscosity self-consistently into the equations, and we derive a dispersion relation by adopting a slab model, where the density inside the slab is greater than that outside. Such a configuration supports two types of modes: surface waves and trapped body waves. In order to determine under what conditions these waves may contribute to the heating of active regions, we solve our dispersion relation for a range of densities, temperatures, magnetic field strengths, density ratios, wavevector magnitudes, wavevector ratios, and slab widths. We find that surface waves exhibit very small damping, but body waves can potentially damp at rates needed to balance radiative losses. However, the required frequencies of these body waves are very high. For example, the wave frequency must be at least 5.0/s for a slab density of 10(exp 9,5)/cc, a slab temperature of 10(exp 6,5) K, a field strength of 100 G, and a density ratio of 5. For a slab density of 10(exp 10)/cc, this frequency increases to 8.8/s. Although these frequencies are very high, there in no observational evidence to rule out their existence, and they may be generated both below the corona and at magnetic reconnection sites in the corona. However, we do find that, for slab densities of 10(exp 10)/cc or less, the dissipation of high-frequency waves will be insufficient to balance the radiative losses if the magnetic field strength exceeds roughly 200 G. Because the magnetic field is known to exceed 200 G in many active region loops, particularly low-lying loops and loops emanating from sunspots, it is unlikely that high-frequency waves can provide sufficient heating in these regions.

Porter, Lisa J.↗

Development of Numerical Tools for the Investigation of Plasma Detachment from Magnetic Nozzles

A multidimensional numerical simulation framework aimed at investigating the process of plasma detachment from a magnetic nozzle is introduced. An existing numerical code based on a magnetohydrodynamic formulation of the plasma flow equations that accounts for various dispersive and dissipative processes in plasmas was significantly enhanced to allow for the modeling of axisymmetric domains containing three.dimensiunai momentum and magnetic flux vectors. A separate magnetostatic solver was used to simulate the applied magnetic field topologies found in various nozzle experiments. Numerical results from a magnetic diffusion test problem in which all three components of the magnetic field were present exhibit excellent quantitative agreement with the analytical solution, and the lack of numerical instabilities due to fluctuations in the value of del(raised dot)B indicate that the conservative MHD framework with dissipative effects is well-suited for multi-dimensional analysis of magnetic nozzles. Further studies will focus on modeling literature experiments both for the purpose of code validation and to extract physical insight regarding the mechanisms driving detachment.

Sankaran, Kamesh↗

Analysis of Kelvin Helmholtz Instabilities of Plasma Jets

Ulysses data indicates density fluctuations which are 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.↗