Researcher profile

L. Chacon

L. Chacon contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

physics.plasm-ph physics.space-ph physics.comp-ph astro-ph.SR

Trust snapshot

Quick read

Trust 17 - Baseline

4works

0followers

4topics

4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Open graph Browse works

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

Building this map preview

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

preprint2021arXiv

A conservative implicit-PIC scheme for the hybrid kinetic-ion fluid-electron plasma model on curvilinear meshes

The hybrid kinetic-ion fluid-electron plasma model is widely used to study challenging multi-scale problems in space and laboratory plasma physics. Here, a novel conservative scheme for this model employing implicit particle-in-cell techniques is extended to arbitrary coordinate systems via curvilinear maps from logical to physical space. The scheme features a fully non-linear electromagnetic formulation with a multi-rate time advance - including sub-cycling and orbit-averaging for the kinetic ions. By careful choice of compatible particle-based kinetic-ion and mesh-based fluid-electron discretizations in curvilinear coordinates, as well as particle-mesh interpolations and implicit midpoint time advance, the scheme is proven to conserve total energy for arbitrary curvilinear meshes. In the electrostatic limit, the method is also proven to conserve total momentum for arbitrary curvilinear meshes. Although momentum is not conserved for arbitrary curvilinear meshes in the electromagnetic case, it is for an important subset of Cartesian tensor-packed meshes. The scheme and its novel conservation properties are demonstrated for several challenging numerical problems using different curvilinear meshes, including a merging flux-rope simulation for a space weather application, and a helical $m=1$ mode simulation for magnetic fusion energy application.

preprint2019arXiv

A cancellation problem in hybrid particle-in-cell schemes due to finite particle size

The quasi-neutral hybrid particle-in-cell algorithm with kinetic ions and fluid electrons is a popular model to study multi-scale problems in laboratory, space, and astrophysical plasmas. Here, it is shown that the treatment of ions as finite-size particles and electrons as a grid-based fluid can cause significant numerical wave dispersion errors in the magnetohydrodynamic limit ($kd_i \ll 1$, where $d_i$ is the ion skin-depth). Practical requirements on the mesh spacing $Δx/d_i$ are suggested to bound these errors from above.

preprint2015arXiv

Fluid vs. kinetic magnetic reconnection with strong guide-fields

The fast rates of magnetic reconnection found in both nature and experiments are important to understand theoretically. Recently, it was demonstrated that two-fluid magnetic reconnection remains fast in the strong guide field regime, regardless of the presence of fast-dispersive waves. This conclusion is in agreement with recent results from kinetic simulations, and is in contradiction to the findings in an earlier two-fluid study, where it was suggested that fast-dispersive waves are necessary for fast reconnection. In this paper, we give a more detailed derivation of the analytic model presented in a recent letter, and present additional simulation results to support the conclusions that the magnetic reconnection rate in this regime is independent of both collisional dissipation and system-size. In particular, we present a detailed comparison between fluid and kinetic simulations, finding good agreement in both the reconnection rate and overall length of the current layer. Finally, we revisit the earlier two-fluid study, which arrived at different conclusions, and suggest an alternative interpretation for the numerical results presented therein.

preprint2009arXiv

Current sheet bifurcation and collapse in electron magnetohydrodynamics

Inertial effects in nonlinear magnetic reconnection are studied within the context of 2D electron magnetohydrodynamics (EMHD) with resistive and viscous dissipation. Families of nonlinear solutions for relevant current sheet parameters are predicted and confirmed numerically in all regimes of interest. Electron inertia becomes important for current sheet thicknesses $δ$ below the inertial length $d_{e}$. In this case, in the absence of electron viscosity, the sheet thickness experiences a nonlinear collapse. Viscosity regularizes solutions at small scales. Transition from resistive to viscous regimes shows a nontrivial dependence on resistivity and viscosity, featuring a hysteresis bifurcation. In all accessible regimes, the nonlinear reconnection rate is found to be explicitly independent of the electron inertia and dissipation coefficients.