UNIDAD DOCENTE 1

Introducción y cinemática

Tema 3

Flujo convectivo y ecuación de la

continuidad.

3.1 Sistemas abiertos y sistemas cerrados.

Las ecuaciones que gobiernan el movimiento de un fluido se obtienen por aplicación
de las leyes de la Mecánica y la Termodinámica a algún volumen fluido, que por su
definición no permite el paso de masa pero si el de energía (los llamados sistemas
cerrados).
Sin embargo, con frecuencia resulta más útil aplicar las leyes a sistemas abiertos, y
por tanto a volúmenes que no son fluidos, esto es: volúmenes que permiten el paso de
masa a través de alguna o todas sus paredes, y que llamaremos volúmenes de control.
Estos se definen como se considere conveniente y sin restricciones, una vez establecido
el sistema de referencia en el que se estudia el flujo.
Necesitamos tener para ello dos conceptos claros que veremos en los apartados
siguientes:
- Flujo convectivo de una magnitud fluida.
- Derivación de magnitudes integrales extendidas a volúmenes fluidos.
Estos conceptos nos permitirá generalizar las leyes obtenidas originalmente para
sistemas cerrados (volúmenes fluidos) para que sean aplicables a sistemas abiertos
(volúmenes de control).

- Comprueba que has comprendido la relación sistema abierto-volumen de control,

sistema cerrado-volumen fluido. Pon ejemplos de unos y otros hasta entender mejor
ambos tipos de análisis.

UD1 – Tema 3. Flujo convectivo y ecuación de la continuidad.

1
 Mecánica de Fluidos.
2º curso del Grado en Ingeniería en Tecnologías Industriales.

3.2 Flujo convectivo de una magnitud fluida.

El flujo convectivo de una magnitud fluida extensiva (siendo φ su valor intensivo por
unidad de volumen) a través de una superficie Σ fija al sistema de referencia, es la
cantidad de esa magnitud que atraviesa con el fluido dicha superficie en la unidad de
tiempo.
Si φ es la magnitud fluida por unidad de
 
volumen, (φ dσ )( v ⋅ n )dt es la cantidad de φ
que atraviesa la superficie dσ durante el
tiempo dt. Por unidad de tiempo será:
 
 
v ⋅ dt φ ( v ⋅ n)dσ , y para una superficie finita es la
n integral:
dσ
∫
 
Flujo de φ = φ(v ⋅ n )dσ
Σ

Si φ es un escalar, la cantidad φ v se llama vector flujo de φ.

Si φ es un vector, la cantidad φ v se llama tensor flujo de φ.
 
Así por ejemplo, el flujo volumétrico es v , el másico ρv y el de cantidad de

movimiento ρv v (un tensor de componentes ρvi v j ).

Si la superficie Σ es cerrada y φ v es continuo en el 
v
interior y en la superficie, el teorema de Gauss-Ostrogradsky
nos permite escribir,
Σ
∫ ∫
  
Flujo de φ = φ (v ⋅ n )dσ = ∇ ⋅ (φ v )dV , (3.1)
Σ V

 n
con lo cual ∇ ⋅ (φ v ) , divergencia del vector flujo, es el
resultado neto del flujo de φ hacia el exterior en cada diferencial de volumen, y por
integración se obtiene el flujo hacia el exterior de todo el volumen de la magnitud φ.
 
Si la superficie Σ se mueve con una velocidad normal de avance vn tal que vn = vc ⋅ n

, donde vc es la velocidad del punto de la superficie que estamos considerando, el flujo
de φ está dado por

∫
  
Flujo de φ = φ (v -vc ) ⋅ n dσ . (3.2)
Σ

De esta expresión se pueden obtener dos conclusiones inmediatas:


- El flujo de φ a través de una superficie de corriente es nulo puesto que vc = 0 y
 
v ⋅n = 0.
 
- El flujo de φ a través de una superficie fluida es también nulo puesto que v = vc .

UD1 – Tema 3. Flujo convectivo y ecuación de la continuidad.

2
 Mecánica de Fluidos.
2º curso del Grado en Ingeniería en Tecnologías Industriales.

- Asegúrate que entiendes bien el significado de la nomenclatura utilizada en las

