APÉNDICE B: DERIVACIÓN DEL COEFICIENTE DE DIFUSIÓN NUMÉRICA
DEL MÉTODO MUSKINGUM-CUNGE
|
Figure B-1 Space-time discretization of kinematic wave equation.
|
Expandiendo la función de la malla Q( jΔx,nΔt ) (Fig. B-1) en serie de Taylor acerca del punto
(jΔx,nΔt) lleva a:
∂Q
1 ∂2Q
Q j n+1 =
Q j n + [ _____ ] j Δt + ___ [ ______ ] j Δt 2 + o (Δt 3)
∂t
2 ∂t 2
| (B.1) |
∂Q
1 ∂2Q
Q j+1n+1 =
Q j+1 n + [ _____ ] j+1 Δt + ___ [ ______ ] j+1 Δt 2 + o (Δt 3)
∂t
2 ∂t 2
| (B.2) |
∂Q
1 ∂2Q
Q j+1n =
Q j n + [ _____ ] n Δx + ___ [ ______ ] n Δx 2 + o (Δx 3)
∂x
2 ∂x 2
| (B.3) |
∂Q
1 ∂2Q
Q j+1n+1 =
Q j n+1 + [ _____ ] n+1 Δx + ___ [ ______ ] n+1 Δx 2 + o (Δx 3)
∂x
2 ∂x 2
| (B.4) |
Substituyendo las ecuaciones B.1 a B.4 en la ecuación 9-61 y despreciando los términos de tercer orden lleva a:
∂Q
1
∂2Q
X
{ [ ____ ] j
Δt +
____
[ ______ ] j Δt 2
}
∂t
2
∂t 2
| |
∂Q
1
∂2Q
+ (1 - X )
{ [ _____ ] j+1
Δt +
____
[ _____ ] j+1 Δt 2
}
∂t
2
∂t 2
| |
C
∂Q
1
∂2Q
+ ____
{ [ _____ ] n
Δx +
____
[ ______ ] n Δx 2
}
2
∂x
2
∂x 2
| |
C
∂Q
1
∂2Q
+ ____
{ [ _____ ] n+1
Δx +
____
[ ______ ] n+1 Δx 2
} = 0
2
∂x
2
∂x 2
| (B.5) |
en la cual C = c (Δt /Δx) es el número de Courant.
Expresando las derivadas en el punto de malla [( j + 1)Δx, (n + 1)Δt ]
en términos de las derivadas en el punto de malla ( jΔx, nΔt ) por medio de la serie de Taylor:
∂Q
∂Q
∂2Q
[ _____ ] j+1 =
[ _____ ] j + [ ______ ] j,n Δx + o (Δx 2)
∂t
∂t
∂x ∂t
| (B.6) |
∂Q
∂Q
∂2Q
[ _____ ] n+1 =
[ _____ ] n + [ _______ ] j,n Δt + o (Δt 2)
∂x
∂x
∂x ∂t
| (B.7) |
∂2Q
∂2Q
∂3Q
[ ______ ] j+1 =
[ ______ ] j + [ _______ ] j Δx + o (Δx
2)
∂t 2
∂t 2
∂t 2∂x
| (B.8) |
∂2Q
∂2Q
∂3Q
[ ______ ] n+1 =
[ ______ ] n + [ _______ ] n Δt + o (Δt 2)
∂x 2
∂x 2
∂x 2∂t
| (B.9) |
Substituyendo las ecuaciones B.6 a B.9 en la ecuación B.5 y despreciando los términos de tercer orden:
∂Q
1
∂2Q
X
{ [ _____ ] j
Δt +
____
[ ______ ] j Δt 2
}
∂t
2
∂t 2
| |
∂Q
∂2Q
1
∂2Q
+ (1 - X )
{ [ _____ ] j
Δt +
[ ______ ] j,n Δx Δt +
____
[ ______ ] j Δt 2
}
∂t
∂x ∂t
2
∂t 2
| |
C
∂Q
1
∂2Q
+ ____
{ [ _____ ] n
Δx +
____
[ ______ ] n Δx 2
}
2
∂x
2
∂x 2
| |
C
∂Q
∂2Q
1
∂2Q
+ ____
{ [ _____ ] n
Δx +
[ ______ ] j,n Δx Δt +
____
[ ______ ] n Δx 2
} = 0
2
∂x
∂x ∂t
2
∂x 2
| (B.10) |
En la ecuación B.10, dividiendo por Δt y simplificando:
∂Q
∂Q
Δt
∂2Q
c Δx
∂2Q
[ _____ ] j +
c [ _____ ] n + ____ [ ______ ] j + ______ [ ______ ] n
∂t
∂x
2
∂t 2
2
∂x 2
| |
C
∂2Q
+ Δx
{ ( 1 - X ) + ____ }
[ ______ ] j,n
= 0
2
∂x ∂t
| (B.11) |
Los dos primeros términos de la ecuación B.11 constituyen la ecuación de la onda cinemática, Ecuación 9-18.
Los términos remanentes son el error R del esquema númerico de primer orden:
Δt
∂2Q
c Δx
∂2Q
C
∂2Q
R = ____ [ ______ ] j + ______ [ ______ ] n + Δx
{ ( 1 - X ) + _____ }
[ ______ ] j,n
= 0
2
∂t 2
2
∂x 2
2
∂x ∂t
| (B.12) |
De la ecuación 9-18:
∂Q
∂Q
____ = - c ____
∂t
∂x
| (B.13) |
Por lo tanto:
∂2Q
∂2Q
______ = - c ______
∂x ∂t
∂x 2
| (B.14) |
∂2Q
∂2Q
______ = c 2 ______
∂t 2
∂x 2
| (B.15) |
Substituyendo las ecuaciones B.14 y B.15 en la ecuación B.12 y simplificando:
1 ∂2Q
R = c Δx ( X - ___ ) _____
2 ∂x 2
| (B.16) |
Comparando la ecuación B.16 con el lado derecho de la ecuación de la onda difusiva, que se repite aquí:
∂Q
∂Q
∂2Q
____ + c _____ = νh _______
∂t
∂x
∂x 2
| (B.17) |
se demuestra que el coeficiente de difusión numérica del método Muskingum-Cunge method es:
1
νh = c Δx ( ___ - X )
2
| (B.18) |
|
