Search NASA⌕ Search

SEARCH · Search NASA

Results for “Lorentz force”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Friction in a strongly magnetized neutral plasma

A charged projectile generates an electrostatic wake as it moves through a background plasma. This wake normally gives rise to a drag force that acts antiparallel to the velocity vector of the projectile, and is commonly referred to as the stopping power. A recent theory predicted that an additional component of the friction force that is perpendicular to both the projectile velocity and the Lorentz force arises when a plasma is strongly magnetized. We extend the previous analysis of this transverse friction force, which was based on the one-component plasma model, to treat a two-component charge-neutral electron-ion plasma. The direction and magnitude of the transverse force is found to depend on the projectile speed relative to the electron and ion thermal velocities, and to change sign three times. For projectile speeds below a certain threshold, the transverse force is dominated by collisions with ions, while above this threshold, only electron collisions contribute. Here, the average trajectory of the projectile is significantly altered by the transverse force as it slows.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Efficient calculation of self magnetic field, self-force, and self-inductance for electromagnetic coils with rectangular cross-section

Abstract For designing high-field electromagnets, the Lorentz force on coils needs to be computed in order to design suitable support structures, and the inductance should be computed to evaluate the stored energy and dynamics. Also, the magnetic field and its variation inside the conductor is of interest for computing stress and strain, and due to superconducting quench limits. For these force, inductance, energy, and internal field calculations, the coils cannot be naively approximated as infinitesimally thin filaments due to divergences when the source and evaluation points coincide, so more computationally demanding calculations are usually required, resolving the finite cross-section of the conductors. Here, we present a new alternative method that enables the internal magnetic field vector, self-force, and self-inductance to be computed rapidly and accurately within a 1D filament model. The method is applicable to coils for which the curve center-line can have general noncircular shape, as long as the conductor width is small compared to the radius of curvature. This paper extends a previous calculation for circular-cross-section conductors (Hurwitz et al 2024 IEEE Trans. Magn. ) to consider the case of rectangular cross-section. The reduced model is derived by rigorous analysis of the singularity, regularizing the filament integrals such that they match the true high-dimensional integrals at high coil aspect ratio. The new filament model exactly recovers analytic results for a circular coil, and is shown to accurately reproduce full finite-cross-section calculations for a non-planar coil of a stellarator magnetic fusion device. Due to the efficiency of the model here, it is well suited for use inside design optimization.

Landreman, Matt (ORCID:000000027233577X)↗
DOE OSTI2025

Study of ponderomotive instability in the FRIB beta=0.53 half-wave resonator

Superconducting radio-frequency niobium cavities are susceptible to deformations caused by external or internal forces, leading to shifts in the cavity resonant frequency. One source of deformation comes from the radiation pressure of the cavity fields, producing the so-called Lorentz force detuning effect. This effect can couple to the cavity mechanical modes in generator driven mode, leading to ponderomotive instability. In the FRIB 322 MHz, β=0.53 Half-wave Resonators (HWR), the instability appeared when the cavity was detuned, with thresholds depending on low-level RF control parameters, such as closed loop gain, as well as accelerating gradient. Using a measured Lorentz Transfer Function, Simulink simulations were conducted to predict the instability thresholds, which were then compared with the experimental results. We will discuss how these thresholds can be broadened to enable stable operation at higher gradients.

Accelerator Physics↗
DOE OSTI2020

A generalized Boltzmann kinetic theory for strongly magnetized plasmas with application to friction

Coulomb collisions in plasmas are typically modeled using the Boltzmann collision operator, or its variants, which apply to weakly magnetized plasmas in which the typical gyroradius of particles significantly exceeds the Debye length. Conversely, O’Neil has developed a kinetic theory to treat plasmas that are so strongly magnetized that the typical gyroradius of particles is much smaller than the distance of the closest approach in a binary collision. Here, we develop a generalized collision operator that applies across the full range of magnetization strength. Since there is no closed-form solution for the scattering cross section when plasma is strongly magnetized, the input to the collision operator is obtained by numerically calculating particle trajectories. To demonstrate novel physics associated with strong magnetization, it is used to compute the friction force on a massive test charge. In addition to the traditional stopping power component, this is found to exhibit a transverse component that is perpendicular to both the velocity and Lorentz force vectors in the strongly magnetized regime, as was predicted recently using linear response theory. As a result, good agreement is found between the collision theory and linear response theory in the regime in which both apply, but the new collision theory also applies to stronger magnetization strength regimes than the linear response theory is expected to apply in.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Friction force in strongly magnetized plasmas

