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

y+ and u+ values with low-Re RANS turbulence models: utility + testcase

Register Blogs Community New Posts Updated Threads Search

Like Tree81Likes

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   November 24, 2012, 11:45
Default
  #61
atg
Member
 
Karl
Join Date: Jan 2011
Posts: 36
Rep Power: 15
atg is on a distinguished road
Thank you Aljazari for your reply.

I followed your instructions. Different error this time. I am on the latest SGI version of OF; perhaps that is the issue.

Here is the terminal output:

karl@Karl-Xi:~/Downloads/yPlus$ wmake
Making dependency list for source file yPlus.C
SOURCE=yPlus.C ; g++ -m64 -Dlinux64 -DWM_DP -Wall -Wextra -Wno-unused-parameter -Wold-style-cast -Wnon-virtual-dtor -O3 -DNoRepository -ftemplate-depth-100 -I/home/karl/OpenFOAM/OpenFOAM-2.1.1/src/meshTools/lnInclude -I/home/karl/OpenFOAM/OpenFOAM-2.1.1/src/turbulenceModels -I/home/karl/OpenFOAM/OpenFOAM-2.1.1/src/transportModels -I/home/karl/OpenFOAM/OpenFOAM-2.1.1/src/thermophysicalModels/basic/lnInclude -I/home/karl/OpenFOAM/OpenFOAM-2.1.1/src/finiteVolume/lnInclude -IlnInclude -I. -I/home/karl/OpenFOAM/OpenFOAM-2.1.1/src/OpenFOAM/lnInclude -I/home/karl/OpenFOAM/OpenFOAM-2.1.1/src/OSspecific/POSIX/lnInclude -fPIC -c $SOURCE -o Make/linux64GccDPOpt/yPlus.o
g++ -m64 -Dlinux64 -DWM_DP -Wall -Wextra -Wno-unused-parameter -Wold-style-cast -Wnon-virtual-dtor -O3 -DNoRepository -ftemplate-depth-100 -I/home/karl/OpenFOAM/OpenFOAM-2.1.1/src/meshTools/lnInclude -I/home/karl/OpenFOAM/OpenFOAM-2.1.1/src/turbulenceModels -I/home/karl/OpenFOAM/OpenFOAM-2.1.1/src/transportModels -I/home/karl/OpenFOAM/OpenFOAM-2.1.1/src/thermophysicalModels/basic/lnInclude -I/home/karl/OpenFOAM/OpenFOAM-2.1.1/src/finiteVolume/lnInclude -IlnInclude -I. -I/home/karl/OpenFOAM/OpenFOAM-2.1.1/src/OpenFOAM/lnInclude -I/home/karl/OpenFOAM/OpenFOAM-2.1.1/src/OSspecific/POSIX/lnInclude -fPIC -Xlinker --add-needed -Xlinker --no-as-needed Make/linux64GccDPOpt/yPlus.o -L/home/karl/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib \
-lcompressibleTurbulenceModel -lcompressibleRASModels -lcompressibleLESModels -lincompressibleTurbulenceModel -lincompressibleTransportModels -lincompressibleRASModels -lincompressibleLESModels -lfiniteVolume -lgenericPatchFields -lOpenFOAM -ldl -lm -o /home/karl/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/bin/yPlus
/usr/bin/ld: error: --add-needed is not supported but is required for libbasicSolidThermo.so in /home/karl/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libcompressibleTurbulenceModel.so
Make/linux64GccDPOpt/yPlus.o:yPlus.C:function main: error: undefined reference to 'Foam::basicThermo::New(Foam::fvMesh const&)'
collect2: ld returned 1 exit status
make: *** [/home/karl/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/bin/yPlus] Error 1
karl@Karl-Xi:~/Downloads/yPlus$ sudo bash
[sudo] password for karl:
root@Karl-Xi:~/Downloads/yPlus# wmake
g++ -m64 -Dlinux64 -DWM_DP -Wall -Wextra -Wno-unused-parameter -Wold-style-cast -Wnon-virtual-dtor -O3 -DNoRepository -ftemplate-depth-100 -I/home/karl/OpenFOAM/OpenFOAM-2.1.1/src/meshTools/lnInclude -I/home/karl/OpenFOAM/OpenFOAM-2.1.1/src/turbulenceModels -I/home/karl/OpenFOAM/OpenFOAM-2.1.1/src/transportModels -I/home/karl/OpenFOAM/OpenFOAM-2.1.1/src/thermophysicalModels/basic/lnInclude -I/home/karl/OpenFOAM/OpenFOAM-2.1.1/src/finiteVolume/lnInclude -IlnInclude -I. -I/home/karl/OpenFOAM/OpenFOAM-2.1.1/src/OpenFOAM/lnInclude -I/home/karl/OpenFOAM/OpenFOAM-2.1.1/src/OSspecific/POSIX/lnInclude -fPIC -Xlinker --add-needed -Xlinker --no-as-needed Make/linux64GccDPOpt/yPlus.o -L/home/karl/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib \
-lcompressibleTurbulenceModel -lcompressibleRASModels -lcompressibleLESModels -lincompressibleTurbulenceModel -lincompressibleTransportModels -lincompressibleRASModels -lincompressibleLESModels -lfiniteVolume -lgenericPatchFields -lOpenFOAM -ldl -lm -o /home/karl/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/bin/yPlus
/usr/bin/ld: error: --add-needed is not supported but is required for libbasicSolidThermo.so in /home/karl/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/lib/libcompressibleTurbulenceModel.so
Make/linux64GccDPOpt/yPlus.o:yPlus.C:function main: error: undefined reference to 'Foam::basicThermo::New(Foam::fvMesh const&)'
collect2: ld returned 1 exit status
make: *** [/home/karl/OpenFOAM/OpenFOAM-2.1.1/platforms/linux64GccDPOpt/bin/yPlus] Error 1
atg is offline   Reply With Quote

