Conocimiento Educación en Ingeniería Química ¿Por qué es fundamental elegir correctamente el modelo de propiedades termodinámicas? Dimensione columnas de destilación no ideales con confianza
Avatar del autor

Equipo técnico · LABPARK

Actualizado hace 2 meses

¿Por qué es fundamental elegir correctamente el modelo de propiedades termodinámicas? Dimensione columnas de destilación no ideales con confianza


El modelo de propiedades termodinámicas es fundamental porque determina los datos de equilibrio de los que depende cada decisión importante del diseño de destilación. Para un sistema binario no ideal, el modelo debe predecir con precisión el equilibrio vapor-líquido (VLE), incluidos los coeficientes de actividad, la volatilidad relativa y los posibles azeótropos. Un modelo inadecuado puede producir requisitos incorrectos de etapas, relaciones de reflujo, presiones de operación, perfiles de temperatura, cargas térmicas e incluso la falsa impresión de que una separación es factible.

El diseño de una destilación solo es tan fiable como sus predicciones del equilibrio de fases. Para sistemas binarios no ideales, los modelos de coeficientes de actividad, como NRTL o UNIQUAC, respaldados por parámetros de interacción binaria verificados, suelen ser más adecuados que la suposición de fluido ideal o una ecuación de estado no validada.

Por qué las predicciones de equilibrio controlan el diseño de la columna

El VLE define las composiciones de vapor y líquido

Una columna de destilación separa los componentes porque las fases vapor y líquido tienen composiciones diferentes en el equilibrio. El modelo de propiedades calcula esta relación, comúnmente representada mediante la relación de equilibrio:

[ K_i = \frac{y_i}{x_i} ]

donde (x_i) es la composición de la fase líquida y (y_i) es la composición de la fase vapor.

Para líquidos no ideales, (K_i) depende en gran medida del coeficiente de actividad (\gamma_i), la presión de vapor, la presión y la fugacidad de la fase vapor. Un error en (\gamma_i) modifica la composición de vapor predicha y, por tanto, cambia el rendimiento de separación calculado.

La volatilidad relativa determina la dificultad de separación

La volatilidad relativa indica con qué facilidad se separa un componente de otro. Cuando la volatilidad relativa predicha es demasiado alta, la columna puede parecer más fácil de operar de lo que realmente es.

Ese error puede dar lugar a columnas subdimensionadas, reflujo insuficiente, pureza inadecuada del producto y baja recuperación durante la operación piloto.

Los azeótropos establecen los límites físicos de la separación

Un azeótropo se produce cuando las composiciones del vapor y del líquido se vuelven iguales a una composición y presión determinadas. En ese punto, la destilación de equilibrio convencional no puede seguir enriqueciendo un componente a través de dicha composición.

Un modelo de fluido ideal no puede predecir azeótropos porque supone que los coeficientes de actividad son iguales a uno. Para un sistema binario azeotrópico o fuertemente no ideal, esa suposición elimina precisamente el comportamiento que determina si la separación deseada es posible.

Por qué los modelos ideales fallan para sistemas binarios no ideales

El comportamiento ideal elimina las interacciones moleculares

Un modelo de líquido ideal supone que cada componente experimenta esencialmente el mismo entorno molecular en la mezcla que en el líquido puro. Sin embargo, las mezclas reales pueden presentar interacciones fuertes o débiles entre moléculas diferentes, enlaces de hidrógeno, efectos de polaridad o asociación.

Estas interacciones modifican la tendencia de escape de cada componente. Las desviaciones resultantes respecto a la idealidad pueden ser lo bastante grandes como para crear azeótropos o invertir el orden de volatilidad esperado.

Las curvas de equilibrio incorrectas distorsionan los cálculos de etapas

El análisis de McCabe-Thiele, los cálculos rigurosos etapa por etapa y los métodos abreviados dependen de la curva de equilibrio. Si la curva es inexacta, también lo serán el número estimado de etapas teóricas y el requisito de reflujo.

El problema no se limita a una pequeña corrección numérica. Una curva incorrecta puede cambiar el reflujo mínimo predicho, las ubicaciones de los puntos de pinzamiento, el intervalo operativo factible y la configuración de la columna.

La factibilidad falsa es más peligrosa que un fallo evidente

