Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “fields”

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

Role of perturbed parallel magnetic field effects in predicting turbulent transport in NSTX

This study presents analysis of gyrokinetic simulations on the National Spherical Torus Experiment (NSTX) to investigate the effects of electromagnetic fields on plasma turbulence and transport. The simulations, performed with varying levels of fidelity using the gyrokinetic CGYRO code, include electrostatic (ES), single-field electromagnetic (EM1), and two-field electromagnetic (EM2) models. A detailed comparison across the simulation database reveals that electromagnetic effects increase both predicted growth rates and quasilinear fluxes, with EM2 simulations producing stronger turbulence than ES and EM1 cases. Quasilinear modeling using QLGYRO demonstrates that while the perturbed parallel magnetic field (δB ∥ ) does not drastically affect the total flux at experimental gradients, it leads to a shift in the dominant instability, altering mode structures from microtearing to kinetic ballooning modes (KBMs). The proximity of the plasma profiles to the KBM threshold is explored, with the experimental conditions being near the onset of KBM-driven transport. The KBM, with its large growth rates, is identified as a potential driver of electron temperature flattening, as it can rapidly transport heat across flux surfaces. Performing stability analysis shows core-localized unstable a low- mode that could contribute to the flattening at the early times of the discharge. TGYRO predictive modeling, incorporating both TGLF and QLGYRO, indicates that the inclusion of δB ∥ significantly improves the accuracy of temperature profile predictions in NSTX high-beta plasmas, although challenges remain in modeling the sharp flux discontinuities caused by KBM-driven instabilities.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Error field identification through torque balance on a saturated island in DIII-D

Measurement of the electromagnetic torque on a magnetic island could be an attractive method for error field identification in the early phase of ITER operation. Previous DIII-D experiments (Strait 2014 Nucl. Fusion 54 073004; Shiraki et al 2015 Plasma Phys. Control. Fusion 57 025016) have demonstrated the principle of this approach using a stationary or slowly rotating island, while recent developments in magnetic data analysis (Sweeney and Strait 2019 Phys. Plasmas 26 012509) allow the field of a rapidly rotating island to be readily distinguished from that of the wall currents induced by its rotation. In a recent experiment, a rotating n = 1 magnetic perturbation forced a saturated magnetic island to rotate, thus sampling all toroidal phases periodically in a single discharge. The phase and amplitude of the error field are inferred from analysis of the time-dependent torque balance on the island, including torques from the error field, the applied magnetic perturbation, and the wall currents induced by rotation of the applied perturbation and the island. Furthermore, results agree well with those from more conventional methods.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Phase-controlled direct laser acceleration enabled by longitudinal variation of the laser-driven quasi-static plasma magnetic field

Direct laser acceleration (DLA) enables energy transfer from an ultra-high-intensity laser to plasma electrons and underpins many laser-driven particle and radiation-source concepts. A laser-driven azimuthal plasma magnetic field is a key player in this process: it confines energetic electrons, induces betatron oscillations, and makes possible a resonant interaction between the betatron motion and the laser field. While this betatron resonance can enhance electron energy gain, the gain itself generally drives frequency detuning and promotes largely reversible energy exchange that limits net acceleration. Here we show, using a test-electron model with prescribed fields, that a slow longitudinal increase of the quasi-static plasma magnetic field qualitatively changes DLA by introducing hysteresis in the ratio of the betatron frequency to the laser frequency experienced by the electron, so that this ratio depends on the prior evolution of the electron even at the same energy. This hysteresis enables phase control of the electron-laser energy exchange and suppresses the usual reversibility of DLA, allowing electrons to retain the acquired energy and sustain energy gain without intermittent losses.

direct laser acceleration↗
DOE OSTI2023

MHD simulation on magnetic compression of field reversed configurations with NIMROD

Abstract Magnetic compression has long been proposed a promising method for plasma heating in a field reversed configuration (FRC). However, it remains a challenge to fully understand the physical mechanisms underlying the compression process, due to its highly dynamic nature beyond the one-dimensional (1D) adiabatic theory model (Spencer et al 1983 Phys. Fluids 26 1564). In this work, magnetohydrodynamics simulations on the magnetic compression of FRCs using the NIMROD code (Sovinec et al 2004 J. Comput. Phys. 195 355) and their comparisons with the 1D theory have been performed. The effects of the assumptions of the theory on the compression process have been explored, and the detailed profiles of the FRC during compression have been investigated. The pressure evolution agrees with the theoretical prediction under various initial conditions. The axial contraction of the FRC can be affected by the initial density profile and the ramping rate of the compression magnetic field, but the theoretical predictions on the FRC’s length in general and the relation r s = 2 r o in particular hold approximately well during the whole compression process, where r s is the major radius of FRC separatrix and r o is that of the magnetic axis. The evolutions of the density and temperature can be affected significantly by the initial equilibrium profile and the ramping rate of the compression magnetic field. During the compression, the major radius of the FRC is another parameter that is susceptible to the ramping rate of the compression field. Basically, for the same magnetic compression ratio, the peak density is higher and the FRC’s radius r s is smaller than the theoretical predictions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