expresiones presentadas en el apartado.
- ¿El concepto definido para flujo convectivo coincide con el concepto intuitivo que
podías tener?.
- Distingue claramente la velocidad de la superficie en un punto de la misma, de la
velocidad del fluido en ese punto de la superficie en la ecuación (3.2). Busca ejemplos
de situaciones en que no coincida.

3.3 Derivación de integrales extendidas a volúmenes fluidos. Teorema del

transporte de Reynolds.

La derivada sustancial nos permite calcular la variación de una magnitud fluida

intensiva siguiendo a una partícula fluida. Cuando lo que se pretende es seguir la
evolución de una magnitud determinada integrada en un volumen fluido, pero el análisis
se hace en un volumen de control que no lo es, es necesario establecer la relación entre
la derivada de la magnitud fluida en uno y otro.
Sea φ una magnitud fluida extensiva por unidad de volumen ligada al fluido. La
cantidad de esa magnitud asociada a un volumen de control cualquiera Vc (t) variable
con el tiempo está dado por,

∫ φ( x, t )dV

Vc ( t )

y la variación en la unidad de tiempo de esa cantidad es,

 
∫
d   
φ ( x, t )dV .
dt  
Vc ( t ) 
Para determinar esta derivada escribiremos,

∫ ∫
  
φ ( x, t + ∆t )dV − φ ( x, t )dV 
   
∫
d     Vc ( t + ∆t )  ,
φ ( x, t )dV = lím Vc ( t )

dt   ∆t →0  ∆t 
Vc ( t )   
 
donde el volumen Vc (t + ∆t ) = Vc (t ) + δV2 − δV1 (véase figura adjunta). Por lo tanto
escribiremos,

∫ φ( x, t + ∆t)dV − ∫ φ( x, t)dV = ∫ [φ( x, t + ∆t) − φ( x, t)]dV +

   
Vc ( t + ∆t ) Vc ( t ) Vc ( t )
(3.3)
∫ ∫
 
+ φ ( x, t + ∆t )dV − φ ( x, t + ∆t )dV
δV2 δV1

UD1 – Tema 3. Flujo convectivo y ecuación de la continuidad.

3
 Mecánica de Fluidos.
2º curso del Grado en Ingeniería en Tecnologías Industriales.

dado que para los volúmenes infinitesimales δ V1

Vc(t+dt)
δV2 y δ V2 , el elemento diferencial de volumen dV se
puede escribir como
 
vnΔt dV = ± v n ∆t ⋅ dσ = ±v ⋅ n ⋅ ∆t ⋅ dσ ,

δV1 donde el signo "+" es para δ V2 y el "-" para δ V1

(debido a que las normales están cambiadas de
signo) y por lo tanto estas integrales de volumen
se convierten en integrales de superficie. A su vez,
Vc(t) -vnΔt 
desarrollando φ ( x, t + ∆t ) en serie de Taylor
∂φ
∫ ∫ ∫ ∫
  
φ ( x, t + ∆t )dV = φ ( x, t + ∆t ) v n dσ ∆t = ∆t φ ( x, t )v n dσ + (∆t )2 v dσ (3.4)
∂t n
δV Σ Σ Σ

y despreciando términos del orden ( ∆t ) 2 frente a ∆t cuando éste es pequeño, podemos

escribir,
 
∫ ∫ ∫ ∫
   
φ ( x, t + ∆t )dV − φ ( x, t + ∆t )dV = ( ∆t )  φ ( x, t )v n dσ + φ ( x, t )v n dσ  =
 
δV2 δV1 Σ1 ( t ) Σ2 ( t )  (3.5)

∫

= (∆t ) φ ( x, t )v n dσ
Σc ( t )

donde Σ c (t ) es la superficie que encierra el volumen Vc ( t ) . Esta superficie se ha

descompuesto en dos partes, Σ c ( t ) = Σ1 ( t ) + Σ 2 ( t ) , la primera es por donde el volumen
disminuye y la segunda por donde aumenta. La velocidad vn es en este caso la
 
componente normal de la velocidad de avance de la superficie de control: vn = vc ⋅ n .
Sustituyendo esta ecuación en la (3.4), dividiendo por ∆t y tomando límites, podemos
concluir que,

∂ φ( x , t )
∫ ∫ ∫
   
dV + φ ( x, t ) vc ⋅ n dσ
d
φ ( x, t )dV = (3.6)
dt ∂ t
Vc ( t ) Vc ( t ) Σc ( t )

