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

Magnetic Helicity of Alfven Simple Waves

The magnetic helicity of fully nonlinear, multi-dimensional Alfven simple waves are investigated, by using relative helicity formulae and also by using an approach involving poloidal and toroidal decomposition of the magnetic field and magnetic vector potential. Different methods to calculate the magnetic vector potential are used, including the homotopy and Biot-Savart formulas. Two basic Alfven modes are identified: (a) the plane 1D Alfven simple wave given in standard texts, in which the Alfven wave propagates along the z-axis, with wave phase varphi=k_0(z-lambda t), where k_0 is the wave number and lambda is the group velocity of the wave, and (b)\ the generalized Barnes (1976) simple Alfven wave in which the wave normal {bf n} moves in a circle in the xy-plane perpendicular to the mean field, which is directed along the z-axis. The plane Alfven wave (a) is analogous to the slab Alfven mode and the generalized Barnes solution (b) is analogous to the 2D mode in Alfvenic, incompressible turbulence. The helicity characteristics of these two basic Alfven modes are distinct. The helicity characteristics of more general multi-dimensional simple Alfven waves are also investigated. Applications to nonlinear Aifvenic fluctuations and structures observed in the solar wind are discussed.

Webb, Gary M.↗
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↗
NASA NTRS1977

Symbolic computer vector analysis

A MACSYMA program is described which performs symbolic vector algebra and vector calculus. The program can combine and simplify symbolic expressions including dot products and cross products, together with the gradient, divergence, curl, and Laplacian operators. The distribution of these operators over sums or products is under user control, as are various other expansions, including expansion into components in any specific orthogonal coordinate system. There is also a capability for deriving the scalar or vector potential of a vector field. Examples include derivation of the partial differential equations describing fluid flow and magnetohydrodynamics, for 12 different classic orthogonal curvilinear coordinate systems.

Stoutemyer, D. R.↗
NASA NTRS1986

Wall reflection of a viscous vortex ring

The behavior of a viscous axisymmetric vortex ring being reflected from a wall is investigated. The incompressible Navier-Stokes equations formulated in terms of the vorticity function and vector potential are numerically integrated by implicit finite difference methods. To specify the vector potential at a far boundary from the wall, the existing integral method used so far only for an unbounded domain is modified by a kind of image method. The trajectory of the vortex ring calcualted as a result closely resembles that observable from the experiment.

Sa, J. Y.↗
NASA NTRS1982

An application of the Braunbeck method to the Maggi-Rubinowicz field representation

The Braunbek method is applied to the generalized vector potential associated with the Maggi-rubinowicz representation. Under certain approximations, an asymptotic evaluation of the vector potential is obtained. For observation points away from caustics or shadow boundaries, the field derived from this quantity is the same as that determined from the geometrical theory of diffraction on a singly diffracted edge ray. An evaluation of the field for the simple case of a plane wave normally incident on a circular aperture is presented showing that the field predicted by the Maggi-Rubinowicz theory is continuous across the shadow boundary.

Meneghini, R.↗
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↗
NASA NTRS1989

A dual potential formulation of the Navier-Stokes equations

A dual potential formulation for numerically solving the Navier-Stokes equations is developed and presented. The velocity field is decomposed using a scalar and vector potential. Vorticity and dilatation are used as the dependent variables in the momentum equations. Test cases in two dimensions verify the capability to solve flows using approximations from potential flow to full Navier-Stokes simulations. A three-dimensional incompressible flow formulation is also described. An interesting feature of this approach to solving the Navier-Stokes equations is the decomposition of the velocity field into a rotational part (vector potential) and an irrotational part (scalar potential). The Helmholtz decomposition theorem allows this splitting of the velocity field. This approach has had only limited use since it increases the number of dependent variables in the solution. However, it has often been used for incompressible flows where the solution scheme is known to be fast and accurate. This research extends the usage of this method to fully compressible Navier-Stokes simulations by using the dilatation variable along with vorticity. A time-accurate, iterative algorithm is used for the uncoupled solution of the governing equations. Several levels of flow approximation are available within the framework of this method. Potential flow, Euler and full Navier-Stokes solutions are possible using the dual potential formulation. Solution efficiency can be enhanced in a straightforward way. For some flows, the vorticity and/or dilatation may be negligible in certain regions (e.g., far from a viscous boundary in an external flow). It is possible to drop the calculation of these variables then and optimize the solution speed. Also, efficient Poisson solvers are available for the potentials. The relative merits of non-primitive variables versus primitive variables for solution of the Navier-Stokes equations are also discussed.

Gegg, S. G.↗
NASA NTRS1997

Analysis of Discontinuities in a Rectangular Waveguide Using Dyadic Green's Function Approach in Conjunction with Method of Moments

The dyadic Green's function for an electric current source placed in a rectangular waveguide is derived using a magnetic vector potential approach. A complete solution for the electric and magnetic fields including the source location is obtained by simple differentiation of the vector potential around the source location. The simple differentiation approach which gives electric and magnetic fields identical to an earlier derivation is overlooked by the earlier workers in the derivation of the dyadic Green's function particularly around the source location. Numerical results obtained using the Green's function approach are compared with the results obtained using the Finite Element Method (FEM).

