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 181 records · Page 10

Decoupling the Role of Inertia and Gravity on Particle Dispersion

Turbulent gas flows laden with small, dense particles are encountered in a wide number of important applications in both industrial settings and aerodynamics applications. Particle interactions with the underlying turbulent flow are exceedingly complex and, consequently, difficult to accurately model. The difficulty arises primarily due to the fact that response of a particle to the local environment is dictated by turbulence properties in the reference frame moving with the particle (particle-Lagrangian). The particle-Lagrangian reference frame is in turn dependent upon the particle relaxation time (time constant) as well as gravitational drift. The combination of inertial and gravitational effects in this frame complicates our ability to accurately predict particle-laden flows since measurements in the particle-Lagrangian reference frame are difficult to obtain. Therefore, in this work we will examine separately the effects of inertia and gravitational drift on particle dispersion through a combination of physical and numerical experiments. In this study, particle-Lagrangian measurements will be obtained in physical experiments using stereo image velocimetry. Gravitational drift will be varied in the variable-g environments of the NASA DC-9 and in the zero-g environment at the drop tower at NASA-Lewis. Direct numerical simulations will be used to corroborate the measurements from the variable-g experiments. We expect that this work will generate new insight into the underlying physics of particle dispersion and will, in turn, lead to more accurate models of particle transport in turbulent flows.

Rogers, Chris↗

Tipping points emerge from weak mutualism in metacommunities

The coexistence of obligate mutualists is often precariously close to tipping points where small environmental changes can drive catastrophic shifts in species composition. For example, microbial ecosystems can collapse by the decline of a strain that provides an essential resource on which other strains cross-feed. Here, we show that tipping points, ecosystem collapse, bistability and hysteresis arise even with very weak (non-obligate) mutualism provided the population is spatially structured. Based on numeric solutions of a metacommunity model and mean-field analyses, we demonstrate that weak mutualism lowers the minimal dispersal rate necessary to avoid stochastic extinction, while species need to overcome a mean threshold density to survive in this low dispersal rate regime. Our results allow us to make numerous predictions for mutualistic metacommunities regarding tipping points, hysteresis effects, and recovery from external perturbations, and let us draw general conclusions for ecosystems even with random, not necessarily mutualistic, interactions and systems with density-dependent dispersal rather than direct mutualistic interactions.

59 BASIC BIOLOGICAL SCIENCES↗

A scalable multidimensional fully implicit solver for Hall magnetohydrodynamics

We propose an optimally performant fully implicit algorithm for the Hall magnetohydrodynamics (HMHD) equations based on multigrid-preconditioned Jacobian-free Newton-Krylov methods. HMHD is a challenging system to solve numerically because it supports stiff fast dispersive waves. The preconditioner is formulated using an operator-split approximate block factorization (Schur complement), informed by physics insight. We use a vector-potential formulation (instead of a magnetic field one) to allow a clean segregation of the problematic $\nabla$ x $\nabla$ x operator in the electron Ohm's law subsystem. This segregation allows the formulation of an effective damped block-Jacobi smoother for multigrid. We demonstrate by analysis that our proposed block-Jacobi iteration is convergent and has the smoothing property. The resulting HMHD solver is verified linearly with wave propagation examples, and nonlinearly with the GEM challenge reconnection problem by comparison against another HMHD code. We demonstrate the excellent algorithmic and parallel performance of the algorithm up to 16384 MPI tasks in two dimensions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Structure and overstability of resistive modes with runaway electrons

In this paper, we investigate the effects of runaway electron current on the dispersion relation of resistive magnetohydrodynamic modes in tokamaks. We present a new theoretical model to derive the dispersion relation, which is based on the asymptotic analysis of the resistive layer structure of the modes. It is found that in addition to the conventional resistive layer, a new runaway current layer can emerge whose properties depend on the ratio of the Alfvén velocity to the runaway electron convection speed. Due to the contribution from this layer, both the tearing mode and kink mode will have a real frequency in addition to a growth rate. The derived dispersion relation has been compared with numerical results using both a simplified eigenvalue calculation and a M3D-C 1 linear simulation, and good agreement is found in both cases.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Verification of a fully implicit particle-in-cell method for the <!--${MathJax: TeX-AMS-MML_HTMLorMML}--> v &#x2225; -formalism of electromagnetic gyrokinetics in the XGC code