Old   November 25, 2012, 01:21
Default
  #62
New Member
 
Abdulhaq
Join Date: Oct 2012
Posts: 27
Rep Power: 14
aljazari is on a distinguished road
Quote:
Originally Posted by atg View Post
/usr/bin/ld: error: --add-needed is not supported but is required for libbasicSolidThermo.so in /home/karl/OpenFOAM/OpenFOAM-
Check this thread, it has an error similar to the one above although that person was compiling everything.

Make sure you have all the required build tools:

Quote:
Originally Posted by wyldckat View Post
Code:
sudo apt-get install build-essential flex bison cmake zlib1g-dev qt4-dev-tools libqt4-dev gnuplot libreadline-dev libncurses-dev libxt-dev
aljazari is offline   Reply With Quote

Old   November 25, 2012, 18:53
Default
  #63
atg
Member
 
Karl
Join Date: Jan 2011
Posts: 36
Rep Power: 15
atg is on a distinguished road
Quote:
Originally Posted by aljazari View Post
Check this thread, it has an error similar to the one above although that person was compiling everything.

Make sure you have all the required build tools:
I do have the required tools thanks.

That makes sense I did add the "gold" package discussed in the thread you referenced, in order to build the latest CUDA. OF has been giving me problems ever since!

I just reinstalled Ubuntu and am recompiling 2.1.x. We will see how it goes; I do not need CUDA for the moment.

Regards and Thanks again,
Karl
atg is offline   Reply With Quote

Old   November 26, 2012, 10:54
Default
  #64
Senior Member
 
Bernhard
Join Date: Sep 2009
Location: Delft
Posts: 790
Rep Power: 22
Bernhard is on a distinguished road
Quote:
Originally Posted by lakeat View Post
Okay, FOAM community,

Since I feel there is a little bit confusion here, and I agree with Jose's concern in this post, so I just rewrote them based on OpenCFD's version of OpenFOAM. You can download them here,
Thanks, I think it is good to have a topic which makes these thing clear.

I have some question on your implementation. What is the reason for using nearWallDist instead of wallDist? Am I correct in concluding that nearWallDist is demanding less memory?
Bernhard is offline   Reply With Quote

Old   November 26, 2012, 11:25
Default
  #65
Senior Member
 
lakeat's Avatar
 
Daniel WEI (老魏)
Join Date: Mar 2009
Location: Beijing, China
Posts: 689
Blog Entries: 9
Rep Power: 21
lakeat is on a distinguished road
Send a message via Skype™ to lakeat
Good catch, that line of include nearwalldist is redundant. Thanks
__________________
~
Daniel WEI
-------------
Boeing Research & Technology - China
Beijing, China
Email
lakeat is offline   Reply With Quote

Old   November 27, 2012, 09:38
Default
  #66
New Member
 
Join Date: Jun 2012
Posts: 19
Rep Power: 14
asonda is on a distinguished road
Hi all,

I am trying to understand the utilities kindly shared by Daniel Wei (thank you so much!).

In the yPlus utility, the implemented equation for the compressible case is presented at the end of this post.

Please double-check this simple computation: if we combine eqs 1-4 presented below ('w' represents wall values), we should get :

y+ = sqrt(muEff * dUmag/dn) * y / mut_w * sqrt(rho_w)

I am not very practical with object-oriented languages, but to me it seems like the factor (sqrt(rho_w)) is not included in the current code implementation. . If I am right, then one easy modification would fix it; if I am wrong, please forgive me.

Alberto

ps: The incompressible formula looks good to me (the density term coming from [3] simplifies with [1] leaving nuEff)
  1. muEff = mu + mut; // effective viscosity
  2. Tau_w = (muEff * dUmag/dn)_w; // wall shear stress
  3. u_t = sqrt(Tau_w / rho_w); // wall friction velocity
  4. y+ = y*u_t / nu = y*u_t*rho_w / mut_w; // y+

if (compressible)
{
yPlus.boundaryField()[patchi] =
d[patchi]
*sqrt
(
muEff.boundaryField()[patchi]
*mag(U.boundaryField()[patchi].snGrad())
)
/muLam.boundaryField()[patchi];
}
asonda is offline   Reply With Quote

Old   November 27, 2012, 11:34
Default
  #67
Senior Member
 
lakeat's Avatar
 
Daniel WEI (老魏)
Join Date: Mar 2009
Location: Beijing, China
Posts: 689
Blog Entries: 9
Rep Power: 21
lakeat is on a distinguished road
Send a message via Skype™ to lakeat
Quote:
Originally Posted by asonda View Post
Hi all,

I am trying to understand the utilities kindly shared by Daniel Wei (thank you so much!).

In the yPlus utility, the implemented equation for the compressible case is presented at the end of this post.

Please double-check this simple computation: if we combine eqs 1-4 presented below ('w' represents wall values), we should get :

y+ = sqrt(muEff * dUmag/dn) * y / mut_w * sqrt(rho_w)

I am not very practical with object-oriented languages, but to me it seems like the factor (sqrt(rho_w)) is not included in the current code implementation. . If I am right, then one easy modification would fix it; if I am wrong, please forgive me.

Alberto

ps: The incompressible formula looks good to me (the density term coming from [3] simplifies with [1] leaving nuEff)
  1. muEff = mu + mut; // effective viscosity
  2. Tau_w = (muEff * dUmag/dn)_w; // wall shear stress
  3. u_t = sqrt(Tau_w / rho_w); // wall friction velocity
  4. y+ = y*u_t / nu = y*u_t*rho_w / mut_w; // y+

