La transformación canónica de una ecuación de

segundo grado como método de análisis de

superficies de respuesta
José R. Guerrero Haber*, Juan M. Mesa Pérez**, Luis E. Brossard Pérez*,
Eleazar Henry Izquierdo*, Luis Brossard González*
*Universidad de Oriente, **Centro Nacional de Electromagnetismo Aplicado (CNEA)

El artículo presenta una metodología para transformar una ecuación de segundo orden a su forma
canónica. Se expone un esquema generalizado de aprendizaje (EGA) para su transformación, así como
la interpretación de la ecuación canónica obtenida.
A manera de ejemplo se muestra la utilización de la metodología en la optimización de la etapa de
extracción básica de alginato de sodio, profundizando en la interpretación del modelo matemático y la
superficie de respuesta.
Palabras claves: ecuación de 2do grado, ecuación canónica.
_____________________

In this paper a methodology for the transformation of a second degree equation to its canonical
expression is presented. A generalized chart for knowledge (EGA) used for this transformation, as well
as the interpretation of the presented canonical equation are exposed.
As an example, the utilization of this methodology in the optimization of basic leaching step of sodium
alginate is showed, making emphasis in the mathematical model and response surface interpretation.
Key words: 2nd degree equation, canonical equation.

Introducción Esta expresión es la llamada ecuación canónica

de la ecuación (1). La interpretación analítica de la
Muchos autores /1, 2, 3, ... 7/ han definido respuesta depende de los valores y signos de los
estrategias para la búsqueda de la mejor respuesta coeficientes J1, J2 y J3 y se realiza a partir de
(máximo o mínimo) a partir de un modelo matemá- conceptos geométricos.
tico. Para tres variables independientes un polinomio La forma de una ecuación cuadrática con dos
de segundo orden presenta la siguiente forma: variables independientes es la siguiente:

y = b0 + b1 x1 + b2 x2 + b3 x 3 + b12 x1 x2 + b13 x1 x 3 + y = b0 + b1 x1 + b2 x 2 + b12 x1 x 2 + b11 x12 + b22 x 22 (3)

+ b23 x 2 x3 + b11 x + b22 x + b33 x
2
1
2
2
2
3
(1)
Esta ecuación puede representarse en el plano
Para este caso es difícil definir el efecto de cada por cualesquiera de los cuatro contornos represen-
factor por simple inspección del modelo matemáti- tados en la figura 1. A partir de los valores de los
co. Para dos factores la representación de la res- coeficientes b0, b1, b2, b12, b11 y b22 no es posible
puesta en un plano brinda la información necesaria, dilucidar a cuál contorno pertenece la ecuación.
sin embargo, para tres o más factores su representa- La transformación canónica de la ecuación (3)
ción resulta imposible. Haciendo transformaciones presenta la siguiente forma:
al modelo se puede llegar a una expresión que
resuelve esta problemática, para este caso el modelo y = y c = J 1 X 12 + J 2 X 22 (4)
transformado será el siguiente:
¿Cómo se interpreta la ecuación en función de
y = yc = J1 X + J 2 X + J 3 X
2
1
2
2
2
3
(2) conceptos geométricos?:

TECNOLOGÍA QUÍMICA Vol. 19, No. 1, 1999

83
 Fig. 1 Diferentes formas geométricas de la superficie de respuesta para dos variables.

84 TECNOLOGÍA QUÍMICA Vol. 19, No. 1, 1999

