Changes between Version 52 and Version 53 of u/erica/GudonovMethodEuler
- Timestamp:
- 05/07/13 10:52:32 (12 years ago)
Legend:
- Unmodified
- Added
- Removed
- Modified
-
u/erica/GudonovMethodEuler
v52 v53 117 117 i. Call Bcond -- Boundary condition SR that sets left and right ghost zones to adjacent physical cells (recall boundary conditions are set to transmissive). 118 118 ii. Call CFLcond -- Computes timestep using cfl coeff and max speed, where the max speed is also computed in this SR. 119 iii. Call GetFluxes -- This SR sets up the local Riemann Problem (LRP) and solves it exactly to get rho(x/t), u(x/t), P(x/t) at x/t = 0 (see above explanation of sampling routine). This is the most complicated structure of the program. 119 iii. Call GetFluxes -- This SR sets up the local Riemann Problem (LRP) and solves it exactly to get rho(x/t), u(x/t), P(x/t) at x/t = 0 (see above explanation of sampling routine). This is the most complicated structure of the program. A schematic of the loop is as follows: 120 120 121 121 [[Image(LoopSR.png, 35%)]] 122 122 123 i v. Call Update -- updates cells using computed flux124 v. Call Output -- output new values of cells125 vi. Begin loop again for n+1 step.123 i.i. Call Update -- updates cells using computed flux 124 ii.i. Call Output -- output new values of cells 125 iii.i. Begin loop again for n+1 step. 126 126 ''' 127 127