if (compressible)
{
yPlus.boundaryField()[patchi] =
d[patchi]
*sqrt
(
muEff.boundaryField()[patchi]
*mag(U.boundaryField()[patchi].snGrad())
)
/muLam.boundaryField()[patchi];
}
Thanks! I don't know much about the compressible flow, sorry for that. Please let me know what equation should I use.
__________________
~
Daniel WEI
-------------
Boeing Research & Technology - China
Beijing, China
Email
lakeat is offline   Reply With Quote

Old   November 28, 2012, 04:12
Default
  #68
Senior Member
 
Bernhard
Join Date: Sep 2009
Location: Delft
Posts: 790
Rep Power: 22
Bernhard is on a distinguished road
I think some other people are using the utility to refine the mesh based on the value of y+ on a patch. I have been using refineWallLayer for this. But, as you sometimes have a patch with a large range of y+ values, it is not always desired to refine all those cells.

I figured out the following procedure to refine wall-adjacent cells based on the value of y+.

1.
Code:
                forAll(yPlus.boundaryField()[patchi], faceI)
                {   
                    label faceCellI = currPatch.faceCells()[faceI];
                    yPlus[faceCellI]=yPlus.boundaryField()[patchi][faceI];
                }
Add these lines in the above utility by Lakeat, just below the line where yPlus is calculated. It just copies the patch-values of yPlus to the internal field.

2.
Put the following in system/topoSetDict, to select the cells that you want to keep unchanged. (and execute it with e.g topoSet -latestTime )
Code:
actions
(
    {
        name    highYPlus;
        type    cellSet;
        action  new;
        source  fieldToCell;
        sourceInfo
        {
            fieldName   yPlus;
            min         0.0;
            max         3.0;

        }
    }
);
3.
Then, use refineWallLayer refine the cells with largest y+ values calculated by the utility.
refineWallLayer -useSet highYPlus -overwrite upperWall 0.5

I hope it is of good use to anyone. Kudo's for the one to put it into a single utility
Bernhard is offline   Reply With Quote

Old   November 28, 2012, 14:51
Default
  #69
New Member
 
Abdulhaq
Join Date: Oct 2012
Posts: 27
Rep Power: 14
aljazari is on a distinguished road
Quote:
Originally Posted by asonda View Post
Hi all,

I am trying to understand the utilities kindly shared by Daniel Wei (thank you so much!).

In the yPlus utility, the implemented equation for the compressible case is presented at the end of this post.

Please double-check this simple computation: if we combine eqs 1-4 presented below ('w' represents wall values), we should get :

y+ = sqrt(muEff * dUmag/dn) * y / mut_w * sqrt(rho_w)

I am not very practical with object-oriented languages, but to me it seems like the factor (sqrt(rho_w)) is not included in the current code implementation. . If I am right, then one easy modification would fix it; if I am wrong, please forgive me.

Alberto

ps: The incompressible formula looks good to me (the density term coming from [3] simplifies with [1] leaving nuEff)
  1. muEff = mu + mut; // effective viscosity
  2. Tau_w = (muEff * dUmag/dn)_w; // wall shear stress
  3. u_t = sqrt(Tau_w / rho_w); // wall friction velocity
  4. y+ = y*u_t / nu = y*u_t*rho_w / mut_w; // y+

if (compressible)
{
yPlus.boundaryField()[patchi] =
d[patchi]
*sqrt
(
muEff.boundaryField()[patchi]
*mag(U.boundaryField()[patchi].snGrad())
)
/muLam.boundaryField()[patchi];
}
I think you're correct, there's a rho term missing in the denominator.



I would appreciate it lakeat if you included it.
aljazari is offline   Reply With Quote

Old   November 28, 2012, 14:53
Default
  #70
Senior Member
 
lakeat's Avatar
 
Daniel WEI (老魏)
Join Date: Mar 2009
Location: Beijing, China
Posts: 689
Blog Entries: 9
Rep Power: 21
lakeat is on a distinguished road
Send a message via Skype™ to lakeat
That has been done already, I will update it today.

Right now, I am considering the suggestion that Bernhard has made. Even though I am not sure if it is a good idea to make yPlus too lengthy.
aljazari likes this.
__________________
~
Daniel WEI
-------------
Boeing Research & Technology - China
Beijing, China
Email
lakeat is offline   Reply With Quote

Old   November 28, 2012, 19:22
Default
  #71
Senior Member
 
lakeat's Avatar
 
Daniel WEI (老魏)
Join Date: Mar 2009
Location: Beijing, China
Posts: 689
Blog Entries: 9
Rep Power: 21
lakeat is on a distinguished road
Send a message via Skype™ to lakeat
Hi again,

I think I have fixed the bug.

yPlus utility: yPlus.tar.gz
yStar utility: yStar.tar.gz

Equations used for incompressible flow:

Code:
y^+=\frac{y\sqrt{\frac{\mu_{eff}\left(\frac{\partial U}{\partial y}\right)_{y=0}}{\rho}}}{\frac{\mu_{lam}}{\rho}}=\frac{y\sqrt{{\nu_{eff}\left(\frac{\partial U}{\partial y}\right)_{y=0}}}}{\nu_{lam}}
Equations used for compressible flow:


Code:
y^+=\frac{y\sqrt{\frac{\mu_{eff}\left(\frac{\partial U}{\partial y}\right)_{y=0}}{\rho}}}{\frac{\mu_{lam}}{\rho}}=\frac{y\sqrt{\rho{\mu_{eff}\left(\frac{\partial U}{\partial y}\right)_{y=0}}}}{\mu_{lam}}

