stardis-solver

commit 4388b4a81332cb506d8c71d2c337ed14d761733f
parent 5bbd028ae62a70ee28a55af9d6b979ab27e0db7c
Author: Vincent Forest <vincent.forest@meso-star.com>
Date:   Tue, 22 May 2018 16:14:30 +0200

Add comments to the volumic_power<3|4>_2d tests

Diffstat:
M
src/test_sdis_volumic_power3_2d.c
 | 
25
+++++++++++++++++++++++++
M
src/test_sdis_volumic_power4_2d.c
 | 
27
+++++++++++++++++++++++++--

2 files changed, 50 insertions(+), 2 deletions(-)
diff --git a/src/test_sdis_volumic_power3_2d.c b/src/test_sdis_volumic_power3_2d.c
@@ -26,6 +26,31 @@
 #define H1 5.0
 #define H2 10.0
 
+/*
+ * The 2D scene is composed of 3 stacked solid slabs whose middle slab has a
+ * volumic power. The +/-X sides of the slabs are stretched far away to
+ * simulate a 1D case. The upper and lower bounds of the "sandwich" has a
+ * convective exchange with the surrounding fluid whose temperature is known.
+ *
+ *           _\  T1
+ *          / /
+ *          \__/
+ * ... -----H1------ ...
+ *       LAMBDA1
+ *
+ * ... ------------- ...
+ *       LAMBDA, Pw
+ * ... ------------- ...
+ *
+ *       LAMBDA2
+ *
+ *
+ * ... -----H2------ ...
+ *            _\  T2
+ *           / /
+ *           \__/
+ */
+
 static const double vertices[8/*#vertices*/*2/*#coords per vertex*/] = {
  -100000.5, 0.0,
  -100000.5, 1.4,
diff --git a/src/test_sdis_volumic_power4_2d.c b/src/test_sdis_volumic_power4_2d.c
@@ -25,6 +25,29 @@
 #define Delta (1.0/20.0)
 #define Nrealisations 10000
 
+/*
+ * The 2D scene is a solid slabs stretched along the X dimension to simulate a
+ * 1D case. The slab has a volumic power and has a convective exchange with the
+ * surrounding fluid whose temperature is fixed to Tfluid.
+ *
+ *
+ *           _\  TFluid
+ *          / /
+ *          \__/
+ *
+ * ... -----Hboundary----- ...
+ *
+ *        Lambda, Power
+ *
+ * ... -----Hboundary----- ...
+ *
+ *           _\  TFluid
+ *          / /
+ *          \__/
+ *
+ */
+
+
 static const double vertices[4/*#vertices*/*2/*#coords per vertex*/] = {
  -10000.5,-0.5,
  -10000.5, 0.5,
@@ -301,8 +324,8 @@ main(int argc, char** argv)
   } else {
     Tinterf = Power*L / (2*Hboundary) + Tfluid;
   }
-  Tref = 
-    Tinterf 
+  Tref =
+    Tinterf
   + Power / (2*Lambda) * ((L*L)/4.0 - pos[1]*pos[1]);
 
   CHK(sdis_solve_probe(scn, Nrealisations, pos, INF, 1.f, -1, 0, &estimator) == RES_OK);