MINISTERE DE L'INDUSTRIE

BUREAU DE RECHERCHES GOLOGIQUES ET MINIRES

SERVICE GOLOGIQUE NATIONAL
B. P. 6009 - 45018 Orlans Cedex -Tl.

(38) 63.80.01

MODLISATION DU COMPORTEMENT
HYDRAULIQUE ET THERMIQUE
DES RSEAUX MAILLS

APPLICATION AU RSEAU D'EAU THERMALE

DE LA VILLE DE DAX

-2.V- ':.} 33A

SERVICE GEOLOGIQUE REGIONAL AQUITAINE

Avenue Docteur-Albert-Schweitzer, 3 3 6 0 0 PESSAC - Tl. (56) 8 0 . 6 9 . 0 0

Rapport du B . R . G . M .
84 AGI 164 AQI

Pessac, le 29 juin 1 9 8 4

MINISTRE DE L'INDUSTRIE ET DE LA RECHERCHE

BUREAU DE RECHERCHES GOLOGIQUES ET MINIRES
SERVICE GOLOGIQUE NATIONAL
B.P. 6009 - 45060 Orlans Cedex - Tl.: (38) 63.80.01

MODELISATION DU COMPORTEMENT
HYDRAULIQUE ET THERMIQUE
DES RESEAUX MAILLES

Ph.DAUGUET
JC.MARTIN
84 AGI 164 AQI \

Pessac le 29 juin 1984

SERVICE GOLOGIQUE RGIONAL AQUITAINE

Avenue Docteur Albert Schweitzer - 33600 PESSAC - Tl.: (56) 80.69.00

S O M M A I R E

1 - INTRODUCTION

2 - DESCRIPTION D'UN RESEAU,GENERALITES

3 - ETUDE HYDRAULIQUE
3.1 - MODELISATION DES PERTES DE CHARGES
a) Pertes de charges dans une conduite simple
b) Calcul de.'la rsistance l'coulement R
c ) Rcapitulatif

7
8
10

3.2 - MODE DE CALCUL DE LA DISTRIBUTION DES DEBITS

a) Etablissement du systme d'quations
b} Exemple
c ) Rsolution du systme non linaire
d) Amlioration de la convergence

16
19
21
23

3.3 - METHODE DE CALCUL DE LA PIEZOMETRIE

a) Units et terminologie employes
b) Calcul des cotes pizometriques aux noeuds
c ) Calcul des pressions effectives
3.4 - CAS DES CONDUITES

24
25
2.6

SINGULIERES

a) Dbit "en route"

28

b ) Vannes et surpresseurs

30

3.5 - VALIDITE DES RESULTATS

32

3. - Extension possible aux rseaux de vapeur

32

4 - ETUDE THERMIQUE
4.1 - MODELISATION DES DEPERDITIONS THERMIQUES
a) Dperdition dans une conduite enterre
b ) Evolution de la temprature
c) Rcapitulatif

33
36
38

4.2 - METHODE DE CALCUL DE LA DISTRIBUTION DES TEMPERATURES

a ) Exemple
b) Gnralisation
c) Validit des rsultats

39
40
42

5 - PROGRAMMATION
5 . 1 - GENERALITES
a)
b)
c)
d)
e)

Utilisation des fichiers disques

Manipulation des diffrents fichiers
Structure du logiciel
Synoptique
Description des variables communes

43
44
47
48
49
50

5.2 - SAISIE ET STOCKAGE DES DONNEES

Programme "SAISIE"
Sous programme "RECDAT"
Fichier "DONNEES-D"

51
54
55

5.3 - MODIFICATION DES DONNEES

Programme "MODIF"

50

5.4 - ETABLISSEMENT DE LA MATRICE HYDRAULIQUE

Programme "MATRIX"

57

5.5 - CALCUL DE LA DISTRIBUTION DES DEBITS

Programme "DEBIX"

59

5.6 - CALCUL DE LA DISTRIBUTION DES PRESSIONS

Programme "PIEZO"

61

5.7 - CALCUL DE LA DISTRIBUTION DES TEMPERATURES

Programme "THERMO"

63

5.8 - INVERSION ET PRODUIT DE MATRICES

Sous programmes " INV" et "PROD"

<54

BIBLIOGRAPHIE

65

CONCLUSIONS

66

ANNEXE 1 : Logiciel LAUREL, Notice d'utilisation.

ANNEXE 2 : Exploitation sur le rseau gothermal de la ville de DAX.

1 - INTRODUCTION

1 - INTRODUCTION

Dans le cadre d'un contrat d'tude du B.R.G.M. avec la ville

de DAX (Landes) portant sur le rseau gothermal de cette ville, il est
apparu l'importance de possder un outil de calcul systmatique, d'utilisation locale (sur micro-ordinateur) permettant la simulation des
comportements hydraulique et thermique d'un tel rseau, ceci en vue
d'optimiser le rseau existant et surtout son extension future.

Le problme de dtermination des dbits dans un rseau est

d'une formulation thorique relativement aise mais sa rsolution se pose
de manire trs complexe ds l'instant o l'on a affaire un rseau maill,
ce qui est gnralement le cas dans la plupart des rseaux de distribution.

Ce problme, au moins pour sa partie hydraulique, a donn lieu

quelques mthodes de rsolution dont la plus connue est sans doute celle
mise au point par

HARDY-CROSS en 1936 et qui est en fait une mthode de

rsolution de systme non linaire par relaxation. Celle-ci connat

d'ailleurs quelques applications sur ordinateur (notamment le logiciel
HARDY dvelopp par la CISI).

L'analyse des rseaux maills s'inscrit en effet parfaitement

dans le champ d'action des techniques informatiques.

Le logiciel LAUREL, issu de.la prsente tude, est bas sur la

linarisation d'un systme par la mthode de Newton applique au calcul
matriciel. La partie thermique prsente l'originalit de faire galement
appel au calcul matriciel.

2 - DESCRIPTION D'UN RESEAU - GENERALITES

2 - DESCRIPTION D'UN RESEAU - GENERALITES

On appellera rseau un ensemble de conduites relies entre elles

et formant une srie de noeuds et de mailles.

- Un noeud est, soit un point de jonction de plusieurs conduites,

soit un point d'extrmit d'une ramification.

Les noeuds comprennent les points d'injection de dbits et les

points de soutirage (consommations ou rejets).

- On dira de manire simple qu'une conduite est une partie du

rseau comprise entre deux noeuds.

- Un ensemble de conduites formant une boucle ferme constitue

une maille du rseau (une maille ne comportera pas de vanne ferme).

Les rseaux ne comportant pas de mailles sont dits ramifis.

Il faut noter l'existence possible de vannes et de surpresseurs.

d<

i tes

p
~f

GENERALITES D'une manire gnrale, les hypothses de travail sont les

suivantes :
. Rgime hydraulique et rgime thermique tablis.
. Sauf dans le cas des vannes, les pertes de charge singulires
ne sont pas tudies spcifiquement mais sont incluses dans le
calcul des pertes de charge linaires.
. On nglige les dperditions thermiques singulires face aux
dperditions linaires.
. Le calcul des dbits est effectu sur la base des dbits aux noeuds
imposs (dbits extrieurs).

3 - ETUDE HYDRAULIQUE

3 - ETUDE HYDRAULIQUE

3.1 - MODELISATION DES PERTES DE CHARGE

a) Pertes de charge dans une conduite simple

jfjl_.
A' / / / r/- /S / S f r

Le rgime est suppos turbulent (1), la perte de charge par unit

de longueur est donne par la formule universelle :

A
p

: masse volumique du fluide

: vitesse de dbit (

: diamtre

( ** /

: est un coefficient

adimensionnel fonction de la rugosit.

~est calcul partir de la formule de COLEBROOK :

JL_ .

.o < 3 7 ^

2,51

( *"
avec k t coefficient de rugosit
VD
On rappelle que Re = (nombre de Reynolds)

Dans la prsente tude, la valeur de i est fournie par un sous-programm

utilisant un calcul itratif dont l'algorithme est le suivant :

On pose

X =

(1) En pratique} te rg-Lme laminaire se limite aux vitesses trs faibles3 petits
diamtres et viscosit forte.

Re, k, D

Xo

Valeur d'initialisation

X = - 2 log

( -

1U

3,7 D

+ -2li
Re

Xo = X

oui-

non
I

b) Calcul de la

rsistance 1'coulement R

On se donne les hypothses suivantes :

La perte de charge est exprime en hauteur d'eau h (unit m . C . E . > mtre de
colonne d'eau) :
4

g h ;

(L : longueur de la conduite)
Le fluide vhicul a une masse volumique voisine de l^eau (il s'agira d'eau
chaude) :

- La perte de charge est calcule partir du dbit

D'o la formule :

- (

rr g
h

= 0,0826. ^L 5 . q 2
D

ou encore :
h - Rq

(m.CE.)

avec :

(SI)

. Adaptation de la formule thorique

Pratiquement, afin de tenir compte des pertes de charge singulires
telles que les coudes

ou les raccords entre tronons, on adapte la formule

en diminuant l'exposant du dbit (les dbits considrs sont infrieurs

l'unit (m3/s)J.

avec

1,7<

^ 2

(Voir figure 3 . 1 . 1 l'influence du choix de e^ sur le calcul des pertes de charge)

Cette adaptation pratique constituera une adaptation globale pour

l'ensemble des conduites du rseau.

On aura une adaptation spcifique chaque conduite par la modulation

du coefficient de rugosit k de la conduite considre.
. k est fonction du matri'au de la conduite et de son ge (dpts sur
parois), sa valeur est difficilement

les

apprciable.

(Voir figure 3 . 1 . 2 l'influence de k sur le calcul de la rsistance)

REMARQUE
Le calcul de la distribution des dbits dans le rseau ncessite la
connaissance des rsistances.
0v3 on sait que :

R = J(Dt k, Re)
en toute rigueur3 R est donc fonction de la vitesse (Re =

YD

).

On montre qu'en faitt dans un domaine restreint des valeurs de 7 , la

variation de Re n'a que trs peu d'influence sur la valeur de la rsistance.
On se fixe donc pralablement une vitesse Vo qui correspond la vitesse la
plus probable dans le rseau.
Cette approximation est justifie dans la figure 3.1.3.

10

D'autre part, on a affaire une eau de temprature qui n'est pas

constante pour l'ensemble du rseau.
Or y (viscosit cinmatique) dpend de la temprature.
On fait alors la mme remarque que prcdemment en se fisant une
temprature To qui correspond la temprature la plus probable dans le
rseau. (Justification en figure 3.1.4.).
Finalement :

Rcapitulatif ;
La dmarche est la suivante :
Pour une conduite donne :

. Saisie du diamtre D
de la longueur L
de la rugosit k
Vo
. Calcul de Re avec Re = 7-7

(Vo et I/o vitesse et viscosit les plus probables dans le rseau)

. Calcul itratif de / \ (Re, D , k)

. Calcul de R ,

R = O,O826

-*~D

. Si ncessaire, on modifie la valeur de la rugosit k (valeur difficilement

apprciable) pour se rapprocher des rsultats exprimentaux (opration de calage)

clu choix

ck

P(

s<sr Its

ptrc/ de. cUa.tr

I\H

r.v*n\

Fiytrs
-C*T

0loo * = <?/
0 30

">

Vo

TQ

o,}-

-4-

15

Proprits physiques de l'eau douce la pression atmosphrique (g = 9 , 8 1 m / s s )

Temperaturei

Poids

Tension
superficieiie
(eau avec
l'air)
U
Centiatokes
kg m * 1

Mini
spcifique

Viscosit
dynamique

V-

kg/m3

999.9

1OI.93

181 X io-c

i,77 X io- e

1.77

0,00771

0,062

1,99 X 10

iooo,o

101,94

160 X 10-

1.57 Xi-

1.57

0,00763

0,083

IO

999.7

1OI.9I

1 3 4 X 1 0 c 1,31X10-

1,31

0,00757

0,125

2,09 x io"

20

998,a

ior,76

103 X IO-6

1,01

0,00743

0,239

2,18 X 10*

30

995,7

101.50

84 X io- 0,83X10-

0,82

0,00726

O.433

2,S0 X 10*

4
5

992,2

101,14

67 X 10- 6

0,66 X io-s

0,66

0,00710

O.753

2,21 X IO

983,1

100,72

56X10

0,55X10-6

O.55

0,00690

1,258

2,22 X 1OS

60

983,2

100,23

47 X 10-

0,46X10-0

0.46

0,00676

2,O33

2,23 X io"

80

97r,8

99,06

37 X IO 6 o , 3 7 X i o - n

o,37

0,00638

4.831

100

958,4

97.70

28 X 10- 0,29X10-6

0,29

O,O06o0

10,333

k m -

s-2

kg m '

Visco sit* cinmatique!

Tension
de vapeur

p ici (1 q u e

m*/.

I,OIXIO

Dans lea calculs hydrauliques ordinaires, on prend UT= 1000 k e / m 3 ;

kg/m3

m
(valeurs approches)
d'eau 4 * C

p 102 kg m " 4 s2; V = 1,01 X 10*' m ^ / s .

Remarque : Indiquons quelques valeurs du poids spcifique de la glace i 0* C (E 916,7 k g / m 3 ;

k g / m 3 ; - 20 C - w 920,3 k g / m 3 .

Module d* lasticit

Vf

a 10a C W " 918,6

16

3.2 - MODE DE CALCUL DE LA DISTRIBUTION DES DEBITS

a) Etablissement du systme d'quation

Les quations caractrisant la distribution des dbits dans un rseau
sont issues des deux principes suivants :

1 - Prin&tpe de onsPVatiort des dbits :

En un noeud quelconque, la somme des dbits qui arrivent au noeud est gale
la somme des dbits qui en partent.
De manire vidente :

q2 + ql = q3 + Q

ou

ql + q2 - q3 = Q

Si NNOEUD est le nombre de noeuds du rseau, on aura un systme de NNOEUD

quations (dites "quations aux noeuds") de la forme :

<+ 1) q (i) = Q (j)

(somme sur les conduites touchant le noeud j)

o q(i) reprsente le dbit dans la conduite i (le sens tant choisi de manire
arbitraire) et Q (j) le dbit extrieur au rseau au noeud j.
La somme des NNOEUD quations de ce systme fournit l'quation de conservation des dbits extrieurs au rseau :

NNOEUD
Q (j) - O
j - 1
Ces dbits extrieurs (dbits sources/rejets/consownation) tant supposs
connus, on retire du systme cette dernire quation.
On obtient alors un systme de (NNOEUD - 1) quations indpendantes faisant
intervenir les N inconnues q(i).
(N : Nombre de conduites du rseau)

17

2 - Prinoipe de ontinuit de la ohcge (ou "onsevvation de l'nergie?)

La chute de pression entre deux points tant indpendante du chemin
suivi, on crit que sur le contour ferm d'une maille quelconque, la somme
algbrique des pertes de charge est gale O ou la valeur de la surpression
(^ ) s'il existe un surpresseur dans la maille.

On appelle J (i) la perte de charge sur la conduite i. J(i)

est compt positivement si le sens du dbit dans i correspond au sens de parcoun

de la maille, ngativement sinon.

est compt positivement, si le surpresseur est dans le sens de la maille,

ngativement sinon.

Dans ce cas:

+ J(2) + J(3) - J(4) - J(5)

A chaque maille du rseau correspond donc une quation du type :

O
(+ 1) J(i)

OU

: valeur du surpresseur

i maille
,

La perte de charge est une fonction du dbit du type :

J = Rq

o R correspond la rsistance l'coulement de la conduite et o o( est

un rel voisin de 2 .
(voir chapitre "Pertes de charge").

18

On obtient un systme de N quations N inconnues en ajoutant

aux(NOEUD - 1) quations de noeuds, NMAILLE quations de mailles.
avec :
NMAILLE = N - (NNOEUD - 1)

NMAILLE est le nombre de mailles indpendante^ ;du rseau

(on choisira les mailles lmentaires).

On a donc le systme suivant :

(+ 1) R(i) q(i)

quations de mailles

= ^ (NMAILL) '

NMAILL
(+ 1) q(i) = Q (1)
quations aux noeuds

(+ 1) q (i) = Q (NNOEUD - 1)
(NNOEUD - 1)

S'il n'existe pas de surpresseur dans la maille j on crit que

Ad) = o.
N.B. Il s'agit d'un systme non linaire .

19

b) Exemple :
Etablissement du systme d'quations de ce rseau :

= 5

NNOEUD

NMAILL

N - (NNOEUD - 1)

2 mailles indpendantes

- quations de mailles :
(maille I)

* +

K_

(maille II)

- quations de noeuds :
Noeud 1
Noeud 2

q3 - q 5 =

d'o

le

Noeud 3

q4 = - Q3

- q2

systme non

linaire :

O -

-Rlql

+0

y,
+ 0

+ O

+ O

+o
+o
+

-q5
+ O

D'autre part l'quation : Q. +Q 2 = Q

+ Q. doit tre respecte.

20

On appellera matrice caractristique du rseau la matrice du systme.

M=

"R2
0

-1

-1

0
R

3 ~ 4
0
0

-1

Le vecteur (O, -&, Q., Q 2 - Q 3 ) constitue le vecteur des termes

indpendants (vect.)

On notera qu' ce stade, il est ncessaire de s'tre fix un sens

arbitraire pour chaque conduite, le sens effectif sera donn par le signe du
dbit dans le rsultat.

21

c) Rsolution du systme non-linaire

La mthode utilise pour rsoudre ce systme non linaire est une

mthode de calcul numrique itrative qui fait appel au procd de linarisation de la mthode de NEWTON :

soit rsoudre F (X) = O o F est une fonction non linaire.

On se donne Xo suffisamment proche de la solution suppose (pour

une fonction F donne, la convergence du procd itratif dpend de Xo) et
_ la prcision du rsultat.
On dtermine X. tel que :

F (Xx) fifi O = F (Xo) +

-j Xo) F(Xo)

(1)

(dveloppement au 1er ordre)

On pose X

= X4 et on reitre le calcul jusqu' ce que

X. - X

Pour le systme tudi, les quations de maille sont du type :

( i d . . M NMAILL))

ri
o les Mij sont les composantes de la matrice caractristique.

22

La relation

entre deux itrations q et q. est alors donne

par la formule itrative :

+

'J
Gnralisation : (o^Mij q.

?L (^j " ^J>( "iJ qj

^ ~

- O (2)

J sont les composantes du Jacobien J de la

fonction F .

On a

avec :

F (qo) + (qj - q ) J (q ) = O

(il

J (q) =
Fn

^ Fn

De (2) on tire la relation linaire :

/
(0-1) \
<

( Mij qo.

I ql

^Zs*M i 3 *

On est donc ramen chaque itration la rsolution d'un

systme linaire de la forme :

avec
et

Aij = Mij qo.

i d r NMAILL)

Aij = Mij

i ^ (NMAILL + 1, N)

Vi

L'algorithme est le suivant :

qo = distribution arbitraire

23

d) Amlioration de la convergence

II est ncessaire, pour dbuter le procd dcrit prcdemment, de

se fixer une distribution arbitraire des dbits (qui doit satisfaire les
quations de continuit aux noeuds). Il est possible de lever cette contrainte
et d'amliorer en mme temps la convergence en faisant la remarque suivante :

Les exemples de calcul de rseau maills effectus sur la base

d'un rgime laminaire, c'est--dire une loi de perte de charge linaire
(J = R .q) ont montr des rsultats relativement proche des rsultats en
rgime turbulent.

La distribution initiale go peut donc tre calcul par la rsolution

du systme linaire : (M) go = V traduisant ainsi un rgime laminaire (au
rapport prs des rsistances laminaire/turbulent).

24

3.3 - Mthode de calcul de la pizomtrie

a) Units et terminologie employes
Que ce soit au niveau de la saisie des donnes, au niveau
des calculs ou au niveau de l'dition des rsultats, toutes les pressions
caractristiques du rseau sont exprimes en hauteur de colonne d'eau
(mtre C E . )
2
avec 10 m CE = 1 kg/cm

1 bar.

(on ne tient pas compte ici, dans le domaine des pressions considr, de
la variation de masse volumique de l'eau avec la temprature).

On travaille alors sur les grandeurs suivantes :

- P

: pression effective au point considr (m.CE)

(pression absolue = P + P. atmosphrique)
Pour un point la pression atmosphrique, on aura donc : P = 0

- Z

: cote au sol (en m NGF)

(N.G.F. : nivellement gographique de la France)
correspond en fait l'altitude du point considr.

- H

: cote pizomtrique (m NGF)

- J

: perte de charge sur la conduite AB (m CE)

- /\

: surpression ventuelle (m CE)

Avec les relations suivantes :

H =Z+ P
J

AB

= H

A - HB + A

Voir figure 331 le graphe caractristique des pressions dans une conduite.
REMARQUE : Seule la cote pizomtrique H est significative du sens de
l'coulement de l'eau dans un rseau ou dans un milieu quelconque (Analogie
Electrique Potentiel/courant).

25

b) Calcul des cotes pizomtriques aux noeuds

Le calcul des cotes pizomtriques s'effectue partir de la

formule :

H. -

J -

(1)

o la perte de charge est donne par :

J =
Ce calcul intervient donc aprs le calcul de la distribution
des dbits.
La mthode utilise pour connatre la cote pizomtrique en
chaque noeud du rseau, consiste venir appliquer la relation (1) sur
toutes les conduites dont on connat dj la cote pizomtrique d'un
des deux noeuds extrmits. On boucle sur ce procd jusqu' la connaissance totale des cotes pizomtriques.

Il faut donc noter/ pour dmarrer le processus, la donne

ncessaire d'un noeud pression dtermine (pression de rfrence).

26

c) Calcul des_pressions effectives

Les cotes pizomtriques tant connues, les pressions effectives
en chaque noeud s'obtiennent simplement par la relation :

P = H - Z
o Z est la cote au sol du point considr.

27

Ceracte.rs/'iot/es
h-

Fiat/re.

d<*s

28

3.4 - CAS DES CONDUITES SINGULIERES

a) Dbit de consommation rparti sur une conduite ("Dbit en route")

