Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Thermal fluid physics”

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

Preliminary Plan to Inform Testing of a Heat Exchanger Test Article

This report presents a preliminary plan to guide the qualification testing of advanced heat exchanger (HX) components for nuclear-to-industrial heat transfer applications. The objective is to establish a defensible, physics-based methodology that integrates computational modeling, targeted experimentation, and in-service inspection considerations to demonstrate component performance and reliability under representative reactor conditions. The analysis identifies Sodium-cooled Fast Reactor (SFR) and High-Temperature Gas-cooled Reactor (HTGR) systems as reference configurations in terms of temperature, pressure, and chemical environment. Within these operating envelopes, dominant degradation mechanisms— including creep–fatigue interaction, flow-induced vibration, corrosion, and diffusion-bond deterioration—were evaluated to define test requirements. A comprehensive computationalexperimental framework is proposed to support life prediction and qualification activities. The framework couples high-fidelity structural-mechanics, thermal-hydraulic, and fluid-structure interaction models with accelerated degradation testing to produce a traceable linkage between microstructural evolution, mechanical performance, and remaining useful life (RUL). The approach adheres to established Verification, Validation, and Uncertainty Quantification (VVUQ) standards (ASME V&V 10/20; NUREG-2152) and incorporates a digital-twin architecture for continuous model refinement through data assimilation. The plan further outlines testing methodologies, including pre-test analyses, test-loop design parameters, and sensor placement strategies that maximize information yield while maintaining mechanistic fidelity. Complementary sections describe in-service inspection (ISI), on-line monitoring (OLM), and structural-health-monitoring (SHM) techniques applicable to compact HX geometries typical of advanced reactors. Collectively, these activities establish the technical foundation for demonstrating 40-60-year equivalent service life of advanced heat exchangers in support of the U.S. Department of Energy’s Advanced Reactor and Integrated Energy Systems programs. The forthcoming phase will execute the defined pre-test analyses, initiate hardware fabrication, and implement the integrated testing campaign to validate the proposed qualification methodology.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

A simple and effective five-equation two-phase numerical model for liquid-vapor phase transition in cavitating flows

Numerical difficulties, notably the non-monotonic behavior of the Wood speed of sound and the volume fraction positivity, associated with the reduced five-equation two-phase flow model of Kapila et al. (2001) [A.K. Kapila, R. Menikoff, J.B. Bdzil, S.F. Son, D.S. Stewart, 2001. Two-phase modeling of deflagration-to-detonation transition in granular materials: reduced equations, Physics of Fluids 13(10), 3002–3024] have been resolved in the past through the introduction of a frozen speed of sound and an algebraic approach for mechanical relaxation afforded by a pressure non-equilibrium six-equation model proposed by [R. Saurel, F. Petitpas, R.A. Berry, 2009. Simple and efficient relaxation methods for interfaces separating compressible fluids, cavitating flows and shocks in multiphase mixture, J. Comput. Phys. 228, 1678–1712]. By contrast, it is explored and demonstrated in this work that these difficulties can in fact still be resolved within the numerical scheme for solving the reduced five-equation model by numerically replacing the Wood speed of sound for the estimates of wave speeds in the approximate Riemann solver HLLC with the monotonic mixture speed of sound for a transport five-equation model. For shock interface (artificial mixture separating pure or nearly pure fluids) interaction problems, with the apparent advantage of monotonic behavior of the speed of sound in the interface, the effect of the numerical replacement is also confined to the interface. Differences other than the behavior of the speed of sound within the interface in the solutions due to the replacement diminish with increasing resolution when reasonable solution can be obtained with Wood speed of sound. For cavitating/expansion problems in physical fluid mixture, it is pointed out and explained why acoustics in the numerical solutions still propagate at the Wood speed of sound (therefore consistent with the reduced five-equation model) even though in some cases a much higher speed of sound like the numerical replacement above for solving the reduced five-equation model or the frozen speed of sound for solving a six-equation model is used for the estimates of wave speeds in the HLLC scheme. A variant of the five-equation two-phase flow model by Saurel et al. (2008) [R. Saurel, F. Petitpas, R. Abgrall, 2008. Modelling phase transition in metastable liquids: application to cavitating and flashing flows, J. Fluid Mech. 607, 313–350] is then constructed for liquid-vapor phase transition in cavitating flows. The relaxation toward thermo-chemical equilibrium during phase transition is achieved by solving a simple system of algebraic equations for the equilibrium state variables for better efficiency, following Pelanti and Shyue (2014) [M. Pelanti, K.-M. Shyue, 2014. A mixture-energy-consistent six-equation two-phase numerical model for fluid with interfaces, cavitation and evaporation waves. J. Comput. Phys. 259, 331–357]. Therefore, the current model retains both the simplicity afforded by the five-equation model and the efficiency of the algebraic relaxation solver. An alternative algebraic approach for handling the non-conservative term (the so-called K∇ · u term) in the reduced five-equation model for mechanical equilibrium of a liquid-vapor mixture is also explored by enforcing the thermal equilibrium at the same time. Finally, numerical results of sample tests in both one and two dimensions in the literature as well as that in three dimensions demonstrate the effectiveness and ability of the proposed model to simulate cavitating flows. An interesting mechanism of shock generation by acoustics in water due to phase transition is then found by the numerical simulations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Machine learning of consistent thermodynamic models using automatic differentiation