A charged particle moving through a plasma experiences a friction force that commonly acts antiparallel to its velocity. It was recently predicted that in strongly magnetized plasmas, in which the plasma particle gyrofrequency exceeds the plasma frequency, the friction also includes a transverse component that is perpendicular to both the velocity and Lorentz force. Here, this prediction is confirmed using molecular-dynamics simulations, and it is shown that the relative magnitude of the transverse component increases with plasma coupling strength. Furthermore, this result influences single-particle motion and macroscopic transport in strongly magnetized plasmas found in a broad range of applications.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Global stellarator coil optimization with quadratic constraints and objectives

Most present stellarator designs are produced by costly two-stage optimization: the first for an optimized equilibrium, and the second for a coil design reproducing its magnetic configuration. Few proxies for coil complexity and forces exist at the equilibrium stage. Rapid initial state finding for both stages is a topic of active research. Most present convex coil optimization codes use the least square winding surface method by Merkel (NESCOIL), with recent improvements in conditioning, regularization, sparsity, and physics objectives. While elegant, the method is limited to modeling the norms of linear functions in coil current. We present QUADCOIL, a global coil optimization method that targets combinations of linear and quadratic functions of the current. It can directly constrain and/or minimize a wide range of physics objectives unavailable in NESCOIL and REGCOIL, including the Lorentz force, magnetic energy, curvature, field-current alignment, and the maximum density of a dipole array. QUADCOIL requires no initial guess and runs nearly $10$ 2 x faster than filament optimization. Integrating it in the equilibrium optimization stage can potentially exclude equilibria with difficult-to-design coils, without significantly increasing the computation time per iteration. QUADCOIL finds the exact, global minimum in a large parameter space when possible, and otherwise finds a well-performing approximate global minimum. It supports most regularization techniques developed for NESCOIL and REGCOIL. We demonstrate QUADCOIL’s effectiveness in coil topology control, minimizing non-convex penalties, and predicting filament coil complexity with three numerical examples.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Coexistence of Merons with Skyrmions in the Centrosymmetric Van Der Waals Ferromagnet Fe 5– x GeTe 2

Fe 5-x GeTe 2 is a centrosymmetric, layered van der Waals (vdW) ferromagnet that displays Curie temperatures T c (270–330 K) that are within the useful range for spintronic applications. However, little is known about the interplay between its topological spin textures (e.g., merons, skyrmions) with technologically relevant transport properties such as the topological Hall effect (THE) or topological thermal transport. Here, via high-resolution Lorentz transmission electron microscopy, it is shown that merons and anti-meron pairs coexist with Néel skyrmions in Fe 5-x GeTe 2 over a wide range of temperatures and probe their effects on thermal and electrical transport. A THE is detected, even at room T, that senses merons at higher T's, as well as their coexistence with skyrmions as T is lowered, indicating an on-demand thermally driven formation of either type of spin texture. Remarkably, an unconventional THE is also observed in absence of Lorentz force, and it is attributed to the interaction between charge carriers and magnetic field-induced chiral spin textures. Here these results expose Fe 5-x GeTe 2 as a promising candidate for the development of applications in skyrmionics/meronics due to the interplay between distinct but coexisting topological magnetic textures and unconventional transport of charge/heat carriers.

two-dimensional magnets↗
DOE OSTI2023

Synergizing a Large Ordinary Nernst Effect and Axis‐Dependent Conduction Polarity in Flat Band KMgBi Crystals

The exploration of quantum materials in which an applied thermo/electrical/magnetic field along one crystallographic direction produces an anisotropic response has led to unique functionalities. Along these lines, KMgBi is a layered, narrow gap semiconductor near a critical state between multiple Dirac phases due to the presence of a flat band near the Fermi level. The valence band is highly anisotropic with minimal cross-plane dispersion, which, in combination with an isotropic conduction band, enables axis-dependent conduction polarity. Thermopower and Hall measurements indicate dominant p-type conduction along the cross-plane direction, and n-type conduction along the in-plane direction, leading to a significant zero-field transverse thermoelectric response when the heat flux is at an angle to the principal crystallographic directions. Additionally, a large Ordinary Nernst effect (ONE) is observed with an applied field. It arises from the ambipolar term in the Nernst effect, whereby the Lorentz force on electrons and holes makes them drift in opposite directions so that the resulting Nernst voltage becomes a function of the difference between their partial thermopowers, greatly enhancing the ONE. It is proven that axis-dependent polarity can synergistically enhance the ONE, in addition to leading to a zero-field transverse thermoelectric performance.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Interactions of particles with “continuous spin” fields

