Ecuación de La Continuidad para Presión Constante

Descargar como docx, pdf o txt
Descargar como docx, pdf o txt
Está en la página 1de 3

Ecuación de la continuidad para presión constante ( ρ ) y coeficiente de difusividad constante

( D¿¿ AB)¿
Demostrar que ara coordenadas cilíndricas:

∂C A ∂ CA 1 ∂ CA ∂C A 1 ∂ ∂ C A 1 ∂ 2 ∂2 C A
∂t (
+ vr +
∂r
+ vθ
r ∂θ
+ vz
∂z
=D AB
r ∂r
r
∂r )
+ 2 2+
r ∂θ ( ( )
∂Z
2
+RA )
También se puede escribir de la siguiente manera:

∂C A 2
+( v •∇ C A )=D AB ( ∇ C A ) + R A
∂t
(1)

A partir del análisis tensorial sabemos que el divergente de φ y su laplaciano son descritas con la
ecuación 2 y 3 respectivamente
3
1 ∂φ
∇ φ=∑ g
i=1 √ gii ∂ x i ⟨ i⟩
(2)
3
2 1 ∂ ∂φ
∇ φ=∑ ( g gij j )
i √
i=1 √g ∂ x ∂x
(3)
Y el tensor métrico para las coordenadas cilíndricas es representado por la ec. 4

1 0 0
ij
g=0
[ ]
0
1
( x 1)
0
2
[
1 0 0
0 g ij = 0 ( x1 )2 0

1
0 0 1 ]
(4)

Y g es el determinante de gij
2
g= ( x1 )

1. Cálculo de la divergente de la función C A

Desarrollando la sumatorio para i= 1, 2 y 3


3
1 ∂ CA 1 ∂ CA 1 ∂ CA 1 ∂CA
∇ C A =∑ g⟨ i⟩ =¿ g⟨ 1⟩ + g⟨ 2 ⟩ + g ¿
i=1 √ gii ∂ x i
√ g11 ∂ x 1
√ g 22 ∂ x 2
√ g 33 ∂ x3 ⟨3 ⟩
Sustituyendo las componentes del tensor métrico:
3
1 ∂CA 1 ∂ CA 1 ∂ CA 1 ∂CA
∇ C A =∑ g i =¿ g 1 + g 2 + g ¿
i=1 √1 ∂ xi ⟨ ⟩
√1 ∂ x 1 ⟨ ⟩ 2
( x1 ) ∂ x
2 ⟨ ⟩
√ 1 ∂ x 3 ⟨3 ⟩

∂C A 1 ∂ CA ∂C A
¿ 1
g ⟨1 ⟩ + 1 2
g⟨ 2 ⟩ + 3 g ⟨ 3 ⟩
∂x x ∂x ∂x
Recordando que

x 1=r ; x 2=θ ; x3 =z
Finamente la ecuación 5 es igual al divergente de la concentración de A

∂C A 1 ∂ C A ∂ C A
∇ C A= + +
∂r r ∂θ ∂z
(5)

2. Cálculo del Laplaciano de la concentración A:

Primero, se desarrolla la ecuación (3) para coordenadas cilíndricas


3
1 ∂ ∂ CA
∇ 2 C A =∑ ( g gij
i √
)
i=1 √g ∂ x ∂x j

1 ∂ ∂C A 1 ∂ ∂CA 1 ∂ 33 ∂C A
2
∇ C A=
√ g ∂ x1 ( √ g g 11
∂x j ) +
√ g ∂ x2 ( √ g g22
∂x 2 ) +
√g ∂ x
( gg
3 √
∂ x3
)

Se sustituyen los componentes del tensor métrico y g

1 ∂ ∂C 1 ∂ ∂C A ∂ CA
2
∇ C A=
1 2
1 ( √ ( x
1 2
) (1)
A
j ) +
1 2
2
( √( x ) ( 1 )
1 2
1 2 2
)√
+
1
2

3√ 1 2
( ( x ) (1) 3 )
√ (x ) ∂x ∂x √ (x ) ∂x x ∂x ( x 1) ∂ x ∂x

2
1 ∂ 1 ∂ CA 1 ∂ x1 ∂ C A ( x1 ) ∂ ∂ C A √
¿
x1 ∂ x1
x
( +
)
∂ x j x 1 ∂ x 2 ( x 1) 2 ∂ x 2
+
(2 3
( x1 ) ∂ x ∂ x
3
) √ ( )
1 2
+√
1 ∂ 1 ∂ CA 1 x1 ∂ ∂ C A (x ) ∂ ∂ CA
¿ 1
x ∂ x1
x
(+
)
∂ x j x 1 ( x 1 )2 ∂ x 2 ∂ x 2 ( ) √(x ) 1 2 ∂ x 3 ∂ x3( )
1 ∂ 1 ∂ CA 1 ∂ CA ∂ CA
¿
x1 ∂ x1
x ( + )
+
∂ x j ( x1 )2 ∂ 2 x 2 ∂ 2 x 3
Recordando que

x 1=r ; x 2=θ ; x3 =z

1 ∂ ∂C A 1 ∂ C A ∂ C A
∇2 C A=
r ∂r(r
∂r )
+ 2 2 + 2
r ∂θ ∂ z
(6)

Sustituyendo ecuación 5 y 6 en 1:

∂C A ∂ CA 1 ∂ CA ∂C A
+ vr + vθ + vz =D AB ¿) + R A
∂t ∂r r ∂θ ∂z

También podría gustarte