In this study, we propose a data-driven method to describe consistent equations of state (EOS) for arbitrary systems. Complex EOS are traditionally obtained by fitting suitable analytical expressions to thermophysical data. A key aspect of EOS is that the relationships between state variables are given by derivatives of the system free energy. In this work, we model the free energy with an artificial neural network and utilize automatic differentiation to directly learn the derivatives of the free energy. We demonstrate this approach on two different systems, the analytic van der Waals EOS and published data for the Lennard-Jones fluid, and we show that it is advantageous over direct learning of thermodynamic properties (i.e., not as derivatives of the free energy but as independent properties), in terms of both accuracy and the exact preservation of the Maxwell relations. Furthermore, the method implicitly provides the free energy of a system without explicit integration.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Impact of moment-based, energy integrated neutrino transport on microphysics and ejecta in binary neutron star mergers

We present an extensive study of the effects of neutrino transport in three-dimensional general relativistic radiation hydrodynamics (GRHD) simulations of binary neutron star (BNS) mergers using our moment-based, energy-integrated neutrino radiation transport (M1) scheme. Here, we consider a total of eight BNS configurations, while varying equation of state models, mass ratios, and grid resolutions, for a total of 16 simulations. We find that M1 neutrino transport is crucial in modeling the local absorption of neutrinos and the deposition of lepton number throughout the medium. We provide an in-depth look at the effects of neutrinos on the fluid dynamics and luminosity during the late inspiral and postmerger phases, the properties of ejecta and outflow, and the postmerger nucleosynthesis. The simulations presented in this work comprise an extensive study of the combined effect of the equation of state and M1 neutrino transport in GRHD simulations of BNS mergers, and establish that the solution provided by our M1 scheme is robust across system properties.

150 ≤ A ≤ 18959 ≤ A ≤ 8990 ≤ A ≤ 149↗

Modeling Geologic Waste Repository Systems Below Residual Saturation

The heat generated by high-level radioactive waste can pose numerical and physical challenges to subsurface flow and transport simulators if the liquid water content in a region near the waste package approaches residual saturation due to evaporation. Here, residual saturation is the fraction of the pore space occupied by liquid water when the hydraulic connectivity through a porous medium is lost, preventing the flow of liquid water. While conventional capillary pressure models represent residual saturation using asymptotically large values of capillary pressure, here, residual saturation is effectively modeled as a tortuosity effect alone. Treating the residual fluid as primarily dead-end pores and adsorbed films, relative permeability is independent of capillary pressure below residual saturation. To test this approach, PFLOTRAN is then used to simulate thermal-hydrological conditions resulting from direct disposal of a dual-purpose canister in unsaturated alluvium using both conventional asymptotic and revised, smooth models. Importantly, while the two models have comparable results over 100 000 years, the number of flow steps required is reduced by approximately 94%.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Multiscale Interaction of Turbulence, Temperature, and Soot Formation: Measurements for Critical Assessments of Chemical Kinetics and Mechanisms