1. Si J1 y J2 son negativos (positivos), y el centro mo) sea el realmente hallado en el punto C; se
del sistema (x1c, x2c) está dentro de la zona aconseja realizar experimentos adicionales que
estudiada, se puede decir que un alejamiento del lo conduzcan hacia la mejor respuesta y luego
centro "C" en cualquier dirección representa en la nueva zona hallar otra correlación de 2do
una pérdida (ganancia) cuadrática de la res- orden y realizarle el mismo análisis.
puesta; por lo tanto, la superficie ajustada tiene 3. Si J1 es negativo y J2 positivo (o viceversa),
un máximo (mínimo) en "C" [(Fig. 1a )]. Pueden estaremos frente a un minimax [(fig. 1 c)]; se
aparecer dos casos : espera que la respuesta se incremente en los
1.1.|J1| ≈ |J2| : Da como resultado circunferencias sentidos donde Ji sea positivo y lo hará mas
concéntricas, por lo que la variación en cual- rápidamente a medida que Ji sea mayor.
quier dirección de x1 o x2, produce la misma 4. Si J1 o J2 es cero o aproximadamente igual a cero,
pérdida (ganancia) en la respuesta. estaremos en presencia de una teja horizontal
1.2.|J1| < |J2| : lo que nos indica que la variación en donde no hay centro definido, sino que hay una
la dirección de x2, produce mayor pérdida (ga- recta de máximos para la misma [(fig. 1 d)].
nancia) en la respuesta que en la dirección x 1. En
este caso serían elipses concéntricas. Desarrollo
2. Si J1 y J2 son negativos (positivos), y el centro
del centro del sistema (x1c, x2c) está fuera de la Para llevar una ecuación de segundo orden a su
región estudiada [(fig. 1 b)], la superficie será forma canónica se sugieren los pasos planteados en
la llamada loma ascendente (descendente), de la el siguiente esquema generalizado de aprendizaje
que no podemos asegurar que el máximo (míni- (EGA) :

Fig. 2 EGA para la transformación de una ecuación cuadrática en la forma canónica.

TECNOLOGÍA QUÍMICA Vol. 19, No. 1, 1999

85
I. Expresión de la ecuación de 2do orden II Determinación de las coordenadas del
en forma matricial: centro: [XC]
Ecuación cuadrática para tres factores: El primer paso en la búsqueda de las coordena-
das del centro es comprobar que existe la matriz
y = b0 + b1 x1 + b2 x 2 + b3 x 3 + b12 x1 x2 + b13 x1 x 3 + inversa B 11-1 *. Esta prueba se reduce a calcular si el
+ b23 x2 x 3 + b11 x12 + b22 x 22 + b33 x 32 (5) determinante de la matriz B11 es distinto de cero.

se expresa en su forma matricial: det B11 ≠ 0 (7)

y = b0 + B1 X + X ′B11 X (6) Si el determinante es cero la superficie no tiene

centro, constituyendo lo que se conoce como super-
donde : ficie de segundo grado degenerada.
B1 :es la matriz fila compuesta por los coeficientes Luego de comprobar la existencia de B11-1 se
que acompañan los términos de primer orden. calculan las coordenadas del centro [XC], igualando
a cero todas las derivadas parciales con respecto a
cada xi o resolviendo el sistema de ecuaciones
B1 = b1 b2 b3
resultante que expresado en forma matricial:

X: es la matriz columna que representa las variables

independientes. XC = − b g⋅ B
1
2 11
−1
⋅ B1 ' (8)

Lx O donde el subíndice "c" nos indica que se trata de las

X = MM x PP
1
coordenadas del centro o punto estacionario.
MN x PQ
2
A partir de las coordenadas del centro XC
3
se calcula la respuesta predicha en el punto
estacionario (yc) sustituyendo X=XC en la ecuación
B11 : es una matriz cuadrada y simétrica que está (6):
compuesta por los coeficientes que acompañan
los términos de 2do orden (incluyendo las y C = b0 + B1 X C + X C ' B11 X C (9)
interacciones).
- para i ≠ j:
Con el centro de la superficie ubicado en un
bij = b ji = 12 bij
máximo, un mínimo o un punto estacionario, el
próximo paso será:
donde:
bij ; b ji → elementos de la matriz; III Trasladar el origen de coordenadas al
bij → coeficientes que acompañan las interac- centro del sistema (superficie)
ciones de 2do orden en la ecuación (1).
El objetivo de este paso es eliminar los términos
- para i = j: (diagonal principal de la matriz B11) de primer orden de la ecuación de 2do grado que
bii = bii define la superficie.
El nuevo valor de X& i referido a las coordenadas
LMb11
b12 b13OP del centro X& iC se determina como:
B11 = Mb b22 b P
MNb b PQ
21 23

31 b32 33
x&i = xi − xiC → xi = x&i + xiC

* La condición necesaria y suficiente para que exista la inversa de la matriz "A" es que sea NO SINGULAR. Una matriz es no singular
cuando su determinante es desigual de cero.

86 TECNOLOGÍA QUÍMICA Vol. 19, No. 1, 1999

y matricialmente se representa de la siguiente ma- entre las variables en el modelo cuadrático (ec. 15),
nera: o sea, para eliminar los bij (donde i ≠ j).
¿En qué consiste el problema ahora?
X = X& + X C (10) Consiste en hallar una matriz real única, la
matriz U, de forma tal que el cambio de variable:
Sustituyendo la ecuación (10) en la (6) se obtie- *
ne que: X =U X (16)

d i d i d
yC = b0 + B1 X& + X C + X& + X C B11 X& + X C (11) i al ser introducido en la ecuación (15) produzca el
resultado siguiente:

Luego restando (9) a (11): d i d i

y = yC = UX B11 XU
|

y − y C = B1 X& + X& ' B11 X& + X& ' B11 X C + X C ' B11 X&
d i
sabiendo que UX = X |U |
|

d i
y − y C = B1 X& + X& ' B11 X& C + X C ' B11 X& X& ' B11 X& (12) c
y = yC = X | U | B11U X h
De la ecuación (8): de la cual queda:

y − y C = J1 X 12 + J 2 X 22 + J 3 X 32 (17)
bg
B11 X C = 12 ⋅ B1 '
bg
X& ' B11 X C = 21 ⋅ X& ' B1 ' (13) que es la ecuación canónica correspondiente a la
ecuación (5) donde los términos Ji fueron designa-
dos como parámetros de transformación, y no son
Como la matriz B11 es simétrica, se cumple que: más que los valores propios (eigenvalues) de la
matriz (B11), determinados a partir de la expresión
X C ' B11 X& = X& ' B11 X C (14) matricial:

Usando las ecuaciones (13) y (14) y conociendo

que:
b
det B11 − J ⋅ I = 0 g (18)

. la cual puede expandirse en un polinomio como

B1 X& = X& ' B1 ' (porque B 1 y X son vectores fila y sigue:
columna respectivamente)

quedará una ecuación sin términos de primer orden:

bg
P J = J n + P1 J n−1 +...+ Pn−1 J + Pn = 0 (19)

El escalar "J" es una de las n raíces de la ecua-

y − y C = X& ' B11 X& (15)
ción; cada una de estas raíces se conoce como un
valor propio. Por lo tanto, existirán "n" valores de Ji
En esta ecuación los ejes Xi presentan la misma
y asociados a cada uno existirá un vector columna Ui
dirección que los ejes originales Xi.
tal que satisfaga la ecuación:

IV. Rotación de los ejes para que la superficie

sea simétrica con respecto a los nuevos ejes X i
bB 11 g
− J1 I U i = 0 (19.1)

designados como ejes principales

Estos vectores se llaman vectores propios de la
El objetivo de hacer la rotación de los ejes es matriz B11 y uniendo los vectores propios Ui se
anular los términos que representan las interacciones forma la matriz U.

TECNOLOGÍA QUÍMICA Vol. 19, No. 1, 1999

87
 Luego, a partir de las ecuaciones (10) y (11) se cesos que requieren agentes suspensor, gelificantes
obtiene: o emulsionante en ciertos compuestos.
La viscosidad de una solución de alginato al 1 %
b
X = U −1 X − X C g (20) constituye uno de los parámetros de calidad más
importante usado en las distintas industrias. Desde
Para tres variables y sabiendo que U-1 = U’, el punto de vista económico es de interés extraer
queda: más alginato con la menor cantidad de algas posible
(altos rendimientos).
b g b g b g
X 1 = u11 x1 − x1 C + u21 x 2 − x2 C + u31 x 3 − x 3 C Se han definido dos variables de respuesta, la
X2 = u bx − x g + u bx − x g + u bx − x g
12 1 1C 22 2 2C 32 3 3C (21) viscosidad de la solución de alginato al 1 % y el
= u bx − x g + u bx − x g + u bx − x g
rendimiento en la extracción. El ejemplo que se
X3 13 1 1C 23 2 2C 33 3 3C presenta se dedicará a estudiar la más importante,
que en este caso resulta ser la viscosidad, para la
Si se quiere explorar una superficie a lo largo de
cual se ajustó la siguiente ecuación de 2do orden
uno de los ejes, por ejemplo, el eje X1, se hace el otro
/7/, y cuya representación gráfica aparece en la
igual a cero y se toman los valores de x1 y x2 que
fig. 1 a).
satisfagan esta restricción, valiéndonos de las
ecuaciones generales para los ejes originales des-
pués de trasladado su origen al centro de la superfi- y = 1 444 − 60 x1 + 59 x 2 + 79 x1 x2 − 362 x12 − 324 x22 (23)
cie, en función de los ejes principales [de la ecua-
ción (16)] : la cual se expresa en su forma matricial (paso I):