FLARE: field line analysis and reconstruction for 3D boundary plasma modeling

The FLARE code is a magnetic mesh generator that is integrated within a suite of tools for the analysis of the magnetic geometry in toroidal fusion devices. A magnetic mesh is constructed from field line segments and permits fast reconstruction of field lines in 3D boundary plasma codes such as EMC3-EIRENE. Both intrinsically non-axisymmetric configurations (stellarators) and those with symmetry breaking perturbations of an axisymmetric equilibrium (tokamaks) are supported. The code itself is written in Modern Fortran with MPI support for parallel computing, and it incorporates object-oriented programming for the definition of the magnetic field and the material surface geometry. Extended derived types for a number of different magnetohydrodynamic equilibrium and plasma response models are implemented. The core element of FLARE is a field line tracer with adaptive step-size control, and this is integrated into tools for the construction of Poincaré maps and invariant manifolds of X-points. A collection of high-level procedures that generate output files for visualization is build on top of that. The analysis modules are build with Python frontends that facilitate customization of tasks and/or scripting of parameter scans.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

A hybrid 3D/2D field response calculation for liquid argon detectors with PCB based anode plane

Liquid Argon Time Projection Chamber (LArTPC) technology is commonly utilized in neutrino detector designs. It enables detailed reconstruction of neutrino events with high spatial precision and low energy threshold. Its field response (FR) model describes the time-dependent electric currents induced in the anode-plane electrodes when ionization electrons drift nearby. An accurate and precise FR is a crucial input to LArTPC detector simulations and charge reconstruction. Established LArTPC designs have been based on parallel wire planes. It allows accurate and computationally economic two-dimensional (2D) FR models utilizing the translational symmetry along the direction of the wires. Recently, novel LArTPC designs utilize electrodes formed on printed circuit board (PCB) in the shape of strips with through holes. The translational symmetry is no longer a good approximation near the electrodes and a new FR calculation that employs regions with three dimensions (3D) has been developed. Extending the 2D models to 3D would be computationally expensive. Fortuitously, the nature of strips with through holes allows for a computationally economic approach based on the finite-difference method (FDM). In this paper, we present a new software package pochoir that calculates LArTPC field response for these new strip-based anode designs. This package combines 3D calculations in the volume near the electrodes with 2D far-field solutions to achieve fast and precise field response computation. We apply the resulting FR to simulate and reconstruct samples of cosmic-ray muons and 39 Ar decays from a Vertical Drift (VD) detector prototype operated at CERN. We find the difference between real and simulated data within 5%. Current state-of-the-art LArTPC software requires a 2D FR which we provide by averaging over one dimension and estimate that variations lost in this average are smaller than 7%.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Cosmic ray transport in mixed magnetic fields and their role on the observed anisotropies

ABSTRACT There is a growing set of observational data demonstrating that cosmic rays exhibit small-scale anisotropies (5°–30°) with amplitude deviations lying between 0.01–0.1 per cent that of the average cosmic ray flux. A broad range of models have been proposed to explain these anisotropies ranging from finite-scale magnetic field structures to dark matter annihilation. The standard diffusion transport methods used in cosmic ray propagation do not capture the transport physics in a medium with finite-scale or coherent magnetic field structures. Here, we present a Monte Carlo transport method, applying it to a series of finite-scale magnetic field structures to determine the requirements of such fields in explaining the observed cosmic ray, small-scale anisotropies.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Nonadiabatic Strong-Field Photoionization Revisited

We measure strong field ionization of cesium atoms, observing a robust feature near 2⁢𝑈p in the photoelectron spectrum, which we call the intermediate energy structure (IES). Using a Coulomb-corrected strong-field approximation, we show it arises from electrons born with large inward velocities that rapidly undergo forward scattering off the Coulomb potential. The IES is similar to the previously identified low energy structure in that they are both due to forward scattering. However, it is different in that IES requires both nonadiabatic ionization conditions (as defined by the Keldysh parameter) and a more weakly bound initial state. Furthermore, our joint experimental and theoretical study supports the need for an additional parameter in characterization of strong field ionization, which depends both on nonadiabaticity of the process and the principal quantum number of the target. Together with Reiss and Keldysh parameters, this dimensionless parameter expands the description of strong field ionization across a broad range of atomic targets and incorporates ionization from excited states.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

