|March 13, 2007, 03:56||
Hi everybody, I'm trying to
Join Date: Mar 2009
Location: Belgrade, Serbia
Posts: 1Rep Power: 0
I'm trying to simulate axisymmetric flow in a pipe with contraction (ERCOFTAC database, Case 58). I've used simpleFOAM, and got very nice results with RNG k-epsilon. I've tried some other models offered in simpleFOAM, and with RSTM I had some problems. Here they are:
1) When I use RSTM I have to specify Reynolds stress tensor (R) at the inlet, am I right? Should it be isotropic, or not? What about the outlet - is it pressure outlet?
2) I've tried with isotropic and non-isotropic R at the inlet, but after about 500 time-steps (iterations) the solution singularity is written in my shell :-). I suppose it's because thr time-step continuity error is increasing during iterations. Any solution?
Thaks in advance,
|March 13, 2007, 05:49||
Heya Aleksandar, Doing axi-
Join Date: Mar 2009
Location: London, England
Posts: 1,785Rep Power: 22
Doing axi-symmetric flows with RSTM is known to be more demanding form the point of view of convergence and stability than the k-epsilon turbulence variants. I would advise running a turbulent viscosity model first, using the application called R which will calculate the Reynolds stress and starting from that as the initial condition. You may need to change the boundary conditions on R before restarting - I don't think they get picked up.
There are no limitations on what your R should look like at the inlet. It would be nice if it were in equilibrium with the flow + satisfying thr Schwartz inequality but the code will run anyway.
I would guess you will need to lower the under-relaxation factors to get the code behaving in a stable manner, especially if you are starting from uniform R field. You can control the continuity error by tightening the pressure solver tolerance - a combination of the above should be sufficient to get the things working.
Providing commercial FOAM/OpenFOAM and CFD Consulting: http://wikki.co.uk
|Thread||Thread Starter||Forum||Replies||Last Post|
|Burgerbs equation non constant Boundary Conditions Initial Conditions||arkangel||OpenFOAM Running, Solving & CFD||1||October 2, 2008 14:48|
|Integral boundary conditions turbulent intensitylength boundary conditions||olesen||OpenFOAM Running, Solving & CFD||0||July 27, 2006 07:18|
|RSTM convergence issues||nazario||Main CFD Forum||0||June 23, 2006 07:23|
|RSTM Implimentation||M. Choi||Main CFD Forum||2||June 15, 2005 11:36|
|use of the RSTM||Olfa Daboussi||Phoenics||1||May 19, 2005 12:42|