Supongamos ahora que el volumen Vc ( t ) fuese un volumen fluido Vf ( t ) , en ese

  
caso vn = v ⋅ n siendo v la velocidad del fluido. De acuerdo con la Ec.(3.6) tendremos,

∂ φ( x, t )
∫ ∫ ∫
d    
φ ( x, t )dV = dV + φ ( x, t )v ⋅ ndσ (3.7)
dt ∂ t
V f (t ) V f (t ) Σ f (t )

Si ahora nos fijamos en un volumen Vc ( t ) , que en el instante t = t 1 coincide con

Vf ( t ) , es decir Vc ( t 1 ) ≡ Vf ( t 1 ) , (y también Σ c ( t 1 ) ≡ Σ f ( t 1 ) ); las primeras integrales del
segundo miembro de las Ecs (3.6 y 3.7) coinciden puesto que los integrandos y los
recintos de integración de cada una de estas integrales coinciden en ese instante. Las
segundas integrales del segundo miembro difieren porque los integrandos son diferentes
   
ya que vc ⋅ n es distinto de v ⋅ n , pues son las componentes normales de dos velocidades

UD1 – Tema 3. Flujo convectivo y ecuación de la continuidad.

4
 Mecánica de Fluidos.
2º curso del Grado en Ingeniería en Tecnologías Industriales.

distintas: la primera de la velocidad de la superficie de control, y la segunda de la

velocidad del fluido. Restando a la Ec.(3.7) la Ec.(3.6) se obtiene,

∫ ∫ ∫
d  d     
φ ( x, t )dV = φ ( x, t )dV + φ ( x, t )(v -vc ) ⋅ ndσ (3.8)
dt dt
V f (t ) Vc ( t ) Σc ( t )

La Ec.(3.8) nos permite calcular la variación en la unidad de tiempo de una magnitud

fluida extensiva ligada a un volumen fluido como la variación de esa magnitud fluida
dentro de un volumen de control más su flujo convectivo a través de la superficie que
limita dicho volumen de control.
- Estudia la demostración hasta poder reproducirla entendiendo cada paso y el objetivo
general del teorema. ¿Está clara toda la nomenclatura utilizada?.

3.4 Principio de conservación de la masa. Ecuación de continuidad en forma

integral y diferencial.
El requerimiento de conservación de la masa de un fluido impone ciertas
restricciones al campo de velocidades, que por tanto no puede ser cualquiera. Estas
restricciones cinemáticas estarán ligadas a la distribución de masa, y por tanto al campo
de densidades.
El principio de conservación de la masa nos dice que la masa de un volumen fluido
(sistema cerrado) no cambia con el tiempo, esto es

∫
d 
ρ ( x , t )dV = 0 . (3.9)
dt
V f (t )

Mediante la ecuación del transporte de Reynolds podemos generalizar la ecuación

anterior para volúmenes de control, es decir para sistemas abiertos,

∫ ∫
d    
ρ ( x , t )dV + ρ (v − vc ) ⋅ n dσ = 0, (3.10)
dt
Vc ( t ) Σc ( t )

indicándonos que la variación en la unidad de tiempo de la masa contenida en un

volumen de control se debe al flujo (convectivo) de masa a través de sus paredes. Ésta
es la ecuación de la continuidad en forma integral.
Si el volumen de control es fijo respecto al sistema de referencia, la Ec.(3.10) puede
escribirse en la forma,
∂ρ
∫ ∂ t dV + ∫ ρ(v ⋅ n )dσ = 0.
 
(3.11)
Vo Σo

Si ρv es continua en el interior de Vo y sobre la superficie Σ o , el teorema de
Gauss-Ostrogradsky permite escribir,
 ∂ρ 
∫  ∂t + ∇ ⋅ (ρ v ) dV = 0.
Vo
(3.12)

UD1 – Tema 3. Flujo convectivo y ecuación de la continuidad.

5
 Mecánica de Fluidos.
2º curso del Grado en Ingeniería en Tecnologías Industriales.

