Papers › On the Kelvin-Helmholtz instability with smooth initial conditions -- Linear theory...

On the Kelvin-Helmholtz instability with smooth initial conditions -- Linear theory and simulations

4 Feb 2019arXiv:1902.01403links table onlyarchive 2025-07-28

Thomas Berlok, Christoph Pfrommer

The archive published only this paper's code-link row. Authors, date and abstract are from arXiv's metadata (CC0), read from the Kaggle arXiv metadata snapshot of 2026-09-12 where its title matched the archive's; the title is the archive's.

The Kelvin-Helmholtz instability (KHI) is a standard test of hydrodynamic and magnetohydrodynamic (MHD) simulation codes and finds many applications in astrophysics. The classic linear theory considers a discontinuity in density and velocity at the interface of two fluids. However, for numerical simulations of the KHI such initial conditions do not yield converged results even at the linear stage of the instability. Instead, smooth profiles of velocity and density are required for convergence. This renders analytical theory to be only approximately valid and hinders quantitative comparisons between the classical theory and simulations. In this paper we derive a linear theory for the KHI with smooth profiles and illustrate code testing with the MHD code Athena. We provide the linear solution for the KHI with smooth initial conditions in three different limits: inviscid hydrodynamics, ideal MHD and Braginskii-MHD. These linear solutions are obtained numerically with the framework Psecas (Pseudo-Spectral Eigenvalue Calculator with an Automated Solver), which generates and solves numerical eigenvalue problems using an equation-parser and pseudo-spectral methods. The Athena simulations are carried out on a periodic, Cartesian domain which is useful for code testing purposes. Using Psecas and analytic theory, we outline the differences between this artificial numerical setup and the KHI on an infinite Cartesian domain and the KHI in cylindrical geometry. We discuss several astrophysical applications, such as cold flows in galaxy formation and cold fronts in galaxy cluster mergers. Psecas, and the linear solutions used for code testing, are publicly available and can be downloaded from the web.

PaperPDFCodeCode Syntology ran

In Syntology Open this paper in Syntology's Atlas, the map of the papers in Syntology's graph and their citations.

For agents, Syntology's MCP tool lists every function and class Syntology harvested from this paper and whether it ran (how to connect): get_harvested_code_for_paper(arxiv_id="1902.01403")

Code

Syntology Ran 0 of 6 code samples harvested from 1 repository linked to this paper; 6 have no recorded run.

By repository: official repository: 6 samples from 1 repository, 0 ran. The run record, sample by sample. “Ran” means executed on a synthesized input, not that the code is correct or reproduces the paper.

tberlok/psecas officialmentioned on GitHubBSD-3-Clause report

Repository list and official/mentioned flags are the archive's, frozen 2025-07-28. Reachability, where shown, is from one Syntology probe window (2026-09-16 to 2026-09-18); repositories not probed show nothing. GitHub stars are not tracked.

Code Syntology ran Syntology

6 samples harvested; 0 ran; 0 honoured the contract we drafted; 6 have no recorded run. Read from Syntology's graph 2026-09-24; that is when this build read the record, not when the samples ran.

6unverified

Licence: 0 of the 6 samples are pointer only, meaning Syntology does not serve that copy's text. This page shows no code text for any sample; each one links to its file in the repository.

Harvested from tberlok/psecas. “Ran” means the sample executed on a synthesized input. It does not mean the output is correct, and nothing here reproduces the paper's results. “Honoured” and “violated” refer to a contract Syntology drafted from the code itself; “our draft was wrong” and “fixture could not drive it” are failures of Syntology's instrument, not of the code.

Each sample ends with its code_sha256, Syntology's identity for that exact code. An agent fetches the stored sample with Syntology's MCP tool get_code(code_sha256="…") (how to connect); click an identity to copy that call.

Repository labels, per sample. official repository: The archive marks this repository official for the paper. named in the paper: The archive records that the paper mentions this repository; it is not marked official. community (archive-listed): In the archive's code links for this paper, not marked official and not recorded as mentioned in the paper. found in paper text by Syntology: Syntology found this repository in the paper's own text; whether it is the authors' implementation is not asserted. community: Not in the archive's code links for this paper; a community repository Syntology harvested. Samples from a repository marked official are listed first. Licence labels name the repository's licence as recorded at harvest. “Pointer only” means Syntology does not serve that copy's text, for one of four reasons: no licence file was found; the licence was not identified; the licence is recorded as permissive but that copy's record is not marked cleared; or the licence is outside the permissive list Syntology serves text under (MIT, Apache-2.0, BSD and similar). Some licences outside that list permit redistribution, such as WTFPL, and GPL-3.0 under its conditions; they are simply not on the list. Hover a licence label for the reason. File links open the file on GitHub at the default branch, which may have changed since the harvest.

get_2D_cylindrical_map tberlok/psecas/psecas/plotting.py official repository unverified BSD-3-Clause (permissive) · 54164a91f74cce6d · report
get_2D_cylindrical_map_in_cylindrical_coords tberlok/psecas/psecas/plotting.py official repository unverified BSD-3-Clause (permissive) · c02ef653e9847056 · report
get_2Dmap tberlok/psecas/psecas/plotting.py official repository unverified BSD-3-Clause (permissive) · a0c64bbdfcd4c675 · report
golden_section tberlok/psecas/psecas/utils.py official repository unverified BSD-3-Clause (permissive) · a02023397920ae8a · report
load_system tberlok/psecas/psecas/serial_io.py official repository unverified BSD-3-Clause (permissive) · afecd0a258a990a5 · report
var_replace tberlok/psecas/psecas/string_methods.py official repository unverified BSD-3-Clause (permissive) · a46ca7fa7814294f · report

Results from the paper archive 2025-07-28

No leaderboard rows for this paper in the archive.

Report a problem or propose a change · a person checks every report against the paper or source before anything changes; decisions are listed on /corrections