x1 − x1C = U 11 X 1 + U 12 X 2 + U 13 X 3 y = b0 + B1 X + X ' B11 X (24)

x2 − x2 C = U 21 X 1 + U 22 X 2 + U 23 X 3 (22)
x3 − x 3C = U 31 X 1 + U 32 X 2 + U 33 X 3 B1 = [ −60 + 59]

X=
LM x OP
1

Estas expresiones resultan más cómodas para Nx Q

2
hacer exploraciones a lo largo de algunos de los ejes
principales. Basta hacer X2=0, por ejemplo, y se B11 =M
L−362 39,5OP = LM b 11 b12
PQO
pueden calcular x2 y x1 que satisfagan esta restric- N 39,5 −324 Q N b 21
b 22
ción.
coordenadas del centro y la respuesta predicha en el
Ejemplo de obtención de la ecuación canónica mismo (Paso II)*:
a partir de una ecuación de 2do orden
δy
= −60 + 79 x2 − 724 x1 = 0
El alginato de sodio es un polímero que se δy1
obtiene a partir de las algas pardas. En el caso que (25)
δy
se estudia estarán representadas por los sargazos. = −59 + 79 x1 − 648 x2 = 0
δy 2
El alginato de sodio ofrece propiedades particu-
lares que permiten ser utilizado en diversas ramas Del sistema (25) se obtiene la solución siguien-
de la industria, fundamentalmente, en aquellos pro- te:

