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

strange processor boundary behavior with linearUpwindV

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

Like Tree1Likes
  • 1 Post By ma-tri-x

Reply
 
LinkBack Thread Tools Display Modes
Old   November 4, 2011, 20:42
Default strange processor boundary behavior with linearUpwindV
  #1
New Member
 
Austin Kimbrell
Join Date: Feb 2011
Location: Tennessee, USA
Posts: 8
Rep Power: 6
akimbrell is on a distinguished road
Code used is a recent version of OpenFOAM-1.6-ext.

I observed the following behavior when running a test case in parallel - see attached image. The case is a point vortex with P = 0 boundaries and U boundaries set to zeroGradient. I am running standard icoFoam solver with 3 processors. For discretization schemes I am using the default Gauss linear everywhere except div(phi,U), for that one I use Gauss linearUpwindV Gauss linear. My max Courant number is ~0.3, grid resolution is reasonable.

The quantity shown is vorticity magnitude - you can see that when the scale is minimized, data is not being shared properly across the processor boundaries - they effectively show up in the solution when they should be invisible. I have run this case previously in serial and have not seen anything like this before. Furthermore, most other schemes for the convection term do not show this behavior. I have observed something similar also using SFCDV but it dissipates shortly after initialization.

Is this a known problem with the linearUpwindV scheme? I have also tried pure linearUpwind and did not see this problem at all.
Attached Images
File Type: png linearUpwindV.png (20.4 KB, 107 views)

Last edited by akimbrell; November 4, 2011 at 20:42. Reason: identified OpenFOAM version
akimbrell is offline   Reply With Quote

Old   November 5, 2011, 13:45
Default
  #2
Senior Member
 
Alberto Passalacqua
Join Date: Mar 2009
Location: Ames, Iowa, United States
Posts: 1,894
Rep Power: 26
alberto will become famous soon enoughalberto will become famous soon enough
This is typically due to a missing update of the BC's. If you can reproduce it in 2.0.x, please report it as a bug on http://www.openfoam.com/mantisbt/main_page.php .
__________________
Alberto Passalacqua

GeekoCFD - A free distribution based on openSUSE 64 bit with CFD tools, including OpenFOAM. Available as live DVD/USB, hard drive image and virtual image.
OpenQBMM - An open-source implementation of quadrature-based moment methods
alberto is offline   Reply With Quote

Old   February 11, 2013, 15:56
Default
  #3
Member
 
David Hora
Join Date: Mar 2009
Location: Zürich, Switzerland
Posts: 63
Rep Power: 8
david is on a distinguished road
http://www.openfoam.org/mantisbt/view.php?id=676

https://github.com/OpenFOAM/OpenFOAM...01782e30bd4291

Regards
David
david is offline   Reply With Quote

Old   November 6, 2013, 04:53
Default
  #4
s.m
Senior Member
 
saeideh mohamadi
Join Date: Aug 2012
Posts: 229
Rep Power: 5
s.m is on a distinguished road
Hi
would you please tell me that how linearUpwind can be define in OpenFOAM-1.6-ext?
e.g linearUpwind in openFoam 2.2.0 is defind in this way:
div(phi,U) Gauss linearUpwindV grad(U);

i want to know how should i define it in OpenFOAM-1.6-ext?

thank you very much
s.m is offline   Reply With Quote

Old   November 9, 2013, 14:16
Default
  #5
Super Moderator
 
Bruno Santos
Join Date: Mar 2009
Location: Lisbon, Portugal
Posts: 8,251
Blog Entries: 34
Rep Power: 84
wyldckat is just really nicewyldckat is just really nicewyldckat is just really nicewyldckat is just really nice
Greetings s.m,

You can find tutorials that use this scheme by running:
Code:
grep -R "linearUpwindV" $FOAM_TUTORIALS/
Best regards,
Bruno
wyldckat is offline   Reply With Quote

Old   May 18, 2014, 20:04
Default compressible InterFoam
  #6
New Member
 
Join Date: Sep 2013
Posts: 9
Rep Power: 3
ma-tri-x is on a distinguished road
Hi !

I worked several months on simualtion of cavitation bubbles. I started with the compressibleInterFoam solver and modified it, but the problem still also remains with the original one(s) (versions: 2.1.1, 2.2.1, 2.3.0):

