The main program sets up geometry of solution domain for 2-D advective diffusion similar to SS version.
- time-stepping parameters nt and dt are specified here. Initial condition is also set here.
- The boundaries of this solution domain are on finite-volume faces, but additional nodes are added on the boundaries (Patankar's Practice B, p.69).
- Boundary conditions can be either fixed phi or fixed phi gradient at each face on the boundary. For example, in bc_e, first column is numerical value of BC, second column specifies BC type - 1 for value, 2 for gradient.
- u.m and w.m find velocities on appropriate volume boundaries
- A.m and F_upwind.m implement power-law advection scheme
Finite-Volume Formulator routine fvm_transient_adv_diff_2D.m now includes time step loop.
- Execution hits "pause" after each step. Hit any key to resume execution when ready.
- solve.m solves the resulting system of linear equations.