CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Bugs

RasInterFoam and LRR turbulence model

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   November 26, 2007, 06:05
Default Description: Although LRR is
  #1
Senior Member
 
dmoroian's Avatar
 
Dragos
Join Date: Mar 2009
Posts: 648
Rep Power: 20
dmoroian is on a distinguished road
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/dragosm-1.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.
dmoroian is offline   Reply With Quote

Old   November 26, 2007, 19:00
Default Try replacing volTensorField w
  #2
Senior Member
 
Join Date: Mar 2009
Posts: 854
Rep Power: 22
henry is on a distinguished road
Try replacing volTensorField with volSymmTensorField in the file /tirian10/lokal/OpenFOAM/dragosm-1.4.1/run/tutorials/rasInterFoam//damBreak/0/R

I will fix the file for the next release.
henry is offline   Reply With Quote

Old   November 27, 2007, 02:59
Default Ok, it was that simple...http:
  #3
Senior Member
 
dmoroian's Avatar
 
Dragos
Join Date: Mar 2009
Posts: 648
Rep Power: 20
dmoroian is on a distinguished road
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!
dmoroian is offline   Reply With Quote

Old   November 27, 2007, 03:48
Default ...but the problems are still
  #4
Senior Member
 
dmoroian's Avatar
 
Dragos
Join Date: Mar 2009
Posts: 648
Rep Power: 20
dmoroian is on a distinguished road
...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
dmoroian is offline   Reply With Quote

Old   July 25, 2009, 12:25
Default
  #5
New Member
 
Barath Ezhilan
Join Date: Jun 2009
Posts: 20
Rep Power: 17
barath.ezhilan is on a distinguished road
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!!
barath.ezhilan is offline   Reply With Quote

Old   November 23, 2009, 12:47
Default
  #6
New Member
 
Martin Romagnoli
Join Date: Mar 2009
Location: Rosario, Santa Fe, Argentina
Posts: 22
Rep Power: 17
martinr is on a distinguished road
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.
martinr is offline   Reply With Quote

Old   November 24, 2009, 03:57
Default
  #7
Senior Member
 
Join Date: Mar 2009
Posts: 854
Rep Power: 22
henry is on a distinguished road
Try it with OpenFOAM-1.6.x.

H
henry is offline   Reply With Quote

Old   March 28, 2011, 18:28
Default R field in dambrak tutorial openfoam
  #8
Member
 
Arnout
Join Date: Nov 2010
Posts: 46
Rep Power: 16
The King is on a distinguished road
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!!
The King 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
Turbulence model for rasInterFoam holger_marschall OpenFOAM Running, Solving & CFD 0 February 8, 2008 13:40
changing model constants in k-e turbulence model Sunil CFX 3 October 3, 2006 13:12
ERRORS in rasInterFoam Turbulence kumar2 OpenFOAM Running, Solving & CFD 0 June 9, 2006 16:06
RasInterFoam cavitation maritozzo OpenFOAM Running, Solving & CFD 2 December 6, 2005 15:09
HELP! TURBULENCE k-e OR k-omega TURBULENCE MODEL? Mirek Kabacinski FLUENT 5 August 24, 2003 23:31


All times are GMT -4. The time now is 01:06.