Search NASA⌕ Search

SEARCH · Search NASA

Results for “motion study”

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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Numerical simulations of three-dimensional ion crystal dynamics in a Penning trap using the fast multipole method

We simulate the dynamics, including laser cooling, of three-dimensional (3-D) ion crystals confined in a Penning trap using a newly developed molecular dynamics-like code. The numerical integration of the ions’ equations of motion is accelerated using the fast multipole method to calculate the Coulomb interaction between ions, which allows us to efficiently study large ion crystals with thousands of ions. In particular, we show that the simulation time scales linearly with ion number, rather than with the square of the ion number. By treating the ions’ absorption of photons as a Poisson process, we simulate individual photon scattering events to study laser cooling of 3-D ellipsoidal ion crystals. Initial simulations suggest that these crystals can be efficiently cooled to ultracold temperatures, aided by the mixing of the easily cooled axial motional modes with the low frequency planar modes. In our simulations of a spherical crystal of 1000 ions, the planar kinetic energy is cooled to several millikelvin in a few milliseconds while the axial kinetic energy and total potential energy are cooled even further. This suggests that 3-D ion crystals could be well suited as platforms for future quantum science experiments.

Zaris, John (ORCID:0009000196476323)↗
DOE OSTI2024

Tropical Cyclone Wind Shear-Relative Asymmetry in Reanalyses

Abstract While tropical cyclones (TCs) are axisymmetric vortices to the first order, they often exhibit noteworthy structural asymmetries. These often result from environmental vertical wind shear, which tilts the vortex and induces a wavenumber 1 pattern in the circulation and precipitation fields. Reanalyses and climate models have improved in representing the TC structure and climatology, but their relatively coarse resolution and dependence on parameterized physics cast doubt on their ability to capture the asymmetric TC structure. We perform the most comprehensive process-oriented assessment of TC asymmetry to date in reanalyses. Specifically, we analyze the composite shear-relative TC structure in ERA5 and Climate Forecast System Reanalysis (CFSR), which vary in their resolutions, physical parameterization suites, and data assimilation techniques. These structures are compared with aircraft reconnaissance radar observations. In agreement with the observations, the strongest tangential winds are usually found left-of-shear, while inner core rainfall, ascent, vortex tilt, and low-level inflow are favored directly downshear or in the downshear-left quadrant. Outer rainband convection generally peaks in the downshear-right quadrant. Thermodynamic asymmetries are also apparent, with anomalous low-level moisture right-of-shear, midlevel warmth in the upshear-right quadrant (uptilt), and cloud properties suggestive of a realistic precipitation life cycle from growth to fallout. We also decompose rainfall contributions from the convective parameterization and large-scale cloud schemes and highlight the roles of vorticity advection, buoyancy advection, and diabatic processes in driving asymmetric vertical motions in the inner core and outer rainband regions. Our results suggest that process-level studies of TC asymmetry and TC–wind shear interaction under future warming are viable using climate models. Significance Statement Asymmetries are common in tropical cyclones (TCs), influencing their intensity, track, and hazards. Vertical wind shear often plays a leading-order role in causing these asymmetries. It is uncertain how well asymmetric structures and processes are captured in reanalyses and global climate models (GCMs) with grid spacings of 0.25° and coarser. In this study, we first evaluate TC asymmetry in reanalyses, which have the benefit of being forced by observations. This helps to assess whether the resolutions associated with GCMs sufficiently capture asymmetric structures and processes and motivates upcoming work with free-running GCMs to study how TC asymmetry may change in a warming climate.

Carstens, Jacob D.↗
DOE OSTI2025

Investigation into the instantaneous centre of rotation for enhanced design of floating offshore wind turbines

The dynamic behaviour of floating offshore wind turbines (FOWTs) involves complex interactions of multivariate loads from wind, waves, and currents, which result in complex motion characteristics. Although methods for analysing global motion responses are well-established, the time- and location-dependent kinematics remain underexplored. This paper investigates the instantaneous centre of rotation (ICR), a point of zero velocity at a time instance of general plane motion. Understanding and strategically positioning the ICR can reduce the dynamic motion in critical structural locations, enhancing the performance and structural robustness of FOWTs. The paper presents a method for computing the ICR using time-domain simulation results and proposes a statistical analysis approach suitable for design studies. Building on prior research, it examines the sensitivity of the ICR to external loading and design features, providing insights into how these factors influence motion response and how the motion response influences the statistics of the ICR, structural loads, and other performance metrics of interest. The study explores two FOWT configurations, a spar and a semisubmersible, identifying design variables that most effectively control the ICR statistics and identifying the ICR statistics most correlated with the responses of interest. Finally, through two case studies, we demonstrate how to apply these new insights in a practical design scenario. By adjusting the design variables most correlated with the ICR (fairlead vertical position and centre of mass for the spar and mooring line length and offset column diameter for the semisubmersible), we successfully modified the designs of the floating support structures to reduce the loads in the mooring lines, tower base, and blade roots, improving the ultimate strength and fatigue characteristics compared to the original designs.

