closed loop under pulsed loads B. Rousset, C.Hoa, B. Lagier, - - PowerPoint PPT Presentation

SLIDE 1

0-D approach to predict P and T along HELIOS SHe closed loop under pulsed loads B. Rousset, 14 October 2011

BASSES TEMPERATURES

0-D thermo hydraulic approach for predicting pressure and temperature along HELIOS SHe closed loop under pulsed loads

• B. Rousset, C.Hoa, B. Lagier, R.Vallcorba

SLIDE 2

2 0-D approach to predict P and T along HELIOS SHe closed loop under pulsed loads B. Rousset, 14 October 2011

1 Build a simple O-D tool able to describe the time dependent thermo physical properties (Pi(t), Ti(t)) of a flow submitted to pulsed loads. Depending on the accuracy of the 0-D model, it can be implemented in open code as EcosimPro, … Furthermore, it will be an help to analyze/predict the response of pulse loads as the CPU time required would be negligible. 2 A such simple model could be used together with iterative calculations to solve (with some additional assumptions) the inverse problem, e.g. find the time dependent power injected from a time dependent pressure evolution. As pressure reacts instantaneously (no delay due to transit time), this can be used to anticipate the pulse loads arrival at the heat exchanger location.

Objectives

SLIDE 3

3 0-D approach to predict P and T along HELIOS SHe closed loop under pulsed loads B. Rousset, 14 October 2011

P t

SLIDE 4

4 0-D approach to predict P and T along HELIOS SHe closed loop under pulsed loads B. Rousset, 14 October 2011

Just remember that pressure variation only depends on energy variation. So during the power injection, pressure increase is independent of mass flow whereas the plateau and the decreasing time only depend on transit time from the start of the heated zone to downstream the heat exchanger, i.e. mass flow.

Time dependent pressure profile

SLIDE 5

5 0-D approach to predict P and T along HELIOS SHe closed loop under pulsed loads B. Rousset, 14 October 2011

Tout t Tout P

SLIDE 6

6 0-D approach to predict P and T along HELIOS SHe closed loop under pulsed loads B. Rousset, 14 October 2011

Simplified HELIOS cooling scheme

T3 T2 T1 31.5 m 25 m 25 m 8 m 8 m 1 2 3 Heat exchanger with saturated bath Circulating pump P

HELIOS is relevant to an isochoric system including a forced flow, some heaters allowing heat pulses and at least one heat exchanger allowing heat pulses to be removed from the loop.

Heat exchanger with saturated bath Heater Heater Heater

SLIDE 7

7 0-D approach to predict P and T along HELIOS SHe closed loop under pulsed loads B. Rousset, 14 October 2011

Problem to be solved

Input data : mass flow rate, loop geometry and time dependent pulse powers Pi(t) Ouput : Find the time dependent pressure and temperature profiles

SLIDE 8

8 0-D approach to predict P and T along HELIOS SHe closed loop under pulsed loads B. Rousset, 14 October 2011

We have seen that transit time seems an important parameter. This transit time depends on considered component geometry. Is there a simple way to express component transit time ? Can we choose length and velocity : ? Not a good choice ! Right choice implies volume and volumetric flow rate : Proceeding so, result becomes quasi universal and can be calculated easily. For example with a density of 132 Kg/m3, 10 l and 100 g/s give 13.2 s (and the same for 100 l and 1 kg/s or 132 s if you consider 1 m3 and 1 kg/s, …) So remind to say this sensor is located 100 liters downstream the inlet and not 12 m for example !

Preliminary remarque

m V Q V t

V

    

v l t 

SLIDE 9

9 0-D approach to predict P and T along HELIOS SHe closed loop under pulsed loads B. Rousset, 14 October 2011

0-D resolution sequence based on superposition principle

Steady state calculation Time dependent pressure profile calculation Time dependent enthalpy change induced by pressure variation (calculation at the exit of each component) Time dependent enthalpy change induced by heat input or removed (calculation at the exit of each component)

T(xi,t) P(x,t)

Split the loop in various 0-D space components (volume and time dependent power for each component)

SLIDE 10

10 0-D approach to predict P and T along HELIOS SHe closed loop under pulsed loads B. Rousset, 14 October 2011

A component is defined by its volume and eventually the time dependent power in case of heated sector. Each heated sector must be considered as a specific component. So the simplest decomposition consists to build components corresponding to heated sectors and non-heated sectors. T emperature must be known at the inlet of one component and will be calculated at its outlet. Consequently, to determine a temperature at a specific abscissa, a component must have its outlet at this abscissa. T

• do

this 1 component can be split in 2 components.

