Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Write meta data to vtk files #199

Merged
merged 24 commits into from
Oct 16, 2023
Merged

Write meta data to vtk files #199

Changes from 1 commit
Commits
Show all changes
24 commits
Select commit Hold shift + click to select a range
47d9a4b
write some meta data
svchb Jul 15, 2023
722e782
Merge branch 'main' into save_meta_data
svchb Jul 16, 2023
91e6b42
add type2string function additional data
svchb Jul 16, 2023
c508cbf
Merge branch 'main' into save_meta_data
svchb Jul 21, 2023
4e59cb8
write the pvd file
svchb Jul 21, 2023
f90f017
Merge remote-tracking branch 'origin/main' into main
svchb Sep 5, 2023
372d67c
Merge branch 'trixi-framework:main' into main
svchb Sep 7, 2023
f564478
Merge branch 'main' of github.com:svchb/TrixiParticles.jl
svchb Sep 20, 2023
9c9e182
Merge branch 'trixi-framework:main' into main
svchb Sep 25, 2023
6b2e62a
Merge branch 'main'
svchb Sep 25, 2023
4a5e0fd
format
svchb Sep 25, 2023
0e95c9d
merge error
svchb Sep 25, 2023
cf07e81
add additional output for viscosity models
svchb Sep 25, 2023
c40fe4e
Merge branch 'main' into save_meta_data
svchb Sep 27, 2023
1951185
Merge branch 'main' into save_meta_data
svchb Oct 5, 2023
fd65721
make meta data optional
svchb Oct 5, 2023
58cd558
fix
svchb Oct 5, 2023
084512b
format
svchb Oct 5, 2023
f80a834
Merge branch 'main' into save_meta_data
svchb Oct 5, 2023
0d6a63b
Ignore markdown file changes
svchb Oct 6, 2023
2861be6
[skip ci] review
svchb Oct 13, 2023
8d99d1a
[skip ci] unpack all models
svchb Oct 14, 2023
51aae5a
[skip ci] add additional meta data
svchb Oct 16, 2023
5ed4caa
Merge branch 'main' into save_meta_data
efaulhaber Oct 16, 2023
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
22 changes: 10 additions & 12 deletions src/visualization/write2vtk.jl
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -72,8 +72,6 @@ function trixi2vtk(v, u, t, system, periodic_box; output_directory="out", prefix

# Store particle index
vtk["index"] = eachparticle(system)

# Write some meta data
vtk["time"] = t

# Extract custom quantities for this system
Expand All @@ -84,7 +82,7 @@ function trixi2vtk(v, u, t, system, periodic_box; output_directory="out", prefix
end
end

# add to collection
# Add to collection
pvd[t] = vtk
svchb marked this conversation as resolved.
Show resolved Hide resolved
end
vtk_save(pvd)
Expand Down Expand Up @@ -150,14 +148,14 @@ function write2vtk!(vtk, v, u, t, system::FluidSystem; write_meta_data=true)

if write_meta_data
vtk["acceleration"] = system.acceleration
vtk["viscosity"] = type2string(system.viscosity)
vtk["viscosity"] = repr(system.viscosity)
vtk["smoothing_kernel"] = type2string(system.smoothing_kernel)
vtk["smoothing_length"] = system.smoothing_length
vtk["density_calculator"] = type2string(system.density_calculator)

if system isa WeaklyCompressibleSPHSystem
vtk["correction_method"] = type2string(system.correction)
vtk["state_equation"] = type2string(system.state_equation)
vtk["correction_method"] = repr(system.correction)
vtk["state_equation"] = repr(system.state_equation)
vtk["solver"] = "WCSPH"
else
vtk["solver"] = "EDAC"
svchb marked this conversation as resolved.
Show resolved Hide resolved
Expand All @@ -170,14 +168,14 @@ function write2vtk!(vtk, v, u, t, system::FluidSystem; write_meta_data=true)
end

function write2vtk!(vtk, viscosity::ViscosityAdami)
vtk["nu"] = viscosity.nu
vtk["epsilon"] = viscosity.epsilon
vtk["viscosity_nu"] = viscosity.nu
vtk["viscosity_epsilon"] = viscosity.epsilon
end

function write2vtk!(vtk, viscosity::ArtificialViscosityMonaghan)
vtk["alpha"] = viscosity.alpha
vtk["beta"] = viscosity.beta
vtk["epsilon"] = viscosity.epsilon
vtk["viscosity_alpha"] = viscosity.alpha
vtk["viscosity_beta"] = viscosity.beta
vtk["viscosity_epsilon"] = viscosity.epsilon
end
svchb marked this conversation as resolved.
Show resolved Hide resolved

function write2vtk!(vtk, v, u, t, system::TotalLagrangianSPHSystem; write_meta_data=true)
svchb marked this conversation as resolved.
Show resolved Hide resolved
Expand Down Expand Up @@ -205,7 +203,7 @@ function write2vtk!(vtk, v, u, t, model::BoundaryModelDummyParticles, system;
vtk["smoothing_kernel"] = type2string(system.boundary_model.smoothing_kernel)
vtk["smoothing_length"] = system.boundary_model.smoothing_length
vtk["density_calculator"] = type2string(system.boundary_model.density_calculator)
vtk["state_equation"] = type2string(system.boundary_model.state_equation)
vtk["state_equation"] = repr(system.boundary_model.state_equation)
end

write2vtk!(vtk, v, u, t, model, model.viscosity, system,
Expand Down