Resolver numéricamente el siguiente problema, además realizar DF y programa en MATLAB.
1.- La entalpia requerida para calentar n moles de un gas desde T 1 a T2, se encuentra integrando
la capacidad calorífica a presión constante sobre el rango de temperaturas según la siguiente
ecuación:
𝑇2
∆𝐻 = 𝑛 ∫ 𝐶𝑝 (𝑇)𝑑𝑇
𝑇1
En la tabla adjunta se muestra datos de capacidad calorífica par el CO 2 a diferentes
temperaturas.
Calcular la entalpia de un mol de CO2 en el rango de temperaturas entre T 1 = 100 y T2 =
800 °C, use la fórmula de Simpson 1/3.
T, °C Cp (J/mol K)
100 40.461
150 42.256
200 43.881
250 45.355
300 46.695
350 47.917
400 49.034
450 50.055
500 50.989
550 51.843
600 52.624
650 53.339
700 53.994
750 54.593
800 55.144
2.-
(16.3 CHAPRA) La razón de crecimiento de una levadura que produce un antibiótico es una
función de la concentración del alimento c,
Como se ilustra en la figura P16.3, el crecimiento parte de cero a muy bajas concentraciones
debido a la limitación de la comida. También parte de cero en altas concentraciones debido a los
efectos de toxicidad. Encuentre el valor de c para el cual el crecimiento es un máximo.
3.-
(16.6 CHAPRA) Hay que separar una mezcla de benceno y tolueno en un reactor flash. ¿A qué
temperatura deberá operarse el reactor para obtener la mayor pureza de tolueno en la fase líquida
(maximizar xT)? La presión en el reactor es de 800 mm Hg. Las unidades en la ecuación de
Antoine son mm Hg y °C para presión y temperatura, respectivamente.
4.- AJUSTE DE CURVAS POR MINIMOS CUADRADOS – REGRESION SIMPLE
Como resultado de experimentos científicos y/o ingenieriles, se obtienen por mediciones
cantidades físico-químicas que son registradas en un conjunto de datos. Es deseable que estos datos
se ajusten a un modelo matemático que describa lo más exactamente posible los hechos
experimentales; para este propósito inicialmente se realiza una gráfica de tendencia de los datos
experimentales obtenidos, luego por analogía con las formas geométricas de funciones estándar
(modelo empírico) o con un modelo matemático obtenido de la abstracción teórica del proceso
(modelo analítico), se elige una función matemática de ajuste a partir de la que se procede a
encontrar los parámetros (coeficientes) de la función de tal manera que a través de un análisis
estadístico se puede confirmar o negar que la función elegida se ajusta bien a los datos
experimentales. El modelo matemático de esta manera encontrado, se puede usar para incluirlo en
cálculos de modelaje y simulación on line, también para estimar propiedades físico-químicas entre
puntos conocidos por interpolación, o para predecir datos fuera del rango de medición por
extrapolación.
TABLA 3.11.- Calor Figura 3.6.- Curva
Específico del Agua Ajustada
Figura 3.4.- Curva Ajustada
T, °C cp, cal/g °C 1,008
10 1,001003 1,006
20 0,998376 1,004
Cp, cal/g °C
f(x)
30 0,997182 1,002
AJUSTE
40 0,996943 1 DATOS
50 0,997659 0,998
60 0,998854 0,996
70 1,000287 0 20 40 60 80 100 120
T, °C
80 1,001959
90 1,004108
100 1,006497
Por el método de mínimos cuadrados plantee un programa capaz de realizar el ajuste curvas de
polinomios de la forma:
𝑐𝑝 = 𝑎1 + 𝑎2 𝑇 + 𝑎3 𝑇 2 + ⋯ + 𝑎𝑛+1 𝑇 𝑛
Encuentre el mejor ajuste polinomial de grado n. Para definir cual es el mejor ajuste, debe
calcular el coeficiente de correlación para cada caso y definir el que más se aproxime a 1.
5.-
REGRESIÓN MÚLTIPLE LINEAL: AJUSTE DE CURVAS POR MINIMOS
CUADRADOS
La ecuación de Antoine permite calcular la presión de vapor (𝑃𝑉 , 𝑚𝑚 𝐻𝑔) de un componente puro en
función de la temperatura (𝑇, °𝐶) :
𝐵
(𝐴+𝑇+𝐶 )
𝑃𝑉 = 10 (𝐸𝑐. 3.91)
La ecuación de Antoine es empírica y sus constantes (𝐴, 𝐵, 𝐶) son características para cada
compuesto puro y su determinación se realiza a partir de las presiones de vapor medidas
experimentalmente.
Calcule las constantes de Antoine del propano a partir de los siguientes datos de presión de vapor
del propano mostrado en la tabla:
TABLA.- Presión de vapor de equilibrio del propano
i Ti, °C Pi, mmHg i Ti, °C Pi, mmHg
1 -70 7.37 11 30 66.3
2 -60 9.72 12 40 78.0
3 -50 12.6 13 50 91.8
4 -40 16.2 14 60 107.1
5 -30 20.3 15 70 124.0
6 -20 25.4 16 80 142.8
7 -10 31.4 17 90 164.0
8 0 38.2 18 100 187.0
9 10 46.0 19 110 213.0
10 20 55.5 20 120 240.0
Estrategia de solución:
Para encontrar las constantes de la ecuación 3.91, debemos linealizarla y utilizando regresión
lineal múltiple, determinar los parámetros A, B y C.
La ecuación de Antoine se puede re-escribir para linealizarla de la siguiente manera:
𝐵 𝐵
𝑃𝑉 = 10(𝐴+𝑇+𝐶 ) → log(𝑃𝑉 ) = 𝐴 + → T log(𝑃𝑉 ) = (𝐴𝐶 + 𝐵) + 𝐴𝑇 − 𝐶 log(𝑃𝑉 )
𝑇+𝐶
(Ec. 3.92)
Esta ecuación es lineal múltiple, donde 𝑦 es una función lineal de 𝑥1 y 𝑥2 de la forma:
𝑦 = 𝑎 + 𝑏𝑥1 + 𝑐𝑥2 (𝐸𝑐. 3.93)
Donde:
𝑦 = T log(𝑃𝑉 )
𝑥1 = 𝑇 ; 𝑥2 = log(𝑃𝑉 )
𝑎 = 𝐴𝐶 + 𝐵; 𝑏 = 𝐴 ; 𝑐 = −𝐶
La ecuación 3.93, es multivariante donde la variable dependiente (𝑦) depende en este caso de dos
variables independientes (𝑥1 , 𝑥2 ). Este caso es de dos dimensiones, por lo que la línea de
regresión se convierte en un plano.
Los mejores valores de los coeficientes (𝐴, 𝐵 𝑦 𝐶) se obtienen mediante el método de mínimos
cuadrados. Primero se formula la sumatoria de los cuadrados de los residuos y luego se minimiza
esta suma de los errores al cuadrado, entonces:
𝑛
2
𝑆 = ∑(𝑦𝑖 − 𝑎 − 𝑏𝑥1,𝑖 − 𝑐𝑥2,𝑖 )
𝑖=1
6.-
La capacidad calorífica del gas propano se puede expresar como un polinomio en función de la
temperatura:
𝑐𝑝 = 𝑎0 + 𝑎1 𝑇 + 𝑎2 𝑇 2 + 𝑎3 𝑇 3 + ⋯
Realice un programa en MATAB para calcular las constantes 𝑎0 , 𝑎1 , 𝑎2 , 𝑎3 , … por el método
de ajuste por mínimos cuadrados de los datos de la siguiente tabla a polinomios de diferentes
grados y elija el mejor ajuste de acuerdo al coeficiente de correlación más cercano a la unidad.
Represente los 𝑐𝑝 calculados y experimentales vs. T en la misma gráfica.
cp, kJ/(kmol
T, K T, K cp, kJ/(kmol K)
K)
50 34.06 700 142.67
100 41.30 800 154.77
150 48.79 900 163.35
200 56.07 1000 174.60
273.16 68.74 1100 182.67
298.15 73.60 1200 189.74
300 73.93 1300 195.85
400 94.01 1400 201.21
500 112.59 1500 205.89
600 128.70
7.-
La potencia generada por un molino de viento varía con la velocidad del viento. En un
experimento, se obtuvieron las cinco mediciones siguientes:
velocidad,
km/h 22.53 35.41 48.28 61.15 74.03
Potencia, W 320 490 540 500 480
a) Determine el polinomio de Lagrange de 4to. orden que pase por los cinco puntos. Usando este
polinomio calcular la potencia a una velocidad de viento de 42 km/hr.
b) Ajustando por mínimos cuadrados estos datos a un polinomio de 4to.grado (calcule el coeficiente
de correlación), que diferencias encuentra con los resultados del anterior inciso.
c) Cual la diferencia de uso entre el método de Lagrange y el de mínimos cuadrados?.
8.-
9.-
.
10.-
24.50 La entalpía de un gas real es función de la presión como se describe a continuación. Los
datos se tomaron para un fluido real. Estime la entalpía del fluido a 400 K y 50 atm (evalúe la
integral de 0.1 atm a 50 atm).