This program targeted an improved understanding of how turbulence affects soot formation chemistry (directly and indirectly) with a particular focus on detailing the multiscale coupling between flow turbulence, mixing, thermal transport, and soot formation kinetics in gas-phase reacting systems. In this project, simultaneous, quantitative multi-dimensional velocity, soot volume fraction, and gas-phase temperature measurements in turbulent non-premixed sooting flames were developed. These new diagnostics tools can be used to elucidate the relative effects of fluid kinematics, mixing, thermal processes, and turbulence/scalar time/length scales on soot formation chemistry and topology. Simultaneous multi-parameter measurements were acquired within a series of flames with different Reynolds numbers and spatial locations to investigate the inter-dependence of soot, temperature, and turbulence within the soot inception, growth, and oxidation region. Preliminary assessment of the results has been completed at the project end and detailed analysis is planned following the completed program.

02 PETROLEUM↗

A High-Fidelity Modeling Tool to Support the Design of Oxy-Combustors for Direct-Fired Supercritical CO2 Cycles

The challenge in the design of oxy-combustors for direct-fired supercritical CO 2 (sCO 2 ) cycles is in addressing disparate performance metrics and objectives. Key design parameters to consider include, among others, injector design for mixing and flame stability, split of recycled CO 2 diluent between injectors and cooling films, target flame temperatures to control noncondensable products, and strategies to inject the diluent CO 2 for film cooling and thermal control. In order to support novel oxy-combustor designs, a high-fidelity yet numerically efficient modeling framework based on the CRUNCH CFD® flow solver is presented, featuring key physics-based submodels relevant in this regime. For computational efficiency in modeling large kinetic sets, a flamelet/progress variable (FPV) based tabulated-chemistry approach is utilized featuring a three-stream extension to allow for the simulation of the CO 2 film cooling stream in addition to the fuel and oxidizer streams. Finite rate chemistry effects are modeled in terms of multiple progress variables for the primary flame as well as for slower-evolving chemical species such as NO x and SO x contaminants. Real fluid effects are modeled using advanced equations of states. The predictive capabilities of this computationally tractable design support tool are demonstrated on a conceptual injector design for an oxy-combustor operating near 30 MPa. Simulations results provide quantitative feedback on the effectiveness of the film cooling as well as the level of contaminants (CO, NO, and N) in the exhaust due to impurities entering from the injectors. These results indicate that this framework would be a useful tool for refining and optimizing the oxy-combustor designs as well as risk mitigation analyses.

Engineering↗

L-H transition trigger physics in ITER-similar plasmas with applied n = 3 magnetic perturbations

