R - A - L - G - 08 - Derivadas Parciales para El Jacobiano
R - A - L - G - 08 - Derivadas Parciales para El Jacobiano
R - A - L - G - 08 - Derivadas Parciales para El Jacobiano
∂ g i ∂ ln f v ∂ ln f l ∂ ln f v ∂ ln f l
( )
J ij =
∂vj
=
∂v j
−
∂v j
=
∂ vj
+
∂ lj
v v vi v
ln f =ln y i + ln φi +ln P=ln + ln φi + ln P
V
Entonces, sea el volumen total una funció n de un pará metro β que es la suma de los
volú menes molares de las especies i de la siguiente forma:
N +C
V ( ¿ β) = ∑ vi
i=1
∂ ln f iv 1 v
1 ∂ ln φi 1 1
∂ vj
= δ ij − +
vi V ∂vj
=
(
δ −1+ Φ́ ij
β y i ij )
En donde Φ está dado para el nú mero de moles de cada especie de la siguiente
manera:
∂ ln φ vi N +C
Φ́ ij=nT ( ∂ nj )
T,P
, nT = ∑ nk
k=1
1 1 1 1
J ij = (
β yi
δ ij −1+ Φ́vij + )
1−β x i (
δ ij −1+ Φ́lij )
Por lo tanto, Expandiendo el modelo del jacobiano:
∂ g1 ∂ g1 ∂ g1 ∂ g1
[ ]
…
∂ ln k 1 ∂ ln k 2 ∂ ln k Nc ∂ ln P o lnT
∂ g2 ∂ g2 ∂ g2 ∂ g1
J= …
∂ ln k 1 ∂ ln k 1 ∂ ln k Nc ∂ ln P o lnT
∂ g N +1 ∂ gN +1 ∂ g N +1 ∂ g1
c c c
…
∂ ln k 1 ∂ ln k 1 ∂ ln k Nc ∂ ln P o lnT
Por otro lado, sea una ecuació n de estado cú bica de dos pará metros a m y bm, (por ejemplo, la
ecuació n S-K-R y la ecuació n de Peng-Robinson). En donde segú n la mecá nica estadística, a m
representa el parámetro de atracció n entre las moléculas y bm, representa el covolumen
ocupado por las moléculas se tiene un modelo como el que se muestra a continuació n:
RT am
P= −
v−bm ( v+ δ 1 bm ) ( v−δ 2 bm )
V
n RT
RTIn ϕi=
−∂
∂ ni (∫[
∞
P−
V ] )
dV −RTIn( Z)
Sea el volumen como V/v, entonces se puede sustituir el término interno de la presió n por la
ecuació n de estado generalizada:
V
nRT am RnT
RTIn ϕi=
−∂
∂ ni ([
∫
∞
− −
( V −b m ) ( V +δ 1 B ) ( V +δ 2 B ) V ] )
dV −RTIn ( Z )
V
nRT am nRT
RTIn ϕi=
−∂
∂ ni (∫[
∞ ( V −b m )
dV −
( V +δ 1 b m )( V + δ 2 bm )
dV −
V ] )
dV −RTIn ( Z )
Realizando las integrales impropias por secciones se tiene que en la primera el límite
V
∫ ( VnRT
[ −bm ) ] ( [ ( ) ])
dV =¿ nRT lim ln
v→ ∞
V −b
V
=0 ¿
∞
a m=κ ( V +δ 2 bm ) + λ ( V +δ 1 b m )
a m=κV +κ δ 2 bm + λV + λ δ 1 bm
a m=V ( κ+ λ ) +b m ( λ δ 1 +κ δ 2 )
V V
am a 1 1
∫ V +δ b V +δ b dV =∫ b δ m−δ
∞ ( 1 m) ( 2 m) ∞ m( 2 1)
( )( −
)
V + δ 1 B V +δ 2 B
dV
V am V + δ 1 bm
RTIn ϕi=
−∂
∂ ni ( (
nRT ln
V −b
+ ( )) (
b m ( δ 2−δ 1 )
ln
V + δ 2 bm )(
−RTIn ( Z ) ))
Derivando término a término:
V
( )
( )
( V −bm )∗∂
V V −bm
( )
∂
∂ ni (
nRT ln
V
V −bm(=−nRT
∂
∂ ni
( ln ( u ) ) ∂ u =nRT
)) ∂ ni
( )
1
V(
V −bm
) ∂
V −bm
∂ ni
=nRT
∂ ni
V
= nRT (
Derivació n del segundo término, sacando constantes de la derivada
∂ am V + δ 1 bm 1 ∂ am V + δ 1 bm 1
(
∂ ni b m ( δ 2−δ 1 )
ln
)( = ) ( )(
V + δ 2 bm ( δ 2 −δ 1 ) ∂ ni bm
ln ) =
V + δ 2 bm ( δ 2−δ 1 )
¿
Entonces, para obtener los términos derivativos correspondientes, se hará uso de las
definiciones de los pará metros am y bm, los cuales será n derivados como funció n de
las reglas de mezclado clá sicas y de Wong y Sandler.
∂ ( n bm )
=bi
∂ ni
2
1 ∂ ( n am )
=2 ∑ x j aij
n ∂ ni j
Mientras que para las reglas de mezclado de Wong y Sandler se tiene que:
∂ ( n bm ) 1 1 ∂2 nQ Q ∂ nD
∂ ni
=
1−D n ∂ ni (
−
( 1−D ))2
1−
∂ ni ( )
2
1 ∂ n am ∂ n bm ∂ nD
n ∂n i
=RT D
∂ ni (+bm
∂ ni )
En donde después de la regla de cadena, se observa una derivada implícita de segundo
orden que se calcula como sigue a continuació n:
1 ∂2 nQ a a
n ∂ ni
=2 ∑ x j b−
j
a
(
RT ij
=2 ∑ x j
j
1
2 ) RT RT [(
bi +b j− i − j ( 1−k ij ) )]
∂ nD a ln γ ∞ i
= i +
∂ ni bi RT C
∑ x j τ ji g ji ∑ x l τ lj glj
( )
E
1 ∂ n A∞ j x j gij l
ln γ ∞ i= = +∑ τ ij −
RT ∂ ni ∑ x k gki j ∑ x k gkj ∑ x k gkj
k k k
En donde los términos derivativos de los pará metros a m y bm son obtenidos de las
reglas de mezclado que se deseen usar para cá lculo.
P ( v−b m ) ∂b
∂
∂ ni( [
−ln
RT
= ])
−RT
RTP ( v −bm )
−P m
∂ ni ( )
∂ 1 Pv ∂ (n bm) Pv ∂ ( n bm ) ∂ 1 1 ∂ (n bm)
((
∂ ni b RT
−1 )
∂ ni
=
RT
−1 )(
∂ ni ∂ ni b
+ )
b ∂ ni ( ()() )
V −b m δ 1 1 ∂ ( n am ) λ ∂ ( n bm )
∂
( 1
∂ ni RT b m ( δ 1−δ 2 )
ln (
V −b m δ 2 )[ n ∂ ni
−
b ∂ ni
=¿ ])
∂ ∂
1 ∂ 1
( ) (
RT ( δ 1−δ 2 ) ∂ ni bm
ln
V −bm δ 1
V −bm δ 2 )[ n ∂ ni
2
1 ∂ ( n am ) λ ∂ ( n bm ) V −b m δ 2
−
b ∂ ni
+
V −b m δ 1 ] ( −δ 1
∂ ni
( bm ) ( V −bm δ 2 ) +δ 2 ( b m )( V −b
( V −b m δ 2 )
∂ ni
2
∂
d
d ni
ln
(V
V −bm )
=
V −bm
V [ (V )
V −b m2
(b )
∂ ni m
]
Reduciendo y derivando el ú ltimo término con la regla de la derivada del liogaritmo
neperiano se obtiene la ecuació n de fugacidad para el componente i.
∂ bm ∂ am ∂a
ln ϕ i=
+
∂n i
RnT
V −b m
+
1
RT ( δ 2−δ 1 ) (
∂ ni
bm − m a m
bm 2
∂ ni
ln
V +δ 1 b m
(
V +δ 2 b m ))
Los términos derivativos corresponden a las reglas de mezclado cuando la ecuació n sea
derivada con respecto al nú mero de componentes:
[( ∂ am ∂a ∂
i
bm − m a m
∂ ln ϕ V −b m ∂ V 1 ∂n i ∂ ni V + δ 2 bm ∂ V +δ 1 b m V + δ 1 bm ∂
∂ ni
=
V ∂ ni( (
ln
V −b m
+
))
RT ( δ 2−δ 1 ) bm 2 ( )
V + δ 1 bm ∂ ni
ln
( ) (
V +δ 2 b m
+ ln
)
V + δ 2 bm