[Nek5000-users] boundary condition issue?

[nek5000-users at lists.mcs.anl.gov]

Hi nek users,

I am encountering a strange issue. Build on the example of cylinder flow, I
write a matlab script to generate the Cartesian mesh and a genbox script to
generate the mesh around the cylinder. When the domain is small, the code
runs. When I increases the domain size and cell numbers, the calculation
blows up after several steps?

See the picture attached here, there is a singular point on the boundary
for the pressure. But the velociy is ok, since I imposed a uniform flow for
inlet and top bottom.

It also sometimes shows another error, the problem size is too large, even
I run them in 32 procs instead of 4 as usual.

Here attached the output of nek.

Does someone have an idea?

Best regards

Tubois
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20140305/ef436b4e/attachment-0001.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: out
Type: application/octet-stream
Size: 41510 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20140305/ef436b4e/attachment-0001.obj>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: Capture du 2014-03-05 19:28:22.png
Type: image/png
Size: 20559 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20140305/ef436b4e/attachment-0001.png>