Choice of volumetric 0-D components

T1 T3 T2 V1 V2 W2(t) V3 V4 V5 W5(t) V6

SLIDE 11

11 0-D approach to predict P and T along HELIOS SHe closed loop under pulsed loads B. Rousset, 14 October 2011

ttr1 ttr2=V5/Q ttr2_hx=V6/Q ttr1_hx This example shows that the knowledge of the pressure increase as a function of the heat injected (W1+W2) is sufficient to determine the time dependent pressure profile !

V is the component volume

Q is the volumetric flow rate

P t

tch1 0 tch2 + tch2 ttr2_hx ttr2_hx + ttr2 + tch1 ttr1_hx ttr1_hx + ttr1 W2 W1 W1 W2

Time dependent pressure profile : example with two heaters

V6 V5

SLIDE 12

12 0-D approach to predict P and T along HELIOS SHe closed loop under pulsed loads B. Rousset, 14 October 2011

Thermal loads due to pump work and heat losses are taken into account during the steady state calculation and are supposed to be not affected by the transient loads. The loops is then considered as a closed volume submitted to an isochoric

• transformation. Applying first principle gives :

Mass internal energy increase is thus equal to : Finally the pressure is calculated using the Hepak code using density and internal energy as input data, the former one being constant in this isochoric process. dQ dW dQ dU   

tot t t t t t t

M Wdt U U



   

 

 

tot t t t t t t t t t t t t t t t

M Wdt U U and with U P



           

        ,

Calculation of the time dependent pressure increase

SLIDE 13

13 0-D approach to predict P and T along HELIOS SHe closed loop under pulsed loads B. Rousset, 14 October 2011

Time dependent pressure profile : calculated and experimental results

3 simultaneous pulses of 333 Watt - 18 seconds and a mass flow rate of 32 g/s

4 4.5 5 5.5 6 6.5 100 150 200 250 300 350 400 450 500 550 600 650 700 Time (s) Pressure (bar) Experimental data 0-D model

SLIDE 14

14 0-D approach to predict P and T along HELIOS SHe closed loop under pulsed loads B. Rousset, 14 October 2011

Time dependent pressure profile : calculated and experimental results

3 simultaneous pulses of 250 Watt - 60 seconds and a mass flow rate of 32 g/s

4 4.5 5 5.5 6 6.5 7 7.5 8 8.5 9 358 408 458 508 558 608 658 708 758 808 858 908 958 time (s) Pressure (Bar) Experimental data 0-D model

SLIDE 15

15 0-D approach to predict P and T along HELIOS SHe closed loop under pulsed loads B. Rousset, 14 October 2011

Time dependent temperature profile will be different along the line and the global approach adopted for the pressure cannot be used here. Furthermore

variation of pressure as well as heat pulse has large impact on temperature evolution and each contribution must be considered.

Appling once again the superposition principle, we will assume that each contribution can be calculated separately and summed afterwards. Finally, temperature evolution will also depend on temperature profile existing upwind previously (convective effect). It is assumed that temperature at the inlet of the line has a constant value (equal to bath temperature + a small delta T). For the point considered, we will first calculate the contribution of the

pressure evolution.

Calculation of the time dependent temperature profile

SLIDE 16

16 0-D approach to predict P and T along HELIOS SHe closed loop under pulsed loads B. Rousset, 14 October 2011

For points located upstream of the heaters (e.g. T1 or T2), the only contribution to temperature change will be pressure change.

Time dependent temperature profile induced by pressure variation

T3 T2 T1 31.5 m 25 m 25 m 8 m 8 m 1 2 3 HX1 Heat exchanger with saturated bath Circulating pump P HX2 Heat exchanger with saturated bath Heater Heater Heater

SLIDE 17

17 0-D approach to predict P and T along HELIOS SHe closed loop under pulsed loads B. Rousset, 14 October 2011

Compression induced by energy injection (resp. pressure discharge induced by energy extracted from the loop) will heat (resp. cool down) fluid inside the loop. Variation of pressure can be thus considered as heater (or cold source) uniformly distributed along the loop. Furthermore, for a constant pressure evolution, temperature change will be limited by the transit time between the inlet heat exchanger (HX1) where outlet temperature is kept constant and the point considered (T2).

Time dependent temperature profile induced by pressure variation

SLIDE 18

18 0-D approach to predict P and T along HELIOS SHe closed loop under pulsed loads B. Rousset, 14 October 2011

Some examples of temperature response to a pressure gradient are shown on following figures. ttransit_hxinlet

P t T2 t Time dependent temperature profile induced by pressure variation

