
[Sponsors] 
July 28, 2011, 04:53 
moving reference frame  weird results!

#1 
Senior Member

Hi experts,
I'm a newbie in OF and I have tried to change simpleFoam solver to achieve oscillatory reference frame, it seems that it's not a big deal but I faced weird results. I did these changes: Code:
scalar Pi = mathematicalConstant::pi; dimensionedVector omega ( "omega", dimensionSet(0, 0, 1, 0, 0, 0, 0), vector (0.0, 0.0, (amp*(2.0*Pi*freq)*Foam::cos((2.0*Pi*freq)*runTime.value()))) ); dimensionedVector alpha ( "alpha", dimensionSet(0, 0, 2, 0, 0, 0, 0), vector (0.0, 0.0, (amp*pow(2.0*Pi*freq,2)*Foam::sin((2.0*Pi*freq)*runTime.value()))) ); volVectorField source = (2.0*omega ^ U) + (omega ^ (omega ^ mesh.C())) + (alpha ^ mesh.C()); Code:
tmp<fvVectorMatrix> UEqn ( fvm::ddt(U) + fvm::div(phi, U)  fvm::laplacian(nu, U) == source ); Any suggestion or experience in this field is really appreciated. Amir Last edited by Amir; July 28, 2011 at 05:45. 

July 30, 2011, 04:21 

