Search NASA⌕ Search

SEARCH · Search NASA

Results for “numerical dispersion”

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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
NASA NTRS2000

Direct Numerical Simulation of Automobile Cavity Tones

The Navier Stokes equation is solved computationally by the Dispersion-Relation-Preserving (DRP) scheme for the flow and acoustic fields associated with a laminar boundary layer flow over an automobile door cavity. In this work, the flow Reynolds number is restricted to R(sub delta*) < 3400; the range of Reynolds number for which laminar flow may be maintained. This investigation focuses on two aspects of the problem, namely, the effect of boundary layer thickness on the cavity tone frequency and intensity and the effect of the size of the computation domain on the accuracy of the numerical simulation. It is found that the tone frequency decreases with an increase in boundary layer thickness. When the boundary layer is thicker than a certain critical value, depending on the flow speed, no tone is emitted by the cavity. Computationally, solutions of aeroacoustics problems are known to be sensitive to the size of the computation domain. Numerical experiments indicate that the use of a small domain could result in normal mode type acoustic oscillations in the entire computation domain leading to an increase in tone frequency and intensity. When the computation domain is expanded so that the boundaries are at least one wavelength away from the noise source, the computed tone frequency and intensity are found to be computation domain size independent.

Kurbatskii, Konstantin↗
DOE OSTI2025

Finite-volume quantization condition from the N/D representation

We propose a new model-independent method for determining hadronic resonances from lattice QCD. The formalism is derived from the general principles of unitarity and analyticity, as encoded in the N/D representation of a partial-wave two-body amplitude. The associated quantization condition relates the finite-volume spectrum to the infinite-volume numerator $\mathscr{N}$, used to reconstruct the scattering amplitude from dispersive relations. Unlike the original Lüscher condition, this new formalism is valid for energies coinciding with the left-hand cuts from arbitrary one- and multiparticle exchanges.

