← Back to branch summary

~ubuntu-branches/ubuntu/hoary/scilab/hoary

« back to all changes in this revision

Viewing changes to macros/scicos_blocks/Hydraulics/VanneReglante.mo

  • Committer: Bazaar Package Importer
  • Author(s): Torsten Werner
  • Date: 2005-01-09 22:58:21 UTC
  • mfrom: (1.1.1 upstream)
  • Revision ID: james.westby@ubuntu.com-20050109225821-473xr8vhgugxxx5j
Tags: 3.0-12
http://bugs.debian.org/255869
changed configure.in to build scilab's own malloc.o, closes: #255869

Show diffs side-by-side

added added

removed removed

Lines of Context:
macros/scicos_blocks/Hydraulics/VanneReglante.mo
 
1
class VanneReglante 
 
2
  parameter Real Cvmax "CV max";
 
3
  parameter Real p_rho "Si > 0, masse volumique impos�e du fluide";
 
4
//  parameter Integer mode=0 "R�gion IF97. 1:liquide - 2:vapeur - 4:saturation - 0:calcul automatique";
 
5
 
 
6
  parameter Real eps=1.e-0 "Limite inf. de la fonction carr�e";
 
7
  Real Cv "Cv";
 
8
  Real Q(start=500) "D�bit massique";
 
9
  Real deltaP "Perte de charge singuli�re";
 
10
  Real rho(start=998) "Masse volumique";
 
11
  Real Tm(start=290) "Temp�rature moyenne";
 
12
  Real Pm(start=1.e5) "Pression moyenne";
 
13
  Real Hm(start=100000) "Enthalpie sp�cifique moyenne";
 
14
  Real region(start=1) "Num�ro de r�gion IF97";
 
15
  Real ThermoCar;
 
16
 
 
17
  Real Ouv ;
 
18
 
 
19
public 
 
20
//  EntreeReelle Ouv ;
 
21
  PortPHQ1 C1 ;
 
22
  PortPHQ2 C2 ;
 
23
 
 
24
equation 
 
25
  
 
26
  C1.H = C2.H;
 
27
  C1.Q = C2.Q;
 
28
  
 
29
  Q = C1.Q;
 
30
  deltaP = C1.P - C2.P;
 
31
  
 
32
  /* Inversions de d�bit */
 
33
  0 = if (Q > 0.0) then C1.H - C1.Hm else C2.H - C2.Hm;
 
34
  
 
35
  /* Calcul de la perte de charge */
 
36
 
 
37
  ThermoCar = if (abs(Q) > eps) then Q*abs(Q) else Q*eps;
 
38
  deltaP*Cv*abs(Cv) = 1.733e12*ThermoCar/rho^2;
 
39
  
 
40
  /* Calcul du Cv en fonction de l'ouverture */
 
41
  Cv = Ouv*Cvmax;
 
42
  
 
43
  /* Calcul des propri�t�s thermodynamiques */
 
44
  Pm = (C1.P + C2.P)/2;
 
45
  Hm = (C1.H + C2.H)/2;
 
46
  
 
47
  Tm = Hm/4187 + 273.15;
 
48
  region = 1;
 
49
  
 
50
  rho = 998;
 
51
 
 
52
end VanneReglante;