Un modelo ideal puede predecir una separación limpia de productos incluso cuando un límite azeotrópico impide alcanzar ese resultado. Esto puede llevar a los ingenieros a proceder con un diseño de destilación convencional que no puede alcanzar la pureza especificada sin métodos suplementarios de separación por cambio de presión, extracción, destilación azeotrópica, membranas u otros.

Por qué suelen preferirse los modelos de coeficientes de actividad

NRTL representa la no aleatoriedad molecular local

El modelo Non-Random Two-Liquid (NRTL) tiene en cuenta las diferencias en las interacciones moleculares locales dentro de la fase líquida. Se utiliza ampliamente para mezclas muy no ideales y polares cuando se dispone de parámetros de interacción binaria fiables.

Para el diseño de columnas, ajustar los parámetros NRTL a datos experimentales de VLE suele ser más preciso que depender de estimaciones genéricas.

UNIQUAC captura el comportamiento del líquido dependiente de la composición

UNIQUAC utiliza información sobre el tamaño molecular y el área superficial junto con parámetros de interacción. Puede representar el comportamiento no ideal del líquido en distintos intervalos de composición y se aplica habitualmente a sistemas polares, asociantes y azeotrópicos.

Su utilidad depende de la calidad de sus parámetros y de la adecuación de la representación subyacente al sistema químico.

UNIFAC ayuda cuando los datos son incompletos

UNIFAC estima los coeficientes de actividad a partir de contribuciones de grupos funcionales. Es útil cuando los datos experimentales de VLE son limitados o cuando deben representarse numerosos pares binarios en un sistema multicomponente.

Sin embargo, UNIFAC es un método predictivo, no un sustituto de la validación. Su precisión depende de la disponibilidad y aplicabilidad de los parámetros de grupo relevantes, especialmente cerca de los azeótropos y en sistemas fuertemente asociantes.

Las fases vapor y líquida pueden necesitar modelos diferentes

Para muchas mezclas líquidas polares, un modelo de coeficientes de actividad describe la fase líquida, mientras que una ecuación de estado estima la fugacidad de la fase vapor. Este enfoque combinado puede ser más adecuado que imponer un único modelo para representar ambas fases.

A presiones moderadas, la no idealidad de la fase vapor puede ser menor en algunos sistemas. A presiones más altas o para vapores condensables y fuertemente no ideales, las correcciones de fugacidad adquieren mayor importancia.

Cómo la elección del modelo modifica el dimensionamiento de la columna

El número de etapas depende de la relación de equilibrio

El número mínimo de etapas y el número real de etapas teóricas dependen del factor de separación predicho por el modelo de propiedades. Un modelo que predice una separación más intensa produce menos etapas calculadas, mientras que uno que predice una separación más débil produce más.

Por ejemplo, ignorar los parámetros de interacción binaria puede cambiar de forma significativa el número calculado de etapas y los requisitos de reflujo. Incluso un cambio aparentemente moderado puede afectar a la altura de la columna, la selección de bandejas o empaques y el coste del proyecto.

La relación de reflujo cambia la demanda energética

La relación de reflujo determina cuánto líquido condensado de cabeza se devuelve a la columna. Un mayor reflujo generalmente mejora la separación, pero aumenta las cargas térmicas del condensador y del rehervidor.

Si el modelo termodinámico subestima la dificultad de separación, puede producir una relación de reflujo irrealmente baja. La columna resultante puede presentar una pureza insuficiente y un consumo real de servicios mayor que el previsto.

La selección de la presión depende del comportamiento de las fases

La presión de operación afecta a las temperaturas de ebullición, las densidades del vapor, la volatilidad relativa y la ubicación de los límites azeotrópicos. Por tanto, se necesita un modelo de propiedades para evaluar las combinaciones de presión y determinar si los cambios de presión mejoran o limitan la separación.

La selección de la presión también afecta a la temperatura del condensador, la temperatura del rehervidor, la compatibilidad con los servicios y los requisitos mecánicos del equipo.

El diseño hidráulico utiliza resultados termodinámicos

Las predicciones termodinámicas determinan los caudales de vapor y líquido a través de la columna. Estos caudales influyen en la carga de las bandejas, los límites de inundación, el riesgo de goteo, la capacidad del empaque, el diámetro y la pérdida de presión.

