3D Plasma Research Articles

FLARE: field line analysis and reconstruction for 3D boundary plasma modeling

Abstract The FLARE code is a magnetic mesh generator that is integrated within a suite of tools for the analysis of the magnetic geometry in toroidal fusion devices. A magnetic mesh is constructed from field line segments and permits fast reconstruction of field lines in 3D boundary plasma codes such as EMC3-EIRENE. Both intrinsically non-axisymmetric configurations (stellarators) and those with symmetry breaking perturbations of an axisymmetric equilibrium (tokamaks) are supported. The code itself is written in Modern Fortran with MPI support for parallel computing, and it incorporates object-oriented programming for the definition of the magnetic field and the material surface geometry. Extended derived types for a number of different magnetohydrodynamic equilibrium and plasma response models are implemented. The core element of FLARE is a field line tracer with adaptive step-size control, and this is integrated into tools for the construction of Poincaré maps and invariant manifolds of X-points. A collection of high-level procedures that generate output files for visualization is build on top of that. The analysis modules are build with Python frontends that facilitate customization of tasks and/or scripting of parameter scans.

Modeling of plasma facing component erosion, impurity migration, dust transport and melting processes at JET-ILW

An overview of the modeling approaches, validation methods and recent main results of analysis and modeling activities related to the plasma-surface interaction (PSI) in JET-ILW experiments, including the recent H/D/T campaigns, is presented in this paper. Code applications to JET experiments improve general erosion/migration/retention prediction capabilities as well as various physics extensions, for instance a treatment of dust particles transport and a detailed description of melting and splashing of PFC induced by transient events at JET. 2D plasma edge transport codes like the SOLPS-ITER code as well as PSI codes are key to realistic description of relevant physical processes in power and particle exhaust. Validation of the PSI and edge transport models across JET experiments considering various effects (isotope effects, first wall geometry, including detailed 3D shaping of plasma-facing components, self-sputtering, thermo-forces, physical and chemically assisted physical sputtering formation of W and Be hydrides) is very important for predictive simulations of W and Be erosion and migration in ITER as well as for increasing quantitative credibility of the models. JET also presents a perfect test-bed for the investigation and modeling of melt material dynamics and its splashing and droplet ejection mechanisms. We attribute the second group of processes rather to transient events as for the steady state and, thus, treat those as independent additions outside the interplay with the first group.

MHD-FiT: MHD-based dynamic reconstruction of tokamak plasma configuration

This paper introduces an innovative method for reconstructing 2D magnetic flux contours and plasma parameters of dynamically moving tokamak plasmas. While conventional methods like EFIT, based on the Grad–Shafranov equation, are suitable for plasma equilibria with a single magnetic axis, our approach utilizes the MHD equations and shows promise for tokamak plasmas in motion or containing multiple magnetic axes, which may not strictly adhere to plasma equilibria. By utilizing limited edge magnetic probe measurements, our developed model successfully reconstructs the time evolution of two merging plasma toroids in the TS-6 experiment. A comparison with direct 2D magnetic probe measurements in a low β regime reveals a reconstruction error of approximately 3%.

Impurity modes in two-dimensional strongly coupled complex plasma crystals

Complex plasma fluctuation processes have been extensively studied in many aspects, especially lattice waves in strongly coupled plasma crystals, which are of great significance for understanding fundamental physical phenomena. A challenge of experimental investigations in two-dimensional strongly coupled complex plasma crystals is to keep the main body and foreign particles of different masses on the same horizontal plane. To solve the problem, we have proposed a potential well formed by two negatively biased grids to bind the negatively charged particles in a two-dimensional (2D) plane, thus achieving a 2D plasma crystal in the microgravity environment. The study of such phenomena in complex plasma crystals under microgravity environment then becomes possible. In this paper, we focus on the continuum spectrum, including both phonon and optic branches of the impurity mode in a 2D system in microgravity environments. The results show the dispersion relation of the longitudinal and transverse impurity oscillation modes and their properties. Considering the macroscopic visibility of complex mesoscopic particle lattices, theoretical and experimental studies on this kind of complex plasma systems will help us further understand the physical nature of a wide range of condensed matters.

Growth mechanism of high‐voltage electric pulse rock breaking 3D plasma channel in drilling fluid environment

