|
[Sponsors] |
rhoSimpleFoam - Residuals do konverge, but oscillate all the way down |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
January 14, 2014, 05:30 |
rhoSimpleFoam - Residuals do konverge, but oscillate all the way down
|
#1 |
New Member
Join Date: Oct 2011
Posts: 19
Rep Power: 14 |
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 |
|
January 14, 2014, 09:57 |
|
#3 | |
New Member
Join Date: Oct 2011
Posts: 19
Rep Power: 14 |
Quote:
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; } } Last edited by okroud; January 16, 2014 at 05:37. |
||
January 28, 2014, 05:05 |
|
#4 |
New Member
Join Date: Oct 2011
Posts: 19
Rep Power: 14 |
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 |
|
Thread Tools | Search this Thread |
Display Modes | |
|
|
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 |