Powerful general arguments allow only a few families of long-range interactions, exemplified by gauge field theories of electromagnetism and gravity. However, all of these arguments presuppose that massless fields have zero spin scale (Casimir invariant) and hence exactly boost invariant helicity. This misses the most general behavior compatible with Lorentz symmetry. We present a Lagrangian formalism describing interactions of matter particles with bosonic “continuous spin” fields with arbitrary spin scale ρ. Remarkably, physical observables are well approximated by familiar theories at frequencies larger than ρ, with calculable deviations at low frequencies and long distances. For example, we predict specific ρ-dependent modifications to the Lorentz force law and the Larmor formula, which lay the foundation for experimental tests of the photon’s spin scale. We also reproduce known soft radiation emission amplitudes for nonzero ρ. The particles’ effective matter currents are not fully localized to their worldlines when ρ ≠ 0, which motivates investigation of manifestly local completions of our theory. Our results also motivate the development of continuous spin analogues of gravity and non-Abelian gauge theories. Given the correspondence with familiar gauge theory in the small ρ limit, we conjecture that continuous spin particles may in fact mediate known long-range forces, with testable consequences.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Twist analysis of the spin-orbit correlation in QCD

We present a QCD analysis of the twist-three parton distribution functions associated with the spin-orbit correlation of quarks and gluons in spin-$\frac{1}{2}$ and spin-0 hadrons. We derive exact non-perturbative identities decomposing the spin-orbit correlations into the Wandzura-Wilczek part and the genuine twist-three part. In the spin-$\frac{1}{2}$ case, the result is partially related to the kinematical twist-three part of the g T (x) distribution familiar in the context of transverse spin physics. We use these identities to obtain a novel longitudinal momentum sum rule which may be regarded as the momentum version of the Jaffe-Manohar spin sum rule. We explore the physical interpretation of the sum rule and make a connection to the color Lorentz forces and their associated potential energies.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Free-Surface Liquid Lithium Flow Modeling and Stability Analysis for Fusion Applications

Liquid metal plasma facing components are considered an attractive design choice for fusion devices including pilot plants. Virtual prototyping of such devices includes modeling of free-surface flow of the electrically conductive liquid, which requires computational fluid dynamics (CFD) and magnetohydrodynamics (MHD) simulations. Numerical tools capable of simulating flows and heat transfer in the free-surface MHD flow were developed at PPPL based on the customized ANSYS CFX. Here, MHD is introduced using a magnetic vector potential approach. Free-surface flow capabilities are available in the code and were tested. Special stabilization procedures were derived and applied to improve convergence of the momentum equations with the source terms due to the Lorentz force and surface tension. Important characteristics of the fusion-relevant liquid metal flow is free surface smoothness and stability. Heat flux from the plasma impacts the liquid surface at a very acute angle, so any change of the free surface from axisymmetry can dramatically increase the local heat flux density and thus create excessive evaporation of liquid lithium into the plasma, which is detrimental to operations. Stability analysis of the liquid metal film flow was performed to determine applicable flow regimes. Thin film flow along horizontal wall is considered. Effects of gravity, magnetic field, and surface tension are included in the analysis.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Effects of an External Magnetic Field on Keyhole Mode Laser Melting of 316 Stainless Steel

Keyhole-mode laser melting is an efficient method for joining or cutting large, thick components, but controlling keyhole depth and fluctuations has remained challenging. Applying an external magnetic field can control melt pool flows and indirectly influence keyhole morphology and dynamics. The induced Lorentz force, comprising Seebeck and damping components, plays a crucial role in the melt pool dynamics, depending on temperature gradient, flow rates, and magnetic field orientation and magnitude. This research investigates the effects of an external magnetic field on keyhole behavior during laser spot melting of 316 stainless steel using synchronized high-speed synchrotron X-ray and thermal imaging. Findings revealed that a longitudinal magnetic field (120 mT) increased keyhole depth but exacerbated lateral fluctuations, resulted in a 20% increase in the melt pool temperature gradient and a 27% decrease in cooling rate. Conversely, a transverse magnetic field (760 mT) reduced keyhole depth and improved porosity formation. The findings suggest that a decrease in keyhole depth correlates with a decrease in fluctuations, and vice versa. These insights enhance understanding of external magnetic fields’ impact on laser melting, with implications for improving part quality.