dado que la integral es nula para cualquier volumen de control Vo fijo (y de paredes Σ o
fijas) al sistema de referencia, esto sólo puede ser posible si,
∂ρ 
+ ∇ ⋅ (ρ v ) = 0, (3.13)
∂t
que constituye la ecuación de conservación de la masa (o ecuación de la continuidad)
en forma diferencial. Formas alternativas de la Ec. (3.13) son,
∂ρ  
+ v ⋅ ∇ρ + ρ∇ ⋅ v = 0, (3.14)
∂t
1 Dρ 
+ ∇ ⋅ v = 0. (3.15)
ρ Dt
La Ec.(3.15) puede interpretarse en términos de cambio en el volumen de un
volumen fluido. En efecto, la variación con el tiempo de un volumen fluido se puede
escribir (con ayuda del teorema del transporte de Reynolds) como,

∫ ∫ ∫
dV f (t )   
d
dV ≡ = 0+ v ⋅ n dσ = (∇ ⋅ v )dV , (3.16)
dt dt
V f (t ) Σ f (t ) V f (t )

y en definitiva,

 1 dV f   
∫
 
  = 1 (∇ ⋅ v )dV  = ∇ ⋅ v . (3.17)
Lím V dt  Lím V
V f →0 
 f  V f →0  f 
Vf 

Por lo tanto ∇ ⋅ v representa la variación de la unidad de volumen en la unidad de
tiempo, por lo que se llama también velocidad de dilatación cúbica unitaria.

- ¿Este principio sigue siendo aplicable si tenemos en el volumen considerado puntos

singulares como manantiales y sumideros? Explica estas aparentes excepciones y
modifica el principio enunciado para tenerlos en cuenta.
- Estudia y reproduce el razonamiento que lleva del concepto general de conservación
de la masa, a la ecuación integral y de ésta a la ecuación diferencial.
- Intenta resolver los problemas 12 (diferencial) y 13 (integral).

3.5 Caso particular de la ecuación de la continuidad: fluidos incompresibles.

Para fluidos incompresibles, aquellos que su densidad permanece constante, la
ecuación de la continuidad en forma integral toma la forma,
  
+ ∫ (v − vc ) ⋅ n dσ = 0,
dVc
(3.18)
dt Σc (t )

UD1 – Tema 3. Flujo convectivo y ecuación de la continuidad.

6
 Mecánica de Fluidos.
2º curso del Grado en Ingeniería en Tecnologías Industriales.

donde la integral de superficie representa el flujo volumétrico‚ Q(t), a través de las

paredes del volumen de control. La ecuación de la continuidad para fluidos
incompresibles se reduce a,

∇ ⋅ v = 0. (3.19)

Si se cumple No se cumple

- Analiza si los campos de velocidades de los problemas 1, 2, y 14 podrían corresponder

a un fluido incompresible.
- ¿Como tendrían que ser las constantes a y b del problema 12 para que el fluido cuyo
movimiento describen sea un líquido perfecto?

3.6 Caso particular de la ecuación de la continuidad: movimientos estacionarios.

En el caso del movimiento estacionario de un fluido compresible la ecuación de la
continuidad se reduce a,

∇ ⋅ (ρ v ) = 0. (3.20)
En este caso, el flujo másico G a través de cualquier
sección de un tubo de corriente es constante. Para ver
esto se toma un volumen de control limitado por la Σ2
superficie del tubo de corriente Σ T y dos secciones Σ1 ΣT
transversales del tubo Σ1 y Σ 2 , y la ecuación de la
continuidad en forma integral nos proporciona,

∫ ∫ ∫
    
ρ v ⋅ n dσ = ∇ ⋅ (ρ v )dV = 0 = ρ v ⋅ n dσ , (3.21)
Σc Vc ΣT + Σ1 + Σ 2

pero por la propia definición de tubo de corriente se tiene,

 
∫ ρ v ⋅ n dσ = 0,
ΣT
(3.22)

 
puesto que en Σ T es v ⋅ n = 0 . Por lo tanto se tiene,

UD1 – Tema 3. Flujo convectivo y ecuación de la continuidad.

7
 Mecánica de Fluidos.
2º curso del Grado en Ingeniería en Tecnologías Industriales.