AbstractHigh‐voltage electric pulse rock breaking has excellent potential for exploiting deep geothermal resources. Numerous researchers have conducted experimental studies on this topic, particularly in rock mechanics, where the breakdown occurs. However, there has been limited scholarly research on drilling fluid. Therefore, the study focuses on the drilling fluid suitable for electric pulse drilling, considering the characteristics of electric pulse rock breaking, which differ from traditional rock breaking. The study focused on the impact of various drilling fluid parameters on the effectiveness of electric impulse rock breaking using red sandstone as the experimental material. This was investigated using the finite element method, and indoor electric rock‐breaking tests were conducted in a drilling fluid environment. The results indicate that the plasma channel mainly grows in the permeable layer of the drilling fluid, resulting in shallow rock breaking depth in the drilling fluid environment. The pore permeated by drilling fluid guides the growth of the plasma channel. The higher the conductivity of the drilling fluid, the closer the ion channel of rock breaking by electric pulse is to the rock surface. This results in a smaller crushing volume and shallower damage depth, which is more detrimental to rock breaking by an electric pulse. The viscosity of drilling fluid can impede the breakdown to some extent. In this paper, the influence of drilling fluid parameters on electro‐pulse rock‐breaking technology is preliminarily studied, which has significant reference value for the selection of actual drilling fluid.

Development of a neural network model for peeling–ballooning stability analysis in the KSTAR tokamak pedestals

The neural network model, MISHKA-NN is developed to mitigate the computational burden associated with the linear ideal magnetohydrodynamic (MHD) stability analysis of the pedestal based on the peeling–ballooning (P–B) model. By utilizing both 1D plasma profiles (current density, pressure gradient, and safety factor) and 0D parameters (plasma geometry, total current, and toroidal mode number), the model predicts linear growth rate of edge-localized ideal MHD instability in a given equilibrium state. By enabling the prediction of each instability within a second, the model reduces the time required for plotting a pedestal P–B stability diagram ( j−α diagram) from approximately 100 CPU hours to a few CPU minutes. Notably, even with the utilization of parametric pressure and current profiles and plasma boundary shapes for the training dataset, the model shows a satisfactory level of performance in benchmarking the j−α diagram for the reconstructed equilibrium from a KSTAR tokamak experiment. We anticipate the model to serve as a versatile alternative to 2D linear MHD stability codes, alleviating numerical costs.

Crossover in densities of confined particles with finite range of interaction

We study a one-dimensional classical system of N particles confined within a harmonic trap. Interactions among these particles are dictated by a pairwise potential V(x), where x is the separation between two particles. Each particle can interact with a maximum of d neighbouring particles on either side (left or right), if available. By adjusting the parameter d, the system can be made nearest neighbour (d=1) to all-to-all (d=N−1) interacting. As suggested by prior studies, the equilibrium density profile of these particles is expected to undergo shape variations as d is changed. In this paper, we investigate this crossover by tuning the parameter f(=d/N) from 1 to 0 in the large N limit for two distinct choices of interaction potentials, V(x)=−|x| and V(x)=−log⁡(|x|) which correspond to 1d one-component plasma and the log-gas model, respectively. For both models, the system size scaling of the density profile for fixed f turns out to be the same as in their respective all-to-all case. However, the scaling function exhibits diverse shapes as f varies. We explicitly compute the average density profile for any f∈(0,1] in the 1d plasma model, while for the log-gas model, we provide approximate calculations for large (close to 1) and small (close to 0) f. Additionally, we present simulation results to numerically demonstrate the crossover and compare these findings with our theoretical results.

Two-Dimensional Electron-Hole Plasma in Colloidal Quantum Shells Enables Integrated Lasing Continuously Tunable in the Red Spectrum.

Combining integrated optical platforms with solution-processable materials offers a clear path toward miniaturized and robust light sources, including lasers. A limiting aspect for red-emitting materials remains the drop in efficiency at high excitation density due to non-radiative quenching pathways, such as Auger recombination. Next to this, lasers based on such materials remain ill characterized, leaving questions about their ultimate performance. Here, we show that colloidal quantum shells (QSs) offer a viable solution for a processable material platform to circumvent these issues. We first show that optical gain in QSs is mediated by a 2D plasma state of unbound electron-hole pairs, opposed to bound excitons, which gives rise to broad-band and sizable gain across the full red spectrum with record gain lifetimes and a low threshold. Moreover, at high excitation density, the emission efficiency of the plasma state does not quench, a feat we can attribute to an increased radiative recombination rate. Finally, QSs are integrated on a silicon nitride platform, enabling high spectral contrast, surface emitting, and TE-polarized lasers with ultranarrow beam divergence across the entire red spectrum from a small surface area. Our results indicate QS materials are an excellent materials platform to realize highly performant and compact on-chip light sources.

Characterization of hydrogen plasma bulk in 2.45 GHz ECR ion source

