Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Transitional flow”

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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

A combined numerical and experimental study to elucidate primary breakup dynamics in liquid metal droplet-on-demand printing

Droplet-on-demand liquid metal jetting is emerging as a powerful technology for the additive manufacturing of metallic parts. The success of this method hinges on overcoming several technological challenges. Additionally, the principal one among these challenges is the controlled repeatable ejection of single uniform droplets. Due to the high density and surface tension of liquid metals, the droplet ejection process occurs near the minimal extremes of the printability phase diagram, defined by acceptable ranges for the Weber (We) and Ohnesorge (Oh) numbers. In this work, we experimentally demonstrate the satellite-free ejection of pneumatically actuated molten tin droplets in this extreme corner of printability and use a combination of high-speed video analysis and volume-of-fluid modeling to elucidate the droplet dynamics. While the simulations at low Oh and We can correctly describe several aspects of the breakup process, such as an increasing tail and pinch-point near the nozzle, no single parameter set can completely capture the droplet shape at breakup. Instead, the experimental droplet dynamics appear to include features from both high and low Oh breakup. This disagreement is ascribed to the incomplete description of the droplet ejection process including wetting and exit effects near the nozzle opening and surface effects such as transient cooling and oxide formation.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Insight into pseudo-film-boiling in supercritical water under simplified interfacial conditions

The focus of this paper is the physics of supercritical water film boiling through the simulation of natural convection of supercritical water in a rectangular cavity. It has been widely believed that above a fluid's critical point, only one phase exists. However, recent research has indicated the existence of distinct gas-like and liquid-like regions separated by the Widom line, which is the locus of the maxima of the specific heat. Along the Widom line, density decreases sixfold, viscosity drops by a factor of two, while specific heat spikes by an order of magnitude. These variations, specifically in density and viscosity, produce a thick pseudo-interface and flow dynamics behavior akin to film-boiling, denoted here as pseudo-two-phase flow dynamics. When the cavity is heated from the bottom, a pseudo-film of gas-like fluid is formed at the heated wall, and the flow complexity evolves, characterized by the formation of multiple vortices. Results demonstrate the ability of modern conjugate heat transfer computational fluid dynamics to capture the pseudo two-phase regime, with the formation of a pseudo-film at the heated wall. Further, the boiling curve is studied, showing an increase in heat flux with an increase in the temperature differential between the wall and the pseudo-critical temperature. Furthermore, results also show a previously uncharacterized instability, where the convection cells in the cavity oscillate and can flip directions under certain conditions.

Computational fluid dynamics↗
DOE OSTI2024

Time-resolved cryogenic electron tomography for the study of transient cellular processes

Cryogenic electron tomography (cryo-ET) is the highest resolution imaging technique applicable to the life sciences, enabling sub-nanometer visualization of specimens preserved in their near native states. The rapid plunge freezing process used to prepare samples lends itself to time-resolved studies, which researchers have pursued for in vitro samples for decades. Here, we focus on developing a freezing apparatus for time-resolved studies in situ. The device mixes cellular samples with solution-phase stimulants before spraying them directly onto an electron microscopy grid that is transiting into cryogenic liquid ethane. By varying the flow rates of cell and stimulant solutions within the device, we can control the reaction time from tens of milliseconds to over a second prior to freezing. In a proof-of-principle demonstration, the freezing method is applied to a model bacterium, Caulobacter crescentus, mixed with an acidic buffer. Through cryo-ET we resolved structural changes throughout the cell, including surface-layer protein dissolution, outer membrane deformation, and cytosolic rearrangement, all within 1.5 seconds of reaction time. This new approach, Time-Resolved cryo-ET (TR-cryo-ET), enhances the capabilities of cryo-ET by incorporating a sub-second temporal axis and enables the visualization of induced structural changes at the molecular, organelle, or cellular level.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

How the birth and death of shear layers determine confinement evolution: from the L → H transition to the density limit

