Parte A Módulo I Puntos 3 y 4-Problema Resuelto

Descargar como pdf o txt
Descargar como pdf o txt
Está en la página 1de 10

MECÁNICA DE FLUIDOS II

1era y 2da Semana II (17 al 22 y 24 al 28 de agosto de 2020)

Dra. Anet de Palma

PARTE A MÓDULO I: Flujo Compresible Isentrópico 1-D con cambio de área y


onda de choque

Formulación de las ecuaciones de:

1. Conservación de masa. 2. Energía

3. Momentum. 4. Balance de entropía en estado estable y en una


dimensión.

 Problema resuelto

3. CANTIDAD DE MOVIMIENTO

¿Qué causa el movimiento de un fluido?


Los fluidos empiezan a moverse cuando sobre ellos se aplica una fuerza resultante distinta
de cero. Ejemplos:
 La P en un lugar es mayor que en otro, el fluido tendera a moverse hacia la región
de menor presión.
 La gravedad, ya que los líquidos fluyen cuesta abajo, convirtiendo su energía
potencial en energía cinética.
 Las diferencias de temperatura causaran que una parte del fluido tenga una densidad
menor a diferencia de otra parte del fluido y el fluido más ligero tendera a subir.
 La fricción, cuando una capa de fluido se mueve con respecto a una capa adyacente,
se desarrolla un esfuerzo viscoso tangente que hace que el flujo se mueva más
rápido o lento.

El principio de conservación de la cantidad de movimiento se deduce de esta segunda ley


de Newton:
𝐹⃗ = 𝑚𝑎⃗
F.- fuerza m.-masa a.- aceleración.
La cantidad de movimiento es el producto de la masa de una partícula y su vector velocidad
o también se le llama a este mismo producto de la masa y de la velocidad de un cuerpo
como momento lineal.

La cantidad de movimiento de un cuerpo rígido de masa “m” que avanza con una velocidad
“𝑉̅ ” es “ 𝑚𝑉̅ ”, entonces, la segunda ley de newton también se puede expresar como que
“la razón de cambio de la cantidad de movimiento de un cuerpo es igual a la fuerza neta
que actúa sobre él”. Por lo tanto, en la mecánica de fluidos suele hacerse referencia a la
segunda ley de Newton como la ecuación de momento lineal.

1
𝑑𝑀
𝐹=
𝑑𝑡

El momento lineal o cantidad de movimiento es el producto de la masa y la velocidad, y su


dirección es la misma que la dirección del vector velocidad

Principio de conservación de la cantidad de movimiento


La cantidad de movimiento de un sistema se mantiene constante cuando la fuerza neta que
actúa sobre él es cero y, por lo tanto, la cantidad de movimiento de esos sistemas de
conserva.

La cantidad de movimiento es una cantidad vectorial que consta de magnitud, dirección y


sentido, y en consecuencia esta tiene la dirección de la velocidad. Cualquier ecuación
vectorial se puede escribir en forma escalar, para una dirección específica, con el uso de
magnitudes; por ejemplo
𝑑(𝑚𝑉𝑥 )
𝐹 = 𝑚𝑎𝑥 =
𝑑𝑡
Si se considera ahora un escurrimiento permanente con gasto Q y se escogen dos secciones,
1 y 2, de dicho escurrimiento, la masa que fluye por cualquiera de ellas en un tiempo ∆t, es:
𝛾𝑄
𝑚= ∆𝑡
𝑔
𝑁 𝑚3 𝑚 𝑚
𝛾 = 3; 𝑄 = ; 𝑡 = 𝑠; 𝑔 = 2 ; 𝑁 = 𝑘𝑔 ∙ 2
𝑚 𝑠 𝑠 𝑠
Y si ∆𝑉⃗⃗ es la diferencia de velocidades medias en ambas secciones, la 2da ley de Newton
puede escribirse:
𝛾𝑄 𝛾𝑄
𝐹⃗ = 𝑔 ∆𝑉
⃗⃗ o ∑2𝑖=1 𝐹⃗ = ⃗⃗1 − 𝑉
∆(𝑉 ⃗⃗2 )
𝑔

El impulso desarrollado por un escurrimiento en una sección de área hidráulica A durante


un tiempo dt, y asumiendo que 𝑑𝑄 = 𝑉 ⃗⃗ ∙ 𝑑𝐴 se puede escribir