A fully implicit particle-in-cell method for handling the v ∥ -formalism of electromagnetic gyrokinetics has been implemented in XGC. By choosing the v ∥ -formalism, here we avoid introducing the nonphysical skin terms in Ampère's law, which are responsible for the well-known “cancellation problem” in the p ∥ -formalism. The v ∥ -formalism, however, is known to suffer from a numerical instability when explicit time integration schemes are used due to the appearance of a time derivative in the particle equations of motion from the inductive component of the electric field. Here, using the conventional δf scheme, we demonstrate that our implicitly discretized algorithm can provide numerically stable simulation results with accurate dispersive properties. We verify the algorithm using a test case for shear Alfvén wave propagation in addition to a case demonstrating the ion temperature gradient-kinetic ballooning mode (ITG-KBM) transition. The ITG-KBM transition case is compared to results obtained from other δf gyrokinetic codes/schemes, whose verification has already been archived in the literature.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Origins of electronic bands in the antiferromagnetic topological insulator MnBi 2 Te 4

Despite the rapid progress in understanding the first intrinsic magnetic topological insulator MnBi 2 Te 4 , its electronic structure remains a topic under debates. In this work, we perform a thorough spectroscopic investigation into the electronic structure of MnBi 2 Te 4 via laser-based angle-resolved photoemission spectroscopy. Through quantitative analysis, we estimate an upper bound of 3 meV for the gap size of the topological surface state. Furthermore, our circular dichroism measurements reveal band chiralities for both the topological surface state and quasi-2D bands, which can be well reproduced in a band hybridization model. A numerical simulation of energy-momentum dispersions based on a four-band model with an additional step potential near the surface provides a promising explanation for the origin of the quasi-2D bands. Our study represents a solid step forward in reconciling the existing controversies in the electronic structure of MnBi 2 Te 4 , and provides an important framework to understand the electronic structures of other relevant topological materials MnBi 2n Te 3n+1 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Modulation of energetic particle fluxes by a mixed mode of transverse and compressional waves

Modulation characteristics of particle fluxes in the presence of a mixed mode of compressional and transverse magnetic waves at hydromagnetic frequencies are investigated through kinetic perturbation of the distribution function. The magnetospheric medium where the particles are modulated contains both the magnetic and pressure gradients. The modulation features are found to be strongly dependent on the energy and pitch angle of the particles. Drifting particles can resonate with waves whose phase velocities are close to their drift velocities. When this occurs, the modulation amplitudes become significantly large and large phase shifts will occur. It is pointed out that resonance is important for particles with mid pitch angles (40-70 deg). The phase shift between the particle modulations and the magnetic field oscillations are strongly controlled by the combined effects of transverse and compressional wave components and/or the occurrence of drift resonance. Numerical calculations are performed using the dispersion relation of drift mirror Alfven waves as an example of waves with both compressional and transverse components.

Lin, C. S.↗

A search for spectroscopic binaries among the runaway O type stars

Numerous radial velocity measurements of medium dispersion were made for the 10 brighter stars given in Stone's list of very probable O type runaways. All plates were measured with the KPNO PDS microdensitometer, and a new iterative reductional analysis was used to derive plate velocities, which are estimated to be 1.6 times more accurate internally than those found by using the traditional method. Of thse stars, psi Per, alpha Cam, HD 188209, and 26 Cep are identified as probable velocity variables, while 9 Sge, lambda Cep, and HD 218915 are classed as possible variables. If the source of this variability is Keplerian rather than atmospheric, which cannot be established unequivocally from the observations of this paper, psi Per could be a spectroscopic binary with a black hole companion, and at least 1.2 solar mass. The detection of runaway binary systems from radial velocity measurements is discussed.

Stone, R. C.↗

Electromagnetic ion beam instabilities - Hot beams at interplanetary shocks

This paper considers Vlasov instabilities driven by a very hot ion beam streaming parallel to a magnetic field. The linear theory of electromagnetic instabilities driven by such a beam in a homogeneous, collisionless, nonrelativistic plasma is used. Numerical solutions of the full linear dispersion equation for bi-Maxwellian distribution functions are presented. If the thermal speed of the beam is greater than its drift speed, the two dominant modes are the right-hand and left-hand resonant ion beam instabilities. The parametric dependencies of the threshold drift speeds and growth rates for both modes are presented. In particular, the thresholds are shown to be sensitive functions of the beam temperature anisotropies. Arguments are presented that these two modes are driven by the suprathermal ion component at interplanetary shocks, and that the growth of these modes is sufficient to account for the amplitudes of MHD-like waves observed at such shocks.

Gary, S. P.↗

Electromagnetic electron temperature anisotropy instabilities

This paper considers electromagnetic Vlasov instabilities driven by electron temperature anisotropies in a homogeneous, nonrelativistic magnetized plasma. Numerical solutions of the full linear dispersion equation for bi-Maxwellian distribution functions and instabilities propagating parallel to the magnetic field are presented. Parametric dependences of the maximum growth rates of the electron fire hose and whistler anisotropy instabilities are given.

Gary, S. P.↗

The ion-ion acoustic instability

