Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ELECTROMAGNETIC 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 271 records · Page 15

Work and energy for particles in electromagnetic field

Defining the energy and work for particles interacting with electromagnetic field (EMF) is an open problem, because—due to the gauge-freedom—there exist various non-equivalent possibilities. It is argued that a consistent definition can be provided via the Lorenz gauge. To this end, I work out a system of two electromagnetically coupled classical particles. One of them is much heavier and models the source of work. The definition of energy in the Lorenz gauge is causal and consistent, because it leads to an approximate conservation law due to which the work done by the heavy particle (source of work) can be defined either via the kinetic energy of the heavy particle, or via the full time-dependent energy (kinetic + potential in the Lorenz gauge) of the light particle.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

High-Gradient Wien Spin Rotators at Jefferson Lab

Nuclear physics experiments performed in the Continuous Electron Beam Accelerator Facility (CEBAF) at Jefferson Laboratory (JLab) require spin manipulation of electron beams. Two Wien spin rotators in the injector keV region are essential at CEBAF to establish longitudinal polarization at the end station target, and to flip the polarization direction by ? rad to rule out false asymmetries. In a Wien filter, the homogeneous and independent electric and magnetic fields, along with the velocity vectors of the electrons that traverse it, form a mutually orthogonal system. The magnitude of the electrostatic field, established by biasing two highly-polished electrodes, defines the desired spin angle at the target yet deviates the beam trajectory due to the Lorentz force. The beam trajectory in the Wien is then re-established by adjusting the magnetic field, induced by an electromagnet encasing the device vacuum chamber. This contribution describes the evolution design and high voltage testing of Wien filters for spin manipulation at increased beam energies in the keV injector region, required by high precision parity violation experiments like MOLLER.

Palacios Serrano, G. G.↗

High-Gradient Wien Spin Rotators at Jefferson Lab

Nuclear physics experiments performed in the Continuous Electron Beam Accelerator Facility (CEBAF) at Jefferson Laboratory (JLab) require spin manipulation of electron beams. Two Wien spin rotators in the injector keV region are essential at CEBAF to establish longitudinal polarization at the end station target, and to flip the polarization direction by ? rad to rule out false asymmetries. In a Wien filter, the homogeneous and independent electric and magnetic fields, along with the velocity vectors of the electrons that traverse it, form a mutually orthogonal system. The magnitude of the electrostatic field, established by biasing two highly-polished electrodes, defines the desired spin angle at the target yet deviates the beam trajectory due to the Lorentz force. The beam trajectory in the Wien is then re-established by adjusting the magnetic field, induced by an electromagnet encasing the device vacuum chamber. This contribution describes the evolution design and high voltage testing of Wien filters for spin manipulation at increased beam energies in the keV injector region, required by high precision parity violation experiments like MOLLER.

Palacios Serrano, G. G.↗

From Bessel beam to complex-source-point cylindrical wave-function

This investigation shows that a scalar Bessel beam can be transformed into the non-paraxial complex-source-point cylindrical wave (CSPCW). High-order CSPCW solutions, termed here high-order quasi-Gaussian cylindrical beams, which exactly satisfy the Helmholtz equation, are derived analytically. Moreover, partial-derivatives of the high-order CSPCW solutions satisfy the Helmholtz equation. In addition, the CSPCW solutions satisfy the nonrelativistic Schrödinger equation within standard quantum mechanics, thus, the results can be used in the description of elementary particle/matter motion and related applications in quantum scattering theory. Furthermore, the analysis is extended to the case of vector beams in which the components of the electromagnetic (EM) field are obtained based on different polarizations of the magnetic and electric vector potentials, which exactly satisfy Maxwell’s vectorial equations and Lorenz’ gauge condition. An attractive feature of the high-order solutions is the rigorous description of strongly focused (or strongly divergent) cylindrical wave-fields without any approximations, nor the need for numerical methods. Possible applications are in beam-forming design using high-aperture or collimated cylindrical laser/electron quasi-Gaussian beams in imaging microscopy, particle manipulation, optical tweezers, and the study of the scattering, and radiation forces on objects. - Highlights: • Bessel beam is transformed into the non-paraxial cylindrical complex-source-point. • Exact high-order tightly focused solutions are derived without any approximations. • The exact solutions also satisfy the nonrelativistic Schrödinger equation. • Electromagnetic beams are obtained as solutions of Maxwell’s vectorial equations. • Applications are in laser/electron beam imaging, tweezers, and radiation force.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Two-temperature-electron Zakharov equations for the dynamics of modulated collapse

