CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM

simulation of irrotational vortex fails

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Display Modes
Old   October 18, 2013, 05:10
Default simulation of irrotational vortex fails
  #1
Member
 
Matthias Walter
Join Date: Mar 2009
Location: Rostock, Germany
Posts: 61
Rep Power: 8
matthias is on a distinguished road
Hello folks,

I would like to address the following OF problem to the board since I got no satisfactory solution for it.

The computational domain is quadratic and consists of 257x257x1 cells. The boundary conditions for U are fixed at all four sides of the domain (bc. type fixedValues). These predefined values have been calculated using the analytical solution of the irrotational vortex equation. The internal velocity field is initialized using also the analytical solution. For pressure p, zero gradient boundary condition has been applied. Since all four "inlets/outlets" are fixed, the viscosity nu is set to zero (inviscid fluid). All numerical and interpolation schemes are second order (linear). The solver settings include PCG for p and PBiCG for U. Turbulence models are deactivated.

After some iterations of icoFoam or pimpleFoam or pisoFoam the solution of U shows unphysically oscillations. Especially the components of the gradient of U (i.e. to see in paraview) oscillate strongly after multiple iterations and therefore blow up the solution. Naturally, the initial analytical solution should not change in time since all sides are fixed (possibly there could be some very small numerical inaccuracies).

Higher order schemes, decreasing of time step, better interpolation schemes, decreasing of solver tolerances up to 10^-16 and application of analytical pressure gradient for p as fixed gradient bc didn't help to solve the problem.

Interestingly, the same case and setup works fine using a simple 2d finite differencing code (http://math.mit.edu/cse/codes/mit18086_navierstokes.pdf). Actually this OF case has no practical background but is an interesting test case for the numerics in OF.


Maybe someone has an idea how to handle this case using i.e. icoFoam?

The case is located at http://www.lemos.uni-rostock.de/uplo...lVortex.tar.gz

Best regards

Matthias
matthias is offline   Reply With Quote

Old   October 21, 2013, 02:19
Default
  #2
Senior Member
 
Jens Höpken
Join Date: Apr 2009
Location: Duisburg, Germany
Posts: 156
Rep Power: 8
jhoepken is on a distinguished road
Send a message via Skype™ to jhoepken
Hi Matthias,

I've tried your case, but it didn't run out of the box. I had to change the solver name for pFinal in system/fvSolution from PCG1 to PCG. An other thing is that I changed pRefValue in the PIMPLE subdict from 1e15 to 0 and now it is running with pimpleFoam, without any problems.

Hope this helps,
Jens
jhoepken is offline   Reply With Quote

Old   October 21, 2013, 03:05
Default
  #3
Member
 
Matthias Walter
Join Date: Mar 2009
Location: Rostock, Germany
Posts: 61
Rep Power: 8
matthias is on a distinguished road
Hi Jens,

how long did you run the case? I tried your settings and still after few iterations using pimpleFoam the derivatives of velocity (grad(U)) show strange oscillations.


Best regards

Matthias
matthias is offline   Reply With Quote

Old   October 21, 2013, 04:49
Default
  #4
Senior Member
 
Jens Höpken
Join Date: Apr 2009
Location: Duisburg, Germany
Posts: 156
Rep Power: 8
jhoepken is on a distinguished road
Send a message via Skype™ to jhoepken
I think I ran it for 20-40 iterations, without having any oscillations and I am using 2.2.x, but not the latest snapshot of the git repo.
jhoepken is offline   Reply With Quote

Reply

Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Differences and functions of Solidworks Simulation and Solidworks Flow Simulation? alpharays Main CFD Forum 0 April 19, 2012 03:13
Boundary condition setting regarding turbine simulation using CFX Lacerlacer CFX 11 March 12, 2012 10:32
Flow simulation of vortex street behind a cylinder Ramnik Main CFD Forum 4 June 24, 2009 09:29
problem with vortex shedding simulation Le Stanc Main CFD Forum 6 November 1, 2006 14:08
3-D Contaminant Dispersal Simulation Apple L S Chan Main CFD Forum 1 December 23, 1998 11:06


All times are GMT -4. The time now is 19:18.