17 WIND ENERGY↗
DOE OSTI2020

Characterizing the Small-Scale Dynamical, Ice Microphysical, and Residual Aerosol Properties of Mid-Latitude Cold Clouds: A Pilot Study (Final Report)

Our understanding of cirrus clouds has been advancing rapidly due to improvements in measurements made inside clouds, computer cloud models, and laboratory experimentation. This advancement in knowledge has made it clear that cirrus are complicated systems where the intricate interactions among dynamic air motions, solar and terrestrial radiation, ice crystal growth, and water and air chemistry can substantially influence cirrus structure and evolution. Measurements made inside clouds have shown that cirrus crystals are often complex with shapes that are often unclassifiable, especially at small sizes. Laboratory measurements have revealed that the ice crystal surface is complex, even at small sizes. Recent laboratory and satellite studies show that this surface complexity modifies the radiative properties of cirrus clouds in important ways. Moreover, observations and computer cloud modeling work show that the larger and smaller scale properties of cirrus are sensitive to the interplay between small-scale air motions and the aerosol articles upon which cirrus ice crystals form. Taken together, these studies indicate that accurately accounting for the radiative properties of cirrus clouds requires knowing the air motions and aerosol properties that influence the number concentration and crystal size, along with an understanding of the surface properties of vapor grown ice. Although this prior work has been illuminating, there have been few studies that provide detailed shape and surface complexity measurements for the small ice crystals that often compose cirrus clouds.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Investigation of Site Amplifications Using Ambient-Noise-Derived Shallow Velocity Structures Under a Dense Array in Oklahoma

The shear-wave velocity (V S ) structure plays an important role in characterizing site amplification. The Large-n Seismic Survey in Oklahoma (LASSO; 1820 stations) revealed large vertical ground-motion variability in a 25 km × 32 km area in northern Oklahoma. The LASSO array has a relatively simple and flat topography, typical in a sedimentary basin environment in the central United States. In this study, we use the dense array to investigate the velocity structure under the LASSO array and how vertical ground motions relate to the shallow-to-deep structures. We extract the fundamental-mode Rayleigh wave by cross-correlating one month of ambient noise (0.7–5 Hz). We use double-beamforming to measure the group and phase velocities and anisotropy. By jointly inverting the group and phase velocities, we obtain the V S structure. Here, we observe correlations between V S at depths of 0.1–1.5 km and vertical ground motions using sites on the stiffer Permian formations. The shallow Quaternary alluvium and terrace deposits can amplify vertical ground motions by a factor of 2–4.5 between 2 and 25 Hz and attenuate signals above 25 Hz. We use 1D V S profiles to simulate the SV-wave transfer functions. An average V S of 250 m/s in the upper 20–40 m may cause the observed amplification between 2 and 40 Hz. V S estimated by topographic slopes cannot predict the relative amplification. Our results highlight the large variability of site-dependent ground motion in a small local region and the importance of characterizing shallow structures to estimate seismic hazards. Small thickness variations of the shallow formation can significantly change the resonance amplitude and frequency, which likely reduces the coherency of the Rayleigh waves extracted from ambient noise.