Electric field profile structure—especially its shear—is a natural order parameter for the edge plasma, and characterizes confinement regimes ranging from the H-mode (Wagner et al. 1982 Phys. Rev. Lett.49, 1408–1412 (doi:10.1103/PhysRevLett.49.1408)) to the density limit (DL) (Greenwald et al. 1988 Nucl. Fusion28, 2199–2207 (doi:10.1088/0029-5515/28/12/009)). The theoretical developments and lessons learned during 40 years of H-mode studies (Connor & Wilson 1999 Plasma Phys. Control. Fusion42, R1–R74 (doi:10.1088/0741-3335/42/1/201); Wagner 2007 Plasma Phys. Control. Fusion49, B1–B33 (doi:10.1088/0741-3335/49/12b/s01)) are applied to the shear layer collapse paradigm (Hong et al. 2017 Nucl. Fusion58, 016041 (doi:10.1088/1741-4326/aa9626)) for the onset of DL phenomena. Results from recent experiments on edge shear layers and DL phenomenology are summarized and discussed in the light of L → H transition physics. The theory of shear layer collapse is then developed. In this work, we demonstrate that shear layer physics captures both the well known current (Greenwald) scaling of the DL (Greenwald 2002 Plasma Phys. Control. Fusion44, R27–R53 (doi:10.1088/0741-3335/44/8/201); Greenwald et al. 2014 Phys. Plasmas21, 110501 (doi:10.1063/1.4901920)), as well as the emerging power scaling (Zanca, Sattin, JET Contributors 2019 Nucl. Fusion59, 126011 (doi:10.1088/1741-4326/ab3b31)). The derivation of the power scaling theory exploits an existing model, originally developed for the L → H transition (Diamond, Liang, Carreras, Terry 1994 Phys. Rev. Lett.72, 2565–2568 (doi:10.1103/PhysRevLett.72.2565); Kim & Diamond 2003 Phys. Rev. Lett.90, 185006 (doi:10.1103/PhysRevLett.90.185006)). We describe the enhanced particle transport events that occur following shear layer collapse. Open problems and future directions are discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Multifidelity validation of digital surrogates using variable-density turbulent mixing models

In this study, ensembles of experimental data are presented and utilized to compare and validate two models used in the simulation of variable-density (Atwood = 0.22), compressible turbulent mixing. Though models of this kind (Reynolds averaged NavierStokes and large-eddy simulations) have been validated extensively with more canonical flows in previous studies, here the present approach offers novelty in the complexity of the geometry, the ensemble-based validation, and the uniformity of the computational framework on which the models are tested. Moreover, all experimental and computational tasks were completed by the authors which has led to a tightly coupled experimental configuration with its “digital twin.” The experimental divergent-shock-tube facility and its data acquisition methods are described and replicated in simulation space. A 2D Euler model which neglects the turbulent mixing at the interface is optimized to experimental data using a Gaussian process. This model then serves as the basis for both the 2D RANS and 3D LES studies that make comparisons to the mixing-layer data from the experiment. A relatively simple RANS model is shown to produce good agreement with experimental data only at late flow development times. The LES ensembles generally show good agreement with experimental data but display sensitivity to the characterization of initial conditions. Resolution-dependent behavior is also observed for certain higher-order statistics of interest. Overall, the LES model successfully captures the effects of divergent geometry, compressibility, and combined nonlinear instabilities inherent to the problem. The successful prediction of mixing width and its growth rate highlight the existence of three distinct regimes in the development of the instability, each with similarities to previously studied instabilities.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Coupling Between Alfvén Wave and Kelvin–Helmholtz Waves in the Low Latitude Boundary Layer

The Kelvin–Helmholtz (KH) instability of magnetohydrodynamic surface waves at the low latitude boundary layer is examined using both an eigenfrequency analysis and a time-dependent wave simulation. The analysis includes the effects of sheared flow and Alfvén velocity gradient. When the magnetosheath flows are perpendicular to the ambient magnetic field direction, unstable KH waves that propagate obliquely to the sheared flow direction occur at the sheared flow surface when the Alfvén Mach number is higher than an instability threshold. Including a shear transition layer between the magnetosphere and magnetosheath leads to secondary KH waves (driven by the sheared flow) that are coupled to the resonant surface Alfvén wave. There are remarkable differences between the primary and the secondary KH waves, including wave frequency, the growth rate, and the ratio between the transverse and compressional components. The secondary KH wave energy is concentrated near the shear Alfvén wave frequency at the magnetosheath with a lower frequency than the primary KH waves. Although the growth rate of the secondary KH waves is lower than the primary KH waves, the threshold condition is lower, so it is expected that these types of waves will dominate at a lower Mach number. Because the transverse component of the secondary KH waves is stronger than that of the primary KH waves, more efficient wave energy transfer from the boundary layer to the inner magnetosphere is also predicted.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Dust-void evolution driven by turbulent dust flux can induce runaway migration of Earth-mass planets