Electron cyclotron resonance ion source (ECRIS) as the front-end for the neutral beam injection (NBI) system has production potential of hydrogen ion beam in tokamaks. In this article, two methods are hired to evaluate the plasma ion fractions: H+, H2+ and H3+ with respect to the ECRIS operating parameters. The methods are matured by utilizing 28 important volume and surface collisions processes. We encrypted a global 0D numerical model (GM) based on solving the nonlinear equations of plasma species as the first method. Next, a 3D plasma simulation (PS) is applied through the finite element method (FEM) by COMSOL Multiphysics. Then, the results are compared to the experimental data.

Toward plasma drifts in EMC3: Implementation of gradient, divergence, and particle tracing schemes

AbstractThis paper presents a first implementation of gradient, divergence, and particle tracing schemes for the EMC3 code, a stochastic 3D plasma fluid code widely employed for edge plasma and impurity transport modeling in tokamaks and stellarators. These schemes are essential to accommodate plasma drift computations, which are currently absent in the code. Plasma drifts have been recognized to significantly influence transport of particles and energy, and their inclusion in future code upgrades will substantially enhance the code's predictive capabilities. For gradient and divergence calculations, we introduce a second‐order least‐squares gradient scheme. We confirm the second‐order convergence properties and assess the accuracy of several analytical test cases in the presence of synthetic noise. In the second part of this paper, we employ the validated gradient scheme in a fourth‐order Runge–Kutta particle tracing scheme to trace a particle through a generic drift velocity field. The impact of synthetic noise on the scheme's performance is investigated by evaluating various error metrics. We find that the implemented schemes function as intended and exhibit sufficient accuracy to enable drift computations.

Development of a toroidal soft x-ray imaging system and application for investigating three-dimensional plasma on J-TEXT

A toroidal soft x-ray imaging (T-SXRI) system has been developed to investigate three-dimensional (3D) plasma physics on J-TEXT. This T-SXRI system consists of three sets of SXR arrays. Two sets are newly developed and located on the vacuum chamber wall at toroidal positions of ° and °, respectively, while one set was established previously at °. Each set of SXR arrays consists of three arrays viewing the plasma poloidally, and hence can be used separately to obtain SXR images via the tomographic method. The sawtooth precursor oscillations are measured by T-SXRI, and the corresponding images of perturbative SXR signals are successfully reconstructed at these three toroidal positions, hence providing measurement of the 3D structure of precursor oscillations. The observed 3D structure is consistent with the helical structure of the m/n = 1/1 mode. The experimental observation confirms that the T-SXRI system is able to observe 3D structures in the J-TEXT plasma.

Mesenchymal Stem Cells Cultured in a 3D Microgel Environment Containing Platelet-Rich Plasma Significantly Modify Their Chondrogenesis-Related miRNA Expression.

The aim of this work is to study the effect of platelet factors on the differentiation of mesenchymal stem cells (MSCs) to hyaline cartilage chondrocytes in a three-dimensional environment. MSCs were cultured in a microgel environment with a chondrogenic medium. The microgel consisted of microspheres that combine gelatin and platelet-rich plasma (PRP). The gelatin/PRP microdroplets were produced by emulsion. The gelatin containing the microdroplets was enzymatically gelled, retaining PRP and, just before seeding the cells, platelets were activated by adding calcium chloride so that platelet growth factors were released into the culture media but not before. Platelet activation was analyzed before activation to rule out the possibility that the gelatin cross-linking process itself activated the platelets. The gene expression of characteristic chondrogenic markers and miRNA expression were analyzed in cells cultured in a differentiation medium and significant differences were found between gelation/PRP microgels and those containing only pure gelatin. In summary, the gelatin microspheres effectively encapsulated platelets that secreted and released factors that significantly contributed to cellular chondrogenic differentiation. At the same time, the microgel constituted a 3D medium that provided the cells with adherent surfaces and the possibility of three-dimensional cell-cell contact.

Excitation of wakefields in carbon nanotubes: a hydrodynamic model approach

The interactions of charged particles with carbon nanotubes (CNTs) may excite electromagnetic modes in the electron gas produced in the cylindrical graphene shell constituting the nanotube wall. This wake effect has recently been proposed as a potential novel method of short-wavelength high-gradient particle acceleration. In this work, the excitation of these wakefields is studied by means of the linearized hydrodynamic model. In this model, the electronic excitations on the nanotube surface are described treating the electron gas as a 2D plasma with additional contributions to the fluid momentum equation from specific solid-state properties of the gas. General expressions are derived for the excited longitudinal and transverse wakefields. Numerical results are obtained for a charged particle moving within a CNT, paraxially to its axis, showing how the wakefield is affected by parameters such as the particle velocity and its radial position, the nanotube radius, and a friction factor, which can be used as a phenomenological parameter to describe effects from the ionic lattice. Assuming a particle driver propagating on axis at a given velocity, optimal parameters were obtained to maximize the longitudinal wakefield amplitude.

