FiniteVolumeFluid

A developing library of finite volume method for compressible inviscous fluid. DO NOT use it in working environment.

Installation

] add https://github.com/HetaoZ/FiniteVolumeFluid.jl.git

Examples

1d SOD problem in 3d space. The results saved into VTK(or VTR/VTI) files can be viewed with Paraview.

using FiniteVolumeFluid

const gamma = 1.4

function initial_condition(x, t)
    if t == 0.0
        if x[1] < 0.5
            rho, u, p = 1, [0,0,0], 1e5
        else
            rho, u, p = 10, [0,0,0], 1e6
        end
        e = p/(gamma - 1)/rho
        return rho, u, e, p
    end
end

f = Fluid(dim = 3, point1 = (0,0,0), point2 = (10,1,1), ncells = (40, 4, 4), nbound = 2);

# set parameters
set_parameters(f, "gamma"=>1.4, "CFL"=>0.5)

set_initial_condition(f, initial_condition)

set_boundaries(f, (FreeBoundary, FreeBoundary), (ReflBoundary, ReflBoundary), (ReflBoundary, ReflBoundary))

set_scheme(f, Muscl())
set_scheme(f, Ausm())
set_scheme(f, RungeKutta(3))

# show(f)

tspan = (0, 1)
nframe = 0
t = 0
N = 10000

save_mesh(f, "../../out/coord")
save_to_vtk(f, ("rho","u","e","p"), (:rho,:u,:e,:p), "../../out/fluid_"*string(N))


for frame = 1:10
    if t > tspan[2]
        break
    end
    
    dt = time_step!(f)
    t += dt

    advance!(f, dt, t)

    save_to_vtk(f, ("rho","u","e","p"), (:rho,:u,:e,:p), "../../out/fluid_"*string(N + frame))

    println((frame, t, dt))
end