Source author record

Xian-Zhu Tang

Xian-Zhu Tang appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

physics.plasm-ph math.NA Numerical Analysis physics.comp-ph cond-mat.mtrl-sci physics.atom-ph

Catalog footprint

What is connected

12works

6topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2022arXiv

An adaptive scalable fully implicit algorithm based on stabilized finite element for reduced visco-resistive MHD

The magnetohydrodynamics (MHD) equations are continuum models used in the study of a wide range of plasma physics systems, including the evolution of complex plasma dynamics in tokamak disruptions. However, efficient numerical solution methods for MHD are extremely challenging due to disparate time and length scales, strong hyperbolic phenomena, and nonlinearity. Therefore the development of scalable, implicit MHD algorithms and high-resolution adaptive mesh refinement strategies is of considerable importance. In this work, we develop a high-order stabilized finite-element algorithm for the reduced visco-resistive MHD equations based on the MFEM finite element library (mfem.org). The scheme is fully implicit, solved with the Jacobian-free Newton-Krylov (JFNK) method with a physics-based preconditioning strategy. Our preconditioning strategy is a generalization of the physics-based preconditioning methods in [Chacon, et al, JCP 2002] to adaptive, stabilized finite elements. Algebraic multigrid methods are used to invert sub-block operators to achieve scalability. A parallel adaptive mesh refinement scheme with dynamic load-balancing is implemented to efficiently resolve the multi-scale spatial features of the system. Our implementation uses the MFEM framework, which provides arbitrary-order polynomials and flexible adaptive conforming and non-conforming meshes capabilities. Results demonstrate the accuracy, efficiency, and scalability of the implicit scheme in the presence of large scale disparity. The potential of the AMR approach is demonstrated on an island coalescence problem in the high Lundquist-number regime ($\ge 10^7$) with the successful resolution of plasmoid instabilities and thin current sheets.

preprint2022arXiv

Bohm criterion of plasma sheaths away from asymptotic limit

The plasma exit flow speed at the sheath entrance is constrained by the Bohm criterion. The so-called Bohm speed regulates the plasma particle and power exhaust fluxes to the wall, and it is commonly deployed as a boundary condition to exclude the sheath region in quasi-neutral plasma modeling. Here the Bohm criterion analysis is performed in the intermediate plasma regime away from the previously known limiting cases of adiabatic laws and the asymptotic limit of infinitesimal Debye length in a finite-size system, using the transport equations of an anisotropic plasma. The resulting Bohm speed has explicit dependence on local plasma heat flux, temperature isotropization, and thermal force. Comparison with kinetic simulations demonstrates its accuracy over the plasma-sheath transition region in which quasineutrality is weakly perturbed and Bohm criterion applies.

preprint2020arXiv

An adaptive discontinuous Petrov-Galerkin method for the Grad-Shafranov equation

In this work, we propose and develop an arbitrary-order adaptive discontinuous Petrov-Galerkin (DPG) method for the nonlinear Grad-Shafranov equation. An ultraweak formulation of the DPG scheme for the equation is given based on a minimal residual method. The DPG scheme has the advantage of providing more accurate gradients compared to conventional finite element methods, which is desired for numerical solutions to the Grad-Shafranov equation. The numerical scheme is augmented with an adaptive mesh refinement approach, and a criterion based on the residual norm in the minimal residual method is developed to achieve dynamic refinement. Nonlinear solvers for the resulting system are explored and a Picard iteration with Anderson acceleration is found to be efficient to solve the system. Finally, the proposed algorithm is implemented in parallel on MFEM using a domain-decomposition approach, and our implementation is general, supporting arbitrary order of accuracy and general meshes. Numerical results are presented to demonstrate the efficiency and accuracy of the proposed algorithm.

preprint2020arXiv

Neural network representability of fully ionized plasma fluid model closures

The closure problem in fluid modeling is a well-known challenge to modelers aiming to accurately describe their system of interest. Over many years, analytic formulations in a wide range of regimes have been presented but a practical, generalized fluid closure for magnetized plasmas remains an elusive goal. In this study, as a first step towards constructing a novel data based approach to this problem, we apply ever-maturing machine learning methods to assess the capability of neural network architectures to reproduce crucial physics inherent in popular magnetized plasma closures. We find encouraging results, indicating the applicability of neural networks to closure physics but also arrive at recommendations on how one should choose appropriate network architectures for given locality properties dictated by underlying physics of the plasma.

preprint2019arXiv

Impact of a minority relativistic electron tail interacting with a thermal plasma containing high-atomic-number impurities

A minority relativistic electron component can arise in both laboratory and naturally-occurring plasmas. In the presence of high-atomic-number ion species, the ion charge state distribution at low bulk electron temperature can be dominated by relativistic electrons, even though their density is orders of magnitude lower. This is due to the relativistic enhancement of the collisional excitation and ionization cross sections. The resulting charge state effect can dramatically impact the radiative power loss rate and the related Bethe stopping power of relativistic electrons in a dilute plasma.

preprint2016arXiv

Comparison of dust charging between Orbital-Motion-Limited theory and Particle-In-Cell simulations

