Skip to content

Latest commit

 

History

History

figure_kelvin_helmholtz

README.md

Kelvin-Helmholtz instability

The files in this folder can be used to reproduce the numerical experiment simulating a Kelvin-Helmholtz instability with Trixi.jl. All results were obtained with Julia v1.6.2 on a computer running Kubuntu 20.04.

To reproduce the numerical experiments and the figure shown in the paper, proceed as follows.

  • Start Julia in this folder and activate the project environment therein (e.g., julia --project=. in this folder). Pass the command line option --check-bounds=no to increase the runtime performance. You can also use multithreading by passing the command line option --threads=XXX to Julia.
  • Execute the following code in Julia. 
    julia> using Trixi

julia> trixi_include("elixir_euler_kelvin_helmholtz_instability_amr.jl", tspan=(0.0, 2.0))
[output of the simulation]

julia> trixi_include("elixir_euler_kelvin_helmholtz_instability_amr.jl", tspan=(0.0, 3.0))
[output of the simulation]

If you want to modify some plots, you can adapt the plotting commands in the elixir or work in the Julia REPL as desired.