Torques from asymmetric dust structures (so-called dust-void and filamentary structures) formed around low-mass planets embedded in a nonturbulent dust-gas disk can exceed the torques produced by the gas disk component and then go on to dominate the planet’s orbital dynamics. Here, we investigate how these structures (hence the dust torque) change when the effect of turbulent dust diffusion and dust feedback are included, along with the direct implications on the migration of Earth-like planets. Using the FARGO3D code, we performed 2D and 3D multifluid hydrodynamic simulations, focusing on a non-migrating planet with a mass of M p = 1.5 M ⊕ in 2D and on migrating planets with M p ∈ [1.5, 12] M ⊕ in 3D. We varied the δ-dimensionless diffusivity parameter in the range [0, 3 × 10 −3 ] and considered three different Stokes numbers, St = {0.04, 0.26, 0.55}, which are representative of the gas-dominated, the transitional, and the gravity-dominated regimes, respectively. In our 2D models, we find that turbulent diffusion of dust prevents the formation of the dust-void and filamentary structures when δ > 3 × 10 −4 . Otherwise, dust structures survive turbulent diffusion flow. However, dust and total torques become positive only in transitional and gravity-dominated regimes. In our 3D models, we find that the dust-void is drastically modified and the high-density ring-shaped barrier delineating the dust-void disappears if δ ≳ 10 −4 , due to the effect of dust turbulent diffusion along with the back-reaction of the dust. For all values of δ, the filament in front of the planet is replaced by a low-density trench. Remarkably, as we have allowed the planets to migrate, the evolving dust-void can drive either runaway migration or outward (inward) oscillatory-torque migration. Our study thus suggests that low-mass Earth-like planets can undergo runaway migration in dusty disks.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

In-situ Observations of Directed Energy Deposition Additive Manufacturing Using High-Speed X-ray Imaging

In laser-based directed energy deposition (DED) additive manufacturing, interactions among the laser beam, particle flow, and melt pool influence the properties of the solidified final part. Two separate DED systems, one with high powder flow rates to represent industrial-scale DED processing and the other with low powder flow rates for individual particle tracking, were synchronized with the high-speed imaging setup at the Advanced Photon Source in Argonne National Laboratory. In-situ x-ray imaging of the DED process using both systems highlighted the influence of powder flow rates. Increased powder flow rates resulted in less laser absorption into the melt pool, leading to a transition from a keyhole mode to a melt pool without a keyhole but with surface fluctuations due to powder flow. Increased velocities of particles during powder flow resulted in a decrease in particle melting times and a greater propensity for porosity formation. Altogether, better understanding of the interactions that occur during various scales of the DED process will enable flexibility, control, and new materials development in DED-based additive manufacturing.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Zonally dominated dynamics and Dimits threshold in curvature-driven ITG turbulence

The saturated state of turbulence driven by the ion-temperature-gradient instability is investigated using a two-dimensional long-wavelength fluid model that describes the perturbed electrostatic potential and perturbed ion temperature in a magnetic field with constant curvature (a $Z$ -pinch) and an equilibrium temperature gradient. Numerical simulations reveal a well-defined transition between a finite-amplitude saturated state dominated by strong zonal-flow and zonal temperature perturbations, and a blow-up state that fails to saturate on a box-independent scale. We argue that this transition is equivalent to the Dimits transition from a low-transport to a high-transport state seen in gyrokinetic numerical simulations (Dimits et al. , Phys. Plasmas , vol. 7, 2000, 969). A quasi-static staircase-like structure of the temperature gradient intertwined with zonal flows, which have patch-wise constant shear, emerges near the Dimits threshold. The turbulent heat flux in the low-collisionality near-marginal state is dominated by turbulent bursts, triggered by coherent long-lived structures closely resembling those found in gyrokinetic simulations with imposed equilibrium flow shear (van Wyk et al. , J. Plasma Phys. , vol. 82, 2016, 905820609). The breakup of the low-transport Dimits regime is linked to a competition between the two different sources of poloidal momentum in the system – the Reynolds stress and the advection of the diamagnetic flow by the $\boldsymbol {E}\times \boldsymbol {B}$ flow. By analysing the linear ion-temperature-gradient modes, we obtain a semi-analytic model for the Dimits threshold at large collisionality.

Physics↗
DOE OSTI2021