Dans un rseau de distribution, certaines conduites prsentent de
nombreux points de soutirage rpartis plus ou moins uniformment et o les
dbits sont tous du mme ordre de grandeur. Il est alors possible de se
ramener au cas d'une conduite simple. (Tous ces points de soutirage devraient
tre logiquement considrs comme autant de noeuds relis par des tronons
simples) :

Pour les conduites de ce type on introduit la notion de densit linique

de dbit q ; si Q est le dbit total distribu de A B et L la longueur :

. Calcul de la perte de charge dans ce cas :

cl*

D : dbit entrant
P : dbit sortant

On a la relation : D = P + O
rsistance linique
On a : dH =-rdx,- dbit 2 (x)
=-rdx (D - qx)
0

=~rdx (D

- 2 qDx + q

x\

29

d'o :
2
HB - HA = - RD + DQR

R
r

HA - HB = R (D - DQ + -*
3

2
Q
)

avec D = P + Q '

HA - HB = R (P + PQ +
2

or

P22 + PQ + 2

2
(P2 + PQ + ) +

12
de mme :

P + PQ + ~ = (P + - )
3
\3
(P

]f

)2

- (2 - )3 ) PQ

73
- o,1547

On a donc :
(P + 0,5 Q ) 2 < (P2 + PQ + - | ) 2 <

(P + 0,577 Q ) 2

II apparat alors que, du point de vue des pertes de charge, on pourra

considrer la conduite comme parcourue par un dbit fictif Qf tel que :
P + 0,5 Q < Q f < P + 0,577 Q
2
On alors J = R Qf (cas d'une conduite simple).

On prend gnralement

I Qf = P + 0,55 Q

On aura compatibilit des dbits aux noeuds en rpartissant le dbit

"en route" O en A et B avec 0,45 Q pour A et 0,55 Q pour B:

" \ y

30

b) Cas des vannes et surpresseurs

- Surpresseurs
Le cas des surpresseurs est envisag de manire rigoureuse
puisque les surpressions ventuelles interviennent dans la mthode de
calcul des dbits et des pressions (voir 3.2 a) et 3.3).

Un surpresseur est caractris par la valeur de sa surpression

(valeur positive). Le sens du dbit dans la conduite concerne doit tre
pralablement connu.

Aprs calcul des dbits, on choisit le surpresseur dont la

courbe caractristique est compatible avec le point de fonctionnement
considr.

31

- Vannes
Dans une vanne, la chute de pression est de la forme :

o C est un coefficient li au degr d'ouverture de la vanne (voir

tableau ci-aprs).

Il est possible et commode de traiter les vannes comme des

surpresseurs valeur ngative. Le calcul des dbits permet alors de
dterminer le degr d'ouverture correspondant la valeur de la chute
de pression dfinie pralablement.

<o
EVOLUTION DU COEFFICIENT C

degr d 'ouverture- %

100

87.5

0.12 0.15

75

62.5

0.26 0.81

50

37.5

25

12.5

2.06 5.52

17

98

32

3.5 - Validit des rsultats

Un exemple de simulation complte de rseau est donn en
annexe 1 (rseau "Exemple"). Voici titre indicatif les rsultats
obtenus par la mthode de Hardy-Cross manuelle applique ce

mme

rseau (source : A. DUPONT "Hydraulique urbaine, exercices et projets")

120 l/s
l= 68,67m

Les rsultats concordent donc parfaitement tant du point de

vue des pressions que des dbits.

3.6 - Extension possible aux rseaux de vapeur

Cette tude peut s'tendre, par analogie, des rseaux de vapeur (ou de
gaz quelconque) moyennant une adaptation adquate de la modlisation des pertes de
charges et la prise en compte de la compressibilit (la convertion des dbits est
appliquer aux dbits massiques).

4 -, ETUDE THERMIQUE

33

4.1

- Modlisation des dperditions thermiques

a

^ PP^^iti2n_tliennique_clans une conduite enterre

On suppose le rgime thermique tabli.

5*/

/ / / / / / / / ss y /

A \ \ \\

\\\\\\\\\\\\\\

tec.

dr

: diamtre intrieur de la conduite

e.

: paisseur du tube

'is

: paisseur de l'isolant

: dbit volumique

: profondeur d'enfouissage

34
On considre un segment de longueur dx, l'nergie calorifique
perdue par le fluide est donne par :

V ~ - O (qdt) C dT = 0 dt

d'o le flux de chaleur :

(1)

D : masse volumique
C

: capacit massique

/^*

- Echange eau/paroi interne

Le flux est donn par :

(//b dx) (T - T)

(1)

o h est le coefficient de convection

2
eau/tube
(W/m C)

dans le tube e t dans

l'isolant

La loi de Fourrier (0 = f\ S grad T)

donne la relation :

Pour le tube :

0=2,

- pour l'isolant :

0=

(2)

(eis
_

, D + 2 eis
T
Ln (

>$> et)
Te)

)
(3)

avec

: coefficient de conduction du tube

: coefficient de conduction de l'isolant (W/m

C)

35

- Echange paroi externe/sol

On a la relation :

sol dx
(4)

o To est la temprature de surface du sol

X sol : coefficient de conductivit moyenne du sol
f est un coefficient fonction de la profondeur et du diamtre extrieur.

On a :

(facteur de forme)/.

36

b) Evolution de la temprature
Des relations prcdentes, on tire la loi de temprature

(T - To) dx
CqdT
(RI + R 2 + R 3 + R4)

avec :
=

Ln

D + 2 et

27\

Ln

D + 2 iso

A sol

'A:

27f/Aiso

Dans la plupart ds cas et dans le domaine des valeurs utilises

ici, les rsistances thermiques de convection et de conduction dans le tube
sont ngligeables devant les rsistances de conduction dans l'isolant et
dans le sol.

On pose R = R. + R,

On a alors :
d (T - To)
(T - To)

dx
C q R

37

Ln (T - To)

+ KO

PC q R

(T - To) =

Ki e

1L
Conditions aux limites :

T (o) = T.
T (L) =

d'o :

(T- - To)

- To)

ou encore :
+ B
en posant :
RP C
A = e
B = (1 - A) To

<

1
RPC

38

c) Dperditions thermiques, rcapitulatif

L
Pour chaque conduite

(i)

-y-

i , on a

=
(proprits phyqiques de l'eau) .
longueur (m)
dbit (m 3 /s)
calcul de f

DONNEESi
(i)

paisseur isolant ou tube

(i)

conductivit isolant ou tube

\ sol
To :

A(i)

temprature de surface

exp ( -

(1 f

e
To
. D + 2e ,
Ln (
)

sol

7T\

Calcul de R

39

4.2 - Mthode de calcul de la distribution des tempratures

a) Exemple
Considrons l'exemple du noeud suivant :

On pose : T.(i) : temprature de dbut de la conduite i

T (i) : temprature de fin de la conduite i

Pendant l'instant dt, la conservation de l'nergie calorifique au noeud
se traduit par l'quation :
T2(l) + q(3)

C (q(l) + q(3) + QJ

dt =

T dt

c'est--dire :
=

q d ) T 2 (D + q(3) T? (3)

Tjnj

+ q(3) + Q

- T

: temprature au noeud (temprature de mlange)

: temprature d'injection au noeud 1

(1)

40

T est aussi la temprature de dbut des conduites qui partent

du noeud.

T^B) = ^(7) = T

' (2)

De plus, on a vu prcdemment que la relation entre T 2 et Tj

est de la forme :

aTj + b

(3)

de (1), (2) et (3), on tire la relation :

q(3) a(3)
q(3)

b(3) q(3)
q(3)

b) Gnralisation
Pour chaque conduite i du rseau, on a donc la relation :

q(k) a(k)
q(k) + Q

b(k) q(k) + Q T i n j
q(k)

La sommation se faisant sur les conduites dbouchant sur le noeud origine

de i
Q et Tj_nj sont respectivement le dbit et la temprature d'injection
sur ce noeud.

On a donc affaire finalement un systme linaire de dimension N

(N : nombre de conduite du rseau) d'inconnue Tj, temprature de
dpart des conduites.

inj

41

Les composantes de

sont donnes par

q(k) a(k)

G (i,k) =
q(k) + Q

si k est une conduite dbouchant sur le noeud origine de i

G (i,k) - 0

sinon

les lments de W sont donns par :

b(k) q(k) + Q T i n j
q(k) + Q

La rsolution de systme fournit les tempratures de dbut de conduites,

donc uniquement les tempratures des noeuds correspondant au moins un
dpart de conduite.

On calcule les tempratures aux noeuds exclusivement d'extrmit

en appliquant la relation entre noeud dpart et noeud extrmit d'une
conduite :

T2

a Tj + b

On a alors la temprature en chaque noeud du rseau.

42

c) Validit des rsultats obtenus

Le calcul des tempratures se fait sur la base du rgime
thermique tabli.

En fait, le temps d'tablissement de ce rgime dans un sol est

suprieur celui des fluctuations de temprature en surface et ce rgime
thermique n'est donc jamais atteint. Compte tenu de cet tat de fait, il
n'est pas utile de chercher affiner au maximum la modlisation thermique.

On considrera donc de manire globale :

- une conductivit moyenne et paisseur moyenne d'isolant

(ou de tube),
- une conductivit moyenne du sol,
- une profondeur moyenne d'enfouissage.

Il n'est alors pas question de considrer les rsultats fournis

exacts moins de 0,1 prs, mais de les interprter comme significatifs
de la distribution des tempratures dans des conditions dfavorables
(faible temprature de surface) et rvlateurs des ventuels dfauts
thermiques du rseau.

5 - PROGRAMMATION

43

5.1 - Gnralits
La programmation s'est effectue essentiellement sur
MICRAL - 21 C dont les caractristiques sont les suivantes :

- microprocesseur
- mmoire vive

8 bits

Z 80

64 K . O .

- 2 units de disquettes 5" 1/4

600 K . O .

Systme d'exploitation PROLOGUE, langage FORTRAN 80 (microsoft)

Principaux utilitaires du systme PROLOGUE :

-

Mode edit :

--> EDIT,Nom de fichier

- Mise a jour (sous edit) :

ESC

; *E

- Compilation :
--> F80,Nom de fichier

