Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Scribd Logo
Cargar

¡Te damos la bienvenida a Scribd!

  • 124 vistas

    Ley de Enfriamiento de Newton

    Cargado por

    Veronica Beatriz
    El documento describe dos métodos para modelar el enfriamiento de una taza de café: el método analítico y el método de espacio de estado. El método analítico aplica la ley de enfriamiento de Newton para derivar una solución analítica a la ecuación diferencial que describe cómo cambia la temperatura con el tiempo. El método de espacio de estado formula un modelo matemático usando un balance de energía térmica y lo resuelve numéricamente usando el método de Euler. Ambos métodos muestran que la temperatura disminuye con el tiempo de forma ex

    Copyright:

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

    Ley de Enfriamiento de Newton

    Cargado por

    Veronica Beatriz
    124 vistas9 páginas
    El documento describe dos métodos para modelar el enfriamiento de una taza de café: el método analítico y el método de espacio de estado. El método analítico aplica la ley de enfriamiento de Newton para derivar una solución analítica a la ecuación diferencial que describe cómo cambia la temperatura con el tiempo. El método de espacio de estado formula un modelo matemático usando un balance de energía térmica y lo resuelve numéricamente usando el método de Euler. Ambos métodos muestran que la temperatura disminuye con el tiempo de forma ex

    Derechos de autor

    © © All Rights Reserved

    Formatos disponibles

    DOCX, PDF, TXT o lea en línea desde Scribd

    ¿Le pareció útil este documento?

    ¿Este contenido es inapropiado?

    El documento describe dos métodos para modelar el enfriamiento de una taza de café: el método analítico y el método de espacio de estado. El método analítico aplica la ley de enfriamiento de Newton para derivar una solución analítica a la ecuación diferencial que describe cómo cambia la temperatura con el tiempo. El método de espacio de estado formula un modelo matemático usando un balance de energía térmica y lo resuelve numéricamente usando el método de Euler. Ambos métodos muestran que la temperatura disminuye con el tiempo de forma ex

    Copyright:

    © All Rights Reserved
    Descargue como DOCX, PDF, TXT o lea en línea desde Scribd
    Descargar como docx, pdf o txt
    124 vistas9 páginas

    Ley de Enfriamiento de Newton

    Cargado por

    Veronica Beatriz
    El documento describe dos métodos para modelar el enfriamiento de una taza de café: el método analítico y el método de espacio de estado. El método analítico aplica la ley de enfriamiento de Newton para derivar una solución analítica a la ecuación diferencial que describe cómo cambia la temperatura con el tiempo. El método de espacio de estado formula un modelo matemático usando un balance de energía térmica y lo resuelve numéricamente usando el método de Euler. Ambos métodos muestran que la temperatura disminuye con el tiempo de forma ex

    Copyright:

    © All Rights Reserved
    Descargue como DOCX, PDF, TXT o lea en línea desde Scribd
    Descargar como docx, pdf o txt
    Está en la página 1de 9

    Mtodo Analtico

    Ley de Enfriamiento de Newton:


    La ley de enfriamiento de Newton establece, que la rapidez de cambio de temperatura
    de un cuerpo en cualquier tiempo t es proporcional a la diferencia de las temperaturas
    del cuerpo y del medio circundante en el tiempo t. Si consideramos a T como la
    temperatura del cuerpo en el tiempo t, T a la temperatura del medio circundante y T 0 la
    temperatura inicial del cuerpo (t = 0).
    Como las variaciones de la temperatura puede ser que aumenta o disminuya. Por
    tanto de acuerdo a la ley de enfriamiento de Newton se expresa mediante la ecuacin
    diferencial:
    (

    Separando variables:
    (

    Integrando esta ecuacin con la condicin inicial de que en el instante t = 0, la


    temperatura del cuerpo es T0:

    (
    (

    )
    )
    )

    Despejando T se obtiene la solucin analtica:


    (

    Donde:
    : Temperatura inicial de la taza.
    : Temperatura del medio que enfra (temperatura del aire).
    : Constante de enfriamiento.
    Datos:
    D (m)
    h (m)
    (kg/m3)
    Cp [J/(kgC)]
    A (m2)
    V (m3)
    U [J/(sm3C)]

    0,08
    0,07
    1000
    4187
    0,028
    3,519E-04
    27,138

    Donde:
    D: Dimetro de la taza de caf.
    h: Altura de lquido en la taza.
    : Densidad del caf, en este caso se consider igual que la del agua.
    CP: Calor especifico del caf, en este caso se consider igual que la del agua.
    A: rea de transferencia

    V: Volumen del lquido en la taza.

    U: Coeficiente de transmisin de calor


    Parmetros del modelo

    t (min)
    0
    2
    4
    6
    8
    10
    12
    14
    16
    18
    20
    22
    24
    26
    28
    30
    32
    34
    36
    38
    40
    42
    44
    46

    T0 C

    82,3

    Ta C
    r (1/s)

    17
    5,093E-04

    t (s) T (C) Calculada Analticamente


    0
    82,300
    120
    78,429
    240
    74,787
    360
    71,362
    480
    68,139
    600
    65,107
    720
    62,256
    840
    59,573
    960
    57,049
    1080
    54,675
    1200
    52,441
    1320
    50,340
    1440
    48,364
    1560
    46,505
    1680
    44,756
    1800
    43,110
    1920
    41,562
    2040
    40,106
    2160
    38,736
    2280
    37,448
    2400
    36,236
    2520
    35,095
    2640
    34,023
    2760
    33,014

    Curva de enfriamiento calculada por el Mtodo


    Analtico

    90
    80
    70

    TEMPERATURA CALCULADA
    ANALTICAMENTE

    T (C)

    60
    50
    40
    30
    20
    10
    0
    0

    10

    20

    30

    40

    50

    t (min)

    Mtodo de espacio de estado


    Definicin del objetivo de estudio: Predecir como vara la temperatura de la taza de
    caf en funcin del tiempo.
    Frontera que encierra el sistema: Volumen del lquido contenido en la taza.
    Formulacin del Modelo conceptual:
    Formulacin del Modelo matemtico:
    Balance de Energa Trmica

    Donde:

    Por lo tanto tenemos:

    La Transferencia de calor (Q) expresado en funcin del coeficiente global de


    transmisin de calor (U) est dado segn la relacin:

    Donde:

    Incgnitas
    Q, T y T

    Parmetros
    , CP, A, V, U, Ta, T0

    Variable de Entrada
    t

    Datos:
    D (m)
    h (m)
    r(1/s)

    0,08
    0,07
    5,093E-04

    Donde:
    D: Dimetro de la taza de caf.
    h: Altura de lquido en la taza.
    Determinacin de parmetros:
    Parmetros: , , A, V, U,
    y
    Donde:
    : Densidad del caf, en este caso se consider igual que la del agua.
    : Calor especifico del caf, en este caso se consider igual que la del agua.
    A: rea de transferencia

    V: Volumen del lquido en la taza.

    U: Coeficiente de transmisin de calor


    : Temperatura inicial de la taza.
    : Temperatura del medio que enfra (temperatura del aire).
    Parmetros del modelo
    (kg/m3)
    Cp [J/(kgC)]
    A (m2)
    V (m3)
    U [J/(sm3C)]

    1000
    4187
    0,028
    3,519E-04
    27,138

    Ta (C)
    T0 C

    17
    82,3

    Simplificacin del Modelo por eliminacin de variables:

    Anlisis de consistencia matemtica


    Grado de libertad del modelo
    Donde:
    : Nmero de incgnitas del modelo
    : Nmero de incgnitas del modelo
    Incgnitas
    T

    Parmetros
    , CP, A, V, U, Ta, T0

    Variable de Entrada
    t

    Resolucin del modelo por el mtodo de Euler en Excel


    t (min) t (s) T (C) Calculada por Euler
    0
    0
    82,300
    2 120
    78,309
    4 240
    74,563
    6 360
    71,045
    8 480
    67,742
    10 600
    64,641
    12 720
    61,730
    14 840
    58,996
    16 960
    56,430
    18 1080
    54,020

    Q
    48,992
    45,998
    43,187
    40,548
    38,070
    35,743
    33,559
    31,508
    29,583
    27,775

    20
    22
    24
    26
    28
    30
    32
    34
    36
    38
    40
    42
    44
    46

    1200
    1320
    1440
    1560
    1680
    1800
    1920
    2040
    2160
    2280
    2400
    2520
    2640
    2760

    51,758
    49,634
    47,640
    45,767
    44,009
    42,359
    40,809
    39,354
    37,988
    36,705
    35,501
    34,370
    33,309
    32,312

    26,077
    24,484
    22,988
    21,583
    20,264
    19,025
    17,863
    16,771
    15,746
    14,784
    13,881
    13,032
    12,236
    11,488

    Resolucin del modelo por el mtodo de Euler en mathcad


    Inicio
    Dat os
    D 0.08m
    h 0.07m

    Parmetros del Modelo


    1000

    C 1K

    kg
    3

    Cp 4.18710

    kg C

    A T 2

    D h

    4
    2

    VT

    D h

    VT 3.519 10

    J
    U 27.138
    2
    s m C

    Ta 17C
    T0 82.30C
    Modelo

    Q(T) U AT T Ta
    dT ( T)

    Q( T)
    Cp VT

    A T 0.028m

    Simulacin

    tf 46min
    NP 24

    tf

    t 2 min

    NP 1

    i 0 NP 1

    t i t
    i

    T T0
    0
    T

    i 1

    T dT T t
    i

    82.3

    78.309

    74.563

    71.045

    67.742

    64.641

    61.73

    58.996

    56.43

    54.02

    10

    51.758

    11

    49.634

    12

    47.64

    13

    45.767

    14

    44.009

    15

    42.359

    16

    40.809

    17

    39.354

    18

    37.988

    19

    36.705

    20

    35.501

    21

    34.37

    22

    33.309

    23

    32.312

    24

    31.376

    90
    80
    70
    60
    50
    40
    30
    3

    110

    210

    310

    ti

    Grfica de la evolucin obtenida en forma analtica y la simulada con Excel para los
    primeros 46 min

    Comparacin de los distintos mtodos

    90

    TEMPERATURA CALCULADA
    ANALTICAMENTE
    TEMPERATURA CALCULADA
    POR EULER EN EXCEL

    80
    70
    60
    T (C)

    Ti

    50
    40
    30
    20
    10
    0
    0

    10

    20

    30

    40

    50

    t (min)

    Determinacin del error porcentual entre la solucin analtica y la solucin numrica


    lograda con Excel.

    t (min)
    0
    2
    4
    6
    8
    10
    12

    T (C)
    T (C)
    Calculada
    Calculada por
    Error
    Analticamente Mtodo de Euler
    porcentual
    82,300
    82,300
    0,00
    78,429
    78,309
    -0,15
    74,787
    74,563
    -0,30
    71,362
    71,045
    -0,44
    68,139
    67,742
    -0,58
    65,107
    64,641
    -0,72
    62,256
    61,730
    -0,84

    14
    16
    18
    20
    22
    24
    26
    28
    30
    32
    34
    36
    38
    40
    42
    44
    46

    59,573
    57,049
    54,675
    52,441
    50,340
    48,364
    46,505
    44,756
    43,110
    41,562
    40,106
    38,736
    37,448
    36,236
    35,095
    34,023
    33,014

    58,996
    56,430
    54,020
    51,758
    49,634
    47,640
    45,767
    44,009
    42,359
    40,809
    39,354
    37,988
    36,705
    35,501
    34,370
    33,309
    32,312

    -0,97
    -1,08
    -1,20
    -1,30
    -1,40
    -1,50
    -1,59
    -1,67
    -1,74
    -1,81
    -1,88
    -1,93
    -1,98
    -2,03
    -2,07
    -2,10
    -2,12

    Determinacin del error porcentual entre la solucin analtica y la solucin numrica


    lograda con Mathcad.

    t (min)
    0
    2
    4
    6
    8
    10
    12
    14
    16
    18
    20
    22
    24
    26
    28
    30
    32
    34
    36
    38
    40
    42
    44
    46

    T (C)
    T (C)
    Calculada
    Calculada por
    Error
    Analticamente Mtodo de Euler porcentual
    82,300
    82,300
    0,00
    78,429
    78,309
    -0,15
    74,787
    74,563
    -0,30
    71,362
    71,045
    -0,44
    68,139
    67,742
    -0,58
    65,107
    64,641
    -0,72
    62,256
    61,730
    -0,84
    59,573
    58,996
    -0,97
    57,049
    56,430
    -1,08
    54,675
    54,020
    -1,20
    52,441
    51,758
    -1,30
    50,340
    49,634
    -1,40
    48,364
    47,640
    -1,50
    46,505
    45,767
    -1,59
    44,756
    44,009
    -1,67
    43,110
    42,359
    -1,74
    41,562
    40,809
    -1,81
    40,106
    39,354
    -1,88
    38,736
    37,988
    -1,93
    37,448
    36,705
    -1,98
    36,236
    35,501
    -2,03
    35,095
    34,370
    -2,07
    34,023
    33,309
    -2,10
    33,014
    32,312
    -2,13

    También podría gustarte