Over the past number of weeks, I have been learning Finite Element Analysis to use for solving the Navier-Stokes equations of fluid dynamics. Finite Element Analysis is used to solve partial differential equations over an unstructured grid. Usually the grid is made up of triangles or quadrilaterals. In my case, I used triangles since any mesh can be constructed from triangles.
Here are some results. The following is a simulation of Vortex Shedding of an incompressible fluid due to a circular obstacle. The simulation was written in Python using Numpy, Scipy, and Matplotlib and employs the pressure correction method.
Properties of the fluid:
Boundary Velocities:
![Rendered by QuickLaTeX.com \[U(0,y) = U(5,y) =U(x,0) = U(x,3) =1.0\]](http://www.rabidgeek.net/wp-content/ql-cache/quicklatex.com-0b3b6d01d2143708f96e512fa6d09b0e_l3.png)
![Rendered by QuickLaTeX.com \[Re = 500\]](http://www.rabidgeek.net/wp-content/ql-cache/quicklatex.com-761ee1efc4bfe82f15129c0926db21e4_l3.png)
Another video with the parameters slightly changed. This uses No-Slip boundary conditions on the top and bottom boundaries as well as a smaller fluid channel
Boundary Velocities:
![Rendered by QuickLaTeX.com \[U(0,y) = U(5,y) = 1.0\]](http://www.rabidgeek.net/wp-content/ql-cache/quicklatex.com-45f07c7de9a0f62cbd069c9d3d1c332c_l3.png)
![Rendered by QuickLaTeX.com \[U(x,0) = U(x,1) = 0.0\]](http://www.rabidgeek.net/wp-content/ql-cache/quicklatex.com-18156d179697267df5c6070f3258cc77_l3.png)
![Rendered by QuickLaTeX.com \[Re = 1000\]](http://www.rabidgeek.net/wp-content/ql-cache/quicklatex.com-cebac4c687b3f21c35549d05e5f027d9_l3.png)
I will release the Python code after I have cleaned it up a bit.


{execi 1200 python /home/commander/scripts/tw.py /home/commander/scripts/tweets.txt}
