SEARCH · Engineering Papers

Results for “Drift vector”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 127 records · Page 7

All sources NASA NTRS DOE OSTI

NASA NTRS2007

The Local Stellar Velocity Field via Vector Spherical Harmonics

We analyze the local field of stellar tangential velocities for a sample of 42,339 nonbinary Hipparcos stars with accurate parallaxes, using a vector spherical harmonic formalism. We derive simple relations between the parameters of the classical linear model (Ogorodnikov-Milne) of the local systemic field and low-degree terms of the general vector harmonic decomposition. Taking advantage of these relationships, we determine the solar velocity with respect to the local stars of (V(sub X), V(sub Y), V(sub Z)) (10.5, 18.5, 7.3) +/- 0.1 km s(exp -1) not corrected for the asymmetric drift with respect to the local standard of rest. If only stars more distant than 100 pc are considered, the peculiar solar motion is (V(sub X), V(sub Y), V(sub Z)) (9.9, 15.6, 6.9) +/- 0.2 km s(exp -1). The adverse effects of harmonic leakage, which occurs between the reflex solar motion represented by the three electric vector harmonics in the velocity space and higher degree harmonics in the proper-motion space, are eliminated in our analysis by direct subtraction of the reflex solar velocity in its tangential components for each star. The Oort parameters determined by a straightforward least-squares adjustment in vector spherical harmonics are A=14.0 +/- 1.4, B=13.1 +/- 1.2, K=1.1 +/- 1.8, and C=2.9 +/- 1.4 km s(exp -1) kpc(exp -1). The physical meaning and the implications of these parameters are discussed in the framework of a general linear model of the velocity field. We find a few statistically significant higher degree harmonic terms that do not correspond to any parameters in the classical linear model. One of them, a third-degree electric harmonic, is tentatively explained as the response to a negative linear gradient of rotation velocity with distance from the Galactic plane, which we estimate at approximately -20 km s(exp -1) kpc(exp -1). A similar vertical gradient of rotation velocity has been detected for more distant stars representing the thick disk (z greater than 1 kpc), but here we surmise its existence in the thin disk at z less than 200 pc. The most unexpected and unexplained term within the Ogorodnikov-Milne model is the first-degree magnetic harmonic, representing a rigid rotation of the stellar field about the axis -Y pointing opposite to the direction of rotation. This harmonic comes out with a statistically robust coefficient of 6.2 +/- 0.9 km s(exp -1) kpc(exp -1) and is also present in the velocity field of more distant stars. The ensuing upward vertical motion of stars in the general direction of the Galactic center and the downward motion in the anticenter direction are opposite to the vector field expected from the stationary Galactic warp model.

Galaxy↗

DOE OSTI2013

On plasma rotation induced by waves in tokamaks

The momentum conservation for resonant wave-particle interactions, now proven rigorously and for general settings, is applied to explain in simple terms how tokamak plasma is spun up by the wave momentum perpendicular to the dc magnetic field. The perpendicular momentum is passed through resonant particles to the dc field and, giving rise to the radial electric field, is accumulated as a Poynting flux; the bulk plasma is then accelerated up to the electric drift velocity proportional to that flux, independently of collisions. The presence of this collisionless acceleration mechanism permits varying the ratio of the average kinetic momentum absorbed by the resonant-particle and bulk distributions depending on the orientation of the wave vector. Both toroidal and poloidal forces are calculated, and a fluid model is presented that yields the plasma velocity at equilibrium.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1972

Direct measurements of plasma drift velocities at high magnetic latitudes.

Description of an incoherent scatter radar experiment performed at the 23-cm radar facility in Chatanika, Alaska. The experiment has provided a direct method for measuring the ionospheric plasma transport velocity vector over long periods with relatively good time resolution. Since the F-region transport at this site is closely associated with magnetospheric convection, particularly at night, the radar can provide important information about the behavior of the magnetosphere.

Doupnik, J. R.↗

NASA NTRS1977

Technique for measuring the differential ion flux vector