Under the framework of the two-timescale method, we derived the two-temperature-electron Zakharov equations which can describe the wave-wave interaction among Langmuir wave, electromagnetic wave, and electron acoustic waves. Numerical results reveal that the modulation instability will lead the wave field collapse into a turbulent state. The spectrum obtained from the Fourier analysis shows a cascade of energy from large to small scales. In the collapse processes, the hot electrons are expelled by the ponderomotive force induced by the high frequency waves, and then density cavitons are formed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Magnetohydrodynamic-particle-in-cell method for coupling cosmic rays with a thermal plasma: application to non-relativistic shocks

We formulate a magnetohydrodynamic-particle-in-cell (MHD-PIC) method for describing the interaction between collisionless cosmic ray (CR) particles and a thermal plasma. The thermal plasma is treated as a fluid, obeying equations of ideal MHD, while CRs are treated as relativistic Lagrangian particles subject to the Lorentz force. Backreaction from CRs to the gas is included in the form of momentum and energy feedback. In addition, we include the electromagnetic feedback due to CR-induced Hall effect that becomes important when the electron–ion drift velocity of the background plasma induced by CRs approaches the Alfvén velocity. Our method is applicable on scales much larger than the ion inertial length, bypassing the microscopic scales that must be resolved in conventional PIC methods, while retaining the full kinetic nature of the CRs. We have implemented and tested this method in the Athena MHD code, where the overall scheme is second-order accurate and fully conservative. As a first application, we describe a numerical experiment to study particle acceleration in non-relativistic shocks. Using a simplified prescription for ion injection, we reproduce the shock structure and the CR energy spectra obtained with more self-consistent hybrid-PIC simulations, but at substantially reduced computational cost. We also show that the CR-induced Hall effect reduces the growth rate of the Bell instability and affects the gas dynamics in the vicinity of the shock front. As a step forward, we are able to capture the transition of particle acceleration from non-relativistic to relativistic regimes, with momentum spectrum f(p)∝p{sup −4} connecting smoothly through the transition, as expected from the theory of Fermi acceleration.

79 ASTRONOMY AND ASTROPHYSICS↗

Force-free motion of a cold plasma during the current quench

Cold disruptive plasma tends to move during the current quench. Its motion is essentially force-free since the current quench timescale is resistive rather than Alfvénic. In contrast to the hot vertical displacement events, the frozen-in condition is violated in the cold plasma case, and the plasma motion is not governed by magnetic flux conservation but rather by its dissipation. In this paper, we present a numerical model of the cold plasma dynamics. This model predicts electromagnetic loads on the vacuum vessel, the plasma flow and density evolution, and the plasma centroid evolution. Our calculations include poloidal wall currents. Finally, we demonstrate their significant contribution to the force acting on the vacuum vessel.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Stochastic regimes in the driven oscillator with a step-like nonlinearity

A nonlinear oscillator with an abruptly inhomogeneous restoring force driven by an uniform oscillating force exhibits stochastic properties under specific resonance conditions. This behaviour elucidates the elementary mechanism of the electron energization in the strong electromagnetic wave interaction with thin targets.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

SU-F-J-10: Sliding Mode Control of a SMA Actuated Active Flexible Needle for Medical Procedures

Purpose: In medical interventional procedures such as brachytherapy, ablative therapies and biopsy precise steering and accurate placement of needles are very important for anatomical obstacle avoidance and accurate targeting. This study presents the efficacy of a sliding mode controller for Shape Memory Alloy (SMA) actuated flexible needle for medical procedures. Methods: Second order system dynamics of the SMA actuated active flexible needle was used for deriving the sliding mode control equations. Both proportional-integral-derivative (PID) and adaptive PID sliding mode control (APIDSMC) algorithms were developed and implemented. The flexible needle was attached at the end of a 6 DOF robotic system. Through LabView programming environment, the control commands were generated using the PID and APIDSMC algorithms. Experiments with artificial tissue mimicking phantom were performed to evaluate the performance of the controller. The actual needle tip position was obtained using an electromagnetic (EM) tracking sensor (Aurora, NDI, waterloo, Canada) at a sampling period of 1ms. During experiment, external disturbances were created applying force and thermal shock to investigate the robustness of the controllers. Results: The root mean square error (RMSE) values for APIDSMC and PID controllers were 0.75 mm and 0.92 mm, respectively, for sinusoidal reference input. In the presence of external disturbances, the APIDSMC controller showed much smoother and less overshooting response compared to that of the PID controller. Conclusion: Performance of the APIDSMC was superior to the PID controller. The APIDSMC was proved to be more effective controller in compensating the SMA uncertainties and external disturbances with clinically acceptable thresholds.

