jinwon park

January 26, 2008 12:10 
Set a set of circular initial conditions on a mesh
I am solving 2D compressible flows. As an initial condition, a highpressure circle is embedded on a cartesian mesh. Currently, I am using the way like sqrt((xxc)^2+(yyc)^2). However, it produced a staircase irregular something which is not a real circle. It may cause instability at later time and lead to the sudden abortion of codes. I am finding the way to get a real circle as possible. Is there a good way to set circular initial condition on a cartesian mesh? I really need your help all. Thanks in advance.
