Search NASA⌕ Search

SEARCH · Search NASA

Results for “LORENTZ FORCE”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 109 records · Page 6

Revealing Color Forces with Transverse Polarized Electron Scattering

The Spin Asymmetries of the Nucleon Experiment measured two double spin asymmetries using a polarized proton target and polarized electron beam at two beam energies, 4.7 and 5.9 GeV. A large acceptance open-configuration detector package identified scattered electrons at 40° and covered a wide range in Bjorken $\mathcal{x}$ (0.3< $\mathcal{x}$ <0.8). Proportional to an average color Lorentz force, the twist-3 matrix element, $\tilde{d}$$p\atop{2}$, was extracted from the measured asymmetries at Q 2 values ranging from 2.0 to 6.0 GeV 2 . The data display the opposite sign compared to most quark models, including the lattice QCD result, and an unexpected scale dependence. Moreover, when combined with the neutron data in the same Q 2 range the results suggest a flavor independent average color Lorentz force.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Transformation of Spinors in Accelerators and Beam Transfer Lines

The stable spin direction and the spin tune are two important concepts used when polarized beams are transferred along beam lines or circulating in an accelerator. One of the methods to calculate the stable spin direction and the spin tune of the reference particle of a polarized beam bunch is by integrating simultaneously the Lorentz Force $\vec{F}$ = q($\vec{E}$ +$\vec{v}$ x $\vec{B}$) and the BMT equation of the reference particle as it is transported by the magnets of a beam line or accelerator. A modified version of the Raytrace code was used to integrate the Lorentz equation and the BMT equation. Using the spin motion which is governed by the BMT equation one can compute both, the stable spin direction of the reference particle at any point along the reference orbit of an accelerator or a beam line, and the spin tune of the closed orbit in an accelerator. The simultaneous integration of the Lorentz force and of the BMT equation is the most accurate method to calculate the trajectory and the spin motion of a charged particle moving in an electromagnetic field, however this method requires the accurate description of the electromagnetic fields of the elements in which the charged particles are moving. An alternative method of calculating the stable spin direction and the spin tune of spin ½ charged particles along the reference trajectory of a beam line or an accelerator is the transformation of the spinor which is the wavefunction of spin ½ particles. This method requires the knowledge of the stable spin direction of each element of the beam line or the accelerator. This technical note describes the “spinor-method” to calculate the stable-spin- direction and the spin-tune of a spin ½ particle moving along the reference trajectory of a beam line or an accelerator.

43 PARTICLE ACCELERATORS↗

Lorentz drift compensation in high harmonic generation in the soft and hard X-ray regions of the spectrum

Here, we present a semi-classical study of the effects of the Lorentz force on electrons during high harmonic generation in the soft and hard X-ray regions driven by near- and mid-infrared lasers with wavelengths from 0.8 to 20 μm, and at intensities below 10 15 W/cm 2 . The transverse extent of the longitudinal Lorentz drift is compared for both Gaussian focus and waveguide geometries. Both geometries exhibit a longitudinal electric field component that cancels the magnetic Lorentz drift in some regions of the focus, once each full optical cycle. We show that the Lorentz force contributes a super-Gaussian scaling which acts in addition to the dominant high harmonic flux scaling of λ -(5-6) due to quantum diffusion. We predict that the high harmonic yield will be reduced for driving wavelengths > 6 μm, and that the presence of dynamic spatial mode asymmetries results in the generation of both even and odd harmonic orders. Remarkably, we show that under realistic conditions, the recollision process can be controlled and does not shut off completely even for wavelengths >10 μm and recollision energies greater than 15 keV.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Comment on “Symplectic integration of magnetic systems”: A proof that the Boris algorithm is not variational

