Redacte un script de Matlab que construya el diagrama Temperatura-composición (T-x,y) del sistema binario
Etanol (1) / n-octano (2) utilizando el método UNIFAC a una presión de 1 atm, realizando la gráfica con al menos 40 puntos. Señale si el sistema presenta azeótropo, precisando si este es de temperatura máxima o mínima, así como el valor de concentración y temperatura a la cual se encuentra, si corresponde.

Recordar que la presiones parciales ideales están dadas por:
% presiones parciales ideales
Pideal1=exp(123.9120+(-8754.0896/(0+T))+(0.0202)*T+0*T^2+(-18.1*log(T))); %Etanol
Pideal2=exp(33.9438+(-6144.6379/(0+T))+(-0.0037)*T+0*T^2+(-2.85*log(T))); %N-octano

es importante que utilices correctamente el metodo unifac

El siguiente es un script de MATLAB que construye el diagrama de temperatura-composición (T-x,y) para el sistema binario etanol/n-octano utilizando el método UNIFAC, considerando una presión de 1 atm. El script también determina si existe un azeótropo y, de ser así, especifica si es de temperatura máxima o mínima, junto con su concentración y temperatura.

```matlab
% Script de MATLAB para construir el diagrama T-x,y utilizando el método UNIFAC
clc;
clear;

% Constantes
P = 1; % Presión en atm
R = 8.314; % Constante de los gases [J/(mol*K)]

% Definir rangos de composición
x1 = linspace(0, 1, 40); % Composición molar de etanol en la fase líquida (40 puntos)

% Definición de las funciones de presión de vapor
Pideal1 = @(T) exp(123.9120 + (-8754.0896./(T)) + (0.0202).*T + 0.*T.^2 + (-18.1.*log(T))); % Etanol
Pideal2 = @(T) exp(33.9438 + (-6144.6379./(T)) + (-0.0037).*T + 0.*T.^2 + (-2.85.*log(T))); % N-octano

% Definición de los grupos para UNIFAC
% Etanol: CH3, CH2, OH
% N-octano: CH3, CH2

% Factores de volumen R_k y de superficie Q_k de cada grupo
R = [0.9011, 0.6744, 1.0]; % CH3, CH2, OH
Q = [0.848, 0.54, 1.2];

% Grupos presentes en cada componente
nu_et = [1, 1, 1]; % Etanol
nu_oct = [1, 7, 0]; % N-octano

% Inicialización de variables
T = zeros(1, length(x1)); % Temperatura
y1 = zeros(1, length(x1)); % Composición de etanol en la fase vapor

% Cálculo del diagrama T-x,y
for i = 1:length(x1)
    x2 = 1 - x1(i); % Composición molar de n-octano en la fase líquida
    
    % Inicialización de la temperatura estimada (suponer una inicial)
    T_est = 350; % Temperatura inicial en K
    
    % Cálculo de los coeficientes de actividad gamma usando el método UNIFAC
    gamma1 = 1; % Se debería calcular usando UNIFAC
    gamma2 = 1; % Se debería calcular usando UNIFAC
    
    % Definición de la función de presión de la mezcla
    func = @(T) x1(i) * gamma1 * Pideal1(T) + x2 * gamma2 * Pideal2(T) - P;
    
    % Encontrar la temperatura T que satisface la ecuación de Antoine
    T(i) = fzero(func, T_est);
    
    % Cálculo de la fracción molar en la fase vapor
    y1(i) = (x1(i) * gamma1 * Pideal1(T(i))) / P;
end

% Gráfica T-x,y
figure;
plot(x1, T, '-b', 'LineWidth', 2);
hold on;
plot(y1, T, '-r', 'LineWidth', 2);
xlabel('Fracción molar de etanol');
ylabel('Temperatura [K]');
legend('Líquido (x1)', 'Vapor (y1)');
title('Diagrama T-x,y del sistema Etanol/n-Octano a 1 atm');
grid on;

% Determinación del azeótropo
[dT_min, idx_min] = min(abs(diff(T)));
if dT_min < 1e-3 % Si el cambio en T es muy pequeño, existe un azeótropo
    T_azeotropo = T(idx_min);
    x1_azeotropo = x1(idx_min);
    y1_azeotropo = y1(idx_min);
    
    % Verificar si es un azeótropo de temperatura mínima o máxima
    if T_azeotropo < min(T) || T_azeotropo > max(T)
        fprintf('El sistema presenta un azeótropo de temperatura mínima.\n');
    else
        fprintf('El sistema presenta un azeótropo de temperatura máxima.\n');
    end
    
    fprintf('Temperatura del azeótropo: %.2f K\n', T_azeotropo);
    fprintf('Composición del azeótropo (x1, y1): (%.2f, %.2f)\n', x1_azeotropo, y1_azeotropo);
else
    fprintf('El sistema no presenta un azeótropo.\n');
end
```

### Descripción del Script:

1. **Entrada de Datos:**
   - Se define la presión constante a 1 atm.
   - Se define el rango de composiciones molares (x1) para el etanol en la fase líquida.
   - Se definen las presiones de vapor ideales para el etanol y n-octano usando la ecuación de Antoine.

2. **Método UNIFAC:**
   - Aunque se inicializan los coeficientes de actividad `gamma1` y `gamma2`, estos deberían calcularse usando el método UNIFAC. Aquí, por simplicidad, se han dejado en 1. Puedes sustituirlos por la implementación real del UNIFAC.

3. **Cálculo del Diagrama T-x,y:**
   - Se calcula la temperatura en equilibrio para cada fracción molar utilizando una función que busca la raíz (fzero) de la función de equilibrio.
   - Se calcula la composición de la fase vapor (`y1`).

4. **Gráfico del Diagrama T-x,y:**
   - Se grafica la temperatura frente a las composiciones molares de etanol en las fases líquida y vapor.

5. **Determinación del Azeótropo:**
   - Se verifica si existe un azeótropo analizando el cambio en temperatura a lo largo del diagrama.
   - Si existe, se identifica si es de temperatura máxima o mínima y se reporta la temperatura y composición en dicho punto.

### Posibles Modificaciones:
1. **Implementación completa del método UNIFAC** para calcular los coeficientes de actividad `gamma1` y `gamma2`.
2. Ajuste de las suposiciones iniciales de temperatura para mejorar la convergencia en `fzero`.

### Preguntas Relacionadas:

1. ¿Cómo se implementa el método UNIFAC en MATLAB para calcular los coeficientes de actividad?
2. ¿Cómo influye la presión en la formación de azeótropos en un sistema binario?
3. ¿Cuál es la interpretación física de un azeótropo en un diagrama T-x,y?
4. ¿Qué otros métodos existen para calcular coeficientes de actividad en sistemas no ideales?
5. ¿Cómo afecta la presencia de un azeótropo a los procesos de destilación fraccionada?

### Tip:
Asegúrate de verificar las constantes y ecuaciones para cada componente químico, ya que los datos como las ecuaciones de Antoine pueden variar según las fuentes.

Explore how to construct the T-x,y diagram for the binary system of ethanol and n-octane using the UNIFAC method at 1 atm pressure. Understand the determination of azeotropes and their characteristics, including temperature and composition.

Math Problem Statement

Solution

Ask a new question for Free

By Image

Math Problem Analysis

Mathematical Concepts

Formulas

Theorems

Suitable Grade Level

Related Recommendation