scalar_wave_1d.lua

simulation = {
    mesh = {
        index_extents = {71},
        domain_bounds = {2}
    },
    domain_boundaries = {
        xmin = "dirichlet"
    },
    shapes = {
        {
            type = "yz_rect",
            psi = 0.1,
            normal = 1,
            boundary_condition = "dirichlet"
        }
    },
    scheme = {
        order = 1,
        type = "E2",
        alpha = {-1.47956280234494, 0.261900367793859, -0.145072532538541, -0.224665713988644}
    },
    system = {
        type = "scalar wave",
        center = {-1},
        radius = 0
        -- diffusivity = 1 / 30
    },
    integrator = {
        type = "rk4"
    },
    step_controller = {
        max_time = 10,
        --max_step = 2,
        cfl = {
            hyperbolic = 0.8,
            parabolic = 0.2
        }
    },
    io = {
        --write_every_step = 1
        write_every_time = 0.8
    }
}