T1 T3 T2 V1 V2 V3 W3(t) V4

ttransit_hxinlet

Pseudo heater due to compression

SLIDE 19

19 0-D approach to predict P and T along HELIOS SHe closed loop under pulsed loads B. Rousset, 14 October 2011

P t P t

tP

P t + =

tP

T t

tP ttransit_hxinlet + tP ttransit_hxinlet

T t

ttransit_hxinlet

T t

tP ttransit_hxinlet + tP

t + =

tP ttransit_hxinlet <

Time dependent temperature profile : pressure effect

SLIDE 20

20 0-D approach to predict P and T along HELIOS SHe closed loop under pulsed loads B. Rousset, 14 October 2011

P t P t

tP

P t + =

tP

T t

tP ttransit_hxinlet + tP ttransit_hxinlet

T t

ttransit_hxinlet

T t

tP ttransit_hxinlet + tP

t + =

tP ttransit_hxinlet >

Time dependent temperature profile : pressure effect

SLIDE 21

21 0-D approach to predict P and T along HELIOS SHe closed loop under pulsed loads B. Rousset, 14 October 2011

T

• calculate the influence of pressure variation on temperature, we will

consider an isentropic evolution. An isentropic evolution means an evolution both adiabatic and reversible. As the superposition principle will be applied, the contribution of the heat on temperature profile is not taken into account at this stage and the assumption of adiabatic evolution is logical. Furthermore, considering the fact that a compression followed by a discharge to obtain same value of pressure will give “in fine” same value of temperature, reversibility is also logical. The code Hepak is then used to calculate the time dependent temperature variation induced by pressure change :

 

t dt t dt t

S P T T ,

  

Time dependent temperature profile induced by pressure variation

SLIDE 22

22 0-D approach to predict P and T along HELIOS SHe closed loop under pulsed loads B. Rousset, 14 October 2011

Comparison between experimental data and 0-D model

4.4 4.45 4.5 4.55 4.6 4.65 4.7 4.75 4.8 4.85 4.9 120 140 160 180 200 220 time (s) Temperature (K) Experimental data 0-D model

T1 (transit time from inlet : 5s) T2 (transit time from inlet : 74s) T3 (transit time from inlet : 480s)

4 4.5 5 5.5 6 120 130 140 150 160 170 180 190 200 210 220 Time (s) Pressure (bar) Experimental data 0-D model

3 simultaneous pulses of 333 Watt - 18 seconds and a mass flow rate of 32 g/s

T3 T2 T1 31.5 m 25 m 25 m 8 m 8 m 1 2 3 Heat exchange saturated ba Circulating pump P Heat exchanger with saturated bath Heater Heater Heater

SLIDE 23

23 0-D approach to predict P and T along HELIOS SHe closed loop under pulsed loads B. Rousset, 14 October 2011

In a 0-D approach, the temperature change will be linked to mass inside the volume and enthalpy change. T

• take into account the convection effect, enthalpy at the outlet of the volume

must be linked to the value of enthalpy at the inlet with a delay equal to the transit time inside the volume. So for a non heated volume, the relation is simply While for a heated volume of mass M with a power W initiated at tinit, the resulting equation becomes :

Time dependent temperature profile induced by heat variation

   

transit inlet

• utlet

t t H t H  

Remarque : in any case enthalpy increase along the volume due to steady state heat losses will be added before converting enthalpy in temperature

   

M Wdxdt L t t H t H

t t L transit inlet

• utlet

init

 

   1

SLIDE 24

24 0-D approach to predict P and T along HELIOS SHe closed loop under pulsed loads B. Rousset, 14 October 2011

Considering a constant spatial repartition and a square temporal power pulse, the previous equation gives :

Time dependent temperature profile induced by heat variation

Remarque 1 : at this stage, we should say that this relation is only valid for a time lower than the transit time. If the power is maintain for a duration larger than the transit time, then the steady state condition is reached for the considered volume and temperature at its outlet becomes constant

     

M t t W t t H t H

init transit inlet

• utlet

   

Remarque 2 : when the heater is turn off, the outlet enthalpy remains constant for a duration equal to the transit time leading to the following equation :

     

M t t W t t H t H

init OFF transit inlet

• utlet

   

Remarque 4 : After the heater is turn off from more than the transit time, the pulse is completely evacuated and the equation is the same as a non heated volume :

   

transit inlet

• utlet

t t H t H  

Remarque 4 : After, the enthalpy decreases with the opposite slope of its increase period leading to the following equation :

       

M t t t t W t t H t H

init transit OFF transit inlet

• utlet

     

OFF init

