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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
NASA NTRS1991

Motion of particles with inertia in a compressible free shear layer

The effects of the inertia of a particle on its flow-tracking accuracy and particle dispersion are studied using direct numerical simulations of 2D compressible free shear layers in convective Mach number (Mc) range of 0.2 to 0.6. The results show that particle response is well characterized by tau, the ratio of particle response time to the flow time scales (Stokes' number). The slip between particle and fluid imposes a fundamental limit on the accuracy of optical measurements such as LDV and PIV. The error is found to grow like tau up to tau = 1 and taper off at higher tau. For tau = 0.2 the error is about 2 percent. In the flow visualizations based on Mie scattering, particles with tau more than 0.05 are found to grossly misrepresent the flow features. These errors are quantified by calculating the dispersion of particles relative to the fluid. Overall, the effect of compressibility does not seem to be significant on the motion of particles in the range of Mc considered here.

Samimy, M.↗
DOE OSTI2023

Three-dimensional coupling of electron cyclotron drift instability and ion–ion two stream instability

Electron cyclotron drift instability (ECDI) and ion–ion two stream instability (IITSI) are both kinetic instabilities that can be present in low-temperature, partially magnetized plasmas. The coupling of instabilities in a three-dimensional configuration leads to the existence of more than one unstable roots to the kinetic dispersion relation. In this paper, a generalized method has been developed for numerically evaluating solutions to the three-dimensional dispersion relation for coupled ECDI and IITSI, assuming cold singly and doubly charged ions and a Maxwellian velocity distribution function for the electrons. The present study demonstrates the coupling between ECDI and IITSI that affects the most unstable mode as a function of the wavenumbers in three dimensions and various plasma properties, including the applied electric field, magnetic field, electron temperature, ion velocities, and plasma density. One of the most notable results is that, while the most unstable mode with the largest growth rate is in the direction of the E×B drift in the two-dimensional cases, the most unstable mode for the three-dimensional configuration occurs in the oblique direction between the applied electric field and the E×B drift. Here, this agrees with experimental observations in cross field plasma sources using coherent Thomson scattering.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1995

Wave propagation and scattering in computational aeroacoustics

This paper describes predictions of model problems in computational aeroacoustics. Two problem classes are considered: The first (Workshop Category 3, Problem 1 and 2) considers two-dimensional wave propagation and non-reflecting boundary conditions in the presence of a mean flow. The second (Workshop Category 4, Problem 1) examines wall boundary conditions. For the last problem we introduce the Impedance Mismatch Method (IMM) to treat the solid wall boundaries. In this method the solid wall is simulated using a wall region in which the characteristic impedance is set to a different value to that in the fluid region. This method has advantages over traditional solid wall boundary conditions including simplicity of coding, speed of computations, and the ability to treat curved boundaries efficiently. Several numerical examples are given in addition to the Workshop Problems. The discretization of the Euler equations is performed in all cases with a Dispersion-Relation-Preserving (DRP) algorithm. The numerical results are compared with either analytical solutions or solutions obtained using traditional solid wall boundary conditions.

Chung, Cathy↗
DOE OSTI2023

Perivascular network segmentations derived from high-field MRI and their implications for perivascular and parenchymal mass transport in the rat brain

A custom segmentation workflow was applied to ex vivo high-field MR images of rat brains acquired following in vivo intraventricular contrast agent infusion to generate maps of the perivascular spaces (PVS). The resulting perivascular network segmentations enabled analysis of perivascular connections to the ventricles, parenchymal solute clearance, and dispersive solute transport within PVS. Numerous perivascular connections between the brain surface and the ventricles suggest the ventricles integrate into a PVS-mediated clearance system and raise the possibility of cerebrospinal fluid (CSF) return from the subarachnoid space to the ventricles via PVS. Assuming rapid solute exchange between the PVS and CSF spaces primarily by advection, the extensive perivascular network decreased the mean clearance distance from parenchyma to the nearest CSF compartment resulting in an over 21-fold reduction in the estimated diffusive clearance time scale, irrespective of solute diffusivity. This corresponds to an estimated diffusive clearance time scale under 10 min for amyloid-beta which suggests that the widespread distribution of PVS may render diffusion an effective parenchymal clearance mechanism. Additional analysis of oscillatory solute dispersion within PVS indicates that advection rather than dispersion is likely the primary transport mechanism for dissolved compounds greater than 66 kDa in the long (> 2 mm) perivascular segments identified here, although dispersion may be significant for smaller compounds in shorter perivascular segments.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS1991

Design of algorithms for a dispersive hyperbolic problem

In order to develop numerical schemes for stiff problems, a model of relaxing heat flow is studied. To isolate those errors unavoidably associated with discretization, a method of characteristics is developed, containing three free parameters depending on the stiffness ratio. It is shown that such 'decoupled' schemes do not take into account the interaction between the wave families, and hence result in incorrect wavespeeds. Schemes can differ by up to two orders of magnitude in their rms errors, even while maintaining second-order accuracy. 'Coupled' schemes which account for the interactions are developed to obtain two additional free parameters. Numerical results are given for several decoupled and coupled schemes.

Roe, Philip L.↗
DOE OSTI2022

Well-Balanced Second-Order Convex Limiting Technique for Solving the Serre–Green–Naghdi Equations

In this article, we introduce a numerical method for approximating the dispersive Serre–Green–Naghdi equations with topography using continuous finite elements. The method is an extension of the hyperbolic relaxation technique introduced in Guermond et al. (J Comput Phys 450:110809, 2022). It is explicit, second-order accurate in space, third-order accurate in time, and is invariant-domain preserving. It is also well balanced and parameter free. Special attention is given to the convex limiting technique when physical source terms are added in the equations. The method is verified with academic benchmarks and validated by comparison with laboratory experimental data.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Crystal structure and polymorphic forms of auranofin revisited

Auranofin, initially developed as a treatment for rheumatoid arthritis, is currently under extensive investigation as a potential drug for various conditions, including cancer, bacterial infections, and parasitic infections. The compound is a known inhibitor of thioredoxin reductase (TXNRD1) and related selenoproteins. Although preliminary studies on the auranofin crystal polymorphism exist, and a low-quality crystal structure has been reported, a comprehensive crystallographic characterization remains unexplored. Utilizing X-ray crystallography techniques, we conducted detailed structural analysis of auranofin and compared our findings with related organogold compounds. Implementation of Hirshfeld atom refinement (HAR) enabled a more accurate hydrogen atom positioning in the structure. The crystal packing reveals a layered arrangement stabilised by numerous weak hydrogen bonds and dispersive interactions. Notably, our attempts to reproduce the previously reported polymorphic form of auranofin, purportedly more water-soluble, were unsuccessful despite following published protocols. To our knowledge, this is the first study describing a “disappearing polymorph” phenomenon of any pharmaceutically relevant transition metal coordination compound. Our findings may have significant implications for medicinal chemistry and pharmacology of coordination complexes, suggesting the need for systematic revision of historical crystallographic data in this field.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Dynamic coupling and spin-wave dispersions in a magnetic hybrid system made of an artificial spin-ice structure and an extended NiFe underlayer

We present a combined experimental and numerical study of the spin-wave dispersion in a NiFe artificial spin-ice (ASI) system consisting of an array of stadium-shaped nanoislands deposited on the top of a continuous NiFe film with non-magnetic spacer layers of varying thickness. The spin-wave dispersion, measured by wavevector resolved Brillouin light scattering spectroscopy in the Damon–Eshbach configuration, consists of a rich number of modes, with either stationary or propagating character. We find that the lowest frequency mode displays a bandwidth of ~0.5 GHz, which is independent of the presence of the film underneath. On the contrary, the Brillouin light scattering intensity of some of the detected modes strongly depends on the presence of the extended thin-film underlayer. Micromagnetic simulations unveil the details of the dynamic coupling between the ASI lattice and film underlayer. Interestingly, the ASI lattice facilitates dynamics of the film either specific wavelengths or intensity modulation peculiar to the modes of the ASI elements imprinted in the film. Our results demonstrate that propagating spin waves can be modulated at the nanometer length scale by harnessing the dynamic mode coupling in the vertical, i.e., the out-of-plane direction of suitably designed magnonic structures.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Oblique propagation and temperature effects on the resonant right-hand ion beam instability

The resonant right-hand instability (RHI) is often the dominant mode driven by reflected ions upstream of Earth’s quasi-parallel bow shock. In the tradition of Peter Gary, this paper further explores the right-hand instability using numerical solutions of the plasma dispersion relation and non-linear kinetic simulations, with parameters inspired by observations from NASA’s Magnetospheric Multiscale (MMS) mission. Agreement is found between the ion distributions in the particle-in-cell simulations and Magnetospheric Multiscale spacecraft data, which show the gyrophase bunching characteristic of the instability. The non-linear structures created by right-hand instability tend to be stronger when the plasma beta is lower. These structures have sizes of around 100 to 200 ion inertial lengths perpendicular to the magnetic field, presenting planet-sized disturbances to the magnetosphere. 2d and 3D hybrid particle-in-cell simulations show that modes with a range of propagation angles oblique to the magnetic field are excited, providing a ground to understand previous statistical studies of observed foreshock waves.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

The HectoMAP Cluster Survey: Spectroscopically Identified Clusters and their Brightest Cluster Galaxies (BCGs)

We apply a friends-of-friends (FoF) algorithm to identify galaxy clusters and we use the catalog to explore the evolutionary synergy between brightest cluster galaxies (BCGs) and their host clusters. We base the cluster catalog on the dense HectoMAP redshift survey (2000 redshifts deg -2 ). The HectoMAP FoF catalog includes 346 clusters with 10 or more spectroscopic members within the range 0.05 < z < 0.55 and with a median z = 0.29. We list these clusters and their members. We also include central velocity dispersions (σ*, BCG ) for the FoF cluster BCGs, a distinctive feature of the HectoMAP FoF catalog. HectoMAP clusters with higher galaxy number density (80 systems) are all genuine clusters with a strong concentration and a prominent BCG in Subaru/Hyper Suprime-Cam images. The phase-space diagrams show the expected elongation along the line of sight. Lower-density systems include some low reliability systems. We establish a connection between BCGs and their host clusters by demonstrating that σ*, BCG /σ cl decreases as a function of cluster velocity dispersion (σ cl ), in contrast, numerical simulations predict a constant σ*, BCG /σ cl . Sets of clusters at two different redshifts show that BCG evolution in massive systems is slow over the redshift range z < 0.4. The data strongly suggest that minor mergers may play an important role in BCG evolution in clusters with σ cl ≳ 300 km s -1 . For lower mass systems (σ cl < 300 km s -1 ), major mergers may play a significant role. The coordinated evolution of BCGs and their host clusters provides an interesting test of simulations in high-density regions of the universe.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1976

Electromagnetic instabilities driven by unequal proton beams in the solar wind

The paper sets forth a numerical investigation of the linear dispersion relation for typical solar wind conditions at 1 AU during those times (high-speed streams) when a secondary beam of protons drifting relative to the main proton component is present. Three beam-driven instabilities were found to occur as the beam drift velocity approaches the Alfven speed: (1) a pure, field-aligned magnetosonic wave that is most important at relatively high beta and/or high beam drift speeds; (2) an oblique magnetosonic wave having highest growth rates 15-30 deg from the magnetic field; and (3) an oblique Alfven wave having maximum growth rates at increasing angle to the magnetic field. The linear growth rates for the field-aligned magnetosonic and the Alfven oblique modes are investigated as a function of relative beam density, varying anisotropic pitch angle distributions for the various components, electron temperature, and electron heat flux.

Montgomery, M. D.↗
NASA NTRS1984

Electromagnetic ion beam instabilities

The linear theory of electromagnetic instabilities driven by an energetic ion beam streaming parallel to a magnetic field in a homogeneous Vlasov plasma is considered. Numerical solutions of the full dispersion equation are presented. At propagation parallel to the magnetic field, there are four distinct instabilities. A sufficiently energetic beam gives rise to two unstable modes with right-hand polarization, one resonant with the beam, the other nonresonant. A beam with sufficiently large T (perpendicular to B)/T (parallel to B) gives rise to the left-hand ion cyclotron anisotropy instability at relatively small beam velocities, and a sufficiently hot beam drives unstable a left-hand beam resonant mode. The parametric dependences of the growth rates for the three high beam velocity instabilities are presented here. In addition, some properties at oblique propagation are examined. It is demonstrated that, as the beam drift velocity is increased, relative maxima in growth rates can arise at harmonics of the ion cyclotron resonance for both right and left elliptically polarized modes.

Gary, S. P.↗
NASA NTRS1985

Whistler damping at oblique propagation - Laminar shock precursors

This paper addresses the collisionless damping of whistlers observed as precursors standing upstream of oblique, low-Mach number terrestrial bow shocks. The linear theory of electromagnetic waves in a homogeneous Vlasov plasma with Maxwellian distribution functions and a magnetic field is considered. Numerical solutions of the full dispersion equation are presented for whistlers propagating at an arbitrary angle with respect to the magnetic field. It is demonstrated that electron Landau damping attenuates oblique whistlers and that the parameter which determines this damping is beta-e. In a well-defined range of parameters, this theory provides damping lengths which are the same order of magnitude as those observed. Thus electron Landau damping is a plausible process in the dissipation of upstream whistlers. Nonlinear plasma processes which may contribute to precursor damping are also discussed, and criteria for distinguishing among these are described.

Gary, S. P.↗
NASA NTRS1985

Coupling of newborn ions to the solar wind by electromagnetic instabilities and their interaction with the bow shock

The process by which the solar wind assimilates newly ionized atoms is important for understanding the presence of planetary or interstellar helium in the solar wind, the dynamics of the Active Magnetospheric Particle Tracer Explorers (AMPTE) lithium releases in front of the earth's bow shock, and the formation of cometary tails. In this paper is examined how newborn ions can be coupled to the solar wind in the direction parallel to the magnetic field by means of electromagnetic instabilities driven by the distribution of newborn ions. The linear properties of three instabilities are analyzed and compared with numerical solutions of the linear dispersion equation, while their nonlinear behavior is followed by means of computer simulation to obtain the characteristic time for the pickup process. With a primary emphasis on the AMPTE lithiuim releases, various degrees of realism are introduced into the calculations to model the upstream conditions and the intersection of the lithium with the bow shock. It is shown that a time-dependent shock model is needed to correctly reproduce the amount of lithium which is transmitted through the shock and that the resulting lithium ion distribution is still likely to be subject to the same type of instabilities in the magnetosheath. Applications of these results to comets, in particular the artificial comet expected to be generated by the AMPTE barium release in the magnetosheath, is also briefly discussed.

Winske, D.↗
NASA NTRS1985

Electromagnetic electron beam instabilities - Hot, isotropic beams

This paper considers the linear theory of electromagnetic instabilities driven by an electron beam in a homogeneous, nonrelativistic, Vlasov plasma. The beam is relatively hot, isotropic in its own frame, and streams parallel or antiparallel to a magnetic field B. Numerical solutions of the full dispersion equation for propagation parallel or antiparallel to B are presented, and the linear properties of the whistler heat flux and electron beam firehose instabilities are exhibited and compared. Under a broad range of parmameters the former mode has the lower beam speed threshold, and the larger maximum growth rate. In addition, it is demonstrated that, for a sufficiently large relative beam density, relative beam temperature, and plasma beta the whistler heat flux instability has a much lower beam speed threshold than the electrostatic electron beam instability. The application of these instabilities to first-order Fermi acceleration of electrons at space plasma shocks is discussed.

Gary, S. P.↗
NASA NTRS1986

Low-frequency waves in a high-beta collisionless plasma Polarization, compressibility and helicity

This paper considers the linear theory of waves near and below the ion cyclotron frequency in an isothermal electron-ion Vlasov plasma which is isotropic, homogeneous and magnetized. Numerical solutions of the full dispersion equation for the magnetosonic/whistler and Alfven/ion cyclotron modes at beta(i) = 1.0 are presented, and the polarizations, compressibilities, helicities, ion Alfven ratios and ion cross-helicities are exhibited and compared. At sufficiently large beta(i) and theta, the angle of propagation with respect to the magnetic field, the real part of the polarization of the Alfven/ion cyclotron wave changes sign, so that, for such parameters, this mode is no longer left-hand polarized. The Alfven/ion cyclotron mode becomes more compressive as the wavenumber increases, whereas the magnetosonic/whistler becomes more compressive with increasing theta.

Gary, S. P.↗
NASA NTRS1986

Oscillations of a viscous compound drop

The effect of viscosity on normal mode oscillations of a compound drop is determined for a range of values of the viscosity, and the effect of a host medium is considered. The general form of the dispersion relation is derived, and numerical values are obtained for the case of an annular liquid region composed of silicon oil and with the core and host regions taken to be air. Results indicate that the sloshing mode of the compound drop is more damped than the bubble mode.

Lyell, M. J.↗