_______________________________________________________________________________________

This code is meant to eventually reach a steady state but the requirements for this to happen will be held in the Reynolds number. That is, the flow into the box must not be too high as compared to the viscosity of the fluid inside the box. In the case that this is not true we will not get enough dissipation in the vorticity function and our solution will eventually diverge. To solve this one must only turn up the viscosity. However we must also note that the time scale at which we advance to the next solution is relative to the size of the box. If we are going to increase the scale at which we are looking to determine a solution we must also decrease the time increment, failure to do so will also result in divergence. This code will determine the stream function and the vorticity function but the often more useful result of the velocity vectors can be easily generated from equations 5. We can see from the next figure an appropriate choice of time step, size, viscosity, and inflow that allows the program to come close to its steady state. These figures show flow lines of constant velocity.