blob: 28219c8d71bf3e292537ab5176e25217c3a5da34 (
plain)
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
|
*----
* Name : TEST Test_SPH.x2m
* Author : Alain Hebert (2018)
*
* 17 X 17 PWR ASSEMBLY WITH NEWTONIAN SPH EQUIVALENCE
*----
* Define STRUCTURES and MODULES used
*----
LINKED_LIST
SPHGEOM MTRACK MACRO MACRO2 SYS OPTIM OUT CALC ;
MODULE
TRIVAT: TRIVAA: FLUD: DELETE: UTL: GREP: END: ;
DOUBLE Fobj ;
SEQ_ASCII _OUT :: FILE './_OUT' ;
PROCEDURE SPH_Optim ;
PROCEDURE assertS ;
INTEGER IterEmax := 1000 ;
*----
* TRANSPORT-DIFFUSION NEWTONIAN SPH EQUIVALENCE
*----
OUT := _OUT ;
SPHGEOM := OUT :: STEP UP 'MACRO-GEOM' ;
MTRACK := TRIVAT: SPHGEOM :: DUAL (*IELEM=*) 1 (*ICOL=*) 1 ;
MACRO := OUT :: STEP UP 'REF-CASE0001' STEP UP 'MACROLIB' ;
MACRO2 OPTIM := SPH_Optim MACRO SPHGEOM MTRACK
:: 0.5 1.5 1.0E-3 <<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 4 ACCE 5 3 ;
assertS CALC :: K-EFFECTIVE 1 0.9554875 ;
ECHO "test Test_SPH completed" ;
END: ;
QUIT "LIST" .
|