Alamdari, Aslan Bafahm [The Ohio State Univ., Colu↗
DOE OSTI2023

Quantitative proton radiography and shadowgraphy for arbitrary intensities

Charged-particle radiography and shadowgraphy data can be directly inverted to obtain a line-integrated transverse Lorentz force or a line-integrated transverse refractive index gradient if intensity modulations due to scattering and absorption are negligible, and angular deflections are small. We develop a new direct-inversion algorithm based on plasma physics and compare it to a new Monge–Ampère code and an existing power diagram code. The measured or source intensity is represented by electrons subject to drag, and the other intensity by fixed ions. The decrease in kinetic plus electrostatic energy determines convergence. The displacement of the electrons from their initial to their equilibrium positions determines the line-integrated force or refractive index gradient. We have implemented two approaches: PIC (particle in cell) and Lagrangian fluid, in 1-D and 2-D. The PIC code works for arbitrary intensities, can work efficiently in parallel, and can make use of existing codes. The Lagrangian code requires less memory and is faster than the PIC code without massively parallel processing, but fails in 2-D for large intensity modulations. The Monge–Ampère code is by far the fastest in 2-D, without massively parallel processing, but fails for intensities with large voids, high contrast ratios and large deflections across the boundaries, and could not obtain the degree of convergence possible with the PIC code. As a result, the power diagram code was by far the slowest and most memory intensive, and failed for large peaks in the measured intensity.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Accurately simulating nine-dimensional phase space of relativistic particles in strong fields

Next-generation high-power laser systems that can be focused to ultra-high intensities exceeding 10 23 W/cm 2 are enabling new physics regimes and applications. The physics of how these lasers interact with matter is highly nonlinear, relativistic, and can involve lowest-order quantum effects. The current tool of choice for modeling these interactions is the particle-in-cell (PIC) method. In the presence of strong electromagnetic fields, the motion of charged particles and their spin is affected by radiation reaction (either the semi-classical or the quantum limit). Standard (PIC) codes usually use Boris or similar operator-splitting methods to advance the particles in standard phase space. These methods have been shown to require very small time steps in the strong-field regime in order to obtain accurate results. In addition, some problems require tracking the spin of particles, which creates a nine-dimensional (9D) particle phase space, i.e., (x, u, s). Therefore, numerical algorithms that enable high-fidelity modeling of the 9D phase space in the strong-field regime (where both the spin and momentum evolution are affected by radiation reaction) are desired. Here we present a new particle pusher that works in 9D and 6D phase space (i.e., with and without spin) based on analytical rather than leapfrog solutions to the momentum and spin advance from the Lorentz force, together with the semi-classical form of radiation reaction in the Landau-Lifshitz equation and spin evolution given by the Bargmann-Michel-Telegdi equation. Analytical solutions for the position advance are also obtained, but these are not amenable to the staggering of space and time in standard PIC codes. These analytical solutions are obtained by assuming a locally uniform and constant electromagnetic field during a time step. The solutions provide the 9D phase space advance in terms of a particle's proper time, and a mapping is used to determine the proper time step duration for each particle as a function of the lab frame time step. Due to the analytical integration of particle trajectory and spin orbit, the constraint on the time step needed to resolve trajectories in ultra-high fields can be greatly reduced. The time step required in a PIC code for accurately advancing the fields may provide additional constraints. We present single-particle simulations to show that the proposed particle pusher can greatly improve the accuracy of particle trajectories in 6D or 9D phase space for given laser fields. We have implemented the new pusher into the PIC code Osiris. Example simulations show that the proposed pusher provides improvement for a given time step. A discussion on the numerical efficiency of the proposed pusher is also provided.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Intense laser interaction with micro-bars

Intense laser fields interact very differently with micrometric rough surfaces than with flat objects. The interaction features high laser energy absorption and increased emission of MeV electrons, ions, and of hard x-rays. In this work, we irradiated isolated, translationally-symmetric objects in the form of micrometric Au bars. The interaction resulted in the emission of two forward-directed electron jets having a small opening angle, a narrow energy spread in the MeV range, and a positive angle to energy correlation. Our numerical simulations show that following ionization, those electrons that are pulled into vacuum near the object’s edge, remain in-phase with the laser pulse for long enough so that the Lorentz force they experience drive them around the object’s edge. After these electrons pass the object, they form attosecond duration bunches and interact with the laser field over large distances in vacuum in confined volumes that trap and accelerate them within a narrow range of momentum. The selectivity in energy of the interaction, its directionality, and the preservation of the attosecond duration of the electron bunches over large distances, offer new means for designing future laser-based light sources.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Plasma facing components with capillary porous system and liquid metal coolant flow

Liquid metal can create a renewable protective surface on plasma facing components (PFC), with an additional advantage of deuterium pumping and the prospect of tritium extraction if liquid lithium (LL) is used and maintained below 450 °C, the temperature above which LL vapor pressure begins to contaminate the plasma. LM can also be utilized as an efficient coolant, driven by the Lorentz force created with the help of the magnetic field in fusion devices. Capillary porous systems can serve as a conduit of LM and simultaneously provide stabilization of the LM flow, protecting against spills into the plasma. Recently, a combination of a fast-flowing LM cooling system with a porous plasma facing wall (CPSF) was investigated [A. Khodak and R. Maingi, Nucl. Mater. Energy 26, 100935 (2021)]. The system takes an advantage of a magnetohydrodynamics velocity profile as well as attractive LM properties to promote efficient heat transfer from the plasma to the LL at low pumping energy cost, relative to the incident heat flux on the PFC. In the case of a disruption leading to excessive heat flux from the plasma to the LM PFCs, LL evaporation can stabilize the PFC surface temperature, due to high evaporation heat and apparent vapor shielding. The proposed CPSF was optimized analytically for the conditions of a fusion nuclear science facility [Kessel et al., Fusion Sci. Technol. 75, 886 (2019)]: 10 T toroidal field and 10 MW/m 2 peak incident heat flux. Finally, computational fluid dynamics analysis confirmed that a CPSF system with 2.5 mm square channels can pump enough LL so that no additional coolant is needed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Contactless excitation of acoustic resonance in insulating wafers

Contactless excitation and detection of high harmonic acoustic overtones in a thin insulator single crystal are described using radio frequency spectroscopy techniques. Single crystal [001] silicon wafer samples were investigated, one side covered with a Nb thin film, the common starting point for the fabrication of quantum devices. Additionally, the coupling between electromagnetic signals and mechanical oscillation is achieved from the Lorentz force generated by an external magnetic field. This method is suitable for any sample with a metallic surface or covered with a thin metal film. High resolution measurements of the temperature dependence of the sound velocity and elastic constants of silicon are reported and compared with known results.

47 OTHER INSTRUMENTATION↗
DOE OSTI2025

On finite-dimensional smoothed-particle Hamiltonian reductions of the Vlasov equation

The inclusion of spatial smoothing in finite-dimensional particle-based Hamiltonian reductions of the Vlasov equation and related models is considered. Here, this work investigates the underlying Hamiltonian structure of such smoothed particle-based methods for Hamiltonian systems and the small-scale regularization such methods implicitly make in approximating the continuum theory. In the context of the Vlasov–Poisson equation and other mean-field Lie–Poisson systems, of which Vlasov–Poisson is a special case, smoothing amounts to a convolutive regularization of the Hamiltonian. This regularization may be interpreted as a change of the inner product structure used to identify the dual space in the Lie–Poisson Hamiltonian formulation. In particular, the shape function used for spatial smoothing may be identified as the kernel function of a reproducing kernel Hilbert space whose inner product is used to define the Lie–Poisson Hamiltonian structure. It is likewise possible to introduce smoothing in the Vlasov–Maxwell system, but in this case the Poisson bracket must be modified rather than the Hamiltonian. The smoothing applied to the Vlasov–Maxwell system is incorporated by inserting smoothing in the map from canonical to kinematic coordinates. In the filtered system, the Lorentz force law and the current, the two terms coupling the Vlasov equation with Maxwell’s equations, are spatially smoothed.

Hamiltonian mechanics↗