The L–H transition power threshold PLH is observed to increase with applied n = 3 resonant magnetic perturbations (RMP) in ITER-similar-shape plasmas with balanced neutral beam torque injection in DIII-D. The increase is most pronounced with added electron–cyclotron heating: PLH increases with decreasing edge plasma collisionality as PLH/PLH-08 ~ (ν*)−0.5, where PLH-08 is the 2008 ITPA multi-machine power threshold scaling. This result raises concerns for H-mode access at low edge collisionality in ITER, where RMP may have to be applied before the L–H transition to safely suppress the first edge-localized mode. Non-axisymmetric modifications with RMP include a simultaneous reduction of the radial electric field (Er) well depth and E × B shear. This can be attributed to increasing edge toroidal co-current rotation, and is consistent with substantially increased local long-wavelength turbulence (measured via beam emission spectroscopy). At high RMP perturbation strength the edge electric field Er reverses sign locally (becomes positive), with changes in dominant turbulence modes. Edge magnetic stochasticity provides an attractive explanation of the observed modifications, and the observed changes in toroidal rotation and Er are consistent with a simple fluid model describing radial electron current flow along stochastic fieldlines. The observed collisionality dependence of the L-mode edge electric field with applied RMP is also qualitatively consistent with this model. Reflectometry data indicate a significant reduction of the normalized L-mode radial density gradient a/Ln at high RMP field with simultaneous increase in radial particle flux and electron thermal flux from power balance analysis. We conjecture that the increase of PLH with RMP results from the combined effects of reduced E × B flow shear (increasing turbulent transport levels) and toroidal/poloidal flow modulation due to edge stochasticity. Initial experiments indicate that non-resonant n = 3 magnetic perturbations lead only to relatively small changes in Er, E × B shear and fluctuation characteristics, and have less impact on the L–H transition power threshold. This motivates further exploration of the RMP spectrum dependence of PLH for possible mitigation of the observed threshold increase.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Non-local transport in radiation-hydrodynamics codes for ICF by efficient coupling to an external Vlasov–Fokker–Planck code

Accurately incorporating non-local transport into radiation-hydrodynamics codes, and indeed any fluid system, has long been elusive. To date, a simplified and accurate theory that can be easily integrated has not been available. This limitation affects modeling in inertial confinement fusion (ICF) and magnetic confinement fusion systems, among others, where non-local transport is well-known to be present. Here, we present a coupling methodology between a full Vlasov–Fokker–Planck (VFP) electron kinetic code and radiation-hydrodynamics (rad-hydro) codes. The VFP code is used to adjust native electron transport in the rad-hydro code, thus enabling improved transport without the need to integrate a full electron VFP solver into the rad-hydro code. This approach necessitates only occasional invocation of the VFP code, reducing computational intensity compared to following the dynamic evolution entirely with the VFP code on fluid time scales. We illustrate that the methodology is more accurate than other simplified methods in thermal decay systems relevant to ICF and can replicate standard theoretical results with high accuracy.

Electronic transport↗

Predicted Heat Flux Performance of Actively Cooled Tungsten-Armored Graphitic Foam Monoblocks

Tungsten (W)–armored graphitic foam monoblocks were developed for applications requiring high-Z plasma-facing material in long-pulse fusion experiments and ultimately deuterium-tritium fusion reactors. The monoblocks are an integrated material system combining the advantages of a chemical vapor deposited (CVD) W coating with a high-conductivity graphitic foam. The W is a high-melting-point, high-Z material with low tritium retention. The graphitic foam coupled to a swirl tube serves as a high-thermal-conductivity heat sink that cannot melt, although it can sublime at much higher temperatures than copper melts. Together, they comprise a robust plasma-facing component (PFC) weighing roughly 5% of an all-W component or 17% of a traditional W-coated copper heat sink.A single-channel mock-up consisting of four graphitic foam monoblocks equipped with a water-cooled swirl tube was fabricated for eventual testing in the 60-kW, EB-60, rastered electron beam at the Applied Research Laboratory of The Pennsylvania State University. Two monoblocks have a thin 50-μm-thick coating of pure W chemically vapor deposited over NbC and pure Nb interlayers. Two others have a 2-mm-thick pure W coating CVD on graphitic monoblocks using the same interlayers. The mock-up will be cooled with available 10 m/s, 0.7 MPa water with a 22°C inlet temperature and subjected to varying uniform heat loads up to 20 MW/m 2 . It is equipped with type-K thermocouples at various depths, and calibrated infrared thermography and spot pyrometry will be used to characterize the heated surface. Real-time water calorimetry will be used to ascertain the absorbed steady-state power and infer the heat flux during testing. Since testing cannot be done under prototypic divertor flow conditions, it is necessary to predict the thermal response of this novel PFC system and investigate the power sharing between radiation and convection at divertor heat flux levels and its inherent ability to avoid critical heat flux. Results are reported for predictions obtained from computational fluid dynamics models up to 30 MW/m 2 of steady-state uniform heat flux. Leading-edge heat loads of 30 MW/m 2 on a 2-mm-wide side strip were also investigated to ascertain if coating delamination is likely.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Effects of non-thermal termolecular reactions on wedge-induced oblique detonation waves

