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
5586 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
                                    boundary_conditions = boundary_condition_periodic)

35


36
###############################################################################

37
# ODE solvers, callbacks etc.

38


39
tspan = (0.0, 2.0)

40
ode = semidiscretize(semi, tspan)

41


42
summary_callback = SummaryCallback()

43


44
analysis_interval = 1000

45
analysis_callback = AnalysisCallback(semi, interval = analysis_interval,

46
                                     extra_analysis_integrals = (energy_total,))

47


48
alive_callback = AliveCallback(analysis_interval = analysis_interval)

49


50
save_solution = SaveSolutionCallback(interval = 1000,

51
                                     save_initial_solution = true,

52
                                     save_final_solution = true,

53
                                     solution_variables = cons2prim)

54


55
callbacks = CallbackSet(summary_callback,

56
                        analysis_callback,

57
                        save_solution,

58
                        alive_callback)

59


60
###############################################################################

61
# run the simulation

62


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

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

65


66