Scribd Logo
Cargar

¡Te damos la bienvenida a Scribd!

  • 3 vistas

    Pra2 21 22

    Cargado por

    Emilio Rodríguez

    Copyright:

    © All Rights Reserved
    Descargue como PDF, TXT o lea en línea desde Scribd

    Pra2 21 22

    Cargado por

    Emilio Rodríguez
    3 vistas3 páginas

    Título original

    pra2_21_22

    Derechos de autor

    © © All Rights Reserved

    Formatos disponibles

    PDF, TXT o lea en línea desde Scribd

    ¿Le pareció útil este documento?

    ¿Este contenido es inapropiado?

    Copyright:

    © All Rights Reserved
    Descargue como PDF, TXT o lea en línea desde Scribd
    Descargar como pdf o txt
    3 vistas3 páginas

    Pra2 21 22

    Cargado por

    Emilio Rodríguez

    Copyright:

    © All Rights Reserved
    Descargue como PDF, TXT o lea en línea desde Scribd
    Descargar como pdf o txt
    de 3

    Departamento de Matemática Aplicada

    Facultade de Matemáticas
    Rúa Lope Gómez de Marzoa, s/n
    15782 Santiago de Compostela

    Ecuaciones Diferenciales. Grado en Ingenierı́a Quı́mica. Curso 2021-22


    Práctica 2: Aplicaciones de las ecuaciones diferenciales

    En esta práctica se repasa la resolución de ecuaciones diferenciales de primer orden mediante el co-
    mando dsolve de Matlab.

    Ejercicios

    1. Se considera la reacción en cadena


    1 k 2 k
    A −→ B −→ C

    donde inicialmente hay Ao gramos del compuesto A. Si y(t) denota la cantidad de sustancia B
    en el instante t, se sabe que la ecuación diferencial que rige el proceso es
    dy
    = k1 Ao e−k1 t − k2 y
    dt
    donde k1 = 1 s−1 , k2 = 1/2 s−1 y Ao = 2 gr.
    Se pide:

    a) Clasificar la ecuación diferencial y plantear el problema de valor inicial que verifica y(t).
    b) Resolver el problema de valor inicial mediante Matlab.
    c) Calcular la cantidad de compuesto A, B y C que hay en cada instante de tiempo y repre-
    sentarlas gráficamente.

    2. Un circuito RL o RC consiste en una fuente de alimentación de fuerza electromotriz E(t), una


    resistencia de valor R y un condensador de capacidad C (circuito RC) o bien una bobina de
    inductancia L en serie (circuito RL), según se indica en la Figura 1. La ecuación diferencial que
    modela la corriente I en estos circuitos viene dada por la segunda ley de Kirchhoff: “la suma
    de las caı́das de tensión o voltajes en los elementos de un circuito debe ser igual a la fuerza
    electromotriz aplicada”.
    En el caso del circuito RL, la EDO resultante es
    dI
    L + RI = E(t).
    dt
    Por otro lado, para un circuito RC, se tiene que
    q
    RI + = E(t),
    C
    donde la corriente en el condensador es el ratio de cambio de su carga, es decir,
    dq
    I(t) = .
    dt

    1
    Figura 1: Circuito RL (izquierda) y circuito RC (derecha). Fuente: Nagle, R., Saff, E. y Snider, A.
    (2017) Fundamentals of differential equations and boundary value problems. Reading, MA: Pearson.

    De este modo, se obtiene la siguiente ecuación:


    dq q
    R + = E(t).
    dt C

    Se pide:

    a) Supón que un circuito RL está formado por una resistencia de 1 ohmio, un inductor de 0.01
    henrios y que la fuerza electromotriz es E(t) = sin(100t) V. Si la corriente inicial es cero,
    obtén la intensidad de corriente para t > 0.
    b) Considera ahora un circuito RC con la misma resistencia, pero con un condensador de
    capacidad 0.2 faradios en lugar del inductor y sujeto a una fuerza electromotriz igual a
    100 voltios. Obtén la intensidad de corriente I(t) y la carga q(t) si la carga inicial es de 0
    culombios. Representa gráficamente ambas funciones en dos figuras diferentes.

    3. Los gases son fluidos compresibles cuya densidad es casi proporcional a la presión. Teniendo en
    cuenta la ley de los gases perfectos
    p = ρRT
    en la EDO
    dp
    = −ρg,
    dz
    se obtiene
    dp g
    =− p,
    dz RT
    siendo p la presión, z la altura, ρ la densidad, g la gravedad, R la constante de los gases perfectos
    y T la temperatura.
    Sabiendo que la temperatura atmosférica disminuye linealmente según una ecuación de la forma

    T (z) = T0 − B z,

    hasta una altura de 11000 m y empleando que la presión a nivel del mar (z = 0) es 101.35 kPa,
    se pide:

    2
    a) Plantear el PVI cuya solución es p(z) y resolverlo con Matlab en función de las variables
    simbólicas T0 , B, g y R, da forma que la variable independiente sea z.
    b) Simplificar la expresión obtenida.
    c) Teniendo en cuenta que T0 = 288.16 K, B = 0.0065 K/m, g = 9.807 m/s2 y R = 287
    m2 /(s2 ·K), sustituir en la expresión anterior las variables simbólicas por su valor.
    d ) Representar las funciones presión y temperatura de 0 a 11000 m en la misma gráfica. Desde
    la ventana de comandos, pon tı́tulo y leyenda a la gráfica y representa la presión con trazo
    rojo y discontinuo.

    También podría gustarte