BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Kai Germaschewski

Kai Germaschewski appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
physics.plasm-phphysics.comp-phphysics.space-phastro-ph.SRMachine Learning

Catalog footprint

What is connected

8works
5topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

8 published item(s)

preprint2021arXiv

Inference and De-Noising of Non-Gaussian Particle Distribution Functions: A Generative Modeling Approach

The particle-in-cell numerical method of plasma physics balances a trade-off between computational cost and intrinsic noise. Inference on data produced by these simulations generally consists of binning the data to recover the particle distribution function, from which physical processes may be investigated. In addition to containing noise, the distribution function is temporally dynamic and can be non-gaussian and multi-modal, making the task of modeling it difficult. Here we demonstrate the use of normalizing flows to learn a smooth, tractable approximation to the noisy particle distribution function. We demonstrate that the resulting data driven likelihood conserves relevant physics and may be extended to encapsulate the temporal evolution of the distribution function.

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
preprint2016arXiv

Effects of electron temperature anisotropy on proton mirror instability evolution

Proton mirror modes are large amplitude nonpropagating structures frequently observed in the magnetosheath. It has been suggested that electron temperature anisotropy can enhance the proton mirror instability growth rate while leaving the proton cyclotron instability largely unaffected, therefore causing the proton mirror instability to dominate the proton cyclotron instability in Earth's magnetosheath. Here, we use particle-in-cell simulations to investigate the electron temperature anisotropy effects on proton mirror instability evolution. Contrary to the hypothesis, electron temperature anisotropy leads to excitement of the electron whistler instability. Our results show that the electron whistler instability grows much faster than the proton mirror instability and quickly consumes the electron free energy, so that there is no electron temperature anisotropy left to significantly impact the evolution of the proton mirror instability.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Effect of electron diamagnetic drifts on cylindrical double-tearing modes

Double-tearing modes (DTMs) have been proposed as a driver of `off-axis sawtooth' crashes in reverse magnetic shear tokamak configurations. Recently differential rotation provided by equilibrium sheared flows has been shown capable of decoupling the two DTM resonant layers, slowing the growth the instability. In this work we instead supply this differential rotation using an electron diamagnetic drift, which emerges in the presence of an equilibrium pressure gradient and finite Larmor radius physics. Diamagnetic drifts have the additional benefit of stabilizing reconnection local to the two tearing layers. Conducting linear and nonlinear simulations with the extended MHD code MRC-3d, we consider an m=2, n=1 cylindrical double-tearing mode. We show that asymmetries between the resonant layers and the emergence of an ideal MHD instability cause the DTM evolution to be highly dependent on the location of the pressure gradient. By locating a strong drift near the outer, dominant resonant surface are we able to saturate the mode and preserve the annular current ring, suggesting that the appearance of DTM activity in advanced tokamaks depends strongly on the details of the plasma pressure profile.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The island coalescence problem: scaling of reconnection in extended fluid models including higher-order moments

As modeling of collisionless magnetic reconnection in most space plasmas with realistic parameters is beyond the capability of today's simulations, due to the separation between global and kinetic length scales, it is important to establish scaling relations in model problems so as to extrapolate to realistic scales. Recently, large scale particle-in-cell (PIC) simulations of island coalescence have shown that the time averaged reconnection rate decreases with system size, while fluid systems at such large scales in the Hall regime have not been studied. Here we perform the complementary resistive MHD, Hall MHD and two fluid simulations using a ten-moment model with the same geometry. In contrast to the standard Harris sheet reconnection problem, Hall MHD is insufficient to capture the physics of the reconnection region. Additionally, motivated by the results of a recent set of hybrid simulations which show the importance of ion kinetics in this geometry, we evaluate the efficacy of the ten-moment model in reproducing such results.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The Plasma Simulation Code: A modern particle-in-cell code with load-balancing and GPU support

Recent increases in supercomputing power, driven by the multi-core revolution and accelerators such as the IBM Cell processor, graphics processing units (GPUs) and Intel's Many Integrated Core (MIC) technology have enabled kinetic simulations of plasmas at unprecedented resolutions, but changing HPC architectures also come with challenges for writing efficient numerical codes. This paper describes the Plasma Simulation Code (PSC), an explicit, electromagnetic particle-in-cell code with support for different order particle shape functions. We focus on two distinguishing feature of the code: patch-based load balancing using space-filling curves, and support for Nvidia GPUs, which achieves substantial speed-up of up to more than 6x on the Cray XK7 architecture compared to a CPU-only implementation.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Instability-Enhanced Friction in the Presheath of Two-Ion-Species Plasmas

The speed at which ions enter a sheath is a fundamental property of a plasma that also provides a useful boundary condition in modeling. A recent theory proposed that this can be significantly influenced by an instability-enhanced friction force arising from two-stream instabilities in the presheath when multiple ion species are present. Although experiments appeared to confirm this theory, recent particle simulations have brought it into question. We reconcile this controversy using direct numerical solutions of the dispersion relation, which show that there is a dependence on the electron-ion temperature ratio that was not considered previously. In addition, particle-in-cell simulations are used to show that ion-ion two-stream instabilities can arise near the sheath edge and generate an enhanced ion-ion friction force. Only by accounting for the instability-enhanced friction force can theory predict the simulated ion speeds at the sheath edge.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Perpendicular Ion Heating by Low-Frequency Alfven-Wave Turbulence in the Solar Wind

We consider ion heating by turbulent Alfven waves (AWs) and kinetic Alfven waves (KAWs) with perpendicular wavelengths comparable to the ion gyroradius and frequencies smaller than the ion cyclotron frequency. When the turbulence amplitude exceeds a certain threshold, an ion's orbit becomes chaotic. The ion then interacts stochastically with the time-varying electrostatic potential, and the ion's energy undergoes a random walk. Using phenomenological arguments, we derive an analytic expression for the rates at which different ion species are heated, which we test by simulating test particles interacting with a spectrum of randomly phased AWs and KAWs. We find that the stochastic heating rate depends sensitively on the quantity epsilon = dv/vperp, where vperp is the component of the ion velocity perpendicular to the background magnetic field B0, and dv (dB) is the rms amplitude of the velocity (magnetic-field) fluctuations at the gyroradius scale. In the case of thermal protons, when epsilon << eps1, where eps1 is a constant, a proton's magnetic moment is nearly conserved and stochastic heating is extremely weak. However, when epsilon > eps1, the proton heating rate exceeds the cascade power that would be present in strong balanced KAW turbulence with the same value of dv, and magnetic-moment conservation is violated. For the random-phase waves in our test-particle simulations, eps1 is approximately 0.2. For protons in low-beta plasmas, epsilon is approximately dB/B0 divided by the square root of beta, and epsilon can exceed eps1 even when dB/B0 << eps1. At comparable temperatures, alpha particles and minor ions have larger values of epsilon than protons and are heated more efficiently as a result. We discuss the implications of our results for ion heating in coronal holes and the solar wind.

0 citations0 reviewsNo code linkedOpen access