
[Sponsors] 
February 8, 2017, 18:32 
Error in LRR (RAS) model in OpenFOAM v4.1

#1 
New Member
Edmond
Join Date: Jan 2017
Posts: 2
Rep Power: 0 
I was trying some simple tutorial about a backwardfacing step in steady state via different turbulence models in RAS using the simpleFoam solver. One of which was the LRR (Reynolds stress model). The tutorials themselves were done so in another version of OpenFOAM and the one I am using was the latest version: 4.1.
All else worked splendidly, however in the case of the LRR, I keep getting an error returned from the fvSchemes file. Basically, the steps in this method that I am following are: after copying the kEpsilon file (by file I mean the whole kEpsilon simulation folder  i.e. the directory containing: 0/ constant/ system/ ), I execute the command: Code:
simpleFoam postProcess func R Next, I change the turbulence model from kEpsilon to LRR as intended. Afterwards, I edit the fvSchemes and include the following for the divergence terms: Code:
divSchemes { default none; div(phi,U) bounded Gauss linearUpwind grad(U); div(phi,k) bounded Gauss linearUpwind grad(k); div(phi,epsilon) bounded Gauss linearUpwind grad(epsilon); div(phi,omega) bounded Gauss limitedLinear 1; div(phi,v2) bounded Gauss limitedLinear 1; div(phi,R) bounded Gauss upwind; div(R) Gauss linear; div((nuEff*dev2(T(grad(U))))) Gauss linear; div(nonlinearStress) Gauss linear; } Code:
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // Create time Create mesh for time = 0 SIMPLE: convergence criteria field p tolerance 0.01 field U tolerance 0.001 field "(kepsilonomegafv2R)" tolerance 0.001 Reading field p Reading field U Reading/calculating face flux field phi Selecting incompressible transport model Newtonian Selecting turbulence model type RAS Selecting RAS turbulence model LRR > FOAM FATAL IO ERROR: cannot find file file: /home/username/OpenFOAM/username4.1/run/tutorials/incompressible/simpleFoam/bfs/BFS_baseCase/0/R at line 0. From function regIOobject::readStream() in file db/regIOobject/regIOobjectRead.C at line 72. FOAM exiting My question hence is why did this happen? is there anything wrong with my implementation or is it a bug ? Thanks ! 

July 29, 2018, 22:28 

#2 
Senior Member
Ruiyan Chen
Join Date: Jul 2016
Location: Hangzhou, China
Posts: 162
Rep Power: 9 
Take a look at this thread
fvSchemes According to Dr. Jasak, "It does not recognise it because this is a straightforward div(thing) term and not convection div(phi, thing). Therefore you cannot do upwinding, since there's no flux to decide the direction from." 

July 2, 2020, 15:09 

#3 
New Member
Giacomo Del Bianco
Join Date: Jul 2020
Location: Venice area  Italy
Posts: 2
Rep Power: 0 
Thank you very much!! You fixed my problem.
I saw browsing from the internet but I didn't realized it was different from the one: Code:
div((nuEff*dev2(T(grad(U))))) Gauss linear; Code:
keyword wallDist is undefined in dictionary Code:
wallDist { method meshWave; } Thank you very much again!! Giacomo 

March 1, 2022, 09:39 

#4  
New Member
SunTime
Join Date: Nov 2020
Posts: 14
Rep Power: 4 
Quote:


Tags 
fvschemes, openfoam4.0, ras 
Thread Tools  Search this Thread 
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
Simulating an aether model of EM using OpenFoam?  lamare  Electromagnetics  2  September 15, 2016 05:08 
Prandtl's oneequation RAS turbulence model  Djub  OpenFOAM Running, Solving & CFD  3  April 9, 2016 06:30 
about Subgridscale model  impecca  OpenFOAM Running, Solving & CFD  4  December 20, 2013 10:36 
OpenFOAM: LES turbulence model names  Ollie  OpenFOAM  5  January 7, 2013 10:24 
Water subcooled boiling  Attesz  CFX  7  January 5, 2013 03:32 