Engineering topics
Zhu, B.
Publications and source records attributed to Zhu, B..
The Compelling Need for a Mid-Scale Stellarator Facility
In the pursuit of the goal of commercial fusion as an abundant and safe source of energy, the stellarator is a leading concept with compelling attractiveness and demonstrated performance. In this white paper we, as a community of US researchers from Universities, National Laboratories, and Private Industry, involved in studying the stellarator concept, lay out the programmatic and technical motivation for a new mid-size stellarator research facility in the US. This contribution is complementary to several other whitepapers authored by members of our community which address different mid-scale stellarator aspects. A community based technical facility proposal has been prepared by F. Parra, et al: Flexible Stellarator Physics Facility. Two private stellarator companies have submitted proposals supporting the development of a mid-scale stellarator: Thea Energy (C.P.S. Swanson, et al.), Type One Energy (W. Guttenfelder, et al.).
The Compelling Need for a Mid-Scale Stellarator Facility
In the pursuit of the goal of commercial fusion as an abundant and safe source of energy, the stellarator is a leading concept with compelling attractiveness and demonstrated performance. A new mid-size stellarator is needed to retire risks and innovate towards a high performance, economically attractive, stellarator Fusion Pilot Plant. In this presentation we, as a community of US researchers from Universities, National Laboratories, and Private Industry, involved in studying the stellarator concept, lay out the programmatic and technical motivation for a new and modern mid-size stellarator research facility. A new mid-scale stellarator is needed to realize the potential predicted by a solid body of theory and simulation along with advances in computational tools for optimization and non-linear turbulence modeling. Notably, it is possible to combine the advantages of the stellarator (steady state, no current drive, no disruptions) with the good confinement regularly achieved in tokamaks. The top priorities for experimental work, and the motivation for a mid-scale stellarator experiment are: turbulence control, non-resonant divertor, MHD stability at large beta, confinement of fast particles, and coil simplification. A new mid-size quasi-symmetric stellarator, built as a user facility, would complement existing research at Wendelstein 7-X and LHD and strongly augment private industry. It would provide a program of innovative research, concept validation, theoretical advancement, and workforce development. Growing support and interest for stellarators by the fusion community and private industry affirms this rationale.
Preface for the Sherwood Fusion Theory 2021–2022 special collection
The Sherwood Fusion Theory conference has a 60-year history and has been held nearly annually for at least five decades. With theoretical and computational research directed toward fusion energy as an organizing principle, the typical Sherwood program encompasses a wide range of plasma physics topics pertaining to a variety of confinement concepts along with advancements in algorithms and numerical methods for fusion plasma simulations. Additionally, the typical annual meeting comprises roughly a dozen invited talks plus a small number of plenary speakers and several poster sessions. The program committee presents several awards each year to graduate (or undergraduate) students with outstanding poster presentations. Beginning in 2021, both invited speakers and student poster award winners were asked to submit manuscripts for this special collection, which combines contributions from the 2021 and 2022 Sherwood meetings.
Singular value decomposition and similarity renormalization group evolution of nuclear interactions
One of the main challenges for ab initio nuclear many-body theory is the growth of computational and storage costs as calculations are extended to heavy, exotic, and structurally complex nuclei. Here, we investigate the factorization of nuclear interactions as a means to address this issue. We perform Singular Value Decompositions of nucleon-nucleon interactions in partial wave representation and study the dependence of the singular value spectrum on interaction characteristics like regularization scheme and resolution scales. We develop and implement the Similarity Renormalization Group (SRG) evolution of the factorized interaction, and demonstrate that this SVD-SRG approach accurately preserves two-nucleon observables. We find that low-resolution interactions allow the truncation of the SVD at low rank, and that a small number of relevant components is sufficient to capture the nuclear interaction and perform an accurate SRG evolution, while the Coulomb interaction requires special consideration. The rank is uniform across all partial waves, and almost independent of the basis choice in the tested cases. This suggests an interpretation of the relevant singular components as mere representations of a small set of abstract operators that can describe the interaction and its SRG flow. Following the traditional workflow for nuclear interactions, we discuss how the transformation between the center-of-mass and laboratory frames creates redundant copies of the partial wave components when implemented in matrix representation, and we discuss strategies for mitigation. Lastly, we test the low-rank approximation to the SRG-evolved interactions in many-body calculations using the In-Medium SRG. By including nuclear radii in our analysis, we verify that the implementation of the SRG using the singular vectors of the interaction does not spoil the evolution of other observables.
Boosting the signal-to-noise of low-field MRI with deep learning image reconstruction
Recent years have seen a resurgence of interest in inexpensive low magnetic field (< 0.3 T) MRI systems mainly due to advances in magnet, coil and gradient set designs. Most of these advances have focused on improving hardware and signal acquisition strategies, and far less on the use of advanced image reconstruction methods to improve attainable image quality at low field. We describe here the use of our end-to-end deep neural network approach (AUTOMAP) to improve the image quality of highly noise-corrupted low-field MRI data. We compare the performance of this approach to two additional state-of-the-art denoising pipelines. We find that AUTOMAP improves image reconstruction of data acquired on two very different low-field MRI systems: human brain data acquired at 6.5 mT, and plant root data acquired at 47 mT, demonstrating SNR gains above Fourier reconstruction by factors of 1.5- to 4.5-fold, and 3-fold, respectively. In these applications, AUTOMAP outperformed two different contemporary image-based denoising algorithms, and suppressed noise-like spike artifacts in the reconstructed images. The impact of domain-specific training corpora on the reconstruction performance is discussed. The AUTOMAP approach to image reconstruction will enable significant image quality improvements at low-field, especially in highly noise-corrupted environments.
Electric field and turbulence in global Braginskii simulations across the ASDEX Upgrade edge and scrape-off layer
Turbulence simulations in diverted geometry across the edge and scrape-off layer (SOL) of ASDEX Upgrade are performed with the GRILLIX code (Stegmeir et al 2019 Phys. Plasmas 26 052517). The underlying global (full-f) drift-reduced Braginskii model allows to concurrently study the self-consistent dynamics of the turbulence and the background as well as the evolution of toroidal and zonal flows. Different contributions to the radial electric field are identified. The dominant contribution on closed flux surfaces comes from the ion pressure gradient, due to the diamagnetic drift in the curved magnetic field. Large deviations can be induced, in particular, by the polarization particle flux, leading to zonal flows. The latter are driven by small-scale eddies, but do not exhibit much impact on the overall transport which is driven by ballooning modes at larger scales. Ion viscosity is found to be important in damping poloidal rotation through adjusting of the parallel velocity profile, but not via direct vorticity damping. The zonal flow drive peaks at the separatrix, where a strong shear layer forms due to the sheath-induced counter-propagating SOL flow, allowing for the formation of a transport barrier. The temperature profile across the separatrix is determined by the competition between cross-field transport and outflow in the SOL, the latter being largely controlled by the parallel heat conductivity.
Pair plasma instability in homogeneous magnetic guide fields
Pair plasmas, collections of both matter and antimatter particles of equal mass, represent a paradigm for the study of basic plasma science, and many open questions exist regarding these unique systems. They are found in many astrophysical settings, such as gamma-ray bursts, and have recently also been produced in carefully designed laboratory experiments. A central research topic in plasma physics is instability; however, unlike their more common ion–electron siblings, pair plasmas are generally thought to be stable to cross field pressure gradients in homogeneous magnetic fields. It is shown here by means of kinetic full-f simulations that, when a pressure gradient is first established, the Gradient-driven Drift Coupling mode is destabilized and becomes turbulent. Force balance is eventually achieved by a combination of flattened pressure profiles due to turbulent transport and establishment of a magnetic field gradient, saturating the growth. During the unstable phase, key physics can be captured by a δf gyrokinetic description, where it is shown analytically and numerically that parallel particle motion results in a coupling of all electromagnetic field components. A fluid model derived therefrom accurately predicts linear eigenmodes and is used to resolve global profile effects. For laser-based electron–positron plasma experiments, prompt instability is predicted with growth times much shorter than plasma lifetimes. Similarly, growth rates are calculated for the planned APEX experiment as well as gamma-ray burst scenarios, suggesting that the instability may contribute to the early evolution of these systems.
Fluid and gyrokinetic turbulence in open field-line, helical plasmas
Two-fluid Braginskii codes have simulated open-field line turbulence for over a decade, and only recently, it has become possible to study these systems with continuum gyrokinetic codes. This work presents a first-of-its-kind comparison between fluid and (long-wavelength) gyrokinetic models in open field-lines, using the GDB and Gkeyll codes to simulate interchange turbulence in the Helimak device at the University of Texas [T. N. Bernard et al., Phys. Plasmas 26, 042301 (2019)]. Partial agreement is attained in a number of diagnostic channels when the GDB sources and sheath boundary conditions (BCs) are selected carefully, especially the heat-flux BCs that can drastically alter the temperature. The radial profile of the fluctuation levels is qualitatively similar and quantitatively comparable on the low-field side, although statistics such as moments of the probability density function and the high-frequency spectrum show greater differences. This comparison indicates areas for future improvement in both simulations, such as sheath BCs, and improvements in GDB like particle conservation and spatially varying thermal conductivity, in order to achieve better fluid-gyrokinetic agreement and increase fidelity when simulating experiments.
Corrigendum: Thermal dynamics in the flux-coordinate independent turbulence code GRILLIX
Here, we have discovered that one of the original articles findings, that the actual choice of the dimensionless parameter χ ∥i0 for the ion heat conductivity seemed to play no role for the simulation results, was not physical. Instead, due to a mistake in the handling of input parameters, χ ∥i0 was actually not used in the code – the dimensionless electron heat conductivity χ ∥e0 was used instead.
Simulations of tokamak boundary plasma turbulence transport in setting the divertor heat flux width
The BOUT++ code has been used to simulate edge plasma electromagnetic (EM) turbulence and transport, and to study the role of EM turbulence in setting the scrape-off layer (SOL) heat flux width λ q . More than a dozen tokamak discharges from C-Mod, DIII-D, EAST, ITER and CFETR have been simulated with encouraging success. The parallel electron heat fluxes onto the target from the BOUT++ simulations of C-Mod, DIII-D, and EAST follow the experimental heat flux width scaling of the inverse dependence on the poloidal magnetic field. Further turbulence statistics analysis shows that the blobs are generated near the pedestal pressure peak gradient region inside the separatrix and contribute to the transport of the particle and heat in the SOL region. Transport simulations indicate two distinct regimes: drift dominant regime and turbulence dominant regime. Goldston's heuristic drift-based (HD) model yields a consistent divertor heat flux width in the drift dominant regime. For C-Mod enhanced D α H-mode discharges, drifts and turbulence are competing in setting the divertor heat flux width, possibly due to its compact machine size and good pedestal confinement. The simulations for ITER and CFETR indicate that divertor heat flux width λ q of the future machines may no longer follows the 1/B pol,OMP HD-based empirical (Eich) scalings and the HD model gives a pessimistic limit of divertor heat flux width. The simulation results show a transition from a drift dominant regime to a turbulence dominant regime from current machines to future machines such as ITER and CFETR for two reasons. (1) The magnetic drift-based radial transport decreases due to large CFETR and ITER machine sizes. (2) The SOL turbulence thermal diffusivity increases due to larger turbulent fluxes ejected from the pedestal into the SOL when operating in a different pedestal structure, from an ELM-free H-mode pedestal regime to a small and grassy ELM regime.
Multigrid treatment of implicit continuum diffusion
Implicit treatment of diffusive terms of various differential orders common in continuum mechanics modeling, such as computational fluid dynamics, is investigated with spectral and multigrid algorithms in non-periodic 2D domains. In doubly periodic time dependent problems these terms are handled efficiently by spectral methods, but in non-periodic systems solved with distributed memory parallel computing and 2D domain decomposition, this efficiency is lost for a large number of processors. Here, we built and present here a multigrid algorithm for these types of problems that outperforms a spectral solution employing the highly optimized FFTW library. This solver is suitable for high performance computing and may be able to efficiently treat implicit diffusion of arbitrary order by introducing auxiliary equations of lower order. We test these solvers for fourth and sixth order diffusion with harmonic test functions as well as turbulent 2D magnetohydrodynamic simulations. It is also shown that an anisotropic operator without mixed-derivative terms improves model accuracy and speed, and we examine the impact that the various diffusion operators have on the energy, the enstrophy, and the qualitative aspect of a simulation.
Global 3D Braginskii simulations of the tokamak edge region of IWL discharges
A study of plasma turbulence and profile evolution in conditions of low (L-mode) and high (H-mode) confinement at the edge of an axisymmetric, nested circular flux-surface approximation to an inner wall limited (IWL) Alcator C-Mod discharge is presented, using numerical simulations with the global drift-ballooning (GDB) code. Further, GDB solves drift-reduced Braginskii two-fluid equations for electromagnetic low-frequency turbulence in a 3D annulus centered on the last closed flux-surface (LCFS). Three simulations that investigate the conditions of a reference L-mode, a high density, and a high temperature (or H-mode-like) shot were performed using realistic parameters. L-mode transport appears to be largely driven by drift resistive ballooning structures. Its pressure profile exhibits a near-SOL breakpoint that Mirror Langmuir Probes (MLP) detect in C-Mod. The high density simulation sees an increase in the size of convective cells and enhanced turbulent transport, while H-mode conditions develop improved confinement, balanced $E$ x $B$ and ion diamagnetic drifts in the closed-flux region, and spontaneous generation of temperature pedestal with a density pedestal remaining absent. A statistical characterization of the turbulence both in the SOL and the closed-flux region is presented.
CFD Analysis of the Aerodynamics of a Business-Jet Airfoil with Leading-Edge Ice Accretion
For rime ice - where the ice buildup has only rough and jagged surfaces but no protruding horns - this study shows two dimensional CFD analysis based on the one-equation Spalart-Almaras (S-A) turbulence model to predict accurately the lift, drag, and pressure coefficients up to near the stall angle. For glaze ice - where the ice buildup has two or more protruding horns near the airfoil's leading edge - CFD predictions were much less satisfactory because of the large separated region produced by the horns even at zero angle of attack. This CFD study, based on the WIND and the Fluent codes, assesses the following turbulence models by comparing predictions with available experimental data: S-A, standard k-epsilon, shear-stress transport, v(exp 2)-f, and differential Reynolds stress.
Computing Aerodynamic Performance of a 2D Iced Airfoil: Blocking Topology and Grid Generation
The ice accrued on airfoils can have enormously complicated shapes with multiple protruded horns and feathers. In this paper, several blocking topologies are proposed and evaluated on their ability to produce high-quality structured multi-block grid systems. A transition layer grid is introduced to ensure that jaggedness on the ice-surface geometry do not to propagate into the domain. This is important for grid-generation methods based on hyperbolic PDEs (Partial Differential Equations) and algebraic transfinite interpolation. A 'thick' wrap-around grid is introduced to ensure that grid lines clustered next to solid walls do not propagate as streaks of tightly packed grid lines into the interior of the domain along block boundaries. For ice shapes that are not too complicated, a method is presented for generating high-quality single-block grids. To demonstrate the usefulness of the methods developed, grids and CFD solutions were generated for two iced airfoils: the NLF0414 airfoil with and without the 623-ice shape and the B575/767 airfoil with and without the 145m-ice shape. To validate the computations, the computed lift coefficients as a function of angle of attack were compared with available experimental data. The ice shapes and the blocking topologies were prepared by NASA Glenn's SmaggIce software. The grid systems were generated by using a four-boundary method based on Hermite interpolation with controls on clustering, orthogonality next to walls, and C continuity across block boundaries. The flow was modeled by the ensemble-averaged compressible Navier-Stokes equations, closed by the shear-stress transport turbulence model in which the integration is to the wall. All solutions were generated by using the NPARC WIND code.