t t t For  

transit OFF

t t t For  

 

init OFF transit transit

t t t t t For    

 

init OFF transit

t t t t For   

SLIDE 25

25 0-D approach to predict P and T along HELIOS SHe closed loop under pulsed loads B. Rousset, 14 October 2011

Time dependent temperature profile induced by heat variation T0

0 tOFF ttransit + tOFF tOFF

t

ttransit

W t

T emperature profile at the outlet of the heated sector (inlet condition assumed to be constant)

T1 T3 T2 V1 V2 V3 W3(t) V4

ttransit

T0

SLIDE 26

26 0-D approach to predict P and T along HELIOS SHe closed loop under pulsed loads B. Rousset, 14 October 2011

Time dependent temperature profile induced by heat variation T0

0 tOFF ttransit + tOFF tOFF

t

ttransit

W t T3

tOFF ttransit + tOFF

t

ttransit + ttransit2 + ttransit2 + ttransit2 ttransit2

T1 T3 T2 V1 V2 V3 W3(t) V4

ttransit

T0

ttransit2

SLIDE 27

27 0-D approach to predict P and T along HELIOS SHe closed loop under pulsed loads B. Rousset, 14 October 2011

Comparison between experimental data and 0-D model

4 4.5 5 5.5 6 6.5 100 150 200 250 300 350 400 450 500 550 600 650 Time (s) Pressure (bar) Experimental data 0-D model

4.5 4.6 4.7 4.8 4.9 5 5.1 5.2 5.3 5.4 100 150 200 250 300 350 400 450 500 550 600 650 Time (s) T° T3 (K) Experimental data 0-D model

SLIDE 28

28 0-D approach to predict P and T along HELIOS SHe closed loop under pulsed loads B. Rousset, 14 October 2011

Summary of 0-D sequence of resolution

Steady state calculation Time dependent pressure profile calculation using internal energy variation (trapezoidal shape due to transit time assumed) at constant density Time dependent enthalpy change induced by pressure variation assuming isentropic evolution Time dependent enthalpy change induced by heat input or removed assuming isobaric evolution

T(x,t) P(x,t)

Split the loop in various components

4 4.5 5 5.5 6 6.5 100 150 200 250 300 350 400 450 500 550 600 650 700 Time (s) Pressure (bar) Experimental data 0-D model

4 4.5 5 5.5 6 6.5 7 7.5 8 8.5 9 358 408 458 508 558 608 658 708 758 808 858 908 958 time (s) Pressure (Bar) Experimental data 0-D model

SLIDE 29

29 0-D approach to predict P and T along HELIOS SHe closed loop under pulsed loads B. Rousset, 14 October 2011

0-D sequence of resolution : first improvement

Choice of loop components + Steady state calculation

T(x,t) P(x,t)

 

 

 ;

t loop

U P

   

 

;

t t

• utlet

S P H

   

 

   

 

 

 

let transitout t t t

• utlet

t t

• utlet

W H S P H S P H



  

'

; ;

 

) ( ) ( ' ) ( ) (

;

t

• utlet

t t

• utlet
• utlet

t

U P H U U   

tot t t t t loop t t loop

M Wdt U U /

) ( ) (



   

 

     

t t loop t loop

U U U   

'

 

 

 ;

' '

t loop

U P

 

) ( ' ) ( ; t t x x

P H S

 

 

) ( ' ; t x x

S P T

t

   

 

   

 

 

 

transitx t t t x t t x

W H S P H S P H



  

'

; ;

SLIDE 30

30 0-D approach to predict P and T along HELIOS SHe closed loop under pulsed loads B. Rousset, 14 October 2011

Comparison between experimental data and improved 0-D model

4 4.5 5 5.5 6 6.5 7 7.5 8 8.5 350 450 550 650 750 850 950 Time (s) Pressure (bar) Experimental data 0-D model first approach 0-D after improvement

3 simultaneous pulses of 250 Watt - 60 seconds and a mass flow rate of 32 g/s

SLIDE 31

31 0-D approach to predict P and T along HELIOS SHe closed loop under pulsed loads B. Rousset, 14 October 2011

Comparison between experimental data and improved 0-D model

T3 temperature profile 4.4 4.6 4.8 5 5.2 5.4 5.6 5.8 6 6.2 6.4 350 450 550 650 750 850 time (s) Temperature (K) experimental data 0-D model after improvement

3 simultaneous pulses of 250 Watt - 60 seconds and a mass flow rate of 32 g/s

SLIDE 32

32 0-D approach to predict P and T along HELIOS SHe closed loop under pulsed loads B. Rousset, 14 October 2011

Conclusions :

A very simple 0-D non iterative approach can be used to understand and reproduce correctly the thermohydraulic behavior of a forced flow submitted to distributed heat pulses. This model may be improved by coupling pressure and temperature effects with an iterative procedure. Among the limitations of this model, it was not able to reproduce the sub-cooling induced by the pressure decrease. Finally, latest improvement shows that it could be easily possible to add the saturated bath component with interaction of energy injected inside this bath on both the loop and the bath in order to analyze thermal buffer effect.

SLIDE 33

33 0-D approach to predict P and T along HELIOS SHe closed loop under pulsed loads B. Rousset, 14 October 2011

Thank you for your attention

SLIDE 34

34 0-D approach to predict P and T along HELIOS SHe closed loop under pulsed loads B. Rousset, 14 October 2011

1) Mass and volumetric flow rates will be considered as constant. 2) Spatial pressure gradient will be considered equal to the steady state case. 3) Power pulses considered here have a constant spatial repartition (W/m) and a square time profile. 4) Heat diffusion through helium or pipe is considered as negligible 5) Equations describing temperature and pressure evolution can be calculated independently, with a superposition principle applied afterward 6) DT between the saturated bath and the outlet of the heat exchanger is supposed to be constant 7) Saturated bath temperature is supposed to be constant

