Book a Demo!
CoCalc Logo Icon
StoreFeaturesDocsShareSupportNewsAboutPoliciesSign UpSign In
trixi-framework
GitHub Repository: trixi-framework/Trixi.jl
 Path: blob/main/examples/structured_1d_dgsem/elixir_advection_nonperiodic.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
boundary_conditions = BoundaryConditionDirichlet(initial_condition)

12


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

14


15
coordinates_min = (0.0,)

16
coordinates_max = (5.0,)

17
mesh = StructuredMesh((16,), coordinates_min, coordinates_max, periodicity = false)

18


19
semi = SemidiscretizationHyperbolic(mesh, equations,

20
                                    initial_condition,

21
                                    solver; boundary_conditions = boundary_conditions)

22


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

24
# ODE solvers, callbacks etc.

25


26
tspan = (0.0, 5.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_restart = SaveRestartCallback(interval = 100,

38
                                   save_final_restart = true)

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.6)

46


47
callbacks = CallbackSet(summary_callback,

48
                        analysis_callback, alive_callback,

49
                        save_restart, save_solution,

50
                        stepsize_callback);

51


52
###############################################################################

53
# run the simulation

54


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

56
            dt = stepsize_callback(ode), # solve needs some value here but it will be overwritten by the stepsize_callback

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

58


59