A diagnostic technique is discussed which gives the angle of incidence of a plasma stream, the energy corresponding to the mean velocity of the ions, and the distribution of the ion thermal motion superimposed on the drift. The technique is shown to be operable for low-energy plasma streams (5-60 eV) and in the presence of multiple plasma streams differing in direction and/or energy. Resolution is better than 3.5 deg with approximately 1% energy spread in the streams. The instrument samples the plasma at a single region in space, independent of the number of plasma streams or their characteristics. The technique was developed to investigate plasma flow interactions in the laboratory. Its capabilities are demonstrated by some preliminary data taken for the case of a long cylindrical body immersed in a drifting, collisionless plasma.

Stone, N. H.↗

DOE OSTI2017

Large scale Brownian dynamics of confined suspensions of rigid particles

We introduce methods for large-scale Brownian Dynamics (BD) simulation of many rigid particles of arbitrary shape suspended in a fluctuating fluid. Our method adds Brownian motion to the rigid multiblob method [F. Balboa Usabiaga et al., Commun. Appl. Math. Comput. Sci. 11(2), 217–296 (2016)] at a cost comparable to the cost of deterministic simulations. We demonstrate that we can efficiently generate deterministic and random displacements for many particles using preconditioned Krylov iterative methods, if kernel methods to efficiently compute the action of the Rotne-Prager-Yamakawa (RPY) mobility matrix and its “square” root are available for the given boundary conditions. These kernel operations can be computed with near linear scaling for periodic domains using the positively split Ewald method. Here we study particles partially confined by gravity above a no-slip bottom wall using a graphical processing unit implementation of the mobility matrix-vector product, combined with a preconditioned Lanczos iteration for generating Brownian displacements. We address a major challenge in large-scale BD simulations, capturing the stochastic drift term that arises because of the configuration-dependent mobility. Unlike the widely used Fixman midpoint scheme, our methods utilize random finite differences and do not require the solution of resistance problems or the computation of the action of the inverse square root of the RPY mobility matrix. We construct two temporal schemes which are viable for large-scale simulations, an Euler-Maruyama traction scheme and a trapezoidal slip scheme, which minimize the number of mobility problems to be solved per time step while capturing the required stochastic drift terms. We validate and compare these schemes numerically by modeling suspensions of boomerang-shaped particles sedimented near a bottom wall. Using the trapezoidal scheme, we investigate the steady-state active motion in dense suspensions of confined microrollers, whose height above the wall is set by a combination of thermal noise and active flows. We find the existence of two populations of active particles, slower ones closer to the bottom and faster ones above them, and demonstrate that our method provides quantitative accuracy even with relatively coarse resolutions of the particle geometry.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1993

Whistlers and plasmaspheric hiss - Wave directions and three-dimensional propagation

Wave propagation directions are determined on the basis of wave data from the DE 1 satellite showing simultaneously nonducted whistlers and hiss. Hiss wave normal angles are determined as about 70 and 77 deg for f = 3.5 and 2.5 kHz, respectively, with the wave vector being almost perpendicular to the meridional plane. A novel approximate analytical formulation of 3D propagation of whistler waves is developed and used to model the drift of magnetospherically reflected whistlers in azimuth. It is shown that depending on initial parameters, the time of arrival of whistler rays at a fixed observation point can differ by 10-20 s, with signals from different magnetospherically reflected whistlers overlapping to evolve into a hisslike signal. The total azimuthal drift of whistler rays is found to not exceed about 30 deg, so that plasmaspheric hiss may be produced by nonducted whistlers at longitudes correlated with the location of thunderstorm activity.

Draganov, A. B.↗

NASA NTRS2003

Measurements of cell-generated deformations on flexible substrata using correlation-based optical flow

The optical flow algorithm presented here is a robust method that rapidly yields a high-density field of substrate displacement vectors based on two optical images. We found that one of the limiting factors, at least for inexperienced experimentalists, is the consistency of focusing or the drift in microscope focus. However, with properly collected images the standard error of the measurement was estimated to be on the order of +/- 0.10 pixels. Finally, although the discussion has been focused on the displacement of flexible substrata, a similar method should be applicable for detecting movements on other types of images, as long as the movement involves a certain degree of local coordination.

Non-NASA Center↗

DOE OSTI2016

Nonlinear electromagnetic formulation for particle-in-cell simulation of lower hybrid waves in toroidal geometry

