Fluid Moment Research Articles

Longitudinal sound waves in a collisionless, quasineutral plasma

The time evolution of slow sound waves in a homogeneous, collisionless and quasineutral plasma, in particular their Landau damping, is investigated using the kinetic-magnetohydrodynamics formulation of Ramos (J. Plasma Phys. vol. 81, 2015 p. 905810325; vol. 82, 2016 p. 905820607). In this approach, the electric field is eliminated from a closed, hybrid fluid-kinetic system that ensures automatically the fulfilment of the charge neutrality condition. Considering the time dependence of a spatial-Fourier-mode linear perturbation with wavevector parallel to the equilibrium magnetic field, this can be cast as a second-order self-adjoint problem with a continuum spectrum of real and positive squared frequencies. Therefore, a conventional resolution of the identity with a continuum basis of singular normal modes is guaranteed, which simplifies significantly a Van Kampen-like treatment of the Landau damping. The explicit form of such singular normal modes is obtained, along with their orthogonality relations. These are used to derive the damped time evolution of the fluid moments of solutions of initial-value problems, for the most general kinds of initial conditions. The non-zero parallel electric field is not used explicitly in this analysis, but it is calculated from any given solution after the later has been obtained.

Bifurcation Analysis of a Rotor Supported by Five-Pad Tilting Pad Journal Bearings Using Numerical Continuation

This paper presents analytical bifurcations analysis of a “Jeffcott” type rigid rotor supported by five-pad tilting pad journal bearings (TPJBs). Numerical techniques such as nonautonomous shooting/arc-length continuation, Floquet theory, and Lyapunov exponents are employed along with direct numerical integration (NI) to analyze nonlinear characteristics of the TPJB-rotor system. A rocker pivot type five-pad TPJB is modeled with finite elements to evaluate the fluid pressure distribution on the pads, and the integrated fluid reaction force and moment are utilized to determine coexistent periodic solutions and bifurcations scenarios. The numerical shooting/continuation algorithms demand significant computational workload when applied to a rotor supported by a finite element bearing model. This bearing model may be significantly more accurate than the simplified infinitely short-/long-bearing approximations. Consequently, the use of efficient computation techniques such as deflation and parallel computing methods is applied to reduce the execution time. Loci of bifurcations of the TPJB-rigid rotor are determined with extensive numerical simulations with respect to both rotor spin speed and unbalance force magnitude. The results show that heavily loaded bearings and/or high unbalance force may induce consecutive transference of response in forms of synchronous to subsynchronous, quasi-periodic responses, and chaotic motions. It is revealed that the coexistent responses and their solution manifolds are obtainable and stretch out with selections of pad preload, pivot offset, and lubricant viscosity so that the periodic doubling bifurcations, saddle node bifurcations, and corresponding local stability are reliably determined by searching parameter sets. In case the system undergoes an aperiodic state, the rate of divergence/convergence of the attractor is examined quantitatively by using the maximum Lyapunov exponent (MLE).

Effects of pivot design on performance of tilting pad journal bearings

Abstract This paper emphasizes on the effects of pivot flexibility upon the performance of a TPJB by making comparison between two model predictions i.e. (a) model with rigid pivots, and (b) model with soft pivots. The latter one uses Hertzian contact theory for the calculation of the pivot stiffness and couples it with the fluid film hydrodynamics. Presently, finite element method solves pressure distribution which integrates to give fluid film forces and moments. Predicted static and dynamic characteristics using two different models illuminate the dominance of pivot flexibility effects on the performance of the TPJB. In general, soft pivot model predictions are closely predicted when compared to the previously published test data.

Simulations of anti-parallel reconnection using a nonlocal heat flux closure

The integration of kinetic effects in fluid models is important for global simulations of the Earth's magnetosphere. In particular, it has been shown that ion kinetics play a crucial role in the dynamics of large reconnecting systems, and that higher-order fluid moment models can account for some of these effects. Here, we use a ten-moment model for electrons and ions, which includes the off diagonal elements of the pressure tensor that are important for magnetic reconnection. Kinetic effects are recovered by using a nonlocal heat flux closure, which approximates linear Landau damping in the fluid framework. The closure is tested using the island coalescence problem, which is sensitive to ion dynamics. We demonstrate that the nonlocal closure is able to self-consistently reproduce the structure of the ion diffusion region, pressure tensor, and ion velocity without the need for fine-tuning of relaxation coefficients present in earlier models.

