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_timeintegration.jl
2055 views
1
using Trixi

2


3
###############################################################################

4
# semidiscretization of the linear advection equation

5


6
advection_velocity = (0.2, -0.7)

7
equations = LinearScalarAdvectionEquation2D(advection_velocity)

8


9
initial_condition = initial_condition_gauss

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

11


12
coordinates_min = (-5.0, -5.0)

13
coordinates_max = (5.0, 5.0)

14
mesh = TreeMesh(coordinates_min, coordinates_max,

15
                initial_refinement_level = 4,

16
                n_cells_max = 30_000)

17


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

19


20
###############################################################################

21
# ODE solvers, callbacks etc.

22


23
tspan = (0.0, 1.0)

24
ode = semidiscretize(semi, tspan)

25


26
summary_callback = SummaryCallback()

27


28
analysis_interval = 100

29
analysis_callback = AnalysisCallback(semi, interval = analysis_interval,

30
                                     extra_analysis_integrals = (entropy,))

31


32
alive_callback = AliveCallback(analysis_interval = analysis_interval)

33


34
save_solution = SaveSolutionCallback(interval = 100,

35
                                     save_initial_solution = true,

36
                                     save_final_solution = true,

37
                                     solution_variables = cons2cons)

38


39
amr_controller = ControllerThreeLevel(semi, IndicatorMax(semi, variable = first),

40
                                      base_level = 4,

41
                                      med_level = 5, med_threshold = 0.1,

42
                                      max_level = 6, max_threshold = 0.6)

43
amr_callback = AMRCallback(semi, amr_controller,

44
                           interval = 5,

45
                           adapt_initial_condition = true,

46
                           adapt_initial_condition_only_refine = true)

47


48
stepsize_callback = StepsizeCallback(cfl = 1.6)

49


50
callbacks = CallbackSet(summary_callback,

51
                        analysis_callback, alive_callback,

52
                        save_solution,

53
                        amr_callback, stepsize_callback)

54


55
###############################################################################

56
# run the simulation

57


58
# sol = solve(ode, CarpenterKennedy2N54(williamson_condition=false);

59
ode_algorithm = Trixi.CarpenterKennedy2N54()

60
sol = Trixi.solve(ode, ode_algorithm;

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

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

63


64