An electromagnetic particle simulation model has been formulated and verified for nonlinear processes of lower hybrid (LH) waves in fusion plasmas. Electron dynamics are described by the drift kinetic equation using either kinetic momentum or canonical momentum. Ion dynamics are treated as the fluid system or by the Vlasov equation. Compressible magnetic perturbation is retained to simulate both the fast and slow LH waves. Numerical properties are greatly improved by using the electron continuity equation to enforce the consistency between electrostatic potential and vector potential, and by using the importance sampling scheme. The simulation model has been implemented in the gyrokinetic toroidal code (GTC), and verified for the dispersion relation and nonlinear particle trapping of the electromagnetic LH waves.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2016

Hydrodynamic description of an unmagnetized plasma with multiple ion species. I. General formulation

A generalization of the Braginskii ion fluid description [S. I. Braginskii, Sov. Phys. JETP 6, 358 (1958)] to the case of an unmagnetized collisional plasma with multiple ion species is presented. An asymptotic expansion in the ion Knudsen number is used to derive the individual ion species continuity, as well as the total ion mass density, momentum, and energy evolution equations accurate through the second order. Expressions for the individual ion species drift velocities with respect to the center of mass reference frame, as well as for the total ion heat flux and viscosity, which are required to close the fluid equations, are evaluated in terms of the first-order corrections to the lowest order Maxwellian ion velocity distribution functions. A variational formulation for evaluating such corrections and its relation to the plasma entropy are presented. Employing trial functions for the corrections, written in terms of expansions in generalized Laguerre polynomials, and maximizing the resulting functionals produces two systems of linear equations (for “vector” and “tensor” portions of the corrections) for the expansion coefficients. A general matrix formulation of the linear systems as well as expressions for the resulting transport fluxes are presented in forms convenient for numerical implementation. The general formulation is employed in the companion paper [A. N. Simakov and K. Molvig, Hydrodynamic description of an unmagnetized plasma with multiple ion species. II. Two and three ion species plasmas, submitted to Phys. Plasmas (2015)] to evaluate the individual ion drift velocities and the total ion heat flux and viscosity for specific cases of two and three ion species plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2018

Symmetries of a reduced fluid-gyrokinetic system

Symmetries of a fluid-gyrokinetic model are investigated using Lie group techniques. Specifically the nonlinear system constructed by Zocco and Schekochihin (Zocco & Schekochihin 2011), which combines nonlinear fluid equations with a drift-kinetic description of parallel electron dynamics, is studied. Significantly, this model is fully gyrokinetic, allowing for arbitrary $k_{\bot \rho I}$ where $k_\bot$ is the perpendicular wave vector of the fluctuations and $\rho I}$ is the ion gyroradius. The model includes integral operators corresponding to gyroaveraging as well as the moment equations relating fluid variables to the kinetic distribution function. A large variety of exact symmetries is uncovered, some of which have unexpected form. Using these results, new nonlinear solutions are constructed, including a helical generalization of the Chapman-Kendall solution for a collapsing current sheet.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2019

Particle filter-based modulation domain infrared targets tracking

Faced with problems of low contrast, poor SNR, and relatively complicated tracking environment, stable infrared target tracking is worth researching for its many potential applications. In this paper, instead of traditional target tracking in the pixel domain, we propose a sampling importance resampling (SIR) particle filter method with indirect velocity measurements to track infrared targets in the modulation domain. The dominant amplitude modulation (AM) features used for tracking is extracted by decomposing the input image using an 18-channel Gabor filter bank followed by the application of the dominant component analysis approach. The dominant AM modulation features provide a significant partial texture characteristic of the target which can be separated from background with better discrimination. To take advantage of observed kinematics, we utilize the augmented state vector with indirect velocity information via combining the measurements of velocity in adjacent frames to the SIR particle filter framework, which weakens weights of particles with bad velocity estimates but still having association with the cluttered background or other moving objects. A dynamic template update strategy is also provided to prevent the tracker from appearance model drift. Experiments indicate that the proposed method is effective for raising the tracking accuracy compared with other tracking methods.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2016

