Libro P1 Traducido
Libro P1 Traducido
Libro P1 Traducido
es de fase en
metales y
aleaciones
TERCERA EDICIÓN
Fase
Transformacion
es en metales y
aleaciones
TERCERA EDICIÓN
DAVID A. PORTER,
KENNETH E. EASTERLING, y
MOHAMED Y. SHERIF
Prensa CRC
Grupo Taylor & Francis
6000 Broken Sound Parkway NW, Suite 300 Boca
Ratón, FL 33487-2742
Este libro contiene información obtenida defuentes autónticas y de gran prestigio. Se han hecho esfuerzos razonables para
publicar datos e información confiables, pero el autor y el editor no pueden asumir la responsabilidad de la validez de todos
los materiales o las consecuencias de su uso. Los autores y publishers han intentado rastrear a los titulares de los derechos
de autor de todo el material reproducido en esta publicación y se disculpan con los titulares de los derechos de autor si no
se ha obtenido permiso para publicar en esta forma. Si algún material con derechos de autor no ha sido reconocido, por
favor escríbanos y háganoslo saber para que podamos rectificar en cualquier reimpresión futura.
Excepto según lo permitido por la Ley de Derechos de Autor de los Estados Unidos, ninguna parte de este libro puede ser
reimpresa, reproducida, transmitida o utilizada de ninguna forma por ningún medio electrónico, mecánico o de otro tipo,
ahora conocido oinventado en el futuro, incluyendo fotocopias, microfilmación y grabación, o en cualquier sistema de
almacenamiento o recuperación de información, sin el permiso por escrito de los editores.
Para obtener permiso para fotocopiar o usar material electrónicamente de este trabajo, acceda a www.copyright.com
(http://www.copyright.com/) o comuníquese con el Copyright Clearance Center, Inc. (CCC), 222 Rosewood Drive,
Danvers, MA 01923, 978-750-8400. CCC es una organización sin fines de lucro que proporciona licencias y registro para
una variedad de sussirvientes. Para las organizaciones a las que el CCC les ha otorgado una licencia de fotocopia, se ha
dispuesto un sistema de pago separado.
Aviso de marca comercial: Los nombres de productos o corporativos pueden ser marcas comerciales o marcas comerciales
registradas, y se utilizan solo para identificación y explicación sin intención de infringir. Visite el sitio web de Taylor &
Francis en http://www.taylorandfrancis.com
en
vi Contenido
Referencias....................................................................................................... 60
Lecturas adicionales ........................
..................................................................... 60
Ejercicios.......................................................................................................... 61
xi
Prefacio a la Segunda Edición
En los diez años transcurridos desde que este libro se publicó por primera vez, ha
habido muchos desarrollos nuevos en el campo metalúrgico. Los metales y vidrios
rápidamente solidificados han alcanzado la mayoría de edad; las nuevas aleaciones
de Al-Li se utilizan ahora en aviones modernos; Los aceros microaleados
(estructurales) y dealta pureza (tuberías) se han vuelto más sofisticados; han
aparecido aceros dispersos por óxido radicalmente nuevos; se han desarrollado una
serie de nuevos metales de memoria; la lista podría continuar. A pesar de esto, los
principios subyacentes que rigen todos estosdesarrollos obviamente no han
cambiado a lo largo de los años. Esta es realmente la fuerza del presente libro de
texto. Desde el principio, nuestro objetivo fue explicar estos principios de una
manera agradable y legible, y una en la que los estudiantes universitarios pudieran
apreciar y ser capaces de desarrollarse por sí mismos. Por lo tanto, el presente texto
ha cambiado deliberadamente poco con respecto al original. Sin embargo, esperamos
haber corregido cualquier error, ampliado las listas de lecturas adicionales y quizás,
lo más importante, incluido un conjunto completo de soluciones a los ejercicios.
Esperamos que la edición revisada continúe siendo disfrutada y apreciada en las
muchas Escuelas de Metalurgia, Ciencia de Materiales y Materiales de Ingeniería
que sabemos que están utilizando nuestro texto en todo el mundo.
Al completar esta edición revisada, estamos agradecidos a las muchas personas,
estudiantes y profesores por igual, que nos han escrito durante la última década. Se
debe agradecer especialmente al Dr. Wen-Bin Li (Universidad de Luleå) por utilizar
un peine de dientes finos para resaltar errores obvios y menos obvios enel texto
original. Quedan, (inevitablemente), algunos "puntos de discordia" con respecto a
nuestra descripción de ciertos fenómenos, como lo plantearon algunos de nuestros
corresponsales, pero no hay nada insalubre en eso. Por último, nos gustaría dar las
gracias al Dr. John Ion (Universidad de Lappeenranta, Finlandia) por su ayuda en la
compilación del capítulo Soluciones a los Ejercicios.
xiii
Prefacio a la Primera Edición
Este libro está escrito como un curso de pregrado en transformaciones de fase para
estudiantes de último año especializados en metalurgia, ciencia de materiales o
materiales de ingeniería. También debería ser útil para los estudiantes de
investigación interesados en revisar susconocimientos sobre el tema. El libro se basa
en conferencias impartidas originalmente por los autores en la Universidad de Luleå
para estudiantes de ingeniería especializados en materiales de ingeniería.
Sorprendentemente, no encontramos tratamientos modernos de este importante tema
en una forma adecuada para un libro de texto, el más reciente probablemente sea
Transformations in Metals de P.G. Shewmon (McGraw-Hill, 1969). Sin embargo,
ha habido algunos desarrollos notables en el tema durante la última década,
particularmente en estudios de interfaces entre fases ymigración de interfaz d, así
como la cinética del crecimiento precipitado y la estabilidad de precipitados.
También ha habido una serie de nuevos desarrollos prácticos importantes basados en
transformaciones de fase, incluida la introducción de aceros TRIP (transformation
inducida por deformación plástica), compuestos eutécticos alineados
direccionalmente y nuevos aceros estructurales sofisticados con propiedades
superiores de soldabilidad y conformado, por mencionar solo algunos. Además, la
fundición continua y la soldadura por fusión de alta velocidad y alta energía han
surgido con fuerza en los últimos años como importantes aplicaciones de producción
de solidificación. El objetivo de este curso era presentar un tratamiento de las
transformaciones de fase en el que estos y otros nuevos desarrollos pudieran
explicarse en términos de los principios básicos de la termodinámica y los
mecanismos atómicos.
El libro está efectivamente en dos partes. Los capítulos 1-3 contienen el material
de fondo necesario para comprender las transformaciones de fase: termodinámica,
cinética, teoría de la difusión y la estructura y propiedades de las interfaces. Los
capítulos 4-6 tratan de las transformaciones específicas: solidificación,
transformaciones de difusión en sólidos y transformaciones sin difusión. Al final de
los capítulos sobre solidificación, transformaciones controladas por difusióny
martensita, ofrecemos algunos estudios de caso seleccionados de aleaciones de
ingeniería para ilustrar algunos de los principios discutidos anteriormente. De esta
manera, esperamos que el texto proporcione un vínculo útil entre la teoría y la
realidad práctica. Cabe señalar que nos pareció necesario impartir este curso junto
con una serie de ejercicios prácticos de laboratorio y ejemplos trabajados. También
se incluyen conjuntos de problemas al final de cada capítulo del libro.
xvi Prefacio a la Primera Edición
1.1 Equilibrio
Es útil comenzar este capítulo sobre termodinámica definiendo algunos de los
términos que se usarán con frecuencia. En el study de transformaciones de fase
trataremos los cambios que pueden ocurrir dentro de un sistema dado, por ejemplo,
una aleación que puede existir como una mezcla de una o más fases. Una fase puede
definirse como una porción del sistema cuyas propiedades y composición
sonhomómicas y que es físicamente distinta de otras partes del sistema. Los
componentes de un sistema dado son los diferentes elementos o compuestos
químicos que componen el sistema, y la composición de una fase o el sistema se
puede describir dandolas cantidades relativas de cada componente.
El estudio de las transformaciones de fase, como su nombre indica, se refiere a
cómo una o más fases en una aleación (el sistema) cambian a una nueva fase o mezcla
de fases. La razón por la que se produce una transformación es porque el estado
inicial de la aleación es inestable en relación con el estado final. Pero, ¿cómo se mide
la estabilidad de fase? La respuesta a esta pregunta es proporcionada por la
termodinámica. Para las transformaciones que ocurren a temperatura y presión
constantes, la estabilidad relative de un sistema está determinada por su energía libre
de Gibbs (G).
1
2 Transformaciones de fase en metales y aleaciones
G = H −TS (1,1)
dG = 0 (1,3)
Se puede ver en la definición de G, Ecuación 1.1, que el estado con mayor estabilidad
seráel mejor compromiso entre baja entalpía y alta entropía. Por lo tanto, a bajas
temperaturas, las fases sólidas son más estables, ya que tienen la unión atómica más
fuerte y, por lo tanto, la energía interna más baja (entalpía). A altas temperaturas, sin
embargo, dominan los 2TS term y las fases con más libertad de movimiento de
átomos, líquidos y gases, se vuelven más estables. Si se consideran los cambios de
presión, se puede ver en la Ecuación 1.2 que las fases con pequeños volúmenes se
ven favorecidas por las altas presiones.
La definición deequilibrium dada por la ecuación 1.3 se puede ilustrar
gráficamente de la siguiente manera. Si fuera posible evaluar la energía libre de un
sistema dado para todas las configuraciones concebibles, se encontraría que la
configuración de equilibrio estable tiene la energía libre más baja. This se ilustra en
Termodinámica y diagramas de fases 3
la Fig. 1.1 donde se imagina que las diversas configuraciones atómicas pueden ser
representadas por puntos a lo largo de la abscisa. La configuración A sería
d =0
d =0
B U
Disposición de los átomosn
Figura 1.1
Una variación esquemática de la energía libre de Gibbs con la disposición de los átomos. La configuración
'A' tiene la energía libre más baja y, por lo tanto, es la disposición cuando el sistema está en equilibrio
estable. La configuración 'B' es unequilibrio metastásico.
∆G = G2 −G1 < 0 (1.4) donde G1 y G2 son las energías libres de los estados inicial y
final respectivamente. La transformación no necesita ir directamente al estado de
equilibrio estable, sino que puede pasar a través de toda una serie de estados
metaestables intermedios.
La respuesta a la pregunta "¿Qué tan rápido ocurre una transformación de fase?"
no es proporcionada por la termodinámica clásica. A veces los estados metaestables
pueden ser muy duros; en otras ocasiones pueden existir casi indefinidamente como
en el caso del diamante a temperatura y presión ambiente. La razón de estas
diferencias es la presencia de la joroba de energía libre entre los estados metaestable
y estable en la Fig. 1.1. El análisisde las tasas de transformación en química física
pertenece al ámbito de la cinética. En general, las jorobas más altas o las barreras de
energía conducen a tasas de transformación más lentas. La cinética obviamente juega
un papel central en el estudio de las transformaciones de fase y muchos exámenesde
los procesos cinéticos se encontrarán a lo largo de este libro.
Las diferentes funciones termodinámicas que se han mencionado en esta sección
se pueden dividir en dos tipos llamados propiedades intensivas y extensivas. Las
propiedades intensivas son aquellas que son independientesdel tamaño del sistema,
como T y P, mientras que las propiedades extensivas son directamente
proporcionales a la cantidad de material en el sistema, por ejemplo, V, E, H, S y G.
La forma habitual de medir el tamaño del sistema es por el número de molesde
material que contiene. Las propiedades extensivas son entonces cantidades molares,
es decir, expresadas en unidades por mol. El número de moles de un componente
dado en el sistema está dado por la masa del componente en gramos dividida por su
peso atómico o molecular.
El número de átomos o moléculas dentro de 1 mol de material está dado por el
número de Avogadro (Na) y es de 6.023 × 1023.
H
∂
Cp = ∂T P (1.5)
Pendiente
=
0 (K)
298
U)
n
0 (K)
0
Figura 1.2
(a) Variación de Cp con temperatura, Cp tiende a un límite de ~3R. (b) Variación de entalpía (H) con
temperatura absoluta para un metal puro, (c) Variación de entropía (S) con temperatura absoluta.
son sólo los cambios en las funciones termodinámicas los que son de interés. En
consecuencia, H se puede medir en relación con cualquier nivel de referencia, lo que
generalmente se hace definiendo H = 0 para un elemento puro en su estado más
estable a 298 K (25 ° C). La variación de H con T se puede calcular integrando la
ecuación 1.5, es decir,
H = ∫ T CpdT (1.6)
6 Transformaciones de fase en metales y aleaciones
298
Cp = ∂∂TS P (1.7)
Tomando la entropía a 0 K como cero, la ecuación 1.7 se puede integrar para dar
S = ∫0T Tp dT
C
(1.8)
=−
Figura 1.3
Variación de la energía libre de Gibbs con la temperatura.
∂G = −S
∂T (1.10)
P
Esto significa que G disminuye con el aumento de T a una tasa dada por 2S. Las
posiciones relativas de las curvas de energía libre de las fases sólida y líquida se
ilustran en la Fig. 1.4. A todas las temperaturas el líquido tiene una entalpía (energía
interna) más alta que el sólido. Por lo tanto a bajas temperaturas GL > GS. Sin
embargo, la fase líquida tiene una entropía más alta que la fase sólida y, por lo tanto,
la energía libre de Gibbs del líquido disminuye más rápidamente con el aumento de
la temperatura que la del sólido. Para temperaturas de hasta Tm, la fase sólida tiene
la energía libre más baja y, por lo tanto, es la fase de equilibrio estable, mientras que
por encima de Tm la fase líquida es el estado de equilibrio del sistema. En Tm ambas
fases tienen el mismo valor de G y tanto sólido como líquido pueden existir en
equilibrio. Tm es, por lo tanto, la temperatura de fusión de equilibrio a la presión en
cuestión.
d
(K)
298
m
U
n
Figura 1.4
Variación de entalpía (H) y energía libre (G) con temperatura para las fases sólida y líquida de un metal
puro. L es el calor latente de la fusión, Tm la temperatura de fusión de equilibrio.
8 Transformaciones de fase en metales y aleaciones
1600
d-Hierr
o
1200
c-Hierr
o
800
400
a- Hierr Hierr
o o
0
25 50 75 100 125 150 175
Presión (kbar)
Figura 1.5
Efecto de la presión sobre el diagrama de fase de equilibrio para el hierro puro.
Termodinámica y diagramas de fases 9
Por ejemplo, la Fig. 1.5 muestra el efecto de la presión sobre las temperaturas de
equilibrio para el hierro puro. El aumento de la presión tiene el efecto de deprimir la
temperatura de α / γ equilibrium y elevar la temperatura de fusión de equilibrio. A
presiones muy altas hcp ε-Fe se vuelve estable. La razón de estos cambios se deriva
de la Ecuación 1.9. A temperatura constante, la energía libre de una fase aumenta
con la presión tal que
∂G
T =V (1.11)
∂P
Si las dos fases en equilibrio tienen diferentes volúmenes molares, sus respectivas
energías libres no aumentarán en la misma cantidad a una temperatura dada y, por lo
tanto, el equilibrio se verá perturbado por los cambios en la presión. La única forma
de mantener el equilibrio a diferentes presiones es variando la temperatura.
Si las dos fases en equilibrio son α y β, aplicación de la ecuación 1.9 a
1 mol de ambos da
dP S Sα S
= ββ −−V α = ∆∆V (1.13) dT V
eq m m
Gα = Hα −TSα Gβ = Hβ −TSβ
GL = H L −TSL
GS = H S −TSS
∆ L
Temperatura
Termodinámica y diagramas de fases 11
Figura 1.6
Diferencia de energía libre entre líquido y sólido cerca del punto de fusión. La curvatura de las líneas GS
y GL ha sido ignorada.
∆G = ∆H −T S∆ (1.15) donde
∆H = H L − H S y ∆S = SL − S
∆G = ∆H −Tm∆S = 0
y por lo tanto en Tm
∆H L
∆S = = (1.16)
Tm T mm
L
∆G L−T
Tm
∆G L TT∆m
(1.17)
12 Transformaciones de fase en metales y aleaciones
Este es un resultado muy útil que se repetirá con frecuencia en capítulos posteriores.
X A+ X B= 1 (1,18)
y X A y X B son las fracciones molares de A y B respectivamente en la aleación. Para
calcular la energía libre de la aleación, la mezcla se puede realizar en dos pasos (ver
Fig. 1.7). Estos son:
Figura 1.7
Energía libre de mezcla.
Después del paso 1 la energía libre del sistema viene dada por
G 1 = X GA + X GB B Jmol−1 (1.19)
U
n
0 B 1
A B
14 Transformaciones de fase en metales y aleaciones
Figura 1.8
Variación de G1 (la energía libre antes de mezclar) con la composición de la aleación (X A o X B).
La energía libre del sistema no permanecerá constante durante la mezcla de los
átomos A y B y después del paso 2 la energía libre de la solución sólida G2 se puede
expresar como
G 1 = H1 −TS1
G2 = H 2 −TS2
Poner
∆H = H 2 − H1
∆S mix = S2 − S1 da
ω
S = kln (1.23) donde k es la constante de Boltzmann y ω es una medida de
aleatoriedad. Hay dos contribuciones a la entropía de una solución sólida: una
contribución térmica Sth y una contribución configuracional Sconfig.
En el caso de la entropía térmica, ω es el número de formas en que la energía
térmica de latapa so puede dividirse entre los átomos, es decir, el número total de
formas en que se pueden establecer vibraciones en el sólido. En las soluciones, existe
aleatoriedad adicional debido a las diferentes formas en que se pueden organizar los
átomos. Esto da una configuración S de entropía adicional para la cual ω es el
número de formas distinguibles de organizar los átomos en la solución.
Si no hay cambio de volumen o cambio de calor durante la mezcla, entonces la
única contribución a la mezcla ΔS es el cambio en la entropía confi guracional. Antes
de mezclar, la A y la B enoms se mantienen por separado en el sistema y solo hay
una forma distinguible en la que se pueden organizar los átomos. En consecuencia
S1 = k ln l = 0 y por lo tanto ΔSmezcla = S2.
Suponiendo que A y B se mezclan para formar una solución sólida sustitutiva y
que todas las confi guraciones de átomos A y B son igualmente probables, el número
de formas distinguibles de organizar los átomos en los sitios atómicos es
N A = X NA y
N B = X NB a
∆S = −R X( A ln X A + X B ln X B) (1.25)
16 Transformaciones de fase en metales y aleaciones
∆G mix = RT X( A ln X A + X B ln X B) (1.26)
A B
Figura 1.9
Energía libre de mezcla para una solución ideal.
Esto se muestra esquemáticamente en la Fig. 1.10. Nótese que, a medida que aumenta
la temperatura, GA y GB disminuyen y las curvas de energía libre asumen un
Termodinámica y diagramas de fases 17
0 B 1
Figura 1.10
La energía libre molar (energía libre por mol de solución) para una solución sólida ideal. Una combinación
de las figs. 1.8 y 1.9.
curvatura. La disminución de GA y GB se debe a la entropía térmica de ambos
componentes y viene dada por la Ecuación 1.10.
Cabe señalar que todos los diagramas de composición de energía libre en este libro
son esencialmente esquemáticos; si se trazan correctamente, las curvas de energía
libre deben terminar asintóticamente en los ejes verticales de los componentes puros,
es decir, tangenciales a los ejes verticales de los diagramas. Esto se puede demostrar
diferenciando la ecuación 1.26 o 1.27.
= G
μA ∂∂n A′ T P n, ,B
(1.29)
µB
µU
n
A XB B
Figura 1.11
La relación entre la curva de energía libre para una solución y los potenciales químicos de los
componentes.
Termodinámica y diagramas de fases 19
Cuando G se conoce como una función de X A y X B, como en la Fig. 1.10 por ejemplo,
(μA y μB se puede obtener extrapolando la tangente a la curva G a los lados del
diagrama de energía libre molar como se muestra en la Fig. 1.11. Esto se puede
observara partir de las ecuaciones 1.30 y 1.31, recordando que X A + X B = 1, es decir,
dX A = –dX B, y esto se deja como un ejercicio para el lector. De la Fig. 1.11 se
desprende claramente que μA y μB varían sistemáticamente con la composición de la
fase.
La comparación de las ecuaciones 1.27 y 1.31 da μA y μB para una solución ideal
como
b B
U U
n
n
– – B
U
n
µU c µB
d
n
U B
XB
n
Figura 1.12
La relación entre la curva de energía libre y los potenciales químicos para una solución ideal.
20 Transformaciones de fase en metales y aleaciones
Figura 1.13
Los diferentes tipos de enlace interatómico en una solución sólida.
Al considerar que la energía cero es el estado en el que los átomos se separan hasta
el infinito εAA, εBB y εAB son cantidades negativas, y se vuelven increíblementemás
negativos a medida que los enlaces se vuelven más fuertes. La energía interna de la
solución E dependerá del número de enlaces de cada tipo PAA, PBB y PAB tal que
ε
∆H mix = PAB (1.33) donde
Si ε < 0 los átomos en la solución preferirán estar rodeados por átomos del tipo
opuesto y esto aumentará PAB, mientras que, si ε > 0, PAB tenderá a ser menor que en
una solución aleatoria. Sin embargo, siempre ε no sea demasiado diferente de cero,
la ecuación 1.35 isigue siendo una buena aproximación, en cuyo caso
∆H = ΩX XA B (1.36)
Dónde
Ω = N zy ε (1.37)
∆G mix = Ω XX A B + RT X ( A ln X A + X B ln X B )
(1.38)
mezcla ∆H −T S∆ mezcla
22 Transformaciones de fase en metales y aleaciones
O
h
B B
Figura 1.14
La variación de la mezcla de ΔH con la composición para una solución regular.
G = X GA + X GB B + ΩX XA B + RT X( A ln X A + X B ln X B) (1.39)
μA = GA + Ω(1− X A )2 + RT ln X A (1.40)
y
Termodinámica y diagramas de fases 23
μB = GB + Ω(1− X B)2 + RT ln X B
+ B + B
0 0
∆
∆
– ∆
– –
– ∆
∆
Mezcla A ∆G B A B
B
0 0
∆ B
– –
– ∆
– ∆
U B U B
n n
(d) Oh>
( c) Oh > O,
alto T
Figura 1.15
El efecto de la mezcla ΔH y T sobre la mezcla ΔG.
1.3.5 Actividad
La expresión 1.32 para el potencial químico de una aleación ideal era simple y es
conveniente retener una expresión similar para cualquier solución. Esto se puede
hacer definiendo la actividad de un componente, a, tal que las distancias ac y bd en
la Fig. 1.16 son 2RT ln a A y 2RT ln aB. En este caso
μA = GA + RT ln a A (1.41) y
μB = GB + RT ln aB
A = Ω (1− X A )2
ln XA RT
Figura 1.16
La relación entre la energía libre molar y la actividad.
aB = Ω (1− X B)2 ln
XB RT
γA = a A / X A (1.43)
1 1
3
B 3 U
n
1 1
2
2
0 0
B 1 1 U 0
B U n B
n
Termodinámica y diagramas de fases 25
(a) (b)
Figura 1.17
La variación de la actividad con la composición (a) a B (b) a A. Línea 1: solución ideal (ley de Raoult). Línea
2: Mezcla ΔH < 0. Línea 3: Mezcla ΔH > 0.
Para una solución diluida de B en A, la ecuación 1.42 se puede simplificar dejando
que X B → 0 en cuyo caso
aB (1.44)
γB = constante (ley de Henry)
XB
La ecuación 1.44 se conoce como la ley de Henry y 1.45 como la ley de Raoult; se
aplican a todas las soluciones cuando se diluyen suficientemente.
Dado que la actividad está simplemente relacionada con el potencial químico a
través de la ecuación 1.41, la actividad de un componente esotro medio para describir
el estado del componente en una solución. No se proporciona información adicional
y su uso es simplemente una cuestión de conveniencia, ya que a menudo conduce a
matemáticas más simples.
La actividad y el potencial químico son simplemente una medida de ladensidad de
un átomo para dejar una solución. Si la actividad o el potencial químico es bajo, los
átomos son reacios a abandonar la solución, lo que significa, por ejemplo, que la
presión de vapor del componente en equilibrio con la solución será relativamente
baja. También será evidente más adelante que la actividad o el potencial químico de
un componente es importante cuando varias fases condensadas están en equilibrio.
para lograr la energía libre mínima. En sistemas con ε < 0 la energía interna del
sistema se reduce aumentando el número de enlaces A-B, es decir, ordenando los
átomos como se muestra en Fig. 1.18a. Si ε > 0, la energía interna puede reducirse
aumentando el número de enlaces A-A y B-B, es decir, mediante la agrupación de
los átomos en grupos ricos en A y B- ricos, Fig. 1.18b. Sin embargo, el grado de
ordenación o agrupación disminuirá a medida que aumente la temperatura debido a
la creciente importancia de la entropía.
Figura 1.18
Representación esquemática de soluciones sólidas: (a) sustitución ordenada, (b) agrupamiento, (c)
intersticial aleatorio.
En sistemas donde hay una diferencia de tamaño entre los átomos, elmodelo emical
cuasicónico subestimará el cambio en la energía interna en la mezcla, ya que no se
tienen en cuenta los campos de deformación elástica que introducen un término de
energía de deformación en la mezcla ΔH. Cuando la diferencia de tamaño es grande,
este efecto puede dominar sobre el ter m químico.
Cuando la diferencia de tamaño entre los átomos es muy grande, entonces las
soluciones sólidas intersticiales son energéticamente más favorables, Fig. 1.18c. Se
necesitan nuevos modelos matemáticos para describir estas soluciones.
En los sistemas donde hay un fuerte enlace químico entrelos átomos hay una
tendencia a la formación de fases intermetálicas. Estos son distintos de las soluciones
basadas en los componentes puros, ya que tienen una estructura cristalina diferente
y también pueden ser muy ordenados. Las fases intermedias y las fases ordenadasse
discuten más a fondo en las siguientes dos secciones.
enlaces A-B es mayor que este, se dice que la solución contiene orden de corto
alcance (SRO). El grado de ordenación se puede cuantificar definiendo un parámetro
SRO tal que
s= PAB − PAB(aleatorio)
PAB(max)− PAB(aleatorio)
(a) (b)
Figura 1.19
(a) Solución aleatoria A-B con un total de 100 átomos y X A = X B = 0.5, PAB ~ 100, S = 0. (b) Misma
aleación con orden de corto alcance PAB = 132, PAB(max) ~ 200, S = (132 –100)/ (200 – 100) = 0.32.
Co
n
Figura 1.20
Estructuras sustitutivas ordenadas en el sistema Cu-Au: (a) estructura desordenada a alta temperatura, (b)
superred CuAu, (c) superred Cu 3Au.
Co
n
0 0.1 0.20.30.40.50.60.70.80.91.0
Cu XAu Au
Figura 1.21
Parte del diagrama de fase Cu-Au que muestra las regiones donde las superredes Cu 3 Au y CuAu son
estables.
Las celosías ordenadas más comunes en otros sistemas se resumen en la Fig. 1.22
junto con su notación Structurbericht y ejemplos de aleaciones en las quese
encuentran t hey. Finalmente, tenga en cuenta que la temperatura crítica para la
pérdida de orden de largo alcance aumenta con el aumento de Ω, o mezcla ΔH, y en
muchos sistemas la fase ordenada es estable hasta el punto de fusión.
Figura 1.22
Las cinco celosías ordenadas comunes, ejemplos de las cuales son: (a) L2 0:CuZn, FeCo, NiAl, FeAl,
AgMg; b) Ll2:Cu3Au, Au 3Cu, Ni3Mn, Ni3Fe, Ni3Al, Pt3Fe; c) Ll0:CuAu, CoPt, FePt; (d) D03:Fe3Al, Fe3Si,
Fe3Be, Cu 3Al; (e) D019:Mg3Cd, Cd 3Mg, Ti3Al, Ni3Sn. (After R.E. Smallman, Modern Physical
Metallurgy, 3ª edición, Butterworths, Londres, 1970.)
30 Transformaciones de fase en metales y aleaciones
B
B
U
n
∆
U B U B
B
n Ideal n
(a) composición (b)
Figura 1.23
Curvas de energía libre para fases intermedias: (a) para un compuesto intermetálico con un rango de
estabilidad muy estrecho, (b) para una fase intermedia con un amplio rango de estabilidad.
Algunas fases intermedias pueden sufrir transformaciones de desorden de orden
en las que unadisposición aleatoria de los átomos es estable a altas temperaturas y
una estructura ordenada es estable por debajo de alguna temperatura crítica. Dicha
transformación ocurre en la fase de β en el sistema Cu-Zn, por ejemplo (ver Sección
5.10).
La estructura de lasfases intervenidas está determinada por tres factores
principales: tamaño atómico relativo, valencia y electronegatividad. Cuando los
átomos componentes difieren en tamaño por un factor de aproximadamente 1.1-1.6,
es posible que los átomos llenen el espacio de manera más eficiente si los átomos los
ordenanen una de las llamadas fases de Laves basadas en MgCu 2, MgZn2 y MgNi2,
Fig. 1.24. Otro ejemplo donde el tamaño atómico determina la estructura es en la
formación de los compuestos intersticiales MX, M2X, MX2 y M6X donde M puede
ser Zr, Ti, V, Cr, etc. y X puede ser H, B, C y N. En este caso, los átomos M forman
una disposición cerrada cúbica o hexagonal y los átomos X son lo suficientemente
pequeños como para entrar en los intersticios entre ellos.
La valencia relativa de los átomos se vuelve importante en los llamados pha ses
de electrones, por ejemplo, α y β latones. La energía libre de estas fases depende del
número de electrones de valencia por unidad de celda, y esto varía con la
composición debido a la diferencia de valencia.
La electronegatividad de un átomo es una medida de la fuerza con la que atrae a
losatrilones e y en sistemas donde los dos componentes tienen electronegatividades
Termodinámica y diagramas de fases 31
Figura 1.24
La estructura de MgCu2 (fase A Laves). (De J.H. Wernick, capítulo 5 en Metalurgia Física,
2ª edn., R.W. Cahn (Ed.) Holanda septentrional, 1974.)
c
f c
U ∆ U
b b
n n
U
y n U B
n
U B B B
(Un) (b )
n
Figura 1.25
a) La curva de energía libre molar para la fase α , b) las curvas de energía libre molar para las fases α y β
.
32 Transformaciones de fase en metales y aleaciones
y
U U
n n
d g
U b
U n U o B B
nB B B
c
n
Figura 1.26
La energía libre molar de una mezcla bifásica (α + β).
La energía libre molar de la mezcla de fase g está dada por el punto en la línea recta
entre α y β como se muestra en la figura. Este resultado se puede probar más
fácilmente utilizando la geometría de la Fig. 1.26. Las longitudes ad y cf representan
Termodinámica y diagramas de fases 33
o
o
U U U
n U no no
n
o o
(a) La aleación X 0 tiene una energía libre G1 como una mezcla de α1 + β1. (b) En equilibrio, la aleación X
0
tiene una energía libre mínima Ge cuando es una mezcla de αe + βe.
α y β fases hasta alcanzar las composiciones αe y βe, Fig. 1.27b. La energía libre del
sistema Ge es ahora un mínimo y no hay deseo de más cambios. En consecuencia, el
sistema está en equilibrio y αe y βe son las composiciones de equilibrio de las fases
α y β.
Este resultado es bastante general y se aplica a cualquier aleación con una
composición global entre αe y βe: solo cambian las cantidades relativas de las dos
fases, según lo dado por la regla de la palanca. Sin embargo, cuando la composición
de la aleación se encuentra fuera de este rango, la energía libre mínima se encuentra
en las curvas Gα o Gβ y el estado de equilibrio de la aleación es unamonofásica
homómica.
De la Fig. 1.27 se puede ver que el equilibrio entre dos fases requiere que las
tangentes a cada curva G en las composiciones de equilibrio se encuentren en una
línea común. En otras palabras, cada componente debe tener el mismo potencial
químico enlas dos fases, es decir, para un equilibrio heterogéneo:
1
.5 Diagramas de fases binarias
En la sección anterior se ha demostrado cómo se puede obtener el estado de
equilibrio de una aleación a partir de las curvas de energía libre a una temperatura
dada.
El siguiente paso es ver cómo el equilibrio se ve afectado por la temperatura.
Termodinámica y diagramas de fases 35
– GU
U
n GB
n
– B
U Uy By B
1n n s
q T
p
0
1
s T q
U
n
0 p
0 1
XB
Figura 1.28
La variación de una A y unaB con composición para un sistema binario que contiene dos soluciones ideales,
α y β.
S S
S
L L
L
U B U B U b c B
(U) n (nb) (nc)
n
m (B) 3 U
m (U)
1
n
c n
2
L b
L
m (B)
S 3 d
S
U B U B U B
(d ) n (ny) (nf)
Figura 1.29
La derivación de un diagrama de fases simple a partir de las curvas de energía libre para el líquido (L) y
el sólido (S).
Fig. 1.30c, y la solución sólida es más estable como una mezcla de dos fases α′ y α"
con composiciones e y f. A temperaturas más altas, cuando − la mezcla de TΔS se
hace más grande, e y f se acercan entre sí y eventualmente desaparecen como se
muestra en el diagrama de fases, Fig. 1.30d. La región de α′ + α" se conoce como
brecha de miscibilidad.
El efecto de una mezcla positiva de ΔH en el sólido ya es evidente a temperaturas
más altas donde da lugar a una mezcla depunto de fusión mini mum. La razón por la
que todas las aleaciones deben fundirse a temperaturas por debajo de los puntos de
fusión
B U c
n
b)
c
2
U
n
y
f 3
U
B
(d n
Figura 1.30
La derivación de un diagrama de fases donde ∆HmezclaS > ∆HmezclaL = 0. Curvas de energía libre v.
composición para (a) T1, (b) T2 y (c) T3.
U
B
n
U C
n a+b
b+c
B B U B
B
n
(a) (b)
Figura 1.31
a) Diagrama de fases cuando ∆HmezclaS < 0; (b) como (a) pero aún más negativo ∆HmixS (After R.A. Swalin,
Thermodynamics of Solids, John Wiley, Nueva York, 1972).
U B
n
B)
Mínimo G
Figura 1.35
Diagrama de energía libre para ilustrar que el rango de composiciones sobre las cuales una fase es estable
depende de las energías libres de las otras fases en equilibrio.
La prueba de esta relación se deja como un ejercicio para el lector (ver Ejercicio
1.10). Una consecuencia deesta condición general es la regla de fase de Gibbs. Esto
establece que si un sistema que contiene componentes C y fases P está en equilibrio,
el número de grados de libertad F viene dado por
P+F=C+2 (1,49)
Un grado de libertad es una variable intensiva como T, P, XA, X B ... que se pueden
variar de forma independiente sin dejar de mantener el equilibrio. Si la presión se
mantiene constante, se pierde un grado de libertad y la regla de fase se convierte en
P+F=C+1 (1.50)
P+F=3
Esto significa que un sistema binario que contiene una fase tiene dos grados de
libertad, es decir, T y X B se pueden variar independientemente. En una región
bifásica de un diagrama de fases P = 2 y por lo tanto F = 1 lo que significa que si la
temperatura se eligeindependientemente las composiciones de las fases son fijas.
Cuando tres fases están en equilibrio, como a una temperatura eutéctica o pertéctica,
44 Transformaciones de fase en metales y aleaciones
∆ B
U B
n
B
B = µB
U B
(b) y B
n B
U B
n
a+b
U B
y B
(U) n B
n
Figura 1.36
Solubilidad de B en A.
μ μαB = Bβ GBβ
Termodinámica y diagramas de fases 45
X Be = exp { −
∆GB + Ω
} (1.52)
RT
Poner
∆GB = ∆HB −T S∆ B da
− Q
X Be = A exp RT (1,53)
Q = ∆HB + Ω (1,54)
∆H ∆H Xvv
Por lo tanto, la energía libre molar del cristal que contiene Xv mol de vacantes está
dada por
G = GA + ∆G = GA + ∆H Xv v −T S X∆ vv v
+ RT(X v ln X v +(1− X v )ln (1− X v )) (1.55)
dG= 0
dX vX =Xv e v
∆H en −T S∆ v + RT ln X ene = 0
y
E
O n
Figura 1.37
Concentración de vacantes de equilibrio.
Gv
e = exp −∆ (1.57)
X en
RT
(U)
n
B
r
∆ C
B
∞
B
∞ r B
(b)
Figura 1.38
El efecto de la energía interfacial sobre la solubilidad de partículas pequeñas.
∆P = 2γ
r
Por definición, la energía libre de Gibbs contiene un término 'PV' y, por lo tanto,
un aumento de la presión P causa un aumento de la energía libre G. De la ecuación
1.9 a temperatura constante
∆G = ∆P V⋅
Por lo tanto, la curva de β en el diagrama de composición de energía libre de molares
en la Fig. 1.38b se elevará en una cantidad
Figura 1.39
Transferencia de dn mol de β de partícula grande a pequeña.
Por ejemplo, considere una derivación alternativa de la ecuación 1.58 que puede
modificarse más fácilmente para tratar casos no esféricos3.
Considere un sistema que contiene dos partículas β, una con una interfaz esférica
de radio r y la otra con una interfaz plana (r = ∞) incrustada en una matriz α como
se muestra en la Fig. 1.39. Si la diferencia de energía libre molar entre las dos
partículas es ΔGγ , la transferencia de una pequeña cantidad (dn mol) de β de la
50 Transformaciones de fase en metales y aleaciones
partícula grande a la pequeña aumentarála energía libre del sistema en una pequeña
cantidad (dG) dada por
dG = ∆Gγdn
dG = γdA
dA
∆Gγ = γ (1,59)
dn
dA = dA/ dr = 2Vm
dn dn/ dr r
X∞
Del mismo modo, X r se puede obtener utilizando (ΔGB − 2γVm / r) en lugar de ΔGB
Por lo tanto
Termodinámica y diagramas de fases 51
Xr X∞ 1+ 2RTrγVm
por ejemplo, para r = 10 nm Xr/ X∞ ~ 1.1. Por lo tanto, se puede ver que pueden surgir
diferencias de solubilidad bastante grandes para partículas en el rango r = 12100 nm.
Sin embargo, para las partículas visibles en el microscopio de luz (r > 1 mm) los
efectos de capilaridad son muy pequeños.
Figura 1.40
2
.7 Equi ternarioLibrium
Dado que la mayoría de las aleaciones comerciales se basan en al menos tres
componentes, la comprensión de los diagramas de fase ternaria es de gran
importancia práctica. Las ideas que se han desarrollado para sistemas binarios se
pueden extender a sistemas con tres o más components4.
52 Transformaciones de fase en metales y aleaciones
El triángulo de Gibbs.
X A+ X B+ X C= 1 (1,62)
La energía libre de Gibbs de cualquier fase ahora se puede representar por una
distancia vertical desde el punto en el triángulo de Gibbs. Si esto se hace para todas
las composiciones posibles, los puntos trazan las superficies de energía libre para
todas las fases posibles, como se muestra en la Fig. 1.41a. Los potenciales químicos
de A, B y C en cualquier fase vienen dados por los puntos donde el plano tangencial
a las superficies de energía libre intersecta los ejes A, B y C. La figura 1.41a se dibuja
para un sistema en el que los tres sistemas binarios AB, BC y CA son eutécticos
simples. Existen superficies de energía libre para tres fases sólidas α, β y γ y la fase
líquida, L. A esta temperatura, la fase líquida es más estable para todas las
composiciones de aleación. A temperaturas más bajas, la superficie de GL se mueve
hacia arriba y eventualmente se cruza con la superficie de Gα como se muestra en la
Fig. 1.41b. Las aleaciones con composiciones en las proximidades de la intersección
de las dos curvas consisten en α + L en equilibrio. Para que los potenciales químicos
sean iguales en ambas fases, las composiciones de los dos phases en equilibrio deben
estar dadas por puntos conectados.
Termodinámica y diagramas de fases 53
U
n
Figura 1.41
a) Energías libres de un líquido y tres fases sólidas de un sistema ternario. (b) Una construcción de plano
tangencial a las superficies de energía libre define el equilibrio entre s y l en el sistema ternario, (c) la
sección isotérmica a través de un diagrama de fase ternario obtade esta manera con una región bifásica
(L1S) y varias líneas de enlace. Las cantidades de l y s en el punto x están determinadas por la regla de la
palanca. (Después de P. Haasen, Metalurgia física, Cambridge University Press, Cambridge, 1978.)
por un plano tangencial común, por ejemplo s y l en la Fig. 1.41b. Estos puntos se
pueden marcar en una sección isotérmica del diagrama de fases de equilibrio como
se muestra en la Fig. 1.41c. Las líneas que unen las composiciones en equilibrio se
conocen como líneas de empate. Al rodar el plano tangencial sobre las dos
superficies de energía libre, se generará toda una serie de líneas de enlace, como pr
y qt, y la región cubierta por estas líneas de enlace pqtr es una región de dos fases en
54 Transformaciones de fase en metales y aleaciones
el diagrama de fases. Por lo tanto, una aleación con composición x en la Fig. 1.41c
minimizará su energía libre al separarse en α sólido con composición s y líquido con
Termodinámica y diagramas de fases 55
d) T= 3 e) T= 2 f) y2 >T >E
g) T= E (h) E>T
Figura 1.42
Secciones isotérmicas a través de la Fig. 1.44. (Después de A. Prince, Alloy Phase Equilibria, Elsevier,
Ámsterdam, 1966.)
1
y
y
Termodinámica y diagramas de fases 57
Figura 1.43
Proyección de las superficies líquidas de la Fig. 1.44 sobre el triángulo de Gibbs.
solo muestra las fases que existen en equilibrio a diferentes temperaturas y no sus
composiciones. Por lo tanto, no se pueden usar como diagramas de fases binarias, a
pesar de la semejanza superficial.
L+a L+b
e L+a+b
n
U
+
n +
L+a+c
C C
a+b+c
1 t
Figura 1.45
Una sección vertical entre los puntos 1, 2 y X de la Fig. 1.44. (Después de A. Prince, Equilibrios de fase
de aleación,
Elsevier, Ámsterdam, 1966.) y que la pendiente de la curva libre de composición
dG = μ μB − A (1,64)
dX B 1
Sustituir esta expresión en la Ecuación 1.63 y multiplicar por X AXB conduce a las
siguientes igualdades:
d G
−X AdμA = X BdμB = X XA B d X2 2 dX B (1.65)
60 Transformaciones de fase en metales y aleaciones
–dm U
n
– µU
dec n
ím
d B
µ –µ = d
d
0 B
1
B
A
Figura 1.46
Evaluación del cambio en el potencial químico debido a un cambio en la composición. (Después de M.
Hillert, en Lectures on the Theory of Phase Transformations, H.I. Aaronson (Ed.), ©The American
Society for Metals y The Metallurgical Society of AIME, Nueva York, 1969.)
que son las ecuaciones requeridas que relacionan dμA, dμB y dX B. La primera igualdad
en esta ecuación se conoce como la relación de Gibbs-Duhem para una solución
binaria. Tenga en cuenta que el subíndice B se ha eliminado de d2G / dX 2 como d 2G
/ dXB2 ≡ d2G / dXA2 . Para una solución regular, la diferenciación de la ecuación 1.39
da
d 2G RT
2 =XX A B − Ω2 (1,66) dX
ddX2G2 = X XRTA B
(1,67)
Por lo tanto
Termodinámica y diagramas de fases 61
Se puede derivar una relación similar para dmA/ dX B. Por lo tanto, la ecuación 1.65
se convierte en
d 2G
= RT
d lnγA d lnγB
X XA B 2 1+ = RT 1+ (1.71)
dX
d ln X A d ln X B
Como se describió anteriormente, si se conoce la energía de Gibbs, otras propiedades
termodinámicas como la entalpía y la entropía podrían derivarse fácilmente. La
información termodinámica recopilada mediante la realización de experimentos
generalmente se compara con los modelos físicos existentes cuyo propósito es
cuantificar y dar predicciones de las propiedades termodinámicas medidas. Sin
embargo, es importante darse cuenta de que tales modelos podrían producir
predicciones poco confiables si se utilizan fuera del rango sobre el cual fueron
verificados.
3
.9 Cálculo de diagramas de fases
La determinación de la energía de Gibbs de fases, con diversas propiedades,
sustancias estequiométricas puras, por ejemplo, en función de la temperatura, la
composición y la presión es necesaria para establecer un banco de datos
termodinámico.
62 Transformaciones de fase en metales y aleaciones
H = H ref
+
∫C p dT (1.72)
Referencia T
S=S
ref + T∫ CTp dT (1.73)
Referencia T
Θ
C=D
T (1.74)
h ωD
Θ= (1.75)
2pk
Dónde
h y k son las constantes de Planck y Boltzmann, respectivamente ω D
se define como la frecuencia de Debye
2 + C42 (1.76) Cp = C1 + C T2
+ C T3 T
G = C1 + C T2+ C T3 ln ( )T + C T4
2 +
+ C T53 CT6 (1.77)
I II III
24
0 1000 6000
Temperatura (K)
Figura 1.47
64 Transformaciones de fase en metales y aleaciones
Capacidad calorífica del gas oxígeno frente a la temperatura. 5 Tenga en cuenta el uso de varios polinomios
en diferentes rangos de temperatura para adaptarse mejor a los datos medidos. (Después de Hack, K., en
The SGTE Casebook: Thermodynamics at Work, Hack, K. (Ed.), The Institute of Materials, Londres,
1996.)
Las contribuciones a la energía de Gibbs de fases estequiométricas puras, como
elementos metálicos u óxidos estequiométricos, SiO2 , por ejemplo, incluyen las que
surgen de la red, las contribuciones magnéticas y de presión.
donde la red G es una función de ΔH y S bajo las condiciones estándar, así como el
polinomio que representa Cp. Otras contribuciones a la energía de Gibbs se discuten
en la siguiente sección.
Al aplicar el esquema anterior a sustancias que exhiben transiciones de fase
magnética de segundo orden, donde la forma de la curva de capacidad calorífica
muestra un cambio bastante repentino,se utilizan muchos polinomios alrededor de la
temperatura de Curie8. Sin embargo, este enfoque tiene limitaciones numéricas dado
el gran número de parámetros de ajuste necesarios, y también sus grandes valores.
Inden sugirió que los cambios excepcionales en Cp deberían tratarse por separado9,10.
Para el caso de una transición de fase magnética, la contribución magnética a la
energía de Gibbs se puede escribir como
T β+ 1) (1.79)
Gmagnético = RTf Tc ln (
Dónde
f es una función de la temperatura dependiente de la estructura
Tc es una temperatura crítica que podría ser la temperatura de Curie o Néel
β es el momento magnético por átomo en la red
Tenga en cuenta que f diferirá en los rangos de temperatura por debajo y por encima
de la temperatura crítica .
Para una descripción completa de la energía libre de Gibbs, también es necesario
considerar la contribución causada por el efecto de la presión sobre el volumen
molar. El consorcio SGTE5 siguió el enfoque sugerido por Murnaghan 11
T
Termodinámica y diagramas de fases 65
Dónde
V0 es el volumen molar a temperatura ambiente α
(T) es la expansión térmica de ese volumen
n es la derivada de presión del módulo a granel que es el inverso de la
compresibilidad
1980 K
GPa
(en
757 K
GPa
(en
500
0 10 20 30 40 50
Presión (GPa)
Figura 1.48
Diagrama de fase presión-temperatura del hierro, calculado. (After Hack, K., en The SGTE Casebook:
Thermodynamics at Work, Hack, K. (Ed.), The Institute of Materials, Londres, 1996.)
A
α( )T = A0 + A T1 + A T2 2 +
3 (1.82)
2
T
Un ejemplo de un diagrama de fase presión-temperatura calculado se muestra en la
Fig. 1.48.
GABex = X XA B
L X −X y (1.84)
Por lo tanto, al vincular los datos experimentales de un binario A-B a esta relación,
se pueden determinar los parámetros. Cuando i = 0, el exceso de energía es X X LA B
(0)
AB , un término que es el del modelo de solución regular, véase la ecuación 1.36.
Aquí el parámetro L(0)AB es independiente de la composición química. Otros
parámetros L( )ABi, cuando i > 0, describen la dependencia de la composición del
parámetro LAB. El parámetro LAB es un parámetro energy que se denomina energía de
interacción. Así que el exceso de energía de Gibbs para una solución ternaria A-B-C
puede escribirse como una primera
aproximación:
GABCex = X XA B L X −X + − i
(1.85)
Termodinámica y diagramas de fases 67
+
X XB Ci
Para los sistemas cuaternarios, a menudo no hay información sobre los cuatro
sistemas ternarios y, a veces, los datos del sistema cuaternario en sí no están
disponibles. Sin embargo, un término cuaternario X X X X LA B C D (0)ABCD puede
añadirse al exceso de Gibbs energy en este caso.
Para otras fases con diferentes propiedades, como las soluciones intersticiales, se
hace referencia al lector, por ejemplo, Hack 5, para obtener más información.
∆ U
n
Figura 1.49
68 Transformaciones de fase en metales y aleaciones
Transformaciones del estado inicial al final a través de un estado activado de mayor energía libre.
∆Ga
∝exp −
tasa kT
exp −∆ Ha
tasa ∝ RT (1.72)
Referencias
1. D.R. Gaskell, Introduction to Metallurgical Thermodynamics, McGraw-Hill, Nueva
York, 1973, p. 342.
2. R.W. Cahn (Ed.), Metalurgia Física, 2ª edn., Capítulos 4 y 5, NorthHolland,
Ámsterdam, 1974.
3. M. Ferrante y R.D. Doherty, Acta Metall. , 27: 1603,1979.
4. M. Hillert, Transformaciones de fase, Capítulo 5, Sociedad Americana de Metales,
Metals Park, OH, 1970.
5. K. Hack (Ed.), The SGTE Casebook: Thermodynamics at Work, Part 1, Institute of
Materials, Londres, 1996.
6. C. Gerthsen y H.O. Kneser, Physik, Springer Verlag, Berlín, Heidelberg, Nueva York,
1969.
7. K.K. Kelley, Boletín 476 de la Oficina de Minas de los Estados Unidos, 1949.
8. I. Barin y O. Knache, Springer-Verlag, Berlín, Heidelberg y Verlag Stahleisen,
Düsseldorf, 1973.
9. G. Inden, Proc. de Calphad V, Düsseldorf, III.4, 1-13, 1976.
10. G. Inden, Proc. de Calphad V, Düsseldorf, IV.1, 1-33, 1976.
11. F.D. Murnaghan, Proc. Natl. Acad. Sci. (Estados Unidos), 30: 244, 1944.
12. O. Redlich y A.T. Kister, Ind. Ing. Chem. , 40: 345-348, 1948.
Lecturas adicionales
A.H. Cottrell, An Introduction to Metallurgy, Capítulo 14 'Aleaciones' y Capítulo 15. 'El
diagrama de fases', Edward Arnold, Londres, 1967.
D.R. Gaskell, Introducción a la termodinámica metalúrgica , McGraw-Hill, Nueva
York, 1973.
P. Gordon, Principios de diagramas de fase en sistemas de materiales, McGraw-Hill,
Nuevo
York, 1968.
M. Hillert, 'Calculation of Phase Equilibria', Capítulo 5 en Phase Transformations,
American Society for Metals, Ohio, 1970.
M. Hillert, 'The uses of the Gibbs free energy-composition diagrams', Capítulo 1 en Lectures
on the Theory of Phase Transformations, H.I. Aaronson (Ed.) publicado por The
Metallurgical Society of AIME, Nueva York, 1975. A. Prince, Alloy Phase Equilibria,
Elsevier, Londres, 1966.
G.V. Raynor, 'Phase diagrams and their determination', Capítulo 7 en R.W. Cahn (Ed.)
Physical Metallurgy, North-Holland, 1970.
70 Transformaciones de fase en metales y aleaciones
Ejercicios
1.1 El calor específico del cobre sólido por encima de 300 K es given por
B
B
1 2
(d
1
µ2
(c )
2
µ1
U B
n
2
U 2 1 B
n (e) f)
Figura 2.1
La energía libre y los cambios de potencial químico durante la difusión, (a) y (b) difusión 'cuesta abajo',
(c) y (d) difusión 'cuesta arriba', (e) μ4A > μ A1 por lo tanto, los átomos A se mueven de (2) a (1), μB2 > μ
4
.1 Mecanismos atómicos de difusión
Hay dos mecanismos comunes por los cuales los átomos pueden difundirse a través
de un sólido y el mecanismo operativo depende del tipo de sitio ocupado en la red.
Sustitutivo átomos generalmente difusos por un vacante mecanismo
B
2
por lo tanto, los átomos B se mueven de (1) a (2). f) μA1 > μA2, por lo tanto, los átomos de E A se mueven
de (1) a (2), μB2 > μB1, por lo tanto, los átomos de B se mueven de (2) a (l).
mucho más fácil de medir que las diferencias de potencial químico, sin embargo, es
más conveniente relacionar la difusión con los gradientes de concentración. Por lo
tanto, el resto de este capítulo se ocupará principalmente de este enfoque de la
difusión.
2
( ) (b)
Figura 2.2
Movimiento de un átomo en una vacante adyacente en una red fcc, (a) Un plano cerrado, (b) Una celda
unitaria que muestra los cuatro átomos (sombreados) que deben moverse antes de que pueda ocurrir el
salto. (Después de P.G. Shewmon, Diffusion in Solids, McGraw-Hill, Nueva York, 1963.)
Losátomos intersticiales más pequeños migran forzando su camino entre los átomos
más grandes, es decir, intersticialmente.
Normalmente, un átomo de sustitución en un cristal oscila alrededor de un sitio
dado y está rodeado por átomos vecinos en sitios similares. La energía racional vib
mediaposeída por cada átomo está dada por 3 kT, y por lo tanto aumenta en
proporción a la temperatura absoluta. Dado que la frecuencia media de vibración es
aproximadamente constante, la energía vibratoria aumenta al aumentar la amplitud
de lasoscilaciones. Normalmente, el movimiento de un átomo sustitutivo está
limitado por sus vecinos y el átomo no puede moverse a otro sitio. Sin embargo, si
un sitio adyacente está vacante, puede suceder que una oscilación particularmente
violenta resulte en que el átomo salte a la vacante. Esto se ilustra en la Fig. 2.2. Tenga
en cuenta que para que se produzca el salto, los átomos sombreados en la Fig. 2.2b
deben separarse para crear suficiente espacio para que el átomo migratorio pase entre
ellos. Por lo tanto, la probabilidad de que cualquier átomo sea unble para saltar a un
sitio vacante depende de la probabilidad de que pueda adquirir suficiente energía
vibratoria. La velocidad a la que cualquier átomo dado es capaz de migrar a través
del sólido estará claramente determinada por la frecuencia con la que se encuentre
con una vacante y esto a su vez depende de la concentración de vacantes en el sólido.
Se demostrará que tanto la probabilidad de saltar como la concentración de vacantes
son extremadamente sensibles a la temperatura.
Cuando un átomo de soluto es apreciablemente más pequeño en diámetro que el
disolvente, ocupa uno de los sitios intersticiales entre los átomos del disolvente. En
los materiales fcc los sitios intersticiales están a medio camino a lo largo de los
bordes del cubo o, equivalentemente, en el centro de la celda unitaria, Fig. 2.3a. Estos
se conocen como sitios octaédricos ya que los seis átomos alrededor del sitio forman
un octaedro. En la red bcc, los átomos intersticiales también ocupan a menudo los
sitios octaédricos que ahora se encuentran en posiciones de centrado de bordes o
centrados en la cara, como se muestra en la Fig. 2.3b.
(a) Intersticios octaédricos (0) en un cristal fcc, (b) intersticios octaédricos en un cristal bcc. (Después
P. Haasen, Metalurgia Física, Cambridge University Press, Cambridge, 1978.)
(U
n
U
n
B
∂ B
U
∂
n
(b
Figura 2.5
Difusión intersticial por saltos aleatorios en un gradiente de concentración.
1 s (J) estará dado por:
1
JB = 6 ΓBn1 átomos m −2 s−1 (2.1)
Durante el mismo tiempo, el número de átomos que saltan del plano (2) al (1),
asumiendo que GB es independiente de la concentración, viene dado por:
Dado que n1 > n2 habrá un flujo neto de átomos de izquierda a derecha dado por:
1 ΓB α 2 ∂∂CxB átomos m
s JB = −
−2 −1
D = 1 C a2 (2.3)
B 6 B rinde:
C
∂ B (2.4)
J=−D
B B ∂x
(U) b
n
(c
Figura 2.6
Átomo intersticial, (a) en posición de equilibrio, (b) en la posición de máxima distorsión de la red, (c)
Variación de la energía libre de la red en función de la posición de intersticial. (Después de P.G. Shewmon,
en Physical Metallurgy, 2ª ed., R.W. Cahn (Ed.), Holanda del Norte, Ámsterdam, 1974.)
Gm
−∆ (2,7)
Γ = zv exp
B RT
ΔGm puede considerarse la suma de una gran entalpía de activación ΔHm y un término
de entropía de activación pequeña –TΔsm.
Combinando esta expresión con la Ecuación 2.3 se obtiene el coeficiente de
difusión como
1 ∆Sm −∆H
DB = α2zv exp
6 R exp RT m (2.8)
Dónde
Q ID = ∆Hm (2.11)
TABle 2.1
Expe rimental Diffusion Data for Interstitials in Ferritic (bcc) Iron
Soluto D0/mm2 s–1 Q/kJ mol−1 Ref.
C 2.0 84.1 2
N 0.3 76.1 3
H 0.1 13.4 4
0
1/T
Figura 2.7
La pendiente del log D v. 1/ T da la energía de activación para la difusión Q.
la difusión aumenta a medida que aumenta el tamaño del átomo intersticial. (Los
diámetros atómicos disminuyen en el orden C, N, H.) Esto es de esperar ya que los
átomos más pequeños causan menos distorsión de la red durantela m igración.
Se puede obtener una representación gráfica conveniente de D en función de la
temperatura escribiendo la Ecuación 2.9 en la forma
Q 1
log D = log D0 −
2 3. T (2.12)
Por lo tanto, si el log D se traza contra (1/ T) se obtiene una línea recta con una
pendiente igual a −(Q/ 2.3 R) y una intercepción en el eje log D en el log D0, ver Fig.
2.7.
(U) 0 +d
n
1
2
(b) 0 1 2
D
1 2
(c)
Figura 2.8
La derivación de la segunda ley de Fick.
donde l es el grosor de la pared. En este básico el flujo a través de la pared viene dado
por
J = D CHH (2.13)
H l
(2.15)
Día2 = D1 + ∂∂xJδx
∂CB = ∂ D ∂CB
∂t ∂x B ∂x (2.17)
que se conoce como la segunda ley de Fick. Si se pueden ignorar las variaciones de
DB con concentración, esta ecuación se puede simplificar a
Figura 2.9
(a) ∂2C/ ∂x2 > 0 todas las concentraciones aumentan con el tiempo, (b) ∂ 2C/ ∂x2 < 0 todas las
concentraciones disminuyen con el tiempo.
=0 l B0
=T
Figura 2.10
El efecto de la difusión sobre una variación sinusoidal de la composición.
Homogeneización
A menudo es de interés poder calcular el tiempo necesario para que una aleación no
homogénea alcance la homogeneidad completa, como por ejemplo en la eliminación
de la segregación en las piezas fundidas.
La variación de composición más simple que se puede resolver matemáticamente
es if CB varía sinusoidalmente con la distancia en una dimensión como se muestra en
la Fig. 2.10. En este caso , los átomos B se difunden por los gradientes de
concentración, y las regiones con curvatura negativa, como entre x = 0 y x = l,
disminuyen en concentración, mientras que las regiones que sereducen a x = l y 2l
aumentan en concentración. La curvatura es cero en x = 0, l, 2l, por lo que las
concentraciones en estos puntos permanecen sin cambios con el tiempo. En
consecuencia, el perfil de concentración después de un cierto tiempo se reduce al
indicado por la línea discontinua in Fig. 2.10. En el tiempo t = 0 la concentración
profi le viene dada por
πx
C = C +β0 sin l (2.19)
_ donde C es la composición media, y β0 es la amplitud de la concentración
inicial profi le. Suponiendo que DB es independiente de la concentración, la solución
de la ecuación 2.18 que satisface esta condición inicial es
πx −t
(2.20)
C = C +β0 sin l
exp τ
τ= 2 2 l (2.21)
π DB
Así, la amplitud del perfil de concentración después de un tiempo t(β) viene dada por
C en x = l/2, es decir.
−t
β β= exp (2.22)
0 t
En otras palabras, la amplitud de la concentración profi le
disminuyeexponencialmente con el tiempo y después de un tiempo suficientemente
largo se acerca a cero de modo que C = C en todas partes. La velocidad a la que esto
ocurre está determinada por el tiempo de relajación t. Después de un tiempo t = t,β
β= 0/ ,e es decir, la amplitud ha disminuido a 1/ 2.72 de su valor en t = 0. Por lo tanto,
la distribución del soluto en esta etapa aparecería como lo muestra la línea
discontinua en la Fig. 2.10. Después de un tiempo t = 2t la amplitud se reduce en un
total de 1/e2, es decir, en aproximadamente un orden de magnitud. De la ecuación
2.21 se puede ver que la tasa de homogeneización aumenta rápidamente a medida
que disminuye la longitud de onda de las fl uctuaciones.
El perfil de concentración inicial no suele ser sinusoidal, pero en general cualquier
perfil de concentración puede considerarse como la suma de una serie infinita de
ondas sinusoidales de longitud de onda y amplitud variables, y cada onda decae a
una velocidad determinada por su propia τ. Por lo tanto, los términos de longitud de
onda corta mueren muy rápidamente y la homogeneización será determinada en
última instancia por τ para el componente de longitud de onda más largo.
C = Cs −(Cs − C0)erf
2√(Dt) (2.23)
s
1
2
½ s + 0)
3
1 0
Figura 2.11
Perfiles de concentración en momentos sucesivos (t3 > t2 > t1) para la difusión en una barra semiinfinita
cuando la concentración superficial Cs se mantiene constante.
donde 'erf' significa función de error que es una integral indefinida definida por la
ecuación
∫
2
P 0
con patio (
)con =exp(−y 2)dy
patio )
( +1
+0.5 =0
+0.5
–1
(U) (b 0
n
Figura 2.12
a) Diagrama esquemático que ilustra las principales características de la función de error. b) Perfiles de
concentración en momentos sucesivos (t2 > t1 > 0) cuando se recocen dos barras semiinfinitas de diferente
composición después de la soldadura.
Por lo tanto, el grosor de la capa carburada es ~ (√ Dt). Nótese también que la
profundidad de cualquier línea de isoconcentración es directamente proporcional a
√(Dt), es decir, para obtener un aumento de dos vecesen la penetración se requiere
un aumento de cuatro veces en el tiempo.
Para el caso de la difusión de carbono en austenita a 1000°C, Do 4 × 10−11 m s5 −1,
lo que significa que una capa carburada de 0,2 mm de espesor requiere un tiempo de
( .0 2 × 10−3 6) / 4 × 10−11, es decir, 1000 s (17 min).
5
.3 Difuso sustitutivoen
La difusión en aleaciones intersticiales diluidas fue relativamente simple porque los
átomos difusores siempre están rodeados de sitios "vacíos" a los que pueden saltar
siempre que tengan suficiente energía para superar la barrera energética para la
migración. En dif sustitutivoSin embargo, un átomo solo puede saltar si hay un sitio
vacante en una de las posiciones de red adyacentes como se muestra en la Fig. 2.2.
El caso más simple de difusión sustitutiva es la autodifusión de átomos en un metal
puro. Esto es susceptible de un simple modelo atómico similar al caso de la difusión
intersticial y será tratado en primer lugar. La difusión sustitutiva en aleaciones
binarias es más compleja y se tratará por separado.
6
.3.1 Autodifusión
Hay otras situaciones en las que la solución a la ecuación de difusión es muy
similar a la ecuación 2.23. Por ejemplo, durante la descarburación del acero, la
concentración superficial se reduce a un valor muy bajo y el carbono se difunde fuera
de la muestra. El perfil de carbono viene dado por
C = C0 erf 2√(Dt)
D
* = DA = 61α2Γ (2,27)
Un
G
−∆ m (2,28)
Γ =υz X v exp RT
G
e = exp −∆ v (2.29)
Xv RT
Combinando estas tres últimas ecuaciones se obtiene
1 −(∆G + ∆G )
DA = 6 α υ2zexp RTm v (2.30)
Sustituyendo ΔG = ΔH – TΔS da
−QSD
DA = D0 exp RT (2,32)
Dónde
R y
La ecuación 2.32 es la misma que se obtuvo para la difusión intersticial, excepto que
la energía de activación para la autodifusión tiene un término adicional (ΔHv). Esto
se debe a que la autodifusión requiere la presencia de vacancies cuya concentración
depende de ΔHv.
Algunos de los datos experimentales sobre la autodifusión sustitutiva se resumen
en la Tabla 2.2. Se puede ver que para una estructura cristalina dada y el tipo de
enlace Q / RTm es aproximadamente constante; es decir, la entalpía de activación para
la autodifusión, Q, es aproximadamente proporcional a la temperatura de fusión de
equilibrio, Tm. Además, dentro de cada clase, la difusividad a la temperatura de
fusión, D(Tm) y D0 son aproximadamente constantes. Por ejemplo, para la mayoría
de los metales empaquetados (fcc y hcp) Q/RTm ~ 18 y D(Tm) 1 μm –2 s–1 (10–12 s–1).
Los datos Q/RTm y D(Tm) también se trazan en la Fig. 2.13 junto con los datos de
otros materiales para la comparación. Una consecuencia inmediata de estas
correlaciones es que los coeficientes de difusión de todos los materiales con una
estructura cristalina y un tipo de enlace dados serán aproximadamenteiguales a la
misma fracción de su temperatura de fusión, es decir, D(T/Tm) = constante. (T/Tm se
conoce como la temperatura homóloga.)
Las correlaciones anteriores se han evaluado para la presión atmosférica. Sin
embargo, hay datos experimentales limitados quesuponen que las mismas
correlaciones se mantienen independientemente de la presión, siempre que, por
supuesto, se tenga en cuenta el efecto de la presión sobre Tm . Dado que el volumen
generalmente aumenta al fundirse, elevar la presión aumenta Tm y, por lo tanto,
disminuye la difusividad en una temperatura de temple dada.
Que exista una correlación aproximada entre Q y Tm no es sorprendente: aumentar
la fuerza del enlace interatómico hace que el proceso de fusión sea más difícil; es
decir, Tm se eleva. También hace que la difusión sea más difícil al aumentar
ΔHv y ΔHm.
Consider el efecto de la temperatura sobre la autodifusión en Cu como ejemplo. A
800°C (1073 K) los datos de la Tabla 2.2 dan Dcu = 5 × 10–9 mm 2 s–1. La distancia
de salto α en Cu es de 0,25 nm y, por lo tanto, la ecuación 2,3 da ΓCu = 5 × 10–5 saltos
s−1. Después de una hora a esta temperatura, √ (Dt) ~ 4μm. Extrapolando el
TABle 2.2
Datos experimentales de presión para la autodifusión sustitutiva en metales puros en la
atmósfera
Tm D0 Q Q A (Lm )
−1 µE 2 s −1
2 s −1 RT
65.7
Metal m
(continuación)
TABle 2.2 (C O NTiNueD )
Datos experimentales de presión para la autodifusión sustitutiva en metales puros en la
atmósfera
Tm D0 Q Q A (Lm )
−1 µE 2 s −1
2 s −1 RTm
Metal β-Co γ-Fe†
Pd
1825
20,5
266,3
17,6
0,49
Th 2023 120
Pt
2046 22
Con
Gráfico
Hie
lo
Figura 2.13
(a) Difusividades del punto de fusión y (b) energías de activación normalizadas
para varias clases de materiales. (Después de A.M. Brown y M.F. Ashby, Acta
Metallurgica 28 (1980) 1085.
Hie
es lo
o
Cc Tr
cco raro
10 15 20 25 30 35
(b) Energía de activación normalizada, Q/RTm
Figura 2.13
(Continuación)
Los datos a 20 ° C, sin embargo, dan DCu, ~ 10–34, es decir, Γ, 10–20 saltos s−1.
Alternativamente, ¡cada átomo daría un salto cada 10 12 años!
Experimentalmente, el método habitual para determinar el coeficiente de
autodifusión es depositar una cantidad conocida (M) de un isótopo radiactivo A* en
los extremos de dos barras de A que luego se unen como se muestra en la Fig. 2.14a.
Después de recocir durante un tiempo conocido a una temperatura fija, A* se habrá
difundido en A y el perfil de concentración se puede determinar mecanizando capas
delgadas de la barra y midiendo la radiactividad en función de la posición. Dado que
A y A* son químicamente idénticos, la difusión de A* en A ocurrirá de acuerdo con
la Ecuación 2.18. La solución de esta ecuación para las condiciones de contorno
actuales es
M exp
=
x
C 2√(πDt) − 4 Dt2 (2.35)
M tiene unidades [cantidad m−2] y C [cantidad m−3]. La Figura 2.14b muestra la forma
de esta ecuación ajustada a puntos experimentales para la autodifusión en oro.
2.3.2 Difusión de vacantes
El salto de átomos a sitios vacantes puede considerarse igualmente como el salto de
vacantes a sitios de átomos. Si se introducen vacantes excesivas en
*
(U)
n
(b) Distancia, x, mm
Figura 2.14
Ilustración del principio de difusión del trazador y del método de la fuente plana para determinar el
coeficiente de autodifusión del oro, (a) acoplamiento de difusión inicial con la fuente plana de oro
radiactivo Au*. b) Distribuciónion de Au* después de la difusión durante 100 h a 920 °C. (Después de
A.G. Guy, Introducción a la Ciencia de los Materiales, McGraw-Hill, Nueva York, 1971.)
D = 1a2C (2.36)
V 6 v
1 ∆ S −∆ H
DV = 6 α υ2zexp R m exp RT m (2.37)
En este caso , ΔHm y ΔSm se aplican a la migración de una vacante y, por lo tanto,
son los mismos que para la migración de un átomo de sustitución. Comparación de
ecuaciones
2.37 y 2.31 se puede ver que
Dv = DA/ X ve (2.38)
CB
∂ (2.40)
J = −D
B B ∂x
donde JA y JB son las fl uxes de los átomos A y B a través de un plano de red dado.
Este punto no necesitaba enfatizarse en el caso de la difusión intersticial porque los
planos de red de los átomos padres no se vieron afectados por el proceso de difusión.
Sinembargo, se verá que la situación es diferente en el caso de la difusión sustitutiva.
Para derivar la segunda ley de Fick consideremos la interdifusión de átomos A y
B en una pareja de difusión que se hace soldando bloques de A yd B puros como se
muestra en la Fig. 2.15a. Si la pareja se recoce a una temperatura lo suficientemente
alta, se desarrollará una producción de concentración como se muestra.
Si hacemos la suposición simplificadora de que el número total de átomos por
unidad de volumen es una constante, C0, independiente de la composición, entonces
C0 = CA + CB (2.41)
A
B
(2.42)
x x
Por lo tanto, en una posición dada, los gradientes de concentración que impulsan la
difusión de los átomos A y B son iguales pero opuestos, y las fluxes de A y B en
relación con
U
n
E B
n P
1 0
U U B B
n n
0 1
(U)
n
+ U
n
0
E
n
B
–
(b
∂ E=– E
∂ n n
+
0
(c Creado
Figura 2.15
Interdifusión y flujo de vacantes. a) Perfil de composición después de la interdifusión de A y B. b) Los
flujos correspondientes de átomos y vacantes en función de la posición x. c) La tasa a la que aumentaría
o disminuiría la concentración de vacantes si no se crearan o destruyeran vacantes por la escalada de
dislocación. la celosía se puede escribir como
C
∂ A (2.43)
J = −D
A A ∂x
∂CA
JB = DB ∂x
Átomos
Figura 2.16
El salto de átomos en una dirección puede considerarse como el salto de vacantes en la otra dirección.
Estos flujos se muestran esquemáticamente en la Fig. 2.15 para el caso DA > DB, es
decir, JA| >| JB| .
Cuando los átomos migran por el proceso de vacante, el salto de un átomo a un
sitio vacante puede considerarse igualmente como el salto de la vacante sobre el
átomo, como se ilustra en la Fig. 2.16. En otras palabras, si hay un flujo neto de
átomos en una dirección, hay un flujo igual de vacantes en ladirección e opuesta. Así,
en la Fig. 2.15a hay una fl ux de vacantes −JA debido a la migración de átomos A
más una fl ux de vacantes −JB debido a la difusión de átomos B. Como JA > JB habrá
un flujo neto de vacantes.
Jv = −JA − JB (2.44)
C
∂ A (2.45)
J = (D −D)
V A B ∂x
Esto conduce a una variación en Jv a través de la pareja de difusión como se ilustra
en la Fig. 2.15b.
⊥
⊥
(a) b) c)
Figura 2.17
(a) antes, (b) después: una vacante se absorbe en un trote en una dislocación de borde (ascenso positivo),
(b) antes, (a) después: una vacante se crea por ascenso negativo de una dislocación de borde, (c) Dibujo
en perspectiva de una dislocación de borde trotado.
=0
Figura 2.18
Un flujo de vacantes hace que los planos atómicos se muevan a través de la muestra.
en la Fig. 2.18. Esto significa que se introducirán planos atómicos adicionales en el
lado rico en B, mientras que planos enteros de átomos serán "comidos" en el lado
rico en A. En consecuencia, los planos de celosía en el medio de la pareja se
desplazarán hacia la izquierda.
La velocidad a la que se mueve cualquier plano de red dado, v, puede estar
relacionada con el flujo de vacantes que lo cruzan. Si el plano tiene un área A, durante
un pequeño intervalo de tiempo δt, el plano barrerá un volumen de Avdt que contiene
átomos de Av dt C0 . Este número de átomos se elimina por el número total de
vacantes que cruzan el plano en el mismo intervalo de tiempo, es decir, J Av ⋅δt, dando
Jv = C en0 (2,46)
Por lo tanto, la velocidad de los planos de la red variará across la pareja de la misma
manera que Jv, ver Fig. 2.15b. La sustitución de la ecuación 2.45 da
X
∂ A (2.47)
v = (D − D )
A B ∂x donde la fracción molar de A, X A =
CA/C0