- Link :
-->L80
* Fichier source,fichiers compils associs + bibliothque
fortran (FORLIB/S),fichier objet/N/E ou G(excution)
- Listing des fichiers (duplication) :
--> CP,DF,Nom de fichier-S,LO (imprimante) .

44
a) Utilisation des fichiers disques
L'utilisation des fichiers permet, pendant l'excution d'un
programme, de venir stocker ou saisir des donnes numriques ou littraires
sur un support magntique de manire squentielle ou de manire slective
(accs direct). Un fichier est caractris par son nom.

- Accs squentiel
Les donnes sont stockes de manire linaire au fur et mesure
de leur criture. Les formats de lecture doivent tre compatibles avec les
formats d'criture, ("tJ

Les fichiers A.S. sont analogues aux fichiers sources, il est

possible de les crer et de les modifier sous le mode EDIT.

L'criture sur un fichier dj existant provoque son effacement.

L'utilisation de ces fichiers se fait de la manire suivante :

- ouverture du fichier nomm

- lecture ou criture des donnes
- fermeture du fichier

La srie des instructions utilises est la suivante :

CALL OPEN (num, NOM-S , O OU 1)

READ ou WRITE (num, 1)
1

FORMAT
ENDFILE

Num

liste

num

: numro logique affect au fichier dans le programme.

O ou 1 : unit de disquette support.

(1) Remarque :
L'intrt de stocker les donnes sous forme formate rside
dans le fait de permettre le contrle, l'dition, voire la modification
de ces donnes.

Si ce contrle ne se justifie pas, on prfrera une criture/

lecture sur fichier des donnes non interprtes (sans format), l'opration
se faisant alors beaucoup plus rapidement.

45

- Accs direct
Le fichier A.D est divis en enregistrements. L'intrt de ce
fichier est qu'il permet l'opration d'entre/sortie un endroit prcis
du fichier correspondant un enregistrement. Chaque enregistrement occupe
une place de longueur prdfinie (ici, chaque enregistrement occupe un
secteur du disque), ces fichiers prsentent donc l'inconvnient d'utiliser
une place plus importante, pour une mme masse d'informations que les
fichiers A.S.

Un fichier existant peut tre prolong par cration de nouveaux

enregistrements.

L'utilisation de ces fichiers se fait de la manire suivante :

- ouverture du fichier nomm

- lecture ou criture sur les enregistrements nomms
- fermeture du fichier.

La srie des instructions est la suivante :

CALL OPEN D (num/ NOM - D, O ou 1, Nrec)

READ (num, f, REC = ni) liste
WRITE (num, f, REC = n2) liste
CALL SVREC (num, n2)
f

FORMAT
CALL CLOSED (num)

Nrec

ni, n2 :

numro du dernier enregistrement effectu

numros d'enregistrement.

OPEN D, CLOSED, SVREC sont des sous-programmes appartenant au

programme de gestion des fichiers accs direct propre au systme
d'exploitation.

L'effacement ou la suppression d'un fichier A.D. se fait par

lancement du programme-systme : CHPURGC.

tt'vHH

150

222

satriow d'vo fkhie1 acc?

se.qventi'e

WRITE(1,150)
F O R M A T I / / . 3 X , ' M i s e en f i c h i e r . ' )
CALL 0 P E N I 4 , ' M A T R I C E - S ' , 1 )
DO 222 1=1,NBRAN

WRITEI4)M(I,J),J=1,NBRAN)
CONTINUE
ENDFILE 4

C
161

CALL 0PEN5,'RESIST-S ,11

WRITE.5,161)R(I),1=1,NBRAN)
FORMAT5F15.2,/)
ENDFILE 5

Exemple d'utilisation d'un fichier accs direct

CALL 0PEND<3,'D0NNIES-D ' , 1 , N R E C 3 >

READ(3^t5UREC=5) NMAX
READ(3,451,REC=) NBRAN
READ(3,/t51,REC=7) NNOEUD
READ<3r451,REC=8) NMAILL

'97
420

C
999

RAD(1r501) MAILL(J)
K1=20+2*NMAX+I
WRITE(3,420,REC=K1) IrNCMA,<MAILL(J),J=1,NCMA)
FOR MAT <1X,13,IX,13,IX,30(14))
CALL SVREC<3,K1>

CONTINUE
CALL CL0SEDC3)

Mar? P </ le fron

/,'

Ficli fe
Sovrct

A5-

oi~i p

ton

A-o

- 0

c -5

C /?

D -
-S

D - fi

-5
dit

1
1

-O

lrht,rc w

48

c) Structure du logiciel
Le logiciel dcrit ici a une structure "chane", il est compos
de plusieurs programmes indpendants lisant ou crivant sur des fichiers
communs (accs squentiel ou direct) et utilisant des sous-programmes
communs.

Outre la souplesse d'emploi du logiciel, une telle structure

apporte une certaine aisance dans la conception et dans la manipulation
des fichiers-source (possibilits de modifications spares, temps de lecture
sur disque et de compilation rduits).

De plus, l'ensemble peut tre facilement complt par l'apport

de nouveaux modules.

Cela permet enfin une optimisation de la place mmoire vive

occupe lors de l'excution.

L AUKi-

50

e) Description des variables communes tous_les programmes

ALPHA

: Exposant du dbit dans la formule des pertes de charge

ERREUR

: Pourcentage d'erreur dans le calcul des dbits

NBRAN

: Nombre de conduites du rseau

NMAILL

: Nombre de mailles

NMAX

: Nombre maximum de conduites

NNOEUD

: Nombre de noeuds du rseau

NREF

: Numro du noeud de rfrence

PREF

: Pression de rfrence

TABLEAUX :
COTE (i)

Cotes au sol des noeuds

DEB (i)

Dbits aux noeuds

D (i)

Diamtre des conduites

K (i)

Rugosit

L (i)

Longueur

MAILL (i,j

Composition de la maille i

NCMA (i)

Nombre de conduite par maille

01 (i)

Noeud origine des conduites

02 (i)

Noeud extrmit des conduites

DELTA (i)

Valc;ur des surpressions

DELTAM (i)

Surpressions dans les mailles

R (i)

Rsistances hydrauliques

Q (i)

Dbit dans la conduite i

M (i, j)

Matrice hydraulique du rseau

VECT (i)

Vecteur des termes indpendants

P (i)

Pression effective au noeud i

H (i)

Cote pizomtrique au noeud i

TI (i)

Temprature d'injection

T (i)

Temprature au noeud

51

5.2. - Saisie et stockage des donnes

PROGRAMME SAISIE
Le but de ce programme est d'tablir le fichier principal
des donnes caractristiques d'un rseau quand celui-ci correspond
une nouvelle tude (nouveau rseau). Le titre de ce fichier correspondra
au titre de l'tude.

Ce fichier est accs direct et a pour nom : DONNEES - D

L'organisation de ce fichier est dcrite plus loin.

Avant cette mise en fichier, le programme "SAISIE" effectue en

mode conversationnel la saisie de toutes les caractristiques gomtriques
et numriques du rseau ncessaire l'ensemble de l'tude.

C'est donc dans le cas d'un nouveau rseau, le premier programme

lancer dans la chane constituant le logiciel.

Particularits

La mise en fichier se fait par appel un sous-programme

extrieur nomm RECDAT.

Les variables communes SAISIE et RECDAT sont dclares

par 1'instruction COMMON.

La saisie est divise en rubriques et se fait sous forme de

questions ou sous forme de tableau pour les donnes rptitives (conduites,
mailles, noeuds). Chaque rubrique est suivie d'une option correction.

Description des variables spcifiques

I, J : test de modification
N, KO, KOND, KONDAB: variables auxiliaires de calcul
NSurp : nombre de surpresseur
REP, NON, OUI : variables caractres de test
ESC, CH, CHO, CHI, CH2 ... : caractres de contrle de l'imprimante
ou de l'cran.

52

Remarque sur la saisie des mailles

Une maille est caractrise de la manire suivante :

NCMA (i)
MAILL (i,j)

: nombre de conduite de la maille

: numro de la jme conduite de la maille i
affect du signe - si le sens de la conduite
est l'inverse du sens de la maille

53
Algorithme du programme SAISIE

Dclaration des variables

Initiali sation

Saisie du titre

et du

nombre de conduites

Saisie des donnes

caractristiques des
conduites

Sortie sur imprimante

Corrections

Saisie du nombre de nombre de noeuds

Calcul du nombre de mailles

Saisie de la composition
des mailles

Corrections

nu
Saisie des dbits
aux noeuds

Corrections

Saisie des cotes au sol

I Surpresseurs

Donnes numriques
complmentaires

1
Appel de RECDAT
(mise en fichier)

I,
STOP

54

SOUS-PROGRAMME RECDAT

Ce sous-programme, appel par SAISIE, cre le fichier accs

direct DONNEES - D et crit les donnes saisies aprs les avoir structures
en bloc, chaque bloc est affect d'un numro d'enregistrement.

L'algorithme est le suivant

Dclaration des
variables

Ouverture du fichier DONNEES-D

Affectation d'un numro

d ' enregistrement

Ecriture du bloc de
donnes correspondant

Fermeture du fichier

STOP

55

ORGANISATION DU FICHIER DONNEES-D

REC

DONNEES

Format d'ecri turo

1 1

1S

TITRE

39X,15A1

13

NSIM

39X.15

14

1 5

1 NMAX

39X,15

1 6

1 NBRAN

39X,I5

17

1 NNOEUD

39X,I5

1 8

1 NMAILL

39X.I5

19

1 ERREUR

39X.F5.2

1 10

1 ALPHA

39X,F5.3

1 11

1 NOREF

39X,I5

1 12

1 PREF

39X,F6.2

13

20
21

Cote(l)

20+NNOEUD

F7.2

Cote(NNOEUD)

F7.2

20+NMAX
20+NMAX+l

1,01(1),0211),D( 1),K(1),

Lll),DELTA! 1)

3(1X,I3),1X,F5.O,1X,F6.3,
F9.2,5X,F9.2

20+NMAX+NBRAN

20+2NMAX
20+2NMAX+1

1,NCMA(1),MAILLI1,1),

IX,13,IX,13,IX,30(14)

DELTAM1)

Fil.5

I20+2NMAX+NMAILL

20+3NMAX
20+3NMAX+1

I20+3NMAX+NMAILL

I21+3NMAX+NMAILL
I .
I .
I20+3NMAX+NMAILL
I +NNOEUD

DEB(1 )

Fil.5

56

5 . 3 ~ Modifications des donnes

PROGRAMME MODIF

Le but de ce programme est de modifier une ou plusieurs

valeurs des donnes prsentes sur le fichier DONNEES-D. Ce programme
est utilis entre chaque simulation d'une mme tude. Il incremente
le numro de simulation.

Algorithme :

Ouverture de DONNEES-D

to

Affichage du Menu des

modifications

Recherche de l'enregistrement
correspondant au choix

Saisie de la modification
Ecriture de la modification

Fermeture du fichier

57

5. i - Etablissement de la matrice caractristique du rseau

PROGRAMME MATRIX

Le but de ce programme est d'tablir la matrice hydraulique

du rseau telle qu'elle a t dfinie dans le paragraphe 3.2. Il tablit
galement le tableau des rsistances dont ie calcul est dcrit au paragraphe
3.1. Il cre les fichiers MATRICE-S et RESIST-S.

En principe, il est utilis chaque simulation aprs SAISIE

ou MODIF.

58
ORGANIGRAMME SIMPLIFIE

Initialisation

////?/~/?7X J

Ouverture de DONNEES-D
Lecture de NMAX, NBRAN,
NNOEUD, NMAILL, ALPHA
i = 1, NBRAN J .

Lecture des caractristiques

de la conduite n i :
01, O2, D, K, L

Calcul de la rsistance
R (i)
M (NMAILL + Ol, i) = 1
M (NMAILL + O2, i) =-1

= -, NMAILL")

Lecture des caractristiques

de la maille N i :
NCMA, MAILL
j = 1, NBRAN

OUI

,j) = signe (Maille(NN)) x

Fermeture de DONNEES-D

Mise en fichier deJMJ

dans MATRICE-S
Mise en fichier de R
dans RESIST-S
STOP

59

5.5 - Calcul de la distribution des dbits

PROGRAMME DEBIX

Ce proframme effectue le calcul itratif sous forme

matricielle des dbits tel qu'il est dcrit au paragraphe 3.2.

L'algorithme est.le suivant :

lecture sur DONNEES-D des variables concernes

lecture sur MATRICE-S de

calcul des dbits

mise en fichier (DEBITS-S)
dition du rsultat

Particularits : l'inversion et le produit des matrices se fait par

appel des sous-programmes "extrieurs contenus dans BIBLIO. Les variables
communes sont dclares sous l'instruction COMMON.

60

ORGANIGRAMME SIMPLIFIE
DU CALCUL DES DEBITS

Lecture de M, VECT
sur fichiers

Appel de INV
M = INV (M)

Appel de PROD
Q = M.VECT

, NBRAN )

s"o\

<

0UI

Signe(j) = - signe(j)

Lecture de M
1, NBRAN J

^T T

1, N B R A N ^

M(i,j) = signe(j)* M(i, j)

= 1, NMAILL )

O<

j -

NBRAN )