60 APPLIED LIFE SCIENCES↗

Numerical modeling of electromagnetic field spatiotemporal evolution to evaluate the effects on calcium carbonate crystallization

Calcium carbonate (CaCO 3 ) scaling is a significant impediment to water systems. Electromagnetic field (EMF) treatment is a promising approach to control scaling owing to its simplicity and low or no energy requirements. However, the underlying mechanisms by which EMF impacts CaCO 3 crystallization remain unclear due to the challenges in measuring the EMFs in feed solutions and the lack of a fundamental understanding of the applied EMFs and the observed physicochemical phenomena. To fill this knowledge gap, a high-fidelity COMSOL model was first developed to simulate EMFs in bulk solutions for three alternating current-induced EMF devices with different configurations and properties. These were then integrated with experimental data to unveil the underlying mechanism by which applied EMFs alter the physicochemical processes. The study revealed that even low-strength EMFs (e.g., electric fields <0.15 V/m and magnetic fields <0.03 mT) promoted CaCO 3 precipitation in bulk solutions. The electric fields created by these EMF devices resulted in higher Lorentz force compared to their induced magnetic fields. In conclusion, the methodology of this study offers the capability to predict the effectiveness of different EMF devices in facilitating crystallization processes, and these mechanistic insights lay the foundation for the smart design of EMF devices for diverse water treatment applications.

COMSOL simulation↗

Electromagnetic and neutral-weak response functions of 4 He and 12 C

A major goal of nuclear theory is to understand the strong interaction in nuclei as it manifests itself in terms of two- and many-body forces among the nuclear constituents, the protons and neutrons, and the interactions of these constituents with external electroweak probes via one- and many-body currents. The objective of the present work is to calculate the quasielastic electroweak response functions in light nuclei within the realistic dynamical framework outlined above. These response functions determine the inclusive cross section as function of the lepton momentum and energy transfers. Their ab initio calculation is a very challenging quantum many-body problem, since it requires summation over the entire excitation spectrum of the nucleus and inclusion in the electroweak currents of one- and many-body terms. Green's functions Monte Carlo methods allow one to circumvent both difficulties by computing the response in imaginary time (the so-called Euclidean response) and hence summing implicitly over the bound and continuum states of the nucleus, and by implementing specific algorithms designed to deal with the complicated spin-isospin structure of nuclear many-body operators. Theoretical predictions for 4 He and 12 C, confirmed by experiment in the electromagnetic case, show that two-body currents generate excess transverse strength from threshold to the quasielastic to the dip region and beyond. In conclusion, these results challenge the conventional picture of quasielastic inclusive scattering as being largely dominated by single-nucleon knockout processes.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Cyclotron resonance in plasma flow

This paper is devoted to the mechanism of particle acceleration via resonant interaction with the electromagnetic circular wave propagating along the inhomogeneous background magnetic field in the presence of a plasma flow. We consider the system where the plasma flow velocity is large enough to change the direction of wave propagation in the rest frame. This system mimics a magnetic field configuration typical for inner structure of a quasi-parallel shock wave. We consider conditions of gyroresonant interaction when the force corresponding to an inhomogeneity of the background magnetic field is compensated by the Lorentz force of the wave-magnetic field. The wave-amplitude is assumed to be about 10% of the background magnetic field. We show that particles can gain energy if kv{sub sw}>ω>kv{sub sw}−Ω{sub c} where k is the wave number, v{sub sw} is a plasma flow velocity, and ω and Ω{sub c} are the wave frequency and the particle gyrofrequency, respectively. This mechanism of acceleration resembles the gyrosurfing mechanism, but the effect of the electrostatic field is replaced by the effect of the magnetic field inhomogeneity.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A per-cent-level determination of the nucleon axial coupling from quantum chromodynamics