Research on microstructure and properties of 3D plasma spray welding repair layer of continuous casting segmented rollers

In this paper, the Fe-Cr-Ni-Si powder was used in a 3D plasma spray welding machine to perform spray welding repair tests on continuous casting segmented rollers for service life improvement. The phase composition, hardness distribution, abrasion resistance and thermal stability of the spray welding layer were analyzed. Results indicated that the spray welding layer produced by 160 A welding current was more uniform and was with finer organization. The main phase zone was a single γ-Fe and the fusion zone was metallurgically combined with the base metal in a planar grain manner. The hardness of the spray welding layer was the maximum when the welding current was 160 A. The abrasion resistance for each position of the spray welding layer was uniform. The γ-Fe phase zone decreased and the α-Fe phase zone increased with reduction of average hardness after the thermal stability test. In industrial tests, the optimum spray welding parameters were verified and the service life of the continuous casting segmented rollers with spray welding was 2 longer than rollers using surfacing welding process.

MAxSIM: multi-angle-crossing structured illumination microscopy with height-controlled mirror for 3D topological mapping of live cells

Mapping 3D plasma membrane topology in live cells can bring unprecedented insights into cell biology. Widefield-based super-resolution methods such as 3D-structured illumination microscopy (3D-SIM) can achieve twice the axial ( ~ 300 nm) and lateral ( ~ 100 nm) resolution of widefield microscopy in real time in live cells. However, twice-resolution enhancement cannot sufficiently visualize nanoscale fine structures of the plasma membrane. Axial interferometry methods including fluorescence light interference contrast microscopy and its derivatives (e.g., scanning angle interference microscopy) can determine nanoscale axial locations of proteins on and near the plasma membrane. Thus, by combining super-resolution lateral imaging of 2D-SIM with axial interferometry, we developed multi-angle-crossing structured illumination microscopy (MAxSIM) to generate multiple incident angles by fast, optoelectronic creation of diffraction patterns. Axial localization accuracy can be enhanced by placing cells on a bottom glass substrate, locating a custom height-controlled mirror (HCM) at a fixed axial position above the glass substrate, and optimizing the height reconstruction algorithm for noisy experimental data. The HCM also enables imaging of both the apical and basal surfaces of a cell. MAxSIM with HCM offers high-fidelity nanoscale 3D topological mapping of cell plasma membranes with near-real-time ( ~ 0.5 Hz) imaging of live cells and 3D single-molecule tracking.

Numerical studies of plasma edge in W7-X with 3D FINDIF code

Abstract Modelling of the plasma transport for inherently three-dimensional (3D) problems as in stellarators requires dedicated complex codes. FINDIF is a 3D multifluid plasma edge transport code that has been previously successfully used for the analysis of energy transport in the TEXTOR-DED tokamak [1], where 3D perturbations led to an ergodic structure of field lines in the plasma edge. The ongoing efforts to apply it meaningfully to Wendelstein 7-X (W7-X) plasma problems resulted in advancements in the main model and accompanying tools for mesh generation and post-processing. In order to verify the applicability of the code and to compare with the reported simulation (EMC3-EIRENE) and experimental (OP1.1) results, a series of simulations for varying plasma density, temperature and anomalous transport coefficients as well as for fixed input power were performed. The connection length pattern of FINDIF traced magnetic field lines on the limiter was reproduced and its impact on heat loads was confirmed. An increase in the peak heat load on the limiter with a rise in plasma density, temperature and anomalous plasma transport coefficients was observed. The decay lengths of density, electron temperature and heat flux did not change with density, and were decreasing with temperature and increasing with anomalous plasma transport coefficient, which was compared to the simple scrape-off layer (SOL) model.

Simulation of a hollow-cathode PECVD process in O2/TMDSO for silicon dioxide deposition – Cross-code validation of 2D plasma model and global plasma model

We investigate a hollow-cathode plasma system operating in the mixture of O2 and TMDSO using physics- and chemistry-based numerical simulation. Two simulation strategies relevant for hollow-cathode PECVD are developed and cross-correlated – a global plasma model and a 2D plasma model. As a part of this effort, a plasma-kinetic system for O2/TMDSO is proposed, based on an existing system for O2/HMDSO. The two simulation approaches are benchmarked against each other, and their respective advantages and disadvantages are discussed. The models are consequently used to provide insight into the plasma dynamics of the hollow-cathode PECVD system and to understand the plasma chemistry of TMDSO. The paper also identifies the appropriate scaling parameter of the investigated HC-PECVD system, which will be practical for its further development. The results are correlated with experimental observations from literature, finding multiple similarities and consistent behavior.