When at least one of the processor patches crosses the interface, the velocity field is not computed correctly. You observe numerical fragments parallel to the patch like in the following picture
https://www.dropbox.com/s/a4cr4qpyoqjhue2/bug.png

The only version I discovered to be able to deal with the decomposition is the openfoam-extend version.

Has anyone observed this? Is this the same error as discussed here?

Furthermore I discovered that even without decomposing, running a 2D axis-symmetric case has some overestimating properties in the cell of the interface directly at the axis. So if you place a bubble that will collapse onto the axis, it will unphysically get thinner at the axis. This also seems to be no problem in the extend version 3.0 (absolutely same input files).

Regards,
Max
ngj likes this.

Last edited by ma-tri-x; May 18, 2014 at 20:08. Reason: forgot sth
ma-tri-x is offline   Reply With Quote

Old   May 20, 2014, 06:04
Default
  #7
Member
 
David Hora
Join Date: Mar 2009
Location: Zürich, Switzerland
Posts: 63
Rep Power: 8
david is on a distinguished road
Hi Max

Do you see the numerical fragments only with linearUpwindV or also with other schemes? I think it would be good to report this bug. Do you have a simple case that demonstrates your problems?

Regards
David
david is offline   Reply With Quote

Old   May 20, 2014, 12:13
Default I don't use linearUpwind
  #8
New Member
 
Join Date: Sep 2013
Posts: 9
Rep Power: 3
ma-tri-x is on a distinguished road
Hi David!

I don't use linearupwind. My fvSchemes looks like:
Code:
ddtSchemes
{
    default         Euler;
}

gradSchemes
{
    default         Gauss linear;
}

divSchemes
{
    default            Gauss vanLeer;
     div(phirb,alpha) Gauss interfaceCompression 1;
}

laplacianSchemes
{
   default              Gauss linear corrected;
}

interpolationSchemes
{
    default         linear;
}

snGradSchemes
{
    default              none;
//     default              Gauss skewCorrected linear;
    snGrad(pd)           limited 0.5;
    snGrad(rho)          limited 0.5;
    snGrad(alpha1)       limited 0.5;
    snGrad(p_rgh)      limited 0.8;
    snGrad(p)          limited 0.8;
}

fluxRequired
{
    default              none;
    p_rgh;
    p;
    pcorr;
    alpha1;
}
A simple example which is reproducable... erm ... I don't know and I'm currently writing my thesis so I have little time to prepare one...

In August, I will be able to prepare one...

I think you will immediately see the symmetry-axis-effect when you set up a bubble with very low pressure on an axis of an axis-symmetric mesh. in the very first timestep, where the bubble wall starts to accelerate, you find that the two cells where the interface hits the axis accelerate faster. It's hidden in the next timesteps because I think you cannot set a threshold for the U-field magnitude in paraview (the two cells are only slightly faster and the range of U-magnitude over the mesh is much larger). I cannot show you this unfortunately, because I didn't save a screenshot. But I can show you another freaky example:

Once I decomposed with scotch and some processor boundaries were parallel to the bubble wall at some time. When the bubble passed the processor patch, droplets were formed. Maybe I can put the screenshot here. A yes, here's an attachment. The left picture is at t=8.18393e-5, the right is at t=8.33299e-5. In between the bubble collapsed further and the interface (="bubble wall") passed the scotch-processor-patch.
Attached Images
File Type: jpg compareDecompositionAndSingleProc.jpg (15.7 KB, 24 views)

Last edited by ma-tri-x; May 20, 2014 at 12:17. Reason: better explanation
ma-tri-x 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
Strange Nut behaviour with K-OmegaSST nicolarre OpenFOAM Running, Solving & CFD 10 January 29, 2015 12:53
Strange Flamelet Generation Problem tar FLUENT 11 March 22, 2014 11:52
strange pressure behaviour with symmetricPlane boudary condition - interFoam duongquaphim OpenFOAM Running, Solving & CFD 10 August 20, 2013 14:00
Coordination Frame Problem - strange Luk_Fiz CFX 2 July 30, 2010 08:20
Strange things in SRF (urgent) Cem FLUENT 0 December 19, 2005 11:37


All times are GMT -4. The time now is 01:23.