Version 16
*
*Solares TWW - System (Hi Flo)
*Schritt 4: KW-Zumischung im TW-Kreis, fiktiver WT im Koll.kreis
*FHTW Berlin / Prof. Dr.-Ing. Friedrich Sick
*
ASSIGN  "DE-Potsdam-103790.tm2"	10
ASSIGN  "WWLast200.dat"		11
ASSIGN  "***.s01"		12
ASSIGN  "***.s02"		22
ASSIGN  "***.o01"		13
ASSIGN  "***.o02"		14
ASSIGN  "***.o03"		15
*
EQUATIONS 17
*
JAN = 0
FEB = 744
MAR = 1416
APR = 2160
MAY = 2880
JUN = 3624
JUL = 4344
AUG = 5088
SEP = 5832
OCT = 6552
NOV = 7296
DEC = 8016
SIL = 8760
*
START  = JAN
STOP   = SIL
STEP   = 0.25
TOL    = 0.001
*
SIMULATION  START   STOP   STEP
TOLERANCES  TOL     TOL
*
EQUATIONS 26
*
* Ausrichtung und Umgebung
RHOG   = 0.2
SLO    = 35
AZI    = 0
*
* Kollektorspezifikationen:
* Fläche eines Moduls (m2), Anzahl Kollektoren, spez. Durchströmung (kg/m2-h)
* eta0, U1 (kJ/h-m2-K) und U2 (kJ/h-m2-K2)
AKOLL  = 2.01
NKOLL  = 3
FLOWSP = 40
ETA0   = 0.837
UWERT1 = 12.312
UWERT2 = 0.05508
AREA   = NKOLL*AKOLL
FLOMAX = FLOWSP*AREA
*
* Spezifikationen zu Regler und Pumpe im Solarkreis
* cp-Koll.Fluid (kJ/kg-K), Delta-T-EIN, Hysterese,
* Max. zulässige Speichertemperatur
CPSOL  = 3.6
DTH    = 15
*DTL    = 5
* sichere untere Schaltdifferenz gem. Duffie/Beckman (siehe TRNSYS Doku):
DTL    = MAX(1,DTH*AREA*UWERT1/FLOMAX/CPSOL)
TAUS   = 95
PUMPOW = 180
*
* Speicher-Spezifikationen
* Vol (m3), U (kJ/h-m2-K), Höhe (m, neg. bei gleich hohen Schichten!)
* Solltemperatur, Hysterese (K), Leistung (kJ/h)
VOLSP  = 0.45
USP    = 2.13
HEIGHT = -1.752
TSPSET = 60
HYSTSP = 5
AUXPOW = 54000
*
* Temperaturspezifikationen (Celsius)
* Kaltwasser, Solltemperatur Warmwasser, Raumtemperaturen im Haus, im Keller
TKW    = 15
TSET   = 50
TRAUM  = 20
TKELL  = 15
*
*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*
*
UNIT 2 TYPE 109	 Weather data and radiation processor combined
PARAMETERS 4
* Mode LUN DiffMode TrackMode
  2    10  4        1
INPUTS 3
* rho_g  slope  azi
  0,0    0,0    0,0
  RHOG   SLO    AZI 
*Outputs
*  1: Tamb,  12: Itot,   13: Ibeam,  14: Idiff_sky,
*            18: IKoll,  19: IbKoll, 20: IdKoll_sky, 21: IdKoll_grd
*
*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*
*
EQUATIONS 8
IBEAM  = [2,13]/3.6
IDIFF  = [2,14]/3.6
IGLOB  = [2,12]/3.6
IKOLL  = [2,18]/3.6
IBKOLL = [2,19]/3.6
IDKOLL = ([2,20]+[2,21])/3.6
TAMB   = [2,1]
*
GAM    = [6,1]
*
*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*
*
UNIT 8 TYPE 9 TRNSYS DATA READER: DHW LOAD
PARAMETERS 10
*MODE SKIP Number of Values  time interval
 5    0    1                 0.25
*Interpol	Mult.factor	Addition term   Average
 -1		1		0		0