\
M(

M(i,j

Qo = Q

( C<- 1)

61

5 . 6 ~ Calcul de la distribution des pressions

PROGRAMME PIEZO

Ce programme effectue le calcul des cotes pizomtriques et

des pressions effectives aux noeuds tel qu'il est dcrit au paragraphe 3.3.

Algorithme :

- Lecture sur DONNEES-D des variables concernes

- Lecture des rsistances sur RESIST-S
- Lecture des dbits sur DEBIT-S
- Calcul des cotes pizomtriques et des pressions effectives
- Mise en fichier (CHARGES-S)
- Edition du rsultat

ORGANIGRAMME SIMPLIFIE DU CALCUL DES PRESSIONS

Lecture de 01, 02, A

, R, q, cote/ sur fichiers

Test (01(j)) = 1

63

5 . 7 - Calcul de la distribution des tempratures

PROGRAMME THERMO

Ce programme effectue le calcul des tempratures tel qu'il

est dcrit au paragraphe 4.2 .

Algorithme :

- Saisie des donnes thermiques

- Lecture sur DONNEES-D des variables concernes
- Lecture des dbits dur DEBITS-S
- Calcul des tempratures
- Mise en fichier (THERM-S)
- Edition du rsultat.

Ce programme fait appel, pour la rsolution du systme

FG"J

T"I = W, aux sous-programmes INV et PROD contenus dans BIBLIO-R.

64

5 . 8 - Inversion d'une matrice, produit d'une matrice par un vecteur

| Sous-programme INV et PROD

Le sous-programme d'inversion utilis reprend la mthode du

Pivot (mthode de GAUSS).

Cette mthode a pour avantage de rduire la taille mmoire

utilise, la matrice initiale est dtruite et remplace par son inverse.
Cette mthode tant numrique, la prcision est excellente mais le temps
d'excution sur micro-ordinateur peut s'avrer trs long dans le cas de
grands rseaux.

Pour les trs grands rseaux, on aura tout intrt chercher

remplacer le sous-programme d'inversion utilis ici par un autre
utilisant une mthode numrique (mthode de GAUSS-SEIDEL).

65

B I B L I O G R A P H I E

CARNAHAN B.
Applied numerical methods.

DUPONT A.
Hydraulique urbaine.

LANCASTRE A.
Manuel d'hydraulique gnrale.

Technique de l'eau (La)/ - n 438/439

VILLEGAS E.
Etude de calcul automatique des rseaux de distribution
d'eaux (thse).

66

- CONCLUSIONS -

L'application de ce logiciel sur un exemple prcis a montr outre des

rsultats trs satisfaisants, une convergence rapide du procd (nombre d'itrations
rduit) et une grande souplesse d'utilisation issue de la structure "chaine" de
ce logiciel. De plus, l'existence de nombreux contrlesde compatibilit des donnes
et la possibilit d'dition de tous les fichiers de travail apportent une certaine
fiabilit d'utilisation.

L'inconvnient majeur reste actuellement, pour de grands rseaux, et pour

une

utilisation sur micro-ordinateur, une dure trop longue pour l'inversion de

matrice due la mthode d'inversion utilise. Il est possible de lever cet inconvnient de 3 manires diffrentes :

- soit utiliser un sous-programme d'inversion plus rapide (mthode numrique)

- soit effectuer l'opration d'inversion sur une unit centrale de classe suprieure
(le SGR Aquitaine est reli au centre de calcul du BRGM d'Orlans (VAX et IBM)
par le rseau TRANSPAC)
- soit si cela est possible, scin'der le rseau en plusieurs sous-rseaux.

La programmation de ce logiciel est adapte son utilisation sur

MICRAL 21-C. Son extension d'autre micro-ordinateur ne doit gure poser de problmes, la principale modification rsidant sans doute dans les commandes d'accs
fichiers qui sont ici des sous-programmespropres au systme. L'intrt d'un tel
logiciel, et notamment dans le cas de DAX o est envisag le recyclage des nombreux
rejets du rseau existant, est de permettre le dimensionnement des conduites de bouclage, des pompes, vannes, circulateurs, nouvelles conduites, e t c . . apportant ainsi
une aide importante la dcision au niveau de l'tude du projet. Son exploitation
permettra galement une optimisation du fonctionnement par des rglages ajusts au
fonctionnement dsir.

Tout cela se traduisant en fait par des conomies au niveau de la ressource en eau, des frais d'exploitation et des investissements.

Le logiciel LAUREL devrait donc se rvler un outil de travail important

pour l'tude des rseaux maills.

MINISTRE DE L'INDUSTRIE ET DE LA RECHERCHE

BUREAU DE RECHERCHES GOLOGIQUES ET MINIRES
SERVICE GOLOGIQUE NATIONAL
B.P. 6009 - 45060 Orlans Cedex - Tl.: (38) 63.80.01

ANNEXE 1
Logiciel LAUREL
Notice d'utilisation
et exemple
I 84 AGI 1<54 AQI \

Pessac le 29 juin 1984

MODELISATION DU COMPORTEMENT HYDRAULIQUE ET

THERMIQUE DES RESEAUX MAILLES

SERVICE GOLOGIQUE RGIONAL AQUITAINE

Avenue Docteur Albert Schweitzer - 33600 PESSAC - Tl.: (56) 80.69.00

Logiciel

L A U R E L

S O M M A I R E

1 - BUT
2 - PRINCIPE

1
1.

3 - DESCRIPTION
- Description

- Composition

4 - DESCRIPTION D'UN RESEAU,TERMINOLOGIE

5 - DESCRIPTION DES DONNEES

6 - EXPLOITATION
- Ut i 1 i sat i on sur MICRAL
- Droulement de l'exploitation
- Ajustement du modle hydraulique
- Dbit "en route"
- Modification du rseau
- Condition d'excution des programmes
7 - EXEMPLE D'UTILISATION
- Rseau "exemple"
- Droulement des programmes et rsultats

6
6
8
8
9
10
11
13

NOTICE D'UTILISATION

oo
oo
oo
oo
L A U R E L
oo
oo
__
_ _ _ __ __ __ _ _=_ _. _ _______

1 - BUT

Ce logiciel est destin a simuler les comportements

hydraulique et thermique en rgime tabli des rseaux maills et/ou
ramifis de conduites vhiculant de l'eau (chaude ou non) sous pression.
Aspect hydraulique
Sont concerns tous les rseaux de distribution correspondants la description prcdente.
Aspect thermique
Sont concerns les rseaux de chaleur monotubes composs
de conduites enterres vhiculant de l'eau chaude sous pression.

Il peut tre utilis lors de la conception d'un rseau ou sur

un rseau existant que l'on souhaite optimiser.

Une simulation complte fournit l'tat des dbits,pressions et

tempratures partir du scnario de fonctionnement qui est appliqu cxu
rseau.Pour un rseau donn,on tabli ce scnario partir de la modification des donnes d'ordre gomtriques,topographiques,hydrauliques
ou thermiques.

2 - PRINCIPE

Tous les rsultats sont issus du calcul des dbits qui constitue
la partie essentielle de 1'exploitt ion.Ce calcul,effectu sous forme
matricielle,uti1ise une mthode de calcul numrique itrative base sur
la mthode de Newton (Linarisation d'un systme).
La matrice caractristique du rseau est issue des principes de
continuit de charge et de conservation des dbits.Cette matrice et la
distribution initiale des dbits sont gnrs automatiquement ce qui re
duit le travail prliminaire de l'utilisateur au minimum.
La convergence du procd est gnralement rapide mais chaque
itration ncessite l'inversion d'une matrice qui,si elle effectue
sur micro-ordinateur,prendra un temps relativement long dans le cas
d'un grand rseau.( 11 est possible de remplacer le sous-programme
d'inversion utilis par un autre plus rapide!mthode numriqueJ).

3 - DESCRIPTION

Le logiciel LAUREL se prsente sous la forme d'une chaine de

programmes et de sous-programmes indpendants exploitant des fichiers
de donnes communs.

Composition :

- Programmes (fichiers excutables)

-

SAISIE
MODIF
MATRIX
DEBIX
P1EZO
THERMO

- Sous-programms (fichiers non excutables)

- RECDAT
- BIBLIO
- MENUX
- Fichiers de donnes (non excutables)
DONNEES
MATRICE
RESIST
DEBITS
CHARGES
THERM

Accs direct)
Accs squentiel)
M

II

4 - DESCRIPTION D'UN RESEAU"TERMINOLOGIE

On adopte la terminologie suivante :

Un NOEUD est,soit un point de jonction de plusieurs conduites soit
un point extrmit d'une ramification.
conduites,soit
de

dbit

C0NDUITE

(une

r m a n t u " e b o u c l e f e r m constitue une MAILLE du rseau

ne comportera pas de vannes fermes).

Il faut noter l'existence possible de vannes et de surpresseurs.

l M

.9
n

5 - DESCRIPTION DES DONNEES

L'tude ncessite :
- La numrotation des conduites
- La numrotation des noeuds
- La numrotation des mailles
- Un sens arbitraire pour chaque dbit
- Le sens arbitraire de contour des mailles
Les donnes ncessaires sont les suivantes :
- Le nombre de conduites maximum
- Le nombre de conduites de l'tude
- Le nombre de noeuds
-

Les noeuds dpart et arrive pour chaque conduite

Le diamtre de chaque conduite
N
La rugosit "
"
(coefficient k)
La longueur

- Le nombre de conduites dans chaque maille

- Le numro des
"
H
M
M
- Le dbit sur chaque noeud
- La cote au sol de chaque noeud
- La donne des surpresseurs ventuels
- La donne des vannes ventuelles
- Une pression effective de rfrence
- Le numro du noeud de rfrence
- La valeur de l'exposant du dbit dans la formule de perte de charge
- Le pourcentage d'erreur dsir dans le calcul des dbits
Remarque
Le programme de calcul de la distribution des dbits travaille
sur la base des dbits aux noeuds imposs et non sur celle des pressions
imposes.Ces pressions sont dtermines par la suite et font partie des
rsultats

================

6 - EXPLOITATION

Utilisation sur MICRAL

Le lancement sur MICRAL se fait de la manire suivante :

- Disquette programmes (LAUREL) dans FLO
- Disquette fichiers dans FL1
{cette disquette doit dj contenir le fichier CATADIR-D)
- Taper RETURN

Edition des fichiers sur imprimante :

- Taper CP,DF,Nom du fichier -S ou -D,LO

Effacement de DONNEES-D (rinitialisation)

- Taper CHPURGC

Droulement de l'exploitation

Une tude est caractrise par son nom (nom du rseau) et

chaque
simulation porte un numro.L'excution de MODIF
incremente le numro de simulation.

L'exploitation s'effectue en plusieurs phases diffrentes

suivant les cas.

Saisie des donnes caractristiques du rseau lors d'une nouvelle

tude ( nouveau rseau ) : excution de SAISIE :
. Le fichier DONNEES est vide
. Saisie des donnes au clavier
. Mise en fichier sur DONNEES
Modification d'un rseau donn en vue d'une nouvelle simulation :
excution de MOD1F
. Saisie des modifications au clavier
. Modification du fichier DONNEES

Simulation du rseau prsent sur DONNEES,1'uti1isateur a les

possibilits suivantes :
Excution de MATRIX :
. Lecture sur DONNEES
. Constitution de la matrice du rseau
. Ecriture sur RESIST et sur MATRICE
Excution de DEBIX
.
.
.
.

Lecture sur DONNEES,MATRICE

Calcul des debits
Ecriture sur DEBITS
Edition

Excution de P1EZO
.
.
.
.

Lecture sur DONNEES,RESIST,DEBITS

Calcul des charges
Ecriture sur CHARGES
Edition

Excution de THERMO
.
.
.
.
.

Saisie des donnes thermiques

Lecture sur DONNEES,DEBITS
Calcul des tempratures
Ecriture sur THERM
Edition

- Ajustement du modle hydraulique

Les pertes de charges sont modelisees par une relation du type :

Alpha
0 = R . Q

ou R est fonction,entre autres,du coefficient de rugosit k (formule

de Colebrook).
L'ajustement du modle hydraulique se fera donc :
- Soit
sur
la modulation de Alpha
(Alpha infrieur
ou gal 2 pour des dbits infrieur au m3/s ).
Ceci constitue une adaptation globale.
- Soit par la modulation de k sur chaque conduite,la valeur
de k tant en gnral difficilement apprciable
(adaptation spcifique).

- Debit "en route" :

Une conduite distribuant un debit Q rpartit sur

toute sa longueur peut tre traite comme une seule conduite
simple en rpartissant le debit Q chacune de ses extrmits,
avec :
- 0.45 Q pour le noeud dpart
- 0.55 Q pour le noeud arriv .
Les debits prlevs doivent tre d'un ordre de grandeur infrieur
au dbit de la conduite principale.

- Modification du rseau
Menu des modifications (programme MODIFJ
.Nombre de conduites
.Nombre de noeuds
.Nombre de mailles
. Prec i s i on
.Alpha
.Noeud de reference
.Charge de reference
.Cotes au sol
Caractristiques des conduites
.Caractristiques des mailles.
.Dbits aux noeuds
.Surpresseurs et vannes

6
7
8
9
10
11
12
13
14
15
16
17

L'ensemble des donnes devant tre cohrent,il est ncessaire

de modifier toutes les donnes concernes par une modification
du rseau
entre autres :
. Les noeuds suplmentaires entranent la donnes de leur debit
extrieur et de leur cote au sol.
. Les conduites suplmentaires peuvent entraner la cration ou
la modification de mailles.
etc. . .

# Suppression de conduites:

Pour des raisons de compatibilit de numrotation d'enregistrement s

fichiers,la suppression d'une conduite se fait par sa mise en "condu
morte" (ce qui pourrait se traduire en fait par la fermeture d'une v
ne ).
Exemple : suppression de 3; :

II est important de noter qu'ainsi,la suppression d'une conde

n'entrane pas la diminution du nombre de conduites.

- Condition d'excution des programmes aprs passage par MODIF

(nouvelle simulation)
* Les modifications numrotes
9,10,11,12,13,10
ncessitent pas l'excution de MATRIX avant DEBIX .

et

17

ne

#
Le calcul des pressions et des tempratures tant issue du
calcul des dbits#PIEZO et THERMO ne s'excutent pas si DEBIX ne
l'a pas t auparavant(apparition d'un message d'erreur)
(Remarque : Au cours d'une mme simulation,le programme THERMO
peut tre r-xcut avec de nouvelles donnes thermiques.)

10

EXEMPLE

EXEMPLE D'UTILISATION

On considere le reseau

"Exemple" de la figure page suivante

- 11 conduites
- 10 noeuds

Les donnes suplementaires sont les suivantes :

Noeud
Debi ts

1
00 1

80 1

100 1

50 1

90 1

120 1

80 1

0 1

45

44. 5

46 .2

45. 3

44. 75

45

50

9
0 1

10
580

Cote
44. 6
sol

- Noeud de reference : 10
- Pression de reference : 0 (p. atmosphrique)
On fixe

alpha

= 1.98

et l'erreur

e=0.1% .

50

105

feSC U

Sens Je

Txe-HsijA It

es

-// / (fiZoo

Coe. Fft'ci'en t
0

cle roo si/r ' :

: iyt\

DEROULEMENT DES PROGRAMMES

CDCDCDCDCDCDCDCDCDCDCDCDCDCDCDCDCDCDCDCDCDCDCDCDCDCDCDCDCDCDCDCDCDCDCD
CD
CD
SIMULATION DU COMPORTEMENT HYDRAULIQUE ET THERMIQUE
CD
CD
DES RESEAUX MAILLES
CD
CD
CD
CD
Los i c i el
CD
CD
CD
CD
L A U R E L
CD
CD
CD
CD
CD

CD
CD
CD
CD
CD
CD
CD
13
CD
CD

CD
CD
CD
CD
CD
CD
CD
CD
CD

Ph.DAUGUET
JC.MARTIN

BORDEAUX, Juin 1984

Bureau de Recherches Gologiques et Minieres r SGR/AQI

Mis en bouteille a la proprit?

CDCDCDCDCDCDCDCDCDCDLDCDCDCDCDCDCDLDCDCDCDCDCDCDCDCDCDCDCDCDCDCDCDCDCD

Menu principal :

L A U R E L
Donnees-D

SAISIE
MODIF

- >

MATRIX

DEBIX
- - >

/
/
/ !

Matrice !\
Resist

-> SAISIE

! Deb its !\
\

!
i

>

.V.
THERMO

\!

> Therm

!\

\.
- ->

PIEZO

>! Charges \
\.

PROGRAMME \" SAISIE""^

* Titre de l'tude [EXEMPLE |

Quel est le nombre maximum suppose de conduites du reseau ?

Nombre de conduites de cette etude

? 11

COMPOSITION DES CONDUITES

!
}

Conduite

Noeud or i s i ne

Noeud extrmit

Cond. 1
8

Cond. 2

Cond. 3

DONNEES DIMENSIONNELLES DES CONDUITES

Diamtre * Rugosit * Longueur

< mm > * < mm )
* < m )
*
Cond

1 * <

3 -

300

* Cond

2 * <

* Cond

3 * <

2 -

200

350

if

1 )

* Ecriture des donnes sur imprimante ? (0/N) : 0

Nombre de conduites :11

Conduite

! N.orisine-N.extreia.

3-8

!
i

Diamtre ( BO )
300.

Rusosite k ( n )
!

2.00

!
!

Lon3ueur ( )

200.00

8-2

350.

2.00

330.00

2 - 1

400.

2.00

400.00

8-4

350.

2.00

220.00

- 5

400.

2.00

400.00

1- ?

400.

2.00

290.00

5- 9

500.

2.00

400.00

5-6

300.

2.00

400.00

7-6

300.

2.00

330.00

10

9-7

400.

2.00

280.00

11

9- 10

800.

2.00

2500.00

Y a t il des corrections a apporter ? (0/N) : N

NOMBRE DE NOEUDS ? 10
II doit y avoir
PAUSE <CR>

mailles

indpendantes

COMPOSITION DES MAILLES

La maille

est forme de combien de conduites ?

Quelles sont ces conduites ?

(Avec + : sens de parcours
- s sens inverse >.
3

-7

La maille

est forme de combien de conduites ?

Y a t 11 des corrections a apporter dans

la composition des mailles ? <.0/N> :|N

Entres les debits sur chaiue noeud: (en m3/s)

< Avec + : debit entrant
- : debit sortant )

Noeud

Debit

0.0

etc.

Le debit sur le no 10 vaut s

-.580

Correction sur les debits ? (0/N) :

Saisie des cotes au sol des noeuds

( en m N.G.F. )

Noeud

Cote(m)

44.6

iuel est le noeud de reference ?

inoeud a pression fixe)

10

ression a ce noeud < en m . C E ) : 0

O
Combien y a f i l

. SURPRESSEURS
de surpresseurs dans ce reseau ? :
VANNES

'"Combien y a t'il de vannes dans ce reseau ? sj 0 I

Precision dsire dans le calcul des debits ? (en %) : 0.1

ALPHA ? s 1.98

EXEMPLE

- 0.

STOCKAGE DES DONNEES DANS LE FICHIER "DONNEES" (Accs direct) .

******K**K#W****tttttt**#*********X**********#**K***************

Veuillez patienter

STOCKAGE TERMINE