The axial coupling of the nucleon, g A , is the strength of its coupling to the weak axial current of the standard model of particle physics, in much the same way as the electric charge is the strength of the coupling to the electromagnetic current. This axial coupling dictates the rate at which neutrons decay to protons, the strength of the attractive long-range force between nucleons and other features of nuclear physics. Precision tests of the standard model in nuclear environments require a quantitative understanding of nuclear physics that is rooted in quantum chromodynamics, a pillar of the standard model. The importance of g A makes it a benchmark quantity to determine theoretically—a difficult task because quantum chromodynamics is non-perturbative, precluding known analytical methods. Lattice quantum chromodynamics provides a rigorous, non-perturbative definition of quantum chromodynamics that can be implemented numerically. It has been estimated that a precision of two per cent would be possible by 2020 if two challenges are overcome: contamination of g A from excited states must be controlled in the calculations and statistical precision must be improved markedly. Here we use an unconventional method inspired by the Feynman–Hellmann theorem that overcomes these challenges. Furthermore, we calculate a g A value of 1.271 ± 0.013, which has a precision of about one per cent.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Recovery of hexagonal Si-IV nanowires from extreme GPa pressure

We use Raman spectroscopy in tandem with transmission electron microscopy and density functional theory simulations to show that extreme (GPa) pressure converts the phase of silicon nanowires from cubic (Si-I) to hexagonal (Si-IV) while preserving the nanowire's cylindrical morphology. In situ Raman scattering of the longitudinal transverse optical (LTO) mode demonstrates the high-pressure Si-I to Si-II phase transition near 9 GPa. Raman signal of the LTO phonon shows a decrease in intensity in the range of 9–14 GPa. Then, at 17 GPa, it is no longer detectable, indicating a second phase change (Si-II to Si-V) in the 14–17 GPa range. Recovery of exotic phases in individual silicon nanowires from diamond anvil cell experiments reaching 17 GPa is also shown. Raman measurements indicate Si-IV as the dominant phase in pressurized nanowires after decompression. Transmission electron microscopy and electron diffraction confirm crystalline Si-IV domains in individual nanowires. Computational electromagnetic simulations suggest that heating from the Raman laser probe is negligible and that near-hydrostatic pressure is the primary driving force for the formation of hexagonal silicon nanowires.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Topics in Theoretical Physics

