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

2
using Trixi

3


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

5
# semidiscretization of the linear advection equation

6


7
advection_velocity = 1.0

8
equations = LinearScalarAdvectionEquation1D(advection_velocity)

9


10
initial_condition = initial_condition_gauss

11


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

13


14
coordinates_min = (-5.0,)

15
coordinates_max = (5.0,)

16
mesh = TreeMesh(coordinates_min, coordinates_max,

17
                initial_refinement_level = 4,

18
                n_cells_max = 30_000, periodicity = true)

19


20
semi = SemidiscretizationHyperbolic(mesh, equations, initial_condition, solver;

21
                                    boundary_conditions = boundary_condition_periodic)

22


23
###############################################################################

24
# ODE solvers, callbacks etc.

25


26
tspan = (0.0, 10.0)

27
ode = semidiscretize(semi, tspan)

28


29
summary_callback = SummaryCallback()

30


31
analysis_interval = 100

32
analysis_callback = AnalysisCallback(semi, interval = analysis_interval,

33
                                     extra_analysis_integrals = (entropy,))

34


35
alive_callback = AliveCallback(analysis_interval = analysis_interval)

36


37
save_solution = SaveSolutionCallback(interval = 100,

38
                                     save_initial_solution = true,

39
                                     save_final_solution = true,

40
                                     solution_variables = cons2prim)

41


42
indicator_function(u, x, t) = u[1]

43


44
amr_controller = ControllerThreeLevel(semi,

45
                                      IndicatorNodalFunction(indicator_function, semi),

46
                                      base_level = 4,

47
                                      med_level = 5, med_threshold = 0.1,

48
                                      max_level = 6, max_threshold = 0.6)

49
amr_callback = AMRCallback(semi, amr_controller,

50
                           interval = 5,

51
                           adapt_initial_condition = true,

52
                           adapt_initial_condition_only_refine = true)

53


54
stepsize_callback = StepsizeCallback(cfl = 1.6)

55


56
callbacks = CallbackSet(summary_callback,

57
                        analysis_callback, alive_callback,

58
                        save_solution,

59
                        amr_callback, stepsize_callback)

60


61
###############################################################################

62
# run the simulation

63


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

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

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

67


68