STOP

.is fry.!
civ A c A/.!

TITRE DE L'ETUDE
NUMERO DE LA SIMULATION

RESEAU EXEMPLE
O

NOMfIRE MAXIMUM DE CONDUITES

NOMBRE DE CONDUITES
NOMBRE DE NOEUDS
NOMBRE DE MAILLES
ERREUR
ALPHA
NOEUD DE REFERENCE
PRESSION DE REFERENCE

15
11
10
2
.10
1.980
10
0.00

COTES AUX NOEUDS (ni.NGF) s

45.00
44.50
46.20
45.30
44.75
45.00
48.00
48.00
105.00

CONDUITES(nuniero 3
8 300.
1
2 350,
2
8
1 400.
2
3
4 350.
4
8
4
5 400.
5
1
9 400.
6
9 500.
5
7
6 300.
8
5
300.
9
7
9
7 400.
10
11
9 10 800.

N.org i ne - N.fi nal - D. - K - L. - Sur p. (niCE)

2.000
200.00
0.00
2.000
330.00
0.00
2.000
400.00
0.00
2.000
0.00
220.00
0.00
2.000
400.00
2.000
290.00
0.00
0.00
2.000
400.00
2.000
400.00
0.00
0.00
2.000
330.00
2.000
280.00
0.00
2.000 2500.00
0.00

MAILLES (numero - ni)rc. conduites- num. conduite:

1 6
2
3
6 -5 -7 -4
2
4
10
9 -8
7

VECTEUR DES TERMES INDEPENDANTS (Surp.- Debits)

0.00000
0.00000
.00000
.08000
.10000
.05000
.09000
.12000
.08000
0.00000
0.00000
-.58000

Retour au menu principal

PROGRAMME " MATRIX

Lecture sur "DONNEES-D" en accs direct

M ise en f ich ier.

Retour au menu principal

Programme DBX

Lecture sur fichiers en cours

Calcul des debits

Numeros des conduites dont vous souhaiter voir

l'volution du debit au cours du calcul <4 numros)
1
3
7
8

! It er.

! Con cl.

! Cond.

! Cond.

! Cond.

! Ecart, max

: Con

.1000 !

.1094

.2671

-.0565

100.00 7. : 1

1
1

.1000 !

.1156

.2568 !

.0524

206.37

.1000 !

.1156

.2568

.0524

* M i se en f i ch i er

* Edition des rsultats

.06 7. :

ESE AU EXEMPLE!

o
G A l_ C

Prec i s i on

O CE S

O E B I T

1 0 'A

RESULTATS s
Noabre d'itrations : 3

Conduite

Debit ( B3/S )

Vitesse (m/s)

Remarque

1
1

.1000

1.41

I
1

.0356

.37

1
1

.1156

.92

.0644

.67

.1144

.91

1
1

.1756

1.40

1
1

.2568

1.31

1
1

.0524

.74

1
1

.0676

.96

1
1

10

.1476

1.17

11

.5800

1.15

ATTENTION : * = Chanseaent de sens .

zi

Retour au menu principal

PROGRAMME: PIEZO

L e c t u r e sur f i c h i e r s en c o u r s

V: M i se en f i en i er
Edition des rsultats

RESEAU EXEMPLE

0
O<*! | _ C U L

oeud de reference s
Pression de reference :

D1~:S

P RESSIOM3

10
0.0Q m

"SULTATS

! Noeud ! Cote (a)

44.60

112.64

45.00

114.01

69.01

116.61

72.11

1 :
2

Cote piezo (s)

P.Effective (a)

1
1

68.04

3 :

44.50

4 :

46.20

113.75

67.55

5 !

45.30

112.41

67.11

i
i

6i

44.75

113.73

68.98

45.00

111.92

66.92

48.00

114.24

66.24

48.00

110.36

62.36

:105.00

105.00

0.00

io

zz

Retour au menu principal

PROGRAMME THERMO

Saisie globale des donnes thermiques

* Epaisseur moyenne de l'isolant
( G si pas d'isolation )

(mm)

* Epaisseur moyenne des tubes (mm) s 25

Conduct vite des tubes ( 1 ambda.? W/m. des ) : 30

Valeur de LAMBDA SOL (W/m.dea)

1 .7

* Profondeur moyenne d'enfouissase

(m) s 0.8

* Temperature de surface du sol

(des C) s 0

Lecture sur fichiers des autres donnes

Quelle est la temperature d'injection

au noeud

3 (des.C) : 80

Quelle est la temperature d'injection au noeud

2 (des.C) s 80

Quelle est la temperature d'injection au noeud

A (des.C) : 80

Mise en fichier du rsultat

Eld i t i on du rsultat .

RESEAU EXEMPLE

0
CALCUL

D d S

T E M R El R A T LJ R 1ST S

Epaisseur moyenne<i solant ou conduites) s

25.0
Lambda (isolant ou conduites, W/m.des) : 30.000000
Lambda sol
<W/m.deg>
:
1.700000
Profondeur moyenne (m)

.80
Temperature de surface
s
0.

Temperatures d" injection s

Ti
Ti
Ti
Ti
Ti
Ti
Ti

<
(
(
(
<
(
<

1 )
2 )
3 )
4 >
5 )
6 )
7 )

=
=
=
=
=
=
=

80.0
80.0
80.0
80.0
80.0
80.0
80.0

RESULTATS :

Noeud

Temperature ( Dea.C )

80.

'80.

80.

80.

so.

80.

so.

i
i

80.

79.

10

79.

i
i

On constate ici des dperditions thermiques infrieures ou

gales I o C sans isolation (celle-ci n'est donc pas ncessaire).
Ceci est d en partie des vitesses d'coulement relativement grandes.
On se propose alors de faire une nouvelle simulation avec des
dbits diviss par dix.

Retour au menu principal

RESEAU EXEMPLE

Mod i -F i cat i on s

Faites votre choix :

A t t e n t i o n a bien m o d i f i e r T O U T E S l e s r u b r i q u e s c o n c e r n e e s
N o m b r e de c o n d u i t e s
Nombre de noeuds
Nombre de mailles
Prec: i s i on
ALPHA
Noeud de reference
Charge de reference
Cotes au sol
Caractristique des conduites
Caractristique des mailles
Debits aux noeuds
Surpresseurs et vannes
F in de mod i f.
Select i on

6
7
B
9
10
11
12
13
14
15
16
17
18

16

Modification des debits aux noeuds - entrez le numero du noeud

Entrez la valeur du debit

I <m3/s)

0.006

. . ,

RESEAU EXEMPLE

O A L C LJ L...
Preci s ion :

DES

O E B X T" S

.10 %

RESULTATS :
Noubre d'itrations 3
Conduite

Debit ( 3/s )

Vitesse ( G / S ) !

Reuariue

.0100

.14

.0036

.04

.OHO

.09

.0064

.07

1
1

1
1

.0114

.09

1
1

l
1

.0176

.14

1
1

1
1

.0257

.13

.0052

.07

i * !

.0068

.10

: :

io

.0148

.12

: :

11

.0580

.12

ATTENTION : * = Chansenent de sens .

RESEAU EXEMPLE

O A L... O U L_

D E S

T E M F* E R A T L J R E S

Epaisseur moyenne(i solant ou conduites) s

25.0
Lambda (isolant ou conduites, W/m.des) : 30.000000
Lambda sol (U/m.des)
" 1.700000
Profondeur moyenne (m)
: ,v .80
Temperature de surface
s
0

Temperatures d'injection
Ti
Ti
Ti
Ti
Ti
Ti
Ti

<
(
(
(
<
(
<

1 )
2 )
3 )
A )
5 )
6 >
7 >

=
=
=
=
=

80.0
80.0
80.0
80.0
80.0
80.0
80.0

RESULTATS :

Noeud

i
i

76.

77.

80.

/
^

1
1

77.

76.

80.

78.

78.

1
1

75.

io

9.

1
1

Temperature ( Des.C )

MINISTRE DE L'INDUSTRIE ET DE LA RECHERCHE

BUREAU DE RECHERCHES GOLOGIQUES ET MINIRES
SERVICE GOLOGIQUE NATIONAL
B.P. 6009 - 45060 Orlans Cedex - Tl.: (38) 63.80.01

ANNEXE 2
Exploitation sur le reseau
geothermal de la ville
de DAX (Landes)
Reseau restreint

I 84 AGI 164 AQI \

Pessac le 29 juin 1984 I

MODELISATION DU COMPORTEMENT HYDRAULIQUE ET

THERMIQUE DES RESEAUX MAILLES

SERVICE GOLOGIQUE RGIONAL AQUITAINE

Avenue Docteur Albert Schweitzer - 33600 PESSAC - Tl.: (56) 80.69.00

Exemple d'application : rseau restreint (Dax).

Cette tude, qui concerne une partie du rseau de

Dax, ne constitue pas une tude complte d'optimisation d'un
rseau
mais
simplement une illustration de ce que
sera
l'exploitation au logiciel LAUREL sur le rseau de Dax.
En fait, elle reprend une tude de ce rseau restreint dj faite
en 1979 et utilisant un logiciel de la CISI (compagnie de service
en informatique).
Pour le rseau restreint considr, le but de l'tude
est
d'envisager le recyclage des rejets existants.

On trouvera dans cette annexe :

- La
numrotation
et
le mai liage
initialement (figure 1 et 2 ) .

du

rseau

- Les donnes initiales pour ce rseau.

- Les fichiers et rsultats issus de la simulation 1
_

H M

M M

M M

J.

- Un plan d'ensemble du rseau de Dax.

considr

sh
(#>v
KK>m

^7

>,

4*

<--

5^^
' ^ ^

f~ora

\ALLAG\

Ta4 / e 8 u

Co* civ It

NOEU

lft I TI AL

d<LS

Co d-vj les (cloytcj)

N O H ' O FINAL

LONGUEUR
<M>

IAHETHE

60.0

.ISO

54.0

.125

50.0

.125

~f

S
5

-il
' C

7?

7
/$
f ^

_.

.12b
.125

60.0

,15

100.0

.125

700.0

.100

90.0

.125

444.0

.100

300.0

.080

260.0

.125

100.0

.125

-fu

260.0

.100

31

1*0.0

.125

100.0

.125

100.0

.ObO

17S.0

.125

?
0

j.

b6.0

fT

Z3

300.0

37

fo

4
7

(M)

*-

Ir

/O

-/
/-/

f f
'

If

1-f

-f-f
fZ.

/3
* ^

Y3
_//

0 X

/^
-?

< ?

J /

zr

-^ r
7^

*/

6
f

<
U
i

-fS

^?

-*?

il

H
S
JT

if

.125

100.0

.050

220.0

.ObO

t>0.0

lb

lbO.O

.ObO

. 3b0.0

*0H0

70.0

.lb

140.0

.Orto

84.0

.UU

1.0

i 0

10.0

Zo

Zl

z-f
ZZ

ij

Zl

.00

bO.O

.ObO

410.0

.060

24.0

OrtO

1*0*0

.060

160.0

.OrtO

100.0

.060

200.0

.000

150.0

.OBO
.00

Z
380.0

.oao

Ceyteli/t c

NOEUO FINAL

LONGUtUH

UIAMETHE
(H)

sr

HO.O

.ObO

Z?

120.0

.oso"

320.0

.080

1.0

.200

1.0

.200

vo.o

.OJO

120.0

.060

1.0

.200

370.0

.2S

1.0.0

.100

}?

100.0

.125

55
5/

be

lbO.O

luo

s?

SO.O

.050

13

*9

SUO.O

.10

lbO.D

.OO

30.0

.060

200.0

.'u*o

NOEUD INITIAL

zs
/

4L

Z
Z

Z7-

J /

Z8
Z3

Z9
BO ,

it

t?

41 '
33

it
5/

4o
4+

Rugosit de conduites
On a retenu, priori, pour ltsnsemble des conduites
un coefficient
k = 0,5 mm
" - ..

rim

0 E u I T

ALTlfUUE

(M)

L/S)

b.50

0.0

.50
' '

.bo

V.00

ZI

1}

V.50

"

10.00

12.SO

Yo

15.00

O'O

*rf

li.,50

0.0

'

"

2J.00

o*

11.00

0.50

"

(M)

0*

6.50

;
,

ALT1TUUE

0
ZY

i/oeuo

lV.bO

lb.50

' Y?

lt.SO

ff

21.50

Y?

23.00

> 0

-.6

2<.00

0.0

22.00

0*0

2<t.Q0

0*0

20.00

: o.o

IS

23.00

2*.00

! o.o

-.5

19 '

23.00

;"

o.o

lb.bO

0.0

0.0

15
3 .
37

28.0

6.00

b*

0.0
0.0

23.00

0.0

0.0

-.9

23.00

35

4*

- . 3

It

0.0

"

0.0

<*.00

H
15

it.SO
.00

24.00

'

0.0

JEOIT
(L/S>

"

3?

20.00

18.50

ALTITUDE
<M)

0.0

6 # 9

OEdlT
(L/S)

6.50

0.0

12.V0

-a.3

32.00

-2.2

b.bO

0.0

La rpartition des dbits aux

noeuds a t dfinie pralablement partir d'une tude sur
les consommations mesures in
si tu.

SIMULATION 1

<n'n

TITRE DE L'ETUDE
NUMERO DE LA SIMULATION
NOMBRE MAXIMUM DE CONDUITES
NOMBRE DE CONDUITES
NOMBRE DE NOEUDS
NOMBRE DE MAILLES
ERREUR
ALPHA
NOEUD DE REFERENCE
PRESSION DE REFERENCE

DAX RESTREINT
1

I
I

