Source: http://slideplayer.es/slide/1716062/
Timestamp: 2017-02-27 06:36:59+00:00

Document:
Modelado y simulación en Ingeniería Química. Manuel Rodríguez Tema 4: Simulación dinámica. - ppt descargar
Modelado y simulación en Ingeniería Química. Manuel Rodríguez Tema 4: Simulación dinámica.
Publicada porCarmelita Morado
Presentación del tema: "Modelado y simulación en Ingeniería Química. Manuel Rodríguez Tema 4: Simulación dinámica."— Transcripción de la presentación:
Modelado y simulación en Ingeniería Química. Manuel Rodríguez Simulación dinámica Simuladores orientados a ecuaciones Métodos de Resolución Introducción a Abacuss ÍNDICE
Modelado y simulación en Ingeniería Química. Manuel Rodríguez Simulación dinámica Operación: JRespuesta de procesos continuos ante perturbaciones JAjuste de controladores JManiobras y desviaciones anormales en el proceso, para estudios de seguridad y de emisiones JAnálisis de operabilidad y riesgo JValidación de procedimientos de emergencia JEntrenamiento de operadores Diseño: JSistema de control y controlabilidad JProcedimientos de puesta en marcha y parada JProcesos discontinuos Estudia el comportamiento transitorio de un sistema. Simula su evolución temporal con distintos fines:
Modelado y simulación en Ingeniería Química. Manuel Rodríguez 1.Las ecuaciones y variables de todos los modelos de unidades se definen individualmente. 2.Las ecuaciones y variables conforman un sistema (grande) de ecuaciones no lineales. 3.Se añaden especificaciones hasta que el sistema tiene grados de libertad CERO. 4.El sistema de ecuaciones se resuelve de forma directa y simultánea. Pasos para establecer la simulación: Simuladores orientados a ecuaciones Es la estrategia empleada para la simulación dinámica.
Modelado y simulación en Ingeniería Química. Manuel Rodríguez Ventajas: Es más eficiente que la estrategia secuencial modular. No distingue entre simulación y diseño de especificación. Es más fácil hacer una librería de modelos extensible y reutilizable. Es utilizable en simulación dinámica, estacionaria y optimización. Es más fácil el diagnóstico de ciertos errores. Como sistemas de ecuaciones mal especificados. Inconvenientes (en la simulación estacionaria) Los algoritmos de resolución de NAEs no son tan robustos y fiables como los del caso secuencial modular. Necesita más recursos de computación (memoria principalmente). Debido al primer inconveniente esta estrategia no es actualmente competitiva frente a la resolución secuencial modular para simulación estacionaria. Orientado a ecuaciones vs. Secuencial modular Simuladores orientados a ecuaciones
Modelado y simulación en Ingeniería Química. Manuel Rodríguez  Seleccionar variables x : Variables de estado u : Variables independientes (generalmente variables de entrada) v : Variables dependientes resultado de ecuaciones algebraicas  Ordenar ecuaciones v 1 =g 1 (x,u,t) v 2 =g 2 (x,u,v 1,t).... v i =g i (x,u,v 1 v 2,...,v i-1,t) x ´=f(x,u,v,t)  Inicializar la integración Elección del paso h Valores iniciales de variables x 0,u 0,v 0 Valores iniciales de f  Bucle de simulación Calcular f r =f(x r,u r,v r,t) Incrementar el tiempo t=t+h Paso de integración x r+1 función de f r,f r-1,... Posible actualización de u r+1 Obtener ordenadamente v i,r+1 para I=1,2,... Cada cierto tiempo imprimir o recolectar resultados Verificar condiciones de fin de simulación Algoritmo de simulación de sistemas continuos Simuladores orientados a ecuaciones
Modelado y simulación en Ingeniería Química. Manuel Rodríguez ECUACIONES DIFERENCIALES ORDINARIAS (ODEs) Métodos de paso único x´(t)= (x(t),t) x(t)=x o + t0 t (x(t),t)dt Son métodos que emplean únicamente el valor calculado en el punto (paso) anterior. Se basan en aproximar la función f(x,t) mediante un desarrollo en serie de Taylor. Métodos de resolución
Modelado y simulación en Ingeniería Química. Manuel Rodríguez Euler implícito Euler explícito x(t+h)=x(t)+h dx(t)/dt dx(t)/dt = (x(t),t) x(t+h)=x(t)+h (x(t),t) x(t+h)=x(t)+h (x(t+h),t) h es el paso de integración No es siempre estable (puede oscilar o divergir para pasos pequeños) Métodos de resolución
Modelado y simulación en Ingeniería Química. Manuel Rodríguez Ejemplo function [x,y] = eulerm(f,tspan,y0,n), % a = tspan(1); b = tspan(2); h = (b-a)/n; t = (a+h:h:b); x(1) = x0 + h*feval(f,a,x0); for i=2:n x(i) = x(i-1) + h*feval(f,t(i-1),x(i-1)); end t = [a t]; x = [x0 x]; Euler explícito Métodos de resolución Modelar mediante Matlab el vaciado de un tanque. Área 100 Flujo de entrada 3 Altura inicial 1m Área salida 1 Simular el comportamiento del tanque empleando el algoritmo de Euler Explícito.
Modelado y simulación en Ingeniería Química. Manuel Rodríguez Runge kutta (orden 45) k 1 =h (t i,x i ) k 2 =h (t i +h/2,x i + k 1 /2) k 3 =h (t i +h/2,x i + k 2 /2) k 4 =h (t i +h/2,x i + k 3 ) x i+1 = x i +k 1 /6 +k 2 /3 +k 3 /3 +k 4 /6 Algoritmo de orden 4 (número de términos) y orden del error 5. Funciones matlab ode23 y ode45 Métodos de resolución
Modelado y simulación en Ingeniería Química. Manuel Rodríguez function [t,x] = rk4m(f,tspan,x0,n), a = tspan(1); b = tspan(2); h = (b-a)/n; t = (a+h:h:b); k1 = h*feval(f,a,x0); k2 = h*feval(f,a+h/2,x0+k1/2); k3 = h*feval(f,a+h/2,x0+k2/2); k4 = h*feval(f,a+h,x0+k3); x(1) = x0 + (k1)/6+ (k2)/3 + (k3)/3 + (k4)/6; for i=1:n-1 k1 = h*feval(f,t(i),x(i)); k2 = h*feval(f,t(i)+h/2,x(i)+k1/2); k3 = h*feval(f,t(i)+h/2,x(i)+k2/2); k4 = h*feval(f,t(i)+h,x(i)+k3); x(i+1) = x(i) + (k1)/6+ (k2)/3 + (k3)/3 + (k4)/6; end t = [a t]; x = [x0 x]; Runge Kutta 45 Ejemplo Métodos de resolución Simular el comportamiento del tanque anterior empleando el algoritmo de Runge Kutta 45. Comparar los resultados con los proporcionados por el algoritmo de Euler según se disminuye el paso de integración. Modelar mediante Matlab dos tanques idénticos en serie. Área 100 Flujo de entrada tanque 1 = 2.2 Flujo de entrada tanque 2 = 1.6 Altura inicial tanque 1 =0.8m Altura inicial tanque 2 =0.2m Simular el comportamiento empleando el algoritmo de Runge Kutta 45.
Modelado y simulación en Ingeniería Química. Manuel Rodríguez Métodos multipaso Adams-Bashford Adams-Moulton x i+1 = x i +h/2 [3 (t i,x i )- (t i-1,x i-1 )] x i+1 = x i +h/24 [9 (t i+1,x i+1 )+19 (t i,x i )- 5 (t i-1,x i-1 )- (t i-2,x i-2 )] Método explícito de 2 pasos Método implícito de 3 pasos Métodos de resolución
Modelado y simulación en Ingeniería Química. Manuel Rodríguez x* i+1 = x i +h/12 [23 (t i,x i )- 16 (t i-1,x i-1 )+ 5 (t i-2,x i-2 )] x i+1 = x i +h/12 [ (t i+1,x* i+1 )+8 (t i,x i )- (t i-1,x i- )] Predictor-Corrector Método explícito para computar x* i+1 Método implícito para computar el valor final x i+1 Evita tener que resolver una ecuación algebraica en el método implícito Métodos de resolución
Modelado y simulación en Ingeniería Química. Manuel Rodríguez ODEs de orden n x (n (t) = f(x(t),x´(t),x´(t),...,x (n-1 (t),t) Para cualquier sistema de ecuaciones diferenciales (de primer orden): RESULTAN VÁLIDOS TODOS LOS MÉTODOS DE INTEGRACIÓN ANTERIORES Un sistema de ODEs de orden n se puede transformar a un sistema de ODEs de primer orden. Métodos de resolución
Modelado y simulación en Ingeniería Química. Manuel Rodríguez ECUACIONES DIFERENCIALES ALGEBRAICAS (DAEs) Stiffness Método de resolución Aparece cuando hay constantes de tiempo muy diferentes en un sistema. Hay fenómenos muy rápidos y fenómenos lentos mezclados. Se presenta en muchos modelos de procesos de la industria química. Se deben emplear métodos implícitos (y de paso variable) DAEs: Conjunto de ecuaciones diferenciales y algebraicas. x' = f(x, u,y, t) g(x, u,y, t) = 0 Semi- explicita Dado x(tn) se resuelve g(x(tn), u(tn)) = 0 ==> se obtiene u(tn) Empleando un método explícito de resolución de ODE a x' = f(x, u, t) se obtiene x(tn+1) Resolver x' = f(x,u, t), g(x, u, t)=0 de forma simultánea mediante métodos implícitos (BDF). Por ejemplo: Método de Gear Forma secuencial Forma simultánea Métodos de resolución
Modelado y simulación en Ingeniería Química. Manuel Rodríguez Métodos implícitos BDF (Backward Differentiation Formulas) DAEs forma Semi-explícita Para cada paso de integración hay que resolver este sistema de ecuaciones algebraicas. Por Newton (u otro método de resolución)
Modelado y simulación en Ingeniería Química. Manuel Rodríguez Introducción a Abacuss II Escribiendo un modelo: Sintaxis # ------------------------ # El símbolo # indica el comienzo de un comentario # ----------------- MODELO TANQUE --------------- Editor de textos (Wordpad, notepad, word,...) DECLARE MODEL SIMULATION TYPE STREAM PARAMETER UNIT VARIABLE STREAM SELECTOR SET EQUATION INEQUALITY SELECTOR INITIAL SCHEDULE SET EQUATION INPUT PRESET OPTIONS UNIT REPORT ASK 3 Secciones principales Subsecciones END Acaban siempre con la palabra END
Modelado y simulación en Ingeniería Química. Manuel Rodríguez Introducción a Abacuss II Nombre = Valor por defecto : valor mín : valor máx [ UNIT = unidades ] Nombre IS variable1 [,variable2,variable3,...] La sección DECLARE Se especifican los tipos de variable a emplear en el modelo Se especifican los tipos de corriente (que tendrá alguna de las variables anteriores)
Modelado y simulación en Ingeniería Química. Manuel Rodríguez Introducción a Abacuss II Inlet: nombre1 [,nombre2,..] AS nombre_tipo_corriente (definido en DECLARE) Nombre1[,nombre2...] AS tipo_variable (definida en DECLARE) Nombre1[,nombre2...] AS REAL, INTEGER O LOGICAL Lista de ecuaciones empleando las variables definidas anteriormente. Las derivadas se denotan añadiendo $ delante de la variable. Admite construcciones con: FOR... TO....END IF... THEN... ELSE..END CASE.... SWITCH La sección MODEL
Modelado y simulación en Ingeniería Química. Manuel Rodríguez Introducción a Abacuss II Instancias de los modelos creados Fija valores de los PARÁMETROS Fija valores de las VARIABLES de entrada Ecuaciones ( o asignaciones) para indicar los valores iniciales para las variables que aparecen derivadas. Indica la evolución de la simulación. RESET cambia los valores de variables de entrada en determinados tiempos. La sección SIMULATION
Modelado y simulación en Ingeniería Química. Manuel Rodríguez Introducción a Abacuss II Ejecutando ABACUSSII Carga un modelo para la simulación Selecciona una simulación para ejecutarla Limpia la pantalla Sale del programa Muestra tipos de variables, modelos, corrientes,... están cargados. Análisis de los bloques que se forman para la resolución matemática
Modelado y simulación en Ingeniería Química. Manuel Rodríguez Introducción a Abacuss II Cargando un modelo...
Modelado y simulación en Ingeniería Química. Manuel Rodríguez Introducción a Abacuss II Ejecutando el modelo...
Modelado y simulación en Ingeniería Química. Manuel Rodríguez Introducción a Abacuss II Resultados en pantalla
Modelado y simulación en Ingeniería Química. Manuel Rodríguez Introducción a Abacuss II Resultados en archivo Excel I
Modelado y simulación en Ingeniería Química. Manuel Rodríguez Introducción a Abacuss II Resultados en archivo Excel y II
Descargar ppt "Modelado y simulación en Ingeniería Química. Manuel Rodríguez Tema 4: Simulación dinámica."
Presentaciones similares Modelado y simulación en Ingeniería Química. Manuel Rodríguez Tema 4: Simulación dinámica.
FORTRAN 77 y Cómputo Científico Guilmer González Facultad de Ciencias – UNAM Programación Avanzada para Ingenieros Petroleros asilados en la Facultad de.
/35 Identificacion de sistemas Simulacion de sistemas continuos 1.
/64 Simulacion de sistemas dinamicos Principios básicos de la integración numérica 1.
ESTABILIDAD EN TENSION. OSCILOGRAMA TIPICO METODOS DE ANALISIS DE LA ESTABILIDAD DE TENSION A) Estáticos - Resolución de sistemas de ecuaciones algebraicas.

References: Resolución 
 resolución 
 resolución 
 resolución

 resolución

 resolución 
 resolución

 resolución 
 resolución

 resolución

 resolución

 resolución 
 resolución 
 resolución

 resolución 
 Resolución