-G1+G2=0; G1=G2 ,
siendo G1 y G2 los flujos másicos a través de las secciones Σ1 y Σ 2 respectivamente (el
signo menos en G1 se debe a que la normal a Σ1 y la velocidad forman un ángulo
superior a 90º si el fluido entra por Σ1 ).
Podemos concluir por tanto que en el movimiento estacionario de un fluido
compresible, el flujo másico a través de cualquier sección transversal al tubo es
constante, si no hay singularidades en el recinto limitado por las paredes del tubo y las
secciones transversales consideradas.
Si el fluido fuese incompresible ( ρ o = ρ ) lo dicho para el flujo másico sería
también cierto para el volumétrico‚ Q = G ρ o , en cada instante, puesto que en este
último caso no sería necesaria la restricción de movimiento estacionario.
Como consecuencia de lo anterior, un tubo de corriente no puede terminar en el
interior de un fluido, ya que al ser G (ó Q) constante, al tender el área de la sección del
 
tubo de corriente a cero, sería a costa de que ρv (ó v ) tendiese a infinito. Por lo tanto
un tubo de corriente en el movimiento estacionario de un fluido compresible (o en el de
un fluido incompresible aun siendo no estacionario) o bien se cierra sobre sí mismo, o
bien termina en los límites del fluido o bien se extiende hasta el infinito.

- Busca ejemplos de tubos de corriente en líquidos y gases.

- Suponiendo condiciones uniformes en cada sección de estos tubos, ¿qué primeras
conclusiones podemos sacar si la sección transversal se va modificando a lo largo del
tubo?
- Intenta hacer el problema 10.

3.7 Movimientos planos y bidimensionales. Función de corriente

En el caso del movimiento de un fluido incompresible o del movimiento
estacionario de un fluido compresible, la ecuación de la continuidad se reduce a la
divergencia nula de un vector.

- Fluido incomprensible: ∇ ⋅ v = 0.

- Movimiento estacionario fluido comprensible: ∇ ⋅ (ρ v ) = 0.
Si además suponemos que el movimiento es bidimensional o axilsimétrico‚ la
divergencia se reduce a la suma de dos términos, en cuyo caso la ecuación de la
continuidad se puede hacer cumplir mediante una función escalar ψ de la que se
 
obtienen las componentes de v o ρ v por derivación.
El procedimiento de obtención de ψ lo haremos para el movimiento bidimensional
de un fluido incompresible, pero de modo análogo se podría hacer para los demás casos
en que la ecuación de la continuidad se pueda escribir como suma de dos derivadas
parciales.

UD1 – Tema 3. Flujo convectivo y ecuación de la continuidad.

8
 Mecánica de Fluidos.
2º curso del Grado en Ingeniería en Tecnologías Industriales.


En un movimiento bidimensional v = (u , v,0 ) donde u y v son independientes de la
coordenada z. La ecuación de la continuidad para un fluido incompresible se reduce a,
 ∂u ∂v
∇⋅v ≡ + = 0, (3.23)
∂x ∂y
de modo que podemos definir una función escalar ψ tal que,
∂ψ ∂ψ
u= ; v=− , (3.24)
∂y ∂x
Utilizando la función ψ, se cumple automáticamente la ecuación de la continuidad y
se reduce el número de incógnitas, de dos variables dependientes: u y v, a una sola
función escalar desconocida.
A la función ψ se le denomina función de corriente porque las líneas ψ=constante
son las líneas de corriente. En efecto,
∂ψ ∂ψ
dψ = dx + dy = −v dx + u dy, (3.25)
∂x ∂y
y para ψ=constante se tiene dψ=0, de modo que,
dx dy
dψ = 0 ⇒ − v dx + u dy = 0 ⇒ = , (3.26)
u v
que es la ecuación de las líneas de corriente.
El flujo volumétrico (por unidad de longitud E a lo largo del eje z) a través de una
línea que una dos puntos O y P, es independiente de la línea elegida para unir dichos
puntos y se tiene,

p p p p p

∫ (u n + v n y )ds = ∫ (u dy − v dx ) = ∫ (ψ dy + ψ dx ) = ∫ dψ = ψ
Qop  
E ∫
= v ⋅ nds =
o o
x
o o
y x
o
p − ψo. (3.27)

Por lo tanto, la diferencia de valores de la función ψ entre dos puntos es igual al

flujo volumétrico (por unidad de envergadura en la dirección Z) que atraviese una línea
cualquiera que una esos dos puntos.
y

ψP
dl
P
ψP > ψO
𝑑𝑑𝑙𝑙⃗ = (𝑑𝑑𝑑𝑑, 𝑑𝑑𝑑𝑑)
ds
v
O
n ψ
O