Por tanto, un error termodinámico puede convertirse en un error de diseño hidráulico, incluso cuando los cálculos hidráulicos se realizan correctamente.

Cómo seleccionar y validar el modelo

Comience por el comportamiento químico

El modelo debe reflejar el comportamiento físico dominante de la mezcla.

  • Sistemas de hidrocarburos a presión elevada: SRK o Peng-Robinson pueden proporcionar predicciones fiables del comportamiento de las fases.
  • Mezclas líquidas polares o asociantes: NRTL, UNIQUAC o Wilson suelen ser más adecuados para la fase líquida.
  • Sistemas azeotrópicos: Utilice un modelo de coeficientes de actividad con parámetros verificados y confirme que reproduce la composición y la temperatura azeotrópicas.
  • Sistemas con datos incompletos: Utilice UNIFAC u otro método predictivo con cautela, seguido de una validación cuando sea posible.
  • Sistemas con separación líquido-líquido: Verifique que el modelo pueda representar el equilibrio líquido-líquido, no solo el VLE.

Ningún modelo es universalmente mejor. La idoneidad depende de los componentes, el intervalo de presiones, el intervalo de temperaturas y el comportamiento de equilibrio que se esté modelando.

Ajuste los parámetros frente a datos experimentales

Cuando se dispone de datos binarios de VLE, los parámetros de interacción deben ajustarse en el intervalo de composición y de temperatura o presión relevante para la columna.

Un conjunto de parámetros que se ajusta a un intervalo estrecho de datos puede funcionar mal en otras condiciones. Por ello, el conjunto de validación debe incluir la región de operación y, en sistemas azeotrópicos, el área cercana al límite azeotrópico.

Compruebe la curva de equilibrio completa

No valide un modelo utilizando únicamente un punto de ebullición o un dato de equilibrio. Compare las composiciones de equilibrio predichas y medidas en todo el intervalo de composición.

Las comprobaciones importantes incluyen:

  • Temperaturas de burbuja y de rocío.
  • Composición del vapor en todo el intervalo binario.
  • Tendencias de la volatilidad relativa.
  • Ubicación y tipo de azeótropo.
  • Separaciones de fases líquido-líquido cuando corresponda.
  • Sensibilidad a la presión.

Realice estudios de sensibilidad

Deben compararse varios modelos termodinámicos plausibles cuando los datos sean inciertos. El diseño debe probarse frente a cambios en el número de etapas, el reflujo, la presión, la carga térmica, la recuperación de producto y la carga hidráulica.

Esto permite identificar si el modelo de propiedades representa una incertidumbre menor o el riesgo dominante del diseño.

Comprensión de las ventajas y desventajas

Los modelos predictivos son convenientes, pero menos seguros

UNIFAC y otros métodos predictivos son valiosos para el diseño preliminar y para sistemas con pocos datos. No obstante, sus predicciones pueden ser poco fiables para grupos funcionales inusuales, asociaciones fuertes, electrolitos o regiones cercanas a los límites de fase.

Deben respaldar el criterio de ingeniería y la recopilación de datos específicos, en lugar de considerarse automáticamente validados.

Los modelos más complejos no son automáticamente mejores

Un modelo sofisticado con parámetros mal ajustados puede ser menos fiable que un modelo más sencillo respaldado por buenos datos experimentales. La complejidad no compensa unos parámetros de interacción incorrectos ni un marco inadecuado de equilibrio de fases.

El modelo adecuado es aquel que reproduce el comportamiento físico relevante con evidencias defendibles.

Los modelos EOS tienen un papel válido, pero limitado

SRK y Peng-Robinson son eficaces para muchos sistemas de hidrocarburos y determinados intervalos de presión en los que el comportamiento de la ecuación de estado está bien representado. Los modelos EOS estándar pueden producir curvas de equilibrio poco fiables para sistemas líquidos polares muy no ideales si no representan adecuadamente las interacciones de la fase líquida.

Esto no constituye un motivo para rechazar universalmente los métodos EOS. Es un motivo para ajustar el modelo y sus reglas de mezcla a la química y a las condiciones de operación.

Una planta piloto no puede corregir un modelo fundamentalmente incorrecto

Los datos de planta piloto pueden revelar discrepancias, pero no pueden hacer físicamente correcto un modelo termodinámico inadecuado. Si la simulación predice que no existe un azeótropo cuando sí existe, o predice una separación de fases incorrecta, ajustar los factores de eficiencia o los parámetros hidráulicos solo puede ocultar el error subyacente.