The shock-induced combustion ramjet (Shcramjet) based on oblique detonation waves (ODWs) is among the promising choices for hypersonic propulsion systems. An understanding of the ignition, propagation, and stability of ODWs is critical to harnessing their propulsive potential. In such high speed reacting flows, there is a high probability of occurrence of non-thermal reactions due to the presence of non-trivial amounts of highly reactive radicals including H, O and OH apart from O 2 as demonstrated recently [M. P. Burke, S. J. Klippenstein, Nat. Chem. 9 (2017) 1078–1082, Y. Tao, A. W. Jasper, Y. Georgievskii, S. J. Klippenstein, R. Sivaramakrishnan, Proc. Combust. Inst. 38 (2021) 515–522]. The present work focuses on examining the initiation, propagation and structure of oblique detonation waves in stoichiometric H 2 -air mixtures through numerical simulations with and without non-thermal reactivity on a two-dimensional adaptive grid. Non-thermal reactions were included in the macroscopic kinetic model as chemically termolecular reactions facilitated by the H + OH radical-radical recombination and the H + O 2 radical-molecule association reactions. Since, the non-thermal reactions result in a corresponding decrease in the reaction fluxes of the incipient recombination/association reactions, an additional simulation was performed by applying corrections to the respective incipient recombination/association rate constants using the methodology demonstrated by Tao et al. [Y. Tao, A. W. Jasper, Y. Georgievskii, S. J. Klippenstein, R. Sivaramakrishnan, Proc. Combust. Inst. 38 (2021) 515–522]. Results show that, under ODWE relevant conditions, non-thermal reactivity fundamentally alters the induction length, intensity as well as the structure of the ODW. Specifically, it is found that non-thermal reactivity leads to a noticeable reduction in initiation length and a simultaneous increase in instantaneous peak heat release rate and the degree of unsteadiness of the ODW. Finally, statistical analysis of key thermo-chemical variables is performed to elucidate the important species as well as reactions responsible for the observed variations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A conservative phase-space moving-grid strategy for a 1D-2V Vlasov–Fokker–Planck Solver

In this work, we develop a conservative configuration- and velocity-space (i.e., phase-space) moving-grid strategy for the Vlasov–Fokker–Planck (VFP) equation in a planar geometry. The velocity-space grid is normalized and shifted in terms of the thermal speed and the bulk-fluid velocity, respectively. The configuration-space grid is moved according to a mesh-motion-partial-differential equation (MMPDE), which equidistributes a monitor function that is inversely proportional to the gradient-length scales of the macroscopic plasma quantities. The resulting inertial terms in the transformed VFP equations are discretized to ensure the discrete conservation of mass, momentum, and energy. To satisfy the discrete conservation theorems in the presence of phase-space mesh motion, we employ the method of discrete nonlinear constraints – explored in previous studies – but the underlying symmetries are determined in a much more efficient manner than before. The conservative grid-adaptivity strategy provides an efficient scheme that resolves important physical structures in the phase-space while controlling the computational complexity at all times. We demonstrate the favorable features of the proposed algorithm through a set of test cases of increasing complexity. The problems test independent components of the algorithms, as well as the integrated capability on settings relevant to inertial confinement fusion.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

BlueCRAB Domain Overlapping Coupling: Theory & Verification