This award supported a broadly based research effort in theoretical particle physics, including research aimed at uncovering the laws of nature at short (subatomic) and long (cosmological) distances. These theoretical developments apply to experiments in laboratories such as CERN, the facility that operates the Large Hadron Collider outside Geneva, as well as to cosmological investigations done using telescopes and satellites. The results reported here apply to physics beyond the so-called Standard Model of particle physics; physics of high energy collisions such as those observed at the Large Hadron Collider; theoretical and mathematical tools and frameworks for describing the laws of nature at short distances; cosmology and astrophysics; and analytic and computational methods to solve theories of short distance physics. Some specific research accomplishments include + Theories of the electroweak interactions, the forces that give rise to many forms of radioactive decay; + Physics of the recently discovered Higgs boson. + Models and phenomenology of dark matter, the mysterious component of the universe, that has so far been detected only by its gravitational effects. + High energy particles in astrophysics and cosmology. + Algorithmic research and Computational methods for physics of and beyond the Standard Model. + Theory and applications of relativity and its possible limitations. + Topological effects in field theory and cosmology. + Conformally invariant systems and AdS/CFT. This award also supported significant training of students and postdoctoral fellows to lead the research effort in particle theory for the coming decades. These students and fellows worked closely with other members of the group as well as theoretical and experimental colleagues throughout the physics community. Many of the research projects funded by this grant arose in response to recently obtained experimental results in the areas of particle physics and cosmology. We describe a few of these below. Relativity is founded on a symmetry property of nature called "Lorentz Invariance". Like all symmetry properties, it is essential to determine precisely how symmetric nature actually is; that is, do the laws of nature fully respect the symmetry or is there room for tiny symmetry violating effects? An important consequence of Lorentz invariance is the existence of a universal limiting velocity for all physical particles. Light travels at this limiting velocity so it is frequently referred to as simply "the speed of light", but relativity requires that ALL particles travel more slowly than this speed. Once the Higgs particle was discovered in 2012 a natural question was whether or not this particle's speed was consistent with relativity. Although the speed of the Higgs particle is not measurable directly, Cohen has shown that, if the maximal speed of the Higgs particle was not precisely the same as the speed of light, then the Higgs would have some unusual properties. In some cases the Higgs would be unstable to some unusual decay modes; in other cases the interactions of the Higgs with other particles would change the properties of these other particles in ways that could be observed in so-called cosmic rays, very energetic particles (such as photons, protons and other atomic nuclei) coming from space. Once these particles hit the upper atmosphere they produce a "shower" of particles that can be seen by ground-based instruments. If the Higgs has a maximal speed that differs even a tiny bit from the speed of light these showers would look quite different from what is observed. In this way Cohen was able to establish that the Higgs travels with a maximal speed that cannot differ from the speed of light by more than one part in a thousand-trillion. This is by far the most precisely determined property of the Higgs particle. Cohen and Schmaltz reviewed evidence from the Large Hadron Collider (LHC), a particle physics experiment operating at the CERN laboratory near Geneva, for a new particle sometimes called a W'. This evidence included certain unexpected by-products in collisions of protons at very high energy. While the evidence was not significant enough to claim a discovery, it was sufficiently intriguing that many particle theorists worked to construct explanations for this signal. Cohen and Schmaltz were able to determine that such explanations are highly constrained by previous experiments involving collisions of very energetic particles. Nevertheless they were able to construct a theory that adequately explains the LHC data and remain consistent with prior experiments. Their explanation predicts the existence of yet another new particle, called a Z', with a mass slightly greater than that of the W'. This additional particle, if it exists, should be seen as more data is collected from the LHC. Amusingly, there is one collision by-product that has already been seen by the CMS experiment at the LHC that supports the existence of this new particle; however, it is not unlikely that this single event is a so-called "background" event, that is a somewhat atypical by-product of a conventional Standard Model process. This theory for the anomalous LHC data will either be confirmed or excluded with further data-taking at the LHC. The ratio of the number of electrons produced in bottom quark decays over the number of muons produced has been measured at the LHC. This ratio is interesting because it can be predicted very precisely from a basic property of the Standard Model: lepton universality. If lepton universality is correct, the ratio of electrons to muons is predicted to be equal to 1. The first measurements of this ratio find a value different from 1 with a statistical significance of about 3 standard deviations. Schmaltz and collaborator proposed a new extension of the Standard Model which can explain the new data. In addition, Schmaltz and collaborators proposed several new measurements of ratios of decay rates which can confirm or rule out the surprising results from the earlier LHC data. The most recent and precise measurements of the cosmic microwave background from the Planck satellite, from a combination of measurements of the dark matter distribution in the universe, and from a measurement of the expansion rate of the universe today show some disagreement when interpreted in terms of the so-called LambdaCDM model. Schmaltz and collaborators proposed an alternative model to LambdaCDM in which the usual cold dark matter is replaced by a new ``dark sector". This sector consists of a cold dark matter particle which interacts with a newly postulated dark radiation component of the universe. The dark radiation can help explain the discrepancy in measurements of the expansion rate, and the dark matter interactions subtly modify the clumping of dark matter at large scales, thus potentially explaining both kinds of tensions in the data. In two publications Schmaltz described the new model and then performed a precision comparison of the predictions of the model with all currently available cosmological data. The results favor the new model at the level of three standard deviations with current data. Quantum Field Theory (QFT) is the language we use to describe quantum systems which are consistent with Einstein’s theory of Special Relativity. In particular, the requirement of Einstein’s theory that signals not travel faster than the speed of light constrains the types of interactions which particles can engage in. One consequence of relativity is that these interactions cannot preserve particle number. The stronger the interactions, the more severe the particle number violation in a given Relativistic QFT. When particle number violation is strong, it becomes very difficult to adequately parameterize the quantum wave function (which characterizes the state of a quantum system). For example, though we can formulate the QFT which describes the strong force as a set of interactions between quark and gluon particles, we have no clear idea how to express the proton state in terms of these quarks and gluons. This is because the proton, though a bound state of quarks and gluons, is not a state of a fixed number of particles due to strong interactions. Yet, understanding the proton state is very important in order to theoretically predict the reaction rates observed at the LHC in Geneva, which is a proton-proton collider. Katz has formulated a new approach to QFT, which among other things offers a way to adequately approximate the quantum wave function of a bound state at strong coupling. The approximation scheme is related to the fact that any sensible QFT (including that of the strong interactions) is at short distances approximately self-similar upon rescaling of space and time. It turns out that keeping track of the response upon this rescaling is important in efficiently parameterizing the state. Katz and collaborators have used this observation to approximate the state of the proton in toy versions of the strong force. In the late 60s Sheldon Glashow, Abdus Salam and Steven Weinberg (1979 Nobel Prize awardees) proposed a theory unifying weak and electromagnetic interaction which assumed the existence of new particles, the W and Z bosons. The W and Z bosons were eventually detected in high-energy collision in a particle accelerator at CERN, and the recent discovery of the Higgs meson at the Large Hadron Collider (LHC), always at CERN, completed the picture. However, deep theoretical considerations indicate that the theory by Glashow, Weinberg and Salam, often referred to as "the standard model" cannot be the whole story: the existence of new particles and new interactions at yet higher energies is widely anticipated. The experiments at the LHC are looking for these, while theorists, like Brower, Rebbi and collaborators, are investigating models for these new interactions. Working in a large national collaboration with access to the most powerful DOE computers Brower, Rebbi and colleagues have been using calculational techniques, similar to those successfully employed for many years to investigate the interactions among quarks in nucleons, to study theories that can describe the expected "beyond the standard model" (BSM) interactions. Their results, which include also a model for dark matter, have been published in several refereed papers in prestigious journals. Various ideas in topologically interesting field theories predict hypothetical objects such as fractional charges and Majorana excitations. However, such fascinating objects have not been seen in particle physics. Nevertheless, these objects demonstrate possible phenomena that quantum field theory can support. Pi used condensed matter physics as a laboratory to study possible realizations and observable effects of these objects predicted by quantum field theory. In recent times there has developed considerable interest among condensed matter field theorists in precisely the same geometrical and topological structures, which were first discovered in particle physics field theories. From particle physicists' point of view, this is an interesting development, since condensed matter provides an arena in which one can concretely realize particle physics ideas. Moreover, particle physicists can learn new ideas from condensed matter physics. Higgs phenomenon is precisely an important particle physics realization of condensed matter ideas. In contrast to the small distance characterizing condensed matter systems, field theory also describes large distance physics characterizing cosmology. Pi worked on various geometrical effects in the standard theory of cosmology, viz general relativity.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