SU-D-207A-07: The Effects of Inter-Cycle Respiratory Motion Variation On Dose Accumulation in Single Fraction MR-Guided SBRT Treatment of Renal Cell Carcinoma

Purpose: To assess the dose deposition in simulated single-fraction MR-Linac treatments of renal cell carcinoma, when inter-cycle respiratory motion variation is taken into account using online MRI. Methods: Three motion characterization methods, with increasing complexity, were compared to evaluate the effect of inter-cycle motion variation and drifts on the accumulated dose for an SBRT kidney MR-Linac treatment: 1) STATIC, in which static anatomy was assumed, 2) AVG-RESP, in which 4D-MRI phase-volumes were time-weighted, based on the respiratory phase and 3) PCA, in which 3D volumes were generated using a PCA-model, enabling the detection of inter-cycle variations and drifts. An experimental ITV-based kidney treatment was simulated in a 1.5T magnetic field on three volunteer datasets. For each volunteer a retrospectively sorted 4D-MRI (ten respiratory phases) and fast 2D cine-MR images (temporal resolution = 476ms) were acquired to simulate MR-imaging during radiation. For each method, the high spatio-temporal resolution 3D volumes were non-rigidly registered to obtain deformation vector fields (DVFs). Using the DVFs, pseudo-CTs (generated from the 4D-MRI) were deformed and the dose was accumulated for the entire treatment. The accuracies of all methods were independently determined using an additional, orthogonal 2D-MRI slice. Results: Motion was most accurately estimated using the PCA method, which correctly estimated drifts and inter-cycle variations (RMSE=3.2, 2.2, 1.1mm on average for STATIC, AVG-RESP and PCA, compared to the 2DMRI slice). Dose-volume parameters on the ITV showed moderate changes (D99=35.2, 32.5, 33.8Gy for STATIC, AVG-RESP and PCA). AVG-RESP showed distinct hot/cold spots outside the ITV margin, which were more distributed for the PCA scenario, since inter-cycle variations were not modeled by the AVG-RESP method. Conclusion: Dose differences were observed when inter-cycle variations were taken into account. The increased inter-cycle randomness in motion as captured by the PCA model mitigates the local (erroneous) hotspots estimated by the AVG-RESP method.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2025

Measurements of the polarization of several instabilities in the DIII-D tokamak

Recently, a method to infer the polarization of modes with frequencies much less than the ion cyclotron frequency was published [X.D. Du et al., Phys. Rev. Lett. 132 (2024) 215101]. The method uses measurements of electron temperature and density fluctuations δTₑ and δnₑ at the same spatial position to infer the local ratio of “acoustic polarization,” |δϕ ∥ |/(|δϕ ∥ |+|δψ|), where δϕ ∥ is the effective parallel potential and δψ is related to the parallel magnetic vector potential A ∥ . This paper summarizes key formulas, with emphasis on their range of validity, and elaborates on the workflow required to infer the acoustic polarization from experimental data. The drift-acoustic polarization of ellipticity-induced, toroidicity-induced, and reversed shear Alfvén eigenmodes is nearly zero, as expected for modes with predominately shear-Alfvénic polarization. The polarization of beta-induced Alfvén eigenmodes contains an acoustic component that increases with poloidal wave number. “Low frequency modes,” (instabilities that appear transiently when the minimum of the safety factor qₘᵢₙ passes through rational values) have large and highly variable acoustic polarization. In both experiment and simulation, fishbones have non-zero acoustic polarization that increases as the mode chirps down in frequency.

Alfven eigenmode↗

NASA NTRS1977

Nonlinear evolution of the collisionless tearing mode

For a broad k vector spectrum, the collisionless tearing mode grows linearly until the Landau electrons become magnetized by the wave field. The electrons then spatially diffuse as a result of turbulent magnetic-gradient drifts which leads to a nonlinearly enhanced growth rate. At large amplitudes the nonlinear growth rate is dominated by ions.

Coroniti, F. V.↗

NASA NTRS1986

Transfer function analysis of thermospheric perturbations