Evaluating the Effects of Transient Purge Flow on Stator-Rotor Seal Performance

Abstract As modern engine designs target higher efficiencies through increased turbine inlet temperatures, critical turbine components are at increased risk of damage from conditions exceeding material melting temperatures. In particular, improperly designed underplatform hardware components are susceptible to damage when hot main gas path flow is ingested into the stator-rotor cavity. While all turbines inherently experience transients during operation, a majority of publicly available turbine studies have been executed using steady operating conditions or inherently transient “blowdown” rigs. For this reason, routine transient events are not well understood. To address this need, the present study utilized a combined experimental and computational approach. The test article is a continuous-duration, one-stage test turbine operating with true-scale engine hardware and seal geometries at engine-representative flow conditions. The nature of the continuous-duration facility uniquely supports the direct assessment of transient events through its ability to transition between steady-state operating conditions. The effects of a transient purge flow were investigated in this study to identify general trends for transient events in a full-scale engine. Results from multiple measurement techniques in the wheelspace region show an interdependence of transient purge flow with a thermal lag of the underplatform hardware. Through experiments conducted at different coolant-to-main gas path temperature ratios, the use of pressure measurements as an indicator of fully purged behavior was introduced, and a thermally driven influence on rim seal performance was quantified. The computational results show good agreement with experimental pressure measurements and provide insight into the physical mechanisms that drive the relationship between pressure and sealing effectiveness measurements observed for the tested geometry.

Engineering↗
DOE OSTI2021

Coupling between Alfven wave and Kelvin-Helmholtz waves in the low latitude boundary layer

The Kelvin-Helmholtz (KH) instability of magnetohydrodynamic surface waves at the low latitude boundary layer is examined using both an eigenfrequency analysis and a time-dependent wave simulation. The analysis includes the effects of sheared flow and Alfven velocity gradient. When the magnetosheath flows are perpendicular to the ambient magnetic field direction, unstable KH waves that propagate obliquely to the sheared flow direction occur at the sheared flow surface when the Alfv\'en Mach number is higher than an instability threshold. Including a shear transition layer between the magnetosphere and magnetosheath leads to secondary KH waves (driven by the sheared flow) that are coupled to the resonant surface Alfven wave. There are remarkable differences between the primary and the secondary KH waves including wave frequency, the growth rate, and the ratio between transverse and the compressional component. The secondary KH wave energy is concentrated near the shear Alfven wave frequency at the magnetosheath with a lower frequency than the primary KH waves. Although the growth rate of the secondary KH waves is lower than the primary KH waves, the threshold condition is lower, so it is expected that these types of waves will dominate at lower Mach number. Because the transverse component of the secondary KH waves is stronger than the primary KH waves, more efficient wave energy transfer from the boundary layer to the inner magnetosphere is also predicted.

Alfven wave↗
DOE OSTI2022

High-energy and low-cost membrane-free chlorine flow battery

Abstract Grid-scale energy storage is essential for reliable electricity transmission and renewable energy integration. Redox flow batteries (RFB) provide affordable and scalable solutions for stationary energy storage. However, most of the current RFB chemistries are based on expensive transition metal ions or synthetic organics. Here, we report a reversible chlorine redox flow battery starting from the electrolysis of aqueous NaCl electrolyte and the as-produced Cl 2 is extracted and stored in the carbon tetrachloride (CCl 4 ) or mineral spirit flow. The immiscibility between the CCl 4 or mineral spirit and NaCl electrolyte enables a membrane-free design with an energy efficiency of >91% at 10 mA/cm 2 and an energy density of 125.7 Wh/L. The chlorine flow battery can meet the stringent price and reliability target for stationary energy storage with the inherently low-cost active materials (~$5/kWh) and the highly reversible Cl 2 /Cl − redox reaction.

25 ENERGY STORAGE↗
DOE OSTI2023

Study of Seamless Microgrid Transition Operation Using Grid-Forming Inverters

