BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Ammar Hakim

Ammar Hakim contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
physics.plasm-phphysics.comp-phMachine Learningmath.NANumerical Analysisphysics.flu-dyn

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
13works
0followers
6topics
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

13 published item(s)

preprint2022arXiv

A kinetic interpretation of the classical Rayleigh-Taylor instability

Rayleigh-Taylor (RT) instabilities are prevalent in many physical regimes ranging from astrophysical to laboratory plasmas and have primarily been studied using fluid models, the majority of which have been ideal fluid models. This work is the first of its kind to present a 5-dimensional (2 spatial dimensions, 3 velocity space dimensions) simulation using the continuum-kinetic model to study the effect of the collisional mean-free-path and transport on the instability growth. The continuum-kinetic model provides noise-free access to the full particle distribution function permitting a detailed investigation of the role of kinetic physics in hydrodynamic phenomena such as the RT instability. For long mean-free-path, there is no RT instability growth, but as collisionality increases, particles relax towards the Maxwellian velocity distribution, and the kinetic simulations reproduce the fluid simulation results. An important and novel contribution of this work is in the intermediate collisional cases that are not accessible with traditional fluid models and require kinetic modeling. Simulations of intermediate collisional cases show that the RT instability evolution is significantly altered compared to the highly collisional fluid-like cases. Specifically, the growth rate of the intermediate collisionality RT instability is lower than the high collisionality case while also producing a significantly more diffused interface. The higher moments of the distribution function play a more significant role relative to inertial terms for intermediate collisionality during the evolution of the RT instability interface. Particle energy-flux is calculated from moments of the distribution and shows that transport is significantly altered in the intermediate collisional case and deviates much more so from the high collisionality limit of the fluid regime.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Electron cyclotron drift instability and anomalous transport: two-fluid moment theory and modeling

In the presence of a strong electric field perpendicular to the magnetic field, the electron cross-field (E$\times$B) flow relative to the unmagnetized ions can cause the Electron Cyclotron Drift Instability (ECDI) due to resonances of the ion acoustic mode and the electron cyclotron harmonics. This occurs in collisionless shock ramps in space, and in $\rm{E \times B}$ discharge devices such as Hall thrusters. ECDI can induce an electron flow parallel to the background E field at a speed greatly exceeding predictions by classical collision theory. Such anomalous transport may lead to particle thermalization at space shocks, and may cause unfavorable plasma flows towards the walls of E$\times$B devices. The development of ECDI and anomalous transport is often considered fully-kinetic. In this work, however, we demonstrate that a reduced variant of this instability, and more importantly, the associated anomalous transport, can be treated self-consistently in a two-fluid framework without any collision. By treating electrons and ions on an equal footing, the free energy allows the growth of an anomalous electron flow parallel to the background E field. We first present linear analyses of the instability in the two-fluid 5- and 10-moment models, and compare them against the fully-kinetic theory. At lower temperatures, the two-fluid fastest-growing mode is in good agreement with the kinetic result. Also, by including more ($>=10$) moments, secondary (and possibly higher) unstable branches can be recovered. The dependence of the instability on various parameters is also explored. We then carry out direct numerical simulations of the cross-field setup using the 5-moment model. The growth of the instability and the anomalous transport is confirmed. Finally, 5-moment and Vlasov simulations using identical parameters in the lower-temperature regime are performed, showing reasonable agreement.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Improved multispecies Dougherty collisions

The Dougherty model Fokker-Planck operator is extended to describe nonlinear full-f collisions between multiple species in plasmas. Simple relations for cross-species interactions are developed which obey conservation laws, and reproduce familiar velocity and temperature relaxation rates. This treatment of multispecies Dougherty collisions, valid for arbitrary mass ratios, satisfies the H-Theorem unlike analogous Bhatnagar-Gross-Krook operators.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Turbulent field fluctuations in gyrokinetic and fluid plasmas

A key uncertainty in the design and development of magnetic confinement fusion energy reactors is predicting edge plasma turbulence. An essential step in overcoming this uncertainty is the validation in accuracy of reduced turbulent transport models. Drift-reduced Braginskii two-fluid theory is one such set of reduced equations that has for decades simulated boundary plasmas in experiment, but significant questions exist regarding its predictive ability. To this end, using a novel physics-informed deep learning framework, we demonstrate the first ever direct quantitative comparisons of turbulent field fluctuations between electrostatic two-fluid theory and electromagnetic gyrokinetic modelling with good overall agreement found in magnetized helical plasmas at low normalized pressure. This framework is readily adaptable to experimental and astrophysical environments, and presents a new technique for the numerical validation and discovery of reduced global plasma turbulence models.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A boundary value "reservoir problem" and boundary conditions for multi-moment multifluid simulations of sheaths

