| Transient geotherms |
||||||||||||||
| Numerical solutions |
||||||||||||||
| •The N equations above can be rearranged into the following matrix equations: | ||||||||||||||
 |
 |
|||||||||||||
| (7) | with | (8) | ||||||||||||
| Thus the RHS of (7) is determined explicitly from the solution at the n-th time step. Note that the boundary conditions are satisfied by ensuring that the coefficients given in (6) are satisfied. Hence in order to find the solution at the time step, one must solve N linear equations. Since the coefficient matrixes in (7) and (8) are tridiagonal, one can make use of efficient algorithms (e.g. Thomas algorithm) to find the solution. | ||||||||||||||
|
|
 |
|
|
|
 |
 |
 |
 |
 |
||