#3 
Senior Member
Alberto Passalacqua
Join Date: Mar 2009
Location: Ames, Iowa, United States
Posts: 1,912
Rep Power: 32 
Oscillatory system (=unsteady) in steady state solver (simpleFoam)?
__________________
Alberto Passalacqua GeekoCFD  A free distribution based on openSUSE 64 bit with CFD tools, including OpenFOAM. Available as in both physical and virtual formats (current status: http://albertopassalacqua.com/?p=1541) OpenQBMM  An opensource implementation of quadraturebased moment methods. To obtain more accurate answers, please specify the version of OpenFOAM you are using. 

July 30, 2011, 05:23 

#4  
Senior Member

Quote:
I have added the transient term to UEqn., fvm::ddt(U); I have used unsteady simpleFoam instead of icoFoam because it's more stable and I could use larger time steps. Do you mean this procedure is incorrect? Amir 

July 30, 2011, 05:34 

#5 
Senior Member
Alberto Passalacqua
Join Date: Mar 2009
Location: Ames, Iowa, United States
Posts: 1,912
Rep Power: 32 
Are you doing subinterations inside each time step? You can find it already done in pimpleFoam.
Best,
__________________
Alberto Passalacqua GeekoCFD  A free distribution based on openSUSE 64 bit with CFD tools, including OpenFOAM. Available as in both physical and virtual formats (current status: http://albertopassalacqua.com/?p=1541) OpenQBMM  An opensource implementation of quadraturebased moment methods. To obtain more accurate answers, please specify the version of OpenFOAM you are using. 

July 30, 2011, 05:44 

#6 
Senior Member


July 30, 2011, 16:32 

#7 
Senior Member
Alberto Passalacqua
Join Date: Mar 2009
Location: Ames, Iowa, United States
Posts: 1,912
Rep Power: 32 
Hi,
to answer your question you should know if your system reaches a steady state or not. Since it is intrinsically unsteady, it might never reach one. Anyway, to use SIMPLE for unsteady calculations you should:
Best,
__________________
Alberto Passalacqua GeekoCFD  A free distribution based on openSUSE 64 bit with CFD tools, including OpenFOAM. Available as in both physical and virtual formats (current status: http://albertopassalacqua.com/?p=1541) OpenQBMM  An opensource implementation of quadraturebased moment methods. To obtain more accurate answers, please specify the version of OpenFOAM you are using. 

August 3, 2011, 08:29 

#9 
Senior Member

Dear Alberto,
I did your recommandations, but now, I have a problem in convergency. I change both pimpleFoam and also pisoFoam but convergency issues exist. I couldn't reduce initail residuals in this case, consequently, there are some differences between results and analytic solution. Do you have any recommendation about another schemes or others? I used different grids and also schemes but I couldn't achieve better convergency. 

August 3, 2011, 11:17 

#10 
Senior Member
Alberto Passalacqua
Join Date: Mar 2009
Location: Ames, Iowa, United States
Posts: 1,912
Rep Power: 32 
The residual behaviour seems to show the solution is oscillating, which might indicate also a problem in the case setup.
However, you are using inconsistent numerical schemes:
Best,
__________________
Alberto Passalacqua GeekoCFD  A free distribution based on openSUSE 64 bit with CFD tools, including OpenFOAM. Available as in both physical and virtual formats (current status: http://albertopassalacqua.com/?p=1541) OpenQBMM  An opensource implementation of quadraturebased moment methods. To obtain more accurate answers, please specify the version of OpenFOAM you are using. 

August 3, 2011, 12:39 

#11  
Senior Member

Quote:
Although the results are better know but there are some differences between this results and analytic one. Anyway, I just want to ensure that whether the solution is converged or not; you said before that in unsteady cases, initial residuals should be around 0.001. Is this criteria acceptable in such cases (oscillatory nature)? (if your answer is No; how can we ensure that the solution is converged?) Thanks again 

August 3, 2011, 12:53 

#12 
Senior Member
Alberto Passalacqua
Join Date: Mar 2009
Location: Ames, Iowa, United States
Posts: 1,912
Rep Power: 32 
Yes, if you perform an unsteady calculation, the residuals should go down at each time step, until they reach the specified tolerance.
Residuals clearly won't do that in a steadystate solver, simply because your system does not have a steady state solution. Best,
__________________
Alberto Passalacqua GeekoCFD  A free distribution based on openSUSE 64 bit with CFD tools, including OpenFOAM. Available as in both physical and virtual formats (current status: http://albertopassalacqua.com/?p=1541) OpenQBMM  An opensource implementation of quadraturebased moment methods. To obtain more accurate answers, please specify the version of OpenFOAM you are using. 

August 3, 2011, 15:35 

#13  
Senior Member

Quote:


August 6, 2011, 17:33 

#15 
Senior Member
Alberto Passalacqua
Join Date: Mar 2009
Location: Ames, Iowa, United States
Posts: 1,912
Rep Power: 32 
Were you having problems with pimpleFoam? If yes, what kind of problems? There should be little difference compared to a "transient SIMPLE".
__________________
Alberto Passalacqua GeekoCFD  A free distribution based on openSUSE 64 bit with CFD tools, including OpenFOAM. Available as in both physical and virtual formats (current status: http://albertopassalacqua.com/?p=1541) OpenQBMM  An opensource implementation of quadraturebased moment methods. To obtain more accurate answers, please specify the version of OpenFOAM you are using. 

August 7, 2011, 04:15 

#16  
Senior Member

Quote:
As I said before, my problem in both pimpleFoam and icoFoam was reducing initial residuals of p and U, I tried different schemes and solution controls and also time steps but I didn't achieve acceptable result particularly peak values in analytic solution; on the other hand, by transientSimpleFoam and with the same time step and schemes I achieved nice results matched with analytic one and both residuals reduced to about 0.01; I know that PISO algorithm is more precise that SIMPLE one but maybe my management for this algorithm was not efficient! Can you explain the reason of that?
__________________
Amir 

August 7, 2011, 04:29 

#17 
Senior Member
Alberto Passalacqua
Join Date: Mar 2009
Location: Ames, Iowa, United States
Posts: 1,912
Rep Power: 32 
Are you using the same relaxation factors in pimpleFoam and transient SIMPLE?
If the setup is the same, they should give the same result.
__________________
Alberto Passalacqua GeekoCFD  A free distribution based on openSUSE 64 bit with CFD tools, including OpenFOAM. Available as in both physical and virtual formats (current status: http://albertopassalacqua.com/?p=1541) OpenQBMM  An opensource implementation of quadraturebased moment methods. To obtain more accurate answers, please specify the version of OpenFOAM you are using. 

August 7, 2011, 05:22 

#18  
Senior Member

Quote:
Thanks,
__________________
Amir 

Tags 
mrf, oscillatory 
Thread Tools  Search this Thread 
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
Enabling Acceleration in Moving Reference Frame  jpo  FLUENT  4  October 8, 2009 14:29 
Moving reference frame  sadanand  FLUENT  0  February 6, 2009 00:18 
G95 + CGNS  Bruno  Main CFD Forum  1  January 30, 2007 01:34 
TwoPhase Buoyant Flow Issue  Miguel Baritto  CFX  4  August 31, 2006 13:02 
Transforming results in rotating reference frame  Mark Render  Main CFD Forum  1  November 27, 2002 08:20 