70
S3
42
I
14
1 .SO
[1.970 [
1
33
1 44.00

COTES AUX NOEUDS (n.NGF) i

8.50
8.50
8.50
9.00
6.50
9.30
VECTEUR DES TERMES INDEPENDANTS
11.00
0.00000
10.00
0.00000
12.50
0.00000
13.00
O.ooooo
$\tr\Are*)iwS
14.50
0.00000

14.50
0.00000
/ -//.c
19.50
0.00000
{ H&.l'tCJ
18.50
0.00000
19.50
0.00000
18.50
0.00000
18.50
0.00000
21.50
0.00000
23.00
0.00000
2 3 . OO
0.00000
20.50
0.00000
24.00
0.00000
22.00
O.ooooo DfLlh
a^x
20.00
o.ooooo
y- 1 * 1 y -*-*
23.00
-.00290
I
24.00
-.00380
n Ot{/(Ai .
24.00
0.00000
23.00
-.00690
23.00
0.00000
18.50
0.00000
20.00
0.00000
14.50
0.00000
6.00
-.00080
4.30
0.00000
12.90
0.00000
32.00
0.00000
4.50
0.00000
18.30
-.00080
0.00000
22.00
0.00000
24.00
24.00
0.00000
23.00
-.00060
0.00000
0.00000
-.00030
0.00000
-.00090
-.00050
0.00000
0.00000
0.00000
0.00000
.02800
0.00000
-.00830
-.00220
0.00000
0.00000
0.00000
0.00000
0.00000
-.00000

-p

CONOUITESCnumero - N.orlslnr - N.final 1

2 150.
.500
60.00
1
2
1
4 123.
.500
34.003
3 123.
.500
50.00
2
4
37
1
2
5
.
.500
3
00.00
2
3
4
1
2
5
.
.500
56.00
s
6
3
6 125.
.500
60.00
7
3 34 1 2 5 .
.500
100.00
8
.500
700.00
S 38 1 0 0 .
9
6
7 125.
.500
90.00
10
8 100.
.500
444.00
4
11
7
8
80.
.500
300.00
12
.500
260.00
7 10 1 2 5 .
13
8
9 125.
.500
100.00
14
9 10 1 0 0 .
.500
260.00
13
9 32,' ' 1 2 5 .
.500
190.00
16 10 13 125.
.500
100.00
50.
.500
100.00
17 11 12
.500
175.00
18 11 14 1 2 5 .
19 11 32 1 2 5 .
.300
10.00
20 12 32
50.
.500
100.00
50.
.500
220.00
21 13 14
.500
60.00
22 13 15 125.
14
30
50.
.500
130.00
23
80.
.500
360.00
24 15 16
.500
70.00
23 15 24 1 2 3 .
80.
.500
140.00
26 14 17
BO.
.300
84.00
27 14 18
.500
1.00
28 17 38 2 0 0 .
29 1B 19
80.
.500
80.00
60.
.500
410.00
30 18 20
.500
80.
24.00
31 19 20
60.
.500
160.00
32 19 21
80.
.500
160.00
33 20 22
60.
.500
34 21 31
100.00
35 22 23
80.
.500
200.00
36 22 24
80.
.500
150.00
37 23 39 2 0 0 .
.500
1.00
38 24 25
80.
.300
380.00
39 24 31
.500
60.
80.00
40 25 27
50.
.500
120.00 *
41 26 28
80.
.500
320.00
.500
1.00
42 26 40 2 0 0 .
.500
1.00
43 27 41 2 0 0 .
44 28 29
30.
.300
90.00
.300
80.
120.00
43 28 30
.500
1.00
44 29 42 2 0 0 .
47 33 34 2 5 0 .
.500
370.00
48 34 35 1 0 0 .
.300
130.00
49 34 37 125.
.500
100.00
50 35 36 1 0 0 .
.500
130.00
50.
.500
31 36 37
230.00
.500
500.00
32 38 39 1 0 0 .
80.
.500
150.00
33 39 40
60.
.500
230.00
54 40 41
.500
40.
200.00
S3 41 42

D. - K - L
0.00
0.00
0.00
0.00
0.00
0.00
0.00
. 0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00

- Suri. (CEJ

MAILLES (nuacro - nbrc". conduites- nu, conduites

1
4 49 -51 -50 -48
2 13
7 49 -4 3 4
9 12 16 22 24 26 28 -8
-1-2-3-3
3
4
4
10 -11 -9
'
3
11 13 14 -12
3
4
6
4 -14 13 -19 18 -21 -14
20 -19 17
3
7
29 31 -30
3
8
9 -28 -26 27 29 31 33 33 37 -32
9
10
7 -24 25 39 -34 -32 -29 -27
11 4 32 34 -39 -34 -33 -31
21 23 -45 -41 42 -53 -37 -33 36 -23 -22
12 11
38 40 43 -54 -53 -37 -35 34 .
8
13
41 44 44 -S3 -34 -42
4
14

DAX

RESTREINT

o <=* i_ o u i~
Precision

RESULTATS

.50

E: S

I:> E :

3 i

Moubre d"itrt ions : 10

1
1

Conduite

Debit ( a3/s )

Vitesse (B/S) !

Remarque

! !

1
1

.0039

.22

.0039

.32

.0004

.52

1
1

.0104

.84

1
1

.0039

.32

! *

1
1

.0104

.84

.0071

.58

1
1

.0042

.54

1
1

.0041

.33

i
i

1
1

10

.0025

.31

1
1

11

.0015

.30

1
1

12

.0026

.21

1
1

13

.0030

.24

1
1

14

.0023

.30

1
1

15

.0006

.os

16

.0003

.02

i
i

1
1

17

.oooo

.01

i
i

18

.OOO

.05

1
1

19

.OOO

.os

20

.oooo

.01

i
i

21

.0004

.19

1
1

22

.0009

.07

1
1

23

.0003

.13

i
i

1
1

24

.0008

.16

25

.0001

.01

1
1

.0021

.42

! *

27

.0013

.26

28

.0021

.07

29

.0004

.09

30

.0001

.03

31

.0002

.05

32

.0002

.07

33

.0003

.07

34

.0002

.07

1
1

35

.0008

: .*

.16

1
1

36

.0005

.11

1
1

37

.0008

.03

1
1

38

.0007

.13

39

.0002

.07

1
1

40

.0004

.18

1
1

41

.0008

.16

42

.0008

.03

43

.0005

.02

44

.0001

.08

i
i

1
1

45

.0003

.os

1
1

46

.0001

.oo

1
1

47

.0280

.57

1
1

48

.0095

1.21

I
1

49

.0113

.92

1
1

1
1

so

.0012

.16

51

.0010

.49

52

.0021

.27

53

.0013

.26

54

.0005

.17

55

.0001

.04

ATTENTION : * = Chansenent de sens .

DAX RESTREINT

CALCUL-

Noeud de reference s
Pression de reference

DES

33
s

44.00 m

RESULTATS :

Noeud

1
1

Cote ()

Cote piezo (i)

19

23.00

44.28

1
I

21.28

20

1
1

23.00

1
I

44.28

t
1

21.28

i
i

21

20.50

44.25

23.75

22

1
1

24.00

1
1

44.26

20.26

23

1
1

22.00

1
1

44.40

22.40

24

1
1

20.00

44.21

1
1

24.21

25

23.00

1
1

44.04

1
1

21.04

26

1
1

24.00

1
1

44.14

t
1

20.14

I
1

27

1
1

24.00

1
1

43.84

1
1

19.84

1
1

28

23.00

1
1

43.93

20.93

t
1

29

23.00

43.86

1
1

20.86

t
1

30

t
1

18.50

43.92

25.42

1
1

31

1
1

20.00

i
i

44.23

24.23

32

1
I

14.50

43.79

29.29

33

1
1

6.00

i
i

50.00

44.00

34

6.50

49.35

42.85

'

PRESSIONS

P.Effective ()

1
1

' 8.50

45.26

36.76

1
1

1
1

8.50

t
1

45.29

36.79

8.50

45.10

36.60

1
1

9.00

1
1

45.18

36.18

6.50

48.88

42.38

1
t

9.50

44.52

35.02

1
1

11.00

44.38

33.38

1
1

1
1

10.00

1
1

43.70

33.70

1
1

12.50

43.78

31.28

'

10

15.00

44.21

29.21

11

t
i

14.50

1
1

43.79

29.29

12

14.50

43.79

29.29

1
1

13

1
1

19.50

44.21

24.71

1
1

14 .

18.50

43.80

25.30

1
1

15

1
1

19.50

44.21

24.71

1
1

16

18.50

44.45

25.95

17

18.50

45.09:

26.59

1
1

18

21.50 :

44.29

'

. !

.'

22.79

35

12.90

45.33

1
1

32.43

1
1

36

32.00

1
1

45.26

1
1

13.26

37

6.50

48.19

41.69

I
t

38

1
1

18.50

45.08

26.58

39

1
1

22.00

44.40

1
1

22.40

40

24.00

44.14

20.14

41

24.00

1
1

43.84

1
1

19.84

42

23.00

f
t

43.86

20.86

DAX

RESTREINT

CALCUL

Noeud de reference :
Pression de reference s

DES

PRESSIONS

e 7

^.

! Cote <)

- -

|l?y

19

23.00

1
1

44.28

21.28

1
1

20

1
1

23.00

1
1

44.28

1
1

21.28

21

20.50

44.25

1
1

23.75

22

24.00

44.26

20.26

23

1
1

22.00

1
1

44.40

1
1

22.40

24

20.00

t
1

44.21

1
1

24.21

1
1

25

23.00

1
t

44.04

1
1

21.04

26

24.00

I
I

44.14

20.14

1
1

27

1
1

24.00

43.84

t
1

19.84

I
t

28

t
1

23.00

1
1

43.93

20.93

29

1
1

23.00

43.86

20.86

30

I
1

18.50

43.92

25.42

1
1

31

20.00

44.23

24.23

i'

43.79

29.29

33
44.00 m

RESULTATS :

S Noeud

! Cote riezo ()

! P.Effective ()

^r^r*ggggggsaaaggra*gggr=gTsg=?gg=?=g==g==

8.50

45.24

36.76

1
1

8.50

45.29

1
1

36.7?

1
1

1
1

1
1

8.50

45.10

I
t

36.60

'.

9.00

45.18

1
1

36.18

6.50

48.88

1
f

42.38

9.50

44.52

I
1

35.02

1
1

11.00

44.38

33.38

1
1

32

I
1

14.50

10.00

43.70

1
1

33.70

33

1
1

6.00

50.00

44.00

1
1

12.50

43.78

1
1

31.28

1
1

34

1
1

6.50

1
1

49.35

42.85

10

15.00

44.21

29.21

35

12.90

45.33

32.43

11

14.50

43.79

29.29

1
1

36

32.00

1
1

45.26

i
i

13.26

1
1

12

14.50

43.79

29.29

1
1

1
1

37

1
1

6.50

1
1

48.19

i
i

41.69

13

19.50

44.21

i
i

24.71

1
1

38

1
1

18.50

1
1

45.08

26.58

14 .

1
1

18.50

43.80

25.30

39

1
1

22.00

44.40

i
i

22.40

15

19.50

44.21

i
i

24.71

1
1

40

24.00

44.14

i
i

20.14

16

18.50

44.45

25.95

1
1

41

24.00

43.84

19.84

17

18.50

45.09

26.59

1
1

42

' 1
1

23.00

1
1

43.86

i
i

20.86

18

21.50

44.29

22.79

'

Avec

DAX RESTREINT

hold ti'oyi .

CALCUL. DES

TEMPERATURES

Epaisseur moyenne(Isolant ou conduites) :

30.0
Lambda (isolant ou conduites* W/m.des) :
.035000
Lambda sol (U/m.des)
:
2.000000
Profondeur moyenne (m)
:
.80
Temperature de surface
.:
0.

Temperatures d'injection s
Tl< 33 )

- 57.0

RESULTATS :

: Noeud

t
t

1
1

56.5

I
1

56.6

1
1

56.4

1
1

1
1

1
1

56.4

I
1

56.8

I
1

56.4

56.3

1
1

1
1

I
1

54.3

1
1

53.0

1
1

10

1
1

55.6

1
1

11

1
1

48.1

1
1

12

25.8

1
1

13

1
l

52.5

1
1

14

49.7

15

Teiperature ( Des.C ) !
'S3SSSSS3

52.6

16

55.5

17

55.8

1
1

18

t
1

55.2

1
1

1?

54.3

t
1

1
1

20

1
1

50.1

21

1
1

51.4

i
i

1
1

22

51.8

23

1
1

54.5

1
1

24

50.0

1
1

25

1
1

47.6

1
1

f
I

26'

54.0

1
1

27

49.1

28

52.1

1
1

1
1

2?

1
1

48.1

1
1

30

50.0

1
I

31

1
1

49.6

t
t

1
1

32

1
1

47.4

1
1

33

57.0

34

56.8

35

f
1

56.8

1
1

36

1
1

56.0

1
1

I
1

37

56.8

1
1

38

i
i

55.8

1
1

39

54.5

40

i
i

54.0

f
1

41

50.8

42

i
i

47.9

(
1

DAX

RESTREINT
CALCUL

DES

16 I

46.0

17 i

48.4

18

, 43.8

19

38.0

t
1

20

26.0

f
1

21

21.2

1
1

22

28.6

1
1

23

39.8

24

19.9

1
1

25 I

12.7

1
1

26

36.4

27

17.7

1
1

28

26.7

1
1

29

8.3

30

18.4

1
1

31

14.7

1
1

1
1

32

12.9

1
t

TEMREIR ATURIS

Epaisseur moyenne(isolant ou conduites) :

10.0
Lambda (isolant ou conduit esr W/m.deg) : 30.000000
Lambda sol
(W/m.deg)
: 1.800000
Profondeur moyenne <m)
s .80
Temperature de surface
s 0.

Temperatures drInjection
Ti< 33 )

= 57.0

hol /

RESULTATS :

1 Noeud !
!

Tewerature ( Des.C )
53.6

1
1

1
1

1
1

54.4

1
1

1
1

53.4

1
1

53.0

1
1

55.5

1
1

53.1

1
1

33

57.0

52.1

1
1

34

56.2

41.0

1
1

t
1

35

55.5

1
1

36.4

36

48.6

1
1

10

I
1

47.8

37

55.8

11

13.3

38

48.5

1
1

12

1
t

39

39.9

1
1

1
1

13

1
1

30.2

1
1

t
1

40

36.5

1
1

14

1
1

16.9

41 I

22.9

1
1

15

.0

30.0

t
I

42

8.2

NOTES

-Pression effective minimale : 13.26 raCE au noeud 36

-Temprature minimale : 25.8 C

au noeud 12 (avec isolation)

La Simulation suivante a pour but d'augmenter la circulation dans les

conduites 18,17,20 et dans la partie sud du rseau afin d'obtenir des
tempratures plus leves dans ces zones.

SIMULATION 2

Les modifications sont les suivantes :

- Conduites 17 et 20
- Conduite 19
- Surpresseur 10 mCE
d*j sens correct
- Surpresseur 30 mCE

: diamtre 100
: diamtre 80
sur la conduite 18 aprs rtablissement
sur la conduite 52

ti

TITRE DE L'ETUDE
NUMERO DE LA SIMULATION

DAX RESTREINT
2

NOMBRE MAXIMUM DE CONDUITES

NOMBRE DE CONDUITES
NOMBRE DE NOEUDS
NOMBRE DE MAILLES
ERREUR
ALPHA
NOEUD DE REFERENCE
PRESSION DE REFERENCE

70
55
42
14
.50
1.970
33
44.00

MAILLES (numero - nbre. conduites- num. conduite:

49 -51 -50 -48
4
7 49 -4
3
13
12 16 22 24 26
3
-1
2 -5 -3
4
4
10 -11 -9
3
5
11 13 14 -12
4
6
6 -14 15 -19 -18 -21 -16
7
3
20 -19 17
8
3
29 31 -30
9
9 -28 -26 27 29 31 33 35 37 -52
10
7 -24 25 39 -34 -32 -29 -27
11
6
32 34 -39 -36 -33 -31.
12 11
21 23 -45 -41 42 -53 -37'-35 36 -25 -22
13
8
38 40 43 -54 -53 -37 -35 36
14
6
41 44 46 -55 -54 -42

28

-8

CONDUITESnumero - N.orisine - N.final - D. - K - L. - Surp. (mCE) s

.500
T5
T"
rno~
.500
0 00
2
1
4 125.
54
VECTEUR DES TERMES INDEPENDANTS (Surp.- Debits) :
.500
0 00
2
3 125.
3
50 ,00
U.O00O
.500
0 00
2 37 125.
4
300 oo
0.00000
.500
0 00
4 125.
3
56 oo
5
0.00000
.500
0 00
6 125.
60 ,00
6
3
0.00000
.500
0. 00
100 oo
7
5 34 125.
0.00000
.500
0 00
700 oo
8
5 38 100.
-10.00000
.500
0 00
90 ,00
6
9
7 125.
0.00000
.500
0 00
444 ,00
10
6
8 100.
COTES AUX NOEUDS (m.NGF) :
.500
0 00
0.00000
300 ,00
11
8
80.
7
8.50
-30.00000
.500
0 00
260 ,00
12
7 10 125.
8.50
.500
0
00
0.00000
,00
100
13
9 125.
8
8.50
.500
0
00
0.00000
260 ,00
14
9 10 100.
9.00
.500
0 00
0.00000
190 ,00
15
9 32 125.
6.50
.500
0.00000
0' 00
100 ,00
16 10 13 125.
9.50
0 00
.500
0.00000
100 ,00
17 11 12 100.
11.00
.500
00
,00
175
0.00000
10
18 14 11 125.
10.00
.500
0.00000
0 00
80.
10 ,00
19 11 32
12.50
.500
0.00000
0 00
100 ,00
20 12 32 100.
15.00
.500
00
,00
0.00000
0
220
50.
21 13 14
14.50
.500
00
,00
-.00290
0
60
22 13 15 125.
14.50
.500
,00
-.00380
0 00
150
50.
23 14 30
19.50
.500
0.00000
0 00
360 .00
80.
24 15 16
18.50
.500
00
.00
-.00690
0
70 .00
25 15 24 125.
19.50
.500
00
0
0.00000
140 ,00
26 16 17
80.
18.50
.500
00
0
0.00000
84 .00
27 16 18
80.
18.50
.500
0 00
0.00000
1
28 17 38 200.
21.50
.500
00
,00
0
0.00000
80 .00
29 18 19
80.
23.00
.500
00
0
-.00080
410 .00
30 18 20
60.
23.00
.500
00
0
0.00000
24 ,00
31 1 9 20
80.
20.50
.500
00
0
0.00000
160 .00
24.00
32 19 21
60.
.500
0 00
0.00000
160 .00
22.00
33 20 22
80.
.500
00
0
0.00000
100 .00
20.00
60.
34 21 31
.500
00
0
-.00080
200 .00
23.00
80.
35 22 23
.500
00
0
0.00000
150 ,00
24.00
36 22 24
80.
.500
00
0
0.00000
1 ,00
24.00
37 23 39 200.
.500
00
0
0.00000
380 ,00
23.00
38 24 25
80.
.500
00
0
-.00060
80 .00
23.00
60.
.500
39 24 31
0 00
0.00000
18.50
50.
.500
120 .00
40 25 27
0 00
0.00000
.500
20.00
80.
320 .00
41 26 28
0 00
-.00030
.500
14.50
1 .00
42 26 40 200.
0 00
0.00000
.500
00
6.00
1 .00
43 27 41 200.
0
-.00090
.500
6.50
30.
90 .00
44 28 29
0 00
-.00050
.500
80.
12.90
120 .00
45 28 30
0 00
0.00000
.500
1 .00
32.00
46 29 42 200.
0 00
0.00000
.500
00
370 .00
6.50
47 33 34 250.
0
0.00000
.500
150 .00
18.50
48 34 35 100.
0 00
.500
0.00000
100 .00
49 34 37 125.
22.00
0 00
.500
.02800
150 .00
50 35 36 100.
24.00
0 00
.500
00
0.00000
250 .00
51 36 37
50.
24.00
.500
-.00830
500 .00
30 00
52 38 39 100.
23.00
.500
-.00220
150 .00
0 00
53 39 40
80.
0.00000
.500
00
250 .00
0
54 40 41
60.
0.00000
200 .00
0
55 41 42
40.
0.00000
0.00000
0.00000
' -.00000

DAX RESTREINT

o A L_ o
Prec i s i on -

E: S

I:> IE: B I T S

- 5 0 "

RESULTATS

"-

Noabre d'itrations 10

Conduite

1
1

.0032

1
1

!"

.0053

1
1

.0086

1
1

1
1

Debit ( H3/S )

Vitesse (/s)

.18

.0032

Renarme

I
1

.43

1
1

.70

.0032

.26

.0086

.70

.0089

.73

i
i

.0060

.77

1
1

.0031

.26

1
1

10

.0016

.21

1
1

11

.0010

.19

i
i

1
1

12

.0022

.18

1
1

1
1

13

.0043

.35

1
1

1
1

14

.0001

.01

1
1

1
1

15

.0042

.35

1
1

16

.0021

.17

1
1

17

.0012

.15

1
1

1
1

18

.0042

.35

i Surpresseur:

19

.0030

.60

20

.0012

.15

1
1
1

21

.0018

.94

22

.0005

.04

1
1

23

.0024

1.22

i
i

1
1

24

.0021

.42

i
i

1
1

25

.0027

.22

i
i

26

.0047

.94

.26

i
i

i
i

10.0 CE !

.52

.0026

27
1
1

28

.0047

.15

1
1

29

.0028

.56

i
i

30

.0000
.0006

.22

31

.0023

.45

32

.0005

33

.0029

34

.0005

35

.0061

!.

36

.0020
.0026

37

1
1

38

i
i

i
i

.57

i
i

i
i

.19

1.20

i
i

.51

.0001
.0061

.19

i
i

.0000
.0006

.13

39

.0005

.19

40

.0009

.47

41

.0030

.60

42

.0030

.10

43

.0018

.06

i
i

44

.0001

.18

i
i

45

.0024

.48

i
i

46

.0001

.00

1
1

47

.0280

.57

48

.0095
.0095

1.21

49
49

.0090
.0096

.78

1
1

50
50

.0012
.0012

.15

1
1

1
1

i
i

51
51

.0010
.0010

.51

52
52

.0108
.0108

1.37

.19

53

.0047

.94

54

.0017

.60

55

.0001

.10

ATTENTION : * = Chanseoent de sens .

i Surpresseur:

30.0 nCE !

i
i

1
1

f
1

DAX RESTREINT

CALCUL.

DES

PRESSIONS
19

Noeud de reference J

23.00

45.33

22.33

20

23.00

45.46

22.46

21

20.50

45.55

25.05

22

24.00

46.79

22.79

33

Pression de reference :

44.00 m

RESULTATS :
f

23

22.00

53.95

31.95

f
1

24

20.00

45.80

25.80

I
1

285 SSSSS SE8S25S"SSSSES5S3SZSSSS

Noeud. ! Cote ()

! Cote piezo <)

P.Effective (>

8.50

4.5o

33.00

1
1

25

23.00

45.95

22.95

1
1

8.50

46.52

38.02

26

24.00

50.66

26.66

1
1

B.SO

46.39

37.89

1
1

27

24.00

47.26

23.26

46.44

37.44

28

23.00

47.74

24.74

1
1

6.50

48.62

42.12

29

23.00

47.41

24.41

9.50

45.99

36.49

30

18.50 , :

47.05

28.55

11.00

45.91

34.91

31

20.00 .^ !

45.68

25.68

10.00

45.62

35.62

32

14.50

46.10

31.60

t
i

9.00

1
1

12.50

45.79

33.29

33

6.00

50.00

44.00

i
i

1
1

10

15.00

45.79

30.79

34

6.50

49.35

42.85

11

14.50

46.24

31.74

35

12.90

45.36

32.46

i
i

12

14.50

46.15

31.65

36

32.00

45.30

13.30

1
1

13

19.50

45.75

26.25

1
1

37

6.50

48.52

42.02

14

18.50

36.53

18.03

1
1

38

18.50

41.02

22.52

15

19.50

45.75

26.25

39

22.00

53.95

31.95

16

18.50

44.12

25.62

40

24.00

50.66

26.66

17

18.50

41.01

22.51

41

24.00

47.26

23.26

1
1

21.50

44.70

23.20

18

42

23.00

47.41

24.41

DAX RESTREINT

CALCUL. DES

TEMRERATURES

Epaisseur moyenne(Isolant ou conduites) :

30.0
Lambda (isolant ou conduites, W/m.des) :
.035000
Lambda sol
<U/m.des)
:
2.000000
Profondeur moyenne (m)
' : .80
Temperature de surface
:
0.

16

52.3

17

52.2

18

1
1

52.6

1
1

19

53.1

t
1

20

53.7

1
1

21

t
1

51.7

22

1
1

54.0

23

1
1

54.1

1
1

1
1

24

1
1

52.9

1
1

25

1
1

52.6

1
1

1
1

26

54.0

1
1

f
1

27

t
1

53.0

1
1

1
I

28

53.5

29

51.8

1
1

30

53.2

31

52.4

1
1

1
1

32

52.9

t
1

1
1

33

57.0

1
1

34

i
i

56.8

1
1

35

i
i

56.8

1
1

1
1

36

i
i

56.0

1
1

37

i
i

56.8

38

i
i

54.4

1
1

39

54.1

1
1

40

i
i

54.0

1
1

41

t
i

53.0

42

51.8

1
1

Temperatures d'injection
Ti( 33 )

57.0

Avec

RESULTATS :

Noeud

!
1
1

1
1

t
t

1
1

1
1

1
1

1
1

t
1

10

1
1

11

1
1

12

I
1

13
14
15

Tearerature ( Oes.C )
56.4
56.5

1
1

56.4

f
1

56.3

56.8

f
1

56.3

56.1

1
1

53.2

52.4

1
1

55.3

f
1

53.2

I
1

52.8

1
1

54.4

I
1

53.5

1
1

52.8

isola ko H

DAX RESTREINT
CALCUL.

DES

TEMPERATURES

Epaisseur moyenne(Isolant ou conduites) : 1 0 . 0

Lambda (isolant ou conduites W / m . d e s ) : 30.000000
1.800000
Lambda sol
<W/m.des)
.80
Profondeur moyenne (m)
0.
Temperature de surface

16

1
1

29.4

1
1

1
1

17

28.7

1
1

18

30.9

1
1

19

33.5

20

1
1

36.8

1
1

21

1
1

25.9

22

1
I

38.4

23

39.4

f
1

24

33.1

1
1

25

29.8

26

38.4

1
1

1
1

27

32.6

t
1

1
1

26

35.4

f
1

1
1

29

22.6

1
1

I
1

30 '

1
1

34.0

I
t

t
1

31-

29.7

32

32.4

1
I

33

1
1

57.0

34

1
1

56.2

35

55.5

1
1

1
1

36

1
1

48.5

1
I

37

55.7

38

t
1

40.9

39

39.4

40 '

38.4

1
1

41

1
1

32.6

1
1

42

22.5

Temperatures d'injection
Tl(

33 )

= 57.0

RESULTATS :

Noeud
1

!
!
!

2
3
4
5
6

!
!
:
:
!
:
!

7 !

Temperature ( Oes.C )
53.1
54.0

1
1

52.8

52.3

1
1

55.6

52.5
51.2

34.3
!
1
1

1
f
I

1
t

10

1
1

11

12
13
14

1
1

15 . i

30.7

46.1

t
1

33.6

1
1

31.4

1
1

41.2
34.9
32.4

1
1

i
i

NOTES

Pression minimale : 13.3 mCE au noeud 36

Tempratures minimum : 51.8 C et 22.5 C au noeud 42 (avec et sans isolation)

La simulation suivante a pour d'envisager une nouvelle alimentation au noeud 41.

SIMULATION 3

Modifications :
- Elimination des surpresseurs
- Alimentation dans la partie sud (conduite 47 : 33-41)

DAX RESTREINT

TITRE DE L'ETUDE
NUMERO DE LA SIMULATION

f~.'

>Le

c** / ^

NOMBRE MAXIMUM DE CONDUITES

NOMBRE DE CONDUITES
NOMBRE DE NOEUDS
NOMBRE DE MAILLES
ERREUR
ALPHA
NOEUD DE REFERENCE
PRESSION DE REFERENCE

3
!
70
t
55
i
42
!
14
!
.50
: 1.970
: .
33

44.00

MAILLES (numero - nbre. conduites- num. conduite:

1
4
49 -51 -50 -48
2 la
/ 49 -4
3
6
9 12 16 22 24 26
4
-1
2 -5 -3
3
3
10 -11 -9
4
4
11 13 14 -12
5
6
6 -14 15 -19 -18 -21 -16
3
20 -19 17
7
8
3
29 31 -30
9
9 -28 -26 27 29 31 33 35 37 -52
7 -24 25 39 -34 -32 -29 -27
10
6
32 34 -39 -36 -33 -31
-s
11
12 11
21 23 -45 -41 42 -53 -37 -35 36 -25 -22
8
38 40 43 -54 -53 -37 -35 36
13
6
41 44 46 -55 -54 -42
14

Ist/A

28

-8

CONDUITES<numero - N.origine - N.final - C . - K - L. - Surp. (mCE) s

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
45
46
47
48
49
SO
51
52
. 53

54
55

1
1
2
2
3
3
5
5
6
6

7
7
8
9
9
10
11
14
11
12
13
13
14
15
15
16
16
17
18
18
19
19
20
21
22
22
23
24
24
25
26
26
27
28
28
29
33
34
34
35
36
38
39
40
41

2
4
3
37
4
6
34
38
7
8
8
10
9
10
32
13
12
11
32
32
14
15
30
16
24
17
18
38
19
20
20
21
22
31
23
24
39
25
31
27
28
40
41
29
30
42
41
35
37
36
37
39
40
41
42

150.
125.
125.
125.
125.
125.
125.
100.
125.
100.
80.
125.
125.
100.
125.
125.
100.
125.
80.
100.
50.
125.
50.
80.
125.
80.
80.
200.
80.
60.
80.
60.
80.
60.
80.
80.
200.
80.
60.
50.
80.
200.
200.
30.
80.
200.
250.
100.
125.
100.
50.
100.
80.

60.
40.

.500
.500
.500
.500
.500
.500
.500
.500
.500
.500
.500
.500
.500
.500
.500
.500
.500
.500
.500
.500
.500
.500
.500
.500
.500
.500
.500
.500
.500
.500
.500
.500
.500
.500
.500
.500
.500
.500
.500
.500
.500
.500
.500
.500
.500
.500
.500
.500
.500
.500
.500
.500
.500
.500
.500

60.00
54.00
50.00
300.00
56.00
60.00
100.00
700.00
90.00
444.00
300.00
260.00
100.00
260.00
190.00
100.00
100.00
175.00
10.00
100.00
220.00
60.00
150.00
360.00
70.00
140.00
84.00
1 .00
80.00
410.00
24.00
160.00
160.00
100.00
200.00
150.00
1.00
380.00
80.00
120.00
320.00
1.00
1.00
90.00
120.00
1.00
370.00
150.00
100.00
150.00
250.00
500.00
150.00
250.00
200.00

0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
v 0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00
0.00

VECTEUR DES TERMES INDEPENDANTS (SurP.- Debits)

0.00000
0.00000
0.00000
0.00000
0.00000
0.00000
0.00000
0.00000
0.00000
0.00000
0.00000
0.00000
0.00000
0.00000
0.00000
0.00000
0.00000
0.00000
-.00290
-.00380
0.00000
-.00690
0.00000
0.00000
0.00000
0.00000
-.00080
0.00000
0.00000
0.00000
0.00000
-.00080
0.00000
0.00000
0.00000
-.00060
0.00000
0.00000
-.00030
0.00000
-.00090
-.00050
0.00000
0.00000
0.00000
0.00000
.02800
0.00000
-.00830
-.00220
0.00000
0.00000
0.00000
0.00000
0.00000
-.00000

,
i
' COTES AUX NOEUDS Ui.NGF) :
8.50
'<
8.50
8.50
9.00
6.50
9.50
11.00
10.00
12.50
15.00
14.50
14.50
19.50
18.50
19.50
18.50
1B.50
21.50
23.00
23.00
20.50
24.00
22.00
20.00
23.00
24.00
24.00
23.00
23.00
18.50
20.00
14.50
6.00
6.50
12.90
32.00
6.50
18.50
22.00
24.00
24.00
23.00

DAX

RESTREINT

o *=* t.... c
P r e c is ion

.50

i:> E;: S

I:> E : B I

Vitesse <B/S) !
.15

.22

7.

RESULTATS

Noabre d'itrations :

Debit ( 3/s )

Reaanue

Conduite

.0027

.0027

.0045

.36

.0072

.58

.0027

.22

!
i

.*

.0072

.58

.0033

.27

.0002

.79

- .0088

io

.0021

11

.0003

i
i

12

.0091

13

.0087

14

15

i
i

i
i

i
i

i
i

.72

.27

.05

.74

.71

.0040

.0047

16

.0132

1.07

17

.0013

1
1

47
.1/

18

.0047

19

.0034

20

21

i
i

i
i

i
i

i
i

.51

i
i

.38

1
1

i
i

1
1

.38

.7

i
i

.0013

.17

i
i

.0008

.43

i
i

22

.0148

1.21

i
i

23

.0039

24

.0015

.29

25

.0133

1.09

26

.0003

.05

1.97

i
i

i
i

i
i

1
1

27

.0012

.24

1
1

28

.0003

.01

1
1

29

i
i

.0016

.33

1
1

30

.0004

.13

1
1

31

.0018

.35

1
1

32

i
i

.0001

.os

33

.0021

.43

34

.0001

.05

1
1

1
1

.0051

1.02

36

.0024

.47

i
i

37

1
1

.0051

.16

1
1

38

1
1

.0108

2.16

1
1

39

.0001

.05

1
1

40

.0111

5.67

1
1

41

.0016

.31

42

1
1

.0016

.os

1
1

43

I
1

.0120

.38

X
*

1
1

44

.0028

3.93

1
1

45

.0039

.77

1
1

46

1
1

.0028

.09

47

1
1

.0280

.57

I
1

1
1

48

1
I

.0093

1.19

1
1

i
i

1
1

49

1
1

.0060

.49

1
1

50

.0010

.13

1
1

51

.0012

.59

<
1

1
1

52

.0065

.83

tf

53

.0116

2.31

54

.0132

4.66

55

.0028

2.21

i
1
1

i
1

-35

ATTENTION : * = Chargement de sens .

1
1

1
t

1
1

1
1

i
i

i
i

DAX

RESTREINT

CALCUL
Noeud de reference s
Pression de reference :

DIZS

PRESSIONS

33
44.00 m

RESULTATS :

Noeud
1

!
!

Cote ()

8.50 !

Cote riezo ()

P.Effective <)

-174.65

-183.15

B.50

-174.67

-183.17

8.50

-174.58

-183.08

-174.62

-183.62

-176.30

-182.80

9.00

6.50

9.50

-174.30

-183.80

1
1

1
1

7 !

11.00

-173.66

-184.66

10.00

-173.68

-183.68

1
1

12.50

-173.55

-186.03

10

15.00

-172.28

-187.28

11

14.50

-173.05

-187.55

1
1

!
1
t
1

1
1

12
13
14

15
16
17

!
!
:

!
:
!

14.50
19.50
18.50
19.50
18.50
18.50

:
!
!
!
!
!

-173.11
-170.73
-172.70
-168.82
-168.03
-168.16

!
!
!
!
i
!
!

-187.61
-190.23
-191.20
-188.32
-186.53
-186.66

19

23.00

-167.69

-190.69

20

23.00

-167.61

' -190.61

21

20.50

-167.70

-188.20

22

24.00

-167.01

-191.01

23

22.00

-161.88

-183.8

24

20.00

-167.71

-187.71

25

23.00

-124.92

-147.92

1
1

26

24.00

-142.50

-166.50

1
1

27

24.00

49.34

25.34

28

23.00

-143.88

-166.88

29

23.00

-12.91

-35.91

30

18.50

-145.65

-164.15

1
1

31

20.00

-167.70

-187.70

32

14.50

-173.17

-187.67

33

6.00

50.00

44.00

34

6.50

-176.40

-182.90

35

12.90

-180.26

-193.16

t
t

36

32.00

-180.31

-212.31

37

6.50

-176.07

-182.57

I
I

1
1

38

18.50

-168.16

-186.66

1
I

39

22.00

-161.88

-183.88

1
I

40

24.00

-142.50

-166.50

1
1

1
f

41

24.00

49.35

25.35

1
1
1
t
1

1
1

42
18

21.50

-167.90

-187.40

23.00

-12.91

-35.91

NOTES

Les tres faibles pressions (ngatives donc sans valeurs pratiques pour la
pression de rfrence donnes) sont dues ici des diamtres trop faibles
pour les conduites proches de la nouvelle alimentation (noeud 4 1 ) .

SIMULATION 4

Modifications :
-

Rtablissement du sens pour les conduites 10,11,12 et 38

Vannes sur 10 (0.2 mCE ) et sur 12 (0.5 mCE )
Surpresseur sur 38 (lOmCE)
Augmentation du diamtre des conduites 52(200mm),53(200mm)
54(125mm),55(125mm) et 40(125mm)

FIGURE

TITRE DE L'ETUDE
,NUMERO DE LA SIMULATION
NOMBRE
NOMBRE
[NOMBRE
NOMBRE
ERREUR
ALPHA

L\ty>a
de
a
S-O

DAX RESTREINT

MAXIMUM DE CONDUITES
DE CONDUITES
DE NOEUDS
DE MAILLES

NOEUD DE REFERENCE
IPRESSION DE REFERENCE

MAILLES (numero - nbre. conduites- num. conclu i t e :

'49 -5-1 -50 -48
4
2 13
7 49 -4
3
9 -12 16 22 24 26
3
-1
2 -5 -3
4
A
3 -10 11 -9
5
4 -11 13 14 12
6
6 -14 15 -19 -18 -21 -16
7
3
20 -19 17
8
3
29 31 -30
9
9 -28 -26 27 29 31 33 35 37 -52
10
7 -24 25 39 -34 -32 -29 -27
6
11
32 34 -39 -36 -33 -31
12 11
21 23 -45 -41 42 -53 -37 -35 36 -25 -22
8 -38 40 43 -54 -53 -37 -35 36
13
6
14
41 44 46 -55 -54 -42

28

-8

CONDUITESCnumero - N.orsine - N.final - D. - K - L. - SurP. <B)CE) S

0.00
60. oT
2 150.
0.00
.500
54. 00
2
A 125.
VECTEUR DES "TERMES"INDEPENDANTS (Surp.- Debits) :
0.00
.500
50. 00
3
3 125.
0.00000
0.00
.500
300. 00
A
2 37 125.
.50000
0.00
.500
56.
5
4 125.
3
0.00000
00
0.00
.500
60. 00

6 125.
3
.20000
0.00
100. 00
.500
7
5 34 125.
-.50000
0.00
700. 00
.500
8
5 38 100.
0.00000
COTES AUX NOEUDS (m.NGF) :
0.00
.500
90. 00
9
6
7 125.
0.00000
8.50
-.20
444. 00
.500
10
8
6 100.
0.00000
8.50
0.00
300. 00
80.
.500
11
8
7
0.00000
8.50
-.50
260. 00
12 10
.500
7 125.
0.00000
9.00
0.00
100,
13
.500
8
9 125.
00
0.00000
6.50
0.00
260. 00
14
.500
9 10 100.
0.00000
9.50
0.00
190, 00
15
.500
9 32 125.
-10.00000
11.00
0.00
100. 00
16 10 13 125.
.500
0.00000
10.00
0.00
100, 00
17 fi 12 100.
.500
0.00000
12.50
0.00
175,
18 14 1 1
.500
00
125.
0.00000
15.00
0.00
19 11 32
10, 00
.500
80.
0.00000
14.50
0.00
20 12 32 100.
100 00
.500
0.00000
14.50
0.00
21 13 14
220, 00
.500
50.
-.00290
19.50
0.00
60 00
.500
22 13 15 125.
-.00380
18.50
0.00
150 00
.500
50.
23 14 30
0.00000
19.50
0.00
360
00
.500
24 15 16
80.
-.00690
0.00
18.50
70 00
.500
25 15 24 125.
0.00000
0.00
18.50
140 00
.500
80.
26 16 17
0.00000
0.00
21.50
84
00
18
.500
80.
27 16
0.00
0.00000
23.00
1 00
.500
28 17 38 200.
0.00
0.00000
23.00
80 00
80.
.500
29 18 19
0.00
-.00080
20.50
410 00
60.
.500
30 18 20
0.00
0.00000
24.00
24 00
80.
.500
31 19 20
0.00
0.00000
22.00
160 00
60.
.500
32 19 21
0.00
0.00000
20.00
160 ,00
.500
80.
33 20 22
0.00
0.00000
23.00
100 00
60.
.500
34 21 31
0.00
-.00080
24.00
200 ,00
.500
80.
35 22 23
0.00
0.00000
24.00
150 ,00
.500
80.
36 22 24
0.00
0.00000
23.00
1 ,00
37 23 39 200.
.500
10.00
0.00000
23.00
380 ,00
.500
38 25 24 80.
0.00
-.00060
18.50
.00
80
.500
39 24 31 60.
0.00
0.00000
20.00
.00
120
.500
40 25 27 125.
0.00
0.00000
14.50
.500
320 .00
41 26 28 80.
0.00
-.00030
6.00
.00
.500
1
0.00
42 26 40 200.
0.00000
6.50
.00
1
.500
0.00
43 27 41 200.
-.00090
12.90
90 .00
0.00
.500
44 28 29 "30.
v
-.00050
32.00
0.00
120 .00
.500
45 28 30 80.
'
0.00000
6.50
.00
0.00
1
4*6 29 42 200.
.500
0.00000
18.50
.00
0.00
370
.500
47 33 41 250.
0.00000
22.00
0.00
150 .00
.500
48 34 35 100.
0.00000
24.00
.00
0.00
100
.500
49 34 37 125.
.02800
24.00
0.00
150 .00
.'500
50 35 36 100.
0.00000
23.00
.00
0.00
250
51 36 37 5 0 .
.500
-.00830
0.00
500 .00
.500
52 38 39 200.
-.00220
.00
0.00
150
.500
53 39 40 200.
0.00000
.00
0.00
250
.500
54 40 41 120.
0.00000
.00
200
.500
55 41 42 120.
0.00000
0.00000
0.00000
-.00000

AX

RESTREINT

O A L O U L_
Precision

O' BZ S

O EL" B I T S

. 5 0 /

RESULTATS s
Moibre d'itrations : 7
!

i
i

Debit ( 6i3/s )

Resanue

Conduite

n -

Vitesse (m/s)

.0022

.12

i
i

.0022

.18

.0036

.29

i
i

.0058

.47

1
t

1
1

.0022

.18

i
i

1
1

.0058

.47

.0047

.39

1
1

.0076

.97

.66

10

.0014

.18

! Vanne

11

.0004

.09

12

.0086

.70

! Vanne

13

.0079

.0031

.*

1
1

1
1

1
1

1
1

1
1

1
1

1
1

-.2 CE !
1

-.5 mCE

1
1

.64

|
1

.60

1
1

.26

i
i

I
1

At

11

.0047

15

.0032

16

.0133

1.08

1
1

17

.0009

.12

1
1

i
i

18

.0032

.26

1
1

19

.0023

.45

i
i

20

.0009

.12

i
i

21

.0010

! .

.52

22

.0151

1.23

i
t

i
i

|
1

23

.0022

1.10

24

.0028

i
i

25

.0123

1.00

1
1

26

.0044

.87

l
1

.56

1
1

27

.0016

.31

1
1

28

.0044

.14

i
i

i
i

29

i
i

.0006

1
1

.13

30

.0001

1
1

.04

1
1

1.
1

31

.0004

1
1

.07

1
1

32

.0010

1
1

.35

1
1

33

1
1

.0002

.05

1
1

34

1
1

.0010

1
t

.35

1
1

35

1
1

.0041

t
1

.82

f
1

1
1

36

1
1

.0033

.66

1
1

37

1
1

.0041

.13

1
1

38

1
1

.0080

1.59

3?

1
1

.0010

.35

1
1

40

.0083

1
1

41

i
i

.0019

1
1

42

.0019

43

1
1

44

>

i
i

1
1

1
1

1
t

1
1

I
1

1
1

10.0 mCE

1
1

.68

.38

i
i

.06

.0092

.29

i
i

1
1

.0008

i
i

1.07

45

.0022

.43

I
1

46

.0008

i
i

.02

I
1

1
1

47

.0280

.57

1
1

1
1

1
1

48

i
i

.0094

1.19

1
1

I
1

I
1

1
1

4?

.0046

.38

1
1

50

.0011

.13

1
1

51

i
i

.0011

.58

I
1

1
1

52

.0120

i
i

.38

1
1

1
1

t
I

1
1

53

.0161

i
i

.51

1
1

1
1

1
1

54

.0180

1.60

1
1

1
1

55

.0008

.07

1
1

i
i

ATTENTION : * = Changenent de sens .

Surpresseur!

1
1

1
1

1
1

I
1

1
1

1
1

I
1

1
I

1
1

DAX RESTREINT

CALCUL.

OES

PRESSIONS

1?

23.00

36.70

20

23.00

36.70

21

20.50

35.95

1
1

22

24.00

34.70

1
I

23

22.00

t
1

40.09

1
1

24

20.00

1
1

25

I
1

23.00

26

27

Noeud de reference s
Pression de reference

13.70
13.70

1
1

15.45

1
1

12.70

18.09

35.09

15.09

48.60

1
1

25.60

1
1

24.00

40.33

1
1

16.38

24.00

49.35

25.35

39.16

1
1

16.16

33
s

44.00 m

RESULTATS :

8.50

1
1

28.38

1
I

19.88

1
1

28

1
1

9.00

1
1

28.34

19.34

1
1

29

. 27.42

t
1

20.92

1
1

I
1

30

23.00

23.00

1
1

49.33

26.33

1
I

18.50

1
1

38.59

20.09

31

f
1

20.00

35.47

15.47

1
1

14.50

29.75

15.25

1
1

i
i

i
i

6.50

1
1

i
i

i
i

9.50

28.57

1
1

19.07

1
1

i
i

11.00

29.11

1
1

18.11

1
1

32

i
i

10.00

28.85

I
1

18.85

1
1

33

6.00

I
1

50.00

t
(

44.00

1
1

1
1

i
i

12.50

i
i

29.57

1
1

17.07

34

t
1

6.50

27.21

t
i

20.71

1
1

10

15.00

i
i

31.28

16.28

1
1

35

1
1

12.90

1
1

23.35

10.45

1
t

11

i
i

14.50

i
i

29.81

15.31

36

32.00

f
1

23.29

i
i

-8.71

f
1

12

i
i

14.50

29.78

I
1

15.28

37

t
t

6.50

t
1

27.41

i
i

20.91

1
1

39.55

i
i

21.05

13

19.50

32.86

13.36

38

I
1

18.50

14

18.50

29.97

11.47

39

1
1

22.00

40.09

18.09

1
1

40

1
1

24.00

40.38

t
t

16.38

25.35

1
1

26.33

15

19.50

34.14

14.44

16

18.50

36.94

18.44

41

I
f

24.00

I
1

49.35

17

18.50

39.55

21.05

1
1

42

23.00

1
1

49.33

18

21.50

36.73

15.23

DAX RESTREINT

CALCUL. DES

TEMPERATURES

Epaisseur moyenne<Isolant ou conduites) s 3 0 . 0

Lambda (isolant ou conduites W/m.des) s
.035000
Lambda sol
(U/m.des)
:
2.000000
Profondeur moyenne
(m)
:
.80
Temperature de surface
:
0.

Temperatures drinjection :
Tl< 33 .)

