Bloc Scicos
fr - eng

##
Perte de charge thermohydraulique (Tuyau)

### Contenu

- Perte de charge thermohydraulique (Tuyau)
  - Description
  - Boîte de dialogue
  - Propriétés par défaut
  - Fonction d'interface
  - Modèle Modelica

### Description

Un bloc de PertDP représente un tuyau hydraulique avec des pertes de charge. Ce composant représente une perte de charge hydraulique où la perte de pression est directement proportionnelle au débit, état laminaire. Ce bloc a une direction; c.-à-d. la direction positive est quand le fluide entre par le port d'entrée (le noir). Ceci signifie que ce bloc suppose que le débit est positif si le fluide découle du port noir vers le port blanc. Les paramètres importants de ce composant sont la longueur, le diamètre du tuyau, les altitudes des ports d'entrée et de sortie, et quelques autres coefficients thermo-hydrauliques.

### Boîte de dialogue

Parametres du tuyau |
Longueur du tube : L (m) | 10
Diamètre interne du tube : D (m) | 0.2
Coefficient de perte de charge-frottement(S.U) : lambda | 0.03
Altitude entrée tuyauterie : z1 (m) | 0
Altitude sortie tuyauterie : z2 (m) | 0
Si 0, masse volumique imposée fu fluide : p_rho (kg/m3) | 0

- Longueur du tube : L (m) : Type 'vec' de taille -1. Longeur du tube
- Diamètre interne du tube : D (m) : Type 'vec' de taille -1. Diamètre interne du tube
- Coefficient de perte de charge-frottement(S.U) : lambda : Type 'vec' de taille -1. Coefficient de perte de charge-frottement(S.U)
- Altitude entrée tuyauterie : z1 (m) : Type 'vec' de taille -1. Altitude entrée tuyauterie (z1)
- Altitude sortie tuyauterie : z2 (m) : Type 'vec' de taille -1. Altitude entrée tuyauterie (z2)
- Si 0, masse volumique imposée fu fluide : p_rho (kg/m3) : Type 'vec' de taille -1. masse volumique imposée fu fluide

### Propriétés par défaut

- toujours actif: non
- direct-feedthrough: oui
- détection de passage à zéro: non
- mode: non
- nombre/taille des entrées régulières: 1 / 1
- nombre/taille des sorties sorties régulières: 1 / 1
- nombre/taille des entrées évènementielles: 0 /
- nombre/taille des sorties évènementielles: 0 /
- possède un état continu: non
- possède un état discret: non
- nom de la fonction modelica: PerteDP

### Fonction d'interface

`PerteDP.sci`

### Modèle Modelica

    class PerteDP
      parameter Real L "Longueur du tube";
      parameter Real D "Diametre interne du tube";
      parameter Real lambda "Coefficient de perte de charge-frottement (s.u.)";
      parameter Real z1 "Altitude entree tuyauterie";
      parameter Real z2 "Altitude sortie tuyauterie";
      parameter Real p_rho "Si > 0, masse volumique imposée du fluide";
    //  parameter Integer mode "Région IF97. 1:liquide - 2:vapeur - 4:saturation - 0:calcul automatique";

      parameter Real g = 9.80665 "Accélération de la pesanteur";
      parameter Real pi = 3.1415927;
      parameter Real eps = 1.e-0 "Limite inf. de la fonction carrée";
      Real khi "Coefficient de perte de charge hydraulique";
      Real deltaPf "Perte de charge par frottement";
      Real deltaP "Perte de charge totale";
      Real Q (start=500) "Débit massique";
      Real rho (start=998) "Masse volumique";
      Real Tm (start=290) "Température moyenne";
      Real Pm (start=1.e5) "Pression moyenne";
      Real Hm (start=100000) "Enthalpie spécifique moyenne";
      Real region (start=1) "Numéro de région IF97";
      Real ThermoCar;

    public
      PortPHQ1 C1 ;
      PortPHQ2 C2 ;

    equation

      C1.P - C2.P = deltaP;
      C2.Q = C1.Q;
      C2.H = C1.H;

      Q = C1.Q;

      /* Inversions de débit */
      0 = if (Q > 0.0) then (C1.H - C1.Hm) else (C2.H - C2.Hm);

      /* Calcul des pertes de charges */
      deltaP = deltaPf + rho*g*(z2 - z1);

      ThermoCar = if (abs(Q) > eps) then Q*abs(Q) else Q*eps;
      deltaPf = 8*khi*ThermoCar/(pi^2*D^4*rho);

        /* Formule de Darcy-Weisbach (Idel'cik p. 55). On suppose qu'on est en régime quadratique. */
      khi = lambda*L/D;

      /* Calcul des propriétés thermodynamiques */
      Pm = (C1.P + C2.P)/2;
      Hm = (C1.H + C2.H)/2;

      Tm = Hm/4187 + 273.15;
      region = 1;

      rho = 998;

    end PerteDP;