𝑛𝑛�⃗ = (𝑑𝑑𝑑𝑑, −𝑑𝑑𝑑𝑑)

x

UD1 – Tema 3. Flujo convectivo y ecuación de la continuidad.

9
 Mecánica de Fluidos.
2º curso del Grado en Ingeniería en Tecnologías Industriales.

Obsérvese‚ que en el movimiento bidimensional estacionario de un fluido

compresible, la función de corriente sería tal que,
∂ψ ∂ψ
ρu = ; ρv = − , (3.28)
∂y ∂x
p

con, ψ p − ψ o = ∫ ρ (u dy − v dx ), (3.29)
o

siendo ψ p − ψ o el flujo másico (por unidad de longitud) a través de una línea

cualquiera que una los puntos O y P.
En el caso más general de movimiento en coordenadas curvilíneas ortogonales se
tiene,
 1  ∂ (h2h3ρ v1 ) ∂ (h1h3ρ v2 ) ∂ (h1h2 ρ v3 )
∇ ⋅ (ρ v ) ≡ + +
h1h2h3  , (3.30)
∂x1 ∂x2 ∂x3 
y si el movimiento es bidimensional nada depende de x3, de modo que la función de
corriente será,
∂ψ ∂ψ
h2 h3 ρv1 = ; h1h3 ρv2 = − . (3.31)
∂x 2 ∂x1
Los factores de escala h1, h2 y h3 y las coordenadas x1 , x2 (y x3 ) en los casos más
corrientes se han visto en el tema anterior.

- Escribe la ecuación de la divergencia del campo de velocidades en cartesianas, polares

cilíndricas y esféricas, y deduce a la vista de estas expresiones la definición de función
de corriente en cada caso.
- Estudia, hasta ser capaz de repetir, la demostración de la relación entre la diferencia
entre el valor de la función de corriente en dos puntos, y el caudal o flujo másico a
través de una línea que los una.
- Intenta hacer los problemas 4 y 9.

Problemas recomendados para este tema: 4, 6, 8, 9, 10, 11, 12, 13 y 25.

Bibliografía básica recomendada del tema 3.
CRESPO, Lec. 8 y 9 (pp. 114-129)
FOX y McDONALD, 4.2 y 4.3 (pp. 109 a 124), 5.1 (pp. 216 a 229)
GERHART et al., 3.3.5 (pp. 153 a 159), 4.2 (pp. 171 a 190) y 5.2.1 (pp. 359 a 367)
SHAMES, 4.7 y 4.8 (pp. 120 a 127), 5.2 (pp. 137 a 141), 7.2 (pp. 238 a 240) y 12.6 a
12.8 (pp. 507 a 510)
WHITE, 3.1 a 3.3 (pp. 130 a 147), 4.2 (pp. 221 a 227) y 4.7 (pp. 243 a 251)

UD1 – Tema 3. Flujo convectivo y ecuación de la continuidad.

10
 Mecánica de Fluidos.
2º curso del Grado en Ingeniería en Tecnologías Industriales.

Otra bibliografía
BARRERO y PÉREZ-SABORID, 3.2 (pp. 44-45), 3.4 (pp.50 a 52), y 5.1.1 (pp. 84 y
85)
BATCHELOR, 2.1 y 2.2 (pp. 73 a 79)
CANDEL, Cap. IV (pp. 77 a 86)
ÇENGEL y CIMBALA, cap. 4.6 (pp 162-170), 5.2 (pp 189-196), 9.2 y 9.3 (pp 442-
463)
FERNÁNDEZ FERIA, cap. 5 (pp. 55-58), cap. 6.1 y 6.2 (pp. 61-65)
HERNÁNDEZ KRAHE, 2.8 a 2.10 (pp. II/21 a 26)
LANDAU y LIFSHITZ, 1 (pp. 1 a 3)
LIÑÁN et al., Lec. 2.8 a 2.9, Lec. 4
LÓPEZ PEÑA, cap. 5.4 (pp. 67-69), cap. 7.1 y 7.2 (pp. 91-95)
POTTER y WIGGERT, 4.3 y4.4 (pp. 118 a 128) y 5.2 (pp. 181 a 186)

UD1 – Tema 3. Flujo convectivo y ecuación de la continuidad.

11