Balances Microscopicos de Calor(1)

download Balances Microscopicos de Calor(1)

of 11

Transcript of Balances Microscopicos de Calor(1)

  • 7/24/2019 Balances Microscopicos de Calor(1)

    1/11TECNOLOGA QUMICAVol. XXVI, No. 2, 2006 53

    BALANCES MICROSCPICOS DE ENERGIA CALORFICAEN COORDENADAS CILNDRICAS APLICANDOECUACIONES DISCRETAS OBTENIDAS EN UN

    VOLUMEN DE CONTROL

    Armando A. Daz Garca

    Universidad de oriente

    En el presente trabajo se expone un mtodo para determinar numricamente perfi les de temperatura

    en un paraleleppedo utilizando diferencias finitas.

    El mtodo se basa en balances de energa calorfica en volmenes de control para los distintos

    sistemas de puntos en que se divide al modelo, utilizando diferencias finitas y sin utilizar ecuaciones

    en derivadas parciales.

    Se toma como modelo fsico para desarrollar el mtodo un cuerpo cilndrico que se encuentra

    inicialmente a una temperatura determinada y se introduce abruptamente en un fluido caliente

    producindose un calentamiento gradual.

    Se explica el mtodo y se obtienen las ecuaciones que pueden ser utilizadas para predecir la

    temperatura en algn punto del cilindro en cualquier tiempo deseado.

    Palabras clave: balances de energa calorfica, ecuaciones discretas.

    In the present work is exposed a method to determine profiles of temperature numerically in a

    cylinder using finite differences.

    The method is based on balances of heating energy in control volumes for the different systems of

    points in that is divided to pattern using finite differences and without using equations in having

    part ial derivatives.

    It is taken as physical model to develop the method a cylindrical body that is initially at a certain

    temperature, and it is introduced abruptly in a hot fluid taking place a gradual heating.

    The method is explained, and are obtained the equations that can be used to predict the temperature

    in any point of the cylinder in any wanted time.

    Key words: heating energy balances, discret equetions.

    _____________________

    Introduccin

    En un artculo de igual nombre, pero aplicado

    a un sistema descrito en coordenadas rectangula-

    res, publicado en esta misma revista, se dieron los

    elementos bsicos indispensables relacionados

    con las diferencias finitas indispensables para

    comprender cada uno de los pasos necesarios

    para llevar a cabo el balance microscpico.

    En este trabajo aplicaremos el mtodo general

    que permite, con conocimientos muy elementales

    de lgebra y geometra, llevar a cabo balances

    microscpicos de calor en un sistema cilndrico

    sometido a calentamiento o enfriamiento por

    conveccin.

    Generalidades

    Para explicar el mtodo utilizaremos un ejem-

    plo basado en un modelo fsico sencillo: un cilindro

    homogneo de longitud finita y propiedades fsi-

    cas constantes, sometido a calentamiento por

    conveccin.

    Supongamos un cuerpo slido homogneo en

    forma cilndrica de longitud 2L y radio R, que seencuentra inicialmente a una temperatura T0 y

    abruptamente (a partir de un tiempo t = 0) se

    introduce en una masa de fluido a la temperatura

    TEX T, constante, intercambiando calor por

    conveccin. Consideremos que los coeficientes

    de transferencia de calor en la pelcula en toda la

    superficie del slido son conocidos y constantes,

  • 7/24/2019 Balances Microscopicos de Calor(1)

    2/11TECNOLOGA QUMICAVol. XXVI, No. 2, 200654

    y se desea conocer la temperatura en cualquier

    punto del cuerpo en un tiempo t.

    Segn las caractersticas especificadas, el ci-

    lindro tendr simetra tanto geomtrica como

    fenomenolgica, y por lo tanto, slo tendremosque desarrollar las ecuaciones para describir la

    mitad del cilindro en cuestin, tal como se muestra

    en la figura 1.

    Puntos bases del eje central: L(0,m).Constituido por todos los puntos bases de la

    lnea central del cuerpo que intercambian calor

    por conduccin con dos puntos vecinos (ver

    esquema en la figura 3).

    Fig. 1 Cilindro slido simtricoexpuesto a calentamiento por

    conveccin.

    Sistemas

    Un sistema es, para los balances microscpi-

    cos, el conjunto de puntos bases de la red que

    pueden ser descritos por una misma ecuacin

    discreta. Para que distintos puntos bases formen

    parte de un sistema, es necesario que tengan la

    misma geometra e igual situacin fenomenolgica.Para el caso que nos ocupa se pueden distinguir

    los siguientes sistemas de puntos bases:

    Puntos bases interiores: V(n,m)

    Constituido por puntos bases que Intercambian

    solo calor por conduccin con cuatro puntos

    vecinos (ver esquema en la figura 2).

    Fig. 2 Esquema del sistema de puntosbases interiores.

    Fig. 3 Esquema del sistema de puntosde bases del eje central.

    Puntos bases de la superficie lateral: S(N,m).

    Constituido por todos los puntos bases de la

    superficie que intercambian calor por conduc-

    cin con tres puntos vecinos y por conveccin

    con un punto vecino ficticio (ver esquema en la

    figura 4).

    Fig. 4 Esquema del sistema de puntosbases de la superficie exterior.

    Puntos bases de la superficie inferior (fon-do):S(N,m).

    Constituido por todos los puntos bases de la

    superficie inferior que intercambia calor por

    conduccin con tres puntos vecinos y por

    conveccin con un punto vecino ficticio (ver

    esquema en la figura 5).

    L

    r

    Z

    TEX

    TEXT

    pT

    0,0

    pT

    0,1

    pT

    0,1

    pm

    T1,0 +

    pmT 1,0

    r

    pmNT ,

    pmNT ,1

    pmNT 1, +

    pmNT 1,

    P(0,0)

    TEXT

    pmnT ,

    p

    mnT

    ,1+

    p

    mnT ,1

    pmnT 1, +

    pmnT 1,

    P(0,0)

  • 7/24/2019 Balances Microscopicos de Calor(1)

    3/11TECNOLOGA QUMICAVol. XXVI, No. 2, 2006 55

    Puntos bases de lnea o borde: L(N,0)

    Constituido por todos los puntos bases del

    borde inferior de la superficie del fondo que

    ocupan la geometra tpica de una circunferen-

    cia y que intercambian calor por conduccin

    con dos puntos vecinos y por conveccin con

    dos puntos vecinos ficticios (ver esquema en la

    figura 6).

    Puntos vecinos

    Puntos vecinos de los puntos bases interiores

    Son los puntos interiores que se relacionan con

    los puntos bases segn se indica en la figura 2.

    Observe que para todos los puntos interiores el

    punto base tiene como puntos vecinos, los puntosp

    mnT ,1+ ,p

    mnT ,1 ,p

    mnT 1, + yp

    mnT 1, todos transfiriendo

    calor por conduccin.

    Teniendo en cuenta que la transmisin de calor

    por conduccin queda definida por las relaciones:

    ( ) ( )

    ( ) ( )pmnpmnpmnpmnZ

    p

    mn

    p

    mn

    p

    mn

    p

    mnr

    TTZ

    kTT

    Z

    kqc

    TTr

    kTT

    r

    kqc

    ,1,,1,

    ,,1,,1

    +

    =

    +

    =

    +

    +

    donde qc es la densidad de flujo de calor en los

    ejes cooordenados especificados, y que el flujo de

    calor es perpendicular a las superficies de trans-

    ferencia, no habr otros puntos vecinos

    involucrados que no sean los geomtricamente

    perpendiculares al punto base.

    Fig. 5 Esquema del sistema de puntosbases de la superficie exterior

    inferior (fondo).

    Fig. 6 Esquema del sistema de puntosbases del borde inferior (fondo).

    Punto base de centro (punto central del fon-do): P(0,0)

    Es similar a los puntos bases del eje central

    pero se diferencia en que presenta intercam-

    bio trmico por conveccin con un punto veci-

    no ficticio (ver esquema en la figura 7).

    Fig. 7 Esquema del sistema de puntosbases del centro en el fondo del eje

    central.

    pnT 0,

    pn

    T0,1

    pn

    T1,

    P(0,0)

    TEX

    pnT 0,1+

    pT

    0,0

    pT

    0,1

    pT

    0,1

    pT

    1,0

    EXTT

    r

    pNT 0,

    pNT 0,1

    pN

    T1,

    P(0,0)

    TEX

    EXTT

  • 7/24/2019 Balances Microscopicos de Calor(1)

    4/11TECNOLOGA QUMICAVol. XXVI, No. 2, 200656

    Puntos vecinos de los puntos bases de lasuperficie lateral

    Para un punto base en la superficie la relacin

    de los puntos vecinos se muestra en la figura 4. Enla misma se observa que a la derecha,en la

    direccin r no existe cuerpo material, sino contac-

    to con el fluido con el cual se establecer la

    tranferencia por conveccin

    Los puntos bases, en este caso para todos los

    puntos de la pared situados en el plano (N,m),

    pueden ser representados de forma general segn

    la figura 4 por p

    mNT , , siendo los puntos vecinos en

    contacto por conduccinp

    mN

    p

    mN TT 1,1, ; + yp

    mNT ,1 y

    en contacto por conveccin con el punto vecino

    ficticio TEXT. Observe que en este caso, todos los

    puntos estn en la posicin (N,m) correspondien-

    te a la pared situada enr

    RN

    = de acuerdo con

    la simetra. .

    Puntos vecinos de los puntos bases de lasuperficies inferior (fondo)

    La relacin con los puntos vecinos para un

    punto base situado en la superficie inferior se

    muestra en la figura 5. En la misma se observa

    que hacia abajo, en la direccin -Z no existe

    cuerpo material sino contacto con el fluido con el

    cual se establecer la tranferencia por conveccin.

    Los puntos bases, en este caso para todos los

    puntos de la pared situados en el plano (n,0),

    pueden ser representados de forma general segn

    la figura 5 porp

    nT 0, , siendo los puntos vecinos en

    contacto por conduccin

    p

    n

    p

    n TT 0,11, ; yp

    nT 0,1+ y elpunto vecino ficticio en contacto por conveccin

    TEXT .

    Puntos vecinos de los puntos bases de linea(bordes)

    Los puntos que conforman el borde de la

    superficie del fondo concuerdan tambin con un

    sistema de puntos pues, presentan transferencia

    de calor por conveccin por dos lados, y a su vez

    presentan conduccin de forma geomtricamente

    diferente de los de la superficie . De este modo,

    teniendo en cuenta que el cuerpo presenta sime-tra, habr slo un sistema de puntos para el borde,

    el correspondiente al borde de la superficie del

    fondo. En la figura 6 se muestra un esquema de la

    distribucin de los puntos vecinos en el borde

    inferior del cilindro, observe que presenta dos

    puntos vecinos por conduccin y dos puntos veci-

    nos ficticios en contacto por conveccin. El punto

    base, en este caso para todos los puntos del borde

    situados en la circunferencia L (n,0), viene dado

    de forma general por el representado en la figura

    6,

    p

    NT 0, siendo los puntos vecinos en contacto porconduccin

    p

    N

    p

    N TyT 1,0,1 .

    Puntos vecinos de los puntos bases de linea(eje central)

    Los puntos bases del eje central se caracteri-

    zan, al igual que los puntos interiores, por inter-

    cambiar calor con cuatro puntos vecinos por

    conduccin (ver figura 3) con la diferencia que

    todos los puntos vecinos en el radio presentan

    simetra fenomenolgica y geomtrica, y por tan-

    to son una misma variable, por lo que la ecuacinque los describe estar en funcin slo de tres

    puntos vecinos.

    Puntos vecinos de centro (punto centraldel fondo): P(0,0)

    Son similares a los puntos vecinos de la super-

    ficie del fondo pero se diferencia en que los

    puntos vecinos radiales son simtricos, corres-

    pondiendo a una misma variable discreta, por lo

    que la ecuacin que lo describe estar en funcin

    slo de tres puntos. El esquema de la distribucinde los puntos se muestra en la figura 7.

    Balances de calor en el volumen de controlBalance de calor no estacionario

    Para llevar a cabo el balance de calor no

    estacionario en el volumen de control, tendremos

    en cuenta que en cualquiera de los sistemas

    analizados se cumple que:

  • 7/24/2019 Balances Microscopicos de Calor(1)

    5/11TECNOLOGA QUMICAVol. XXVI, No. 2, 2006 57

    =

    controldevolumen

    delsaleque

    calordeFlujo

    controldevolumen

    alentraque

    calordeFlujo

    calorde

    flujodel

    Variacion (1)

    La variacin del flujo de calor viene dada por la expresin:

    (2)( )t

    TTViCp

    calorde

    flujodel

    Variacinp

    basepunto

    p

    basepunto

    =

    +1

    2

    y el flujo de calor que entra o sale del volumen de

    control es especfico del punto base de que se trate.

    Se toma la mitad del volumen de control teniendo en

    cuenta que el flujo de calor se produce en todo el

    volumen, y se supone que la acumulacin ocurra en

    la mitad del volumen de control y no en su totalidad.

    Balance de calor en puntos interiores

    Para el sistema de puntos bases interiores en

    que slo tiene lugar intercambio de calor por

    conduccin la ecuacin de balance viene dada

    por:

    ( ) ( )[ ]

    ( ) ( )[ ]pmnp mnpmnp mnZ

    p

    mn

    p

    mnrINT

    p

    mn

    p

    mnrEXT

    TTTTAZ

    k

    TTATTAr

    k

    calorde

    flujodel

    Variacin

    ,,1,,1

    ,,1,,1

    +

    +

    ++

    =

    +

    +

    (3)

    donde ArINT, ArEXTy Azson las reas de transferen-

    cia del volumen de control perpendiculares a la

    direccin del flujo, ya sea en r o en Z. Obsrvese de

    la existencia de dos reas en la direccin r, una

    interna y otra externa. Esto es debido a que el rea

    en la direccin r es variable en un cilindro, siendo una

    funcin del radio, por lo que posteriormente cuando

    se analicen las reas de transferencia habr que

    definir el valor de r donde se efecta el balance.

    Balance de calor en puntos bases de la

    superficie

    Para los puntos de la pared lateral, cuyos

    puntos vecinos se muestran en la figura 4, en

    los que tiene lugar transferencia de calor por

    conveccin y conduccin, la ecuacin de ba-

    lance ser segn el esquema mostrado en la

    figura 4:

    (4)

    ( ) ( )

    ( ) ( )[ ]p mNp mNp mNp mNZ

    p

    mNEXTrEXT

    p

    mN

    p

    mNrINT

    TTTTAZ

    k

    TThATTAr

    k

    calorde

    flujodel

    Variacin

    ,1,,1,

    ,,,1

    +

    +

    ++

    =

    +

    Balance de calor en puntos bases de la

    superficie del fondo

    Para los puntos del fondo del cilindro, cuyos

    puntos vecinos se muestran en la figura 5, en los

    que tiene lugar transferencia de calor por

    conveccin y conduccin, la ecuacin de ba-

    lance ser, segn el esquema mostrado en la

    figura mencionada:

  • 7/24/2019 Balances Microscopicos de Calor(1)

    6/11TECNOLOGA QUMICAVol. XXVI, No. 2, 200658

    Balance de calor en puntos bases del borde

    Los puntos bases del borde se asemejan enmucho a los puntos bases de la superficie delfondo, pero se diferencian en que presentan dos

    puntos vecinos ficticios ocupados por la tempera-tura externa. En la figura 6 se muestra un esque-ma para el sistema de puntos del borde en surelacin con los puntos vecinos.

    ( ) ( )[ ]

    ( ) ( )pnEXTZ

    p

    n

    p

    nZ

    p

    n

    p

    nrINT

    p

    n

    p

    nrEXT

    TThATTAZ

    k

    TTATTAr

    k

    calorde

    flujodel

    Variacin

    0,0,1,

    0,0,10,0,1

    +

    +

    +

    =

    +

    (6)

    Observe que hay dos puntos vecinos ficticiosdescritos por la ecuacin de conveccin.

    Balance de calor en los puntos bases del ejecentral

    Los puntos bases del eje central se diferenciande los puntos bases interiores en que en stosexisten dos puntos simtricos alrededor del ejeque son descritos por una sola variable. De estemodo, segn el esquema mostrado en la figura 3.

    ( ) ( )

    ( ) ( )pNEXTZpNpNZ

    p

    NEXTrEXT

    p

    N

    p

    NrINT

    TThATTAZ

    k

    TThATTAr

    k

    calorde

    flujodel

    Variacin

    0,0,1,

    0,0,0,1

    +

    +

    ++

    =

    (7)

    Como los puntos vecinosp

    m

    p

    m TyT ,1,1 son iguales y las reas rINTArEXTA = la ecuacin se puede escribir

    de la siguiente manera:

    (8)( )[ ] ( ) ( )[ ]pmpmpmpmZpmpmrINT TTTTA

    Z

    kTTA

    r

    k

    calorde

    flujodel

    Variacin

    ,01,0,01,0,0,1

    2+

    +

    =

    +

    Balance de calor en el punto base del

    centro en el fondo

    El balance de calor en el punto del centro del

    rea del fondo tiene mucha similitud con el

    balance en la superficie del fondo, pero en este

    caso existen dos puntos vecinos simtricos que

    son representados por un solo punto. Para este

    punto cuyo esquema se muestra en la figura 7.

    ( ) ( )[ ]

    ( ) ( )pEXTZppZ

    pp

    rINT

    pp

    rEXT

    TThATTAZ

    k

    TTATTAr

    k

    calorde

    flujodel

    Variacin

    0,00,01,0

    0,00,10,00,1

    +

    +

    =

  • 7/24/2019 Balances Microscopicos de Calor(1)

    7/11TECNOLOGA QUMICAVol. XXVI, No. 2, 2006 59

    Comopp TT 0,10,1 = por simetra y las reas son iguales A rINT= ArEXT, la ecuacin se puede escribir de

    la forma:

    (9)( ) ( ) ++=

    pEXT

    ppZ

    pprINT TThTT

    ZkATTA

    rk

    calorde

    flujodel

    Variacin

    0,00,01,00,00,12

    Determinacin de las reas y elvolumen de control

    El rea de transferencia de calor del volumen

    de control en el caso de las superficies cilndricas,

    tiene caractersticas singulares.

    El rea seccional AZ es constante siempre,para el caso de la transferencia de calor cuando

    se utilizan espacios discretos igualmente espacia-

    dos, como es el ejemplo que tratamos.

    El rea lateral del cilindro Ar es funcin

    del radio y se incrementa en la medida

    que el radio crece segn la ecuacin:

    2422 rnZrAZ == , por lo que al

    efectuar el balance de calor debe elegirse el rea

    que mejor represente el flujo global de transferen-

    cia de calor, en el trabajo ( ) desarrollado para

    discutir este aspecto se demostr que para coor-

    denadas cilndricas el balance debe desarrollarse

    en el rea anular intermedia del incremento, esto

    es en el punto del ngulo2

    r.

    reas para el balance en los puntos interiores

    Para el caso de los puntos interiores el rea de

    transferencia se determina de acuerdo con el

    diagrama de la figura 2. El rea de transferencia

    se determina calculando el rea del volumen de

    control perpendicular al flujo de calor para todos

    los casos es:

    Teniendo en cuenta que rnr = y haciendo

    El volumen de control ser:

    382 rnZAVi Z ==

    reas para el balance en puntos bases de la superficie

    Para la superficie lateral segn el esquema mostrado en la figura 4:

    ( )[ ]2222;22

    2 rRRAyZRAZr

    RA ZrEXTrINT ==

    =

    Como los incrementos son iguales y rNR =

    El volumen de control ser: ( )1222 == NZAViZ

    +=

    +=

    =

    22

    222

    22;2

    22

    rr

    rrAyZ

    rrAZ

    rrA

    ZrEXTrINT

    ( ) ( ) 222 4122;122 rnAyrnArnA ZrEXTrINT =+==

    ( ) [ ] 222 124;12 rNAyrNArNA ZrEXTrINT ===

  • 7/24/2019 Balances Microscopicos de Calor(1)

    8/11TECNOLOGA QUMICAVol. XXVI, No. 2, 200660

    rea para el balance en puntos bases de la superficie del fondo

    Para los bordes, mostrado el ejemplo en la figura 5:

    las cuales se pueden escribir de forma discreta:

    El volumen de control ser:

    Para el resto de los puntos, las reas se dan en la tabla 1.

    ( ) ( )[ ]222;2

    2;

    2

    2 rrrrAZr

    rAZr

    rA ZrEXTrINT +=

    +=

    =

    [ ] ( ) 222 4;12;12 rnArnArnA ZrEXTrINT =+==

    Tabla 1reas de transferencia y volumen de los volmenes

    de control de los puntos bases en un cilindro 32 rVdyrAd ==

    Puntos bases ArINT/Ad ArINT/Ad AZ/Ad Vi/Vd

    T(n,m) 2(2n-1) 2(2n+1) 4n 8n

    T(0,m) 2 2 1 2

    T(N,m) 2(2N-1) 4N (2N-1) 2(2N-1)

    T(n,0) (2n-1) (2n+1) 4n 4n

    T(N,0) (2N-1) 2N 2N-1 2 N-1

    T(0,0) 1 1 1 1

    Obtencin de las ecuacionesdiscretas

    Ecuacin discreta para los puntos baseinteriores

    Obtenidas las reas de transferencia sesustituyen en la ecuacin (3) y este resultadose sust i tuye luego en la ecuacin (2) ,obtenindose:

    la cual se puede escribir de la forma:

    (10)( )

    ++

    +

    ++= +++ p

    qmn

    p

    mn

    p

    qmn

    p

    qmn

    p

    mn

    p

    qmn TTTnTnFoTFoT ,1,1,,,1,,1,1

    ,, 2

    1

    12

    1

    141

    donde Fo es el nmero adimensional de Fourier yviene dado por:

    (11)

    Ecuacin discreta para los puntos bases desuperficie

    Sustituyendo (2) en (4) y sustituyendo lasreas por su valor dado en la tabla 1, se tiene:

    2rCp

    tkFo

    =

    ( )

    ++

    +

    += ++

    + pmn

    p

    mn

    p

    mn

    p

    mn

    p

    mn

    p

    mqn

    p

    qmn TTTT

    n

    nT

    n

    nFoTT ,1,1,,1,1,

    1

    ,, 42

    12

    2

    12)

    34 rnZAVi

    Z ==

  • 7/24/2019 Balances Microscopicos de Calor(1)

    9/11TECNOLOGA QUMICAVol. XXVI, No. 2, 2006 61

    la cual se puede escribir de la forma:

    donde Bio es el nmero adimensional de Biot y viene dado por:

    (13)

    Balance de calor en puntos bases de la superficie del fondo

    Sustituyendo (5) en (2) y sustituyendo las reas por su valor dado en la tabla 1, se tiene:

    la cual se escribe definitivamente de la forma:

    (14)

    Procediendo de igual forma en todos los casos se obtiene:

    Ecuacin discreta para los puntos bases del borde

    (15)

    Ecuacin discreta para los puntos bases del eje central

    (16)

    Ecuacin discreta para el punto base del centro en el fondo

    (17)

    ( )

    +

    +++= ++

    + 112

    412

    42 1,01,1,,1,

    1,

    p

    mEXT

    p

    mN

    p

    mN

    p

    mN

    p

    mN

    p

    mN TN

    NBioFoT

    N

    NBioTTTFoTT

    +++++

    += +

    +EXT

    p

    qm

    p

    mN

    p

    mN

    p

    mN

    p

    mN

    p

    mN TN

    NBioTTTTFoT

    N

    NBioFoT

    122

    12141 ,,01,1,,1,

    1,

    k

    rhBio

    =

    ( ) ( )BioFoTBioTnn

    BioTTn

    n

    FoTT EXT

    p

    n

    p

    N

    p

    np

    n

    p

    n +

    ++++= ++ 22

    22

    122

    1211,0,1

    0,1

    0,

    +

    +++

    +

    +=

    +EXT

    p

    N

    p

    N

    p

    N

    p

    N TN

    NBioTTFoTN

    NBioFoT 112

    221`12

    2221 1,0,10,1

    0,

    ( ) ( )11,011,0,1,01,0 461 ++++ +++= pmpmpmpmpm TTTFoTFoT

    ( )( ) ( )EXTpppp BioTTTFoTBioFoT ++++=+ 1,00,10,010,0 2321

    Las seis ecuaciones discretas obtenidas paralos diferentes sistemas de puntos bases debernser aplicados en los puntos correspondientes a

    Variable Mnimo Mximo

    n 1 N-1

    .m 1 M

    sus lmites superior en inferior. En el caso expues-to de acuerdo con la simetra, todos los sistemasde puntos bases tendrn las siguientes lmites:

    ( )[ ]

    ++

    +

    +++= +

    +EXT

    p

    n

    p

    n

    p

    n

    p

    n

    p

    n TBioTTn

    nT

    n

    nFoTBioFoT 22

    2

    12

    2

    12221

    1,0,10,10,

    1

    0,

  • 7/24/2019 Balances Microscopicos de Calor(1)

    10/11TECNOLOGA QUMICAVol. XXVI, No. 2, 200662

    Las ecuaciones se evalan paso a paso par-tiendo del valor de temperatura inicial conocida, yactualizando las variables cada vez que se requie-ra un nuevo clculo.

    Convergencia y estabilidad de lasecuaciones discretas

    Las ecuaciones (10), (12), (14), (15), (16) y(17) son ecuaciones en diferencias finitas explci-tas, las cuales deben cumplir condiciones de con-vergencia y estabilidad.

    Todas las ecuaciones obtenidas son del tipo:

    y para que sean convergentes y estables los

    coeficientes a1deben cumplir las siguientes con-diciones:

    0.1 >ia Todos los coeficientes deben serpositivos.

    1.2

  • 7/24/2019 Balances Microscopicos de Calor(1)

    11/11TECNOLOGA QUMICA V l XXVI N 2 2006 63

    Bibliografa

    Barba, D., Calcolo Elettronico Nell IngegneriaChimica, Roma, Edicin Siderea, 1971.

    Bird, R. B., W. E. Steward, E. N. Lightfoot , Fenme-nos de Transporte, La Habana, Editorial Puebloy Educacin, 1979.

    Carnahan , B., Luthers, H. A., James, O. Wilker,pplied Numerical Methods, New York, Edito-rial John Wiley & Sons, 1969.

    Daz Garca, Armando, Mtodos numricos aplica-dos a la Ingeniera Qumica, Ediciones Internas,Facultad de Ingeniera Qumica, Universidad deOriente, 1998.

    ______, Soluc in de problemas de conduccin enms de una dimensin, en rgimen no estaciona-rio por mtodos de discretizacin, Revista Tec-nologa Qumica, Vol. XVII, No. 2, 1 997, pgs. 1-

    5._____, Apl icacin de mtodos de discre tizacin

    para la so lucin de problemas de conduccin decalor en coordenadas esfricas , Revista Tec-nologa Qumica, Vol. XVIII, Nos. 1 y 2, 1998,

    pgs . 75-77.

    _____, Apl icac in de mtodos discre tos para lasolucin de problemas de conduccin de caloren coordenadas cilndricas, Revista Tecnolo-ga Qumica, Vol. XVIII, No. 3, 1998, pgs. 45-48 .

    Daz, H., Meja, E., DazA., Estudio de la sol ucin de

    las ecuaciones diferenciales parciales aplicadasa Ingeniera Qumica, Trabajo de Diploma, Facul-tad de Ingeniera Qumica, Universidad de Orien-te, Santiago de Cuba, 2003.

    Incropera,B., Fundamentals of Heat and Transfer,La Habana, Editorial Pueblo y Educacin, 1973.

    Jenson, V. G., Jeffreys, G. V.,Mathemat ical Methodsin Chemical Engineering, La Habana, Institutodel libro, 1967.

    Lapidus, L., Digi ta l Computat ion for ChemicalEngineers, New York, Editorial McGraw Hill,1962.

    Mickley, H. S., Sherwood, T. K., Reed, Ch.,AppliedMa th ema ti cs in Ch em ic al En gi nee rin g , NewYork, Edicin Mc Graw Hill, 1957.

    Perry, J. H., Chemicals Engineers Handbook, NewYork, 5ta Edicin, Mc Graw Hill, 1973.

    Robsenow, W. H. , Choi , H. , Heat Mass andMomentum Transfer, New Jersey, Prentice Hall,1961.