One popular technique for the numerical time advance of charged particles interacting with electric and magnetic fields according to the Lorentz force law [1], [2], [3] and [4] is the Boris algorithm. Its popularity stems from simple implementation, rapid iteration, and excellent long-term numerical fidelity [1] and [5]. Excellent long-term behavior strongly suggests the numerical dynamics exhibit conservation laws analogous to those governing the continuous Lorentz force system [6]. Moreover, without conserved quantities to constrain the numerical dynamics, algorithms typically dissipate or accumulate important observables such as energy and momentum over long periods of simulated time [6]. Identification of the conservative properties of an algorithm is important for establishing rigorous expectations on the long-term behavior; energy-preserving, symplectic, and volume-preserving methods each have particular implications for the qualitative numerical behavior [6], [7], [8], [9], [10] and [11].

97 MATHEMATICS AND COMPUTING↗

Novel Applications of Magnetic Fields for Fluid Flow Control and for Simulating Variable Gravity Conditions

Static and dynamic magnetic fields have been used to control convection in many materials processing applications. In most of the applications, convection control (damping or enhancement) is achieved through the Lorentz force that can be tailored to counteract/assist dominant system flows. This technique has been successfully applied to liquids that are electrically conducting, such as high temperature melts of semiconductors, metals and alloys, etc. In liquids with low electrical conductivity such as ionic solutions of salts in water, the Lorentz force is weak and hence not very effective and alternate ways of flow control are necessary. If the salt in solution is paramagnetic then the variation of magnetic susceptibility with temperature and/or concentration can be used for flow control. For thermal buoyancy driven flows this can be accomplished in a temperature range below the Curie point of the salt. The magnetic force is proportional to the magnetic susceptibility and the product of the magnetic field and its gradient. By suitably positioning the experiment cell in the magnet, system flows can be assisted or countered, as desired. A similar approach can be extended to diamagnetic substances and fluids but the required magnetic force is considerably larger than that required for paramagnetic substances. The presentation will provide an overview of work to date on a NASA fluid physics sponsored project that aims to test the hypothesis of convective flow control using strong magnetic fields in protein crystal growth. The objective is to understand the nature of the various forces that come into play, delineate causative factors for fluid flow and to quantify them through experiments, analysis, and numerical modeling. The seminar will report specifically on the experimental results using paramagnetic salts and solutions in magnetic fields and compare them to analytical predictions. Applications of the concept to protein crystallization studies will be discussed. The use of strong magnetic fields for terrestrially simulating variable gravity environments and applications supporting the NASA Exploration Initiative will also be briefly discussed.

Ramachandran, N.↗

Canted–Cosine–Theta Magnet (CCT)—A Concept for High Field Accelerator Magnets

Canted-Cosine-Theta (CCT) magnet is an accelerator magnet that superposes fields of nested and tilted solenoids that are oppositely canted. The current distribution of any canted layer generates a pure harmonic field as well as a solenoid field that can be cancelled with a similar but oppositely canted layer. The concept places windings within mandrel's ribs and spars that simultaneously intercept and guide Lorentz forces of each turn to prevent stress accumulation. With respect to other designs, the need for pre-stress in this concept is reduced by an order of magnitude making it highly compatible with the use of strain sensitive superconductors such as Nb 3 Sn or HTS. Intercepting large Lorentz forces is of particular interest in magnets with large bores and high field accelerator magnets like the one foreseen in the future high energy upgrade of the LHC. This paper describes the CCT concept and reports on the construction of CCT1 a "proof of principle" dipole.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Dynamical evolution of dust particles in the Kuiper disk

A Kuiper belt dust disk will have a resonant structure, arising because the Plutinos are in the 3:2 mean motion resonance with Neptune. We run numerical integrations of particles originating from Plutinos to determine what percentage of particles remain in the resonance for a variety of particle and source body sizes. The dynamical evolution of the particles is followed from source to sink with Poynting-Robertson light drag, solar wind drag, radiation pressure, the Lorentz force, neutral interstellar gas drag, and the effects of planetary gravitational perturbations included. The number of particles in the 3:2 resonance increases with decreasing p for the cases where the initial source bodies are small and the percentage of particles in resonance is not significantly changed by either the addition of the Lorentz force, as long as the potential of the particles is small (U = 5 V), or the effect of neutral interstellar gas drag.