Exact solutions for the electromagnetic fields of a flying focus

The intensity peak of a “flying” focus travels at a programmable velocity over many Rayleigh ranges while maintaining a near-constant profile. Assessing the extent to which these features can enhance laser-based applications requires an accurate description of the electromagnetic fields. Here we present exact analytical solutions to Maxwell’s equations for the electromagnetic fields of a constant-velocity flying focus, generalized for arbitrary polarization and orbital angular momentum. Furthermore, the approach combines the complex source-point method, which transforms multipole solutions into beam-like solutions, with the Lorentz invariance of Maxwell’s equations. Propagating the fields backward in space reveals the space-time profile that an optical assembly must produce to realize these fields in the laboratory. Comparisons with simpler paraxial solutions provide conditions for their reliable use when modeling a flying focus.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Strategies for simulating the time evolution of Hamiltonian lattice field theories

Simulating the time evolution of quantum field theories given some Hamiltonian H requires developing algorithms for implementing the unitary operator e -iHt . A variety of techniques exist that accomplish this task, with the most common technique used so far being Trotterization, which is a special case of the application of a product formula. However, other techniques exist that promise better asymptotic scaling in certain parameters of the theory being simulated, the most efficient of which are based on the concept of block encoding. In this work we study the performance of such algorithms in simulating lattice field theories. We derive and compare the asymptotic gate complexities of several commonly used simulation techniques in application to Hamiltonian lattice field theories. Using the scalar $\hat{φ}$ 4 theory as a test, we also perform numerical studies and compare the gate costs required by product formulas and signal-processing-based techniques to simulate time evolution. For the latter, we use the linear combination of unitaries (LCU) construction augmented with the quantum Fourier transform circuit to switch between the field and momentum eigenbases, which leads to immediate order-of-magnitude improvement in the cost of preparing the block encoding. Further, this paper also includes a pedagogical review of the techniques used, in particular product formulas, LCU, qubitization, quantum signal processing, as well as the technique for simulating geometrically-local Hamiltonians developed by Haah, Hastings, Kothari, and Low.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Theory and phase-field simulations of electrical control of spin cycloids in a multiferroic

Cycloidal spin orders are common in multiferroics. One of the prototypical examples is BiFeO 3 (BFO) which shows a large polarization and a cycloidal antiferromagnetic (AFM) order at room temperature. Here, we employ Landau theory and phase-field simulations to analyze the coupled switching dynamics of polarization and cycloidal AFM orders in BFO. We identify 14 types of transitional spin structures between two cycloids and nine electric-field-induced spin switching paths. We demonstrate the electric-field-induced rotation of wave vectors of the cycloidal spins and discover two types of cycloidal spin switching dynamics: fast local spin flips and slow rotation of wave vectors. Also, we construct roadmaps to achieve the switching between any two spin cycloids through multistep applications of electric fields. Furthermore, this paper provides a theoretical framework for the phenomenological description of spin cycloids and a fundamental understanding of the switching mechanisms to achieve electrical control of magnetic orders.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Protected fermionic zero modes in periodic gauge fields

It is well known that macroscopically normalizable zero-energy wave functions of spin-$\frac{1}{2}$ particles in a two-dimensional inhomogeneous magnetic field are spin-polarized and exactly calculable with degeneracy equaling the number of flux quanta linking the whole system. Here, extending this argument to massless Dirac fermions subjected to magnetic fields that have zero net flux but are doubly periodic in real space, we show that there exist only two Bloch-normalizable zero-energy eigenstates, one for each spin flavor. This result is immediately relevant to graphene multilayer systems subjected to doubly periodic strain fields, which at low energies enter the Hamiltonian as periodic pseudogauge vector potentials. Furthermore, we explore various related settings including nonlinearly dispersing band structure models and systems with singly periodic magnetic fields.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Constructing chiral effective field theory around unnatural leading-order interactions

A momentum-dependent formulation based on a stationary spin-0 and isospin-1 dibaryon field is proposed to improve convergence of chiral effective field theory in the 1 S 0 channel of NN scattering. Although the two-parameter leading-order interaction appears to be unnatural, it nevertheless has the necessary features of an effective field theory. A rapid order-by-order convergence is found in 1 S 0 . As an application beyond the two-body level, the triton binding energy is studied and compared to standard chiral effective field theory with partly perturbative pions. Here, the consistency of the chiral Lagrangian for the new formulation is examined by working out the pionic radiative corrections, and consequences of nontrivial chiral-connection terms are discussed.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Deformations, relaxation, and broken symmetries in liquids, solids, and glasses: A unified topological field theory