Y si se calcula con la velocidad media para toda la sección A, habrá que corregirlo con coeficiente β
𝛾
de manera que: 𝐼 = 𝛽 𝑔 𝑉𝑚2 𝐴
⃗⃗ 2 𝑑𝐴
ʃ𝐴 𝑉
𝛽=
𝑉𝑚2 𝐴

Flujo de cantidad de movimiento a través de un diferencial de área ρV2dA

2
Flujo total de cantidad de movimiento a través de la sección ʃA ρV2dA
La velocidad varía en los diferentes puntos de la sección transversal, y el ʃA ρV2dA=𝛽𝜌𝑉𝑚2 𝐴
resultado del integral requiere un ajuste para poderlo expresar en términos de la
velocidad media en la sección. El coeficiente que permite igualar las
expresiones,β , se conoce como coeficiente de Boussinesq para la corrección de
la cantidad de movimiento.
Y la expresión para el coeficiente de corrección de Boussinesq es ʃ𝐴 𝑉 2 𝑑𝐴
𝛽=
𝑉𝑚2 𝐴
A partir de este coeficiente se redefine la ecuación de fuerza dinámica F= 𝛽𝜌𝑉𝑚2 𝐴

4) ENTROPÍA Y EL PRINCIPIO DE INCREMENTO DE ENTROPÍA

CONCEPTO DE ENTROPÍA
La desigualdad de Clausius es la base para la definición de una nueva propiedad
llamada entropía. Se considera un ciclo reversible formado por dos procesos internamente
reversibles A y B como se muestra en la figura. Aplicando la
desigualdad de Clausius, se tiene lo siguiente:
𝛿𝑄
∮𝛿 𝑇
=0
2𝐴 𝛿𝑄 1𝐵 𝛿𝑄 2𝐴 𝛿𝑄 2𝐵 𝛿𝑄
∮1𝐴 + ∮2𝐵 =0 ∮1𝐴 − ∮1𝐵 =0
𝑇 𝑇 𝑇 𝑇

2𝐴 𝛿𝑄 2𝐵 𝛿𝑄 Ciclo Reversible
∮1𝐴 = ∮1𝐵
𝑇 𝑇

𝛿𝑄
Puesto que ∫ 𝑇 es la misma para ambas trayectorias entre los estados 1 y 2, el valor de
esta integral depende sólo de los estados extremos y no de la trayectoria seguida. En
consecuencia debe representar el cambio de una propiedad ya que es independiente de la
trayectoria. A esta propiedad se denomina entropía y se designa por S.
𝛿𝑄 𝐾𝐽
𝑑𝑆 = ( ) ( )
𝑇 𝑅𝑒𝑣 𝐾

La entropía por unidad de masa, denominada s, es una propiedad intensiva y se mide con la
unidad kJ/(kg∙K). Si integramos esta ecuación entre los estados inicial y final, obtenemos el
cambio de entropía de un sistema durante un proceso:
2 2
𝛿𝑄
∫ 𝑑𝑆 = ∫ ( )
1 1 𝑇 𝑅𝑒𝑣

3
2 𝛿𝑄
𝑆2 − 𝑆1 = = ∫1 ( 𝑇 )
𝑅𝑒𝑣

Si vamos a realizar la integración de la ecuación debemos conocer la relación entre Q y T


durante un proceso.
CAMBIO DE ENTROPÍA EN EL CICLO DE CARNOT
En la figura se representa el ciclo de Carnot en un
diagrama T-s. Para hacer la evaluación de la entropía Ciclo de Carnot
para cada proceso se comienza por el proceso de
compresión adiabático de 4 a 1; luego con el proceso de
calentamiento isotérmico de 1 a 2; después con una
expansión adiabática de 2 a 3 y finalmente con un
proceso de enfriamiento isotérmico de 3 a 4.
Proceso isotérmico reversible de 1 a 2
1 2 𝑄𝐻
𝑆2 − 𝑆1 = ∫ 𝛿𝑄 =
𝑇𝐻 1 𝑇𝐻

Proceso adiabático reversible de 2 a 3 𝑆3 − 𝑆2 = 0

1 4 𝑄𝐿
Proceso isotérmico reversible de 3 a 4 𝑆4 − 𝑆3 = ∫ 𝛿𝑄 =
𝑇𝐿 3 𝑇𝐿

