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

Usage of R and LRR

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

Like Tree14Likes
  • 2 Post By GGerber
  • 4 Post By GGerber
  • 3 Post By s.m
  • 5 Post By amuzeshi

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   March 27, 2011, 05:11
Default Usage of R and LRR
  #1
New Member
 
George
Join Date: Oct 2010
Posts: 18
Rep Power: 15
GGerber is on a distinguished road
Hi Foamers,

I have a wall-jet simulation, which uses the kOmegaSST turbulence model. The model performs ok and I would now like to investigate the Reynolds stresses for this flow. Therefore I have modified my model to use the LLR turbulence model. The LRR simulation blows up when I run it from time 0.

From the forums I gather that I can use the R utility to derive the Reynolds stresses fields from the kOmegaSST results. This will set up the 'initial conditions' for the LRR model.

However, when I run the command "R" from the casedir the utility simply zip's the R file that I placed in the directory of the last time-step. I would have expected that the utility would modify the 'internalField' entry in the R-file, but it is left unchanged...

If someone can tell me how to properly use R and setup LRR it would be much appreciated. Below is my R file and controlDict file.

Thanks,
George

~~~~~~~~
FoamFile
{
version 2.0;
format ascii;
class volSymmTensorField;
object R;
}
// * * * * * * * * * * * * * * * * * * * * * //

