Researcher profile

E. Khomenko

E. Khomenko contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
7works
0followers
3topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

7 published item(s)

preprint2022arXiv

Generalized Fluid Models of the Braginskii Type

Several generalizations of the well-known fluid model of Braginskii (Rev. of Plasma Phys., 1965) are considered. We use the Landau collisional operator and the moment method of Grad. We focus on the 21-moment model that is analogous to the Braginskii model, and we also consider a 22-moment model. Both models are formulated for general multi-species plasmas with arbitrary masses and temperatures, where all the fluid moments are described by their evolution equations. The 21-moment model contains two "heat flux vectors" (3rd and 5th-order moments) and two "viscosity-tensors" (2nd and 4th-order moments). The Braginskii model is then obtained as a particular case of a one ion-electron plasma with similar temperatures, with de-coupled heat fluxes and viscosity-tensors expressed in a quasi-static approximation. We provide all the numerical values of the Braginskii model in a fully analytic form (together with the 4th and 5th-order moments). For multi-species plasmas, the model makes calculation of transport coefficients straightforward. Formulation in fluid moments (instead of Hermite moments) is also suitable for implementation into existing numerical codes. It is emphasized that it is the quasi-static approximation which makes some Braginskii coefficients divergent in a weakly-collisional regime. Importantly, we show that the heat fluxes and viscosity-tensors are coupled even in the linear approximation, and that the fully contracted (scalar) perturbations of the 4th-order moment, which are accounted for in the 22-moment model, modify the energy exchange rates. We also provide several Appendices, which can be useful as a guide for deriving the Braginskii model with the moment method of Grad.

preprint2022arXiv

Modeling the thermal conduction in the solar atmosphere with the code MANCHA3D

Thermal conductivity is one of the important mechanisms of heat transfer in the solar corona. In the limit of strongly magnetized plasma, it is typically modeled by Spitzer's expression where the heat flux is aligned with the magnetic field. This paper describes the implementation of the heat conduction into the code MANCHA3D with an aim of extending single-fluid MHD simulations from the upper convection zone into the solar corona. Two different schemes to model heat conduction are implemented: (1) a standard scheme where a parabolic term is added to the energy equation, and (2) a scheme where the hyperbolic heat flux equation is solved. The first scheme limits the time step due to the explicit integration of a parabolic term, which makes the simulations computationally expensive. The second scheme solves the limitations on the time step by artificially limiting the heat conduction speed to computationally manageable values. The validation of both schemes is carried out with standard tests in one, two, and three spatial dimensions. Furthermore, we implement the model for heat flux derived by Braginskii (1965) in its most general form, when the expression for the heat flux depends on the ratio of the collisional to cyclotron frequencies of the plasma, and, therefore on the magnetic field strength. Additionally, our implementation takes into account the heat conduction in parallel, perpendicular, and transverse directions, and provides the contributions from ions and electrons separately. The model also transitions smoothly between field-aligned conductivity and isotropic conductivity for regions with a low or null magnetic field. Finally, we present a two-dimensional test for heat conduction using realistic values of the solar atmosphere where we prove the robustness of the two schemes implemented.

preprint2020arXiv

Acoustic-gravity wave propagation characteristics in 3D radiation hydrodynamic simulations of the solar atmosphere

There has been tremendous progress in the degree of realism of three-dimensional radiation magneto-hydrodynamic simulations of the solar atmosphere in the past decades. Four of the most frequently used numerical codes are Bifrost, CO5BOLD, MANCHA3D, and MURaM. Here we test and compare the wave propagation characteristics in model runs from these four codes by measuring the dispersion relation of acoustic-gravity waves at various heights. We find considerable differences between the various models. The height dependence of wave power, in particular of high-frequency waves, varies by up to two orders of magnitude between the models, and the phase difference spectra of several models show unexpected features, including $\pm180^\circ$ phase jumps.

preprint2020arXiv

Influence of ambipolar and Hall effects on vorticity in 3D simulations of magneto-convection

This paper presents the results of the analysis of 3D simulations of solar magneto-convection that include the joint action of the ambipolar diffusion and the Hall effect. Three simulation-runs are compared: one including both ambipolar diffusion and Hall effect; one including only ambipolar diffusion; and one without any of these two effects. The magnetic field is amplified from initial field to saturation level by the action of turbulent local dynamo. In each of these cases, we study 2 hours of simulated solar time after the local dynamo reaches the saturation regime. We analyze the power spectra of vorticity, of magnetic field fluctuations and of the different components of the magnetic Poynting flux responsible for the transport of vertical or horizontal perturbations. Our preliminary results show that the ambipolar diffusion produces a strong reduction of vorticity in the upper chromospheric layers and that it dissipates the vortical perturbations converting them into thermal energy. The Hall effect acts in the opposite way, strongly enhancing the vorticity. When the Hall effect is included, the magnetic field in the simulations becomes, on average, more vertical and long-lived flux tube-like structures are produced. We trace a single magnetic structure to study its evolution pattern and the magnetic field intensification, and their possible relation to the Hall effect.