57.0

RESULTATS :

Noeud

Teverature ( Des.C )

1 !

54.8

54.9

1
1

55.1

1
1

1
1

55.0

t
1

55.7

1
1

1
1

55.2

1
1

55.7

t
I

54.7

1
1

1
1

54.9

1
1

10

1
1

55.9

1
1

11

I
1

54.8

t
t

12

54.2

1
1

13

1
1

55.9

1
1

14

55.1

1
1

15

55.9

16

56.1

17 i

56.3

1
1

18

55.8

I
1

f
1

19

53.3

1
1

20

50.2

1
1

21

52.7

1
1

22

56.4

1
1

23

56.7

1
1

24

56.1

1
1

25

56.7

1
1

26 i

56.7

t
I

1
1

27 !

56.8

1
1

28

55.6

1
I

29

55.1

30

55.4

31- !

52.4

1
f

32 !

54.4

33

57.0

34

55.0

35

54.9

1
1

36 ;

54.0

37 !

54.5

38

56.3

39

56.7

1
1

40

56.8

1
1

41

56.8

t
1

42 ;

55.1

DAX RESTREINT
DEIS

TEMP ER ATURES

Epaisseur moyenne<Isolant ou conduites) s

10.0
Lambda (isolant ou conduites, W/m.des) : 30.000000
2.000000
Lambda sol
(U/m.des)
.80
Profondeur moyenne (m)
0.
Temperature de surface

