Write meta data to vtk files

src/util.jl

@@ -299,3 +299,7 @@ macro autoinfiltrate(condition=true)
                 lnn,
                 esc(condition))
 end
+
+function type2string(type)
+    return string(nameof(typeof(type)))
+end

src/visualization/write2vtk.jl

@@ -41,6 +41,9 @@ function trixi2vtk(v, u, t, system; output_directory="out", prefix="", iter=noth
         # Store particle index
         vtk["index"] = eachparticle(system)
 
+        # Write some meta data
+        vtk["time"] = t
+
         # Extract custom quantities for this system
         for (key, func) in custom_quantities
             value = func(v, u, t, system)
@@ -75,6 +78,16 @@ function write2vtk!(vtk, v, u, t, system::WeaklyCompressibleSPHSystem)
                       for particle in eachparticle(system)]
     vtk["pressure"] = system.pressure
 
+    # write meta data
+    vtk["solver"] = "WCSPH"
+    vtk["correction_method"] = type2string(system.correction)
+    vtk["acceleration"] = system.acceleration
+    vtk["viscosity"] = type2string(system.viscosity)
+    vtk["smoothing_kernel"] = type2string(system.smoothing_kernel)
+    vtk["smoothing_length"] = system.smoothing_length
+    vtk["density_calculator"] = type2string(system.density_calculator)
+    vtk["state_equation"] = type2string(system.state_equation)
+
     return vtk
 end
 
@@ -97,6 +110,15 @@ function write2vtk!(vtk, v, u, t, model, system)
 end
 
 function write2vtk!(vtk, v, u, t, model::BoundaryModelDummyParticles, system)
+    @unpack boundary_model = system
+
+    # write meta data
+    vtk["boundary_model"] = "BoundaryModelDummyParticles"
+    vtk["smoothing_kernel"] = type2string(boundary_model.smoothing_kernel)
+    vtk["smoothing_length"] = boundary_model.smoothing_length
+    vtk["density_calculator"] = type2string(boundary_model.density_calculator)
+    vtk["state_equation"] = type2string(boundary_model.state_equation)
+
     write2vtk!(vtk, v, u, t, model, model.viscosity, system)
 end
 
@@ -105,6 +127,8 @@ function write2vtk!(vtk, v, u, t, model::BoundaryModelDummyParticles, viscosity,
                                    for particle in eachparticle(system)]
     vtk["pressure"] = model.pressure
 
+    # write meta data
+    vtk["viscosity_model"] = type2string(viscosity)
     return vtk
 end
 
@@ -115,5 +139,8 @@ function write2vtk!(vtk, v, u, t, model::BoundaryModelDummyParticles,
     vtk["pressure"] = model.pressure
     vtk["wall_velocity"] = view(model.cache.wall_velocity, 1:ndims(system), :)
 
+    # write meta data
+    vtk["viscosity_model"] = "ViscosityAdami"
+
     return vtk
 end