La validación termodinámica debe realizarse antes de utilizar los resultados de la planta piloto para ajustar el equipo o los parámetros de operación.

Cómo tomar la decisión adecuada para su objetivo

Elija el método de propiedades de acuerdo con el riesgo dominante del diseño y valídelo frente al comportamiento de equilibrio que controla la columna.

  • Si su prioridad es obtener una pureza de producto precisa: Utilice NRTL o UNIQUAC con parámetros de interacción binaria verificados y confirme la curva de VLE predicha y el comportamiento azeotrópico frente a datos experimentales.
  • Si su prioridad es realizar una evaluación preliminar con datos limitados: Utilice UNIFAC u otro método predictivo de coeficientes de actividad, pero considere provisionales los resultados de número de etapas, reflujo y factibilidad.
  • Si su prioridad es el procesamiento de hidrocarburos a presión elevada: Evalúe SRK o Peng-Robinson y verifique sus predicciones en el intervalo real de presión y temperatura.
  • Si su prioridad es el escalado de una planta piloto: Compare los modelos termodinámicos mediante estudios de sensibilidad que incluyan etapas, reflujo, cargas térmicas, temperaturas, caudales y límites hidráulicos.
  • Si su prioridad es un sistema con separación de fases líquido-líquido: Seleccione y valide un modelo capaz de representar el LLE, incluidos los datos de distribución o los parámetros de interacción respaldados por mediciones experimentales.

Un modelo termodinámico defendible convierte el dimensionamiento de una columna de un ejercicio matemático en una predicción creíble de cómo se comportará la separación real.

Tabla resumen:

Factor de diseño Efecto del modelo de propiedades Necesidad principal de validación
VLE y coeficientes de actividad Determina las composiciones de equilibrio y la precisión de la separación Comparar los datos de VLE predichos y medidos
Volatilidad relativa Controla el número de etapas y la dificultad de separación Verificar las tendencias de volatilidad en todo el intervalo de composición
Azeótropos Define si la destilación convencional puede alcanzar la pureza objetivo Confirmar la ubicación y el tipo de azeótropo
Reflujo y carga térmica Influye en los requisitos energéticos del condensador y del rehervidor Realizar estudios de sensibilidad del modelo
Presión de operación Afecta a las temperaturas, la volatilidad, el comportamiento de las fases y los servicios Validar las predicciones en todo el intervalo de presión
Diseño hidráulico Determina los caudales de vapor y líquido, el diámetro, la inundación y la pérdida de presión Comprobar conjuntamente los resultados termodinámicos e hidráulicos
Selección del modelo NRTL o UNIQUAC suelen ser adecuados para líquidos polares no ideales; SRK o Peng-Robinson son apropiados para muchos hidrocarburos a presión elevada Adaptar el modelo a la química y a las condiciones de operación

¿Necesita un diseño de destilación fiable para un sistema binario no ideal? Contacte con LABPARK para analizar su aplicación. LABPARK proporciona plantas piloto de operaciones unitarias con fines educativos y de formación profesional para ingeniería química, bioprocesos y biotecnología, así como para el tratamiento medioambiental y de aguas. Nuestros sistemas de destilación a escala piloto ayudan a universidades, institutos de investigación y empresas a validar el comportamiento del VLE, comparar modelos termodinámicos, evaluar etapas y reflujo, medir las cargas energéticas y aumentar la confianza antes del escalado. Trabaje con LABPARK para seleccionar equipos que respalden la enseñanza práctica, la investigación, el desarrollo de procesos y la formación industrial.

Productos relacionados

La gente también pregunta

Productos relacionados

Planta Piloto de Destilación Continua de Platos de Cribado para la Educación de Laboratorio de Operaciones Unitarias

Planta Piloto de Destilación Continua de Platos de Cribado para la Educación de Laboratorio de Operaciones Unitarias

Sistema de enseñanza integrado a escala piloto para estudios de destilación continua de platos de cribado. Demostración visual de la hidráulica de los platos, posiciones de alimentación flexibles y control automático de reflujo para la educación práctica de operaciones unitarias en laboratorios de ingeniería. Diseñado para laboratorios de ingeniería de educación superior.

