stardis-solver

commit a0c93d9ef2a914a72239be53aa812645b2795897
parent ec7f4632e459f3787a961bc0ca789f0874422e34
Author: Vincent Forest <vincent.forest@meso-star.com>
Date:   Fri, 19 Jan 2024 14:56:07 +0100

First test of WoS diffusion

Using the WoS algorithm on some pre-existing tests calculate the probe
temperature in a simple stationary configuration with dirichlet boundary
conditions. This mainly serves to check that the current implementation
does not present any significant numerical problems (in 2D and 3D).

Diffstat:
M
src/test_sdis_solve_probe.c
 | 
3
+++
M
src/test_sdis_solve_probe2.c
 | 
2
+-
M
src/test_sdis_solve_probe2_2d.c
 | 
1
+

3 files changed, 5 insertions(+), 1 deletion(-)
diff --git a/src/test_sdis_solve_probe.c b/src/test_sdis_solve_probe.c
@@ -397,6 +397,9 @@ main(int argc, char** argv)
   solve_args.picard_order = 0;
   BA(sdis_solve_probe(scn, &solve_args, &estimator));
   solve_args.picard_order = 1;
+  solve_args.diff_algo = SDIS_DIFFUSION_NONE;
+  BA(sdis_solve_probe(scn, &solve_args, &estimator));
+  solve_args.diff_algo = SDIS_DIFFUSION_DELTA_SPHERE;
   OK(sdis_solve_probe(scn, &solve_args, &estimator));
 
   BA(sdis_estimator_get_type(estimator, NULL));
diff --git a/src/test_sdis_solve_probe2.c b/src/test_sdis_solve_probe2.c
@@ -255,7 +255,7 @@ main(int argc, char** argv)
   solve_args.position[2] = 0.5;
   solve_args.time_range[0] = INF;
   solve_args.time_range[1] = INF;
-
+  solve_args.diff_algo = SDIS_DIFFUSION_WOS;
   OK(sdis_solve_probe(scn, &solve_args, &estimator));
 
   ref = 350 * solve_args.position[2] + (1-solve_args.position[2]) * 300;
diff --git a/src/test_sdis_solve_probe2_2d.c b/src/test_sdis_solve_probe2_2d.c
@@ -246,6 +246,7 @@ main(int argc, char** argv)
   solve_args.position[1] = 0.5;
   solve_args.time_range[0] = INF;
   solve_args.time_range[1] = INF;
+  solve_args.diff_algo = SDIS_DIFFUSION_WOS;
 
   OK(sdis_solve_probe(scn, &solve_args, &estimator));