The Orbital-Motion-Limited (OML) theory has been modified to predict the dust charge and the results were contrasted with the Whipple approximation [Tang and Delzanno, Phys. Plasmas 21, 123708 (2014)]. To further establish its regime of applicability, in this paper the OML predictions (for a non-electron-emitting, spherical dust grain at rest in a collisionless, unmagnetized plasma) are compared with Particle-In-Cell simulations that retain the absorption radius effect. It is found that for large dust grain radius $r_d$ relative to the plasma Debye length $λ_D$, the revised OML theory remains a very good approximation as, for the parameters considered ($r_d/λ_D\le10$, equal electron and ion temperatures), it yields the dust charge to within $20\%$ accuracy. This is a substantial improvement over the Whipple approximation. The dust collected currents and energy fluxes, which remain the same in the revised and standard OML theories, are accurate to within $15-30\%$.

preprint2015arXiv

Orbital-motion-limited theory of dust charging and plasma response

The foundational theory for dusty plasmas is the dust charging theory that provides the dust potential and charge arising from the dust interaction with a plasma. The most widely used dust charging theory for negatively charged dust particles is the so-called orbital motion limited (OML) theory, which predicts the dust potential and heat collection accurately for a variety of applications, but was previously found to be incapable of evaluating the dust charge and plasma response in any situation. Here we report a revised OML formulation that is able to predict the plasma response and hence the dust charge. Numerical solutions of the new OML model show that the widely-used Whipple approximation of dust charge-potential relationship agrees with OML theory in the limit of small dust radius compared with plasma Debye length, but incurs large (order-unity) deviation from the OML prediction when the dust size becomes comparable with or larger than plasma Debye length. This latter case is expected for the important application of dust particles in a tokamak plasma.

preprint2014arXiv

Charging and heat collection by a positively charged dust grain in a plasma

Dust particulates immersed in a quasineutral plasma can emit electrons in several important applications. Once electron emission becomes strong enough, the dust enters the positively charged regime where the conventional Orbital-Motion-Limited (OML) theory can break down due to potential well effects on trapped electrons. A minimal modification of the trapped-passing boundary approximation in the so-called OML$^+$ approach is shown to accurately predict the dust charge and heat collection flux for a wide range of dust size and temperature.

preprint2014arXiv

Reflection and implantation of low energy helium with tungsten surfaces

Reflection and implantation of low energy helium (He) ions by tungsten (W) substrate are studied using molecular dynamics (MD) simulations. Motivated by the ITER divertor design, our study considers a range of W substrate temperatures (300 K, 1000 K, 1500 K), a range of He atom incidence energies ($\le$100 eV) and a range of angles of incidence ($0^{\circ}$-$75^{\circ}$) with respect to substrate normal. The MD simulations quantify the reflection and implantation function, the integrated moments such as the particle/energy reflection coefficients and average implantation depths. Distributions of implantation depths, reflected energy, polar and azimuthal angles of reflection are obtained, as functions of simulation parameters, such as W substrate temperature, polar angle of incidence, the energy of incident He, and the type of W substrate surface. Comparison between the MD simulation results, the results obtained using SRIM simulation package, and the existing experimental and theoretical results is provided.

preprint2013arXiv

Thermo-diffusion in inertially confined plasmas

In a plasma of multiple ion species, thermodynamic forces such as pressure and temperature gradients can drive ion species separation via inter-species diffusion. Unlike its neutral mix counterpart, plasma thermo-diffusion is found comparable to, or even much larger than, baro-diffusion. It is shown that such a strong effect is due to the long-range nature of the Coulomb potential, as opposed to short-range interactions in neutral gases. A special composition of the tritium and 3He fuel is identified to have vanishing net diffusion during adiabatic compression, and hence provides an experimental test in which yield degradation is minimized during ICF implosions.

preprint2012arXiv

Electro-diffusion in a plasma with two ion species

Electric field is a thermodynamic force that can drive collisional inter-ion-species transport in a multicomponent plasma. In an inertial confinement fusion (ICF) capsule, such transport causes fuel ion separation even with a target initially prepared to have equal number densities for the two fuel ion species. Unlike the baro-diffusion driven by ion pressure gradient and the thermo-diffusion driven by ion and electron temperature gradients, electro-diffusion has a critical dependence on the charge-to-mass ratio of the ion species. Specifically, it is shown here that electro-diffusion vanishes if the ion species have the same charge-to-mass ratio. An explicit expression for the electro-diffusion ratio is obtained and used to investigate the relative importance of electro- and baro-diffusion mechanisms. In particular, it is found that electro-diffusion reinforces baro-diffusion in the deuterium and tritium mix, but tends to cancel it in the deuterium and helium-3 mix.

preprint2011arXiv

Bounce-free spherical hydrodynamic implosion

In a bounce-free spherical hydrodynamic implosion, the post-stagnation hot core plasma does not expand against the imploding flow. Such an implosion scheme has the advantage of improving the dwell time of the burning fuel, resulting in a higher fusion burn-up fraction. The existence of bounce-free spherical implosions is demonstrated by explicitly constructing a family of self-similar solutions to the spherically symmetric ideal hydrodynamic equations. When applied to a specific example of plasma liner driven magneto-inertial fusion, the bounce-free solution is found to produce at least a factor of four improvement in dwell time and fusion energy gain.