Dynamic Coupling of Boltzmann Plasma Model to Surface Erosion Model for Kinetic Treatment of Plasma-Material Interactions

We present an innovative coupled Boltzmann–binary collision approximation (BCA) method for the simulation of the near-wall plasma in the presence of a material-releasing wall. The method is based on a full-f multispecies Boltzmann solver for the plasma (charged and neutral species) coupled to a modification of the classical BCA code TRIDYN. Both the plasma ions and the impurities are treated as Boltzmann kinetic species, allowing high resolution even at very disparate densities, particle fluxes, drift velocities, and energy fluxes. From the distribution functions, all the fluid moments (density, heat flux, etc.) and the net and gross erosion rates are derived. An example of calculation of a helium plasma facing a beryllium wall is reported, showing the evolution of the phase-spaces of ions, neutrals, and material impurities in the near-wall region at nominal ITER conditions.

Heat transfer analysis in an annular cone subjected to power law variations

Present study deals with the analysis of heat transfer and fluid flow behavior in an annular cone fixed with saturated porous medium. The inner surface of the cone is assumed to have power law variable wall temperature. The governing partial differential equations are solved using well known Finite Element Method (FEM). The coupled nonlinear differential equations are converted into the algebraic equations by using Galerkin method. A 3 noded triangular element is used to divide the porous domain into smaller segments. The effects of various geometrical parameters on the cone angle are presented. It is found that the effect of cone angle on the heat transfer characteristics and fluid flow behavior is considerably significant. The fluid moment is found to shift towards the upper side of cone with increase in the power law coefficient. The fluid velocity decreases with increase in the power law coefficient.

Numerical investigation of flow-induced rotary oscillation of circular cylinder with rigid splitter plate

Numerical results of fluid flow over a rotationally oscillating circular cylinder with splitter plate are presented here. Different from the previous examinations with freely rotatable assembly, the fluid and structure interactions are treated as a coupled dynamic system by fully considering the structural inertia, stiffness, and damping. The hydrodynamic characteristics are examined in terms of reduced velocity Ur at a relatively low Reynolds number Re = 100 for different plate lengths of L/D = 0.5, 1.0, and 1.5, where Ur = U/(Dfn), Re = UD/υ and fn = (κ/J)0.5/2π with U the free stream velocity, D the diameter of the circular cylinder, υ the fluid kinematic viscosity, fn the natural frequency, J the inertial moment, κ the torsional stiffness, and L the plate length. Contrast to the freely rotating cylinder/plate body, that is, in the limit of κ → 0 or Ur →∞, remarkable rotary oscillation is observed at relatively low reduced velocities. For the typical case with L/D = 1.0, the maximum amplitude may reach five times that at the highest reduced velocity of Ur = 15.0 considered in this work. At the critical reduced velocity Ur = 4.2, notable hydrodynamic jumps are identified for the rotation amplitude, response frequency, mean drag coefficient, lift amplitude, and vortex shedding frequency. Moreover, the phase angle between the fluid moment and rotary oscillation abruptly changes from 0 to π at Ur = 6.5. Due to the combined effect of fluid moment, rotation response, and phase difference, the natural frequency of the rotating body varies in flow, leading to a wide regime of lock-in/synchronization (Ur ≥4.2, for L/D = 1.0). The phenomenon of rotation bifurcation, i.e., the equilibrium position of the rotary oscillation deflects to a position which is not parallel to the free stream, is found to only occur at higher reduced velocities. The longer splitter plate has the lower critical reduced velocity. The occurrence of bifurcation is attributed to the anti-symmetry breaking of the wake flow evolution. The resultant asymmetric mean pressure distribution on the splitter plate gives rise to the net lift force and the deviated moment on the assembly, leading to the offset mean position of splitter plate. The global vortex shedding is identified to be the classic 2S mode for both cases with and without the bifurcation, although the second vortex formation and the shedding pattern in the near wake for the bifurcate case are different from the non-bifurcate case with lower reduced velocities.

Development of the fluid-type transport code on the flux coordinates in a tokamak

