Diferencias Finitas 1D PDF

6.
Diferencias finitas 1D
6. DIFERENCIAS FINITAS 1D
Existen dos métodos para generar un sistema de ecuaciones de diferencias

finitas para problemas de transferencia de calor:
• Comenzar con las ecuaciones diferenciales parciales del calor y expresarlas en

forma de diferencias finitas.
• Realizar un balance de energía para regiones específicas de interés.
En este proyecto se ha utilizado el primero de los métodos que a continuación

se va a detallar.
6.1. APROXIMACIÓN POR DIFERENCIAS FINITAS PARA DERIVADAS

La función T(x) puede ser aproximada por una serie de Taylor en x como:
dTx ∆x
d
Tx ∆x d Tx
Tx ∆x Tx ∆x · · ·
dx 2 dx
6 dx
dTx ∆x
d
Tx ∆x d Tx
Tx ∆x Tx ∆x · · ·
dx 2 dx
6 dx
Truncando la primera ecuación después de los dos primeros términos y

despejando la primera derivada de dicha ecuación obtenemos la aproximación en
diferencias finitas en forma implícita:
dTx Tx ∆x Tx

error∆x
dx ∆x
El error debido al truncamiento de la serie de Taylor es del orden de ∆x.
Adoración Cerezuela Parish Estudio térmico de las fachadas

Grupo de Termotecnia 55 ventiladas con aletas de PCM
6. Diferencias finitas 1D
Truncando la segunda ecuación después de los dos primeros términos y

despejando la primera derivada de dicha ecuación obtenemos la aproximación en
diferencias finitas en forma explícita:
dTx Tx Tx ∆x

error∆x
dx ∆x
Si truncando ambas ecuaciones en el tercer término, le restamos la segunda de las

ecuaciones a la primera y despejamos la primera derivada de la función, obtenemos la
aproximación de la primera derivada en forma centrada:
dTx Tx ∆x Tx ∆x

error∆x

dx 2 · ∆x
Para obtener la aproximación en diferencias finitas de la segunda derivada truncamos

las ecuaciones en el tercer término y las sumamos, para posteriormente despejar la
segunda derivada, obteniendo la siguiente relación, que es también una relación de
diferencias centrada:
d
Tx Tx ∆x Tx ∆x 2 · Tx
error∆x

dx
∆x
Tanto a la fase en estado sólido como a la fase en estado gaseoso se les puede
aplicar la ecuación de conducción de calor ya que por simplificación se suele asumir
que no existe convección en la fase líquida. Utilizaremos el subíndice a para
referirnos tanto al estado sólido, líquido o al cambio de fase, siendo sus propiedades
en cada caso distintas y que serán definidas en el próximo apartado. De manera que la
ecuación de conducción de calor que se utiliza es:
∂ ∂T ∂T
K · ρ · C ·
∂x ∂x ∂t

Aplicándole a esta ecuación el método de diferencias finitas explícitas

obtenemos:
Tx ∆x, t Tx ∆x, t 2 · Tx, t Tx, t ∆t Tx, t

K · Cp · ρ ·
∆x
∆t
Despejando el término de la temperatura en el tiempo (t+∆t) obtenemos la

ecuación que se utilizará para el cálculo del campo de temperaturas.
K ∆t
Tx, t ∆t ·
· Tx ∆x, t Tx ∆x, t 2 · Tx, t Tx, t
Cp · ρ ∆x
α · ∆t
M
∆x
K
α
ρ · Cp
6.2. ESTUDIO DE LAS HIPÓTESIS DE SIMULACIÓN
Debido a que es necesario establecer un rango en el cual ocurre el cambio de

fase, la primera variable en estudio es justamente este rango.
El segundo criterio a estudiar es el nivel de mallado necesario, dado que

aunque es bien conocido que cuanto más fina sea la malla, menor será el error,
consideraciones relacionadas con la acumulación de errores y el tiempo de cálculo y
los recursos computacionales, hacen recomendable que el afinamiento de la malla
tenga un límite.
Finalmente la tercera variable estudiada es el paso de tiempo, siendo válidas

las mismas consideraciones realizadas para el nivel de mallado.

6.2.1. Intervalo de temperaturas en el que ocurre el cambio

de fase
La siguiente figura muestra a nivel cualitativo la forma en la que se ha

