Page 894 - Chapra y Canale. Metodos Numericos para Ingenieros 5edición_Neat
P. 894

870                     DIFERENCIAS FINITAS: ECUACIONES ELÍPTICAS

                                                                      100C


                                                               (1, 3)  (2, 3)  (3, 3)


                                                               (1, 2)  (2, 2)  (3, 2)
                                                       75C                           50C


                                                               (1, 1)  (2, 1)  (3, 1)


                                                                       0C

                                      FIGURA 29.4
                                      Una placa calentada donde las temperaturas frontera se mantienen a niveles constantes. Este
                                      caso se denomina condición de frontera de Dirichlet.



                                      En la malla cuadrada de la figura 29.3, ∆x = ∆y, y reagrupando términos, la ecuación se
                                      convierte en
                                         T   + T  + T   + T  – 4T  = 0                                 (29.8)
                                          i+1,j  i–1,j  i,j+1  i,j–1  i,j
                                      Esta relación, que se satisface por todos los puntos interiores de la placa, se conoce como
                                      ecuación laplaciana en diferencias.
                                         Además, se deben especificar las condiciones de frontera en los extremos de la
                                      placa para obtener una solución única. El caso más simple es aquel donde la temperatu-
                                      ra en la frontera es un valor fijo. Ésta se conoce como condición de frontera de Dirichlet.
                                      Tal es el caso de la figura 29.4, donde los extremos se mantienen a temperaturas cons-
                                      tantes. En el caso ilustrado en la figura 29.4, un balance en el nodo (1, 1) es, de acuerdo
                                      con la ecuación (29.8),
                                         T  + T  + T  + T  – 4T  = 0                                   (29.9)
                                          21   01  12   10   11
                                      Sin embargo, T  = 75 y T  = 0, y, por lo tanto, la ecuación (29.9) se expresa como
                                                  01      10
                                         –4T  + T  + T  = –75
                                            11   12  21
                                         Ecuaciones similares se pueden desarrollar para los otros puntos interiores. El resul-
                                      tado es el siguiente conjunto de nueve ecuaciones simultáneas con nueve incógnitas:

                                       4T     –T             –T                                      =  75
                                         11     21             12
                                       –T     +4T     –T            –T                               =   0
                                         11      21    31             22
                                              –T     +4T                   –T                        =   50
                                                21      31                   32
                                       –T                   +4T     –T            –T                 =  75
                                         11                    12     22            13
                                              –T             –T    +4T     –T            –T          =     0
                                                21             12     22     32            23
                                                      –T            –T     +4T                   –T  =   50
                                                       31             22      32                   33
                                                             –T                   +4T    –T          = 175
                                                               12                    13    23
                                                                    –T            –T     +4T     –T  = 100
                                                                      22            13      23     33
                                                                           –T            –T     +4T  = 150
                                                                             32            23      33
                                                                                                      (29.10)
                                                                                                         6/12/06   14:04:14
          Chapra-29.indd   870
          Chapra-29.indd   870                                                                           6/12/06   14:04:14
   889   890   891   892   893   894   895   896   897   898   899