Computational study on microstructure evolution and magnetic property of laser additively manufactured magnetic materials

Additive manufacturing (AM) offers an unprecedented opportunity for the quick production of complex shaped parts directly from a powder precursor. But its application to functional materials in general and magnetic materials in particular is still at the very beginning. Here we present the first attempt to computationally study the microstructure evolution and magnetic properties of magnetic materials (e.g. Fe-Ni alloys) processed by selective laser melting (SLM). SLM process induced thermal history and thus the residual stress distribution in Fe-Ni alloys are calculated by finite element analysis (FEA). The evolution and distribution of the $\gamma$-Fe-Ni and FeNi$_3$ phase fractions were predicted by using the temperature information from FEA and the output from CALculation of PHAse Diagrams (CALPHAD). Based on the relation between residual stress and magnetoelastic energy, magnetic properties of SLM processed Fe-Ni alloys (magnetic coercivity, remanent magnetization, and magnetic domain structure) are examined by micromagnetic simulations. The calculated coercivity is found to be in line with the experimentally measured values of SLM-processed Fe-Ni alloys. This computation study demonstrates a feasible approach for the simulation of additively manufactured magnetic materials by integrating FEA, CALPHAD, and micromagnetics.Comment: 20 pages, 15 figure

Numerical study of tearing mode seeding in tokamak X-point plasma

A detailed understanding of island seeding is crucial to avoid (N)TMs and their negative consequences like confinement degradation and disruptions. In the present work, we investigate the growth of 2/1 islands in response to magnetic perturbations. Although we use externally applied perturbations produced by resonant magnetic perturbation (RMP) coils for this study, results are directly transferable to island seeding by other MHD instabilities creating a resonant magnetic field component at the rational surface. Experimental results for 2/1 island penetration from ASDEX Upgrade are presented extending previous studies. Simulations are based on an ASDEX Upgrade L-mode discharge with low collisionality and active RMP coils. Our numerical studies are performed with the 3D, two fluid, non-linear MHD code JOREK. All three phases of mode seeding observed in the experiment are also seen in the simulations: first a weak response phase characterized by large perpendicular electron flow velocities followed by a fast growth of the magnetic island size accompanied by a reduction of the perpendicular electron velocity, and finally the saturation to a fully formed island state with perpendicular electron velocity close to zero. Thresholds for mode penetration are observed in the plasma rotation as well as in the RMP coil current. A hysteresis of the island size and electron perpendicular velocity is observed between the ramping up and down of the RMP amplitude consistent with an analytically predicted bifurcation. The transition from dominant kink/bending to tearing parity during the penetration is investigated

Trans-Debye Scale Plasma Modeling & Stochastic GRB Wakefield Plasma Processes

Modeling plasma physical processes in astrophysical context demands for both detailed kinetics and large scale development of the electromagnetic field densities. We present a new framework for modeling plasma physics of hot tenuous plasmas by a two-split scheme, in which the large scale fields are modeled by means of a particle-in-cell (PIC) code, and in which binary collision processes and single-particle processes are modeled through a Monte-Carlo approach. Our novel simulation tool -- the PhotonPlasma code -- is a unique hybrid model; it combines a highly parallelized (Vlasov) particle-in-cell approach with continuous weighting of particles and a sub-Debye Monte-Carlo binary particle interaction framework. As an illustration of the capabilities we present results from a numerical study of Gamma-Ray Burst - Circumburst Medium interaction and plasma preconditioning via Compton scattering. We argue that important microphysical processes can only viably be investigated by means of hybrid codes such as the PhotonPlasma code. Our first results from 3D simulations with this new simulation tool suggest that magnetic fields and plasma filaments are created in the wakefield of prompt gamma-ray bursts. Furthermore, the photon flux density gradient impacts on particle acceleration in the burst head and wakefield. We discuss some possible implications of the circumburst medium being preconditioned for a trailing afterglow shock front. We also discuss important improvements for future studies of GRB wakefields processes, using the PhotonPlasma code.Comment: 11 pages, 6 figures, to appear in proceedings of "Cool discs, hot flows: The varying faces of accreting compact objects", Sweden, March 2008, AIP Conference Proc. Rev1: Minor changes, in compliance with referee comments. Rev2: corrected some citation formattin

Langevin Simulation of Thermally Activated Magnetization Reversal in Nanoscale Pillars

Numerical solutions of the Landau-Lifshitz-Gilbert micromagnetic model incorporating thermal fluctuations and dipole-dipole interactions (calculated by the Fast Multipole Method) are presented for systems composed of nanoscale iron pillars of dimension 9 nm x 9 nm x 150 nm. Hysteresis loops generated under sinusoidally varying fields are obtained, while the coercive field is estimated to be 1979 $\pm$ 14 Oe using linear field sweeps at T=0 K. Thermal effects are essential to the relaxation of magnetization trapped in a metastable orientation, such as happens after a rapid reversal of an external magnetic field less than the coercive value. The distribution of switching times is compared to a simple analytic theory that describes reversal with nucleation at the ends of the nanomagnets. Results are also presented for arrays of nanomagnets oriented perpendicular to a flat substrate. Even at a separation of 300 nm, where the field from neighboring pillars is only $\sim$ 1 Oe, the interactions have a significant effect on the switching of the magnets.Comment: 19 pages RevTeX, including 12 figures, clarified discussion of numerical technique