Proceso adiabático reversible de 4 a 1 𝑆1 − 𝑆4 = 0


Cuando la entropía permanece constante durante el proceso, es decir ΔS=0, el proceso se
llama proceso isoentrópico.

CAMBIO DE ENTROPÍA EN PROCESOS


IRREVERSIBLES
En la figura se considera un ciclo irreversible formado
por dos procesos, uno reversible y otro irreversible

Aplicando la desigualdad de Clausius


𝛿𝑄
∮𝛿 <0
𝑇
Transformando la integral cíclica en la suma de integrales de línea:
2𝐴 𝛿𝑄 1𝐵 𝛿𝑄 2𝐴 𝛿𝑄 2𝐵 𝛿𝑄
∮1𝐴 + ∮2𝐵 <0 ∮1𝐴 < ∮1𝐵
𝑇 𝑇 𝑇 𝑇

4
Tenemos la expresión, corresponde a la variación de entropía para un proceso irreversible.
𝛿𝑄
𝑑𝑆 > ( )
𝑇 𝐼𝑟𝑟𝑒𝑣

PRINCIPIO DE INCREMENTO DE ENTROPÍA


El principio de incremento de entropía se expresa para cualquier proceso como:
ΔS Total = ΔS sistema + ΔS medio exterior ≥ 0

Este principio es aplicable tanto para sistemas cerrados como para volúmenes de control.
La igualdad se cumple para procesos reversibles y la desigualdad para procesos
irreversibles. La variación de entropía para el medio exterior se expresa como:
𝑄𝑖
∆𝑆𝑚𝑒𝑑𝑖𝑜 𝑒𝑥𝑡𝑒𝑟𝑖𝑜𝑟 = ∑
𝑇𝑖

Se tiene Qi es la transferencia de calor a través de la frontera a la temperatura Ti.


Si el sistema se encuentra aislado térmicamente, es decir no hay transferencia de calor con
el medio exterior, el principio de incremento de entropía se expresa como:
ΔS sistema ≥ 0

RELACIONES TdS
Las relaciones TdS permiten determinar la variación de la entropía a través de trayectorias
reversibles. Pero los resultados obtenidos son válidos tanto para procesos reversibles como
irreversibles. Para obtenerlas, se parte de la primera ley y se considera que el proceso es
reversible:
𝛿𝑄𝑖𝑛𝑡 𝑟𝑒𝑣 − 𝛿𝑊𝑖𝑛𝑡 𝑟𝑒𝑣,𝑠𝑎𝑙𝑒 = 𝑑𝑈

𝛿𝑄𝑖𝑛𝑡 𝑟𝑒𝑣 = 𝑇𝑑𝑠; 𝛿𝑊𝑖𝑛𝑡 𝑟𝑒𝑣,𝑠𝑎𝑙𝑒 = 𝑃𝑑𝑉


Sustituyendo en la ecuación de la Primera Ley, tenemos:
TdS = dU + PdV (kJ) (ecuación Tds o Gibbs)
o Tds = du + Pdv (kJ/kg)
Reemplazamos:
h = u +pv → dh = du + Pdv + vdp
Tds = du + Pdv → Tds = dh –vdP
Tds = du + Pdv Tds = dh –vdP Estas ecuaciones son válidas para procesos
reversibles o irreversibles. Se aplican tanto a sistema cerrado o sistema abierto.

Para cambios diferenciales se obtienen relaciones explicitas en la entropía al resolver s las


ecuaciones:

5
𝑑𝑢 𝑃𝑑𝑣 𝑑ℎ 𝑣𝑑𝑃
𝑑𝑠 = + 𝑑𝑠 = −
𝑇 𝑇 𝑇 𝑇

CAMBIO DE ENTROPÍA DE GASES IDEALES.


Vamos a obtener dos expresiones para el cambio de entropía de un gas ideal mediante las
𝑅𝑇
relaciones Tds, sustituyendo: du= cvdT y 𝑃 = 𝑣

𝑑𝑇 𝑑𝑣
𝑑𝑠 = 𝑐𝑉 +𝑅
𝑇 𝑣
El cambio de entropía para un proceso se obtiene integrando esta relación entre los estados
extremos:
2
𝑑𝑇 𝑣2
𝑠2 − 𝑠1 = ∫ 𝑐𝑣 (𝑇) +𝑅
1 𝑇 𝑣1