De aquí se infiere que el centro de la superficie se

x1C = −0,074
→ encuentra dentro de la región estudiada (-1,1) para
x2 C = 0,081
ambas variables.
y C = 1 448,64

* Por ser más sencillo en este caso, se hallarán las derivadas parciales y se igualan a cero, para encontrar el punto estacionario buscado.

88 TECNOLOGÍA QUÍMICA Vol. 19, No. 1, 1999

 Trasladando el origen de coordenadas al centro
det
LM−362 − J 39,5 OP = 0 (26)
de la superficie (Paso III):
N 39,5 −324 − J Q
*
*
x = x1 + 0,074
x i = xi − xi → *
donde queda:
x = x 2 − 0,081
C

Determinación de los valores propios (27)

J 2 + 686 J + 115 726,75 = 0
(eigenvalues) de la matriz (B11) (Paso IV):

b
det B11 − J ⋅ I = 0 g Resolviendo la ecuación 27 se obtiene:

Como J 1 y J 2 son negativas, en el centro de

J 1 = −299 ,17 la superficie se encuentra un máximo, y como
→
J 2 = −386,83 J 1 ≠ J 2 la superficie está constituida por
elipses concéntricas.

A partir de los valores propios y utilizando la Resolviendo 19 se obtiene el sistema de

ecuación (19.1), se determinan los vectores propios ecuaciones siguientes:
de la matriz B11:
24,83u1 + 39,50u2 = 0 (29.1)
Para J=-299,17:
39,50u1 + 62,83u2 = 0 (29.2)
LM−362 + 299,17 39,5 OP LMu OP = 0
1
(28)
N 39,5 −324 + 299,17Q Nu Q 2
Tomando u1=1 y sustituyendo se halla u2:
donde se llega al siguiente sistema de ecuaciones:
de (29.1): u2 = -0,628 6
de (29.2): u2 = -0,628 7 → u2=-0,628 7
−62,83u1 + 39,50u2 = 0 (28.1)

−39,50u1 + 24 ,83u2 = 0 (28.2)

l= b1g + b0,628 7g
2 2
= 1181
, 2
Tomando u1=1 y sustituyendo se halla u2:
de (28.1): u2 = 1,591 1
→ u1 = = 0,847
de (28.2): u2 = 1,591 u2=1,591
1181
, 2
→ U2 =
LM 0,847 OP
l= b1g + b1,519g
2 2
= 1,879 17 u2 =
1,591
1181
, 2
= 0,532 N−0,532Q
1
u1 =
1,879 17
= 0,532
→ U1 =
LM0,532OP LM0,532 0,847OP
u2 =
1,591
1,879 17
= 0,847 N0,847Q U =
N0,847 −0,532Q