Multifluid simulations of plasma sheaths are increasingly used to model a wide variety of problems in plasma physics ranging from global magnetospheric flows around celestial bodies to plasma-wall interactions in thrusters and fusion devices. For multifluid problems, accurate boundary conditions to model an absorbing wall that resolves a classical sheath remains an open research area. This work justifies the use of vacuum boundary conditions for absorbing walls to show comparable accuracy between a multifluid sheath and lower moments of a continuum-kinetic sheath.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Alias-free, matrix-free, and quadrature-free discontinuous Galerkin algorithms for (plasma) kinetic equations

Understanding fundamental kinetic processes is important for many problems, from plasma physics to gas dynamics. A first-principles approach to these problems requires a statistical description via the Boltzmann equation, coupled to appropriate field equations. In this paper we present a novel version of the discontinuous Galerkin (DG) algorithm to solve such kinetic equations. Unlike Monte-Carlo methods we use a continuum scheme in which we directly discretize the 6D phase-space using discontinuous basis functions. Our DG scheme eliminates counting noise and aliasing errors that would otherwise contaminate the delicate field-particle interactions. We use modal basis functions with reduced degrees of freedom to improve efficiency while retaining a high formal order of convergence. Our implementation incorporates a number of software innovations: use of JIT compiled top-level language, automatically generated computational kernels and a sophisticated shared-memory MPI implementation to handle velocity space parallelization.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Conservative discontinuous Galerkin scheme of a gyro-averaged Dougherty collision operator

A conservative discontinuous Galerkin scheme for a nonlinear Dougherty collision operator in full-f long-wavelength gyrokinetics is presented. Analytically this model operator has the advective-diffusive form of Fokker-Planck operators, it has a non-decreasing entropy functional, and conserves particles, momentum and energy. Discretely these conservative properties are maintained exactly as well, independent of numerical resolution. In this work the phase space discretization is performed using a novel version of the discontinuous Galerkin scheme, carefully constructed using concepts of weak equality and recovery. Discrete time advancement is carried out with an explicit time-stepping algorithm, whose stability limits we explore. The formulation and implementation within the long-wavelength gyrokinetic solver of Gkeyll are validated with relaxation tests, collisional Landau-damping benchmarks and the study of 5D gyrokinetic turbulence on helical, open field lines.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Dynamo in weakly collisional nonmagnetized plasmas impeded by Landau damping of magnetic fields

We perform fully kinetic simulations of flows known to produce dynamo in magnetohydrodynamics (MHD), considering scenarios with low Reynolds number and high magnetic Prandtl number, relevant for galaxy cluster scale fluctuation dynamos. We find that Landau damping on the electrons leads to a rapid decay of magnetic perturbations, impeding the dynamo. This collisionless damping process operates on spatial scales where electrons are nonmagnetized, reducing the range of scales where the magnetic field grows in high magnetic Prandtl number fluctuation dynamos. When electrons are not magnetized down to the resistive scale, the magnetic energy spectrum is expected to be limited by the scale corresponding to magnetic Landau damping or, if smaller, the electron gyroradius scale, instead of the resistive scale. In simulations we thus observe decaying magnetic fields where resistive MHD would predict a dynamo.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Fluid & Gyrokinetic turbulence in open field-line, helical plasmas

Two-fluid Braginskii codes have simulated open-field line turbulence for over a decade, and only recently has it become possible to study these systems with continuum gyrokinetic codes. This work presents a first-of-its-kind comparison between fluid and (long-wavelength) gyrokinetic models in open field-lines, using the GDB and Gkeyll codes to simulate interchange turbulence in the Helimak device at the University of Texas (T. N. Bernard, et. al., Phys. of Plasmas 26, 042301 (2019)). Partial agreement is attained in a number of diagnostic channels when the GDB sources and sheath boundary conditions (BCs) are selected carefully, especially the heat-flux BCs which can drastically alter the temperature. The radial profile of the fluctuation levels is qualitatively similar and quantitatively comparable on the low-field side, although statistics such as moments of the probability density function and the high-frequency spectrum show greater differences. This comparison indicates areas for future improvement in both simulations, such as sheath BCs, as well as improvements in GDB like particle conservation and spatially varying thermal conductivity, in order to achieve better fluid-gyrokinetic agreement and increase fidelity when simulating experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Investigating shear flow through continuum gyrokinetic simulations of limiter biasing in the Texas Helimak