*LUN	FRMT
 11     -1
* Output: 	1
*		Liter WW
*
EQUATIONS 1
MWW  = [8,1]/0.25
*
*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*
*
UNIT 3 TYPE 3 PUMP
PARAMETERS 4
*MMAX   CP    PMAX   FPAR
FLOMAX  CPSOL PUMPOW 0.0
INPUTS 3
33,1  33,2   GAM
20.   0      0
*
*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*
*
UNIT 4 TYPE 1 SOLAR COLLECTOR
PARAMETERS 13
*NS A     CP    EffMode Gtest  eta0   k1     k2     optM alpha NG etaR  KL
 1  AREA  CPSOL 1       60     ETA0   UWERT1 UWERT2 4    0.95  1  1.526 0.016
INPUTS 9
* Ti   mdot Tamb  Ikoll ItotHor IdifHor albedo theta beta
  3,1  3,2  TAMB  2,18  2,12    2,14    0,0    2,22  0,0
  20   200  10    0     0       0       RHOG   0     SLO
* Output: 	1	2	3
*               To      mdot    Quse
*
*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*
*
UNIT 33 TYPE 91 Heat Exchanger
PARAMETERS 3
* eff    Cph   CpC
  1      CPSOL 4.2
INPUTS 4
* Thi  mdot Tci  mdot
  4,1  4,2  7,1  4,2
  20   0    TKW  0
*
*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*
*
UNIT 7 TYPE 4 STRATIFIED FLUID STORAGE TANK
PARAMETERS 22
* MODE Vol   Cp  rho  Ut  H1     AuxM L1  LT1 Tset1  DT1    Qdot1  
  3    VOLSP 4.2 1000 USP HEIGHT 1    2   1   TSPSET HYSTSP AUXPOW
* L2  LT2 Tset2 DT2 Qdot2 UAf Tf  Tboil HotIn KWin
  4   3   15    3   0     0   20  99    4     5
INPUTS 5
* TH    MH    TL   ML   Tenv
  33,3  33,4  30,1 30,2 0,0
  20    0     TKW  0    TKELL
DERIVATIVES 5
15 15 15 15 15
*
*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*
*
UNIT 6 TYPE 2 CONTROLLER
PARAMETERS 2
* NSTK TMAX
  5    TAUS
INPUTS 6
* TH  TL   TIN  gam  DTH  DTL
  4,1 7,15 7,3  6,1  0,0  0,0
  20  20   20   0    DTH  DTL
*
*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*
*
UNIT 30 TYPE 11 Temperature controlled Flow diverter
PARAMETERS 2
*mode NSTK
 4    5
INPUTS 4
* Ti   mdot Thot Tset
  39,1 39,2 7,3  0,0
  20   0    20   TSET
*
*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*
*
UNIT 31 TYPE 11 T piece
PARAMETERS 1
*mode
 1
INPUTS 4
* T1   mdot1 T2   mdot2
  30,3 30,4  7,3  7,4
  20   0     20   50
*
*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*
*
EQUATIONS 1
MENTN = MWW
*
UNIT 32 TYPE 91 Heat Exchanger
PARAMETERS 3
* eff  Cph  CpC
  0.9  4.2  4.2
INPUTS 4
* Thi  mdot Tci  mdot
  31,1 MENTN  0,0  MWW
  20   0    TKW  0
*
*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*
*
EQUATION 2
DINS = 0.016
LMIN = 4*100*STEP/3.141592/1000/DINS/DINS
*
* Ohne Rohr bekommt der Flow Diverter zu Zeiten ohne Massenstrom
* die Temperatur am Eingang des HX. Obwohl dann ohne Bedeutung,
* werden lästige Warnungen geschrieben, ggf. Abbruch.
* Daher Einbau eines Rohres mit Volumen > Massenstrom*dt
* ohne thermische Verluste.
UNIT 39 TYPE 31 Rohr WT-FlowDiverter
PARAMETERS 6
*di    l     U   rho  cp   T_initial
DINS   LMIN  .1  1000 4.2  TKELL
INPUTS 3
*Ti    mdot  Tenv
 32,1  32,2  0,0
 20.0  0.0   TKELL*