preprint2020arXiv

Two-fluid simulations of Rayleigh-Taylor instability in a magnetized solar prominence thread. I. Effects of prominence magnetization and mass loading

Solar prominences are formed by partially ionized plasma with inter-particle collision frequencies generally warranting magnetohydrodynamic treatment. In this work, we explore the dynamical impacts and observable signatures of two-fluid effects in the parameter regimes when ion-neutral collisions do not fully couple the neutral and charged fluids. We perform 2.5D two-fluid (charges - neutrals) simulations of the Rayleigh-Taylor instability (RTI) at a smoothly changing interface between a solar prominence thread and the corona. The purpose of this study is to deepen our understanding of the RTI and the effects of the partial ionization on the development of RTI using non-linear two-fluid numerical simulations. Our two-fluid model takes into account viscosity, thermal conductivity, and collisional interaction between neutrals and charges: ionization/recombination, energy and momentum transfer, and frictional heating. In this paper I, the sensitivity of the RTI dynamics to the prominence equilibrium configuration, including the impact of the magnetic field strength and shear supporting the prominence thread, and the amount of prominence mass-loading is explored. We show that, at small scales, a realistically smooth prominence-corona interface leads to qualitatively different linear RTI evolution than that expected for a discontinuous interface, while magnetic field shear has the stabilizing effect of reducing the growth rate or eliminating the instability. In the non-linear phase, we observe that in the presence of field shear the development of the instability leads to formation of coherent and interacting 2.5D magnetic structures, which, in turn, can lead to substantial plasma flow across magnetic field lines and associated decoupling of the fluid velocities of charges and neutrals.

preprint2019arXiv

An introductory guide to fluid models with anisotropic temperatures Part 1 -- CGL description and collisionless fluid hierarchy

We present a detailed guide to advanced collisionless fluid models that incorporate kinetic effects into the fluid framework, and that are much closer to the collisionless kinetic description than traditional magnetohydrodynamics. Such fluid models are directly applicable to modeling turbulent evolution of a vast array of astrophysical plasmas, such as the solar corona and the solar wind, the interstellar medium, as well as accretion disks and galaxy clusters. The text can be viewed as a detailed guide to Landau fluid models and it is divided into two parts. Part 1 is dedicated to fluid models that are obtained by closing the fluid hierarchy with simple (non Landau fluid) closures. Part 2 is dedicated to Landau fluid closures. Here in Part 1, we discuss the CGL fluid model in great detail, together with fluid models that contain dispersive effects introduced by the Hall term and by the finite Larmor radius (FLR) corrections to the pressure tensor. We consider dispersive effects introduced by the non-gyrotropic heat flux vectors. We investigate the parallel and oblique firehose instability, and show that the non-gyrotropic heat flux strongly influences the maximum growth rate of these instabilities. Furthermore, we discuss fluid models that contain evolution equations for the gyrotropic heat flux fluctuations and that are closed at the 4th-moment level by prescribing a specific form for the distribution function. For the bi-Maxwellian distribution, such a closure is known as the "normal" closure. We also discuss a fluid closure for the bi-kappa distribution. Finally, by considering one-dimensional Maxwellian fluid closures at higher-order moments, we show that such fluid models are always unstable. The last possible non Landau fluid closure is therefore the "normal" closure, and beyond the 4th-order moment, Landau fluid closures are required.

preprint2019arXiv

An introductory guide to fluid models with anisotropic temperatures Part 2 -- Kinetic theory, Padé approximants and Landau fluid closures

In Part 2 of our guide to collisionless fluid models, we concentrate on Landau fluid closures. These closures were pioneered by Hammett and Perkins and allow for the rigorous incorporation of collisionless Landau damping into a fluid framework. It is Landau damping that sharply separates traditional fluid models and collisionless kinetic theory, and is the main reason why the usual fluid models do not converge to the kinetic description, even in the long-wavelength low-frequency limit. We start with a brief introduction to kinetic theory, where we discuss in detail the plasma dispersion function $Z(ζ)$, and the associated plasma response function $R(ζ)=1+ζZ(ζ)=-Z'(ζ)/2$. We then consider a 1D (electrostatic) geometry and make a significant effort to map all possible Landau fluid closures that can be constructed at the 4th-order moment level. These closures for parallel moments have general validity from the largest astrophysical scales down to the Debye length, and we verify their validity by considering examples of the (proton and electron) Landau damping of the ion-acoustic mode, and the electron Landau damping of the Langmuir mode. We proceed by considering 1D closures at higher-order moments than the 4th-order, and as was concluded in Part 1, this is not possible without Landau fluid closures. We show that it is possible to reproduce linear Landau damping in the fluid framework to any desired precision, thus showing the convergence of the fluid and collisionless kinetic descriptions. We then consider a 3D (electromagnetic) geometry in the gyrotropic (long-wavelength low-frequency) limit and map all closures that are available at the 4th-order moment level. In the Appendix A, we provide comprehensive tables with Padé approximants of $R(ζ)$ up to the 8th-pole order, with many given in an analytic form.