Here are a few ideas concerning combining yPlus and refine wall layer mesh together:

1. First, Thanks, Bernhard. I have written another utility to refine the wall layer mesh based on yPlus/yStar value. And I have posted to another thread, since this is just a thread discussing yPlus. Let's make yPlus light-weighted.

2. Due to my limited knowledge, I don't see a need to use cellSet to choose which wall cells to be refined. IMHO, all cells on the wall should be simultaneously refined.

3. Again, I think one should better rebuild his mesh if he found y+ too high. That is to say, from turb modeling point of view, I think refineWallLayer utility is not the best method to refine the mesh. Because the abrupt change due to this utility, and this resulting abrupt change would most probably in the boundary layer. Some turb models will even have troubles with this.


I'll post the utility here too, but more details or discussion I think another thread would be a better fit.
refineWallLayerSmart utility: refineWallLayerSmart.tar.gz
__________________
~
Daniel WEI
-------------
Boeing Research & Technology - China
Beijing, China
Email

Last edited by lakeat; November 28, 2012 at 19:48.
lakeat is offline   Reply With Quote

Old   December 12, 2012, 21:45
Default
  #72
Member
 
sqing
Join Date: Sep 2012
Location: Dalian
Posts: 77
Rep Power: 14
Sunxing is on a distinguished road
Quote:
Originally Posted by jona View Post
Hi *,

for converting the code of Florian to OF 2.0, "RASModel->nu()" had to be explicitely treated as "volScalarField&".

The corrected code is attached.

Have fun,

Jona
Hi Jona,

I have downloaded the utility you posted and it works well for my OF(version 2.1.1).
Now I want to put out the wall shear stress so I have added a the following sentence to the plusPostRANS.C

Info << " avg. friction velocity uTau is: "
<< uTauAvg.value() << " (averaged over " << nPatch << " wall(s))" << nl <<endl;

//adding*******************************
Info << " wall shear stress is: "
<< uTauAvg.value()*uTauAvg.value() << " (averaged over " << nPatch << " wall(s))" << nl <<endl;
//adding********************************

yPlus = y.y() * uTauAvg / RASModel->nu();
uPlus = U / uTauAvg;


Did I modify it correctly?

Regards,
Sunxing
Sunxing is offline   Reply With Quote

Old   January 5, 2013, 04:03
Default
  #73
Member
 
Dan Kokron
Join Date: Dec 2012
Posts: 33
Rep Power: 14
dkokron is on a distinguished road
Daniel,

Is it safe to run your yPlus utility on a case that is decomposed for parallel execution?

Dan
dkokron is offline   Reply With Quote

Old   January 5, 2013, 08:48
Default
  #74
New Member
 
Abdulhaq
Join Date: Oct 2012
Posts: 27
Rep Power: 14
aljazari is on a distinguished road
Quote:
Originally Posted by dkokron View Post
Daniel,

Is it safe to run your yPlus utility on a case that is decomposed for parallel execution?

Dan
If I may answer, yes you can run it on a parallel case. If you have doubts, run it on a parallel case then reconstruct the latest time and run it again to check if you obtain similar results

1-
Code:
mpirun -np 8 yPlus -parallel
2-
Code:
reconstructPar -latestTime
yPlus -latestTime
aljazari is offline   Reply With Quote

Old   January 5, 2013, 12:23
Default
  #75
Member
 
Dan Kokron
Join Date: Dec 2012
Posts: 33
Rep Power: 14
dkokron is on a distinguished road
Abdulhaq,

Okay, I just wanted to make sure it was sane to even attempt before I dig into the code. I am getting a failure when running yPlus with 8 processes on the simpleFoam output from the incompressible/pisoFoam/les/motorBike tutorial.

Build : 2.1.x-45bb1ecfe0c7
Exec : yPlus -parallel -latestTime
Date : Jan 04 2013
Time : 23:57:01
Host : "#####"
PID : 42369
Case : /OpenFOAM/dkokron-2.1.x/tutorials/incompressible/pisoFoam/les/motorBike/motorBike_8p_scotch12b
nProcs : 8
Slaves : 7

Runs fin until ..................

Reading/calculating face flux field phi

Selecting incompressible transport model Newtonian
[0] #0 Foam::error:rintStack(Foam::Ostream&)--------------------------------------------------------------------------
An MPI process has executed an operation involving a call to the
"fork()" system call to create a child process. Open MPI is currently
operating in a condition that could result in memory corruption or
other system errors; your MPI job may hang, crash, or produce silent
data corruption. The use of fork() (or system() or other calls that
create child processes) is strongly discouraged.
dkokron is offline   Reply With Quote

Old   January 6, 2013, 10:42
Default
  #76
New Member
 
Abdulhaq
Join Date: Oct 2012
Posts: 27
Rep Power: 14
aljazari is on a distinguished road
Quote:
Originally Posted by dkokron View Post
Abdulhaq,

Okay, I just wanted to make sure it was sane to even attempt before I dig into the code. I am getting a failure when running yPlus with 8 processes on the simpleFoam output from the incompressible/pisoFoam/les/motorBike tutorial.

Build : 2.1.x-45bb1ecfe0c7
Exec : yPlus -parallel -latestTime
Date : Jan 04 2013
Time : 23:57:01
Host : "#####"
PID : 42369
Case : /OpenFOAM/dkokron-2.1.x/tutorials/incompressible/pisoFoam/les/motorBike/motorBike_8p_scotch12b
nProcs : 8
Slaves : 7

Runs fin until ..................

Reading/calculating face flux field phi