𝑇 𝑃
Segunda relación de entropía: 𝑠2 − 𝑠1 = 𝑐𝑃 𝐼𝑛 𝑇2 − 𝑅 𝐼𝑛 𝑃2
1 1

CALORES ESPECÍFICOS CONSTANTES (análisis aproximado)


Esta forma simplifica el análisis, pero hay una pérdida de exactitud. Ejemplo: 1) Para los
gases monoatómicos ideales como el helio, los calores
específicos son independientes de la temperatura y por
consiguiente suponer el calor constante específico no provoca
error; 2) Para gases ideales cuyos calores específicos varían
casi linealmente en el intervalo en que se hallan las
temperaturas de interés, el error posible se minimiza usando
los valores de calores específicos evaluados a temperatura
promedio.
𝑇 𝑣 𝑘𝐽
𝑠2 − 𝑠1 = 𝑐𝑣,𝑝𝑟𝑜 𝐼𝑛 𝑇2 + 𝑅 𝐼𝑛 𝑣2 (𝑘𝑔∙𝐾)
1 1

𝑇2 𝑃2 𝑘𝐽
𝑠2 − 𝑠1 = 𝑐𝑃,𝑝𝑟𝑜 𝐼𝑛 − 𝑅 𝐼𝑛 ( )
𝑇1 𝑃1 𝑘𝑔 ∙ 𝐾

Los cambios de entropía también pueden expresarse por unidad de mol, al multiplicar estas
relaciones por la masa molar;
𝑇2 𝑣2 𝑘𝐽
𝑠̅2 − 𝑠̅1 = 𝑐̅𝑣,𝑝𝑟𝑜 𝐼𝑛 + 𝑅𝑢 𝐼𝑛 ( )
𝑇1 𝑣1 𝑘𝑚𝑜𝑙∙𝐾

𝑇 𝑃 𝑘𝐽
𝑠̅2 − 𝑠̅1 = 𝑐̅𝑝,𝑝𝑟𝑜 𝐼𝑛 𝑇2 − 𝑅𝑢 𝐼𝑛 𝑃2 (𝑘𝑚𝑜𝑙∙𝐾)
1 1

6
Calores específicos variables (análisis exacto)
Durante un proceso, cuando el cambio de temperatura es grande, no podemos suponer
calores específicos constantes, ya que nos llevaría a errores considerables al calcular el
cambio de entropía. En estos casos, debemos utilizar las relaciones exactas para estos
calores como una función de temperatura.

Para esto, se elige el cero absoluto como la temperatura de referencia y se define una
función s° como:
𝑇 𝑑𝑇 2 𝑑𝑇
𝑠 0 = ∫0 𝑐𝑝 (𝑇) 𝑠20 − 𝑠10 = ∫1 𝑐𝑝 (𝑇)
𝑇 𝑇

𝑠20 es el valor de s° en T2 y 𝑠10 es el valor en T1.


𝑃2 𝑘𝐽
𝑠2 − 𝑠1 = 𝑠20 − 𝑠10 − 𝑅 𝐼𝑛 ( )
𝑃1 𝑘𝑔 ∙ 𝐾
0 0 𝑃2 𝑘𝐽
𝑠̅2 − 𝑠̅1 = 𝑠̅ 2 − 𝑠̅ 1 − 𝑅𝑢 𝐼𝑛 ( )
𝑃 1 𝑘𝑚𝑜𝑙 ∙ 𝐾

PROCESO ISOENTRÓPICO DE GASES IDEALES


Es esta sección se van a obtener varias relaciones para los procesos isentrópicos de gases
ideales al igualar a cero las relaciones del cambio de entropía desarrolladas previamente.
 Calores específicos constantes (análisis aproximado)

R= Cp-CV; K=CP/CV en consecuencia, R/Cv = k-1


𝑅
𝑇 𝑅 𝑣 𝑇 𝑣
𝐼𝑛 𝑇2 = − 𝑐 𝐼𝑛 𝑣2 𝐼𝑛 𝑇2 = ( 𝐼𝑛 𝑣1 ) 𝑐𝑣
1 𝑣 1 1 2