Kuiper↗

Magnetic Susceptibility Effects and Lorentz Damping in Diamagnetic Fluids

A great number of crystals (semi-conductor and protein) grown in space are plagued by convective motions which contribute to structural flaws. The character of these instabilities is not well understood but is associated with density variations in the presence of residual gravity and g-jitter. Both static and dynamic (rotating or travelling wave) magnetic fields can be used to reduce the effects of convection in materials processing. In semi-conductor melts, due to their relatively high electrical conductivity, the induced Lorentz force can be effectively used to curtail convective effects. In melts/solutions with reduced electrical conductivity, such as aqueous solutions used in solution crystal growth, protein crystal growth and/or model fluid experiments for simulating melt growth, however, the variation of the magnetic susceptibility with temperature and/or concentration can be utilized to better damp fluid convection than the Lorentz force method. This paper presents a comprehensive, comparative numerical study of the relative damping effects using static magnetic fields and gradients in a simple geometry subjected to a thermal gradient. The governing equations are formulated in general terms and then simplified for the numerical calculations. Operational regimes, based on the best damping technique for different melts/solutions are identified based on fluid properties. Comparisons are provided between the numerical results and available results from experiments in surveyed literature.

Ramachandran, Narayanan↗

Demonstrating electromagnetic control of free-surface, liquid-metal flows relevant to fusion reactors