Planta Piloto Educativa de Destilación Especial Multifuncional

Planta Piloto Educativa de Destilación Especial Multifuncional

Planta piloto de destilación especial multifuncional versátil para la educación en ingeniería química. Admite destilación continua, al vacío, azeotrópica, reactiva y extractiva. Las columnas de vidrio transparente permiten la observación visual en tiempo real de la hidrodinámica y los procesos de separación.

Planta Piloto de Entrenamiento en Operaciones de Unidad de Destilación Multimodal

Planta Piloto de Entrenamiento en Operaciones de Unidad de Destilación Multimodal

Planta piloto de destilación multimodal para entrenamiento práctico en operaciones de unidad para la enseñanza de ingeniería química. Cuenta con modos de simulación real, analógica y semifísica, construcción industrial y es personalizable para laboratorios universitarios. Incluye columnas de fraccionamiento prácticas, control SCADA y sistemas de seguridad. Incorpora mirillas de inspección y puertos de muestreo con reciclaje en circuito cerrado.

Planta Piloto Educativa de Destilación Extractiva por Lotes Continua

Planta Piloto Educativa de Destilación Extractiva por Lotes Continua

Planta piloto versátil para la formación en destilación continua, por lotes y extractiva. Columna de vidrio de alto borosilicato para visualizar la hidráulica, pantalla táctil de 15,6 pulgadas con registro de datos, control preciso de la relación de reflujo 1-99 y marco resistente a la corrosión. Ideal para la educación en ingeniería química y la investigación de procesos.

Planta Piloto de Rectificación de Doble Modo para la Formación Práctica en Operaciones Unitarias

Planta Piloto de Rectificación de Doble Modo para la Formación Práctica en Operaciones Unitarias

Planta piloto de rectificación de doble modo a escala industrial para la formación práctica en ingeniería química. Cuenta con modos de operación de material real y material simulado, columna de platos perforados con visores para la observación visual de la hidrodinámica y control SCADA personalizable para un aprendizaje seguro y práctico de operaciones unitarias y transferencia de materia.

Planta Piloto Educativa de Extracción Líquido-Líquido con Disco Rotatorio

Planta Piloto Educativa de Extracción Líquido-Líquido con Disco Rotatorio

Una columna transparente de discos rotatorios para experimentos educativos de extracción líquido-líquido. Esta planta piloto permite a los estudiantes estudiar la transferencia de masa, la dinámica de gotas y el comportamiento de inundación, uniendo la teoría y la práctica en la educación de operaciones unitarias de ingeniería química. Cuenta con agitación de velocidad variable y control por PLC.

Planta Piloto Integral de Extracción Líquido-Líquido para Educación en Ingeniería

Planta Piloto Integral de Extracción Líquido-Líquido para Educación en Ingeniería

Planta piloto integral de extracción líquido-líquido para educación en ingeniería, que integra columnas rotativas y vibratorias para la observación práctica del comportamiento de fases, límites de inundación y eficiencia de transferencia de masa, permitiendo cálculos precisos de Unidad de Transferencia de Altura (HTU) y coeficientes de transferencia de masa.

Planta Piloto Educativa de Purificación de Agua por Intercambio Iónico para Operaciones Unitarias de Ingeniería

Planta Piloto Educativa de Purificación de Agua por Intercambio Iónico para Operaciones Unitarias de Ingeniería

Esta planta piloto de intercambio iónico a escala de banco entrena a estudiantes de ingeniería en la purificación de agua. Dos columnas transparentes duales simulan el ablandamiento y desmineralización industriales. Los estudiantes observan la dinámica de fluidos, realizan la regeneración de resinas y analizan curvas de ruptura. El marco resistente a la corrosión garantiza durabilidad en los experimentos de operaciones unitarias.

Planta Piloto Educativa de Operaciones Unitarias de Absorción en Lecho Empaquetado

Planta Piloto Educativa de Operaciones Unitarias de Absorción en Lecho Empaquetado

Estudie la absorción gas-líquido, la caída de presión, el inundamiento y los coeficientes de transferencia de materia con esta planta piloto. Columna empacada transparente, pantalla táctil industrial, datos de sensores en tiempo real, análisis automatizado. Investigue el flujo bifásico, los puntos de carga y la eficiencia de la columna. Incluye software integral de registro de datos y evaluación.