The one-dimensional fluid-type transport code, TASK/TX, is developed compatible with the flux coordinates in a tokamak. Unlike diffusive transport equations usually adopted in conventional transport codes, the governing equations conform to a two-fluid model consisting of Maxwell’s equations and the multiple fluid moment equations for each species.Quasi-neutrality and ambipolar flux conditions are not imposed, which are inherently satisfied as a consequence of the equation system solved. The neoclassical particle flux is not approximated by the flux–gradient relationship, and the total particle flux composed of the neoclassical and turbulent contributions is directly treated as the dependent variable. The quantities related to neoclassical transport are intrinsically calculated without external neoclassical transport modules. In other words, TASK/TX by itself has the function of a neoclassical transport solver based on the moment approach as well. Several numerical tests clearly reveal the unique features of TASK/TX not possessed by conventional transport codes.

Fluid moments of the nonlinear Landau collision operator

An important problem in plasma physics is the lack of an accurate and complete description of Coulomb collisions in associated fluid models. To shed light on the problem, this Letter introduces an integral identity involving the multivariate Hermite tensor polynomials and presents a method for computing exact expressions for the fluid moments of the nonlinear Landau collision operator. The proposed methodology provides a systematic and rigorous means of extending the validity of fluid models that have an underlying inverse-square force particle dynamics to arbitrary collisionality and flow.

Connecting Collisionless Landau Fluid Closures to Collisional Plasma Physics Models

AbstractIn this work, collisionless Landau fluid models are combined with collisional fluid models using a Padé approximation that is accurate in both collisionality limits. The collisionless models capture Landau damping through a nonlocal integro‐differential sub‐diffusion operator with ballistic characteristics. Collisional extensions of Landau fluid models are derived by analyzing the higher order moment equations which combine a Landau closure with collisional friction forces. The model derived here evolves fluid moments for density, parallel velocity, and anisotropic pressure and includes the frictional heat flux, the parallel thermal force and anisotropic electrical conductivity. Since anisotropies must vanish in the collisional limit, a simple closure can be derived if the friction force neglects pitch‐angle scattering for the closure moments themselves. The resulting plasma physics model is potentially quite useful for applications in magnetic fusion and astrophysics. (© 2016 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)

Suppression of phase mixing in drift-kinetic plasma turbulence

Transfer of free energy from large to small velocity-space scales by phase mixing leads to Landau damping in a linear plasma. In a turbulent drift-kinetic plasma, this transfer is statistically nearly canceled by an inverse transfer from small to large velocity-space scales due to “anti-phase-mixing” modes excited by a stochastic form of plasma echo. Fluid moments (density, velocity, and temperature) are thus approximately energetically isolated from the higher moments of the distribution function, so phase mixing is ineffective as a dissipation mechanism when the plasma collisionality is small.

Moment tensors of a dislocation in a porous medium

A dislocation can be represented by a moment tensor for calculating seismic waves. However, the moment tensor expression was derived in an elastic medium and cannot completely describe a dislocation in a porous medium. In this paper, effective moment tensors of a dislocation in a porous medium are derived. It is found that the dislocation is equivalent to two independent moment tensors, i.e., the bulk moment tensor acting on the bulk of the porous medium and the isotropic fluid moment tensor acting on the pore fluid. Both of them are caused by the solid dislocation as well as the fluid–solid relative motion corresponding to fluid injection towards the surrounding rocks (or fluid outflow) through the fault plane. For a shear dislocation, the fluid moment tensor is zero, and the dislocation is equivalent to a double couple acting on the bulk; for an opening dislocation or fluid injection, the two moment tensors are needed to describe the source. The fluid moment tensor only affects the radiated compressional waves. By calculating the ratio of the radiation fields generated by unit fluid moment tensor and bulk moment tensor, it is found that the fast compressional wave radiated by the bulk moment tensor is much stronger than that radiated by the fluid moment tensor, while the slow compressional wave radiated by the fluid moment tensor is several times stronger than that radiated by the bulk moment tensor.

Analysis of heat transfer due to stretching cylinder with partial slip and prescribed heat flux: A Chebyshev Spectral Newton Iterative Scheme