Selecting incompressible transport model Newtonian
[0] #0 Foam::error:rintStack(Foam::Ostream&)--------------------------------------------------------------------------
An MPI process has executed an operation involving a call to the
"fork()" system call to create a child process. Open MPI is currently
operating in a condition that could result in memory corruption or
other system errors; your MPI job may hang, crash, or produce silent
data corruption. The use of fork() (or system() or other calls that
create child processes) is strongly discouraged.
What happens if you reconstructed it and ran it on a single processor?
aljazari is offline   Reply With Quote

Old   January 6, 2013, 16:14
Default
  #77
Member
 
Dan Kokron
Join Date: Dec 2012
Posts: 33
Rep Power: 14
dkokron is on a distinguished road
Abdulhaq,

reconstructPar actually failed too.

I'm suspicious that already having output from yPlusRAS is getting in the way. I'll check this and get back.

Dan
dkokron is offline   Reply With Quote

Old   January 6, 2013, 17:10
Default
  #78
Member
 
Dan Kokron
Join Date: Dec 2012
Posts: 33
Rep Power: 14
dkokron is on a distinguished road
Nope, having output from yPlusRAS already in the case directory isn't the cause. I was able to reproduce the failure of yPlus using the tutorials/incompressible/simpleFoam/motorBike tutorial. The only modification was to add a the following to the Allrun script.

runApplication yPlus -latestTime

This is a 1 process case, so that removes parallelizm as the cause.

Can anyone reproduce this failure in their build? Here is the stack trace.

Program received signal SIGSEGV
raise () in /lib64/libc-2.11.1.so
(idb) where
#0 0x00002aaaafc58945 in raise () in /lib64/libc-2.11.1.so
#1 0x00002aaaaeed258a in _ZN4Foam7sigSegv10sigHandlerEi () in /OpenFOAM/OpenFOAM-2.1.x/platforms/linux64IccDPOpt/lib/libOpenFOAM.so
#2 0x00002aaaafc589e0 in __restore_rt () in /lib64/libc-2.11.1.so
#3 0x00002aaaaea7b5eb in _ZN4Foam7IOerrorclEPKcS2_iRKNS_8IOstreamE () in /OpenFOAM/OpenFOAM-2.1.x/platforms/linux64IccDPOpt/lib/libOpenFOAM.so
#4 0x00002aaaaeb69efd in _ZN4Foam11regIOobject10readStreamERKNS_4wordE () in /OpenFOAM/OpenFOAM-2.1.x/platforms/linux64IccDPOpt/lib/libOpenFOAM.so
#5 0x00002aaaaeb6411e in _ZN4Foam12IOdictionary8readFileEb () in /nobackupp1/dkokron/Projects/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64IccDPOpt/lib/libOpenFOAM.so
#6 0x00002aaaaeb634e3 in _ZN4Foam12IOdictionaryC1ERKNS_8IOobjectE () in /OpenFOAM/OpenFOAM-2.1.x/platforms/linux64IccDPOpt/lib/libOpenFOAM.so
#7 0x00002aaaab999ebe in _ZN4Foam14incompressible15turbulenceModel3NewERKNS _14GeometricFieldINS_6VectorIdEENS_12fvPatchFieldE NS_7volMeshEEERKNS2_IdNS_13fvsPatchFieldENS_11surf aceMeshEEERNS_14transportModelERKNS_4wordE () in /OpenFOAM/OpenFOAM-2.1.x/platforms/linux64IccDPOpt/lib/libincompressibleTurbulenceModel.so
#8 0x0000000000413a49 in main () in OpenFOAM/OpenFOAM-2.1.x/platforms/linux64IccDPOpt/bin/yPlus


Dan
dkokron is offline   Reply With Quote

Old   January 7, 2013, 15:18
Default
  #79
New Member
 
Abdulhaq
Join Date: Oct 2012
Posts: 27
Rep Power: 14
aljazari is on a distinguished road
Quote:
Originally Posted by dkokron View Post
Nope, having output from yPlusRAS already in the case directory isn't the cause. I was able to reproduce the failure of yPlus using the tutorials/incompressible/simpleFoam/motorBike tutorial. The only modification was to add a the following to the Allrun script.

runApplication yPlus -latestTime

This is a 1 process case, so that removes parallelizm as the cause.

Can anyone reproduce this failure in their build? Here is the stack trace.

Program received signal SIGSEGV
raise () in /lib64/libc-2.11.1.so
(idb) where
#0 0x00002aaaafc58945 in raise () in /lib64/libc-2.11.1.so
#1 0x00002aaaaeed258a in _ZN4Foam7sigSegv10sigHandlerEi () in /OpenFOAM/OpenFOAM-2.1.x/platforms/linux64IccDPOpt/lib/libOpenFOAM.so
#2 0x00002aaaafc589e0 in __restore_rt () in /lib64/libc-2.11.1.so
#3 0x00002aaaaea7b5eb in _ZN4Foam7IOerrorclEPKcS2_iRKNS_8IOstreamE () in /OpenFOAM/OpenFOAM-2.1.x/platforms/linux64IccDPOpt/lib/libOpenFOAM.so
#4 0x00002aaaaeb69efd in _ZN4Foam11regIOobject10readStreamERKNS_4wordE () in /OpenFOAM/OpenFOAM-2.1.x/platforms/linux64IccDPOpt/lib/libOpenFOAM.so
#5 0x00002aaaaeb6411e in _ZN4Foam12IOdictionary8readFileEb () in /nobackupp1/dkokron/Projects/OpenFOAM/OpenFOAM-2.1.x/platforms/linux64IccDPOpt/lib/libOpenFOAM.so
#6 0x00002aaaaeb634e3 in _ZN4Foam12IOdictionaryC1ERKNS_8IOobjectE () in /OpenFOAM/OpenFOAM-2.1.x/platforms/linux64IccDPOpt/lib/libOpenFOAM.so
#7 0x00002aaaab999ebe in _ZN4Foam14incompressible15turbulenceModel3NewERKNS _14GeometricFieldINS_6VectorIdEENS_12fvPatchFieldE NS_7volMeshEEERKNS2_IdNS_13fvsPatchFieldENS_11surf aceMeshEEERNS_14transportModelERKNS_4wordE () in /OpenFOAM/OpenFOAM-2.1.x/platforms/linux64IccDPOpt/lib/libincompressibleTurbulenceModel.so
#8 0x0000000000413a49 in main () in OpenFOAM/OpenFOAM-2.1.x/platforms/linux64IccDPOpt/bin/yPlus