Applying perturbation theory, a spectral model in terms of vectors spherical harmonics (Legendre polynomials) is used to describe the short term thermospheric perturbations originating in the auroral regions. The source may be Joule heating, particle precipitation or ExB ion drift-momentum coupling. A multiconstituent atmosphere is considered, allowing for the collisional momentum exchange between species including Ar, O2, N2, O, He and H. The coupled equations of energy, mass and momentum conservation are solved simultaneously for the major species N2 and O. Applying homogeneous boundary conditions, the integration is carred out from the Earth's surface up to 700 km. In the analysis, the spherical harmonics are treated as eigenfunctions, assuming that the Earth's rotation (and prevailing circulation) do not significantly affect perturbations with periods which are typically much less than one day. Under these simplifying assumptions, and given a particular source distribution in the vertical, a two dimensional transfer function is constructed to describe the three dimensional response of the atmosphere. In the order of increasing horizontal wave numbers (order of polynomials), this transfer function reveals five components. To compile the transfer function, the numerical computations are very time consuming (about 100 hours on a VAX for one particular vertical source distribution). However, given the transfer function, the atmospheric response in space and time (using Fourier integral representation) can be constructed with a few seconds of a central processing unit. This model is applied in a case study of wind and temperature measurements on the Dynamics Explorer B, which show features characteristic of a ringlike excitation source in the auroral oval. The data can be interpreted as gravity waves which are focused (and amplified) in the polar region and then are reflected to propagate toward lower latitudes.

Mayr, H. G.↗

NASA NTRS1972

Correlation techniques and measurements of wave-height statistics

Statistical measurements of wave height fluctuations have been made in a wind wave tank. The power spectral density function of temporal wave height fluctuations evidenced second-harmonic components and an f to the minus 5th power law decay beyond the second harmonic. The observations of second harmonic effects agreed very well with a theoretical prediction. From the wave statistics, surface drift currents were inferred and compared to experimental measurements with satisfactory agreement. Measurements were made of the two dimensional correlation coefficient at 15 deg increments in angle with respect to the wind vector. An estimate of the two-dimensional spatial power spectral density function was also made.

Guthart, H.↗

DOE OSTI2022

Primary track recovery in high-definition gas time projection chambers

Abstract We develop and validate a new algorithm called primary track recovery () that effectively deconvolves known physics and detector effects from nuclear recoil tracks in gas time projection chambers (TPCs) with high-resolution readout. This gives access to the primary track charge, length, and vector direction (helping to resolve the “head-tail” ambiguity). Additionally, provides a measurement of the transverse and longitudinal diffusion widths, which can be used to determine the absolute position of tracks in the drift direction for detector fiducialization. Using simulated helium recoils in an atmospheric pressure TPC with a 70:30 mixture of $$\hbox {He:CO}_2$$ He:CO 2 we compare the performance of to traditional methods for all key track variables. We find that the algorithm reduces reconstruction errors, including those caused by charge integration, for tracks with mean length-to-width ratios 1.4 and above, corresponding to recoil energies of 20 keV and above in the studied TPCs. We show that improves on existing methods for head-tail disambiguation, particularly for highly inclined tracks, and improves the determination of the absolute position of recoils on the drift axis via transverse diffusion. We find that can partially recover charge structure integrated out by the detector in the z direction, but that its determination of energy and length have worse resolution compared to existing methods. We use experimental data to qualitatively verify these findings and discuss implications for future directional detectors at the low-energy frontier.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2014

Matching to Drift Spaces Containing Acceleration

The transfer matrix of a drift space of length L is shown in the equation below: [$^x_{x^{\prime}}$] $(L)$ = [$^1_0$ $^L_1$][$^x_{x^{\prime}}$]$(0$), with $y$ and $y^{\prime}$ behaving behaving similarly. Writing the phase space vector as $s$ = [$^x_{x^{\prime}}$], it can be seen that under a linear mapping $s \mapsto Ms$. the matrix $\langle{ss^T}\rangle$, where the average is taken over a beam of particles, transforms as $\langle{ss^T}\rangle$ $\mapsto$ $\langle{Mss^TM^T}\rangle$ = $M\langle{ss^T}\rangle$, assuming all particles have the same transfer matrix $M$ (this is true if, for example, there is no momentum spread).

43 PARTICLE ACCELERATORS↗