38,859 research outputs found
Comparisons and Applications of Four Independent Numerical Approaches for Linear Gyrokinetic Drift Modes
To help reveal the complete picture of linear kinetic drift modes, four
independent numerical approaches, based on integral equation, Euler initial
value simulation, Euler matrix eigenvalue solution and Lagrangian particle
simulation, respectively, are used to solve the linear gyrokinetic
electrostatic drift modes equation in Z-pinch with slab simplification and in
tokamak with ballooning space coordinate. We identify that these approaches can
yield the same solution with the difference smaller than 1\%, and the
discrepancies mainly come from the numerical convergence, which is the first
detailed benchmark of four independent numerical approaches for gyrokinetic
linear drift modes. Using these approaches, we find that the entropy mode and
interchange mode are on the same branch in Z-pinch, and the entropy mode can
have both electron and ion branches. And, at strong gradient, more than one
eigenstate of the ion temperature gradient mode (ITG) can be unstable and the
most unstable one can be on non-ground eigenstates. The propagation of ITGs
from ion to electron diamagnetic direction at strong gradient is also observed,
which implies that the propagation direction is not a decisive criterion for
the experimental diagnosis of turbulent mode at the edge plasmas.Comment: 12 pages, 10 figures, accept by Physics of Plasma
Numerical Simulation of the Interaction between a Planar Shock Wave and a Cylindrical Bubble
Three-dimensional (3D) computational fluid dynamics (CFD) simulations have been carried out to investigate the complex interaction of a planar shock wave (Ma = 1.22) with a cylindrical bubble. The unsteady Reynolds-averaged Navier–Stokes (URANS) approach with a level set coupled with volume of fluid (LSVOF) method has been applied in the present study. The predicted velocities of refracted wave, transmitted wave, upstream interface, downstream interface, jet, and vortex filaments are in very good agreement with the experimental data. The predicted non-dimensional bubble and vortex velocities also have great concordance with the experimental data compared with a simple model of shock-induced Rayleigh–Taylor instability (i.e., Richtmyer–Meshkov instability) and other theoretical models. The simulated changes in the bubble shape and size (length and width) against time agree very well with the experimental results. Comprehensive flow analysis has shown the shock–bubble interaction (SBI) process clearly from the onset of bubble compression up to the formation of vortex filaments, especially elucidating the mechanism on the air–jet formation and its development. It is demonstrated for the first time that turbulence is generated at the early phase of the shock cylindrical bubble interaction process, with the maximum turbulence intensity reaching about 20% around the vortex filament regions at the later phase of the interaction process
Planetesimal collisions in binary systems
We study the collisional evolution of km-sized planetesimals in tight binary
star systems to investigate whether accretion towards protoplanets can proceed
despite the strong gravitational perturbations from the secondary star. The
orbits of planetesimals are numerically integrated in two dimensions under the
influence of the two stars and gas drag. The masses and orbits of the
planetesimals are allowed to evolve due to collisions with other planetesimals
and accretion of collisional debris. In addition, the mass in debris can evolve
due to planetesimal-planetesimal collisions and the creation of new
planetesimals. We show that it is possible in principle for km-sized
planetesimals to grow by two orders of magnitude in size if the efficiency of
planetesimal formation is relatively low. We discuss the limitations of our
two-dimensional approach.Comment: 5 pages, 5 figures, accepted for publication in MNRA
- …