Deshpande, M. D.↗
NASA NTRS1986

Mutual coupling effects in antenna arrays, volume 1

Mutual coupling between rectangular apertures in a finite antenna array, in an infinite ground plane, is analyzed using the vector potential approach. The method of moments is used to solve the equations that result from setting the tangential magnetic fields across each aperture equal. The approximation uses a set of vector potential model functions to solve for equivalent magnetic currents. A computer program was written to carry out this analysis and the resulting currents were used to determine the co- and cross-polarized far zone radiation patterns. Numerical results for various arrays using several modes in the approximation are presented. Results for one and two aperture arrays are compared against published data to check on the agreement of this model with previous work. Computer derived results are also compared against experimental results to test the accuracy of the model. These tests of the accuracy of the program showed that it yields valid data.

Collin, R. E.↗
NASA NTRS1985

Collisionless tearing in a field-reversed sheet pinch assuming nonparallel propagation

The problem of collisionless linear tearing is examined assuming a wave vector with a component normal to the equilibrium field. The geometry is defined and the general form of the linear dispersion equation is calculated. The linear theory results when k is parallel to B are reviewed, and Ampere's law is calculated for the external adiabatic region when k times B does not equal zero, using two-fluid theory. A solution is obtained for the approximate form of the perturbed currents and vector potential assuming quasi-parallel k. The resonant current contributions within the singular layer are calculated, obtaining an estimate of the dispersion equation. The form of the adiabatic currents within the singular layer is calculated, showing that an x-z current system persists even in the limit k perpendicular to B goes to zero. Finally, the perturbed vector potential solutions across the singular layer are matched to obtain the shape of the complete eigenfunction.

Quest, K. B.↗
NASA NTRS1979

Solving the current element problem over lossy half-space without Sommerfeld integrals

An approach is presented for efficient computation of the vector potentials arising in the problem of a current element radiating over a lossy half-space. The present approach departs from the conventional ones in that it works primarily with the transform domain representations rather than with the Sommerfeld integrals which are the corresponding spatial domain counterparts. The key step in the present method is to approximate the transforms using a suitable approximation which is valid for a wide range of parameters of practical interest. The approximated transforms can be inverted in a closed form for the horizontal component of the vector potentials and can be expressed in a computationally efficient form for the vertical component. Numerical results illustrating the accuracy of the method are presented and some estimates of comparative computational times are also included.

Mittra, R.↗
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 ↗
NASA NTRS1988

A generalization of Petterssen's frontogenesis function and its relation to the forcing of vertical motion

The classic Peterssen frontogenesis function, defined as the Lagrangian rate of change of the magnitude of the potential temperature gradient, is generalized to apply to the vector potential temperature gradient. The derivation of vector frontogenesis in natural coordinates is shown, and analytical examples of frontal evolution in nondivergent horizontal velocity fields are given which suggest that both the frontogenetic and rotational components of the vector frontogenesis function F may be comparable in developing frontal zones observed in nature. The relative importance of the magnitude and direction contributions to F is quantitatively investigated, and it is found that the Lagrangian rates of change of the magnitude and direction of the potential temperature gradient are comparable. The frontal circulation is found to be related to the magnitude component of the Q vector, whereas the background circulation is related to the direction component.

Keyser, Daniel↗
NASA NTRS1989

Off disk-center potential field calculations using vector magnetograms

A potential field calculation for off disk-center vector magnetograms that uses all the three components of the measured field is investigated. There is neither any need for interpolation of grid points between the image plane and the heliographic plane nor for an extension or a truncation to a heliographic rectangle. Hence, the method provides the maximum information content from the photospheric field as well as the most consistent potential field independent of the viewing angle. The introduction of polarimetric noise produces a less tolerant extrapolation procedure than using the line-of-sight extrapolation, but the resultant standard deviation is still small enough for the practical utility of this method.

Venkatakrishnan, P.↗
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↗
NASA NTRS1992

Estimating high mosquito-producing rice fields using spectral and spatial data

The cultivation of irrigated rice provides ideal larval habitat for a number of anopheline vectors of malaria throughout the world. Anopheles freeborni, a potential vector of human malaria, is associated with the nearly 240,000 hectares of irrigated rice grown annually in Northern and Central California; therefore, this species can serve as a model for the study of rice field anopheline population dynamics. Analysis of field data revealed that rice fields with early season canopy development, that are located near bloodmeal sources (i.e., pastures with livestock) were more likely to produce anopheline larvae than fields with less developed canopies located further from pastures. Remote sensing reflectance measurements of early-season canopy development and geographic information system (GIS) measurements of distanes between rice fields and pastures with livestock were combined to distinguish between high and low mosquito-producing rice fields. Using spectral and distance measures in either a discriminant or Bayesian analysis, the identification of high mosquito-producing fields was made with 85 percent accuracy nearly two months before anopheline larval populations peaked. Since omission errors were also minimized by these approaches, they could provide a new basis for directing abatement techniques for the control of malaria vectors.

Wood, B. L.↗