From 7dfcc480ba1e19bd3232349fc733caef94034292 Mon Sep 17 00:00:00 2001
From: stainer_t <thomas.stainer@oecd-nea.org>
Date: Mon, 8 Sep 2025 13:48:49 +0200
Subject: Initial commit from Polytechnique Montreal

---
 Donjon/data/Test_SPH_void.x2m | 66 +++++++++++++++++++++++++++++++++++++++++++
 1 file changed, 66 insertions(+)
 create mode 100644 Donjon/data/Test_SPH_void.x2m

(limited to 'Donjon/data/Test_SPH_void.x2m')

diff --git a/Donjon/data/Test_SPH_void.x2m b/Donjon/data/Test_SPH_void.x2m
new file mode 100644
index 0000000..f1a6b60
--- /dev/null
+++ b/Donjon/data/Test_SPH_void.x2m
@@ -0,0 +1,66 @@
+*----
+*  Name          : TEST Test_SPH_void.x2m
+*  Author        : Alain Hebert (2018)
+*
+*  17 X 17 PWR VOIDED ASSEMBLY WITH NEWTONIAN SPH EQUIVALENCE
+*  Apply SPH correction on albedo factors
+*----
+*  Define STRUCTURES and MODULES used
+*----
+LINKED_LIST
+  SPHGEOM MTRACK MACRO MACRO2 SYS OPTIM OUT CALC ;
+MODULE
+  GEO: MAC: TRIVAT: TRIVAA: FLUD: DELETE: UTL: GREP: END: ;
+DOUBLE Fobj ;
+SEQ_ASCII _OUTV :: FILE './_OUTV' ;
+PROCEDURE SPH_Optim ;
+PROCEDURE assertS ;
+INTEGER IterEmax := 3000 ;
+
+*----
+* TRANSPORT-DIFFUSION NEWTONIAN SPH EQUIVALENCE
+*----
+OUT := _OUTV ;
+SPHGEOM := GEO: :: CAR2D 9 9
+  X- DIAG X+ ALBE 1 Y- SYME Y+ DIAG
+  MIX   1  2  3  4  5  6  7  8  9
+          10 11 12 13 14 15 16 17
+             18 19 20 21 22 23 24
+                25 26 27 28 29 30
+                   31 32 33 34 35
+                      36 37 38 39
+                         40 41 42
+                            43 44
+                               45
+  MESHX 0.0 1.26471996 2.52943993 3.79415989 5.05887985
+        6.32359982 7.58831978 8.85303974 10.1177597 11.4324799
+  MESHY 0.0 1.26471996 2.52943993 3.79415989 5.05887985
+        6.32359982 7.58831978 8.85303974 10.1177597 11.4324799
+;
+MTRACK := TRIVAT: SPHGEOM :: DUAL (*IELEM=*) 1 (*ICOL=*) 1
+                  SPN 3 SCAT 1 ;
+MACRO := OUT :: STEP UP 'REF-CASE0001' STEP UP 'MACROLIB' ;
+*----
+* Set physical albedos
+*----
+MACRO := MAC: MACRO :: ALBP 1 0.0 0.0 ;
+*----
+* Perform Newtonian SPH equivalence
+*----
+MACRO2 OPTIM := SPH_Optim MACRO SPHGEOM MTRACK
+   :: 0.5 2.5 1.0E-4 <<IterEmax>> ;
+ECHO "control variables at convergence" ;
+UTL: OPTIM :: IMPR 'VAR-VALUE' * ;
+
+GREP: OPTIM :: GETVAL 'FOBJ-CST-VAL' 1 >>Fobj<< ;
+ECHO "optimal objective function=" Fobj ;
+*----
+* Verification calculation
+*----
+SYS := TRIVAA: MACRO2 MTRACK ;
+CALC := FLUD: SYS MTRACK :: EDIT 2 ADI 5 ACCE 5 3 ;
+assertS CALC :: K-EFFECTIVE 1 0.4555094 ;
+
+ECHO "test Test_SPH_void completed" ;
+END: ; 
+QUIT "LIST" .
-- 
cgit v1.2.3