Book a Demo!
CoCalc Logo Icon
StoreFeaturesDocsShareSupportNewsAboutPoliciesSign UpSign In
trixi-framework
GitHub Repository: trixi-framework/Trixi.jl
 Path: blob/main/examples/tree_2d_dgsem/elixir_advection_mortar.jl
2055 views
1
using OrdinaryDiffEqLowStorageRK

2
using Trixi

3


4
###############################################################################

5
# semidiscretization of the linear advection equation

6


7
advection_velocity = (0.2, -0.7)

8
equations = LinearScalarAdvectionEquation2D(advection_velocity)

9


10
initial_condition = initial_condition_convergence_test

11
solver = DGSEM(polydeg = 3, surface_flux = flux_lax_friedrichs)

12


13
coordinates_min = (-1.0, -1.0)

14
coordinates_max = (1.0, 1.0)

15
refinement_patches = ((type = "box", coordinates_min = (0.0, -1.0),

16
                       coordinates_max = (1.0, 1.0)),)

17
mesh = TreeMesh(coordinates_min, coordinates_max,

18
                initial_refinement_level = 2,

19
                refinement_patches = refinement_patches,

20
                n_cells_max = 10_000)

21


22
semi = SemidiscretizationHyperbolic(mesh, equations, initial_condition, solver)

23


24
###############################################################################

25
# ODE solvers, callbacks etc.

26


27
tspan = (0.0, 1.0)

28
ode = semidiscretize(semi, tspan)

29


30
summary_callback = SummaryCallback()

31


32
analysis_interval = 100

33
analysis_callback = AnalysisCallback(semi, interval = analysis_interval,

34
                                     extra_analysis_integrals = (entropy,))

35


36
alive_callback = AliveCallback(analysis_interval = analysis_interval)

37


38
save_solution = SaveSolutionCallback(interval = 100,

39
                                     save_initial_solution = true,

40
                                     save_final_solution = true,

41
                                     solution_variables = cons2prim)

42


43
stepsize_callback = StepsizeCallback(cfl = 1.6)

44


45
callbacks = CallbackSet(summary_callback,

46
                        analysis_callback, alive_callback,

47
                        save_solution,

48
                        stepsize_callback)

49


50
###############################################################################

51
# run the simulation

52


53
sol = solve(ode, CarpenterKennedy2N54(williamson_condition = false);

54
            dt = 1.0, # solve needs some value here but it will be overwritten by the stepsize_callback

55
            ode_default_options()..., callback = callbacks);

56


57