16

50.8

17

52.3

18

1
1

48.5

1
1

19

34.9

1
1

I
1

20

20.7

1
1

t
1

21

30.7

22

I
1

52.5

1
1

23

54.8

'

24

50.2

I
t

1
1

25

55.3

1
1

26

55.3

1
1

27

56.1

28

1
f

45.6

29

1
1

43.5

1
1

30

1
1

43.5

1
1

31-,

I
1

28.3

1
1

32

1
1

37.0

1
1

f
1

33

1
1

57.0

34

43.3

1
1

35

1
1

42.6

1
1

1
1

36

35.6

1
1

37

1
1

40.1

38

1
1

52.3

39

1
1

54.8

1
1

1
1

40

t
1

55.4

1
1

t
1

41

1
1

56.1

42

43.6

Temperatures d'Injection :
Ti< 33 )

}5olfiTfCA/\

= 57.0

RESULTATS :

Noeud

Tewerature ( Des.C )

41.7

42.2

'

43.8
42J
48.1
44.2

1
I

41.6

1
1

10

1
1

48.6

11

I
1

39.1

12

1
1

35.4

1
1
1

! .

1
1

47.2
40.3

!
!

t
1

13

48.9

t
1

14

41.2

15

49.1

NOTES

La pression encore trop faible au noeud 36 ncessiterait la mise en place

d'un surpresseur sur cette drivation.

Tempratures minimales : 50.2 C et 20.7 C au noeud 20 (avec et sans isolation)

On note donc dans ce cas une remarquable homognit des tempratures.