Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Vector Potential”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Vector potential and surface magnetic field in magnetoelectric antiferromagnetic materials

A general formula for the average vector potential of bulk periodic systems is proposed and shown to set the boundary conditions at magnetic interfaces. Here, for antiferromagnetic materials, the study reveals a unique relation between the macroscopic potential and the orientation-dependent magnetic quadrupole as a result of the different crystalline and magnetic symmetries. In particular, at surfaces and interfaces of a truncated bulk without inversion and time-reversal symmetries, the average vector potential exhibits a discontinuity which results in an interfacial magnetic field. In general, however, due to the surface and interface electronic and atomic relaxations, additional magnetization may result. For the experimentally observed magnetoelectric antiferromagnets, in particular, our symmetry analysis suggests that the relaxation effects could well be a system response to the presence of such a potential discontinuity.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

High order interpolation of magnetic fields with vector potential reconstruction for particle simulations

We propose a method for interpolating divergence-free continuous magnetic fields via vector potential reconstruction using Hermite interpolation, which ensures high-order continuity for applications requiring adaptive, high-order ordinary differential equation (ODE) integrators, such as the Dormand-Prince method. The method provides C(m) continuity and achieves high-order accuracy, making it particularly suited for particle trajectory integration and Poincaré section analysis under optimal integration order and timestep adjustments. Through numerical experiments, we demonstrate that the Hermite interpolation method preserves volume and continuity, which are critical for conserving toroidal canonical momentum and magnetic moment in guiding center simulations, especially over long-term trajectory integration. Furthermore, we analyze the impact of insufficient derivative continuity on Runge-Kutta schemes and show how it degrades accuracy at low error tolerances, introducing discontinuity-induced truncation errors. Lastly, we demonstrate performant Poincaré section analysis in two relevant settings of field data collocated from finite element meshes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Comparative investigation of iterative solutions of the all-frequency stable formulation and its vector-potential-only variation

An all-frequency stable formulation has been proposed and applied to solve low-frequency and multiscale electromagnetic problems to avoid the low-frequency breakdown catastrophe of the traditional vector wave equation. Based on the potential representation of the fields, this formulation involves vector and scalar potentials, as well as an auxiliary potential to enforce the gauge condition. In this paper, such a formulation is first simplified to involve vector potential only. Their iterative solutions are then sought by using an incomplete LU preconditioner and an iterative solver. The iterative performance of both formulations are investigated in a comparative study.

Mekonnen, Minyechil↗
DOE OSTI2025

Superfluid weight of strongly inhomogeneous superconductors

In this study, we obtain the superfluid weight of strongly inhomogeneous superconductors. We find that, in general, the response of the superconductor’s pairing potential to the perturbing vector potential is essential to obtain the correct value of the superfluid weight. We consider two exemplary cases: the case when strong inhomogeneities in the pairing potential are induced by a periodic potential, and the case when superconducting vortices are induced by an external magnetic field. For both cases we show that the correction to the superfluid weight due to the response of the paring potential to the perturbing vector potential can be significant, it must be included to obtain quantitatively correct results, and that for the case when vortices are present the expression of the superfluid weight that does not include such correction returns qualitatively wrong results. The results are directly relevant for superconducting states in moiré multilayer where the almost inevitable spatial variations of the magic angle and the proximity of superconductivity to anisotropic phases induce inhomogeneities in the pairing potential, and to the efforts to optimally pin superconducting vortices to achieve large superconducting critical currents.

Bogoliubov-de Gennes equations↗
DOE OSTI2024

A Particle-in-Cell Method for Plasmas with a Generalized Momentum Formulation, Part II: Enforcing the Lorenz Gauge Condition

