[TRNSYS-users] controlling the passage to the next time steps
Michaël Kummert
michael.kummert at polymtl.ca
Mon Jul 12 05:51:52 PDT 2021
Hello,
The co-simulator with ESP-r available in TRNSYS 17 and 18 implements
something that is similar to what you are trying to do, if I understand
correctly.
Once TRNSYS has converged for a time step, Type 130 checks that the
iterations between the 2 programs have converged, and if that is not the
case, Type 130 “forces” TRNSYS to continue iterating at the current time
step.
Perhaps you can take a look at Type 130 and how it interacts with the Kernel
to see whether you can use that for what you are trying to achieve.
The part I am not sure to understand is that in the co-simulation case, it
makes sense to force additional iterations because the inputs coming from
ESP-r will be different so the TRNSYS results after new iterations will be
different. In your case, I am not sure what would cause TRNSYS to find
different results, but I am probably missing something.
In any case, I hope this helps.
Michaël Kummert
Polytechnique Montréal
From: Pierre Jaboyedoff via TRNSYS-users
<trnsys-users at lists.onebuilding.org>
Sent: July 9, 2021 08:35
To: trnsys-users at lists.onebuilding.org
Cc: pierre.jaboyedoff at effinart.ch
Subject: [TRNSYS-users] controlling the passage to the next time steps
In trnsys 13.1, i had (long time ago ..) modified the kernel of TRNSYS (see
below).
The reason was to be able to check if the control routine developed was
converging to the correct physical case (e.g. solar, then auxiliary etc
).
If not then i was getting TRNSYS to restart the iteration on the same time
step until the correct control case was obtained.
Has anybody done something similar for TRNSYS 17.2 ?
I have started looking for this, it seems to be located between, exec.for,
trnsys.for, etc
Thanks in advance,
Pierre Jaboyedoff
Head of BEEP PMTU Switzerland (www.beepindia.org <http://www.beepindia.org>
)
Effin’Art
Associate
Chemin de Pré-Fleuri 6
Entrance at the back, 3rd floor on the left
CH-1006 Lausanne
Tel : 00 41 21 616 11 00
Mobile : 00 41 79 668 50 88
pierre.jaboyedoff at effinart.ch <mailto:pierre.jaboyedoff at effinart.ch>
-------------------------------------------------------------------------
-------------------------------------------------------------------------
in blockdat
C MODIF PJ BEGIN
COMMON /CHKSTEP/ ireiter, icas, icaschk
COMMON /sunsp/ tsunsp, tset, tmax, tload, tsto, text
C MODIF PJ END
-------------------------------------------------------------------------
-------------------------------------------------------------------------
-------------------------------------------------------------------------
in main
-------------------------------------------------------------------------
c modif pj begin
COMMON /CHKSTEP/ IREITER , icas, icaschk
COMMON /sunsp/ tsunsp, tset, tmax, tload, tsto, text
c modif pj end
-------------------------------------------------------------------------
C ********************** MAIN SIMULATION LOOP **************************
100 IF (DONE) GO TO 999
CALL CLOCK(TIME,DONE,2)
C modif pj begin
c added label 10001 used for reiter of the time step
10001 IF (NDERIV.EQ.0) THEN
c modif pj end
-------------------------------------------------------------------------
370 CONTINUE
ENDIF
C SWITCH PREDICTED AND CORRECTED VALUES
I = PRED
PRED = CORR
CORR = I
INTG = 2
IF (IRECUR.GT.0) INTG = 3
IF ((IRECUR.GT.0.OR.NCALLS.GT.0).AND.(ICT.EQ.0)) GO TO 200
C ************************************************************
C MODIF PJ BEGIN (LABELED COMMON ADDED IN BLOCKDAT_M.FOR, )
c control_iter retuen a variable reiter (in common labelled CHKSTEP)
C REITER =0 GO ON TO NEXT TIME STEP
C =1 REITER ON THE SAME TIME STEP
call control_iter(time)
if (ireiter .gt. 0) then
ireiter=ireiter+1
icasmain=icas
icaschkm=icaschk
timereit=max(0.,time)
c WRITE(*,*) "reiter ", timereit, icasmain, icaschkm
tsunspm=tsunsp
tsetm=tset
tloadm=tload
tmaxm=tmax
stom=tsto
textm=text
c write (*,*) " tsunsp=", tsunspm, " tset=",tsetm,"
text=",textm
c write (*,*) " tload=", tloadm, " tmax=",tmaxm," tsto=",tsto
c write (*,*) icasmain
c write (*,*) icaschkm
c CALL PRINT(TIME,DONE)
goto 10001
endif
C MODIF PJ END
-------------------------------------------------------------------------
-------------------------------------------------------------------------
-------------------------------------------------------------------------
-------------------------------------------------------------------------
in exec
-------------------------------------------------------------------------
c modif pj begin
COMMON /CHKSTEP/ IREITER, icas, icaschk
c modif pj end
-------------------------------------------------------------------------
110 INFO(7,J) = -1
ELSE
DO 111 J = 1, NUNITS
CALL(J) = CALL1(J)
c modif pj begin
c added info(7,j) =0 for the case where one wants to reiter on the same
c time step
if (ireiter .le. 0) then
INFO(7,J) = -1
else
CALL(J) = .TRUE.
INFO(7,J) = 0
endif
111 continue
c modif pj end
-------------------------------------------------------------------------
-------------------------------------------------------------------------
-------------------------------------------------------------------------
-------------------------------------------------------------------------
In my own control routine
c modif pj begin
COMMON /CHKSTEP/ IREITER, icas, icaschk
COMMON /sunsp/ tsunsp, tset, tmax, tload, tsto, text
SAVE /sunsp/
c modif pj end
Pierre Jaboyedoff
Head of BEEP PMTU Switzerland ( <http://www.beepindia.org>
www.beepindia.org)
Effin’Art
Associate
Chemin de Pré-Fleuri 6
Entrance at the back, 3rd floor on the left
CH-1006 Lausanne
Tel : 00 41 21 616 11 00
Mobile : 00 41 79 668 50 88
pierre.jaboyedoff at effinart.ch <mailto:pierre.jaboyedoff at effinart.ch>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.onebuilding.org/pipermail/trnsys-users-onebuilding.org/attachments/20210712/8c910a56/attachment-0001.html>
More information about the TRNSYS-users
mailing list