Previous limiter-biasing experiments on the Texas Helimak, a simple magnetized torus, have been inconclusive on the effect of flow shear on turbulence levels. To investigate this, the first gyrokinetic simulations of limiter biasing in the Helimak using the plasma physics code Gkeyll have been carried out, and results are presented here. For the scenarios considered, turbulence is mostly driven by the interchange instability, which depends on gradients of equilibrium density profiles. An analysis of both experimental and simulation data demonstrates that shear rates are mostly less than than local linear growth rates, and not all requirements for shear stabilization are met. Rather, the mostly vertical shear flow has an important effect on bulk transport and experimental equilibrium density profiles, and changes to the gradients correspond to changes in turbulence levels.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Measurement of temperature of a dusty plasma from configuration

A new method called `Configurational Temperature' is introduced in the context of dusty plasma, where the temperature of the dust particles, submerged in the plasma, can be measured directly from the positional information of the individual dust particles and the interaction potential between the dust grains. This method does not require the velocity information of individual particles which is a key parameter to measure the dust temperature in the conventional method. The technique is initially tested using two dimensional OpenMP parallel Molecular Dynamics and Monte-Carlo simulation and then compared with the temperature evaluating from the experimental data. The experiments have been carried out in Dusty plasma experimental (DPEx) device where a two dimensional stationary plasma crystal of melamine formaldehyde particles is formed in the cathode sheath of a DC glow discharge argon plasma. The dust kinetic temperature is calculated using standard PIV technique at different pressures. The simulation results matches well with the experimental data at relatively higher pressures where the dust particles arranged into crystalline state or in a strongly coupled fluid state. An extended simulation results for three dimensional case is also presented which can be employed for the temperature measurement of three dimensional dust crystal in laboratory devices.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Conservative Discontinuous Galerkin Schemes for Nonlinear Fokker-Planck Collision Operators

We present a novel discontinuous Galerkin algorithm for the solution of a class of Fokker-Planck collision operators. These operators arise in many fields of physics, and our particular application is for kinetic plasma simulations. In particular, we focus on an operator often known as the `Lenard-Bernstein,' or `Dougherty,' operator. Several novel algorithmic innovations are reported. The concept of weak-equality is introduced and used to define weak-operators to compute primitive moments needed in the updates. Weak-equality is also used to determine a reconstruction procedure that allows an efficient and accurate discretization of the diffusion term. We show that when two integration by parts are used to construct the discrete weak-form, and finite velocity-space extents are accounted for, a scheme that conserves density, momentum and energy exactly is obtained. One novel feature is that the requirements of momentum and energy conservation lead to unique formulas to compute primitive moments. Careful definition of discretized moments also ensure that energy is conserved in the piecewise linear case, even though the $v^2$ term is not included in the basis-set used in the discretization. A series of benchmark problems are presented and show that the scheme conserves momentum and energy to machine precision. Empirical evidence also indicates that entropy is a non-decreasing function. The collision terms are combined with the Vlasov equation to study collisional Landau damping and plasma heating via magnetic pumping. We conclude with an outline of future work, in particular with some indications of how the algorithms presented here can be extended to use the Rosenbluth potentials to compute the drag and diffusion coefficients.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Exact and Locally Implicit Source Term Solvers for Multifluid-Maxwell Systems

Recently, a family of models that couple multifluid systems to the full Maxwell equations draw a lot of attention in laboratory, space, and astrophysical plasma modeling. These models are more complete descriptions of the plasma than reduced models like magnetohydrodynamic (MHD) since they naturally retain non-ideal effects like electron inertia, Hall term, pressure anisotropy/nongyrotropy, etc. One obstacle to broader application of these model is that an explicit treatment of their source terms leads to the need to resolve rapid kinetic processes like plasma oscillation and electron cyclotron motion, even when they are not important. In this paper, we suggest two ways to address this issue. First, we derive the analytic forms solutions to the source update equations, which can be implemented as a practical, but less generic solver. We then develop a time-centered, locally implicit algorithm to update the source terms, allowing stepping over the fast kinetic time-scales. For a plasma with $S$ species, the locally implict algorithm involves inverting a local $3S+3$ matrix only, thus is very efficient. The performance can be further elevated by using the direct update formulas to skip null calculations. Benchmarks illustrated the exact energy-conservation of the locally implicit solver, as well as its efficiency and robustness for both small-scale, idealized problems and large-scale, complex systems. The locally implicit algorithm can be also easily extended to include other local sources, like collisions and ionization, which are difficult to solve analytically.

0 citations0 reviewsNo code linkedOpen access