In a previous paper Christlieb et al. (A particle-in-cell method for plasmas with a generalized momentum formulation, part I: Model formulation, 2024), we developed a new particle-in-cell (PIC) method for the relativistic Vlasov–Maxwell system in which the electromagnetic fields and the equations of motion for the particles were cast in terms of scalar and vector potentials through a Hamiltonian formulation. This new method evolved the potentials under the Lorenz gauge using integral equation methods. New methods to construct spatial derivatives of the potentials that converge at the same rates as the fields were also presented. The new particle method was compared against standard explicit discretizations, including the well-known FDTD-PIC method, for a range of applications involving sheaths and particle beams. Here, this paper extends this new class of methods by focusing on the enforcement the Lorenz gauge condition in both exact and approximate forms using co-located meshes. A time-consistency property of the proposed field solver for the vector potential form of Maxwell’s equations is established, which is shown to preserve the equivalence between the semi-discrete Lorenz gauge condition and the analogous semi-discrete continuity equation. Using this property, we present three methods to enforce a semi-discrete gauge condition. The first method introduces an update for the continuity equation that is consistent with the discretization of the Lorenz gauge condition. Both the finite difference and spectral implementations satisfy this discrete gauge condition to machine precision. The second approach we propose enforces a semi-discrete continuity equation using the boundary integral solution to the field equations. The potential benefit of this approach is that it eliminates spatial derivatives that appear on the particle data, namely the current density, which is often calculated by linear combinations of low-order spline basis functions. This method is ideally suited to boundary integral equation methods that invert multi-dimensional operators without dimensional splitting techniques and will be the subject of future work. The third approach introduces a gauge correcting method that makes direct use of the gauge condition to modify the scalar potential and uses local maps for both the charge and current densities. This results in a gauge error, as the maps do not enforce the continuity equation. The vector potential coming from the current density is taken to be exact, and using the Lorenz gauge, we compute a correction to the scalar potential that makes the two potentials satisfy the gauge condition. This method also enforces the gauge condition to machine precision. We demonstrate two of the proposed methods in the context of periodic domains. Problems defined on bounded domains, including those with complex geometric features remain an ongoing effort. However, this work shows that it is possible to design computationally efficient methods that can effectively enforce the Lorenz gauge condition in a non-staggered PIC formulation.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Engineering and revealing Dirac strings in spinor condensates

Artificial monopoles have been engineered in various systems, yet there has been no systematic study of the singular vector potentials associated with the monopole field. We show that the Dirac string, the line singularity of the vector potential, can be engineered, manipulated, and made manifest in a spinor atomic condensate. We elucidate the connection among spin, orbital degrees of freedom, and the artificial gauge, and show that there exists a mapping between the vortex filament and the Dirac string. We also devise a proposal where preparing initial spin states with relevant symmetries can result in different vortex patterns, revealing an underlying correspondence between the internal spin states and the spherical vortex structures. Such a mapping also leads to a new way of constructing spherical Landau levels, and monopole harmonics. Our observation provides insights into the behavior of quantum matter possessing internal symmetries in curved spaces. Published by the American Physical Society 2024

Xu, Gui-Sheng↗
DOE OSTI2025

Field–potential finite-difference time-domain (FiPo FDTD) technique for computational electromagnetics

Modeling light–matter interactions at the nanoscale requires accurate handling of coupled quantum and electromagnetic systems. This coupling requires information about the electric scalar potential Φ and the magnetic vector potential A, which are not typically calculated in standard computational electromagnetics implementations. To that end, we have developed a field–potential finite-difference time-domain (FiPo FDTD) algorithm, which solves a set of first-order equations for Φ and A alongside equations for the electric and magnetic fields E and H. The FiPo Basic code is essentially conventional FDTD, but with an added module that calculates the potentials. The FiPo Hybrid code self-consistently calculates both fields and potentials and is particularly suitable for coupling with quantum electronic transport solvers because it can be sourced by the potentials themselves. To terminate the domain and mimic infinite space, we have derived and implemented a convolutional perfectly matched layer (CPML) absorbing boundary condition for FiPo FDTD whose performance is on par with state-of-the-art CPMLs for standard FDTD. We present FiPo simulation results on several example systems.