Coupling low- and high-fidelity codes is a useful way to model complex engineering systems. Thanks to the high-fidelity code, complex phenomena can be resolved in areas of the system where this is required, yet the efficiency of the lower-fidelity code is still retained in modeling the rest of the system. This document details the theory and implementation of coupling two different thermal-hydraulic codes: the system thermal-hydraulics (STH) code System Analysis Module (SAM) and the coarse-mesh computational fluid dynamics (CFD) code Pronghorn. Both applications are included in the Comprehensive Reactor Analysis Bundle (BlueCRAB) code suite and are based on the Multiphysics Object-Oriented Simulation Environment (MOOSE) framework. The domain overlapping (DO) coupling approach was adopted, as it offers proven advantages over more conventional domain decomposition methods. In the DO coupling, SAM provides Pronghorn with boundary conditions that depend on the system-level simulation of the entire plant. In return, the overlapping coupled SAM components–termed “surrogate components”–receive friction factors and source terms computed online based on the Pronghorn simulation. The framework is designed to be generic and enable coupling regardless of geometry and the number of inlet/outlet boundaries in the DO coupled domain. The developed method leads to consistent pressure drops, enthalpies, and scalar concentrations when comparing the coupled SAM and Pronghorn simulations. This document presents the DO coupling approach, along with two verification and two demonstration cases. The proposed problems explore different physical aspects relevant to nuclear reactor analysis, including buoyancy-driven flows, complex flow patterns, and multiple inlets and outlets. Periodically, new versions of this “BlueCRAB Domain Overlapping Coupling: Theory & Verification” report will be issued to reflect future developments and verification tests.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗

Aqueous Diels–Alder reactions for thermochemical storage and heat transfer fluids identified using density functional theory

Thermal storage and transfer fluids have important applications in industrial, transportation, and domestic settings. Current thermal fluids have relatively low specific heats, often significantly below that of water. However, by introducing a thermochemical reaction to a base fluid, it is possible to enhance the fluid's thermal properties. Here, density functional theory (DFT) is used to screen Diels–Alder reactions for use in aqueous thermal fluids. From an initial set of 52 reactions, four are identified with moderate aqueous solubility and predicted turning temperature near the liquid region of water. These reactions are selectively modified through 60 total functional group substitutions to produce novel reactions with improved solubility and thermal properties. Among the reactions generated by functional group substitution, seven have promising predicted thermal properties, significantly improving specific heat (by as much as 30.5%) and energy storage density (by as much as 4.9%) compared to pure water.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Spin polarization induced by magnetic field and the relativistic Barnett effect

First, I study the analogy between the magnetization of a material and the spin polarization of particles in a fluid. Using the relativistic version of the Barnett effect, i.e. the magnetization of a material induced by mechanical rotation, the spin polarization induced by thermal vorticity is obtained within a purely classical model, where spin is treated as an intrinsic magnetic moment and rotation is included as a non-inertial effect. I argue that since spin polarization induced by thermal vorticity can be obtained in a classical theory, it can not be dominated by quantum anomalies. Second, the spin polarization induced by magnetic field is obtained for a fluid at local thermal equilibrium using statistical quantum field theory. The obtained formula is valid beyond the weak field approximation and when contributions from the non-homogeneity of the magnetic field are small. The exact form of spin polarization is studied for a free Dirac field at global equilibrium, and, like magnetic susceptibility, it oscillates according to the de Haas - van Alphen effect. Finally, I briefly review how magnetic field contributes to the difference between the spin polarization of Λ and $\barΛ$ observed in heavy-ion collisions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Collaborative Research: Explosive reconnection in relativistic magnetically-dominated plasmas (Final Report)

