diff options
| author | stainer_t <thomas.stainer@oecd-nea.org> | 2025-09-08 13:48:49 +0200 |
|---|---|---|
| committer | stainer_t <thomas.stainer@oecd-nea.org> | 2025-09-08 13:48:49 +0200 |
| commit | 7dfcc480ba1e19bd3232349fc733caef94034292 (patch) | |
| tree | 03ee104eb8846d5cc1a981d267687a729185d3f3 /Donjon/data/Test_SPH_void.x2m | |
Initial commit from Polytechnique Montreal
Diffstat (limited to 'Donjon/data/Test_SPH_void.x2m')
| -rw-r--r-- | Donjon/data/Test_SPH_void.x2m | 66 |
1 files changed, 66 insertions, 0 deletions
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" . |