Book a Demo!
CoCalc Logo Icon
StoreFeaturesDocsShareSupportNewsAboutPoliciesSign UpSign In
trixi-framework
GitHub Repository: trixi-framework/Trixi.jl
 Path: blob/main/examples/tree_2d_fdsbp/elixir_euler_convergence.jl
2055 views
1
# !!! warning "Experimental implementation (upwind SBP)"

2
#     This is an experimental feature and may change in future releases.

3


4
using OrdinaryDiffEqSSPRK

5
using Trixi

6


7
###############################################################################

8
# semidiscretization of the compressible Euler equations

9


10
equations = CompressibleEulerEquations2D(1.4)

11


12
initial_condition = initial_condition_convergence_test

13
source_terms = source_terms_convergence_test

14


15
D_upw = upwind_operators(SummationByPartsOperators.Mattsson2017,

16
                         derivative_order = 1,

17
                         accuracy_order = 4,

18
                         xmin = -1.0, xmax = 1.0,

19
                         N = 16)

20
flux_splitting = splitting_steger_warming

21
solver = FDSBP(D_upw,

22
               surface_integral = SurfaceIntegralUpwind(flux_splitting),

23
               volume_integral = VolumeIntegralUpwind(flux_splitting))

24


25
coordinates_min = (-1.0, -1.0)

26
coordinates_max = (1.0, 1.0)

27
mesh = TreeMesh(coordinates_min, coordinates_max,

28
                initial_refinement_level = 3,

29
                n_cells_max = 30_000,

30
                periodicity = true)

31


32
semi = SemidiscretizationHyperbolic(mesh, equations, initial_condition, solver,

33
                                    source_terms = source_terms_convergence_test)

34


35
###############################################################################

36
# ODE solvers, callbacks etc.

37


38
tspan = (0.0, 2.0)

39
ode = semidiscretize(semi, tspan)

40


41
summary_callback = SummaryCallback()

42


43
analysis_interval = 1000

44
analysis_callback = AnalysisCallback(semi, interval = analysis_interval,

45
                                     extra_analysis_integrals = (energy_total,))

46


47
alive_callback = AliveCallback(analysis_interval = analysis_interval)

48


49
save_solution = SaveSolutionCallback(interval = 1000,

50
                                     save_initial_solution = true,

51
                                     save_final_solution = true,

52
                                     solution_variables = cons2prim)

53


54
callbacks = CallbackSet(summary_callback,

55
                        analysis_callback,

56
                        save_solution,

57
                        alive_callback)

58


59
###############################################################################

60
# run the simulation

61


62
sol = solve(ode, SSPRK43(); abstol = 1.0e-9, reltol = 1.0e-9,

63
            ode_default_options()..., callback = callbacks)

64


65