This paper investigates operational techniques to achieve seamless (smooth) microgrid (MG) transitions by dispatching a grid-forming (GFM) inverter. In traditional approaches, the GFM inverter must switch between grid-following (GFL) and GFM control modes during MG transition operation. Today's inverter technology allows GFM inverters to always operate in GFM control mode, so it is worth exploring how to use them to achieve smooth MG transition operation. This paper proposes three operational techniques: a traditional scheme of switching between GFL and GFM control; a new scheme of consistent GFM control and shifting the droop intercept up before islanding operation; and a new scheme of consistent GFM control and shifting the droop intercept up before synchronization operation. A full hardware setup is established to compare the three techniques and showcase their implementations in real-world applications. The results show that the third technique outperforms the others and exhibits the best transition performance because the GFM inverter maintains the same operating points during the transition operation. Therefore, we conclude that ensuring smooth MG transition operation requires that the GFM inverter(s) maintain the same operating points (v, f, P, Q, and phase angle) during the transition operation in addition to minimizes the point of common coupling power flow.

grid-forming↗
DOE OSTI2022

Full Length Assembly Testing in PELICAN (Final Report)

In support of the development of the U.S. Department of Energy (DOE) Versatile Test Reactor (VTR), a thermal hydraulics test facility was constructed to generate experimental measurement of the pressure drop across a single full-scale assembly containing prototypic axial reflectors, fuel, and plena components. Constructed and operated at Argonne National Laboratory, the Pressure drop Experimental Loop for Investigations of Core Assemblies in Nuclear reactors (PELICAN) facility was designed to achieve hydraulic conditions identical to those anticipated for a full-scale fuel assembly located in the VTR core in the region with the highest flow rate. Using water as surrogate for liquid sodium, the flow loop was operated at elevated temperatures and pressures to match the thermophysical properties of liquid sodium and ensure matching Reynolds and Euler numbers. The measurement objectives for data generated from this test facility was driven primarily by the validation needs for code calculations and simulations of the reference VTR core. These objectives focused on the need to validate pressure drop results across the various segments of the fuel assembly as they relate directly to the pumping power and safety of the reactor. Presented in this report are experimental results and analytical comparisons based on testing of a full-length assembly in PELICAN. Housed within a hexagonal test section extending 3.4 m in length, the tested assembly features a prototypic lower reflector, grid plates, wire-wrapped rod bundle, upper reflector, and exit region. The rod bundle extends over 1.5 m in length and contains 217 individual wire-wrapped rods with dimensions that best reflect the reference VTR design. The as-tested bundle assembly was fabricated using 316 stainless steel 0.25-inch (6.35-mm) diameter rods wrapped with 0.04-inch (1.016-mm) diameter wire at a helical pitch of 10.51 inch (26.6 cm). Details of the method for in-house wire-wrapping, assembly, and installation are provided later in this report. Experimental measurements of pressure drop at 19 positions along the test assembly were recorded for a range of flow conditions, with special attention paid to key locations within the assembly, including component inlet and outlet, transition, and wire-wrapped rod bundle regions. Testing conditions were based on 110°C water with flow rates ranging from 50 to 450 GPM (3 to 27 kg/s) at the inlet of the test assembly generating Reynolds numbers and velocities up to ~8.0×10 4 and ~7.8 m/s, respectively, within the rod bundle region. Non-dimensional values for the friction factor were then calculated based on these experimental measurements and compared against those predicted by various analytical correlations available from open literature. Predictions by the upgraded Cheng and Todreas, Rehme, and Novendstern correlations fell within 4% to those values measured experimentally.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Surface-to-Interior Transport Timescales and Ventilation Patterns in a Time-Dependent Circulation Driven by Sustained Climate Warming

