Book a Demo!
CoCalc Logo Icon
StoreFeaturesDocsShareSupportNewsAboutPoliciesSign UpSign In
trixi-framework
GitHub Repository: trixi-framework/Trixi.jl
 Path: blob/main/examples/tree_3d_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, 0.5)

8
equations = LinearScalarAdvectionEquation3D(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, -1.0)

14
coordinates_max = (1.0, 1.0, 1.0)

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

16
                       coordinates_max = (1.0, 1.0, 1.0)),

17
                      (type = "box", coordinates_min = (0.0, -0.5, -0.5),

18
                       coordinates_max = (0.5, 0.5, 0.5)))

19
mesh = TreeMesh(coordinates_min, coordinates_max,

20
                initial_refinement_level = 2,

21
                refinement_patches = refinement_patches,

22
                n_cells_max = 10_000)

23


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

25


26
###############################################################################

27
# ODE solvers, callbacks etc.

28


29
tspan = (0.0, 5.0)

30
ode = semidiscretize(semi, tspan)

31


32
summary_callback = SummaryCallback()

33


34
analysis_interval = 100

35
analysis_callback = AnalysisCallback(semi, interval = analysis_interval,

36
                                     extra_analysis_integrals = (entropy,))

37


38
alive_callback = AliveCallback(analysis_interval = analysis_interval)

39


40
save_solution = SaveSolutionCallback(interval = 100,

41
                                     save_initial_solution = true,

42
                                     save_final_solution = true,

43
                                     solution_variables = cons2prim)

44


45
stepsize_callback = StepsizeCallback(cfl = 1.2)

46


47
callbacks = CallbackSet(summary_callback,

48
                        analysis_callback,

49
                        alive_callback,

50
                        save_solution,

51
                        stepsize_callback)

52


53
###############################################################################

54
# run the simulation

55


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

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

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

59


60