Avazpour, L. [University of Wisconsin-Madison, WI ↗
DOE OSTI2024

NSTXU Diagnostic Disruption Dynamic Loading Represented by Response Spectra

This article presents the results of transient dynamic simulations of loads due to disruption eddy currents on the NSTXU vacuum vessel. Dynamic loading at diagnostic mounting locations is expressed as response spectra derived from the time history results of the dynamic structural simulations of a variety of disruption scenarios. The disruption simulations draw on a history of the project assessments of worst case disruptions for specific components. Major efforts to assess disruption loading have included the vacuum vessel which is the major structural support for the machine, as well as the passive plates (PPs), high harmonic fast wave (HHFW) antenna, and centerstack casing. Each one of these efforts included transient electromagnetic simulations producing time-dependent eddy current Lorentz loads (and in some cases halo loads) which then were applied to time-dependent structural dynamic analyses intended to obtain the proper dynamic amplification factors. In some instances, the EM model and structural model were identical allowing direct transfer of EM forces to the structural model. In other cases, the EM and structural model were not identical and the vector potential (VP) transfer method was used. The results files from these analyses were available (or re-run) to post process in ANSYS Classic time history postprocessor. In conclusion, the ANSYS command is used to create response spectra from time history data at desired points on the vessel.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Chaotic and integrable magnetic fields in one-dimensional hybrid Vlasov–Maxwell equilibria

The construction of kinetic equilibrium states is important for studying stability and wave propagation in collisionless plasmas. Thus, many studies over the past decades have been focused on calculating Vlasov–Maxwell equilibria using analytical and numerical methods. However, the problem of kinetic equilibrium of hybrid models is less studied, and self-consistent treatments often adopt restrictive assumptions ruling out cases with irregular and chaotic behaviour, although such behaviour is observed in spacecraft observations of space plasmas. In this paper, we develop a one-dimensional (1-D), quasineutral, hybrid Vlasov–Maxwell equilibrium model with kinetic ions and massless fluid electrons and derive associated solutions. The model allows for an electrostatic potential that is expressed in terms of the vector potential components through the quasineutrality condition. The equilibrium states are calculated upon solving an inhomogeneous Beltrami equation that determines the magnetic field, where the inhomogeneous term is the current density of the kinetic ions and the homogeneous term represents the electron current density. We show that the corresponding 1-D system is Hamiltonian, with position playing the role of time, and its trajectories have a regular, periodic behaviour for ion distribution functions that are symmetric in the two conserved particle canonical momenta. For asymmetric distribution functions, the system is nonintegrable, resulting in irregular and chaotic behaviour of the fields. The electron current density can modify the magnetic field phase space structure, inducing orbit trapping and the organization of orbits into large islands of stability. Thus, the electron contribution can be responsible for the emergence of localized electric field structures that induce ion trapping. We also provide a paradigm for the analytical construction of hybrid equilibria using a rotating two-dimensional harmonic oscillator Hamiltonian, enabling the calculation of analytic magnetic fields and the construction of the corresponding distribution functions in terms of Hermite polynomials.

Kaltsas, Dimitrios A. (ORCID:0000000300769015)↗
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↗
DOE OSTI2025

First Principles Cable Braid Electromagnetic Penetration Model - Potential Impedance Clarification

In this report the lack of assumed connectivity between braid carriers imposes the need to use the original magnetic vector potential formulas for transfer and self impedances rather than the final simplified results in the magnetic coupling and propagation formulas. Each of the magnetic sections in the previous journal article and Sandia report are addressed with references to the appropriate formulas.

42 ENGINEERING↗
DOE OSTI2023

Simulations Supporting the Development of Northstar's Indirect Beam Parameters Monitoring System

NorthStar Medical Radioisotopes, LLC is planning to produce the important medical radioisotope molybdenum-99 (Mo-99), the parent of technetium-99m (Tc-99m), through photonuclear reactions in molybdenum-100 (Mo-100). In this approach, a target comprising multiple thin disks of enriched molybdenum metal is bombarded with a 40-MeV electron beam. Electrons impinged on the molybdenum target produce bremsstrahlung X-rays that cause the nuclear reaction. Because enriched Mo-100 is expensive, there is a desire to utilize as much beam power as possible to achieve maximum production yield and minimize the size of the target. This requirement leads to very high beam power density (and heat deposition in the target), which creates challenging requirements for the cooling of the target. The critical part of the target is the target window. It separates the high-pressure helium cooled target from the vacuum beamline and the subject of structural and thermal stress. The temperature of the target window is proportional to the energy density deposited by the beam, so it is critical to maintain the desired beam profile on the target window. The feasibility of indirectly monitoring the maximum energy density of the beam on the beam window through beam losses at the main collimator (Collimator) before the production target was verified. A model of the NorthStar beam transport line was constructed for this purpose using MAD-X and Tao/Bmad codes. Beam optics were computed for the standard operational scenario, followed by an investigation involving approximately 400 cases with parameter variations in the last tuning quadrupoles. This was done to assess the correlation between losses in the collimator and the peak energy density on the target. We developed a model to explore the potential application of Optical transition radiation (OTR) for controlling beam parameters in the NorthStar beam delivery system. This model was based on a generic formula derived from the fundamental solution of the inhomogeneous wave equation of the vector potential, and allowed us to consider various surfaces, even those with irregular or random features, using numerical integration. We applied the model to OTR generated by relativistic electrons impacting an Inconel® 718 beam window. We examined cases with different levels of the window’s surface roughness, ranging from 0.5 to 3.0 microns of root square mean (RMS) deviation. The results of the OTR simulations provided distributions of OTR photons that can be used to study the limitations of optical systems for controlling beam parameters.

43 PARTICLE ACCELERATORS↗
DOE OSTI2021

Robust and ultrafast state preparation by ramping artificial gauge potentials

Abstract The implementation of static artificial magnetic fields in ultracold atomic systems has become a powerful tool, e.g. for simulating quantum-Hall physics with charge-neutral atoms. Taking an interacting bosonic flux ladder as a minimal model, we investigate protocols for adiabatic state preparation via magnetic flux ramps. Considering the fact that it is actually the artificial vector potential (in the form of Peierls phases) that can be experimentally engineered in optical lattices, rather than the magnetic field, we find that the time required for adiabatic state preparation dramatically depends on which pattern of Peierls phases is used. This can be understood intuitively by noting that different patterns of time-dependent Peierls phases that all give rise to the same magnetic field ramp, generally lead to different artificial electric fields during the ramp. As an intriguing result, we find that an optimal choice allows for preparing the ground state almost instantaneously in the non-interacting system, which can be related to the concept of counterdiabatic driving. Remarkably, we find extremely short preparation times also in the strongly-interacting regime. Our findings open new possibilities for robust state preparation in atomic quantum simulators.

Wang, Botao (ORCID:0000000282202452)↗
DOE OSTI2022

Topological monopole's gauge field‐induced anomalous Hall effect in artificial honeycomb lattice

Abstract Vortex magnetic structure in artificial honeycomb lattice provides a unique platform to explore emergent properties due to the additional Berry phase curvature imparted by chiral magnetization to circulating electrons via direct interaction. We argue that while the perpendicularly aligned magnetic component leads to the quantized flux of monopole at the center of the Berry sphere, the in‐plane vortex circulation of magnetization gives rise to unexpected nontrivial topological Berry phase due to the gauge field transformation. The unprecedented effect signifies the importance of vector potential in multiply connected geometrical systems. Experimental confirmations to proposed hypotheses are obtained from Hall resistance measurements on permalloy honeycomb lattice. Investigation of the topological gauge transformation due to the in‐plane chirality reveals anomalous quasi‐oscillatory behavior in Hall resistance as a function of perpendicular field. The oscillatory nature of is owed to the fluctuation in equilibrium current as a function of Fermi wave‐vector , envisaged under the proposed new formulation in this article. Our synergistic approach suggests that artificially tunable nanostructured material provides new vista to the exploration of topological phenomena of strong fundamental importance. Key Points Demonstration of gauge field flux due to vortex magnetism. Anomalous Hall effect in magnetic honeycomb lattice. Topological monopole‐induced magnetic flux.

Guo, Jiasen↗
DOE OSTI2024

Charge Conservation beyond Uniformity: Spatially Inhomogeneous Electromagnetic Response in Periodic Solids

Nonlinear electromagnetic response functions have reemerged as a crucial tool for studying quantum materials, due to recently appreciated connections between optical response functions, quantum geometry, and band topology. Most attention has been paid to responses to spatially uniform electric fields, relevant to low-energy optical experiments in conventional solid state materials. However, magnetic and magnetoelectric phenomena are naturally connected by responses to spatially varying electric fields due to Maxwell’s equations. Furthermore, in the emerging field of moiré materials, characteristic lattice scales are much longer, allowing spatial variation of optical electric fields to potentially have a measurable effect in experiments. In order to address these issues, we develop a formalism for computing linear and nonlinear responses to spatially inhomogeneous electromagnetic fields. Starting with the continuity equation, we derive an expression for the second-quantized current operator that is manifestly conserved and model independent. Crucially, our formalism makes no assumptions on the form of the microscopic Hamiltonian and so is applicable to model Hamiltonians derived from tight-binding or calculations. We then develop a diagrammatic Kubo formalism for computing the wave vector dependence of linear and nonlinear conductivities, using Ward identities to fix the value of the diamagnetic current order by order in the vector potential. We apply our formula to compute the magnitude of the Kerr effect at oblique incidence for a model of a moiré-Chern insulator and demonstrate the experimental relevance of spatially inhomogeneous fields in these systems. We further show how our formalism allows us to compute the (orbital) magnetic multipole moments and magnetic susceptibilities in insulators. Turning to nonlinear response, we use our formalism to compute the second-order transverse response to spatially varying transverse electric fields in our moiré-Chern insulator model, with an eye toward the next generation of experiments in these systems. Published by the American Physical Society 2024

Physics↗
DOE OSTI2020

Meshfree methods on manifolds for hydrodynamic flows on curved surfaces: A Generalized Moving Least-Squares (GMLS) approach

We utilize generalized moving least squares (GMLS) to develop meshfree techniques for discretizing hydrodynamic flow problems on manifolds. We use exterior calculus to formulate incompressible hydrodynamic equations in the Stokesian regime and handle the divergence-free constraints via a generalized vector potential. This provides less coordinate-centric descriptions and enables the development of efficient numerical methods and splitting schemes for the fourth-order governing equations in terms of a system of second-order elliptic operators. Using a Hodge decomposition, we develop methods for manifolds having spherical topology. We show the methods exhibit high-order convergence rates for solving hydrodynamic flows on curved surfaces. The methods also provide general high-order approximations for the metric, curvature, and other geometric quantities of the manifold and associated exterior calculus operators. The approaches also can be utilized to develop high-order solvers for other scalar-valued and vector-valued problems on manifolds.

97 MATHEMATICS AND COMPUTING↗