Para J=-386,83: Finalmente, la ecuación canónica será:

LM−362 + 386,83 39,5 OP LMu OP = 0
1
(29) y − 1 448,64 = −299,17 X 12 − 386,83 X 12
N 39,5 −324 + 386,83Q Nu Q 2

TECNOLOGÍA QUÍMICA Vol. 19, No. 1, 1999

89
 Si se desea conocer el valor de las viejas coorde- de donde:
nadas donde se obtiene una respuesta determinada,
conociendo las coordenadas de los ejes canónicos o x1 + 0,074 = 0,532 X 1 + 0,847 X 2 (31.1)
viceversa se puede emplear el procedimiento analí-
tico que a continuación se ilustra: x2 − 0,081 = 0,847 X 1 − 0,532 X 2 (31.2)

Como: que son las ecuaciones que definen una superficie

*
elíptica.
X =U⋅X Estas ecuaciones brindan la posibilidad de trans-
formar el sistema de coordenadas original en el
se puede definir: descrito por la ecuación canónica y viceversa.

LM X OP = L0,532 0,847 OPLM x + 0,074OP ¿Cuándo es necesario realizar una

N X Q MN0,847
1 1
(30)
−0,532Q N x − 0,081Q
transformación canónica?
2 2
Hasta dos variables independientes se puede
prescindir de este método, ya que el análisis gráfico
b g
X 1 = 0,532 x1 + 0,074 + 0,847 x 2 − 0,081 b g (30.1) es posible, pero: ¿Cómo analizar gráficamente 3
variables o más?
Muchos investigadores resuelven este problema
b g
X 2 = 0,847 x1 + 0,074 − 0,532 x 2 − 0,081 b g (30.2) dejando 2 variables libres y "fijan" las restantes a
determinado nivel (generalmente -1, 0, ó +1) aten-
diendo a razones generalmente económicas o de
y además, se puede decir que: diseño. Tomando como ejemplo lo planteado por
* Mesa y colaboradores en su artículo /7/; donde se
X = U −1 X , y como: U -1 = U ' optimizó la etapa de extracción de alginato de sodio,
y en la cual se identificaron como variables de
luego, posible influencia X1: temperatura, X2: concentra-
LM x + 0,074 OP = LM0,532
1
0,847 OPLM X OP
1
(31)
ción y X3: relación líquido/sólido y como respuesta
y: la viscosidad de una solución de alginato de sodio
N x = −0,081Q N0,847
2 −0,532Q N X Q2 al 1 %. La ecuación de 2do grado obtenida fue:

y = 1 763,58 − 121,9 X 1 + 85,08 X 1 − 95,03 X 1 − 31,87 X 1 X 2 + 115,87 X 1 X 3 − 90,87 X 2 X 3 −

−496,55 X 12 − 445,06 X 22 − 484,35 X 32 (32)

Las coordenadas y la respuesta predicha en el Luego se fijó x3 en el nivel inferior (-1) y

centro son: basándose en ello se analizó la superficie de
respuesta.
x1c = −0,1411
¿Qué dificultades se presentan cuando el proble-
x2 c = 0,11 3 5 ma se ataca de esta manera?
x3c = −0,126 2 Según la tabla 1 y las superficies de respuesta
y e = 1 783,04 generadas por los modelos, para cada valor de x3 se
Tabla 1
Optimos locales condicionados por X3
Valor en que se fija Coordenadas del centro Respuesta predicha Observaciones
X3 X1c X2c
+1 -0,005 9 -0,006 3 1 279,27 fig. 3 a)
0 -0,126 0 0,100 1 1 775,52 fig. 3 b)
-1 -0,246 1 0,206 5 1 469,61 fig. 3 c)

90 TECNOLOGÍA QUÍMICA Vol. 19, No. 1, 1999

 Fig.Fig.
3 a)
3 a)
1
1500
0,6
1200
0,2
900 X2
Y
-0,2
600

300 -0,6
1
0 0.20.6
-0.2 -1
-1 -0.6
-1 -0.6 -0.2 0.2 0.6 X2 -1 -0,6 -0,2 0,2 0,6 1
1
X1 X1