This study is dedicated to analyze the combined effects of partial slip and prescribed surface heat flux when the fluid is moving due to stretching cylinder. A very moderate and powerful technique Chebyshev Spectral Newton Iterative Scheme is used to determine the solution of the present mathematical model. Involved physical parameters, namely the slip parameter, Casson fluid parameter, curvature parameter and Prandtl number are utilized to control the fluid moments and temperature distribution. The results show that the fluid velocity and the skin friction coefficient on the stretching cylinder are strongly influenced by the slip parameter. It is further analyzed that hydrodynamic boundary layer decreases and thermal boundary layer increases with the slip parameter. Influence of Casson fluid parameter on temperature profile provides the opposite behavior as compared to the slip parameter. The comparison of numerical values of skin friction coefficient and the local Nusselt number is made with the results available in the literature. The accuracy and convergence of Chebyshev Spectral Newton Iterative Scheme is compared with finite difference scheme (Keller box method) through tables. The CPU time is calculated for both schemes. It is observed that CSNIS is efficient, less time consuming, stable and rapid convergent.

Hamiltonian fluid reductions of electromagnetic drift-kinetic equations for an arbitrary number of moments

We present an infinite family of Hamiltonian electromagnetic fluid models for plasmas, derived from drift-kinetic equations. An infinite hierarchy of fluid equations is obtained from a Hamiltonian drift-kinetic system by taking moments of a generalized distribution function and using Hermite polynomials as weight functions of the velocity coordinate along the magnetic guide field. Each fluid model is then obtained by truncating the hierarchy to a finite number N+1 of equations by means of a closure relation. We show that, for any positive N, a linear closure relation between the moment of order N+1 and the moment of order N guarantees that the resulting fluid model possesses a Hamiltonian structure, thus respecting the Hamiltonian character of the parent drift-kinetic model. An orthogonal transformation is identified which maps the fluid moments to a new set of dynamical variables in terms of which the Poisson brackets of the fluid models become a direct sum and which unveils remarkable dynamical properties of the models in the two-dimensional (2D) limit. Indeed, when imposing translational symmetry with respect to the direction of the magnetic guide field, all models belonging to the infinite family can be reformulated as systems of advection equations for Lagrangian invariants transported by incompressible generalized velocities. These are reminiscent of the advection properties of the parent drift-kinetic model in the 2D limit and are related to the Casimirs of the Poisson brackets of the fluid models. The Hamiltonian structure of the generic fluid model belonging to the infinite family is illustrated treating a specific example of a fluid model retaining five moments in the electron dynamics and two in the ion dynamics. We also clarify the connection existing between the fluid models of this infinite family and some fluid models already present in the literature.

Extension of the flow-rate-of-strain tensor formulation of plasma rotation theory to non-axisymmetric tokamaks

A systematic formalism for the calculation of rotation in non-axisymmetric tokamaks with 3D magnetic fields is described. The Braginskii Ωτ-ordered viscous stress tensor formalism, generalized to accommodate non-axisymmetric 3D magnetic fields in general toroidal flux surface geometry, and the resulting fluid moment equations provide a systematic formalism for the calculation of toroidal and poloidal rotation and radial ion flow in tokamaks in the presence of various non-axisymmetric “neoclassical toroidal viscosity” mechanisms. The relation among rotation velocities, radial ion particle flux, ion orbit loss, and radial electric field is discussed, and the possibility of controlling these quantities by producing externally controllable toroidal and/or poloidal currents in the edge plasma for this purpose is suggested for future investigation.

Extensions and applications of the Bohm criterion

The generalized Bohm criterion is revisited in the context of incorporating kinetic effects of the electron and ion distribution functions into the theory. The underlying assumptions and results of two different approaches are compared: the conventional ‘kinetic Bohm criterion’ and a fluid-moment hierarchy approach. The former is based on the asymptotic limit of an infinitely thin sheath (λD/l = 0), whereas the latter is based on a perturbative expansion of a sheath that is thin compared to the plasma (λD/l ≪ 1). Here λD is the Debye length, which characterizes the sheath length scale, and l is a measure of the plasma or presheath length scale. The consequences of these assumptions are discussed in terms of how they restrict the class of distribution functions to which the resulting criteria can be applied. Two examples are considered to provide concrete comparisons between the two approaches. The first is a Tonks–Langmuir model including a warm ion source (Robertson 2009 Phys. Plasmas 16 103503). This highlights a substantial difference between the conventional kinetic theory, which predicts slow ions dominate at the sheath edge, and the fluid moment approach, which predicts slow ions have little influence. The second example considers planar electrostatic probes biased near the plasma potential using model equations and particle-in-cell simulations. This demonstrates a situation where electron kinetic effects alter the Bohm criterion, leading to a subsonic ion flow at the sheath edge.