A second-order semi-implicit δf method for hybrid simulation

A second-order accurate semi-implicit Lorentz force ions, fluid electrons δf hybrid model has been developed using a “current closure” scheme. The model assumes quasi-neutrality and is fully electromagnetic. The implicit field solver improves numerical accuracy by separating the equilibrium terms in the presence of small perturbations. The equilibrium part of the generalized Ohm’s law is solved by direct matrix inversion along the direction of gradients for every Fourier mode in the other two directions, while the nonlinear part is solved iteratively. The simulation has been benchmarked on Alfvén waves, ion sound waves and whistler waves against analytical dispersion relation in a slab. In particular, the first-order and second-order schemes are compared by studying the numerical damping of whistler waves. The full evolution of the resistive tearing mode using the Harris sheet equilibrium is also investigated. The linear growth rate and mode structure are compared with the resistive MHD theory. Important tearing mode nonlinear phenomena such as the Rutherford regime and saturation are demonstrated. We also presented systematic study of Rutherford growth rates and saturation island width, which is consistent with previous MHD studies.

97 MATHEMATICS AND COMPUTING↗

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↗

Effect of interparticle fields and radiation reaction on beam dynamics

The dynamics of relativistic particles in an intense electromagnetic field can be described by the Landau-Lifshitz (LL) equation, where the radiation reaction (RR) is accounted for via a self-force, and interparticle fields are often neglected as an approximation. However, the inclusion of interparticle fields is necessary to ensure energy-momentum conservation, particularly during coherent emission. Here we present (i) an analytical proof showing that the energy-momentum conservation law of the Hamilton-Rohrlich-Dirac action, which is divergence free and describes a generic system of interacting charges, respects causality and provides physically sensible results; (ii) a simple generalization of the LL equation for many particles evaluated as a function of the total field, i.e., the sum of the external and interparticle fields. By performing first-principles numerical simulations of a neutral, relativistic bunch of electrons and positrons (e − /e + ) colliding with a laser pulse, this theory is shown to satisfy energy-momentum conservation when interparticle fields and RR are simultaneously taken into account; and (iii) the combined effect of interparticle fields and RR primarily affects the tail of the particle energy distribution. Additionally, our first-principles simulations show that the effect of interparticle fields on beam energy loss becomes smaller when most of the radiated energy is incoherent.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