
[Sponsors] 
November 26, 2007, 06:05 
Description:
Although LRR is

#1 
Senior Member
Dragos
Join Date: Mar 2009
Posts: 648
Rep Power: 12 
Description:
Although LRR is a valid choice in const/turbulenceProperties, the case doesn't run. Solver/Application: rasInterFoam Source file: (if possible name of OpenFOAM library and file in the library) Testcase: Case : damBreak Platform: SLED 10 (SuSE) amd64 Version: 1.4.1 Notes: Nprocs : 1 Create time Create mesh for time = 0 Reading environmentalProperties Reading field pd Reading field gamma Reading field U Reading/calculating face flux field phi Reading transportProperties Selecting incompressible transport model Newtonian Selecting incompressible transport model Newtonian Calculating field g.h Selecting turbulence model LRR > FOAM FATAL IO ERROR : unexpected class name volTensorField expected volSymmTensorField while reading object R file: /tirian10/lokal/OpenFOAM/dragosm1.4.1/run/tutorials/rasInterFoam//damBreak/0/R at line 21. From function regIOobject::readStream(const word&) in file db/regIOobject/regIOobjectRead.C at line 113. FOAM exiting I tried kEpsilon, RNGkEpsilon, NonlinearKEShih models, and they run well. 

November 26, 2007, 19:00 
Try replacing volTensorField w

#2 
Senior Member
Join Date: Mar 2009
Posts: 854
Rep Power: 14 
Try replacing volTensorField with volSymmTensorField in the file /tirian10/lokal/OpenFOAM/dragosm1.4.1/run/tutorials/rasInterFoam//damBreak/0/R
I will fix the file for the next release. 

November 27, 2007, 02:59 
Ok, it was that simple...http:

#3 
Senior Member
Dragos
Join Date: Mar 2009
Posts: 648
Rep Power: 12 
Ok, it was that simple..., and it wasn't actually a bug.
Well, for the information to be completed, it is necessary to change the specification of the tensor, too. In the same file the following occurences (0 0 0 0 0 0 0 0 0) have to be replaced with (0 0 0 0 0 0). Since it is a symmetric tensor, it needs only 6 value to be completely defined. And by the way, thank you Henry! 

November 27, 2007, 03:48 
...but the problems are still

#4 
Senior Member
Dragos
Join Date: Mar 2009
Posts: 648
Rep Power: 12 
...but the problems are still there .
If I keep the default settings (which work for the 2 equation models) I still get a solution, but not a physical one. Lowering the MaxCo from 0.2 to 0.02 or even lower, makes the things even worse: the computation ends with a floating point exception. Switching the discretization schemes from Gauss linear to Gauss upwind doesn't help, the same floating point exception occurs at the end. Anybody could give me a hint on how to tackle this thing? Dragos 

July 25, 2009, 11:25 

#5 
New Member
Barath Ezhilan
Join Date: Jun 2009
Posts: 20
Rep Power: 9 
Hey.. Am having a problem with simulating using LRR turbulence model in transientSimpleFoam!!
The timeStepContinuityError blows out of proportion and the simulation stops! Kindly tell me if you have any solution to this problem!! 

November 23, 2009, 12:47 

#6 
New Member
Martin Romagnoli
Join Date: Mar 2009
Location: Rosario, Santa Fe, Argentina
Posts: 22
Rep Power: 9 
Hi users, I am using OF 1.4.1 rasInterFoam solver with LRR turbulence model on damBreak case. I replaced volTensorField with volSymmTensorField in the R file but simulation stops and floating exception occurs.
Could be due to a bug in LRR.C file? There are some commented lines in Dissipation and Reynolds stress equations in LRR.C. Do they have to be commented? // Dissipation equation tmp<fvScalarMatrix> epsEqn ( fvm::ddt(epsilon_) + fvm::div(phi_, epsilon_) // fvm::laplacian(Ceps*(K/epsilon_)*R, epsilon_)  fvm::laplacian(DepsilonEff(), epsilon_) == C1*G*epsilon_/k_  fvm::Sp(C2*epsilon_/k_, epsilon_) ); epsEqn().relax(); .... // Reynolds stress equation const fvPatchList& patches = mesh_.boundary(); forAll(patches, patchi) { const fvPatch& curPatch = patches[patchi]; if (typeid(curPatch) == typeid(wallFvPatch)) { forAll(curPatch, facei) { label faceCelli = curPatch.faceCells()[facei]; P[faceCelli] *= min(G[faceCelli]/(0.5*tr(P[faceCelli]) + SMALL), 1.0); } } } tmp<fvSymmTensorMatrix> REqn ( fvm::ddt(R_) + fvm::div(phi_, R_) // fvm::laplacian(Cs*(k_/epsilon_)*R_, R_)  fvm::laplacian(DREff(), R_) + fvm::Sp(Clrr1*epsilon_/k_, R_) == P  (2.0/3.0*(1  Clrr1)*I)*epsilon_  Clrr2*dev(P) ); REqn().relax(); Thanks in advance Martín. 

November 24, 2009, 03:57 

#7 
Senior Member
Join Date: Mar 2009
Posts: 854
Rep Power: 14 
Try it with OpenFOAM1.6.x.
H 

March 28, 2011, 17:28 
R field in dambrak tutorial openfoam

#8 
New Member
Arnout
Join Date: Nov 2010
Posts: 23
Rep Power: 7 
Hi,
I'm was looking on the ras dambreak tutorial. In the /0 dir, there are a lot of boundary conditions set. Alpha1, k, epsilon, U and p_rgh are logical. But what are nut (normaly calculated from k and epsilon), mut, nuTilda and R doing there? What is the physical meaning of R [m2/s2] and mud [kg/m/s]? Is there amy documentation where I can find what is solved there? Thx!! 

Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
Turbulence model for rasInterFoam  holger_marschall  OpenFOAM Running, Solving & CFD  0  February 8, 2008 13:40 
changing model constants in ke turbulence model  Sunil  CFX  3  October 3, 2006 12:12 
ERRORS in rasInterFoam Turbulence  kumar2  OpenFOAM Running, Solving & CFD  0  June 9, 2006 15:06 
RasInterFoam cavitation  maritozzo  OpenFOAM Running, Solving & CFD  2  December 6, 2005 15:09 
HELP! TURBULENCE ke OR komega TURBULENCE MODEL?  Mirek Kabacinski  FLUENT  5  August 24, 2003 22:31 