realizado
ealizado la selección de rango de temperaturas en la que ocurre el cambio de fase.
En la figura resulta evidente el efecto de aumentar el rango de temperaturas.

Si se selecciona un rango estrecho de temperaturas para el cambio de fase, es
necesario considerar
rar una capacidad térmica muy grande. A medida que dicho rango
se amplia, la capacidad térmica equivalente se reduce.
Figura 38. Evolución de la capacidad térmica para distintos rangos de temperatura

de cambio de fase
El problema de seleccionar un rango

rango de temperaturas muy estrecho, es que
dificulta la convergencia del método numérico y aumenta el error que éste comete.
Por otro lado, un rango demasiado amplio, se aleja
leja de la situación real y finalmente el
resultado puede guardar poca relación con el caso
caso que se está estudiando.
A continuación se presentan los resultados obtenidos para diferentes rangos de

temperatura en los cuales se ha asumido que ocurre el cambio de fase y se compara
con la solución analítica exacta.

6.2.2. Nivel de mallado y paso de tiempo
Las dos últimas variables en estudio es el nivel de mallado y los pasos de

tiempo. La siguiente figura esquematiza lo que implica variar el paso de tiempo de
cálculo:
La ecuación del calor discretizada en diferencias finitas mediante el método

explícito solo converge si (Carnahan 1969)
0 $ % & 1(2
Esta restricción significa si se quiere que el error de truncar la serie de Taylor

sea reducido debe reducirse ∆x, por lo que al deber cumplirse la restricción el paso de
tiempo también debe reducirse:
∆x
∆t &
2·α
Esta puede ser una restricción severa en términos de tiempo computacional,

ya que requiere muchas iteraciones para alcanzar un tiempo determinado.
Sin embargo, este método presenta una gran ventaja que, como su propio
nombre indica, las ecuaciones de cada nodo se resuelven una a una y por estricto
orden.
Las limitaciones del método explícito se resuelven considerando un diferente

esquema de paso de tiempo, es decir, utilizando el método implícito o el centrado.
Además, el método centrado tiene la ventaja adicional de tener un menor error de
truncamiento que el método explícito y el implícito. Sin embargo, ambos métodos
requieren la resolución de un sistema de ecuaciones algebraicas simultáneo.
A partir de simulaciones y análisis, se encontró que se puede definir un

número adimensional que permite definir un paso de tiempo adecuado para un
problema particular.
∆t
N α · St ·
∆x

En las pruebas realizadas se ha encontrado que la convergencia y el grado de

error son buenos si N & 1.
Cuanto menores sean ∆x y ∆t el resultado será más aproximado a la solución

exacta.
6.3. CAMBIO DE FASE EN DIFERENCIAS FINITAS

El área encerrada por la curva del calor específico debe ser igual que el efecto
del calor latente en el incremento de temperatura ∆T, pero Bonacina et al.(1974) notó
que la forma que asumía la curva no parecía tener importancia para el cálculo del
campo de temperaturas. Bonacina et al. (1973) sugirió que el efecto del calor latente
puede asumirse que ocurre en un rango de temperaturas Tf ± ∆T. Entonces el calor
específico aparente puede escribirse como:
12 123∆1
2 · C · ∆T - Ci · dT L - Cs · dT
124∆1 12
Si los calores específicos son constantes con la temperatura, el calor

específico aparente es:
Ci si T $ :; ∆:
9
7 Ci Cs L
Ca si Tf ∆T $ : $ :; ∆: >
8 2 2 · ∆T
7
6 Cs si T = :; ∆:
Bonacina et al. (1973) también propuso que la conductividad térmica podía

ser aproximada como:
Ki si T $ :; ∆:
9
7 Ks Ki
Ka Ki · T Tf ∆T si Tf ∆T $ : $ :; ∆: >
8 2 · ∆T
7
6 Ks si T = :; ∆:

6.4. CONDICIONES DE TEMPERATURA IMPUESTA

Para comprobar la validez de las soluciones numéricas propuestas las vamos a
comparar con las soluciones analíticas. También comprobaremos cómo varían las
respuestas variando el intervalo de temperatura en que ocurre el cambio de fase en la
solución numérica modelada. Se parte de las propiedades del material cambio de fase
con 23ºC de temperatura de fusión PLUSICE® , sacadas del catálogo de la página
web de dicha empresa, siendo sus propiedades las siguientes:
- Densidad: 785 Kg/m3

- Calor latente: 170 KJ/Kg
- Conductividad térmica: 0.18 W/(m⋅K)
- Calor específico: 2.22 KJ/(Kg⋅K)
- Espesor del bloque: 0.5 m
- Tiempo final: 100000 s
- Temperatura superficial constante a lo largo del tiempo: -100ºC
- Intervalo en que se produce el cambio de fase: ∆T = 0.02ºC
- Incremento de tiempo: ∆t = 80.71 s
- Incremento de espacio: ∆x = 0.005 m
- Número de nodos espaciales: N=100
- Número de nodos temporales: M=1239
La temperatura inicial del bloque para la solución numérica será de23.01ºC,

pero para la solución analítica será se 23ºC.Esto se debe a que si la temperatura
inicial de la solución numérica coincidiera con 23ºC, en el cambio de fase no se
utilizaría todo el calor disponible en el material.
El número de nodos temporales se calcula de modo que el número de Fourier

sea 3 para de esta forma asegurar la convergencia.
Se van a calcular varios tipos de errores, el error máximo relativo en

temperaturas, el error medio temporal máximo relativo en temperaturas, el error
medio espacial máximo relativo en temperaturas, el error medio total relativo en
temperaturas y el error máximo relativo en calor transferido.

EVOLUCIÓN A LO LARGO DEL EJE
20
-20
anal_5000
anal_50000s
anal_100000s
-40 df_5000s
df_50000s
df_100000s
-60
-80
-100
0 0.02 0.04 0.06 0.08 0.1 0.12 0.14 0.16 0.18 0.2
Figura 39. Evolución de la temperatura a lo largo del eje para distintos instantes de
tiempo para condición de temperatura impuesta
Se calcula el error relativo en temperatura para cada instante de tiempo y para

cada nodo, de manera que tenemos 100x1239=123900 errores relativos:
T@AíCDEF TGD2.HD@DCI
error1 t, x abs
T@AíCDEF
De todos estos errores el valor máximo es de 5.73%.
Calculamos a continuación la media de los errores en cada instante de tiempo

para todos los nodos, de manera que obtendremos 1239 errores relativos medios en
temperatura para cada instante de tiempo:
∑L error1 t, x
error
JJJJJJJJt
1
N
De estos errores el valor máximo es de 0.073%.

Calculamos a continuación la media de los errores en cada nodo para todos los
instantes de tiempo, de manera que obtendremos 100 errores relativos medios en
temperatura para cada nodo:
∑C error1 t, x
error1
JJJJJJJJx
M
De estos errores el valor máximo es de 0.16%.
Calculamos a continuación la media de todos los errores para todos los nodos
y para todos los instantes de tiempo, de manera que obtendremos un error relativo
medio en temperatura:
∑C ∑L error1 t, x
error1
MMMMMMMM
M·N
El error relativo medio es de 0.03%
Esta diferencia en la magnitud de los errores en temperatura se debe a que se

producen errores mayores en los primeros instantes de tiempo y en los primeros
nodos espaciales, sin embargo, a medida que se avanza en el tiempo y en los nodos
los errores se reducen drásticamente.
Para calcular el error relativo del calor transferido, primero calcularemos el

calor transferido en cada instante de tiempo tomando la diferencia de temperaturas
entre los dos primeros nodos más cercanos a la superficie para a continuación sumar
todos estos calores obteniendo el calor transferido total.
∆Tt T
t TO t
qt K · · ∆t K · · ∆t
∆x ∆x
q P qt
C
El error relativo en calor transferido queda por tanto:
∑C q @AíCDEF t ∑C q GD2.HD@DCI t

errorQ 0.78%
∑C q @AíCDEF t

Con las mismas especificaciones y propiedades vamos a ver cómo influye el

intervalo de temperaturas en el que se efectúa el cambio de fase. Vamos a estudiar 3
casos:
• ∆T=0.2 K
• ∆T=0.02 K
• ∆T=0.002 K
EVOLUCIÓN A LO LARGO DEL EJE X
23.1
23.05
23
AT=0.002ºC
22.95 AT=0.02ºC
AT=0.2ºC
22.9
T em peratura [ºC]
solución analítica
22.85
22.8
22.75
22.7
22.65
22.6
0.1346 0.1348 0.135 0.1352 0.1354 0.1356 0.1358 0.136 0.1362 0.1364
x[m]
Figura 40. Evolución de la temperatura a lo largo del eje x para distintos ∆T
Aunque la temperatura inicial en cada caso es de 23.001ºC, 23.01ºC y 23,1ºC

se observa claramente como al aumentar ∆T la curva de la solución numérica se aleja
de la solución analítica, cuya temperatura inicial se ha cogido de 23ºC.

6.5. CONDICIONES CONVECTIVAS

Para problemas transitorios 1-D con propiedades constantes ρ, CU y K sujetos
a condiciones de contorno convectivas, para un pequeño intervalo de tiempo ∆t el
método del balance de energía da la siguiente ecuación para el nodo del contorno
(D.R.Croft y D.G.Lilley):
∆TX K
ρ · CU · VX · · TO TX h · TZ TX
∆t ∆x
donde VX ∆x(2.
0 1 2 3
∆x ∆x ∆x
Figura 41. Esquema de los nodos más superficiales de material
Expresando la ecuación anterior de forma explícita obtenemos:
1
TX t ∆t 2 · Fo · \TO t Bi · TZ TX t · 1 Bi^
2 · Fo
Donde
K h · ∆x
Fo Bi
ρ · CU · ∆x
K

La comprobación de la validez de las soluciones numéricas propuestas, se

hace mediante la comparación con las soluciones analíticas con condiciones
convectivas; de forma análoga a lo realizado con las soluciones en condiciones de
temperatura impuesta.
El material y las condiciones seleccionadas para la simulación son las mismas

que para el caso de condiciones de temperatura impuesta, únicamente que la
temperatura superficial no se sabe a priori, siendo la temperatura conocida la
ambiental y teniendo que añadir el valor del coeficiente de película:
- Temperatura ambiental constante a lo largo del tiempo: 278.15 K

- Coeficiente de película: h=12 W/(m2⋅K)
Figura 42. Evolución de la temperatura superficial para condiciones convectivas

para 300 y 350 nodos

Figura 43. Evolución de la temperatura a lo largo del eje con condición convectiva
para distintos instantes de tiempo
La evolución de la temperatura a lo largo del bloque se ha simulado con

únicamente 100 nodos, que aunque se ha visto que se comete mayor error en la
temperatura superficial, el error cometido a lo largo del bloque es bastante pequeño.
Vemos como la exactitud se incrementa a medida que se avanza en el tiempo.
Aunque el bloque mide medio metro, para los tiempos simulados, el resto del bloque
se queda en la temperatura inicial, es decir en la temperatura de cambio de fase.
Analizando los resultados obtenidos por el programa diseñado, se observa

como la evolución de la temperatura superficial con el tiempo coincide con la
solución analítica. Al haberse obtenido la temperatura superficial de un balance de
energía discretizando el bloque, a medida que aumentamos el número de nodos se
observa como la solución obtenida por diferencias finitas se acerca más a la analítica
y es más continua.
Vamos a calcular los mismos tipos de errores que se calcularon en la solución

de condiciones de temperatura impuesta: el error máximo relativo en temperaturas, el
error medio temporal máximo relativo en temperaturas, el error medio espacial

máximo relativo en temperaturas, el error medio total relativo en temperaturas y el

error máximo relativo en calor transferido.
El error relativo en temperatura para cada instante de tiempo y para cada nodo
es:
T@AíCDEF TGD2.HD@DCI
error1 t, x abs
T@AíCDEF
El valor máximo es de 0.82%.
El error relativo medio en temperatura para cada instante de tiempo es:
∑L error1 t, x
error
JJJJJJJJt
1
N
El error relativo medio en temperatura para cada nodo es:
∑C error1 t, x
error1
JJJJJJJJx
M
El error relativo medio en temperatura es:
∑C ∑L error1 t, x
error1
MMMMMMMM
M·N
El error relativo medio es de 0.008%
El error relativo en calor transferido queda es:
∑C q @AíCDEF t ∑C q GD2.HD@DCI t

errorQ 0.26%
∑C q @AíCDEF t
Todos los errores son menores del 1%, por lo que es aceptable considerar que
COMSOL es un buen método para calibrar el método de diferencias finitas para
condiciones convectivas.

6.6. VALIDACIÓN CON COMSOL EN CONDICIONES DISTINTAS A

LAS ANALÍTICAS
A continuación utilizando COMSOL como programa numérico de referencia

vamos a comparar una solución numérica con nuestro modelo propio con temperatura
impuesta y una con condición convectiva con condiciones distintas a las analíticas. Se
realizarán las simulaciones con el mismo material A_23.
6.6.1. Condición de temperatura impuesta
La temperatura superficial impuesta será de -100ºC, mientras que la

temperatura inicial será de 50ºC.El resto de especificaciones y condiciones serán las
mismas que en los apartados anteriores.

60
40
20
0
COMSOL_5000s
Temperatura[ºC]
COMSOL_20000s
COMSOL_50000s
-20
COMSOL_100000s
DF_5000s
-40 DF_20000s
DF_50000s
DF_100000s
-60
-80
-100
0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5
x[m]
Figura 44. Evolución a lo largo del eje para condiciones de temperatura impuesta
distintas a las analíticas para distintos instantes de tiempo

Se calculará el máximo error en los cuatro tiempos simulados. El error en

cada uno de los instantes de tiempo se va a calcular cómo:
T_àb`c Tde
error abs
T_àb`c
El error máximo en cada instante de tiempo es:
• t=5000s → 1.3%
• t=20000s→1.11%
• t=50000s→1.35%
• t=100000s→0.76%
El error es un poco mayor que la comparación con las soluciones analíticas,

sin embargo, sigue siendo muy bajo, de modo que la evolución seguida por el método
de diferencias finitas tiene una gran semejanza con los resultados obtenidos con
COMSOL.
6.6.2. Condición convectiva
La temperatura ambiente será de 5ºC, mientras que la temperatura inicial será

de 50ºC.El resto de especificaciones y condiciones serán las mismas que en los
apartados anteriores.


55
50
45
40
35 COMSOL_5000s
Temperatura[ºC]
COMSOL_20000s
COMSOL_50000s
30
COMSOL_100000s
DF_5000s
25
DF_20000s
DF_50000S
20
DF_100000s
15
10
5
0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5
x[m]
Figura 45. Evolución a lo largo del eje para condiciones de temperatura impuesta
distintas a las analíticas para distintos instantes de tiempo
Al observar la evolución en varios instantes de tiempo (5000s, 20000s,

50000s, 100000s) se ve como la temperatura superficial influye de gran manera en
dicha evolución. Se aprecia cómo en la temperatura del cambio de fase, en el
programa COMSOL existe una gran diferencia de pendientes, que el programa en
diferencias finitas no es tan pronunciado. Esto se debe a que en COMSOL el
intervalo de temperaturas en que se produce el cambio de fase es 1ºC, mientras que
en programa de diferencias finitas es de 0.02ºC.
En la siguiente figura se representa la evolución de la temperatura superficial

con respecto al tiempo. Se destaca cómo las curvas se distancian a partir de que se
alcanza la temperatura de cambio de fase, a partir de la cual la evolución es fluctuante
hasta que se consigue una evolución más suave donde ambas curvas tienden a volver
a acercarse.

Figura 46. Evolución de la temperatura superficial a lo largo del tiempo para

condiciones convectivas distintas a las analíticas
Se calculará el máximo error en los cuatro tiempos simulados. El error en

cada uno de los instantes de tiempo se va a calcular cómo:
T_àb`c Tde
error abs
T_àb`c
El error máximo en cada instante de tiempo es:
• t=5000s → 1.07%
• t=20000s→0.43%
• t=50000s→1.14%
• t=100000s→0.87%
El máximo error en la evolución de la temperatura superficial a lo largo del

tiempo es 1.3%.Es lógico que el error sea mayor en la superficie, ya que como la
solución tiende a la temperatura inicial, el error tiende a disminuir.


Diferencias Finitas 1D PDF

Cargado por

Copyright:

Formatos disponibles

Diferencias Finitas 1D PDF

Cargado por

Información del documento

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Diferencias Finitas 1D PDF

Cargado por

Copyright:

Formatos disponibles

6.

Existen dos métodos para generar un sistema de ecuaciones de diferencias

• Comenzar con las ecuaciones diferenciales parciales del calor y expresarlas en

• Realizar un balance de energía para regiones específicas de interés.

En este proyecto se ha utilizado el primero de los métodos que a continuación

6.1. APROXIMACIÓN POR DIFERENCIAS FINITAS PARA DERIVADAS

Truncando la primera ecuación después de los dos primeros términos y

dTx Tx ∆x  Tx

El error debido al truncamiento de la serie de Taylor es del orden de ∆x.

Adoración Cerezuela Parish Estudio térmico de las fachadas

Truncando la segunda ecuación después de los dos primeros términos y

dTx Tx  Tx  ∆x

Si truncando ambas ecuaciones en el tercer término, le restamos la segunda de las

dTx Tx ∆x  Tx  ∆x

Para obtener la aproximación en diferencias finitas de la segunda derivada truncamos

Adoración Cerezuela Parish Estudio térmico de las fachadas

Aplicándole a esta ecuación el método de diferencias finitas explícitas

Tx ∆x, t Tx  ∆x, t  2 · Tx, t Tx, t ∆t  Tx, t

Despejando el término de la temperatura en el tiempo (t+∆t) obtenemos la

6.2. ESTUDIO DE LAS HIPÓTESIS DE SIMULACIÓN

Debido a que es necesario establecer un rango en el cual ocurre el cambio de

El segundo criterio a estudiar es el nivel de mallado necesario, dado que

Finalmente la tercera variable estudiada es el paso de tiempo, siendo válidas

Adoración Cerezuela Parish Estudio térmico de las fachadas

6.2.1. Intervalo de temperaturas en el que ocurre el cambio

La siguiente figura muestra a nivel cualitativo la forma en la que se ha

En la figura resulta evidente el efecto de aumentar el rango de temperaturas.

Figura 38. Evolución de la capacidad térmica para distintos rangos de temperatura

El problema de seleccionar un rango

A continuación se presentan los resultados obtenidos para diferentes rangos de

Adoración Cerezuela Parish Estudio térmico de las fachadas

6.2.2. Nivel de mallado y paso de tiempo

Las dos últimas variables en estudio es el nivel de mallado y los pasos de

La ecuación del calor discretizada en diferencias finitas mediante el método

Esta restricción significa si se quiere que el error de truncar la serie de Taylor

Esta puede ser una restricción severa en términos de tiempo computacional,

Las limitaciones del método explícito se resuelven considerando un diferente

A partir de simulaciones y análisis, se encontró que se puede definir un

Adoración Cerezuela Parish Estudio térmico de las fachadas

En las pruebas realizadas se ha encontrado que la convergencia y el grado de

Cuanto menores sean ∆x y ∆t el resultado será más aproximado a la solución

6.3. CAMBIO DE FASE EN DIFERENCIAS FINITAS

Si los calores específicos son constantes con la temperatura, el calor

Bonacina et al. (1973) también propuso que la conductividad térmica podía

Adoración Cerezuela Parish Estudio térmico de las fachadas

6.4. CONDICIONES DE TEMPERATURA IMPUESTA

- Densidad: 785 Kg/m3

La temperatura inicial del bloque para la solución numérica será de23.01ºC,

El número de nodos temporales se calcula de modo que el número de Fourier

Se van a calcular varios tipos de errores, el error máximo relativo en

Adoración Cerezuela Parish Estudio térmico de las fachadas

EVOLUCIÓN A LO LARGO DEL EJE

Se calcula el error relativo en temperatura para cada instante de tiempo y para

De todos estos errores el valor máximo es de 5.73%.

Calculamos a continuación la media de los errores en cada instante de tiempo

De estos errores el valor máximo es de 0.073%.

Adoración Cerezuela Parish Estudio térmico de las fachadas

De estos errores el valor máximo es de 0.16%.

El error relativo medio es de 0.03%

dTx Tx ∆x Tx

dTx Tx Tx ∆x

dTx Tx ∆x Tx ∆x

Tx ∆x, t Tx ∆x, t 2 · Tx, t Tx, t ∆t Tx, t

∑C q @AíCDEF t ∑C q GD2.HD@DCI t

∑C q @AíCDEF t ∑C q GD2.HD@DCI t