Here, we combine hydrodynamic and field theoretic methods to develop a general theory of phonons as Goldstone bosons in crystals, glasses, and liquids based on nonaffine displacements and the consequent Goldstone phase relaxation. We relate the conservation, or lack thereof, of specific higher-form currents with properties of the underlying deformation field—nonaffinity—which dictates how molecules move under an applied stress or deformation. In particular, the single-valuedness of the deformation field is associated with conservation of higher-form charges that count the number of topological defects. Our formalism predicts, from first principles, the presence of propagating shear waves above a critical wave vector in liquids, thus giving a formal derivation of the phenomenon in terms of fundamental symmetries. The same picture provides also a theoretical explanation of the corresponding “positive sound dispersion” phenomenon for longitudinal sound. Importantly, accordingly to our theory, the main collective relaxation timescale of a liquid or a glass (known as the α relaxation for the latter) is given by the phase relaxation time, which is not necessarily related to the Maxwell time. Finally, we build a nonequilibrium effective action using the in-in formalism defined on the Schwinger-Keldysh contour, that further supports the emerging picture. In summary, our work suggests that the fundamental difference between solids, fluids, and glasses has to be identified with the associated generalized higher-form global symmetries and their topological structure, and that the Burgers vector for the displacement fields serves as a suitable topological order parameter distinguishing the solid (ordered) phase and the amorphous ones (fluids, glasses).

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
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 OSTI2023

Nonlinear ablative Rayleigh-Taylor instability: Increased growth due to self-generated magnetic fields

The growth rate of the nonlinear ablative Rayleigh-Taylor (RT) instability is enhanced by magnetic fields self-generated by the Biermann battery mechanism; a scaling for this effect with perturbation height and wavelength is proposed and validated with extended-magnetohydrodynamic simulations. The magnetic flux generation rate around a single RT spike is found to scale with the spike height. The Hall parameter, which quantifies electron magnetization, is found to be strongly enhanced for short-wavelength spikes due to Nernst compression of the magnetic field at the spike tip. Furthermore, the impact of the magnetic field on spike growth is through both the suppressed thermal conduction into the unstable spike and the Righi-Leduc heat flow deflecting heat from the spike tip to the base. Righi-Leduc is found to be the dominant effect for small Hall parameters, while suppressed thermal conduction dominates for large Hall parameters. These results demonstrate the importance of considering magnetic fields in all perturbed inertial confinement fusion hot spots.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Magnetoelectric effect arising from a field-induced pseudo Jahn-Teller distortion in a rare-earth magnet

Magnetoelectric materials are attractive for several applications, including actuators, switches, and magnetic field sensors. Typical mechanisms for achieving a strong magnetoelectric coupling are rooted in transition metal magnetism. In sharp contrast, here we identify CsEr(MoO 4 ) 2 as a magnetoelectric material without magnetic transition metal ions, thus ensuring that the Er ions play a key role in achieving this interesting property. Our detailed study includes measurements of the structural, magnetic, and magnetoelectric properties of this material. Bulk characterization and neutron powder diffraction show no evidence for structural phase transitions down to 0.3 K and therefore CsEr(MoO 4 ) 2 maintains the room temperature P2/c space group over a wide temperature range without external magnetic field. These same measurements also identify collinear antiferromagnetic ordering of the Er 3+ moments below T N =0.87K. Complementary dielectric constant and pyroelectric current measurements reveal that a ferroelectric phase with a maximum polarization P~0.6 nC/cm 2 emerges when applying a modest external magnetic field, which indicates that this material has a strong magnetoelectric coupling. Finally, we argue that the magnetoelectric coupling in this system arises from a pseudo Jahn-Teller distortion induced by the magnetic field.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Towards Feature-Preserving 2D and 3D Vector Field Compression

The objective of this work is to develop error-bounded lossy compression methods to preserve topological features in 2D and 3D vector fields. Specifically, we explore the preservation of critical points in piecewise linear vector fields. We define the preservation of critical points as, without any false positive, false negative, or false type change in the decompressed data, (1) keeping each critical point in its original cell and (2) retaining the type of each critical point (e.g., saddle and attracting node). The key to our method is to adapt a vertex-wise error bound for each grid point and to compress input data together with the error bound field using a modified lossy compressor. Our compression algorithm can be also embarrassingly parallelized for large data handling and in situ processing. We benchmark our method by comparing it with existing lossy compressors in terms of false positive/negative/type rates, compression ratio, and various vector field visualizations with several scientific applications.

critical points↗