Assumptions :

SLIDE 35

35 0-D approach to predict P and T along HELIOS SHe closed loop under pulsed loads B. Rousset, 14 October 2011

Steady state calculation

Inlet line temperature, mass flow rate, heat losses and geometry of the line are used in the momentum balance and energy balance equation to determine pressure and temperature in each abscissa. Typically, pressure drop is calculated using some correlations (colebrook correlation for regular pressure drop and adapted correlation for singular pressure drop as elbows or valves), while energy balance is considered as enthalpy balance. The resolution of pressure and enthalpy in each point is sufficient to solve the thermodynamic problem as these variables are independent.

SLIDE 36

36 0-D approach to predict P and T along HELIOS SHe closed loop under pulsed loads B. Rousset, 14 October 2011

P t P t

tP1

P t + =

tP2

T t

tP2 ttr +tP2 ttr

T t

ttransit_hxinlet

T t

tP1 ttr +tP1

+ = Time dependent temperature profile : pressure effect P t

tP2 tP1

T t

ttr +tP2

t + t +

tP1 tP2 ttr +tP1

P t

tP3 tP3

+ T

tP3 ttr + tP3 = ttr

+

tP3 ttr +tP3

SLIDE 37

37 0-D approach to predict P and T along HELIOS SHe closed loop under pulsed loads B. Rousset, 14 October 2011

P t

SLIDE 38

38 0-D approach to predict P and T along HELIOS SHe closed loop under pulsed loads B. Rousset, 14 October 2011

Tout t Tout P

SLIDE 39

39 0-D approach to predict P and T along HELIOS SHe closed loop under pulsed loads B. Rousset, 14 October 2011

Parmi les nombreuses remarques énergétiques que l’on peut faire, on notera :

• On obtient des extrema de pression (minima ou maxima) à chaque fois que

l’énergie extraite dans la boucle est égale à l’énergie extraite en régime stationnaire (ce qui se calcule facilement en additionnant les bilans mDH aux bornes des 2 échangeurs).

• Pendant le pulse de puissance, la température en entrée des échangeurs

augmente par effet de montée de pression, ce qui revient à extraire plus de puissance de la ligne et donc à diminuer un peu l’impact du pulse de chaleur en terme d’accroissement de pression.

Conclusions :

SLIDE 40

40 0-D approach to predict P and T along HELIOS SHe closed loop under pulsed loads B. Rousset, 14 October 2011

Rmq1 : la connaissance de la montée en pression en fonction du temps peut donc dans la plupart des cas être un indicateur de l’énergie déposée sur la ligne en fonction du temps, par contre cela n’est pas suffisant pour remonter à la répartition spatio-temporelle de la puissance injectée au cours du temps (deux secteurs de longueurs différentes et éventuellement disposés à des endroits différents de la ligne) pouvant ainsi induire la même montée en pression si la puissance totale déposée sur chaque secteur est identique et synchrone.

Evolution temporelle de la pression: remarque finale

Rmq2 : Pour retrouver le profil spatio-temporelle de puissance, on peut s’aider du temps de plateau et de la courbe de descente, qui, avec la connaissance du débit devraient permettre de localiser le début et la fin de la zone chauffée, ce qui permettrait ensuite de fixer de manière non ambigüe la valeur de celle-ci au cours du temps.