Plasma-facing components (PFC's) made from solid materials may not be able to withstand the large heat and particle fluxes that will be produced within next-generation fusion reactors. To address the shortcomings of solid PFC's, a variety of liquid-metal (LM) PFC concepts have been proposed. Many of the suggested LM-PFC designs rely on electromagnetic restraint (Lorentz force) to keep free-surface, liquid-metal flows adhered to the interior surfaces of a fusion reactor. However, there is very little, if any, experimental data demonstrating that free-surface, LM-PFC's can actually be electromagnetically controlled. Therefore, in this study, electrical currents were injected into a free-surface liquid-metal that was flowing through a uniform magnetic field. The resultant Lorentz force generated within the liquid-metal affected the velocity and depth of the flow in a controllable manner that closely matched theoretical predictions. Furthermore, these results show the promise of electromagnetic control for LM-PFC's and suggest that electromagnetic control could be further developed to adjust liquid-metal nozzle output, prevent splashing within a tokamak, and alter heat transfer properties for a wide-range of liquid-metal systems.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Demonstrating electromagnetic control of free-surface, liquid-metal flows relevant to fusion reactors

Plasma-facing components (PFC's) made from solid materials may not be able to withstand the large heat and particle fluxes that will be produced within next-generation fusion reactors. To address the shortcomings of solid PFC's, a variety of liquid-metal (LM) PFC concepts have been proposed. Many of the suggested LM-PFC designs rely on electromagnetic restraint (Lorentz force) to keep free-surface, liquid-metal flows adhered to the interior surfaces of a fusion reactor. However, there is very little, if any, experimental data demonstrating that free-surface, LM-PFC's can actually be electromagnetically controlled. Therefore, in this study, electrical currents were injected into a free-surface liquid-metal that was flowing through a uniform magnetic field. The resultant Lorentz force generated within the liquid-metal affected the velocity and depth of the flow in a controllable manner that closely matched theoretical predictions. These results show the promise of electromagnetic control for LM-PFC's and suggest that electromagnetic control could be further developed to adjust liquid-metal nozzle output, prevent splashing within a tokamak, and alter heat transfer properties for a wide-range of liquid-metal systems.

Plasma-Facing Component↗

Magnetohydrodynamics (MHD) Aerocapture System for Enabling Faster-Larger Planetary Science & Human Exploration Missions

Since our completing the NIAC Phase I NIAC on this Advanced Aerocapture System, NASA Langley Research Center has funded or supported a number of studies and code enhancements through its Center Innovation Fund (CIF) and NASA’s NSTGRO and Internship Programs to mature the analysis capabilities and quantify the merits of the MHD Aerocapture System technology. These efforts have resulted in a plug and play analysis capability for assessing MHD aerocapture system performance for arrival at many planetary bodies of interest. Our efforts have especially focused on the potential mass savings for improving the capacity for science observations at Neptune and Triton. A re-cent Forbes article published “‘Orbital mechanics is probably going to decide for us whether we go to Uranus or Neptune because we need to flyby Jupiter,’ said Kunio Sayanagi at Hampton University, Virginia, who also worked on the Neptune Odyssey proposal…. Exactly when a mission can be sent to Uranus, or Neptune, depends on the relative position of Jupiter, which can help give a spacecraft a gravitational slingshot. That drastically shortens the cruise phase.” [1] Since shortening the cruise phase is important for these science missions, any mass savings enabled by the MHD Aerocapture System could be reallocated to increasing Thermal Protection System mass to allow faster arrival speeds and/or for onboarding additional payloads such as science instruments, batteries, or propellant for conducting more science for longer durations in the desired orbits. The analysis steps and codes for conducting trades and sizing vehicles for aerocapture are as follows: Step 1 is to conduct aeroheating analysis using LAURA of the selected entry vehicle shape to identify locations on the forebody where ionization and flow velocity are sufficient for producing Lo-rentz forces. LAURA is a multiblock structured grid finite-volume CFD solver developed at the NASA Langley Research Center. [2] LAURA has been used for aerothermal analysis support of the entry, de-scent and landing (EDL) phase of interplanetary missions over the last three decades [3-7]. Step 2 is to port the LAURA results into CFDWARP to calcu-late electrical and thermal conductivities of ionized flow for sizing MHD patch system and calculating Lorentz forces needed for controls analysis. CFDWARP is a CFD code that uses advanced nu-merical methods that enable the simulation of the full coupling between the aerodynamics, the magne-tohydrodynamics, and the non-neutral plasma sheaths. CFDWARP has the unique capability to simulate efficiently the non-neutral sheaths (near the electrodes) in coupled form with the quasi-neutral bulk MHD flow [8-11]. Step 3 is to link re-sults from LAURA and CFDWARP into POST2 for calculating entry trajectories and comparing MHD control results with other aerodynamic control strategies. The Program to Optimize Simulated Tra-jectories II (POST2) is a generalized point mass, discrete parameter targeting and optimization pro-gram. POST2 provides the capability to target and optimize point mass trajectories for multiple pow-ered or un-powered vehicles near an arbitrary rotat-ing, oblate planet [12]. Step 4: TPS sizing was per-formed using the Fully Implicit Ablation and Ther-mal-response code (FIAT) tool which computes the transient one-dimensional thermal response and surface thermochemistry of a multilayer stackup of thermal protection, bonding, and structural materi-als subject to aeroheating on one surface [13]. The sizing and margining methodology used was based on the approach documented by Mahzari and Milos [14] for the dual-layered heatshield for extreme entry environment technology (DL-HEEET) TPS concept. TPS analysis utilizes trajectory information from POST2. Using this step-wise plug and play MHD Aerocapture performance assessment process, our analysis targets a Neptune aerocapture trajectory that will place the spacecraft in an observation orbit for Triton. [15]. Magnetohydrodynamic (MHD) control of a 4.5-meter diameter MSL-style capsule resulted in TPS mass savings of nearly 2000 kg when using an MHD system mass of under 200 kg. The flight path for a vehicle using the MHD control strategy has a much lower heat rate and heat load compared to the conventional aerodynamic aerocapture strategies known as bank angle con-trolled (BAC) and direct force controlled (DFC). Both BAC and DFC have heat rates significantly greater than 1500 W/cm2 typically used as an upper limit for PICA. Thus, DL-HEEET TPS concept was required for the BAC and DFC control strategies. However, considering the more benign environ-ments for the MHD case, additional TPS concepts with improved mass efficiency were also assessed. PICA was considered for the MHD controlled strat-egy since the maximum heat rate was well within the limits (<1500 W/cm2) of PICA. TPS sizing re-sulted in a significant mass reduction. The PICA layer for this sizing case was about 7.8 cm. As a point of reference, the Mars 2020 mission, which used this same PICA concept, had a PICA thickness of 3.18 cm [16]. The trajectories used for the TPS sizing originat-ed from the POST2 simulations. The current, I, to an electromagnet configuration can be manipulated to allow for active control of the vehicle. Manipula-tion of the current, I, changes the magnetic field, B, which affects the Lorentz force and therefore the MHD drag force on the vehicle. Our analysis in-cluded both open-loop and close-loop control. Closed-loop control will enable improved overall performance when taking into account mission level uncertainties, such as interplanetary delivery errors and atmospheric modeling uncertainties. The open-loop and closed-loop MHD control cases do not dip as deep into the atmosphere as the aerodynamic cases. Three types of aerodynamic-only approach-es are investigated: bank angle modulation (BAM), director force control (DFC), and Drag Modulated. BAM and DFC make use of vehicle aerodynamic angles to steer the vehicle. Thus, changing the aer-odynamic forces acting on the vehicle for control, aerodynamic drag modulated case requires a vary-ing drag area to modulate the drag force. The MHD drag modulated case modulates MHD generated drag force that adds to the aerodynamic drag. This higher atmospheric activation of drag forces by the MHD patch results in significantly less heat flux on the vehicle. The MHD technology will enable shorter cruise times and deceleration of larger payloads for increasing the capacity for science at the Ice Giants or for returning astronauts to Earth from cislunar space or from Mars. The purpose of this presentation is to provide more details about this work and to highlight plans for further research and development including a flight demonstration.

R. W. Moses↗

A Computational Procedure for Assessing I$_{c}(\varepsilon$) in Nb$_{3}$Sn/Bi-2212 Hybrid Magnets

The critical current of superconductors is commonly measured by testing unloaded wires under an external magnetic field. While stressed by intense Lorentz forces, the existing HTS/LTS superconductors are prone to a reduction in critical current before reaching their structural mechanical limit. Here, in this work, the magnetic and mechanical analysis of the FNAL 4-layer Bi-2212/Nb$_{3}$Sn hybrid dipole magnet is reported, aimed at predicting the critical current degradation for both the superconductors during powering at 16T. All the Rutherford cables in the coils of the hybrid magnet were modeled at the strand level in Ansys APDL with the heterogeneous cable model. Utilizing this detailed geometry, it was possible to evaluate the effects of strain on the critical current degradation for both the Nb$_{3}$Sn and Bi-2212 superconductors under the intense Lorentz forces. The analysis presented in this paper integrates strain-dependent critical current laws, with parameters derived from experimental data, to simulate the hybrid magnet's performance for all possible current-powering configurations. The proposed methodology enables a detailed assessment of conductor integrity and I$_{C}(\varepsilon$) reduction in existing hybrid magnet designs, providing a versatile and rigorous framework for optimizing future high-field hybrid magnets.

D'Agliano, A. [Lawrence Berkeley National Laborato↗

Torus CLAS12-Superconducting Magnet Quench Analysis

The JLAB Torus magnet system consists of six superconducting trapezoidal racetrack-type coils assembled in a toroidal configuration. These coils are wound with SSC-36 Nb-Ti superconductor and have the peak magnetic field of 3.6 T. The first coil manufacturing based on the JLAB design began at FNAL. The large magnet system dimensions (8 m diameter and 14 MJ of stored energy) dictate the need for quench protection. Each coil is placed in an aluminum case mounted inside a cryostat and cooled by 4.6 K supercritical helium gas flowing through a copper tube attached to the coil ID. The large coil dimensions and small cryostat thickness drove the design to challenging technical solutions, suggesting that Lorentz forces due to transport currents and eddy currents during quench and various failure scenarios are analyzed. The paper covers the magnet system quench analysis using the OPERA3d Quench code.

Kashikhin, V S↗

Static current-sheet models of quiescent prominences

A particular class of theoretical models idealize the prominence to be a discrete flat electric-current sheet suspended vertically in a potential magnetic field. The weight of the prominence is supported by the Lorentz force in the current sheet. These models can be extended to have curved electric-current sheets and to vary three-dimensionally. The equation for force balance is 1 over 4 pi (del times B) times Bdel p- p9 z=zero. Using Cartesian coordinates we take, for simplicity, a uniform gravity with constant acceleration g in the direction -z. If we are interested not in the detailed internal structure of the prominence, but in the global magnetic configuration around the prominence, we may take prominence plasma to be cold. Consideration is given to how such equilibrium states can be constructed. To simplify the mathematical problem, suppose there is no electric current in the atmosphere except for the discrete currents in the cold prominence sheet. Let us take the plane z =0 to be the base of the atmosphere and restrict our attention to the domain z greater than 0. The task we have is to solve for a magnetic field which is everywhere potential except on some free surface S, subject to suit able to boundary conditions. The surface S is determined by requiring that it possesses a discrete electric current density such that the Lorentz force on it is everywhere vertically upward to balance the weight of the material m(S). Since the magnetic field is potential in the external atmosphere, the latter is decoupled from the magnetic field and its plane parallel hydrostatic pressure and density can be prescribed.

Wu, F.↗

Test Results of the CERN HL-LHC low- Quadrupole Short Models MQXFS3c and MQXFS4

For the High Luminosity Upgrade of the CERN Large Hadron Collider, lower * quadrupole magnets based on advanced Nb3Sn conductors will be installed on each side of the ATLAS and CMS experiment insertion zones. As part of the technological developments needed to achieve the required field gradient of 132.6 T/m within a 150-mm aperture, short length model magnets, named MQXFS, are tested both at the CERN SM18 and FERMI Lab test facilities. The model magnets rely on two types of Nb3Sn conductors (RRP and PIT) and on an innovative bladders and keys design to provide mechanical support against the Lorentz forces. In 2016 and 2017, the powering tests of the first two models MQXFS3 (RRP) and MQXFS5 (PIT) proved that nominal performance (16.5 kA) could be reached with excellent memory of the quench current after thermal cycle. However both magnets showed a slow training behavior with clear observations of voltage disturbances before the quench. Besides MQXFS5 only could reach ultimate current (17.9 kA) whereas erratic behavior was observed on MQXFS3 due to conductor local degradation at the head of one of the coil. In 2018, this limiting coil was changed and the applied azimuthal pre-stress increased. If ultimate current could then be reached, no stable current could be maintained due to identified defect on the outer layer of the new coil. Finally the outcome of the test of the new model MQXFS4, featuring the final RRP conductors that will be used for the series production and variation on the inner layer quench heater designs are here reported in details.

43 PARTICLE ACCELERATORS↗

Magnetoacoustic tomography with magnetic induction for high-resolution bioimepedance imaging through vector source reconstruction under the static field of MRI magnet

Purpose: Magnetoacoustic tomography with magnetic induction (MAT-MI) is an imaging modality to reconstruct the electrical conductivity of biological tissue based on the acoustic measurements of Lorentz force induced tissue vibration. This study presents the feasibility of the authors' new MAT-MI system and vector source imaging algorithm to perform a complete reconstruction of the conductivity distribution of real biological tissues with ultrasound spatial resolution. Methods: In the present study, using ultrasound beamformation, imaging point spread functions are designed to reconstruct the induced vector source in the object which is used to estimate the object conductivity distribution. Both numerical studies and phantom experiments are performed to demonstrate the merits of the proposed method. Also, through the numerical simulations, the full width half maximum of the imaging point spread function is calculated to estimate of the spatial resolution. The tissue phantom experiments are performed with a MAT-MI imaging system in the static field of a 9.4 T magnetic resonance imaging magnet. Results: The image reconstruction through vector beamformation in the numerical and experimental studies gives a reliable estimate of the conductivity distribution in the object with a ∼1.5 mm spatial resolution corresponding to the imaging system frequency of 500 kHz ultrasound. In addition, the experiment results suggest that MAT-MI under high static magnetic field environment is able to reconstruct images of tissue-mimicking gel phantoms and real tissue samples with reliable conductivity contrast. Conclusions: The results demonstrate that MAT-MI is able to image the electrical conductivity properties of biological tissues with better than 2 mm spatial resolution at 500 kHz, and the imaging with MAT-MI under a high static magnetic field environment is able to provide improved imaging contrast for biological tissue conductivity reconstruction.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

An Adaptive Newton-Based Free-Boundary Grad–Shafranov Solver

Equilibria in magnetic confinement devices result from force balancing between the Lorentz force and the plasma pressure gradient. In an axisymmetric configuration like a tokamak, such an equilibrium is described by an elliptic equation for the poloidal magnetic flux, commonly known as the Grad–Shafranov equation. It is challenging to develop a scalable and accurate free-boundary Grad–Shafranov solver, since it is a fully nonlinear optimization problem that simultaneously solves for the magnetic field coil current outside the plasma to control the plasma shape. In this work, we develop a Newton-based free-boundary Grad–Shafranov solver using adaptive finite elements and preconditioning strategies. The free-boundary interaction leads to the evaluation of a domain-dependent nonlinear form of which its contribution to the Jacobian matrix is achieved through shape calculus. The optimization problem aims to minimize the distance between the plasma boundary and specified control points while satisfying two nontrivial constraints, which correspond to the nonlinear finite element discretization of the Grad–Shafranov equation and a constraint on the total plasma current involving a nonlocal coupling term. The linear system is solved by a block factorization, and AMG is called for subblock elliptic operators. The unique contributions of this work include the treatment of a global constraint, preconditioning strategies, nonlocal reformulation, and the implementation of adaptive finite elements. Furthermore, it is found that the resulting Newton solver is robust, successfully reducing the nonlinear residual to 1e-6 and lower in a small handful of iterations while addressing the challenging case to find a Taylor state equilibrium where conventional Picard-based solvers fail to converge.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

EMPHASIS/Nevada UTDEM user guide. Version 2.0.

The Unstructured Time-Domain ElectroMagnetics (UTDEM) portion of the EMPHASIS suite solves Maxwell's equations using finite-element techniques on unstructured meshes. This document provides user-specific information to facilitate the use of the code for applications of interest. UTDEM is a general-purpose code for solving Maxwell's equations on arbitrary, unstructured tetrahedral meshes. The geometries and the meshes thereof are limited only by the patience of the user in meshing and by the available computing resources for the solution. UTDEM solves Maxwell's equations using finite-element method (FEM) techniques on tetrahedral elements using vector, edge-conforming basis functions. EMPHASIS/Nevada Unstructured Time-Domain ElectroMagnetic Particle-In-Cell (UTDEM PIC) is a superset of the capabilities found in UTDEM. It adds the capability to simulate systems in which the effects of free charge are important and need to be treated in a self-consistent manner. This is done by integrating the equations of motion for macroparticles (a macroparticle is an object that represents a large number of real physical particles, all with the same position and momentum) being accelerated by the electromagnetic forces upon the particle (Lorentz force). The motion of these particles results in a current, which is a source for the fields in Maxwell's equations.

99 GENERAL AND MISCELLANEOUS↗