𝑅
𝑇2 𝑣1 𝑐𝑣
( ) =( )
𝑇1 𝑠=𝑐𝑜𝑛𝑠𝑡𝑎𝑛𝑡𝑒 𝑣2
𝑇 𝑅 𝑃
𝐼𝑛 𝑇2 = 𝑐 𝐼𝑛 𝑃2
1 𝑃 1

(𝑘−1)
𝑐𝑃 −𝑐𝑣 𝑅 𝑘
𝑇2 𝑃2 𝑐𝑃 𝑃2 𝑘𝑐𝑣 𝑃2
𝐼𝑛 = 𝐼𝑛 ( ) = 𝐼𝑛 ( ) = 𝐼𝑛 ( )
𝑇1 𝑃1 𝑃1 𝑃1
𝑘−1
𝑇 𝑃 𝐾
(𝑇2 ) = (𝑃2 ) (Gas ideal)
1 𝑠,𝑐𝑜𝑛𝑠𝑡𝑎𝑛𝑡𝑒 1

 Calores específicos variables (tratamiento exacto)


𝑃 𝑃
0 = 𝑠20 − 𝑠10 − 𝑅 𝐼𝑛 𝑃2 𝑠20 = 𝑠10 + 𝑅 𝐼𝑛 𝑃2
1 1

7
Donde 𝑠20 es el valor de s° al final del proceso isoentrópico.

Presión relativa y volumen específico relativo


La ecuación anterior proporciona una manera exacta de evaluar cambios de propiedades de
gases ideales durante los procesos isentrópicos, ya que toma en cuenta la variación de
calores específicos debido a la temperatura. Sin embargo, involucra iteraciones tediosas
cuando la razón de volúmenes se da en lugar de la presiones. Ésta es una molestia real en
los estudios de optimización porque normalmente requiere numerosos cálculos repetitivos.
Para remediar esta deficiencia, se definen dos nuevas cantidades adimensionales asociadas
con los procesos isentrópicos, estas son:
𝑃2 𝑠20 −𝑠10 𝑃2 𝑒𝑥𝑝(𝑠0 /𝑅) 𝑠20 −𝑠10 𝑃2 𝑠20 −𝑠10
= 𝑒𝑥𝑝 = 𝑒𝑥𝑝(𝑠20 /𝑅) 𝑒𝑥𝑝 = 𝑒𝑥𝑝
𝑃1 𝑅 𝑃1 1 𝑅 𝑃1 𝑅

La cantidad exp(s°/R) se define como la presión relativa Pr, definición con la que la última
relación se vuelve:
𝑃2 𝑃𝑟2
( ) =
𝑃1 𝑠=𝑐𝑜𝑛𝑠𝑡𝑎𝑛𝑡𝑒 𝑃𝑟1

𝑃1 𝑣1 𝑃2 𝑣2 𝑣2 𝑇2 𝑃1 𝑇2 𝑃𝑟1 𝑇2 /𝑃𝑟2
= → = = =
𝑇1 𝑇2 𝑣1 𝑇1 𝑃2 𝑇1 𝑃𝑟2 𝑇1 /𝑃𝑟1

𝑣2 𝑣𝑟2
( ) =
𝑣1 𝑠=𝑐𝑜𝑛𝑠𝑡𝑎𝑛𝑡𝑒 𝑣𝑟1

8
PROBLEMA RESUELTO
2) Un comprensor de aire que maneja 300 m3/ min aumenta la presión desde 0.1 hasta 0.23
MP, y se extrae calor con una rapidez de 1700 kJ/min. La temperatura de entrada y el área
son 17°C y 280 cm2, respectivamente, mientras que en la salida los valores
correspondientes son 137 °C y 200 cm2. Determine:
1) Tasas de flujo másico (𝒎̇), en kg/s.
2) Las velocidades de entrada y salida, en m/s.
3) La potencia necesaria, en kW.

Solución:
1𝑚 2
Datos: T1=17+273=290 K P1= 0.1 MPa 𝐴𝑒𝑛𝑡 = 280 𝑐𝑚2 (100 𝑐𝑚) = 0.028 𝑚2
1𝑚 2
T2=137+273=410 K P2= 0.23 MPa 𝐴𝑠𝑎𝑙 = 200 𝑐𝑚2 (100 𝑐𝑚) = 0.020 𝑚2