Dan
I just ran this case using Allrun, then ran yPlus -latestTime, it went fine.

It didn't work the first time because the file "turbulenceProperties" was not found. I have no idea how the simulation worked without it. By including this file, yPlus worked with no problems.

Code:
/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.1.1                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    location    "constant";
    object      turbulenceProperties;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

simulationType  RASModel;


// ************************************************************************* //
yPlus-latestTime

Code:
/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.1.1                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
Build  : 2.1.1-221db2718bbb
Exec   : yPlus -latestTime
Date   : Jan 07 2013
Time   : 17:12:58
Host   : "p150"
PID    : 11118
Case   : /home/abdulhaq/OpenFOAM/abdulhaq-2.1.1/run/tutorials/incompressible/simpleFoam/motorBike
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Disallowing user-supplied system call operations

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create mesh for time = 500

Time = 500
Calculating wall distance

Writing wall distance to field y

Reading field U

Reading/calculating face flux field phi

Selecting incompressible transport model Newtonian
Selecting turbulence model type RASModel
Selecting RAS turbulence model kOmegaSST
kOmegaSSTCoeffs
{
    alphaK1         0.85034;
    alphaK2         1;
    alphaOmega1     0.5;
    alphaOmega2     0.85616;
    gamma1          0.5532;
    gamma2          0.4403;
    beta1           0.075;
    beta2           0.0828;
    betaStar        0.09;
    a1              0.31;
    c1              10;
}

Patch 3 named lowerWall
 d  : min: 0.000397575 max: 0.112139 average: 0.0112107
 y+ : min: 3.1981 max: 1512.07 average: 120.613

Patch 5 named motorBike_frt-fairing:001%1
 d  : min: 0.000268932 max: 0.0140758 average: 0.00364481
 y+ : min: 1.38902 max: 742.058 average: 128.173

Patch 6 named motorBike_windshield:002%2
 d  : min: 0.000430623 max: 0.00514934 average: 0.00220496
 y+ : min: 6.71703 max: 181.397 average: 94.8434

Patch 7 named motorBike_rr-wh-rim:005%5
 d  : min: 0.000326252 max: 0.0118033 average: 0.00415699
 y+ : min: 2.63874 max: 500.176 average: 115.064

Patch 8 named motorBike_rr-wh-rim:010%10
 d  : min: 0.000327294 max: 0.0115265 average: 0.00661764
 y+ : min: 2.88247 max: 599.602 average: 190.942

Patch 9 named motorBike_fr-wh-rim:011%11
 d  : min: 0.0023297 max: 0.0119296 average: 0.00677875
 y+ : min: 36.5513 max: 545 average: 223.732

Patch 10 named motorBike_fr-wh-brake-disk:012%12
 d  : min: 0.00284218 max: 0.0120569 average: 0.00762753
 y+ : min: 28.1762 max: 504.851 average: 243.077

Patch 11 named motorBike_frame:016-shadow%13
 d  : min: 0.000390573 max: 0.014721 average: 0.00620935
 y+ : min: 6.03975 max: 680.548 average: 223.03

Patch 12 named motorBike_rear-susp:014%14
 d  : min: 0.000280177 max: 0.012389 average: 0.00399244
 y+ : min: 1.01297 max: 528.538 average: 59.9596

Patch 13 named motorBike_rear-susp:014-shadow%15
 d  : min: 0.000306696 max: 0.0128792 average: 0.00373295
 y+ : min: 1.12489 max: 417.635 average: 57.7303

Patch 14 named motorBike_frame:016%16
 d  : min: 0.000362031 max: 0.0151413 average: 0.00555179
 y+ : min: 5.3446 max: 648.579 average: 202.398

Patch 15 named motorBike_rr-wh-rim:005-shadow%17
 d  : min: 0.000328156 max: 0.0118508 average: 0.0043253
 y+ : min: 3.76975 max: 495.027 average: 120.949

Patch 16 named motorBike_rr-wh-chain-hub:022%22
 d  : min: 0.0026924 max: 0.0112263 average: 0.00619253
 y+ : min: 33.3899 max: 528.787 average: 169.589

Patch 17 named motorBike_rearseat%24
 d  : min: 0.000355204 max: 0.0129121 average: 0.00357538
 y+ : min: 4.37037 max: 592.86 average: 125.88

Patch 18 named motorBike_frt-fairing%25
 d  : min: 0.000288805 max: 0.0124963 average: 0.0032173
 y+ : min: 1.6081 max: 559.374 average: 120.952

Patch 19 named motorBike_windshield%26
 d  : min: 0.000327636 max: 0.0121903 average: 0.00318627
 y+ : min: 2.67853 max: 491.891 average: 109.613

Patch 20 named motorBike_headlights%27
 d  : min: 0.00033471 max: 0.00865439 average: 0.00208301
 y+ : min: 4.35855 max: 252.623 average: 85.2525

