* fichier : pressuhx2.dgibi * * Containment pressurization (Phebus size) with heat losses. * 3D mesh of a 14.5 m3 cylindrical containment with a 10cm in depth vertical wall. * * Initial pressure and temperature are 1bar and 40oC with the wall at 60oC. * Steam is injected at a constant 50g/s mass flow rate and a temperature of 150oC. * * Transient is computed for 1000s. We check the condensation mass flow * rate evolution (that starts at about 20s). * * Compare to pressu.dgibi, hereafter heat losses are considered : the * containment is in contact with air or water. The water table elevation * is moving with time ; the convective heat transfer coefficient is * computed by the Mc Adams correlation (by use of a personnel procedur) ; * air and water characteristic lenght and temperature are known with * time for both. * * Auteurs : E. Studer, J.P. Magnaud Novembre 1999 * 'OPTI' 'DIME' 3 'ELEM' 'CU20' 'DENS' 1. 'TRAC' 'PSC' ; * COMPLET = FAUX ; GRAPH = faux ; 'SI' COMPLET ; nbit = 100 ; DT0 = 1. ; n1 = 1 ; n2 = 4 ; n3 = 4 ; nn = 2 ; 'SINO' ; nbit = 100 ; DT0 = 10. ; n1 = 1 ; n2 = 2 ; n3 = 4 ; nn = 1 ; 'FINS' ; epsi = 1.D-2 ; * * * Mesh * * * ri=cavity radius, h1=cavity height, sp=wall depth, * fg1=break to cavity radius ratio ri = 1.052 ; h1 = 4.163 ; sp = 0.10 ; fg1 = 0.25 ; fg2 = fg1 * (2.0 ** 0.5) / 2. ; * p0 = 0.000 0.000 0.000 ; p1 = (ri*fg1) 0.000 0.000 ; p2 = (ri*fg2) (ri*fg2) 0.000 ; p3 = 0.000 (ri*fg1) 0.000 ; p4 = ri 0.000 0.000 ; p5 = 0.000 ri 0.000 ; p6 = (ri+sp) 0.000 0.000 ; p7 = 0.000 (ri+sp) 0.000 ; * * Fluid and structure basement surfaces * (built by symetry) l1 = 'DROI' p0 p1 n1 ; l2 = 'DROI' p1 p2 n1 ; l3 = 'DROI' p2 p3 n1 ; l4 = 'DROI' p3 p0 n1 ; l5 = 'CERC' p4 p0 p5 (2*n1) ; l6 = 'CERC' p6 p0 p7 (2*n1) ; basf0 = 'DALL' l1 l2 l3 l4 'PLAN' ; basf1 = ('REGL' (l2 'ET' l3) l5 n2) ; * l44 = 'COTE' 2 basf1; ax4 = ('INVE' l4) 'ET' l44 ; l11 = 'COTE' 4 basf1; ax1 = l11 'ET' ('INVE' l1) ; basf = basf0 'ET' ('REGL' (l2 'ET' l3) l5 n2) ; 'ELIM' basf epsi ; basf = basf 'ET' ('SYME' basf 'DROI' p0 p3) ; ax11 = ('SYME' ax1 'DROI' p0 p3) 'ET' ('INVE' ax1) ; 'ELIM' basf epsi ; basf = basf 'ET' ('SYME' basf 'DROI' p0 p1) ; ax44 = ('INVE' ax4) 'ET' ('SYME' ax4 'DROI' p0 p4) ; 'ELIM' basf epsi ; basm = 'REGL' l5 l6 n3 ; basm = basm 'ET' ('SYME' basm 'DROI' p0 p3) ; 'ELIM' basm epsi ; basm = basm 'ET' ('SYME' basm 'DROI' p0 p1) ; 'ELIM' basm epsi ; l5 = l5 'ET' ('SYME' l5 'DROI' p0 p1) ; l5 = l5 'ET' ('SYME' l5 'DROI' p0 p3) ; l6 = l6 'ET' ('SYME' l6 'DROI' p0 p1) ; l6 = l6 'ET' ('SYME' l6 'DROI' p0 p3) ; * * Fluid and structure volumes * (built by translation) nz1 = ('ENTI' (h1 '/' (ri '/' 2.))) '*' nn ; v1 = 0. 0. h1 ; mt = basf 'VOLU' nz1 'TRAN' v1 ; wall = (basm 'VOLU' nz1 'TRAN' v1) 'COUL' 'ROUG' ; * plan1 = ax11 'TRAN' nz1 v1 ; plan4 = ax44 'TRAN' nz1 v1 ; pari = l5 'TRAN' nz1 v1 ; pare = l6 'TRAN' nz1 v1 ; pti = 'POIN' 1 l5 ; pte = 'POIN' 1 l6 ; lzi = pti d nz1 (pti 'PLUS' v1) ; lze = pte d nz1 (pte 'PLUS' v1) ; * 'ELIM' (mt 'ET' wall 'ET' plan1 'ET' plan4 'ET' pari 'ET' pare 'ET' lzi 'ET' lze) epsi ; * * Break at the basement if any pjg = 'POIN' basf 'PROC' (0. 0. 0.) ; breche = ('ELEM' basf 'APPUIE' 'LARGEMENT' pjg) 'COUL' 'VERT' ; * * * Data for execrxt.procedur * * rxt = 'TABLE' ; rxt . 'VERSION' = 'V0' ; rxt . 'vtf' = mt ; rxt . 'epsi' = epsi ; rxt . 'pi' = 0. 0. 0.5 ; * rxt . 'DISCR' = 'MACRO'; rxt . 'KPRE' = 'CENTRE'; rxt . 'DT0' = DT0 ; * rxt . 'MODTURB' = 'NUTURB' ; rxt . 'NUT' = 1.D-2 ; * rxt . 'THERMP' = VRAI ; rxt . 'vtp' = wall ; rxt . 'LAMBDA' = 15. ; rxt . 'ROCP' = 3.9E6 ; rxt . 'Tp0' = 60. ; rxt . 'ECHAN' = 10. ; * * rxt . 'VAPEUR' = VRAI ; rxt . 'TF0' = 40.0 ; rxt . 'PT0' = 1.0e5 ; rxt . 'Yvap0' = 0.0023 ; * rxt . 'Breches' = 'TABLE' ; rxt . 'Breches' . 'A' = 'TABLE' ; rxt . 'Breches' . 'A' . 'scenario' = 'TABLE' ; rxt . 'Breches' . 'A' . 'Maillage' = breche ; rxt . 'Breches' . 'A' . 'diru' = (0. 0. 1.) ; rxt . 'Breches' . 'A' . 'scenario' . 't' = 'PROG' 0.0 1000.0 ; rxt . 'Breches' . 'A' . 'scenario' . 'qeau' = 'PROG' 0.050 0.050 ; rxt . 'Breches' . 'A' . 'scenario' . 'qair' = 'PROG' 0.000 0.000 ; rxt . 'Breches' . 'A' . 'scenario' . 'tinj' = 'PROG' 150.0 150.0 ; * rxt . 'GRAPH' = GRAPH ; rxt . 'DETMAT' = VRAI ; rxt . 'FRPREC' = 5 ; rxt . 'RENU' = 'RIEN' ; * *=================================================================== * Use of a personnal procedur to impose external temperature and heat * transfer coefficient * rxt . 'ECHEXT' = vrai ; rxt . 'parext' = pare ; rxt . 'HEXT' = 0. ; rxt . 'TPEXT' = 0. ; rxt . 'PERSO' = vrai ; rxt . 'PRCPERSO' = 'MAPROCX' ; rxt . 'TABPERSO' = 'TABLE' ; * 'DEBPROC' MAPROCX TPS*'FLOTTANT' rxt*'TABLE' ; * * Heat transfer coeff given by Mc Adams correlation with h=hair if * z>z(t) and h=hwater otherwise (case KAS2 of the PROCHEXT procedur). zt = 'EVOL' 'MANU' ('PROG' 0. 1.e3) ('PROG' 2. 2.) ; Lair = 3. ; tair = 'EVOL' 'MANU' ('PROG' 0. 1.e3) ('PROG' 45. 45.) ; Lmer = 2. ; tmer = 'EVOL' 'MANU' ('PROG' 0. 1.e3) ('PROG' 25. 25.) ; PROCHEXT TPS RXT 'KAS2' ZT TAIR LAIR TMER LMER ; 'FINP' ; * *=================================================================== * * * Transient (with restart after 2 time steps) * * EXECRXT 2 rxt ; EXECRXT (nbit - 2) rxt ; * * * Tests * * 'SI' ('NON' COMPLET) ; ERR1 = 0 ; list rxt . 'TIC' . 'Qc' ; lqc = 'PROG' 0.0000 0.0000 0.0000 4.06187E-04 2.65273E-03 4.88913E-03 9.33217E-03 1.21696E-02 1.62809E-02 1.83569E-02 2.15760E-02 2.30394E-02 2.49146E-02 2.71021E-02 2.88218E-02 3.11268E-02 3.25368E-02 3.41987E-02 3.55088E-02 3.65889E-02 3.75529E-02 3.84019E-02 3.91586E-02 3.98414E-02 4.04639E-02 4.10347E-02 4.15599E-02 4.20449E-02 4.24941E-02 4.29111E-02 4.32985E-02 4.36591E-02 4.39950E-02 4.43082E-02 4.46005E-02 ; lqc = lqc 'ET' ('PROG' 4.48734E-02 4.51285E-02 4.53669E-02 4.55900E-02 4.57989E-02 4.59945E-02 4.61778E-02 4.63498E-02 4.65111E-02 4.66625E-02 4.68047E-02 4.69383E-02 4.70639E-02 4.71820E-02 4.72932E-02 4.73978E-02 4.74964E-02 4.75892E-02 4.76767E-02 4.77591E-02 4.78369E-02 4.79102E-02 4.79795E-02 4.80448E-02 4.81065E-02 4.81648E-02 4.82198E-02 4.82719E-02 4.83210E-02 4.83675E-02 4.84115E-02 4.84531E-02 4.84924E-02 4.85296E-02 4.85648E-02) ; lqc = lqc 'ET' ('PROG' 4.85982E-02 4.86297E-02 4.86596E-02 4.86879E-02 4.87147E-02 4.87400E-02 4.87640E-02 4.87868E-02 4.88084E-02 4.88288E-02 4.88481E-02 4.88664E-02 4.88838E-02 4.89003E-02 4.89158E-02 4.89306E-02 4.89446E-02 4.89578E-02 4.89704E-02 4.89822E-02 4.89935E-02 4.90041E-02 4.90142E-02 4.90237E-02 4.90327E-02 4.90413E-02 4.90493E-02 4.90569E-02 4.90641E-02 4.90709E-02 4.90774E-02) ; tic = rxt . 'TIC' ; ERQc = 'SOMM' ('ABS'(lqc - tic . 'Qc' )) ; 'MESS' 'ERQc=' ERQc ; 'SI' (ERQc '>' 1.e-4) ; err1 = err1 '+' 1 ; 'FINS' ; 'SI' ('NEG' ERR1 0) ; 'ERRE' 5 ; 'FINS' ; 'FINS' ; * * * Plots devoted to pressu... cases * * 'SI' GRAPH ; tbt = rxt . 'TBT' ; tic = rxt . 'TIC' ; * $vtf = rxt . 'GEO' . '$vtf' ; vtf = 'DOMA' $vtf 'MAILLAGE' ; * Mpl1 = 'CHAN' 'QUAF' plan1 ; Mpl4 = 'CHAN' 'QUAF' plan4 ; 'ELIM' (vtf 'ET' Mpl1 'ET' Mpl4) epsi ; $mpl1 = 'MODE' Mpl1 'NAVIER_STOKES' 'MACRO' ; $mpl4 = 'MODE' Mpl4 'NAVIER_STOKES' 'MACRO' ; plan1 = 'DOMA' $mpl1 'MAILLAGE' ; plan4 = 'DOMA' $mpl4 'MAILLAGE' ; plan = plan1 'ET' plan4 ; cplan = 'CONT' plan ; * 'SI' ('EXIS' tic 'TP') ; $vtp = rxt . 'GEO' . '$vtp' ; vtp = 'DOMA' $vtp 'MAILLAGE' ; 'FINS' ; paroif = rxt . 'GEO' . 'paroif'; cparoif = 'CONT' paroif ; * axe = p0 d nz1 (p0 plus v1) ; axe = 'CHAN' axe 'QUAF' ; 'ELIM' (axe 'ET' mt) epsi ; * un = tic . 'UN'; unp = 'REDU' un plan ; ung = 'VECT' un 0.5 'UX' 'UY' 'UZ' 'JAUN' ; ungp = 'VECT' unp 0.5 'UX' 'UY' 'UZ' 'JAUN' ; tf = tic . 'TF' ; rho = tic . 'RHO' ; rair = tic . 'RAIR' ; 'SI' tbt . 'THE' ; rhe = tic . 'RHE' ; 'FINS' ; 'SI' tbt . 'TH2' ; rh2 = tic . 'RH2' ; 'FINS' ; 'SI' tbt . 'TCO' ; rco = tic . 'RCO' ; 'FINS' ; 'SI' tbt . 'TCO2' ; rco2 = tic . 'RCO2' ; 'FINS' ; 'SI' tbt . 'VAPEUR' ; rvp = tic . 'RVP' ; 'FINS' ; * evauz = 'EVOL' 'CHPO' ('EXCO' un 'UZ') axe ; 'DESS' evauz 'TITR' 'Velocity with the z axis' 'MIMA' 'GRIL' 'POIN' 'GRIS' 'TITX' 'z' 'TITY' ' m/s' ; 'TRAC' ung plan ('CONT' plan) 'TITR' ' Velocity' ; * evatf = 'EVOL' 'CHPO' tf axe ; 'DESS' evatf 'TITR' 'Gas temperature with the z axis' 'MIMA' 'GRIL' 'POIN' 'GRIS' 'TITX' 'z' 'TITY' ' C' ; 'TRAC' tf plan cplan 'TITR' ' Temperature' ; 'TRAC' tf paroif cparoif 'TITR' ' Temperature' ; * evarh = 'EVOL' 'CHPO' rho axe ; 'DESS' evarh 'TITR' 'Gas density with the z axis' 'MIMA' 'GRIL' 'POIN' 'GRIS' 'TITX' 'z' 'TITY' ' kg/m3' ; 'TRAC' rho plan ('CONT' plan) ungp 'TITR' ' Density & velocity' ; 'TRAC' rho paroif cparoif 'TITR' ' Density' ; * evavap = 'EVOL' 'CHPO' rvp axe ; 'DESS' evavap 'TITR' 'Steam density with the z axis' 'MIMA' 'GRIL' 'POIN' 'GRIS' 'TITX' 'z' 'TITY' ' kg/m3' ; 'TRAC' rvp plan cplan 'TITR' ' Steam density' ; * Fcond = rxt . 'TIC' . 'Fcondw'; 'TRAC' fcond paroif cparoif 'TITR' ' Fcond kg/m**2' ; * * Wall temperature 'SI' ('EXIS' tic 'TP') ; 'TRAC' tic . 'TP' vtp 'TITR' ' Wall temperature' ; 'FINS' ; * * * Specific plots * * TAB1 = TABLE ; TAB1 . 1 = 'MARQ LOSA REGU' ; TAB1 . 2 = 'MARQ TRIU REGU' ; TAB1 . 3 = 'MARQ PLUS TIRM' ; TAB1 . 4 = 'MARQ CARR TIRM' ; TAB1 . 'TITRE' = TABLE ; * *---------------------------------------> Hext, Text , ... TAB1 . 'TITRE' . 1 = 'MOT' 'TP int' ; TAB1 . 'TITRE' . 2 = 'MOT' 'TP ext' ; TAB1 . 'TITRE' . 3 = 'MOT' 'Text' ; TAB1 . 'TITRE' . 4 = 'MOT' 'TF int' ; evtpi = ('EVOL' 'CHPO' tic . 'TP' lzi) 'COUL' 'ROUG' ; evtpe = ('EVOL' 'CHPO' tic . 'TP' lze) 'COUL' 'BLEU' ; evtpxt = ('EVOL' 'CHPO' tic . 'TPEXT' lze) 'COUL' 'BLEU' ; evtfi = ('EVOL' 'CHPO' tic . 'TF' lzi) 'COUL' 'ROUG' ; 'DESS' (evtpi 'ET' evtpe 'ET' evtpxt 'ET' evtfi) 'TITR' 'Gas, wall and external temperatures (KAS2)' 'MIMA' 'GRIL' 'POIN' 'GRIS' 'TITX' 'z(m)' 'TITY' ' C' 'LEGE' tab1 ; TAB1 . 'TITRE' . 1 = 'MOT' 'Hext ' ; hext = 'ELNO' tic . 'HEXTc' rxt . 'GEO' . '$parext' ; evhxt = 'EVOL' 'CHPO' hext lze ; 'DESS' evhxt 'TITR' 'Hext (KAS2)' 'MIMA' 'GRIL' 'POIN' 'GRIS' 'TITX' 'z(m)' 'TITY' 'W/m^2/K' 'LEGE' tab1 ; 'TRAC' hext pare 'TITR' 'Hext (KAS2)' ; 'TRAC' tic . 'TPEXT' pare 'TITR' 'Text (KAS2)' ; * *---------------------------------------> Pressure * Why TPi and not TFi for air ans steam ? TAB1 . 'TITRE' . 1 = 'MOT' 'PT' ; TAB1 . 'TITRE' . 2 = 'MOT' 'Psat(TPi)' ; TAB1 . 'TITRE' . 3 = 'MOT' 'Pair(TPi)' ; TAB1 . 'TITRE' . 4 = 'MOT' 'Pvap(TPi)' ; Psattp = PSATT (tic . 'TP' + 273.15) ; PT = 'EXTR' ('DIME' tic . 'PT') (tic . 'PT') ; Pair = ('EXTR' tic . 'Rhomn' 1) '*' 287.1 '*' tic . 'TP' ; pvi = PT - Pair ; evpsattpi = 'EVOL' 'CHPO' Psattp lzi ; evpt = 'EVOL' 'MANU' ('PROG' 0. h1) ('PROG' PT PT) ; evpair = 'EVOL' 'CHPO' pair lzi ; evpvi = 'EVOL' 'CHPO' pvi lzi ; 'DESS' (evPT et evpsattpi et evpair et evpvi) 'TITR' 'Pressure' 'MIMA' 'GRIL' 'POIN' 'GRIS' 'TITX' 'z(m)' 'TITY' ' Pa' 'LEGE' tab1 ; * *---------------------------------------> Density TAB1 . 'TITRE' . 1 = 'MOT' 'Rho' ; TAB1 . 'TITRE' . 2 = 'MOT' 'Rvp' ; TAB1 . 'TITRE' . 3 = 'MOT' 'Rnc' ; TAB1 . 'TITRE' . 4 = 'MOT' 'Rsat(TPi)' ; Rsattp = Psattp '/' 461.73 '/' (tic . 'TP' '+' 273.15) ; evrhpi = 'EVOL' 'CHPO' tic . 'RHO' lzi ; evrvpi = 'EVOL' 'CHPO' tic . 'RVP' lzi ; evrinci = 'EVOL' 'CHPO' (tic . 'RHO' '-' tic . 'RVP') lzi ; evrsati = 'EVOL' 'CHPO' Rsattp lzi ; 'DESS' (evrhpi 'ET' evrvpi 'ET' evrinci 'ET' evrsati) 'TITR' 'Density at the wall' 'MIMA' 'GRIL' 'POIN' 'GRIS' 'TITX' 'z(m)' 'TITY' 'kg/m3' 'LEGE' tab1 ; * *---------------------------------------> Mass fraction TAB1 . 'TITRE' . 1 = 'MOT' 'Yvap' ; TAB1 . 'TITRE' . 2 = 'MOT' 'Ysat(TPi)' ; TAB1 . 'TITRE' . 3 = 'MOT' 'Yvap - Ysat(TPi)' ; yvi = tic . 'RVP' '/' tic . 'RHO' ; ysati = ('REDU' Rsattp lzi) '/' ('REDU' tic . 'RHO' lzi) ; evysati = 'EVOL' 'CHPO' ysati lzi ; evyvi = 'EVOL' 'CHPO' yvi lzi ; evdyv = evyvi '-' evysati ; 'DESS' (evyvi 'ET' evysati 'ET' evdyv) 'TITR' 'Steam mass fraction at wall' 'MIMA' 'YBOR' 0. 0.4 'GRIL' 'POIN' 'GRIS' 'TITX' 'z(m)' 'TITY' ' Mass Fraction' 'LEGE' tab1 ; * *- Fin GRAPH * 'FINSI' ; 'FIN' ;