En este trabajo se presenta un modelado termodinámicamente consistente del proceso de destilación multicomponente, considerando el volumen de fases constante en cada etapa. Se generalizan los balances de moles totales, balances por componente y balances de energía, formulando las matrices de masa para sistemas ideales que cumplen con la ley de Raoult. Las simulaciones implementadas en MATLAB muestran resultados coherentes, con estados estacionarios que convergen a valores cercanos a los obtenidos mediante métodos clásicos como Fenske-Underwood-Gilliland. El objetivo principal es analizar el comportamiento dinámico y optimizar el desempeño energético de una columna de destilación, partiendo de un modelo detallado para mezclas binarias en condiciones adiabáticas. Además, se realiza un estudio de las propiedades termodinámicas, tales como la energía total, entropía y energía libre de Gibbs, para identificar los procesos fenomenológicos relevantes. En etapas futuras, se planea implementar una optimización basada en la minimización de la producción de entropía, con el fin de mejorar la eficiencia energética del proceso. Este enfoque evita las simplificaciones comunes, tales como flujos internos constantes, desprecio de la dinámica en la fase gaseosa o entalpías constantes, las cuales limitan la modelación a sistemas de dos componentes, permitiendo así abordar sistemas multicomponentes más complejos. El trabajo constituye una contribución significativa hacia la modelación precisa y optimizada de procesos de destilación multicomponente.