dimensions [0 2 -2 0 0 0 0];
internalField uniform (0 0 0 0 0 0);
boundaryField
{
inlet {
type fixedValue;
value uniform (0.0015325 0 0 0.0015325 0 0.00306504); // (k/2 0 0 k/2 0 k);
}
walls
{
type kqRWallFunction;
value uniform ( 0 0 0 0 0 0 );
}
bed
{
type kqRWallFunction;
value uniform ( 0 0 0 0 0 0 );
}
top
{
type symmetryPlane;
}
outlet
{
type zeroGradient;
}
}
~~~~~~~~~~~~
FoamFile
{
version 2.0;
format ascii;
class dictionary;
object controlDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * //
application simpleFoam;
startFrom latestTime;
startTime 0;
stopAt endTime;
endTime 2820;
deltaT 1;
writeControl timeStep;
writeInterval 20;//100;
purgeWrite 0;
writeFormat ascii;
writePrecision 6;
writeCompression compressed;
timeFormat general;
timePrecision 6;
runTimeModifiable yes;
~~~~~~~~~~~~~`
GGerber is offline   Reply With Quote

Old   March 29, 2011, 03:47
Default
  #2
New Member
 
George
Join Date: Oct 2010
Posts: 18
Rep Power: 15
GGerber is on a distinguished road
Hi Foamers,

I have solved my problem.
Issue closed.

Best regards,
George
s.m and lpz456 like this.
GGerber is offline   Reply With Quote

Old   April 8, 2011, 18:13
Default About R
  #3
New Member
 
Sergio
Join Date: Apr 2011
Posts: 8
Rep Power: 14
Sergio13 is on a distinguished road
Hi George
I am quite new with the CFD and OpenFoam, and I need to calculate the Reynolds Stresses, and if you can help me with the following questions, I would really appreciate it.

1. Is necessary to create a R file?
2. If so, should the location for this file be the system directory?
3. How will the output file be created? (I mean, is it required to run again the simulation?)

I would also appreciate if you can give me some instructions about the usage of R.

Thanks
Sergio13 is offline   Reply With Quote

Old   April 10, 2011, 17:09
Default
  #4
New Member
 
George
Join Date: Oct 2010
Posts: 18
Rep Power: 15
GGerber is on a distinguished road
Quote:
Originally Posted by Sergio13 View Post
Hi George
I am quite new with the CFD and OpenFoam, and I need to calculate the Reynolds Stresses, and if you can help me with the following questions, I would really appreciate it.

1. Is necessary to create a R file?
2. If so, should the location for this file be the system directory?
3. How will the output file be created? (I mean, is it required to run again the simulation?)

I would also appreciate if you can give me some instructions about the usage of R.

Thanks
Hi Sergio,
1) No, I dont think its necessary to create a R file. Just run 'R' from the command line after you finished a k-epsilon simulation. I think you also need to set the ControlDict keyword endTime (or lastTime?) to latestTime.
2) I have not investigated how the Reynolds stresses are computed from the k-epsilon results. You will probably see in the 'R' application code how it is computed. The Boussinesq hypothesis is probably used (using the effective viscosity and mean velocity gradients) to compute the individual Reynolds stresses.

HTH,
George
GGerber is offline   Reply With Quote

Old   April 10, 2011, 17:33
Default
  #5
New Member
 
Sergio
Join Date: Apr 2011
Posts: 8
Rep Power: 14
Sergio13 is on a distinguished road
Hi George.
First than all, i want to thank you for take time to read my question, and answer it.

On the other hand I have realized that I am using a LES turbulence model, and thus my simulation is not steady state. This explain why after I typed R an error message, asking for the RASProperties.

Have you ever tried to calculate Reynolds stresses in such a simulation (LES)?

Thanks again for your help
Sergio13 is offline   Reply With Quote

Old   April 12, 2011, 03:39
Default
  #6
New Member
 
George
Join Date: Oct 2010
Posts: 18
Rep Power: 15
GGerber is on a distinguished road
Hi Sergio,

Sorry I do not have any experience with LES.

Good luck,
George
GGerber is offline   Reply With Quote

Old   April 12, 2011, 11:36
Default
  #7
New Member
 
Sergio
Join Date: Apr 2011
Posts: 8
Rep Power: 14
Sergio13 is on a distinguished road
Thanks for your help!
Sergio13 is offline   Reply With Quote

Old   May 7, 2012, 11:49
Default
  #8
Member
 
Join Date: Jun 2011
Posts: 42
Rep Power: 14
mikeP is on a distinguished road
Hi,

I have the same problem as George. I have a case solved with k-epsilon. When I run the utility R, it does not give any warnings or errors and completes successfully. Afterwards when I check the field R, it remains unchanged.
It seems that George has found the answer to it, but how is it?
mikeP is offline   Reply With Quote

Old   October 16, 2013, 04:36
Default
  #9
s.m
Senior Member
 
saeideh mohamadi
Join Date: Aug 2012
Posts: 229
Rep Power: 14
s.m is on a distinguished road
Quote:
Originally Posted by GGerber View Post
Hi Foamers,

I have solved my problem.
Issue closed.

Best regards,
George
Hi GGerber,
What did you do to solve your promblem, could you please tell me, i really need some guidance to solve my problem,
thank you very much
i simulate the flow over an airfoil, using the simpleFoam solver, my solution be complete with kOmegaSST, but it become diverged when i change the kOmegaSST to LRR, i don't know what should i do.
This is the divergence error:
Time = 14

smoothSolver: Solving for Ux, Initial residual = 0.348073, Final residual = 6.66759e-10, No Iterations 69
smoothSolver: Solving for Uy, Initial residual = 0.520936, Final residual = 9.19728e-10, No Iterations 68
GAMG: Solving for p, Initial residual = 1, Final residual = 0.256866, No Iterations 100
GAMG: Solving for p, Initial residual = 1, Final residual = 8.38804e-11, No Iterations 25
time step continuity errors : sum local = 5.57903e+32, global = -3.39894e+28, cumulative = -3.39894e+28
#0 Foam::error:rintStack(Foam::Ostream&) at ??:?
#1 Foam::sigFpe::sigHandler(int) at ??:?
#2 in "/lib64/libc.so.6"
#3 Foam::GaussSeidelSmoother::smooth(Foam::word const&, Foam::Field<double>&, Foam::lduMatrix const&, Foam::Field<double> const&, Foam::FieldField<Foam::Field, double> const&, Foam::UPtrList<Foam::lduInterfaceField const> const&, unsigned char, int) at ??:?
#4 Foam::GaussSeidelSmoother::smooth(Foam::Field<doub le>&, Foam::Field<double> const&, unsigned char, int) const at ??:?
#5 Foam::smoothSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
#6 Foam::fvMatrix<double>::solveSegregated(Foam::dict ionary const&) at ??:?
#7 Foam::fvMatrix<double>::solve(Foam::dictionary const&) at ??:?
#8 Foam::fvMatrix<double>::solve() at ??:?
#9 Foam::SolverPerformance<double> Foam::solve<double>(Foam::tmp<Foam::fvMatrix<doubl e> > const&) at ??:?
#10 Foam::incompressible::RASModels::LRR::correct() at ??:?
#11
at ??:?
#12 __libc_start_main in "/lib64/libc.so.6"
#13
at /home/abuild/rpmbuild/BUILD/glibc-2.17/csu/../sysdeps/x86_64/start.S:126
Floating point exception
s.m is offline   Reply With Quote

Old   November 23, 2016, 04:38
Default
  #10
Senior Member
 
ali
Join Date: Jul 2016
Posts: 147
Rep Power: 9
ebtedaei is on a distinguished road
Quote:
Originally Posted by GGerber View Post
Hi Foamers,

I have solved my problem.
Issue closed.

Best regards,
George
Dear GGerber,
I need your help.
Please send me your Email address.

Thanks
Ali
ebtedaei is offline   Reply With Quote

Old   May 25, 2018, 13:15
Default Lrr
  #11
Senior Member
 
Ali Shayegh
Join Date: Oct 2015
Posts: 130
Rep Power: 10
amuzeshi is on a distinguished road
This is what you need.
Error in LRR (RAS) model in OpenFOAM v4.1
follow the steps; it works. Note that you must apply the following change (RED)also in fvSolution:


Code:
solvers
{
    p
    {
        solver          GAMG;
        tolerance       1e-06;
        relTol          0.1;
        smoother        GaussSeidel;
    }

    "(U|k|epsilon|omega|f|v2|R)"
    {
        solver          smoothSolver;
        smoother        symGaussSeidel;
        tolerance       1e-05;
        relTol          0.1;
    }
}

SIMPLE
{
    nNonOrthogonalCorrectors 0;
    consistent      yes;

    residualControl
    {
        p               1e-2;
        U               1e-3;
        "(k|epsilon|omega|f|v2|R)" 1e-3;
    }
}
faiazk, samuel_rff, Maahs and 2 others like this.
amuzeshi is offline   Reply With Quote

Old   January 29, 2024, 08:33
Default
  #12
New Member
 
Ellen Caroline
Join Date: Nov 2023
Posts: 2
Rep Power: 0
Ellen Caroline is on a distinguished road
HI! I'm new in openfoam. I ran "R" in terminal and I got this message instead an R file... can you help me, please?



openfoam2312:~/Documents/Ellen/testesKepsilon/teste6/
n003$ R

R version 4.1.2 (2021-11-01) -- "Bird Hippie"
Copyright (C) 2021 The R Foundation for Statistical Computing
Platform: x86_64-pc-linux-gnu (64-bit)

R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.

Natural language support but running in an English locale

R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.

Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.
Ellen Caroline is offline   Reply With Quote

Old   January 29, 2024, 09:55
Default
  #13
Senior Member
 
Join Date: Apr 2020
Location: UK
Posts: 663
Rep Power: 13
Tobermory will become famous soon enough
You typed R in a terminal window? It looks like you have the R programming language (https://en.wikipedia.org/wiki/R_(programming_language)) installed, and so the OS tried to run this. Clearly not what you wanted!

Try running postProcess with the R functionObject:

Code:
postProcess -func R
or you may need to do this through the solver, eg:

Code:
<solverName> -postProcess -func R
Tobermory is offline   Reply With Quote

Reply

Tags
lrr, rasmodel

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
RasInterFoam and LRR turbulence model dmoroian OpenFOAM Bugs 7 March 28, 2011 18:28
LRR Rstm pipe flow takes too long to get fully developed kjetil OpenFOAM Running, Solving & CFD 0 February 13, 2010 12:26
what is the special need for LRR RSM model? yuhai OpenFOAM Running, Solving & CFD 0 October 9, 2009 07:05
Problems with Unsteady simulations using LRR Model in transientSimpleFoam barath.ezhilan OpenFOAM Running, Solving & CFD 1 July 25, 2009 14:58
LRR model andimb OpenFOAM Running, Solving & CFD 0 March 17, 2006 05:29


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