|
[Sponsors] |
November 26, 2007, 06:05 |
Description:
Although LRR is
|
#1 |
Senior Member
Dragos
Join Date: Mar 2009
Posts: 648
Rep Power: 20 |
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. |
|
November 26, 2007, 19:00 |
Try replacing volTensorField w
|
#2 |
Senior Member
Join Date: Mar 2009
Posts: 854
Rep Power: 22 |
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. |
|
November 27, 2007, 02:59 |
Ok, it was that simple...http:
|
#3 |
Senior Member
Dragos
Join Date: Mar 2009
Posts: 648
Rep Power: 20 |
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: 20 |
...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, 12:25 |
|
#5 |
New Member
Barath Ezhilan
Join Date: Jun 2009
Posts: 20
Rep Power: 17 |
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: 17 |
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: 22 |
Try it with OpenFOAM-1.6.x.
H |
|
March 28, 2011, 18:28 |
R field in dambrak tutorial openfoam
|
#8 |
Member
Arnout
Join Date: Nov 2010
Posts: 46
Rep Power: 16 |
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 | Search this Thread |
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 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 |