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
5586 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, periodicity = true)

23


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

25
                                    boundary_conditions = boundary_condition_periodic)

26


27
###############################################################################

28
# ODE solvers, callbacks etc.

29


30
tspan = (0.0, 5.0)

31
ode = semidiscretize(semi, tspan)

32


33
summary_callback = SummaryCallback()

34


35
analysis_interval = 100

36
analysis_callback = AnalysisCallback(semi, interval = analysis_interval,

37
                                     extra_analysis_integrals = (entropy,))

38


39
alive_callback = AliveCallback(analysis_interval = analysis_interval)

40


41
save_solution = SaveSolutionCallback(interval = 100,

42
                                     save_initial_solution = true,

43
                                     save_final_solution = true,

44
                                     solution_variables = cons2prim)

45


46
stepsize_callback = StepsizeCallback(cfl = 1.2)

47


48
callbacks = CallbackSet(summary_callback,

49
                        analysis_callback,

50
                        alive_callback,

51
                        save_solution,

52
                        stepsize_callback)

53


54
###############################################################################

55
# run the simulation

56


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

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

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

60


61