Tomographic optical emission spectroscopy of an atmospheric pressure plasma jet and surface ionization waves on planar and structured surfaces

In this paper, an approach for 3D plasma structure diagnostics using tomographic optical emission spectroscopy (Tomo-OES) of a nanosecond pulsed atmospheric pressure plasma jet (APPJ) is presented. In contrast to the well-known Abel inversion, Tomo-OES does not require cylindrical symmetry to recover 3D distributions of plasma light emission. Instead, many 2D angular projections are measured with intensified cameras and the multiplicative algebraic reconstruction technique is used to recover the 3D distribution of light emission. This approach solves the line-of-sight integration problem inherent to optical diagnostics, allowing recovery of localized OES information within the plasma that can be used to better infer plasma parameters within complex plasma structures. Here, Tomo-OES was applied to investigate an APPJ operated with helium in ambient air and impinging on planar and structured dielectric surfaces. Surface charging caused the guided streamer from the APPJ to transition to a surface ionization wave (SIW) that propagated along the surface. The SIW experienced variable geometrical and electrical material properties as it propagated, leading to 3D configurations that were non-symmetric and spatially complex. Light emission from He, N, and N2 were imaged at ten angular projections and the respective time-resolved 3D emission distributions in the plasma were then reconstructed. The spatial resolution of each tomographic reconstruction was 7.4 µm and the temporal resolution was 5 ns, sufficient to observe the guided streamer and the effects of the structured surface on the SIW. Emission from He showed the core of the jet and emission from N and N2 indicated effects of entrainment of ambient air. Penning ionization of N2 created a ring or outer layer of N that spatially converged to form the ‘plasma bullet’ or spatially diverged across a surface as part of a SIW. The SIW entered trenches of size 150 µm, leading to decreases in plasma light emission in regions above the trenches. The plasma light emission was higher in some regions with trenches, possibly due to effects of field enhancement.

Energy transfer and scale dynamics in 2D and 3D laser-driven jets

We demonstrate a methodology for diagnosing the multiscale dynamics and energy transfer in complex HED flows with realistic driving and boundary conditions. The approach separates incompressible, compressible, and baropycnal contributions to energy scale-transfer and quantifies the direction of these transfers in (generalized) wavenumber space. We use this to compare the kinetic energy (KE) transfer across scales in simulations of 2D axisymmetric vs fully 3D laser-driven plasma jets. Using the FLASH code, we model a turbulent jet ablated from an aluminum cone target in the configuration outlined by Liao et al. [Phys. Plasmas, 26 032306 (2019)]. We show that, in addition to its well known bias for underestimating hydrodynamic instability growth, 2D modeling suffers from significant spurious energization of the bulk flow by a turbulent upscale cascade. In 2D, this arises as vorticity and strain from instabilities near the jet's leading edge transfer KE upscale, sustaining a coherent circulation that helps propel the axisymmetric jet farther (≈25% by 3.5 ns) and helps keep it collimated. In 3D, the coherent circulation and upscale KE transfer are absent. The methodology presented here may also help with inter-model comparison and validation, including future modeling efforts to alleviate some of the 2D hydrodynamic artifacts highlighted in this study.

Functionalized 3D-Printed PLA Biomimetic Scaffold for Repairing Critical-Size Bone Defects.

The treatment of critical-size bone defects remains a complicated clinical challenge. Recently, bone tissue engineering has emerged as a potential therapeutic approach for defect repair. This study examined the biocompatibility and repair efficacy of hydroxyapatite-mineralized bionic polylactic acid (PLA) scaffolds, which were prepared through a combination of 3D printing technology, plasma modification, collagen coating, and hydroxyapatite mineralization coating techniques. Physicochemical analysis, mechanical testing, and in vitro and animal experiments were conducted to elucidate the impact of structural design and microenvironment on osteogenesis. Results indicated that the PLA scaffold exhibited a porosity of 84.1% and a pore size of 350 μm, and its macrostructure was maintained following functionalization modification. The functionalized scaffold demonstrated favorable hydrophilicity and biocompatibility and promoted cell adhesion, proliferation, and the expression of osteogenic genes such as ALP, OPN, Col-1, OCN, and RUNX2. Moreover, the scaffold was able to effectively repair critical-size bone defects in the rabbit radius, suggesting a novel strategy for the treatment of critical-size bone defects.