Fig.Fig.
3 b)3 b)
1
1800
1600 0,6
1400 0,2
Y 1200 X2
-0,2
1000
800 -0,6
1
600 0,20,6
-0,2
-1 -0,6 -0,2 0,2 0,6 -1-0,6 X2 -1
1 -1 -0,6 -0,2 0,2 0,6 1
X1 X1
Fig. 3 Representación de las superficies de respuesta de la ecuación 32 con x3 fijo a determinado nivel.

obtiene un máximo condicionado para la viscosi- cualquier dirección que me desplace habrá una
dad, y para cada condición habrá también una caída en la respuesta.
estrategia para la búsqueda de mejores condiciones. 2. Como J 3 es mayor (modularmente), es la que más
Esta situación conduce a una pérdida de tiempo y de influencia ejerce sobre la respuesta y por lo
información sobre el sistema que se estudia, situa- tanto, sobre la que se debe tener mejor control
ción que puede ser resuelta haciendo el análisis de para evitar su variación.
la superficie de respuesta a partir de los modelos 3. Los valores de x1, x2 y x3 deben tratar de mante-
transformados a su forma canónica. nerse cerca del centro de la superficie para
Transformando la ec. 32 a su forma canónica por obtener la mejor respuesta y que su rango de
el procedimiento descrito en el desarrollo se obtie- variación sea mínimo.
ne:

Conclusiones
y − y c = −392,78 X 12 − 480,69 X 22 − 542,49 X 32 (33)
Se ha demostrado que el análisis de superficies
A partir de este modelo se puede llegar a las de respuesta mediante el uso de la ecuación canóni-
siguientes conslusiones: ca de una ecuación de 2do grado es importante e
1. Como J1, J2 y J3 son negativos, el punto central imprescindible cuando existen más de dos varia-
hallado está ubicado sobre el óptimo y hacia bles, ya que permite:

TECNOLOGÍA QUÍMICA Vol. 19, No. 1, 1999

91
1. Analizar el comportamiento de la respuesta a 3. COCHRAN, W.G. y G.M. COX: Diseños experi-
partir de los coeficientes de la ecuación canónica mentales, Ed. F. Trillas S.A., México D.F.,
desde términos geométricos. 1965.
4. GEORGE E.P. BOX, W. HUNTER, G. HUNTER:
2. Determinar el rango de variabilidad con que Introducción al diseño de experimentos, análisis
debe trabajarse los factores para obtener una de datos y construcción de modelos, Editorial
respuesta cercana al deseado. Reverté S.A., España, 1993.
3. Definir qué factor debe ser controlado más es- 5. BOIZÁN JÚSTIZ, M. : Optimización, Ed. Pueblo
trictamente y a cuál con menos precisión. y Educación, Cuba,
6. LÓPEZ PLANES, REINALDO: Diseño estadísti-
co de experimentos, Editorial Científico - Técni-
Bibliografía ca, La Habana , Cuba 1988
7. A R I A S D I O B E R , J . C . B E T A N C O U R T :
1. MESA PEREZ, J.M., L. E. BROSSARD PÉREZ, Optimización del proceso tecnológico de obten-
J. R. GUERRERO HABER, E. HENRY IZQUIER- ción de alginato de sodio para la determinación de
DO: "Estrategia de utilización del diseño de expe- los parámetros de diseño, Trabajo de Diploma,
rimentos". Tecnología Química, 18(3):101-111, Universidad de Oriente, Santiago de Cuba, 1996.
Santiago de Cuba, 1998. 8. MESA PÉREZ J.M, M. VALLE y J. R. GUERRE-
2. BARROS NETO B, I. SPACINO, R. EDWARD: RO: "Optimización de la etapa de extracción bási-
Planejamento E. Optimicao de Experimentos, ca de alginato de sodio", Tecnología Química,
Campinas, SP.Editora de UNICAMP, Brasil, 1995. 18(2):90-96, Santiago de Cuba, 1998.

UNIVERSIDAD DE ORIENTE
Facultad de Ingeniería Química

Cursos de superación profesional

Título: Flujo de fluidos no-newtonianos

Contenido: Estudio de los métodos de cálculo hidráulico en flujo de fluidos no-

newtonianos en régimen laminar y turbulento. Selección de conductos óptimos. Flujo de
suspensiones homogeneas y no homogeneas.

Profesor: Ing. Armando Díaz García, PT

Fecha: Mayo 1999

Total de horas: 50
Centro solicitante: Oferta

92 TECNOLOGÍA QUÍMICA Vol. 19, No. 1, 1999