*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*
*
EQUATIONS 16
TKVL  = [4,1]
TKRL  = [7,1]
TSPO  = [7,3]
TSPM  = [7,13]
TSPU  = [7,15]
TWT   = [32,3]
TWW   = [32,3]
ITOT  = IKOLL*AREA
QKOLL = [4,3]/3.6
QSPIN = [7,11]/3.6
QAUX  = [7,8]/3.6
QSPV  = [7,5]/3.6
QSPLD = [7,6]/3.6
QHXLD = [32,5]/3.6
LOADI = MWW*(TWW-TKW)*4.2/3.6
LOADS = MWW*(TSET-TKW)*4.2/3.6
*
*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*
*
UNIT 65 TYPE 28 SIMULATION SUMMARY RESULTS 1
PARAMETERS 39
-1 START STOP 12 2
-11 -1 1000 2 -4    -12 -1 1000 2 -4    -13 -1 1000 2 -4   -14 -1 1000 2 -4
-15 -1 1000 2 -4    -16 -1 1000 2 -4    -17 -2 2 -4
INPUTS 7
IBEAM IDIFF IGLOB IBKOLL IDKOLL IKOLL TAMB
0 0 0 0 0 0 0
LABELS 7
IBEAM IDIFF IGLOB IBKOLL IDKOLL IKOLL TAMB
*
*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*
*
UNIT 66 TYPE 28 SIMULATION SUMMARY RESULTS 2
PARAMETERS 40
-1 START STOP 22 2
-11 -1 1000 2 -4    -12 -1 1000 2 -4    -13 -1 1000 2 -4   -14 -1 1000 2 -4
-15 -1 1000 2 -4    -16 -1 1000 2 -4    -17 -1 1000 2 -4
INPUTS 7
ITOT QKOLL QSPIN QAUX QSPV LOADI LOADS
0 0 0 0 0 0 0
LABELS 7
ITOT QKOLL QSPIN QAUX QSPV LOADI LOADS
*
*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*
*
UNIT 68 TYPE 65 ONLINE GRAPHICS
PARAMETERS 12
*NLEFT NRIGHT YMIN1 YMAX1 YMIN2 YMAX2 NPIC XGRID on/off LU  units ofd
 6     1      0     1000  -20   80    12   4     0      -13 0     0
INPUTS 7
IBEAM IDIFF IGLOB IKOLL IBKOLL IDKOLL TAMB
IBEAM IDIFF IGLOB IKOLL IBKOLL IDKOLL TAMB
LABELS 3
"Einstrahlung W/m²"
"Temperatur °C"
"Einstrahlung"
*
*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*
*
EQUATIONS 1
GAM300  = 300*GAM
*
UNIT 69 TYPE 65 ONLINE GRAPHICS
PARAMETERS 12
*NLEFT NRIGHT YMIN1 YMAX1 YMIN2 YMAX2 NPIC XGRID on/off LU  units ofd
 7     2      0     100   0     2000  12   4     1      -14 0     0
INPUTS 9
TKVL TSPO TSPM TSPU 32,1 TWW 39,1 IKOLL GAM300
TKVL TSPO TSPM TSPU T321 TWW T391 IKOLL GAMPUMP
LABELS 3
"Temperatur °C"
"Massenstrom kg/h"
"Temperaturen"
*
*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*
*
UNIT 70 TYPE 65 ONLINE GRAPHICS
PARAMETERS 12
*NLEFT NRIGHT YMIN1 YMAX1 YMIN2 YMAX2 NPIC XGRID on/off LU  units ofd
 8     0      0     100   0     2000  12   4     1      -15 0     0
INPUTS 8
39,2 32,4 30,2 30,4 31,2 7,4 33,2 MWW
M392 M324 M302 M304 M312 M74 M332 MWW
LABELS 3
"Massenstrom kg/h"
"-"
"Massenströme"
*
*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*
*
END