software de dimensionamiento y predicción de la caracterÃstica ...
software de dimensionamiento y predicción de la caracterÃstica ...
software de dimensionamiento y predicción de la caracterÃstica ...
Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
SOFTWARE DE DIMENSIONAMIENTO Y PREDICCIÓN DE LA CARACTERÍSTICA<br />
INSTALADA PARA VÁLVULAS DE CONTROL<br />
Ana Milena Borrego, Fabio Castrillón, Rafael Esteban Vásquez<br />
Universidad Pontificia Bolivariana<br />
Grupo <strong>de</strong> Investigación en Automática y Diseño A+D<br />
Circu<strong>la</strong>r primera # 70-01 bloque 11. Me<strong>de</strong>llín, Colombia.<br />
Tel. (+574) 4159020, fax. (+574) 4118779<br />
e-mail: rafavasquez@upb.edu.co<br />
Pa<strong>la</strong>bras C<strong>la</strong>ve: válvu<strong>la</strong>s <strong>de</strong> control, selección y <strong>dimensionamiento</strong> <strong>de</strong> válvu<strong>la</strong>s, características <strong>de</strong> flujo.<br />
Introducción<br />
La característica <strong>de</strong> flujo <strong>de</strong> <strong>la</strong> válvu<strong>la</strong> se <strong>de</strong>fine<br />
como <strong>la</strong> re<strong>la</strong>ción entre el flujo <strong>de</strong> salida <strong>de</strong> <strong>la</strong> válvu<strong>la</strong><br />
y <strong>la</strong> posición <strong>de</strong> apertura <strong>de</strong> ésta (0-100%).<br />
El tipo <strong>de</strong> característica que presente una válvu<strong>la</strong><br />
una vez instada <strong>de</strong>pen<strong>de</strong> exclusivamente <strong>de</strong> <strong>la</strong>s<br />
re<strong>la</strong>ciones <strong>de</strong> caída <strong>de</strong> presión entre ésta y <strong>la</strong> total<br />
disponible en el sistema para un flujo dado. Esta<br />
re<strong>la</strong>ción viene dada por el factor α:<br />
∆P<br />
α= ∆ P<br />
V<br />
T Q<br />
∆ P = ∆ P +∆ P<br />
T V Q L Q<br />
Don<strong>de</strong>:<br />
∆P v ⏐ Q : caída <strong>de</strong> presión a través <strong>de</strong> <strong>la</strong> válvu<strong>la</strong>, para<br />
un caudal Q.<br />
∆P T : caída <strong>de</strong> presión dinámica total <strong>de</strong>l sistema, a<br />
un valor <strong>de</strong> caudal Q.<br />
∆P L ⏐ Q : caída <strong>de</strong> presión dinámica en <strong>la</strong> tubería,<br />
accesorios y equipos <strong>de</strong>ntro <strong>de</strong>l sistema (diferentes<br />
a <strong>la</strong> válvu<strong>la</strong>), para un caudal Q.<br />
Existen dos c<strong>la</strong>ses <strong>de</strong> características <strong>de</strong> flujo según<br />
el valor <strong>de</strong> α:<br />
Característica Inherente <strong>de</strong> Flujo (α=1): se <strong>de</strong>fine<br />
como <strong>la</strong> re<strong>la</strong>ción “teórica” entre el flujo <strong>de</strong> salida <strong>de</strong><br />
<strong>la</strong> válvu<strong>la</strong> y <strong>la</strong> posición <strong>de</strong> apertura <strong>de</strong> ésta, cuando<br />
<strong>la</strong> caída <strong>de</strong> presión a través <strong>de</strong> <strong>la</strong> válvu<strong>la</strong> es<br />
constante e igual a <strong>la</strong> caída <strong>de</strong> presión total <strong>de</strong>l<br />
sistema. Esta característica es única para cada tipo<br />
<strong>de</strong> válvu<strong>la</strong> (lineal, isoporcentual y apertura rápida).<br />
Característica Insta<strong>la</strong>da <strong>de</strong> Flujo<br />
∆P<br />
V<br />
∆P<br />
T Q<br />
< 1 ó α < 1 :<br />
es <strong>la</strong> re<strong>la</strong>ción “real” entre el flujo <strong>de</strong> salida y <strong>la</strong><br />
posición <strong>de</strong> apertura, cuando <strong>la</strong> válvu<strong>la</strong> está en<br />
funcionamiento <strong>de</strong>ntro <strong>de</strong>l sistema y <strong>la</strong> caída <strong>de</strong><br />
presión a través <strong>de</strong> el<strong>la</strong> es menor que <strong>la</strong> caída <strong>de</strong><br />
presión total <strong>de</strong>l sistema.<br />
Esta característica <strong>de</strong>pen<strong>de</strong> <strong>de</strong> <strong>la</strong> característica<br />
inherente y <strong>de</strong> <strong>la</strong>s condiciones <strong>de</strong> operación <strong>de</strong>l<br />
proceso (caídas <strong>de</strong> presión, caudal, propieda<strong>de</strong>s<br />
físicas <strong>de</strong>l fluido, etc.).<br />
La característica inherente <strong>de</strong> una válvu<strong>la</strong> se<br />
modifica en mayor o menor grado, <strong>de</strong>pendiendo <strong>de</strong><br />
<strong>la</strong> razón <strong>de</strong> <strong>la</strong>s caídas <strong>de</strong> presión entre <strong>la</strong> válvu<strong>la</strong> y<br />
<strong>la</strong> total. A medida que α se hace más pequeño, una<br />
válvu<strong>la</strong> <strong>de</strong> característica inherente isoporcentual<br />
tien<strong>de</strong> a ser lineal, y una característica inherente<br />
lineal tien<strong>de</strong> a ser <strong>de</strong> apertura rápida.
Una válvu<strong>la</strong> <strong>de</strong> control <strong>de</strong>ntro <strong>de</strong> un sistema pue<strong>de</strong><br />
mantener, corregir o empeorar <strong>la</strong> linealidad <strong>de</strong>l <strong>la</strong>zo<br />
<strong>de</strong> control, <strong>de</strong>pendiendo <strong>de</strong> <strong>la</strong> característica <strong>de</strong> flujo<br />
que exhiba, una vez insta<strong>la</strong>da en el proceso.<br />
Es por esta razón que el ingeniero <strong>de</strong> control <strong>de</strong>be<br />
analizar en <strong>de</strong>talle cuál será <strong>la</strong> curva <strong>de</strong> ganancia<br />
real <strong>de</strong> <strong>la</strong> válvu<strong>la</strong> y <strong>la</strong> <strong>de</strong>l sistema -antes <strong>de</strong> adquirir<br />
una <strong>de</strong>terminada característica inherente- ,<br />
buscando que esta última “compense” <strong>la</strong>s nolinealida<strong>de</strong>s<br />
<strong>de</strong>l <strong>la</strong>zo.<br />
El programa <strong>de</strong>sarrol<strong>la</strong>do en Visual Basic 6.0<br />
permite comprobar cómo los parámetros <strong>de</strong> diseño<br />
<strong>de</strong> una válvu<strong>la</strong> (caudal máximo, rango <strong>de</strong> operación,<br />
y característica inherente) varían <strong>de</strong> acuerdo a <strong>la</strong>s<br />
condiciones <strong>de</strong> operación <strong>de</strong>l proceso. Un gráfico<br />
muestra cómo <strong>la</strong> característica insta<strong>la</strong>da se <strong>de</strong>svía<br />
<strong>de</strong> <strong>la</strong> inherente a medida que el parámetro α se<br />
hace más pequeño. Así mismo el usuario pue<strong>de</strong><br />
dimensionar cualquier tipo <strong>de</strong> válvu<strong>la</strong> (lineal,<br />
isoporcentual y <strong>de</strong> apertura rápida) ingresando <strong>la</strong>s<br />
variables operativas y un factor <strong>de</strong><br />
sobre<strong>dimensionamiento</strong>.<br />
Metodología<br />
Se tienen <strong>la</strong>s ecuaciones para <strong>la</strong>s condiciones<br />
inherentes y <strong>de</strong> insta<strong>la</strong>ción <strong>de</strong> una válvu<strong>la</strong> <strong>de</strong><br />
control. Ambas parten <strong>de</strong> <strong>la</strong> ecuación general para<br />
<strong>la</strong> característica <strong>de</strong> una válvu<strong>la</strong> <strong>de</strong> control.<br />
Característica Inherente<br />
La ecuación para <strong>la</strong> característica <strong>de</strong> flujo inherente<br />
se <strong>de</strong>riva <strong>de</strong> <strong>la</strong> expresión para el caudal a través <strong>de</strong><br />
<strong>la</strong> válvu<strong>la</strong>, teniendo en cuenta que <strong>la</strong> caída <strong>de</strong><br />
presión a través <strong>de</strong> ésta, para este caso, no<br />
<strong>de</strong>pen<strong>de</strong> <strong>de</strong>l caudal (∆P v constante).<br />
∆Pv<br />
Q(gpm) = CVmaxf ( x *) (1)<br />
G<br />
f<br />
Don<strong>de</strong>:<br />
Cv max : capacidad máxima <strong>de</strong> flujo <strong>de</strong> <strong>la</strong> válvu<strong>la</strong>,<br />
1/ 2<br />
gpm / psi<br />
f(x*): función que <strong>de</strong>termina <strong>la</strong> característica<br />
inherente <strong>de</strong> <strong>la</strong> válvu<strong>la</strong>.<br />
Gf: gravedad específica <strong>de</strong>l fluido a <strong>la</strong> temperatura<br />
<strong>de</strong> flujo.<br />
La Tab<strong>la</strong> 1 muestra <strong>la</strong> función f(x*) para cada tipo <strong>de</strong><br />
válvu<strong>la</strong>.<br />
Tab<strong>la</strong> 1. Tipo <strong>de</strong> función f(x*) según <strong>la</strong> característica<br />
inherente <strong>de</strong> <strong>la</strong> válvu<strong>la</strong><br />
Tipo <strong>de</strong> Válvu<strong>la</strong><br />
f(x*)<br />
Lineal x*<br />
(x* 1)<br />
Isoporcentual R −<br />
Apertura Rápida x*<br />
Don<strong>de</strong>:<br />
R: parámetro <strong>de</strong> rango <strong>de</strong> <strong>la</strong> válvu<strong>la</strong>.<br />
x<br />
x * = : fracción <strong>de</strong> apertura <strong>de</strong> <strong>la</strong> válvu<strong>la</strong>.<br />
x T<br />
x: carrera <strong>de</strong>l vástago.<br />
x T : carrera total <strong>de</strong>l vástago.<br />
Si se evalúa esta expresión para flujo máximo<br />
(Q=Q max , f(x*)=1), se tiene que:<br />
Q (gpm) = C<br />
max<br />
Vmax<br />
∆P<br />
G<br />
f<br />
v<br />
Dividiendo el caudal Q entre Q max se tiene <strong>la</strong><br />
expresión para <strong>la</strong> característica inherente <strong>de</strong> <strong>la</strong><br />
válvu<strong>la</strong>:<br />
Q<br />
Q<br />
max<br />
( )<br />
= f x*<br />
Característica <strong>de</strong> Insta<strong>la</strong>ción<br />
Para obtener <strong>la</strong> expresión matemática <strong>de</strong> <strong>la</strong><br />
característica insta<strong>la</strong>da <strong>de</strong> una válvu<strong>la</strong> <strong>de</strong> control, se<br />
realizaron <strong>la</strong>s siguientes suposiciones:<br />
• Flujo turbulento<br />
• Caída <strong>de</strong> presión total constante (bomba <strong>de</strong><br />
cabeza p<strong>la</strong>na).<br />
De <strong>la</strong> mecánica <strong>de</strong> fluidos, <strong>la</strong> caída <strong>de</strong> presión a<br />
través <strong>de</strong> una línea en régimen turbulento viene<br />
dada por:<br />
Don<strong>de</strong>:<br />
2<br />
∆ PL Q<br />
= KLGfQ<br />
(2)<br />
K L : constante que cuantifica <strong>la</strong>s pérdidas por fricción<br />
<strong>de</strong>bidas a accesorios, equipos y fricción en <strong>la</strong><br />
tubería.
Para régimen turbulento:<br />
K<br />
=<br />
∆P<br />
L<br />
L 2<br />
GQ<br />
f<br />
∆ PT = ∆ PV Q<br />
+∆ PL Q<br />
= cte (3)<br />
Reemp<strong>la</strong>zando (1) y (2) en (3) se obtiene:<br />
GQ<br />
∆ P = + K G Q<br />
Q =<br />
2<br />
f<br />
T 2 2<br />
L f<br />
Cvmax<br />
f (x*)<br />
⎛<br />
Gf<br />
⎜<br />
⎝C<br />
∆P<br />
T<br />
2 2<br />
vmax<br />
2<br />
1 ⎞<br />
+ KL⎟<br />
f (x*) ⎠<br />
(4)<br />
Este es el caudal “insta<strong>la</strong>do” a través <strong>de</strong> una válvu<strong>la</strong><br />
<strong>de</strong> control. La expresión <strong>de</strong>muestra que su valor no<br />
<strong>de</strong>pen<strong>de</strong> únicamente <strong>de</strong> <strong>la</strong> fracción <strong>de</strong> apertura <strong>de</strong><br />
<strong>la</strong> válvu<strong>la</strong> x* (y por lo tanto <strong>de</strong> su caída <strong>de</strong> presión),<br />
sino también <strong>de</strong> <strong>la</strong>s pérdidas dinámicas en <strong>la</strong> línea.<br />
Evaluando el caudal para flujo máximo (Q=Q max ,<br />
f(x*)=1),<br />
Inherente<br />
Evaluando (1) a caudal máximo y mínimo, y<br />
teniendo en cuenta que para condiciones inherentes<br />
<strong>la</strong> caída <strong>de</strong> presión en <strong>la</strong> válvu<strong>la</strong> es constante, se<br />
tiene que:<br />
Q f<br />
max<br />
( x* = 0.95)<br />
=<br />
Q f x* = 0.05<br />
min<br />
( )<br />
Insta<strong>la</strong>do<br />
Evaluando <strong>la</strong> expresión para el flujo <strong>de</strong> insta<strong>la</strong>ción<br />
en x*=0.95 y 0.05 respectivamente, y <strong>de</strong> <strong>la</strong><br />
expresión para K L , se tiene:<br />
2 2<br />
max =<br />
f op L Qop<br />
= ×<br />
2 2<br />
min f op L Qop<br />
( )<br />
( )<br />
Q f(x* 0.95) G Q +∆ P f x* = 0.05<br />
Q f(x* = 0.05) G Q +∆ P f x* = 0.95<br />
Se observa el rango <strong>de</strong> operación <strong>de</strong> insta<strong>la</strong>ción<br />
siempre será menor o igual al inherente; es por esta<br />
razón que una válvu<strong>la</strong> <strong>de</strong> control se <strong>de</strong>be<br />
sobredimensionar, con el fin <strong>de</strong> lograr un control<br />
lineal en un amplio rango <strong>de</strong> valores <strong>de</strong> caudal.<br />
Q<br />
max<br />
=<br />
∆P<br />
⎛ 1<br />
Gf<br />
⎜<br />
⎝C<br />
T<br />
2<br />
vmax<br />
⎞<br />
+ KL⎟<br />
⎠<br />
(5)<br />
Dividiendo (4) entre (5), se llega a:<br />
⎛ 1 ⎞<br />
⎜ + K<br />
2 L ⎟<br />
Q Cvmax<br />
∴ =<br />
⎝<br />
⎠<br />
Qmax<br />
⎛ 1 ⎞<br />
⎜KL +<br />
2 2 ⎟<br />
⎝ Cvmax<br />
f (x*) ⎠<br />
La re<strong>la</strong>ción Q/Q max <strong>de</strong>termina <strong>la</strong> característica <strong>de</strong><br />
insta<strong>la</strong>ción <strong>de</strong> <strong>la</strong> válvu<strong>la</strong>. Se pue<strong>de</strong> observar que<br />
cuado <strong>la</strong>s pérdidas en <strong>la</strong> línea son <strong>de</strong>spreciables<br />
(K L ≈0), <strong>la</strong> característica insta<strong>la</strong>da se convierte en <strong>la</strong><br />
inherente, conservándose <strong>la</strong> curva <strong>de</strong>l fabricante<br />
para una caída <strong>de</strong> presión dada.<br />
Rango <strong>de</strong> operación<br />
El rango <strong>de</strong> operación se <strong>de</strong>fine como <strong>la</strong> re<strong>la</strong>ción<br />
entre el caudal mínimo y máximo contro<strong>la</strong>bles que<br />
tiene una válvu<strong>la</strong>. Usualmente el caudal máximo se<br />
da para un 95% apertura <strong>de</strong>l vástago, y el flujo<br />
mínimo para un 5% respectivamente.<br />
Capacidad <strong>de</strong> operación <strong>de</strong> diseño <strong>de</strong> <strong>la</strong> válvu<strong>la</strong> <strong>de</strong><br />
control<br />
La capacidad <strong>de</strong> operación (C v]op ) se <strong>de</strong>fine como<br />
los galones por minuto <strong>de</strong> agua que pue<strong>de</strong>n pasar a<br />
través <strong>de</strong> una válvu<strong>la</strong> con una fracción <strong>de</strong> apertura<br />
x* y con una caída <strong>de</strong> presión <strong>de</strong> 1 psi.<br />
op<br />
V max<br />
( )<br />
Cv = C f x * (6)<br />
Reemp<strong>la</strong>zando (1) en (6), se llega a:<br />
Cv =Cv *factor<br />
req<br />
op<br />
Para fines <strong>de</strong> <strong>dimensionamiento</strong>, <strong>la</strong> capacidad <strong>de</strong><br />
diseño o capacidad requerida ( Cv ) <strong>de</strong> <strong>la</strong> válvu<strong>la</strong><br />
req<br />
<strong>de</strong>be ser mayor que Cv .<br />
op<br />
El valor para el factor <strong>de</strong> sobre<strong>dimensionamiento</strong> se<br />
pue<strong>de</strong> fijar a partir <strong>de</strong> recomendaciones prácticas.<br />
El comprador <strong>de</strong>be elegir una válvu<strong>la</strong> comercial<br />
cuya capacidad sea mayor o igual al<br />
Cv calcu<strong>la</strong>do.<br />
req
Resultados y Discusión<br />
Uso <strong>de</strong>l programa e ilustración <strong>de</strong> <strong>la</strong>s condiciones<br />
inherentes e insta<strong>la</strong>das <strong>de</strong> una válvu<strong>la</strong> <strong>de</strong> control<br />
El sistema <strong>de</strong> control <strong>de</strong> nivel que se muestra en <strong>la</strong><br />
Figura 1, tiene los siguientes parámetros <strong>de</strong><br />
operación:<br />
Figura 2. Ingreso <strong>de</strong> datos, cálculo <strong>de</strong> condiciones <strong>de</strong><br />
operación y <strong>dimensionamiento</strong> <strong>de</strong> <strong>la</strong> válvu<strong>la</strong> <strong>de</strong> control.<br />
Caudal: 700 gpm.<br />
Caída <strong>de</strong> presión total: 12 psi<br />
Caída <strong>de</strong> presión en <strong>la</strong> línea: 5 psi<br />
Gravedad específica <strong>de</strong>l fluido: 0.94<br />
El tipo <strong>de</strong> válvu<strong>la</strong> es lineal.<br />
Figura 1. Sistema <strong>de</strong> control <strong>de</strong> nivel<br />
Figura 3. Cálculo <strong>de</strong> condiciones inherentes e insta<strong>la</strong>das<br />
<strong>de</strong> <strong>la</strong> válvu<strong>la</strong> <strong>de</strong> control<br />
Ingresando los datos anteriores, el programa calcu<strong>la</strong><br />
los siguientes datos para <strong>la</strong> válvu<strong>la</strong> <strong>de</strong> control:<br />
• Caída <strong>de</strong> presión<br />
• Capacidad <strong>de</strong> operación<br />
• Capacidad <strong>de</strong> operación diseño requerida<br />
(ingresando un factor <strong>de</strong> sobre<strong>dimensionamiento</strong> <strong>de</strong><br />
1.3)<br />
• El caudal máximo y <strong>la</strong> habilidad <strong>de</strong> rango<br />
inherente e insta<strong>la</strong>da son también calcu<strong>la</strong>dos.<br />
Figura 4. Desviación <strong>de</strong> <strong>la</strong> característica inherente e<br />
insta<strong>la</strong>da según <strong>la</strong> re<strong>la</strong>ción <strong>de</strong> caídas <strong>de</strong> presión entre <strong>la</strong><br />
válvu<strong>la</strong> y el sistema.<br />
El usuario pue<strong>de</strong> observar cómo <strong>la</strong> gráfica <strong>de</strong> <strong>la</strong><br />
característica insta<strong>la</strong>da se <strong>de</strong>svía <strong>de</strong> <strong>la</strong> inherente <strong>de</strong><br />
acuerdo a <strong>la</strong> re<strong>la</strong>ción <strong>de</strong> caídas <strong>de</strong> presión entre <strong>la</strong><br />
válvu<strong>la</strong> y <strong>la</strong> total (P v /P T en <strong>la</strong> pantal<strong>la</strong>).
Conclusiones<br />
Las válvu<strong>la</strong>s son los elementos finales <strong>de</strong> control<br />
más comunes que se encuentran en los procesos<br />
industriales, por lo tanto <strong>de</strong> su a<strong>de</strong>cuado<br />
<strong>dimensionamiento</strong>, <strong>de</strong>pen<strong>de</strong> su correcta operación<br />
en el proceso para el cual son seleccionadas.<br />
La característica insta<strong>la</strong>da <strong>de</strong> una válvu<strong>la</strong> <strong>de</strong> control<br />
cambia consi<strong>de</strong>rablemente con <strong>la</strong>s condiciones <strong>de</strong><br />
operación <strong>de</strong>l proceso, por lo tanto su ganancia está<br />
cambiando y <strong>de</strong> esta forma se afectan los<br />
parámetros <strong>de</strong> operación <strong>de</strong> todo el sistema<br />
automático <strong>de</strong> control.<br />
Es necesario pre<strong>de</strong>cir <strong>la</strong> característica insta<strong>la</strong>da <strong>de</strong><br />
<strong>la</strong> válvu<strong>la</strong> a partir <strong>de</strong> <strong>la</strong> característica inherente y <strong>de</strong><br />
<strong>la</strong>s condiciones <strong>de</strong> operación que tiene el proceso.<br />
Al ser un proceso dispendioso, el uso <strong>de</strong><br />
herramientas informáticas para el diseño se hace<br />
<strong>de</strong>seable para disminuir el tiempo <strong>de</strong> selección y<br />
para garantizar una a<strong>de</strong>cuada operación en el<br />
futuro.<br />
Un correcto <strong>dimensionamiento</strong> <strong>de</strong> una válvu<strong>la</strong> <strong>de</strong><br />
control pue<strong>de</strong> garantizar <strong>la</strong>s características lineales<br />
<strong>de</strong> todo el sistema, al hacer uso <strong>de</strong> <strong>la</strong> predicción <strong>de</strong><br />
<strong>la</strong> característica insta<strong>la</strong>da para corregir otras nolinealida<strong>de</strong>s<br />
que existan en los diferentes<br />
componentes.<br />
El trabajo actual fue <strong>de</strong>sarrol<strong>la</strong>do para cuando <strong>la</strong><br />
bomba <strong>de</strong>l proceso es <strong>de</strong> cabeza p<strong>la</strong>na, por lo tanto<br />
se va a <strong>de</strong>sarrol<strong>la</strong>r una ampliación <strong>de</strong>l trabajo sobre<br />
<strong>la</strong> p<strong>la</strong>taforma propuesta, con el fin <strong>de</strong> incluir <strong>la</strong><br />
posibilidad <strong>de</strong> que el usuario ingrese <strong>la</strong> curva <strong>de</strong> <strong>la</strong><br />
bomba, a partir <strong>de</strong> datos experimentales obtenidos<br />
en p<strong>la</strong>nta o que provengan <strong>de</strong>l fabricante.<br />
Adicionalmente se van a implementar <strong>la</strong>s<br />
ecuaciones necesarias para que pueda hacerse <strong>la</strong><br />
predicción <strong>de</strong> <strong>la</strong> característica insta<strong>la</strong>da en válvu<strong>la</strong>s<br />
que manejen fluidos compresibles, esto es, gases y<br />
vapores.<br />
variabilidad <strong>de</strong>l proceso. In: Ingeniería Química.<br />
Madrid. Vol. 31, No. 359; p. 127-134<br />
4. Siler, B. y Spotts, J. (1999). Visual Basic 6.<br />
Madrid: Prentice Hall.<br />
5. Smith, C. y Corripio, A. (1997). Principles and<br />
practice of automatic process control, 2 nd ed. New<br />
York: John Wiley & Sons, Inc.<br />
Bibliografía<br />
1. Anaya, A. (1995). Dimensionamiento <strong>de</strong> válvu<strong>la</strong>s<br />
<strong>de</strong> control. Parte II. In: Ingeniería Química. Madrid.<br />
Vol. 27, No. 309; p. 169-173.<br />
2. Ceballos, F. (1997). Enciclopedia <strong>de</strong> Visual Basic<br />
4. México: Alfaomega.<br />
3. Metcalf, J. (1999). Válvu<strong>la</strong>s <strong>de</strong> control y