This paper considers the linear theory of ion-acoustic-like instabilities in a homogeneous Vlasov plasma with two ion components, a less dense beam and a more dense core, with a relative drift velocity. Numerical solutions of the full electrostatic dispersion equation are presented, and the properties of the ion-ion acoustic instability are studied in detail. At relatively cold beam temperatures, the instability is fluid-like, but it becomes a beam resonant kinetic instability as the beam temperature becomes of the order of the core temperature; if the mode is unstable, its threshold lies well below the threshold of the electron-ion acoustic instability. An electron temperature anisotropy enhances the instability, and at sufficiently large beam-core relative drift speeds, electron magnetization can either detract from or enhance the instability.

Gary, S. Peter↗

Magnetohydrodynamic waves and particle acceleration in a coronal loop

A model is presented for the acceleration of electrons in a flaring solar coronal loop. It is suggested that the primary energy release event in the flare establishes a spectrum of MHD waves within the loop which accelerate electrons by wave-particle resonant interactions. The nature of the MHD oscillations are described and numerical and analytic solution of the dispersion relation are examined. It is found that the behavior of a particle depends on whether or not it is in resonance with the wave spectrum. Resonant particles follow chaotic trajectories and can diffuse in velocity space between wave effective potentials, thereby gaining energy.

De La Beaujardiere, Jean-Francois↗

Theory of circular dielectric waveguide with anisotropic sheet cover

A circular dielectric waveguide consisting of an isotropic core covered with a thin anisotropic sheet is considered. The sheet is represented as a jump immittance and Maxwell's equations are applied. Solution of the boundary value problem yields the characteristic equation, or dispersion relation, which is then solved numerically. The results are verified for the step-index fiber and circular, metallic waveguides. Finally, examples are included to investigate the effects of the anisotropic sheet.

Kish, Jerome Thomas↗

Lagrangian and Eulerian statistics obtained from direct numerical simulations of homogeneous turbulence

Direct numerical simulation is used to study dispersion in decaying isotropic turbulence and homogeneous shear flow. Both Lagrangian and Eulerian data are presented allowing direct comparison, but at fairly low Reynolds number. The quantities presented include properties of the dispersion tensor, isoprobability contours of particle displacement, Lagrangian and Eulerian velocity autocorrelations and time scale ratios, and the eddy diffusivity tensor. The Lagrangian time microscale is found to be consistently larger than the Eulerian microscale, presumably due to the advection of the small scales by the large scales in the Eulerian reference frame.

Squires, Kyle D.↗

Dispersion Caused by the Penetration Effect in X-ray Compressors

Chirped X-ray pulse compression is a promising approach for generating ultra-short X-ray free electron laser (XFEL) pulses. The design of X-ray pulse compressors requires the careful control of group delay dispersion (GDD), which plays a critical role in achieving optimal compression. However, the penetration dispersion of crystals and multilayers can induce an extra GDD, which may result in over-compression or under-compression. In this study, we investigate the penetration dispersion of crystals and multilayers theoretically and numerically. Our results indicate that the extra GDD induced by the penetration effect increases as the bandwidth of the rocking curve decreases. Moreover, the extra GDD is nonlinear and can be mitigated by optimizing the configuration of X-ray pulse compressors. This work provides insights into the dispersion compensation and configuration optimization of X-ray pulse compressors, which are essential for generating ultra-short XFEL pulses.

47 OTHER INSTRUMENTATION↗

Particle dispersion in confined turbulent swirling flows

This paper reports a numerical investigation of confined swirling flows of gas-particle mixtures. A recently developed two-fluid, multiple-scale mixing model is applied to study the influence of particles on the intensity of the reverse flows of the gas phase and the effects of swirl on the particle dispersion in an annular expanding chamber under isothermal condition. The calculations were made for different swirl strength of injection of the annular jet into the mixing chamber. Results agree well qualitatively with experimental information available. It is also found that the calculated flow fields depend heavily on the prescription of the inlet flow conditions.

Chen, C. P.↗

Assessment of Total Variation Diminishing schemes in combustion instability problems

Numerical simulations of acoustic waves in a shear layer and in an idealized combustion chamber using high resolution Total Variation Diminishing (TVD) schemes has been carried out to study the effects of inherent scheme dissipation and dispersion errors of this class of problems. The numerical results are compared against available exact solutions to quantify these errors. Several popular TVD limiters widely used in the Computational Fluid Dynamics (CFD) community have been assessed. Osher-Chakravarthy limiters are modified so that they can be used in explicit schemes. Among all the limiters investigated, Osher-Chakravarthy third-order limiter is identified as having performed the best. It is also found that all TVD schemes have exceptionally small dispersive errors.

Wang, Z. J.↗