Patch 21 named motorBike_driversseat%28
 d  : min: 0.000305473 max: 0.0124167 average: 0.00515002
 y+ : min: 0.914448 max: 386.061 average: 108.287

Patch 22 named motorBike_rear-body%29
 d  : min: 0.000292051 max: 0.0138328 average: 0.00451681
 y+ : min: 1.30529 max: 668.102 average: 109.362

Patch 23 named motorBike_fuel-tank%30
 d  : min: 0.000271334 max: 0.0136267 average: 0.00280589
 y+ : min: 1.14045 max: 348.381 average: 61.4069

Patch 24 named motorBike_exhaust%31
 d  : min: 0.000274933 max: 0.0146777 average: 0.00556152
 y+ : min: 1.52015 max: 693.262 average: 142.121

Patch 25 named motorBike_rr-wh-rim%32
 d  : min: 0.00030843 max: 0.0137104 average: 0.00350059
 y+ : min: 0.650911 max: 510.374 average: 78.4556

Patch 26 named motorBike_fr-mud-guard%33
 d  : min: 0.000303715 max: 0.0147352 average: 0.0047335
 y+ : min: 2.22657 max: 681.515 average: 144.748

Patch 27 named motorBike_fr-wh-rim%34
 d  : min: 0.00207944 max: 0.0125446 average: 0.00642802
 y+ : min: 13.0476 max: 553.624 average: 184.083

Patch 28 named motorBike_fr-wh-brake-disk%35
 d  : min: 0.000289026 max: 0.0130605 average: 0.00431061
 y+ : min: 1.66993 max: 600.062 average: 132.527

Patch 29 named motorBike_fr-brake-caliper%36
 d  : min: 0.00222925 max: 0.0133883 average: 0.00812573
 y+ : min: 23.8092 max: 525.755 average: 230.996

Patch 30 named motorBike_fr-wh-tyre%37
 d  : min: 0.000351147 max: 0.0128806 average: 0.00646112
 y+ : min: 3.27931 max: 699.639 average: 224.714

Patch 31 named motorBike_hbars%38
 d  : min: 0.00224132 max: 0.0150466 average: 0.00752225
 y+ : min: 19.7341 max: 794.713 average: 259.965

Patch 32 named motorBike_fr-forks%39
 d  : min: 0.00210806 max: 0.0161755 average: 0.00715355
 y+ : min: 11.4949 max: 719.727 average: 256.833

Patch 33 named motorBike_chain%40
 d  : min: 0.000361469 max: 0.016935 average: 0.00941252
 y+ : min: 6.89163 max: 705.813 average: 264.784

Patch 34 named motorBike_rr-wh-tyre%41
 d  : min: 0.000326759 max: 0.0145594 average: 0.00315657
 y+ : min: 1.84505 max: 714.16 average: 100.195

Patch 35 named motorBike_square-dial%42
 d  : min: 0.00394236 max: 0.0087235 average: 0.0063372
 y+ : min: 127.826 max: 257.126 average: 192.043

Patch 36 named motorBike_round-dial%43
 d  : min: 0.00475714 max: 0.0119439 average: 0.00785584
 y+ : min: 142.534 max: 345.877 average: 229.736

Patch 37 named motorBike_dial-holder%44
 d  : min: 0.00256609 max: 0.016704 average: 0.00648968
 y+ : min: 9.52715 max: 450.862 average: 123.603

Patch 38 named motorBike_rear-susp%45
 d  : min: 0.000282067 max: 0.0149661 average: 0.00578968
 y+ : min: 0.928264 max: 672.979 average: 161.333

Patch 39 named motorBike_rear-brake-lights%46
 d  : min: 0.00220472 max: 0.0152285 average: 0.00747612
 y+ : min: 8.98794 max: 213.055 average: 92.0237

Patch 40 named motorBike_rear-light-bracket%47
 d  : min: 0.000347537 max: 0.0117151 average: 0.00487275
 y+ : min: 2.40855 max: 269.656 average: 72.6323

Patch 41 named motorBike_frame%48
 d  : min: 0.000277013 max: 0.0141278 average: 0.00438018
 y+ : min: 0.86809 max: 700.369 average: 123.249

Patch 42 named motorBike_rear-mud-guard%49
 d  : min: 0.000320764 max: 0.0136426 average: 0.00384965
 y+ : min: 2.19329 max: 557.792 average: 105.802

Patch 43 named motorBike_rear-susp-spring-damp%50
 d  : min: 0.00198034 max: 0.0125607 average: 0.00641261
 y+ : min: 5.82302 max: 329.074 average: 85.663

Patch 44 named motorBike_fairing-inner-plate%51
 d  : min: 0.000285014 max: 0.0120323 average: 0.0041921
 y+ : min: 2.07677 max: 476.768 average: 97.1625

Patch 45 named motorBike_clutch-housing%52
 d  : min: 0.000271344 max: 0.0130613 average: 0.00312992
 y+ : min: 0.830421 max: 611.69 average: 79.8648

Patch 46 named motorBike_radiator%53
 d  : min: 0.00374399 max: 0.0202953 average: 0.0104274
 y+ : min: 55.7322 max: 709.455 average: 377.254

Patch 47 named motorBike_water-pipe%54
 d  : min: 0.00159212 max: 0.0143308 average: 0.00883537
 y+ : min: 32.0802 max: 762.355 average: 299.535

Patch 48 named motorBike_water-pump%55
 d  : min: 0.00297517 max: 0.012476 average: 0.00692214
 y+ : min: 4.92335 max: 262.53 average: 108.687

