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

rhoSimpleFoam - Residuals do konverge, but oscillate all the way down

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

Like Tree1Likes
  • 1 Post By okroud

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   January 14, 2014, 05:30
Default rhoSimpleFoam - Residuals do konverge, but oscillate all the way down
  #1
New Member
 
Join Date: Oct 2011
Posts: 19
Rep Power: 14
okroud is on a distinguished road
Hello guys,

I am calculating a simple laminar channel flow, with one inlet, one outlet, two walls, and periodic bc's in z-direction. Since i had stability issues, I'm starting with first order Schemes and a relaxation of 0.1 for p and 0.7 for U. Fortunetaly, the residuals do konverge (finally!). But they show oscillations as shown in the attached picture.



Could someone give me a hint about the way I could change my fvSchemes or fvSolution to dampen these oscillations. My fvSchemes and fvSolutions are attached in this Post.

Thanks in advance!

okroud
Attached Images
File Type: jpg Residuals_0_2000.jpg (71.4 KB, 95 views)
Attached Files
File Type: txt fvSolution.txt (1.8 KB, 30 views)
File Type: txt fvSchemes.txt (1.7 KB, 17 views)
okroud is offline   Reply With Quote

Old   January 14, 2014, 07:55
Default
  #2
Super Moderator
 
Tobi's Avatar
 
Tobias Holzmann
Join Date: Oct 2010
Location: Tussenhausen
Posts: 2,708
Blog Entries: 6
Rep Power: 51
Tobi has a spectacular aura aboutTobi has a spectacular aura aboutTobi has a spectacular aura about
Send a message via ICQ to Tobi Send a message via Skype™ to Tobi
Hi,

can you take a picture of your mesh and Initial conditions / velocity?
Additionally a slice of U at 1000, 1010 and 1020 iterations?
Tobi is offline   Reply With Quote

Old   January 14, 2014, 09:57
Default
  #3
New Member
 
Join Date: Oct 2011
Posts: 19
Rep Power: 14
okroud is on a distinguished road
Quote:
Originally Posted by Tobi View Post
can you take a picture of your mesh and Initial conditions / velocity?
Additionally a slice of U at 1000, 1010 and 1020 iterations?
Hi Tobi,
thank you for your interest. I will be able to post these things on Thursday.


Edit: Unfortunately i only have solutions of Timestep 1800,1900 an 2000. There is no visible difference between them. Indeed these timesteps could be interesting regarding the p-residual, I will reproduce them over the weekend. Pictures of the mesh are attached in this post.

Boundary conditions for U:

dimensions [0 1 -1 0 0 0 0];

internalField uniform (0 0 0);

boundaryField
{
in
{
type flowRateInletVelocity;
flowRate 0.17359;
value uniform (0 0 0); // placeholder
}
out
{
type zeroGradient;
}

lid
{
type fixedValue;
value uniform (0 0 0);
}
bottom
{
type fixedValue;
value uniform (0 0 0);
}

PER_MERGED //periodic bc's
{
type cyclic;
}
}

Boundary conditions for p:

dimensions [1 -1 -2 0 0 0 0];

internalField uniform 101325;

boundaryField
{
in
{
type zeroGradient;
}
out
{
type fixedValue;
value uniform 101325;
}

lid
{
type zeroGradient;
}
bottom
{
type zeroGradient;
}

PER_MERGED //periodic bc's
{
type cyclic;
}
}
Attached Images
File Type: jpg mesh_all.jpg (24.4 KB, 27 views)
File Type: jpg mesh_zoom1.jpg (60.5 KB, 22 views)
File Type: jpg mesh_zoom2.jpg (63.5 KB, 22 views)
File Type: jpg Velocity_Step2000.jpg (15.2 KB, 24 views)
File Type: jpg Velocity_Outlet_Step2000.jpg (15.0 KB, 23 views)
Tobi likes this.

Last edited by okroud; January 16, 2014 at 05:37.
okroud is offline   Reply With Quote

Old   January 28, 2014, 05:05
Default
  #4
New Member
 
Join Date: Oct 2011
Posts: 19
Rep Power: 14
okroud is on a distinguished road
I found the Solution to my Problem: The relaxation factor of the oscillating Solution for rho was 0.2. I deleted it completely (which i guess means that it is now 1). Now the residuals konverge without any oscillation. But this rises further questions:

1. Does anyone have a clue why this relaxation of rho is producing such great oscillations?

My guess: Since i use the thermoPhysical Model:

hPsiThermo<pureMixture<sutherlandTransport<specieT hermo<hConstThermo<perfectGas>>>>>;

rho is just calculated with the equation of state for perfect gas. So after rho is underrelaxated, the solution of rho is not anymore in harmony with the equations of momentum and energy, which are solved in the following iteration.

2. But why is the underrelaxation of rho implemented in the solver? What is it's use?

Would be great to hear some opinions
okroud is offline   Reply With Quote

Reply

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
under-relaxation factors -> level of residuals Zigainer FLUENT 19 July 21, 2017 17:53
motorBike Residuals for SST k-omega... and mine JR22 OpenFOAM Running, Solving & CFD 6 August 1, 2013 10:08
what to monitor besides residuals? franzdrs FLUENT 5 March 21, 2013 04:59
judging convergence through residuals MachZero Main CFD Forum 7 December 25, 2012 13:18
Convergence - scaled vs unscaled residuals HS FLUENT 1 November 7, 2005 06:45


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