Dawid, Sebastian M. [Univ. of Washington, Seattle,↗
DOE OSTI2021

X-ray free electron laser heating of water and gold at high static pressure

Probing of reactive materials such as H 2 O ices and fluids at the high pressures and temperatures of planetary interiors is limited by unwanted chemical reactions and confinement failure. Faster experiments can mitigate such issues, but the common approach of adiabatic compression limits the conditions achieved. This study demonstrates a fast experimental strategy for the creation and probing of selected extreme states using static compression coupled with ultrafast X-ray laser heating. Indirect X-ray heating of H 2 O through the use of a gold absorber is evidenced by sample melting inferred from textural changes in the H 2 O diffraction lines and inter-dispersion of gold and H 2 O melts. Coupled with numerical analysis of femtosecond energy absorption, thermal equilibration, and heat transfer, all evidence indicates that temperatures in excess of an electron volt have been reached in the H 2 O at high pressure. Even after repeated heating, samples stayed chemically unchanged from the starting material.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Dispersion and the speed-limited particle-in-cell algorithm

This paper discusses temporally continuous and discrete forms of the speed-limited particle-in-cell (SLPIC) method first treated by Werner et al. [Phys. Plasmas 25, 123512 (2018)]. The dispersion relation for a 1D1V electrostatic plasma whose fast particles are speed-limited is derived and analyzed. By examining the normal modes of this dispersion relation, we show that the imposed speed-limiting substantially reduces the frequency of fast electron plasma oscillations while preserving the correct physics of lower-frequency plasma dynamics (e.g. ion acoustic wave dispersion and damping). We then demonstrate how the timestep constraints of conventional electrostatic particle-in-cell methods are relaxed by the speed-limiting approach, thus enabling larger timesteps and faster simulations. Here, these results indicate that the SLPIC method is a fast, accurate, and powerful technique for modeling plasmas wherein electron kinetic behavior is nontrivial (such that a fluid/Boltzmann representation for electrons is inadequate) but evolution is on ion timescales.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Latest Updates on JSPEC - An IBS and Electron Cooling Simulation Program

JSPEC (JLab Simulation Package on Electron Cooling) is an open-source C++ program developed at Jefferson Lab, which simulates the evolution of the ion beam under the influence of both IBS and electron cooling effects. In this paper, we will report the latest updates to JSPEC. Firstly, we have added theoretical and numerical models that simulate the effect of the electron beam dispersion on non-magnetized cooling. Secondly, the cooler can now be treated as an element with length, rather than a thin lens. This change will impact the modelling of the ions and the electrons in the cooling rate calculation for both magnetized and non-magnetized cooling. Numeri-cal results will be provided to demonstrate the perfor-mance of the new models.

Zhang, H.↗
NASA NTRS1975

Experimental observation of the ion-ion beam instability in a cylindrical geometry

The ion-ion beam instability is investigated in a cylindrical geometry. In the presence of a cylindrically symmetric ion beam a coherent ion-ion beam instability is observed for a well-defined range of beam velocity. This instability gives rise to cylindrical standing waves whose frequency depends upon beam velocity. The beam distribution function is measured to be a ring in velocity space. The dispersion relation appropriate for the cylindrical standing waves is solved numerically and is shown to be in good agreement with the experimentally observed instability.

Romesser, T.↗
NASA NTRS1984

Tearing instability in an anisotropic neutral sheet

A study is made of the collisionless tearing-mode stability properties of a field-reversed plasma layer whose temperature distribution is anisotropic. The plasma is confined by its self magnetic field with no external field. A kinetic description is used for both ions and electrons. The effects of the axis-crossing and nonaxis-crossing orbits are discussed. It is found that the conventional technique of matching the inner and outer asymptotic solutions at the electron inner-region is inadequate for the anisotropic case. An intermediate region in which the axis-crossing ion orbits are important is identified. The eigenvalue equation is solved using both analytic approximations and numerical methods to obtain the eigenmode structure and the linear dispersion relation. Previously announced in STAR as N84-14920

Chen, J.↗
NASA NTRS1985

Perturbation-iteration theory for analyzing microwave striplines

A perturbation-iteration technique is presented for determining the propagation constant and characteristic impedance of an unshielded microstrip transmission line. The method converges to the correct solution with a few iterations at each frequency and is equivalent to a full wave analysis. The perturbation-iteration method gives a direct solution for the propagation constant without having to find the roots of a transcendental dispersion equation. The theory is presented in detail along with numerical results for the effective dielectric constant and characteristic impedance for a wide range of substrate dielectric constants, stripline dimensions, and frequencies.

Kretch, B. E.↗
NASA NTRS1985

Broadband electrostatic noise produced by ion beams in the earth's magnetotail

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.

Omidi, N.↗
NASA NTRS1990

Recurrence in truncated Boussinesq models for nonlinear waves in shallow water

The rapid spatial recurrence of weakly nonlinear and weakly dispersive progressive shallow-water waves is examined using a numerical integration technique on the discretized and truncated form of the Boussinesq equations. This study primarily examines recurrence in wave fields with Ursell number O(1) and characterizes the sensitivity of recurrence to initial spectral shape and number of allowed frequency modes. It is shown that the rapid spatial recurrence is not an inherent property of the considered Boussinesq systems for evolution distances of 10-50 wavelengths. The main result of the study is that highly truncated Boussinesq models of resonant shallow-water ocean surface gravity waves predict rapid multiple recurrence cycles, but that this is an artifact dependent on the number of allowed modes. For initial conditions consisting of essentially all energy concentrated in a single mode, damping of the recurrence cycles increases as the number of low-power background modes increases. When more than 32 modes are allowed, the recurrence behavior is relatively insensitive to the number of allowed modes.

Elgar, Steve↗
NASA NTRS1991

A computer code for multiphase all-speed transient flows in complex geometries. MAST version 1.0

The operation of the MAST code, which computes transient solutions to the multiphase flow equations applicable to all-speed flows, is described. Two-phase flows are formulated based on the Eulerian-Lagrange scheme in which the continuous phase is described by the Navier-Stokes equation (or Reynolds equations for turbulent flows). Dispersed phase is formulated by a Lagrangian tracking scheme. The numerical solution algorithms utilized for fluid flows is a newly developed pressure-implicit algorithm based on the operator-splitting technique in generalized nonorthogonal coordinates. This operator split allows separate operation on each of the variable fields to handle pressure-velocity coupling. The obtained pressure correction equation has the hyperbolic nature and is effective for Mach numbers ranging from the incompressible limit to supersonic flow regimes. The present code adopts a nonstaggered grid arrangement; thus, the velocity components and other dependent variables are collocated at the same grid. A sequence of benchmark-quality problems, including incompressible, subsonic, transonic, supersonic, gas-droplet two-phase flows, as well as spray-combustion problems, were performed to demonstrate the robustness and accuracy of the present code.

Chen, C. P.↗
NASA NTRS1992

Parametric instabilities of large amplitude Alfven waves with obliquely propagating sidebands

This paper presents a brief report on properties of the parametric decay and modulational, filamentation, and magnetoacoustic instabilities of a large amplitude, circularly polarized Alfven wave. We allow the daughter and sideband waves to propagate at an arbitrary angle to the background magnetic field so that the electrostatic and electromagnetic characteristics of these waves are coupled. We investigate the dependance of these instabilities on dispersion, plasma/beta, pump wave amplitude, and propagation angle. Analytical and numerical results are compared with numerical simulations to investigate the full nonlinear evolution of these instabilities.

Vinas, A. F.↗
NASA NTRS1995

Far-ultravoilet and optical spectrscopy of the detached white dwarf-M binary RE J1629+780

We present optical and far-ultraviolet observations, including IUE high- and low-dispersion spectroscopy, of the (DA2+dM4-5) close binary, RE J1629+780, a very bright ROSAT EUV source. The optical data reveal emission reversals in the cores of the Balmer lines, which vary in strength but reveal no evidence of wavelength shifts reflecting orbital motion, suggesting that the red dwarf and the white dwarf are widely separated. Thus it is unlikely that irradiation plays a dominant role in the activity of the red dwarf. One dramatic increase in the intensity of the H-beta emission strength in the course of one hour indicates that the observed intensity variations are most likely related to flare activity on the red dwarf. Our analysis of Balmer and Lyman alpha features yield T(sub eff) = 42,500 +/- 1300 K and log g = 7.6 +/- 0.3 for the DA2 white dwarf, which is consistent with the overall 900 to 3200 A energy distribution derived from Voyager and IUE low-dispersion spectra. Our analysis of three IUE echelle spectra reveal numerous lines of interstellar origin. Using a total of 13 interstellar matter (ISM) features we find a mean radial velocity of -13.57 +/- 1.66 km/s for the interstellar features. However, we found no evidence of trace metal absorbers levitated or accreted in the white-dwarf photosphere nor any evidence of circumstellar or wind features. We estimate a lower limit hydrogen layer mass of 2.1 x 10(exp -15) solar mass for a stratified atmosphere based upon the nondetection of He II (1640). Tighter constraints on the helium abundance and the hydrogen layer mass, from EUV and X-ray data, are presented for both homogeneous and stratified atmospheric configurations.

Sion, Edward M.↗
NASA NTRS2007

Numerical Investigation of Plasma Flows in Magnetic Nozzles

Magnetic nozzles are used in many laboratory experiments in which plasma flows are to be confined, cooled. accelerated, or directed. At present, however, there is no generally accepted theoretical description that explains the phenomena of plasma detachment from an externally-imposed magnetic field. This is an important problem in the field of plasma propulsion, where the ionized gas must detach from the applied, solenoidal magnetic field to realize thrust production. In this paper we simulate a plasma flowing in the presence of an applied magnetic field using a multidimensional numerical simulation tool that includes theoretical models of the various dispersive and dissipative processes present in the plasma. This is an extension of the simulation tool employed in previous work by Sankaran et al. The new tool employs the same formulation of the governing equation set. but retains the axial and radial components of magnetic field and the azimuthal component of velocity that were neglected in other works. We aim to compare the computational results with the various proposed magnetic nozzle detachment theories to develop an understanding of the physical mechanisms that cause detachment. An applied magnetic field topology is obtained using a magnetostatic field solver and this field is superimposed on the time-dependent magnetic field induced in the plasma to provide a self-consistent field description. The applied magnetic field and model geometry match those found in experiments by Kurtki and Okada. We model this geometry because there ts a substantial amount of experimental data that can be compared to our computations, allowing for validation of the model. In addition, comparison of the simulation results with the experimentally obtained plasma parameters will provide insight into the mechanisms that lead to plasma detachment, revealing how the 3 scale with different input parameters.

Polzin, K. A.↗
NASA NTRS2008

Numerical Investigation of Plasma Detachment in Magnetic Nozzle Experiments

At present there exists no generally accepted theoretical model that provides a consistent physical explanation of plasma detachment from an externally-imposed magnetic nozzle. To make progress towards that end, simulation of plasma flow in the magnetic nozzle of an arcjet experiment is performed using a multidimensional numerical simulation tool that includes theoretical models of the various dispersive and dissipative processes present in the plasma. This is an extension of the simulation tool employed in previous work by Sankaran et al. The aim is to compare the computational results with various proposed magnetic nozzle detachment theories to develop an understanding of the physical mechanisms that cause detachment. An applied magnetic field topology is obtained using a magnetostatic field solver (see Fig. I), and this field is superimposed on the time-dependent magnetic field induced in the plasma to provide a self-consistent field description. The applied magnetic field and model geometry match those found in experiments by Kuriki and Okada. This geometry is modeled because there is a substantial amount of experimental data that can be compared to the computational results, allowing for validation of the model. In addition, comparison of the simulation results with the experimentally obtained plasma parameters will provide insight into the mechanisms that lead to plasma detachment, revealing how they scale with different input parameters. Further studies will focus on modeling literature experiments both for the purpose of additional code validation and to extract physical insight regarding the mechanisms driving detachment.

Sankaran, Kamesh↗
NASA NTRS2009

Numerical Investigation of Plasma Flows in Magnetic Nozzles

Magnetic nozzles are used in many laboratory experiments in which plasma flows are to be confined, cooled, accelerated, or directed. At present, however, there is no generally accepted theoretical description that explains the phenomena of plasma expansion in and detachment from an externally-imposed magnetic field. The latter is an especially important problem in the field of plasma propulsion, where the ionized gas must detach from the applied, solenoidal magnetic field to realize thrust production. In this paper we simulate a plasma flowing in the presence of an applied magnetic field using a multidimensional numerical simulation tool that includes theoretical models of the various dispersive and dissipative processes present in the plasma. This is an extension of the simulation tool employed in previous work by Sankaran et al. The new tool employs the same formulation of the governing equation set, but retains the axial and radial components of magnetic field and the azimuthal component of velocity that were neglected. We aim to compare the computational results with the various proposed magnetic nozzle detachment theories to develop an understanding of the physical mechanisms that cause detachment. An applied magnetic field topology is obtained using a magnetostatic field solver, and this field is superimposed on the time-dependent magnetic field induced in the plasma to provide a self-consistent field description. The applied magnetic field and model geometry match those found in experiments by Kuriki and Okada. 4 A schematic showing the setup used in those experiments is shown. We model this geometry because there is a substantial amount of experimental data that can be compared to our computations, allowing for validation of the model. In addition, comparison of the simulation results with the experimentally obtained plasma parameters will provide insight into the mechanisms that lead to plasma detachment, revealing how they scale with different input parameters.

Sankaran, Kamesh↗
NASA NTRS1991

Hydrodynamic behavior and electrochemical impedance of the Hanging Meniscus Rotating Disk (HMRD) electrode. I - Meniscus shape under rotation. II - I-BIEM calculations of frequency dispersion and minimization

The shape equations for an HMRD in static and rotating configurations are developed and solved numerically. A rationale for the applicability of the standard Levich equations to the rotating case is given. The region of stability of the HMRD is examined, and the observed small negative intercept for a Levich plot is explained. The iterative boundary integral equation method is applied to the problem of frequency dispersion at an HMRD electrode. It is shown that a range of disk sizes and heights can be chosen to give almost uniform primary and secondary current distribution and minimal frequency disperison.

Cahan, Boris D.↗
DOE OSTI2025

Ten-moment fluid modeling of the Weibel instability

We investigate the one-dimensional non-relativistic Weibel instability through the capture of anisotropic pressure tensor dynamics using an implicit 10-moment fluid model that employs the electromagnetic Darwin approximation. The results obtained from the 10-moment model are compared with an implicit particle-in-cell simulation. The linear growth rates obtained from the numerical simulations are in good agreement with the theoretical fluid and kinetic dispersion relations. The fluid dispersion relations are derived using Maxwell’s equations and the Darwin approximation. We also show that the magnetohydrodynamic approximation can be used to model the Weibel instability if one accounts for an anisotropic pressure tensor and unsteady terms in the generalised Ohm’s law. In addition, we develop a preliminary theory for the saturation magnetic field strength of the Weibel instability, showing good agreement with the numerical results.

Kuldinow, D. A. (ORCID:0000000319730196)↗