summaryrefslogtreecommitdiff
path: root/examples/one_particle_sedimentation_2d/common.hpp
diff options
context:
space:
mode:
authorClaudius "keldu" Holeksa <mail@keldu.de>2026-08-22 22:24:12 +0200
committerClaudius "keldu" Holeksa <mail@keldu.de>2026-08-22 22:24:25 +0200
commite6168bb3f656ebab9b63be697c79819b08fc1ae0 (patch)
tree0724726c6b16ba27fa7972fffb64501b22e89eda /examples/one_particle_sedimentation_2d/common.hpp
parent47c1f28230ec96ca446a86ce3af0294837359d94 (diff)
downloadlibs-lbm-e6168bb3f656ebab9b63be697c79819b08fc1ae0.tar.gz
Moved folde
Diffstat (limited to 'examples/one_particle_sedimentation_2d/common.hpp')
-rw-r--r--examples/one_particle_sedimentation_2d/common.hpp71
1 files changed, 71 insertions, 0 deletions
diff --git a/examples/one_particle_sedimentation_2d/common.hpp b/examples/one_particle_sedimentation_2d/common.hpp
new file mode 100644
index 0000000..37e0402
--- /dev/null
+++ b/examples/one_particle_sedimentation_2d/common.hpp
@@ -0,0 +1,71 @@
+#pragma once
+
+#include <kel/lbm/sycl/lbm.hpp>
+#include <kel/lbm/lbm.hpp>
+#include <kel/lbm/particle.hpp>
+#include <kel/lbm/particle/particle.hpp>
+
+#include <forstio/io/io.hpp>
+#include <forstio/remote/filesystem/easy.hpp>
+#include <forstio/codec/json/json.hpp>
+#include <forstio/codec/simple.hpp>
+
+namespace kel {
+namespace lbm {
+
+/**
+ * length 2.2m
+ * height 0.41m
+ * viscosity 10^-3
+ * particle diameter 0.1m
+ * position at (0.16m,0.2m)
+ */
+
+constexpr uint64_t dim_y = 512ul;
+constexpr uint64_t dim_x = 256ul;
+
+constexpr uint64_t particle_amount = 1ul;
+
+namespace sch {
+using namespace saw::schema;
+
+using InfoChunk = Chunk<UInt8, 0u, dim_x, dim_y>;
+
+template<typename T, typename Desc>
+using DfChunk = Chunk<FixedArray<T,Desc::Q>, 1u, dim_x, dim_y>;
+
+template<typename T, typename Desc>
+using ScalarChunk = Chunk<Scalar<T>, 0u, dim_x, dim_y>;
+
+template<typename T, typename Desc>
+using VectorChunk = Chunk<Vector<T,Desc::D>, 0u, dim_x, dim_y>;
+
+template<typename T, typename Desc>
+using ChunkStruct = Struct<
+ Member<InfoChunk, "info">,
+ Member<DfChunk<T,Desc>, "dfs">,
+ Member<DfChunk<T,Desc>, "dfs_old">,
+ Member<VectorChunk<T,Desc>, "particle_N">,
+ Member<ScalarChunk<T,Desc>, "particle_D">
+>;
+
+template<typename T, typename Desc>
+using VelChunk = Chunk<Vector<T,Desc::D>, 0u, dim_x, dim_y>;
+
+template<typename T>
+using RhoChunk = Chunk<Scalar<T>, 0u, dim_x, dim_y>;
+
+template<typename T, typename Desc>
+using MacroStruct = Struct<
+ Member<VelChunk<T,Desc>, "momentum">,
+ Member<RhoChunk<T>, "density">,
+ Member<ScalarChunk<T,Desc>, "porosity">,
+ Member<VelChunk<T,Desc>, "force">
+>;
+
+template<typename T, typename Desc>
+using ParticleSpheroidGroup = ParticleGroup<T,Desc::D,1u,coll::Spheroid<T>>;
+}
+
+}
+}