On push-forward representations in the standard gyrokinetic model

Two representations of fluid moments in terms of a gyro-center distribution function and gyro-center coordinates, which are called push-forward representations, are compared in the standard electrostatic gyrokinetic model. In the representation conventionally used to derive the gyrokinetic Poisson equation, the pull-back transformation of the gyro-center distribution function contains effects of the gyro-center transformation and therefore electrostatic potential fluctuations, which is described by the Poisson brackets between the distribution function and scalar functions generating the gyro-center transformation. Usually, only the lowest order solution of the generating function at first order is considered to explicitly derive the gyrokinetic Poisson equation. This is true in explicitly deriving representations of scalar fluid moments with polarization terms. One also recovers the particle diamagnetic flux at this order because it is associated with the guiding-center transformation. However, higher-order solutions are needed to derive finite Larmor radius terms of particle flux including the polarization drift flux from the conventional representation. On the other hand, the lowest order solution is sufficient for the other representation, in which the gyro-center transformation part is combined with the guiding-center one and the pull-back transformation of the distribution function does not appear.

1618 バックシュラウドに作用するロータダイナミック流体モーメントによるふれまわり運動の発生、およびロータダイナミック流体モーメントの発生メカニズム

1618 バックシュラウドに作用するロータダイナミック流体モーメントによるふれまわり運動の発生、およびロータダイナミック流体モーメントの発生メカニズム

Moment-Based Acceleration for Neutral Gas Kinetics with BGK Collision Operator

In this work, we present a moment-based accelerator algorithm for a Picard iteration applied to a neutral gas dynamics Boltzmann transport equation with a Bhatnagar-Gross-Krook collision operator. Traditional approaches relying on either explicit or Picard iteration schemes (i.e., source iteration) are severely limited for investigating time-scales much larger than the collisional relaxation time-scale, τ. We have developed a nonlinear accelerator algorithm that allows one to step over this stiff collision time scale and follow the hydrodynamic time scale of the problem when appropriate. The new algorithm relies on formulating a nonlinear, coupled system comprised of a high-order (HO) kinetic equation and a low-order (LO) fluid moment equation system. The HO equation provides self-consistent closures to the LO fluid equations, while the latter provides the required implicit-moment variables to evaluate the collision operator. We characterize the performance of the new algorithm on a Sod shock tube and a strong shock tube problem with varying Knudsen number.

Hamiltonian structure of a drift-kinetic model and Hamiltonian closures for its two-moment fluid reductions

We address the problem of the existence of the Hamiltonian structure for an electrostatic drift-kinetic model and for the related fluid models describing the evolution of the first two moments of the distribution function with respect to the parallel velocity. The drift-kinetic model, which accounts for background density and temperature gradients as well as polarization effects, is shown to possess a noncanonical Hamiltonian structure. The corresponding Poisson bracket is expressed in terms of the fluid moments and it is found that the set of functionals of the zero order moment forms a sub-algebra, thus automatically leading to a class of one-moment Hamiltonian fluid models. In particular, in the limit of weak spatial variations of the background quantities, the Charney-Hasegawa-Mima equation, with its Hamiltonian structure, is recovered. For the set of functionals of the first two moments, which, unlike the case of the Vlasov equation, turns out not to form a sub-algebra, we look for closures that lead to a closed Poisson bracket restricted to this set of functionals. The constraint of the Jacobi identity turns out to select the adiabatic equation of state for an ideal gas with one-degree-of-freedom molecules, as the only admissible closure in this sense. When the so called δf ordering is applied to the model, on the other hand, a Poisson bracket is obtained if the second order moment is a linear combination of the first two moments of the total distribution function. By means of this procedure, three-dimensional Hamiltonian fluid models that couple a generalized Charney-Hasegawa-Mima equation with an evolution equation for the parallel velocity are derived. Among these, a model adopted by Meiss and Horton [Phys. Fluids 26, 990 (1983)] to describe drift waves coupled to ion-acoustic waves, is obtained and its Hamiltonian structure is provided explicitly.