Planta Piloto de Entrenamiento en Operaciones Unitarias de Destilación Extractiva para la Purificación de Etanol Anhidro Verde

Planta Piloto de Entrenamiento en Operaciones Unitarias de Destilación Extractiva para la Purificación de Etanol Anhidro Verde

Planta piloto modular produce etanol anhidro de alta pureza a partir de etanol crudo mediante destilación extractiva en un proceso de circuito cerrado de cero emisiones que proporciona formación práctica en operaciones unitarias con software SCADA de control basado en PLC y gestión de procesos digitalizada centrándose en principios de ingeniería verde

Planta Piloto Educativa de Destilación, Purificación y Formulación de Electrolitos

Planta Piloto Educativa de Destilación, Purificación y Formulación de Electrolitos

Planta piloto educativa integrada a escala de banco a piloto para destilación, purificación y formulación de electrolitos con construcción de vidrio de borosilicato, automatización PLC, HMI de pantalla táctil y características avanzadas de seguridad industrial para capacitación práctica en procesos químicos, ideal para planes de estudio de ingeniería química y ciencia de materiales.

Planta Piloto de Operaciones Unitarias Educativa de Absorción y Desorción

Planta Piloto de Operaciones Unitarias Educativa de Absorción y Desorción

Planta piloto de absorción y desorción con columna empaquetada dual para formación en ingeniería química, que ofrece medición en tiempo real del coeficiente de transferencia de masa, estructura móvil duradera, interfaz industrial de pantalla táctil y diseño personalizable para diversos currículos de laboratorio, que permite el estudio práctico de la absorción y desorción de gases.

Planta Piloto Educativa de Escala de Banco de Doble Columna para Separación y Captura de Gas

Planta Piloto Educativa de Escala de Banco de Doble Columna para Separación y Captura de Gas

Esta planta piloto educativa de doble columna proporciona enseñanza práctica sobre los procesos de adsorción, separación y captura de gases. Cuenta con columnas de acero inoxidable, regeneración de hasta 400 °C y un PLC con pantalla táctil de 15,6 pulgadas para estudios de TSA y PSA en planes de estudios de ingeniería química y simulación de procesos.

Planta Piloto Educativa para Demostración de Hidrodinámica de Columnas de Platos

Planta Piloto Educativa para Demostración de Hidrodinámica de Columnas de Platos

Planta piloto educativa avanzada transparente para laboratorios de ingeniería química que demuestra la hidrodinámica de columnas de platos con bandejas industriales de tamiz, tapa de burbuja y válvula dentada para la observación visual del contacto gas-líquido, la medición de caída de presión y el análisis de límites operativos, incluyendo inundación, goteo y arrastre

Planta Piloto de Formación Práctica para el Refino de Etanol Anhidro Verde

Planta Piloto de Formación Práctica para el Refino de Etanol Anhidro Verde

Planta piloto integrada avanzada para laboratorios universitarios que demuestra la destilación extractiva para producir etanol absoluto de alta pureza a partir de materias primas crudas, con operación continua de múltiples columnas, reciclaje de disolvente en circuito cerrado y controles personalizables para la educación en ingeniería práctica, ideal para la formación y la investigación en ingeniería química.

Planta Piloto de Operaciones Unitarias para la Síntesis de Acetato de Etilo para Formación Práctica

Planta Piloto de Operaciones Unitarias para la Síntesis de Acetato de Etilo para Formación Práctica

Planta piloto modular y personalizable para la formación práctica en la síntesis de acetato de etilo. Integra operaciones unitarias de reacción de esterificación, extracción líquido-líquido, neutralización y destilación de platos perforados. Conecta la teoría con los procesos industriales reales. Diseñada para laboratorios de ingeniería química universitarios.

Planta Piloto Educativa de Operaciones Unitarias de Cristalización por Membrana Multifuncional

Planta Piloto Educativa de Operaciones Unitarias de Cristalización por Membrana Multifuncional

Planta piloto integrada de cristalización por membrana a escala de banco para educación en ingeniería. Proporciona formación práctica en tecnologías de separación avanzadas, combinando cristalización por destilación de membrana e intensificación de procesos. Cuenta con vasos de escalado variables, control de flujo de grado industrial y adquisición de datos digital interactiva. Personalizable para laboratorios universitarios.


Deja tu mensaje