Patch 49 named motorBike_engine%56
 d  : min: 0.000306856 max: 0.0148045 average: 0.00332984
 y+ : min: 0.82158 max: 772.277 average: 89.4703

Patch 50 named motorBike_rear-shock-link%57
 d  : min: 0.00419348 max: 0.0141306 average: 0.00752708
 y+ : min: 20.9748 max: 254.993 average: 123.762

Patch 51 named motorBike_rear-brake-fluid-pot-bracket%58
 d  : min: 0.000339113 max: 0.0124309 average: 0.00418814
 y+ : min: 2.10608 max: 373.013 average: 65.2629

Patch 52 named motorBike_rear-brake-fluid-pot%59
 d  : min: 0.00273476 max: 0.0141282 average: 0.0083839
 y+ : min: 36.6767 max: 200.626 average: 123.659

Patch 53 named motorBike_footpeg%60
 d  : min: 0.00250189 max: 0.0133891 average: 0.00802141
 y+ : min: 19.6311 max: 670.721 average: 284.097

Patch 54 named motorBike_rr-wh-chain-hub%61
 d  : min: 0.000360962 max: 0.0115398 average: 0.005402
 y+ : min: 3.1775 max: 524.782 average: 146.742

Patch 55 named motorBike_rear-brake-caliper%62
 d  : min: 0.000357519 max: 0.013734 average: 0.00643976
 y+ : min: 5.35867 max: 768.159 average: 294.679

Patch 56 named motorBike_rider-helmet%65
 d  : min: 0.000303999 max: 0.0142727 average: 0.00289369
 y+ : min: 2.85383 max: 744.973 average: 124.842

Patch 57 named motorBike_rider-visor%66
 d  : min: 0.000910223 max: 0.00908061 average: 0.00246434
 y+ : min: 46.2841 max: 365.14 average: 158.312

Patch 58 named motorBike_rider-boots%67
 d  : min: 0.000271891 max: 0.01308 average: 0.00379892
 y+ : min: 1.56582 max: 648.84 average: 154.276

Patch 59 named motorBike_rider-gloves%68
 d  : min: 0.00228877 max: 0.0150876 average: 0.00731412
 y+ : min: 23.4853 max: 618.235 average: 208.937

Patch 60 named motorBike_rider-body%69
 d  : min: 0.000276985 max: 0.0129954 average: 0.00327886
 y+ : min: 0.846533 max: 790.118 average: 102.37

Patch 61 named motorBike_frame:0%70
 d  : min: 0.000351371 max: 0.0105157 average: 0.00345257
 y+ : min: 2.68923 max: 344.934 average: 102.693

Patch 62 named motorBike_frt-fairing:001-shadow%74
 d  : min: 0.000277253 max: 0.0139939 average: 0.00357307
 y+ : min: 1.44804 max: 662.019 average: 122.051

Patch 63 named motorBike_windshield-shadow%75
 d  : min: 0.000342722 max: 0.00956641 average: 0.00318237
 y+ : min: 3.87715 max: 376.309 average: 104.55

Patch 64 named motorBike_fr-mud-guard-shadow%81
 d  : min: 0.000328392 max: 0.012176 average: 0.00430565
 y+ : min: 2.78677 max: 518.462 average: 126.926

Patch 65 named motorBike_fr-wh-brake-disk-shadow%83
 d  : min: 0.000285772 max: 0.0144996 average: 0.00413188
 y+ : min: 2.25531 max: 699.929 average: 135.181

Patch 66 named motorBike_rear-mud-guard-shadow%84
 d  : min: 0.000295653 max: 0.013536 average: 0.00399838
 y+ : min: 2.03781 max: 510.218 average: 104.402

Patch 67 named motorBike_rear-susp-spring-damp-shadow%85
 d  : min: 0.00327178 max: 0.0123216 average: 0.00651848
 y+ : min: 8.7858 max: 331.149 average: 80.401

Patch 68 named motorBike_radiator-shadow%86
 d  : min: 0.00376333 max: 0.0202953 average: 0.0111452
 y+ : min: 51.8406 max: 764.851 average: 436.34

Patch 69 named motorBike_rear-shock-link-shadow%87
 d  : min: 0.00462519 max: 0.0119373 average: 0.0082026
 y+ : min: 75.663 max: 242.782 average: 141.328

Patch 70 named motorBike_rear-brake-fluid-pot-bracket-shadow%88
 d  : min: 0.000350817 max: 0.0120304 average: 0.00447425
 y+ : min: 2.50756 max: 201.903 average: 69.0392

Patch 71 named motorBike_rr-wh-chain-hub-shadow%89
 d  : min: 0.000430666 max: 0.0102251 average: 0.00555544
 y+ : min: 6.13939 max: 479.068 average: 147.676

Writing yPlus to field yPlus

End
aljazari is offline   Reply With Quote

Old   January 8, 2013, 22:36
Default
  #80
Member
 
Dan Kokron
Join Date: Dec 2012
Posts: 33
Rep Power: 14
dkokron is on a distinguished road
Abdulhaq,

That's it! I was missing constant/turbulenceProperties. Wish the code would have been a little more forthcoming about this. Anyway, problem solved.

Dan
dkokron is offline   Reply With Quote

Reply

Tags
low-re rans, y* value, y+ value


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
Low values of Lift force vmlxb6 CFX 1 February 2, 2011 06:13
Incompressible Turbulence models achinta OpenFOAM 4 May 27, 2010 11:35
turbulence models? haider FLUENT 0 March 8, 2006 00:58
Turbulence Models and external flow. Alan FLUENT 3 November 22, 2005 05:46
Turbulence boundary values lego CFX 9 October 25, 2002 12:55


All times are GMT -4. The time now is 15:40.