Abstract The effect of climate warming in response to rising atmospheric CO 2 on the ventilation of the ocean remains uncertain. Here we make theoretical advances in elucidating the relationship between ideal age and transit time distribution (TTD) in a time-dependent flow. Subsequently, we develop an offline tracer-transport model to characterize the ventilation patterns and time scales in the time-evolving circulation for the 1850–2300 period as simulated with the Community Earth System Model version 1 (CESMv1) under a business-as-usual warming scenario. We found that by 2300 2.1% less water originates from the high-latitude deep water formation regions (both hemispheres) compared to 1850. In compensation, there is an increase in the water originating from the subantarctic. We also found that slowing meridional overturning circulation causes a gradual increase in mean age during the 1850–2300 period, with a globally averaged mean-age increase of ∼110 years in 2300. Where and when the water will be re-exposed to the atmosphere depends on the post-2300 circulation. For example, if we assume that the circulation persists in its year-2300 state (scenario 1), the mean interior-to-surface transit time in year 1850 is ∼1140 years. In contrast, if we assume that the circulation abruptly recovers to its year-1850 state (scenario 2), the mean interior-to-surface transit time in 1850 is only ∼740 years. By 2300, these differences become even larger; in scenario 1, the mean interior-to-surface transit time increases by ∼200 years, whereas scenario 2 decreases by ∼80 years. The dependence of interior-to-surface transit time on the future ocean circulation produces an additional unavoidable uncertainty in the long-term durability of marine carbon dioxide removal strategies. Significance Statement The ocean’s circulation, when altered by climate warming, can affect its capacity to absorb heat and CO 2 , which are crucial for the global climate. In our study, we investigated how global warming, caused by rising CO 2 levels, might impact the ocean circulation—the way water moves from deep ocean to the surface and vice versa. We discovered that by 2300, if we continue on our current warming trajectory, the origins of water within the ocean will shift, with less coming from deep, cold zones near the poles and more from subantarctic regions. As a result, deep water will take longer time before it resurfaces than shallow water. How quickly this water travels from deep regions to the surface could change, depending on the state of future ocean circulation. If the circulation remains as predicted in 2300, this journey will take longer. Conversely, if it reverts to the pattern in 1850, the process will be quicker. This variability introduces added uncertainty to strategies aimed at mitigating climate change by storing CO 2 in the ocean. Our work highlights the intricate ways in which climate change can influence our oceans, potentially affecting our plans to mitigate global warming.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Fragile-to-strong transition in liquid As 2 S 3 under pressure: The effect of melt metallization

The well-known classification of glass-forming melts into fragile and strong liquids has several notable exceptions, including water, silica, and certain phase-change materials (PCMs). These exceptional fluid systems exhibit a fragile-to-strong transition (FST) upon cooling: a transformation from a high-temperature liquid with fast atomic dynamics, low viscosity, and low flow activation energy, to a viscous supercooled melt with high energy barriers near the glass transition temperature T g . This behavior is critically important for non-volatile memories, photonic tensor cores, reconfigurable metamaterials, and other devices, that use PCMs, enabling nanosecond-scale crystallization in the fragile regime and long data retention in the strong regime near or below T g . A significant structural transformation is expected between these two viscosity regimes, along with a semiconductor-metal (SC-M) transition upon heating, driven by high internal pressure and associated density increase. By applying high external pressure to the canonical low-conducting chalcogenide melt As 2 S 3 , we observed both the FST and the SC-M transition, occurring simultaneously within the same domain of the P, T−phase space. These findings suggest that the FST is not limited to a few exceptional liquids but is a common phenomenon, at least in systems that exhibit melt metallization within specific regions of their P, T−phase diagrams.

first-principles molecular dynamics↗
DOE OSTI2021

Explicit structure-preserving geometric particle-in-cell algorithm in curvilinear orthogonal coordinate systems and its applications to whole-device 6D kinetic simulations of tokamak physics

Explicit structure-preserving geometric particle-in-cell (PIC) algorithm in curvilinear orthogonal coordinate systems is developed. The work reported represents a further development of the structure-preserving geometric PIC algorithm achieving the goal of practical applications in magnetic fusion research. The algorithm is constructed by discretizing the field theory for the system of charged particles and electromagnetic field using Whitney forms, discrete exterior calculus, and explicit non-canonical symplectic integration. In addition to the truncated infinitely dimensional symplectic structure, the algorithm preserves exactly many important physical symmetries and conservation laws, such as local energy conservation, gauge symmetry and the corresponding local charge conservation. As a result, the algorithm possesses the long-term accuracy and fidelity required for first-principles-based simulations of the multiscale tokamak physics. The algorithm has been implemented in the SymPIC code, which is designed for high-efficiency massively-parallel PIC simulations in modern clusters. The code has been applied to carry out whole-device 6D kinetic simulation studies of tokamak physics. A self-consistent kinetic steady state for fusion plasma in the tokamak geometry is numerically found with a predominately diagonal and anisotropic pressure tensor. The state also admits a steady-state sub-sonic ion flow in the range of 10 km s -1 , agreeing with experimental observations and analytical calculations Kinetic ballooning instability in the self-consistent kinetic steady state is simulated. It is shown that high-n ballooning modes have larger growth rates than low-n global modes, and in the nonlinear phase the modes saturate approximately in 5 ion transit times at the 2% level by the E × B flow generated by the instability. These results are consistent with early and recent electromagnetic gyrokinetic simulations.

43 PARTICLE ACCELERATORS↗