Magnetic fields can be the dominant component of astrophysical plasmas, so that the magnetic energy density might exceed even the rest-mass energy density of matter. In this extreme (and largely unexplored) regime the magnetic field controls the overall plasma evolution, dissipation, and acceleration of non-thermal particles. This plasma regime, applicable to a variety of astrophysical sources - magnetars, pulsars and pulsar wind nebulae (PWNe), jets of Active Galactic Nuclei (AGNs) and Gamma-Ray Bursters (GRBs) - is dramatically different from laboratory plasmas, the magnetospheres of planets, and the interplanetary plasma. Relativistic astrophysical sources then provide an unique opportunity to study the fundamental plasma physics of magnetically-dominated plasmas; a novel and fast-evolving field of theoretical research which, by investigating energy conversion and particle energization processes in plasmas, is of vital importance to the Fusion Energy Sciences DoE program. Data coming from astrophysical high-energy missions, especially the Crab Nebula flares recently observed by the Fermi and AGILE satellites, suggest that the acceleration of non-thermal particles to the highest energies occurs in magnetic reconnection events - a major change of paradigm in high-energy astrophysics. Most importantly, observations demand that particle acceleration should proceed extremely fast (with accelerating electric field of the order of the magnetic field) and on macroscopic scales (much larger, e.g., than the microscopic plasma skin depth). We are conducting studies of the microphysics of magnetically-dominated plasmas focussing in particular on the highly dynamic regime of explosive reconnection and associated particle acceleration in relativistic plasmas. We are studying the stability and explosive plasma dynamics, particle acceleration and radiation production in a number of idealized plasma configurations that approximate relevant astrophysical sources (like the magnetic ABC structures and interacting flux tubes, as well as generalizations of analytical models of X-point collapse to relativistic plasmas). The well-studied case of the Crab Nebula is taken as a prototypical example for the application of the model. We are combining analytical studies of explosive magnetic dynamics and dissipation in relativistic plasmas with particle-in-cell (PIC) simulations and fluid simulations. The theoretical model, fluid and particle-in-cell simulations are cross-checked for agreement and convergence.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Hydrodynamics of thermally driven chiral suspensions

Considerable effort has been directed towards the characterization of chiral mesoscale structures, as shown in chiral protein assemblies and carbon nanotubes. Here, we establish a thermally driven hydrodynamic description for the actuation and separation of mesoscale chiral structures in a fluid medium. Cross-flow of a Newtonian liquid with a thermal gradient gives rise to an effective torque that propels each particle of a chiral suspension according to its handedness. In turn, the chiral suspension alters the liquid flow, which thus acquires a transverse (chiral) velocity component. Since observation of the predicted effects requires a low degree of sophistication, our work provides an efficient and inexpensive approach to test and calibrate chiral particle propulsion and separation strategies.

Mechanics↗

Effect of anisotropic fast ions on internal kink stability in DIII-D negative and positive triangularity plasmas

Recent DIII-D experiments show that sawtooth stability is strongly affected by anisotropic fast ions from neutral beam injection (NBI) in both negative and positive triangularity plasmas. Fast ions from co-current NBI are stabilizing for the sawtooth stability, resulting in longer sawtooth periods. On the other hand, fast ions from counter-current NBI are destabilizing, leading to small and frequent sawteeth. The relative change of sawtooth period and amplitude is more than a factor of two. These observations appear to hold in both plasma shapes. Non-perturbative toroidal modeling, utilizing the magnetohydrodynamic-kinetic hybrid stability code MARS-K (Liu et al 2008 Phys. Plasmas 15 112503), reveals an asymmetric dependence of the stability of the n = 1 (n is the toroidal mode number) internal kink mode on the injection direction of NBI, being qualitatively consistent with the experimentally observed sawtooth behavior. The MARS-K modeling results suggest that anisotropic fast ions affect the mode growth rate and frequency through both adiabatic and non-adiabatic contributions. Here, the asymmetry of the internal kink mode instability relative to the NBI direction is mainly due to the non-adiabatic contribution of passing fast ions, which stabilize (destabilize) the internal kink with the co-(counter-) current NBI as compared to the fluid counterpart. However, finite orbit width (FOW) correction to passing particles partially cancels the asymmetry. Trapped particles are always stabilizing due to precessional drift resonance. Modeling also shows that fast ions affect the internal kink in a similar manner in both negative and positive triangularity plasmas, although being slightly more unstable in the negative triangularity configuration already in the fluid limit. The similarity is mainly attributed to the fact that the mode is localized in the plasma core region, with very similar eigenmode structures in both negative and positive configurations. Furthermore, MARS-K modeling indicates that other factors, such as the plasma rotation and the drift kinetic effects of thermal plasmas, weakly modify the mode stability as compared to the drift kinetic resonance effects and FOW correction of fast ions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