CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Running, Solving & CFD

Stripes fluctuations or clouds when using refinement boxes

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   November 8, 2017, 11:30
Default Stripes fluctuations or clouds when using refinement boxes
  #1
Member
 
hulli graemer
Join Date: Oct 2014
Posts: 48
Rep Power: 12
hulli is on a distinguished road
Dear Foamers,

I was wondering, if anybody could help me with that issue.

I am trying to simulate a flow past an obstacle (could be anything a cube an car etc.).
The model itself runs under pimpleFoam, and is a LES model (Wale).

For benchmarking reasons, and to keep the computational time relatively low, I am only focusing at the flow within the domain or the boxes in first a place. In general there are 3 refinement regions so far, ranging from level 0 – 2. I did that using snappyHexMesh (please see the overview figure). At the left side I have a simple velocity inlet with a uniform velocity fixed value (0.6 cm/s), and at the right side of the domain there is an inlet outlet boundary. Front, back, and bottom of the domain are walls/patches whereas the top is a slip boundary.

So here comes my problem:
I always observe strange “clouds / fluctuations or stripes” in the velocity flow field. As far I can tell these seem the result of nut, and become more pronounced in regions of bordering refinement boxes.

Why do I get these strange fluctuations or clouds of nut that affect the flow field?

Why is this mainly happening between regions of bordering refinement boxes and the inlet?

I attached some more figures of Ux and nut as Z slices cut through the domain

… some more background information:
The model uses only U,p and nut.
checkMesh – Mesh OK
The Courant number is below 1
The pressure seams also fine

Thank you so much.
Attached Images
File Type: png case_boxes_nut_surface.png (29.7 KB, 13 views)
File Type: png case_boxes_overview.png (70.8 KB, 13 views)
File Type: png case_boxes_U.png (129.2 KB, 18 views)
hulli is offline   Reply With Quote

Old   November 8, 2017, 11:43
Default
  #2
Member
 
Join Date: Jun 2017
Posts: 58
Rep Power: 9
sturgeon is on a distinguished road
Disclaimer: I am very new to OpenFOAM and more or less useless at it. With that in mind, have you tried using non-orthogonal correctors?

Like I said, I have very little experience but it's something someone recommended when I was having similar problems.

See the end of this tutorial:

https://www.openfoam.com/documentati...utorialse3.php
sturgeon is offline   Reply With Quote

Old   November 8, 2017, 11:57
Default
  #3
Member
 
hulli graemer
Join Date: Oct 2014
Posts: 48
Rep Power: 12
hulli is on a distinguished road
Dear sturgeon,

thank you so much for your suggestion!!!!
I just double checked that, and I have the following settings:
// Outer Loops (Pressure-Momentum Correction)
nOuterCorrectors 3;
// Inner Loops (Pressure Correction)
nCorrectors 2;
nNonOrthogonalCorrectors 1;

I also have the following relexation factors that should do things any better:

relaxationFactors
{
fields
{
p 0.4;
pFinal 1; //< 1 only if alot of nOuterCorrectors (~ 100) are done: Underrelaxation of the final iteration
nuTilda 0.4;
nuTildaFinal 1;
}

equations
{
U 0.6;
UFinal 1; //< 1 only if alot of nOuterCorrectors (~ 100) are done: Underrelaxation of the final iteration
}
}

so I will give that a try ... Thank you
hulli is offline   Reply With Quote

Old   November 8, 2017, 12:02
Default
  #4
Member
 
Join Date: Jun 2017
Posts: 58
Rep Power: 9
sturgeon is on a distinguished road
You are very welcome, I hope it works! Otherwise I am out of ideas, but I am sure someone else will be able to help.

Cheers
sturgeon
sturgeon is offline   Reply With Quote

Old   November 9, 2017, 04:59
Default
  #5
Member
 
Rodrigo
Join Date: Mar 2010
Posts: 98
Rep Power: 16
guin is on a distinguished road
check your input from fvSchemes dictionary too
guin is offline   Reply With Quote

Old   November 9, 2017, 05:15
Default
  #6
Member
 
hulli graemer
Join Date: Oct 2014
Posts: 48
Rep Power: 12
hulli is on a distinguished road
Dear guin,

thank you so much for your suggestion. I am doing that right now. Just in case I also uploaded my fvSolution as well as the fvScheme, if anyone is keen to double check that..

increasing the number of nNonOrthogonalCorrectors did not really help.

Again thank you so much
Best H
Attached Files
File Type: txt fvSchemes.txt (1.7 KB, 5 views)
File Type: txt fvSolution.txt (2.7 KB, 4 views)
hulli is offline   Reply With Quote

Reply

Tags
les, refinement regions, stripes

Thread Tools Search this Thread
Search this Thread:

Advanced Search
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 Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
[snappyHexMesh] SnappyHexMesh running killed! Mark JIN OpenFOAM Meshing & Mesh Conversion 7 June 14, 2022 02:37
[snappyHexMesh] Error snappyhexmesh - Multiple outside loops avinashjagdale OpenFOAM Meshing & Mesh Conversion 53 March 8, 2019 10:42
[snappyHexMesh] Edge Refinement fracasce OpenFOAM Meshing & Mesh Conversion 3 December 2, 2017 14:30
killed "snappyHexMesh" parkh32 OpenFOAM Pre-Processing 2 April 8, 2012 18:12
[snappyHexMesh] snappyHexMesh aborting Tobi OpenFOAM Meshing & Mesh Conversion 0 November 10, 2010 04:23


All times are GMT -4. The time now is 14:29.