Chang, Hilary [Massachusetts Inst. of Technology (↗
DOE OSTI2024

Joint Inversion of Regional Waveform, First-Motion Polarity, and Synthetic Aperture Radar Surface Displacement for the Fourth and Sixth North Korean Declared Nuclear Explosions

Here, this study analyzed the Democratic People’s Republic of Korea’s (DPRK) fourth (DPRK4, 6 January 2016 M w 4.49) and sixth (DPRK6, 7 September 2017 M w 5.2) declared nuclear tests, employing a joint seismic and Interferometric Synthetic Aperture Radar (InSAR) inversion to improve understanding of these events and enhance moment tensor (MT) inversion capabilities. The recent efforts have focused on employing seismic waveform and InSAR geodetic deformation data separately to analyze these and the previous nuclear tests (e.g., Chiang et al., 2018; Myers et al., 2018; Wang et al., 2018). Building upon our previous work (Chi-Durán et al., 2021), we performed a joint regional waveform, first-motion (FM) polarity, and surface displacement inversion, which demonstrated improved source-type discrimination, a revised MT solution with reduced scalar moment uncertainty, and an independently constrained location. In this article, we build on the previous results for DPRK6 by including an analysis using a four-layered velocity model with free-surface topography to compute the near-source static deformation Green’s functions. The model consists of a 50 m basalt layer (⁠V P = 2.07 km/s, V S = 1.2 km/s⁠), a 250 m stratified volcanic deposit layer (⁠V P = 1.73 km/s, V S = 1.0 km/s⁠), a 700 m weathered granodiorite layer (⁠V P = 2.5 km/s, V S = 1.3 km/s⁠⁠), and a granodiorite half-space (⁠V P = 5.35 km/s, V S = 3.09 km/s⁠⁠). The half-space shares the velocity of the regional MDJ2 velocity model (Ford et al., 2010), which has proven effective for waveform inversion in the region. This model considers the range of reported values for various lithologies and weathering effects. Our findings show that using the layered velocity model enhances the recovery of source location and depth for both the explosions by improving fits and reducing uncertainties. The joint inversion also improves source-type discrimination and better constrains the scalar seismic moment necessary for downstream yield estimation.

58 GEOSCIENCES↗
DOE OSTI2024

Assessing the Expansion of Ground-Motion Sensing Capability in Smart Cities via Internet Fiber-Optic Infrastructure

Monitoring ground motion in smart cities can improve the public safety by providing critical insights on natural and anthropogenic hazards, for example, earthquakes, landslides, explosions, infrastructure failures, and so forth. Although seismic activity is typically measured using dedicated point sensors (e.g., geophones and accelerometers), techniques such as distributed acoustic sensing have demonstrated the utility of using fiber-optic cable to detect seismic activity over comparable distances. In this article, we present the results of a study that quantifies the expansion in an area monitored for low-amplitude ground-motion events by augmenting existing point sensors with the internet fiber-optic cable infrastructure. Here we begin by describing our methodology, which utilizes geospatial data on point sensors and internet optical fiber deployed in metropolitan statistical areas (MSAs) in the United States. We extend these data to identify the area that can be monitored by (1) considering the observed seismic noise data in target locations, (2) applying the model from Wilson et al. (2021) to understand the potential coverage area gains using optical fiber sensing, and (3) optimizing the selection of fiber segments to maximize coverage and minimize deployment costs. We implement our methodology in ArcGIS to assess the additional area that can be monitored for low-amplitude ground-motion events (i.e., magnitude >0.5) by utilizing internet fiber-optic cables in the 100 most populous MSAs in the United States. We find that the addition of internet fiber-based sensors in MSAs would increase the area monitored on average by over an order of magnitude from 1% to 12%, if the subset of fiber cable segments that maximize coverage and minimize deployment costs is chosen even if only 20% of all fibers are used.

58 GEOSCIENCES↗
DOE OSTI2022

Energizing charged particles by an orbit instability in a slowly rotating magnetic field

In this work, the stability of charged particle motion in a uniform magnetic field with an added spatially uniform transverse rotating magnetic field (RMF) is studied analytically. It is found that the stability diagram of a single particle's orbit depends critically on the chosen boundary conditions. We show that for many boundary conditions and wide regions in the parameter space, RMFs oscillating far below the cyclotron frequency can cause linear instabilities in the motion which break μ invariance and energize particles. Such energization may appear at odds with the adiabatic invariance of μ; however, adiabatic invariance is an asymptotic result and does not preclude such heating by magnetic fields oscillating at slow frequencies. This mechanism may contribute to heating in the edge plasma of field-reversed configurations (FRCs) in rotamak-FRC experiments. Furthermore, these RMF-driven instabilities may significantly enhance azimuthal current drive during the formation of FRCs in such devices.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

On the origin of internal obstacles to dislocation glide in single-phase NiFe random alloys

Over the past decade, a new class of unconventional alloys, where several constitutive elements are presented in close to equal concentrations, is considered as promising structural materials in many applications. Properties and behavior of many concentrated solid solution alloys (CSSAs) under different conditions, including mechanical loading, differs from that of conventional alloys and has been the subject of intensive studies by different techniques. A key feature of these materials is that there is no clear distinction between solute and solvent species of atoms, which is the conventional approach to considering dislocation motion in alloys. To understand this regime, we have recently reported on the glide of a screw dislocation glide in a simple equiatomic NiFe alloy and reported the two-mode character. At low applied stress, the dislocation motion is similar to a glide through a field of obstacles. In this work, we report the results of a detailed study on the nature of these internal obstacles. We demonstrate that these are not conventional localized “obstacles” but is the result of significantly different energies of instantaneous dislocation configurations over the alloy bulk. The strong pinning positions are associated with composition changes in the partial dislocation cores, while observing no significant composition changes along the stacking faults. Local energies can vary significantly, as demonstrated using simulations of short-periodicity dislocation lines. Therefore, the critical stress to move a dislocation under such conditions is not a constant Peierls stress but represents a wide spectrum of barriers between different non-straight valleys that the dislocation line accommodates in the bulk at particular stress/strain conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2020

STRUCTURE AND PROPERTIES STUDY ON ENERGY MATERIALS: THERMOELECTRIC MATERIAL TETRAHEDRITE AND LITHIUM ION CONDUCTOR LiPON

Development of efficient energy materials is critical in order to ease the energy demand and reduce our dependence of fossil fuel. Thermoelectric materials are promising due to their capability of generating electrical power by recovering waste heat. The performance of thermoelectric materials is quantified by a dimensionless figure of merit zT, which depends on their properties such as electrical conductivity, Seebeck coefficient and thermal conductivity. Tetrahedrites, a copper antimony sulfosalt mineral, typified by Cu12-xMxSb4S13, where M is a transition metal element such as Ni, Zn, Fe or Mn, have great potential for thermoelectric application due to their relatively high zT (close to 1 at 700 K), earth-abundance, environmental friendliness, favorable electrical properties, and most importantly intrinsic low lattice thermal conductivity (less than 1 W m-1 K-1) in wide temperature. In addition to energy recovery, reliable energy storage devices are also emerging to relieve the energy demand and improve the efficiency of consuming energy resources. Lithium-ion batteries are known to be reliable and successful electrochemical energy storage devices and appliable in various aspects, including laptops, smartphones and electrical vehicles. Lithium phosphorous oxynitride (LiPON) are widely used as thin-film solid-state electrolytes in Li-ion battery, which is the only demonstrated solid-state electrolyte that is quite stable in direct contact with Li metal at potentials from 0-5 V. However, the structure of LiPON, the effects of N doping, and the origin of its good electrochemical stability remains inconclusive.In this thesis, reliable modeling techniques accompanied with experimental tools, are applied to study the thermoelectric material tetrahedrite and the ionic conductor LiPON, in order to study their structural and dynamical properties. Accurate and efficient density-functional theory (DFT) and density-functional tight-binding (DFTB) methods, combined with molecular dynamics (MD) simulations are utilized in order to investigate the structures and properties of these energy materials. The incoherent and coherent atomic dynamics study of tetrahedrite Cu10.5NiZn0.5Sb4S13 provides the origin of softening upon cooling by investigate the motion of Cu12e at different temperatures. The dynamic structure factors in the longitudinal and transverse direction will also be discussed. The Cu movement of Cu-rich tetrahedrite Cu14Sb4S13 is revealed by Cu self-diffusivity, nuclear density map and "nudged elastic band" (NEB). Moreover, we investigate the effect of simulation cell size and basis sets on the DFT-based MD simulation results using tetrahedrite Cu10Zn2Sb4S13 thermoelectric as a model material, showing the advantage of larger cell by accessing smaller Q range. In addition, the low-temperature structural properties of Cu12Sb4S13 is measured by neutron diffraction, which indicates that no cubic to tetragonal transition occurs at metal-semiconductor transition (MST) temperature. Thermoelectric properties such as Seebeck coefficient, electrical resistivity and electrical thermal conductivity will also be investigated. DFTB method is implemented to study the structure and transport properties of Li3PO4 and LiPON, while the exploration of N doping effect is included. Lastly, the LiPON/Li interphase will be revealed in order to study the origin of electrochemical stability

25 ENERGY STORAGE↗
DOE OSTI2021

Sensitive parameter identification and uncertainty quantification for the stability of pipeline conveying fluid

In this study, several uncertainty quantification and sensitivity analysis methods are used to determine the most sensitive geometric and material input parameters of a cantilevered pipeline conveying fluid when uncertainty is introduced to the system at the onset of instability. The full nonlinear equations of motion are modeled using the extended Hamilton’s principle and then discretized using Galerkin’s method. A parametric study is first performed, and the Morris elementary effects are calculated to obtain a preliminary understanding of how the onset speed changes when each parameter is introduced to a ± 5% uncertainty. Then, four different input uncertainty distributions, mainly, uniform and Gaussian distribution, are chosen to investigate how input distributions affect uncertainty in the output. A convergence analysis is used to determine the number of samples needed to maintain simulation accuracy while saving the most computational time. Then, Monte Carlo simulations are run, and the output distributions for each input distribution at ± 1%, ± 3% and ± 5% input uncertainty range are found and discussed. Additionally, the Pearson correlation coefficients are evaluated for different uncertainty ranges. A final Monte Carlo study is performed in which single parameters are held constant while all others still have uncertainty. Overall, the flow speed at the onset of instability is the most sensitive to changes in the outer diameter of the pipe.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Observer-Based Nonlinear Control Scheme to Reduce Oscillations and Zero Crossing in Skid-Steer Vehicles

Motion sickness is a common condition experienced by drivers of skid-steer vehicles, primarily caused by zero crossing and oscillations in undamped systems. This study proposes an observer-based nonlinear control scheme to reduce transient oscillations and zero-crossing phenomena in skid-steer vehicles, thereby potentially alleviating motion sickness. Reducing transient oscillations and zero crossing in the transient response may alleviate motion sickness. A nonlinear damping controller is designed to improve transient response by reducing oscillations and zero-crossing. To design the controller, a reduced-order kinematic model based on coordinate transformation is developed. This transformation not only converts the system modeling into a controllable form but also enhances control performance. Modeling error is addressed by considering the distance between the center of the vehicle and the sensor location. Despite these improvements, model uncertainties and external disturbances remain, which may degrade control performance. To ensure robustness and estimate such disturbances, a high-order sliding mode observer (HOSMO) is incorporated. The effectiveness of the proposed method is validated through MATLAB/Simulink and TruckMaker simulations. From the simulation results, it was shown that the proposed method reduced the mean squared error of the tracking error to within 10 % compared to the state feedback controller with the HOSMO.

Seo, Jiwon [Chung-Ang University, Seoul (Korea, Re↗
DOE OSTI2025

Chemo‐Mechanical Coupling in Hydrogels: Dynamics in the Diffusion‐Limited Regime

Hydrogels are characterized by substantial volume changes in response to external stimuli, making them promising candidates for developing smart materials with enhanced adaptability and responsiveness. By integrating chemical reactions, hydrogels acquire dynamic and tunable responsiveness to external stimuli through chemo‐mechanical coupling, expanding their potential in emerging applications. However, capturing their transient behavior remains challenging due to the complex interplay of chemical reactions, solvent transport, and polymer network deformation. Classical theories capture equilibrium swelling but fail to describe time‐dependent phenomena. To address this, a time‐dependent continuum model is developed that explicitly couples these processes. Volume phase transition in hydrogels with homogeneous chemical reactions is investigated, then the effects of reaction kinetics on these transitions are analyzed. The coupling of these mechanisms is further explored through a study of transient mechanical instabilities. To illustrate the impact of distinct reaction and diffusion timescales, a photo‐active gripper is studied for robotic applications. Finally, a photo‐active microswimmer that exhibits non‐reciprocal motion is proposed to highlight the solvent diffusion for locomotion at the micro‐ and nano‐ scales. The work establishes that transient dynamics of chemo‐mechanical hydrogels generate functions not accessible by steady state models and provides a predictive platform for designing adaptive materials in emerging applications.

chemo-mechanical coupling↗
DOE OSTI2024

An atomistic survey of shear coupling in asymmetric tilt grain boundaries and interpretation using the disconnections framework

Grain Boundaries (GB) play an important role in determining the behavior of polycrystalline materials. While the mechanisms of motion and associated shear response for symmetric tilt grain boundaries (STGBs) are well studied, the same is not true for asymmetric tilt grain boundaries (ATGBs) despite their greater prevalence in polycrystals. Here, this study aims to investigate the shear response of a large collection of asymmetric tilt grain boundaries (ATGBs) using molecular dynamics (MD) simulations and interpret the data using a discrete disconnections model that works remarkably well for STGBs. MD simulations of shear-driven ATGBs show that the plastic shear (shear coupling factor) in the region swept by a GB exhibits a complex dependence on the inclination angle, and this dependence changes with the misorientation of the GB. In addition, the shear response was observed to be highly sensitive to the applied shear rate and temperature. Recognizing the spatial and temporal scale limitations of MD simulations, we extended the discrete disconnections mesoscale model of Khateret al. (2012) to calculate the nucleation barriers of disconnection modes and predict the effective shear coupling of an ATGB. We observed that the mesoscale model’s predictions of the shear coupling factor of ATGBs do not agree with those observed in MD simulations. Finally, we examine the hypotheses of our mesoscale model that contribute to disagreements between MD simulations and the mesoscale model and propose improvements to the mesoscale model for future work.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Generation of straight and curved hollow plasma channels by laser-generated nonlinear wakefields and studies of ultra-intense laser pulse guiding

We report on a study of highly controllable, quasi-static hollow plasma channels generated by ion motion following a laser–plasma interaction via numerical simulations. These channels are generated by ion motion in a nonlinear wakefield, which has an asymmetry between focusing and defocusing periods. Such plasma channels, driven by a fs duration, 100 TW level laser pulse, have the capability to guide more powerful laser pulses. Moreover, we show that curved plasma channels can be generated in a transverse density gradient plasma and alloptical guiding of a second, significantly more powerful laser pulse in such a curved plasma channel. This mechanism may be useful for creating plasma optics for multi-stage TeV laser plasma accelerators and compact synchrotron radiation sources.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Lagrangian conditional statistics and flow topology in edge plasma turbulence

Lagrangian statistics and particle transport in edge plasma turbulence are investigated using the Hasegawa–Wakatani model and its modified version. The latter shows the emergence of pronounced zonal flows. Different values of the adiabaticity parameter are considered. The main goal is to characterize the role of coherent structures, i.e., vortices and zonal flows, and their impact on the Lagrangian statistics of particles. Computationally intensive long time simulations following ensembles of test particles over hundreds of eddy turnover times are considered in statistically stationary turbulent flows. The flow topology is characterized using the Lagrangian Okubo–Weiss criterion in order to split the flow into topologically different domains. In elliptic and hyperbolic regions, the probability density functions (PDFs) of the residence time have self-similar algebraic decaying tails. However, in the intermediate regions, the PDFs exhibit exponentially decaying tails. Topologically conditioned PDFs of the Lagrangian velocity, and acceleration and density fluctuations are likewise computed. The differences between the classical Hasegawa–Wakatani system and its modified version are assessed, and the role of zonal flows is highlighted. Here, the density flux spectrum, which characterizes the contributions of different length scales, is studied, and its inertial scaling is found to be in agreement with predictions based on dimensional arguments. Analyzing the angular change of particle tracers at different time scales, corresponding to coarse grained curvature, completes the study, and these multiscale geometric statistics quantify the directional properties of the particle motion in different flow regimes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Linear analysis of the Kelvin–Helmholtz instability in relativistic magnetized symmetric flows

ABSTRACT We study the linear stability of a planar interface separating two fluids in relative motion, focusing on the symmetric configuration where the two fluids have the same properties (density, temperature, magnetic field strength, and direction). We consider the most general case with arbitrary sound speed cs, Alfvén speed vA, and magnetic field orientation. For the instability associated with the fast mode, we find that the lower bound of unstable shear velocities is set by the requirement that the projection of the velocity on to the fluid-frame wavevector is larger than the projection of the Alfvén speed on to the same direction, i.e. shear should overcome the effect of magnetic tension. In the frame where the two fluids move in opposite directions with equal speed v, the upper bound of unstable velocities corresponds to an effective relativistic Mach number $M_{\rm re}\equiv v/v_{\rm {f}\perp }\sqrt{(1-v_{\rm {f}\perp }^2)/(1-v^2)} \cos \theta =\sqrt{2}$, where $v_{\rm {f}\perp }=[v_{\rm {A}}^2+c_{\rm s}^2(1-v_{\rm {A}}^2)]^{1/2}$ is the fast speed assuming a magnetic field perpendicular to the wavevector (here, all velocities are in units of the speed of light), and θ is the laboratory-frame angle between the flow velocity and the wavevector projection on to the shear interface. Our results have implications for shear flows in the magnetospheres of neutron stars and black holes – both for single objects and for merging binaries – where the Alfvén speed may approach the speed of light.

Astronomy & Astrophysics↗