stardis-solver

commit a58caadfdb60d2c0ca9efd3a5bbb59b6476e4d29
parent 7e0983198ec7ac6653b74d02ae2e04f46ca0349a
Author: Vincent Forest <vincent.forest@meso-star.com>
Date:   Wed, 20 Jun 2018 16:22:31 +0200

Fix how volumic bias on boundary is handled in 3D

Diffstat:
M
src/sdis_solve_Xd.h
 | 
8
+++++++-

1 file changed, 7 insertions(+), 1 deletion(-)
diff --git a/src/sdis_solve_Xd.h b/src/sdis_solve_Xd.h
@@ -993,9 +993,15 @@ XD(solid_temperature)
         /* Add the power corrective term */
         if(h < delta_solid) {
           const double sin_a = h / delta_solid;
+#if DIM==2
           /* tmp1 = sin(2a) / (PI - 2*a) */
           const double tmp1 = sin_a * sqrt(1 - sin_a*sin_a)/acos(sin_a);
-          tmp += -(power*delta_s_in_meter*delta_s_in_meter)/(2.0*DIM*lambda) * tmp1;
+          tmp += -(power*delta_s_in_meter*delta_s_in_meter)/(4.0*lambda) * tmp1;
+#else
+          const double tmp1 = (sin_a*sin_a*sin_a - sin_a)/ (1-sin_a);
+          tmp += (power*delta_s_in_meter*delta_s_in_meter)/(6*lambda) * tmp1;
+#endif
+
         } else if (h == delta_solid) {
           tmp += -(delta_s_in_meter*delta_s_in_meter*power)/(2.0*DIM*lambda);
         }