𝐾𝐽 1 𝐾𝐽
𝑄̇𝑒𝑥𝑡𝑟𝑎í𝑑𝑜 = 1700 ( ) = 28.33
𝑚𝑖𝑛 60 𝑠 𝑠
𝑚3 1 𝑚𝑖𝑛 𝑚3
𝑉̇ = 𝑄 = (300 )( )=5
𝑚𝑖𝑛 60 𝑠 𝑠

Respuesta 1) Tasas de flujo másico (𝒎̇), en kg/s


PV = mRT
1000 𝐾𝑃𝑎 𝑚3
𝑃1 𝑄̇ (0.1 𝑀𝑃𝑎)( )(5 ) 𝐾𝑔
1 𝑀𝑃𝑎 𝑠
𝑚̇ = = 𝐾𝑃𝑎∙𝑚3 𝐾𝑔 𝑚𝑜𝑙
= 6.01
𝑅𝑇1 8.314( )( )(290.15)𝐾 𝑠
𝐾𝑔 𝑚𝑜𝑙∙𝐾 29 𝐾𝑔

Respuesta 2) Las velocidades de entrada y salida 𝑉̇ = 𝑉1 𝐴1

𝑚3
5 𝑠
𝑉1 = = 178.57 𝑚/𝑠
0.028 𝑚2

𝑉̇ = 𝑄 = 𝑉𝐴

𝐾𝑔 𝐾𝑃𝑎 ∙ 𝑚3 𝐾𝑚𝑜𝑙
𝑚̇𝑅𝑇1 (6.01 ) 8.314 ( )( ) (290.15 K)𝐾 𝑚
𝑠 𝐾𝑔 𝑚𝑜𝑙 ∙ 𝐾 29 𝐾𝑔
𝑉1 = = = 178.54
𝐴1 𝑃1 0.028𝑚 (100 𝐾𝑃𝑎)
2 𝑠

𝐾𝑔 𝐾𝑃𝑎 ∙ 𝑚3 𝐾𝑚𝑜𝑙
𝑚̇𝑅𝑇2 (6.01 ) 8.314 ( )( ) (410.15)𝐾 𝑚
𝑠 𝐾𝑔 𝑚𝑜𝑙 ∙ 𝐾 29 𝐾𝑔
𝑉2 = = = 153.63
𝐴2 𝑃2 0.02𝑚 (230 𝐾𝑃𝑎)
2 𝑠

9
Los valores de entalpia se encuentran en la tabla A-17 Propiedades de gas ideal del aire:
página 934.

T (K) h (KJ/kg) Pr u (KJ/kg) Vr s°


(KJ/kg∙K)
290 290.16 1.2311 206.91 676.1 1.66802
410 411.12 4.153 293.43 283.3 2.01699

T1=290.15 K h1=290.16 KJ/kg T2=410.15 K h2=411.12 KJ/kg


𝑉12 𝑉22
ℎ1 + + 𝑔𝑧1 + 𝑤𝑐𝑜𝑚𝑝𝑟𝑒𝑠𝑜𝑟 = ℎ2 + + 𝑔𝑧2 + 𝑞
2 2
𝑉22 − 𝑉12
𝑤̇𝑐𝑜𝑚𝑝 = 𝑚̇(ℎ2 − ℎ1 ) + 𝑚̇ ( ) + 𝑄̇
2

𝐾𝐽
𝐾𝑔 𝑘𝐽 153.632 − 178.542 𝑚2 ⁄𝐾𝑔 ̇ 𝐾𝐽
𝑤𝑐𝑜𝑚𝑝 = 6.01 ((411.12 − 290.16) +( )( 2 )( 2 )) + (28.33 )
𝑠 𝑘𝑔 2 𝑠 1000 𝑚 ⁄ 2 𝑠
𝑠

𝑤̇𝑐𝑜𝑚𝑝 = 6.01 (120.96 − 4.1372)+(28.33 ) = 702.1 + ̇ 28.33 = 730.43 KW

10

También podría gustarte

pFad - Phonifier reborn

Pfad - The Proxy pFad of © 2024 Garber Painting. All rights reserved.

Note: This service is not intended for secure transactions such as banking, social media, email, or purchasing. Use at your own risk. We assume no liability whatsoever for broken pages.


Alternative Proxies:

Alternative Proxy

pFad Proxy

pFad v3 Proxy

pFad v4 Proxy