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 Members List Search Today's Posts Mark Forums Read

Like Tree81Likes

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   February 14, 2013, 09:28
Default
  #81
New Member
 
Sha Huang
Join Date: Dec 2012
Posts: 22
Rep Power: 13
Joanna Huang is on a distinguished road
Quote:
Originally Posted by aljazari View Post
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


Hi:

I saw your results of yPlus, the value seems very big, and how can we judg the y+ is in the proper range?
Joanna Huang is offline   Reply With Quote

Old   April 5, 2013, 13:48
Default
  #82
Member
 
Dan Kokron
Join Date: Dec 2012
Posts: 33
Rep Power: 13
dkokron is on a distinguished road
@post 71.

I posting here as I have not found any other source for the yPlus utility from Daniel WEI.

I have attached a patch to allow yPlus.C to be compiled under OpenFOAM-2.2.x. Runs fine in my environment. Please check my work.
Attached Files
File Type: patch yPlus.C.patch (1.1 KB, 74 views)
dkokron is offline   Reply With Quote

Old   May 27, 2013, 09:33
Default
  #83
New Member
 
bscphil's Avatar
 
Phillip
Join Date: Mar 2012
Location: Germany
Posts: 27
Rep Power: 14
bscphil is on a distinguished road
Quote:
Originally Posted by romant View Post
hej,

unfortunately, one cannot compile this one with OF 2.0.0. It throws the following error message:

Code:
plusPostRANS.C: In function ‘int main(int, char**)’:
plusPostRANS.C:96:40: error: ‘class Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> >’ has no member named ‘boundaryField’
plusPostRANS.C:99:36: error: ‘class Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> >’ has no member named ‘boundaryField’
apparently, the class geometric field was changed.

Hi foamers,

i have change the source code of the "plusPostRANS" tool for the usage in OpenFOAM 2.2.0 (see the attached file)

Just compile it and it works

@romant: i hope this helps..
Attached Files
File Type: gz plusPostRANS.tar.gz (97.4 KB, 115 views)
romant, Uyan and benqing like this.
bscphil is offline   Reply With Quote

Old   May 28, 2013, 12:22
Default
  #84
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 dkokron View Post
@post 71.

I posting here as I have not found any other source for the yPlus utility from Daniel WEI.

I have attached a patch to allow yPlus.C to be compiled under OpenFOAM-2.2.x. Runs fine in my environment. Please check my work.

Sorry, I was too busy. Thanks a lot for your contribution.
__________________
~
Daniel WEI
-------------
Boeing Research & Technology - China
Beijing, China
Email
lakeat is offline   Reply With Quote

Old   May 29, 2013, 07:47
Default
  #85
New Member
 
bscphil's Avatar
 
Phillip
Join Date: Mar 2012
Location: Germany
Posts: 27
Rep Power: 14
bscphil is on a distinguished road
Quote:
Originally Posted by lakeat View Post
Hi again,

I think I have fixed the bug.

yPlus utility: Attachment 17317
yStar utility: Attachment 17318

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:


Hi foamers,

firstly thanks for the nice tool. I want to compile the yPlus.C under OpenFoam 2.2.0 and have some errors due to a new header file: they change the name from "basicThermo.H" to "fluidThermo.H" (OpenFoam 2.1.x --> OpenFoam 2.2.x). This is why i change all entries in the source code of "yPlus.C" and it works fine, now. I have attached the new files for OpenFoam 2.2.x:
bscphil is offline   Reply With Quote

Old   May 29, 2013, 07:56
Default
  #86
New Member
 
bscphil's Avatar
 
Phillip
Join Date: Mar 2012
Location: Germany
Posts: 27
Rep Power: 14
bscphil is on a distinguished road
Quote:
Originally Posted by bscphil View Post
Hi foamers,

firstly thanks for the nice tool. I want to compile the yPlus.C under OpenFoam 2.2.0 and have some errors due to a new header file: they change the name from "basicThermo.H" to "fluidThermo.H" (OpenFoam 2.1.x --> OpenFoam 2.2.x). This is why i change all entries in the source code of "yPlus.C" and it works fine, now. I have attached the new files for OpenFoam 2.2.x:
ok, sorry, the file::
Attached Files
File Type: gz yPlus_openfoam220.tar.gz (2.3 KB, 94 views)
immortality, edooo0 and benqing like this.
bscphil is offline   Reply With Quote

Old   June 5, 2013, 17:44
Default
  #87
Senior Member
 
immortality's Avatar
 
Ehsan
Join Date: Oct 2012
Location: Iran
Posts: 2,208
Rep Power: 26
immortality is on a distinguished road
Hi
is this useful tool appropriate for low-Re turbulent models only or is general?
also I have found another tool for yPlus calculating.
could someone have a look into it to see if its correct too and what are the differences and which is more suitable for compressible unsteady cases?
and also could someone suggest some turbulent methods for compressible unsteady internal flow cases with low-Re mesh(y+<=1) and with wall function(30<=y+<60 or 100 in some unsteady cases)
thanks.
----------------------------------------------------------
complement:
in the file I have attached the fomula of yPlus is:
Code:
yPlus.boundaryField()[patchi] =
	        y[patchi]
	      * sqrt
	        (
		   mut.boundaryField()[patchi]
		 * mag(U.boundaryField()[patchi].snGrad())
		 / rho.boundaryField()[patchi]
                )
	      / (RASModel->mu().boundaryField()[patchi]/rho.boundaryField()[patchi]);
and in previous:
Code:
yPlus.boundaryField()[patchi] =
                        d[patchi]
                       *sqrt
                        (
                            rho.boundaryField()[patchi]
                           *muEff.boundaryField()[patchi]
                           *mag(U.boundaryField()[patchi].snGrad())
                        )
                       /muLam.boundaryField()[patchi];
the difference look to be only in mu's.could someone give more explanation on this subject?
(if there are other differences anyone notices please let me know)
Attached Files
File Type: gz yPlus220.tar.gz (97.1 KB, 30 views)
__________________
Injustice Anywhere is a Threat for Justice Everywhere.Martin Luther King.
To Be or Not To Be,Thats the Question!
The Only Stupid Question Is the One that Goes Unasked.
immortality is offline   Reply With Quote

Old   July 19, 2013, 12:03
Default
  #88
Senior Member
 
Join Date: Mar 2010
Posts: 173
Rep Power: 17
Jonathan is on a distinguished road
Hi Dan,

just a heads up to say thanks for the x2 tools - worked great for me and i thought well distributed functionality (i.e. in terms of separating yStar / yPlus / refinedWallLayer etc).

much appreciated
Jon

Quote:
Originally Posted by lakeat View Post
Sorry, I was too busy. Thanks a lot for your contribution.
Jonathan is offline   Reply With Quote

Old   August 5, 2013, 03:48
Default
  #89
Member
 
Glenn Carlson, PE, PhD (ret)
Join Date: Oct 2012
Location: US
Posts: 49
Rep Power: 13
gcengineer is on a distinguished road
I'm having trouble building yPlus with wmake. I have no problem running OF 2.2.1 in Ubuntu 12.04 64-bit.

I extracted yPlus.c and the Make folder to /opt/openfoam221/src. It seems that, according to

http://www.openfoam.org/docs/user/co...0-720003.2.2.1,

all I need to do is run wmake (no arguments, no directory) from /opt/openfoam221/src/yPlus. But when I run wmake, I get the following errors.

gac@gac-vm01:/opt/openfoam221/src/yPlus$ wmake
mkdir: cannot create directory `linuxGccDPOpt': Permission denied
/bin/sh: 1: cannot create linuxGccDPOpt/options: Directory nonexistent
make: *** [linuxGccDPOpt/options] Error 2
/opt/openfoam221/wmake/MakefileFiles:40: linuxGccDPOpt/options: No such file or directory
make: *** No rule to make target `linuxGccDPOpt/options'. Stop.
wmake error: file 'Make/linuxGccDPOpt/objectFiles' could not be created in /opt/openfoam221/src/yPlus

Thanks for any assistance.

Regards,
Glenn
gcengineer is offline   Reply With Quote

Old   August 5, 2013, 03:58
Default
  #90
Senior Member
 
Bernhard
Join Date: Sep 2009
Location: Delft
Posts: 790
Rep Power: 21
Bernhard is on a distinguished road
Try either "sudo wmake" if you have root privileges, or compile it from a directory where you do have write access, e.g. "cd $WM_PROJECT_USER_DIR"
Bernhard is offline   Reply With Quote

Old   August 5, 2013, 04:27
Default
  #91
Member
 
Glenn Carlson, PE, PhD (ret)
Join Date: Oct 2012
Location: US
Posts: 49
Rep Power: 13
gcengineer is on a distinguished road
Thanks for quick reply.

Extracting yPlus.c and Make folder to $WM_PROJECT_USER_DIR/yPlus helped (proper permissions), but I get another permission error when I run wmake from there:

gac@gac-vm01:~/OpenFOAM/gac-2.2.1/yPlus$ wmake
Making dependency list for source file yPlus.C
SOURCE=yPlus.C ; g++ -m32 -Dlinux -DWM_DP -Wall -Wextra -Wno-unused-parameter -Wold-style-cast -O3 -DNoRepository -ftemplate-depth-100 -I/opt/openfoam221/src/meshTools/lnInclude -I/opt/openfoam221/src/turbulenceModels -I/opt/openfoam221/src/transportModels -I/opt/openfoam221/src/thermophysicalModels/basic/lnInclude -I/opt/openfoam221/src/finiteVolume/lnInclude -IlnInclude -I. -I/opt/openfoam221/src/OpenFOAM/lnInclude -I/opt/openfoam221/src/OSspecific/POSIX/lnInclude -fPIC -c $SOURCE -o Make/linuxGccDPOpt/yPlus.o
g++ -m32 -Dlinux -DWM_DP -Wall -Wextra -Wno-unused-parameter -Wold-style-cast -O3 -DNoRepository -ftemplate-depth-100 -I/opt/openfoam221/src/meshTools/lnInclude -I/opt/openfoam221/src/turbulenceModels -I/opt/openfoam221/src/transportModels -I/opt/openfoam221/src/thermophysicalModels/basic/lnInclude -I/opt/openfoam221/src/finiteVolume/lnInclude -IlnInclude -I. -I/opt/openfoam221/src/OpenFOAM/lnInclude -I/opt/openfoam221/src/OSspecific/POSIX/lnInclude -fPIC -Xlinker --add-needed -Xlinker --no-as-needed Make/linuxGccDPOpt/yPlus.o -L/opt/openfoam221/platforms/linuxGccDPOpt/lib \
-lcompressibleTurbulenceModel -lcompressibleRASModels -lcompressibleLESModels -lincompressibleTurbulenceModel -lincompressibleTransportModels -lincompressibleRASModels -lincompressibleLESModels -lfiniteVolume -lgenericPatchFields -lOpenFOAM -ldl -lm -o /opt/openfoam221/platforms/linuxGccDPOpt/bin/yPlus
/usr/bin/ld: cannot open output file /opt/openfoam221/platforms/linuxGccDPOpt/bin/yPlus: Permission denied
collect2: ld returned 1 exit status

When I run 'sudo wmake', I get the error 'sudo: wmake: command not found'.

Perhaps next time I will not install OpenFOAM as root.

Thanks again.
gcengineer is offline   Reply With Quote

Old   August 5, 2013, 05:08
Default
  #92
Member
 
Glenn Carlson, PE, PhD (ret)
Join Date: Oct 2012
Location: US
Posts: 49
Rep Power: 13
gcengineer is on a distinguished road
I successfully built yPlus by changing 'FOAM_APPBIN' in Make/files to 'FOAM_USER_APPBIN'.
Miguel Ángel likes this.
gcengineer is offline   Reply With Quote

Old   October 15, 2013, 18:12
Default
  #93
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
Hey guys,

1. I found calculate wall distance twice is unnecessary, so I disable it.
2. Here are the two utilities doing the same job as before, with small changes. Tested on OpenFOAM-2.2.x
3. I have also added a few write-out control.

Please use "-help" to find more info.

yPlus: yPlus.tar.gz
yStar: yStar.tar.gz
__________________
~
Daniel WEI
-------------
Boeing Research & Technology - China
Beijing, China
Email
lakeat is offline   Reply With Quote

Old   November 18, 2013, 13:44
Default
  #94
Member
 
Jubayer
Join Date: Oct 2009
Location: The University of Western Ontario, London, Ontario
Posts: 42
Blog Entries: 1
Rep Power: 16
cm_jubayer is on a distinguished road
may be of interest to someone:

http://www.openfoam.org/mantisbt/view.php?id=835

Jubayer
immortality likes this.
cm_jubayer is offline   Reply With Quote

Old   November 24, 2013, 21:42
Default Ask for your help
  #95
New Member
 
borhan
Join Date: May 2011
Posts: 5
Rep Power: 14
borhanalhoseini is on a distinguished road
Dear all,

This issue was already discussed in some length in various threads. Therefore I finally reviewed and attached my plusPostRANS utility.

It calculates y+ and u+ values when using one of the available low-Re RANS turbulence models. I hope the code as well as the output will be self-explaining

I also attached a case. It is a straight pipe (wedge) with periodic inlet/outlet boundary conditions. The Reynolds number based on the mean axial velocity is Re=5300. You can run the case using pisoFoam.

I would appreciate any feedback, it might be still improved at some point!

Best regards,
Florian
borhanalhoseini is offline   Reply With Quote

Old   November 24, 2013, 21:50
Default Ask for your help
  #96
New Member
 
borhan
Join Date: May 2011
Posts: 5
Rep Power: 14
borhanalhoseini is on a distinguished road
Dear all,

This issue was already discussed in some length in various threads. Therefore I finally reviewed and attached my plusPostRANS utility.

It calculates y+ and u+ values when using one of the available low-Re RANS turbulence models. I hope the code as well as the output will be self-explaining

I also attached a case. It is a straight pipe (wedge) with periodic inlet/outlet boundary conditions. The Reynolds number based on the mean axial velocity is Re=5300. You can run the case using pisoFoam.

I would appreciate any feedback, it might be still improved at some point!

Best regards,
Florian


Dear Florian,
I guess you have been discussing about y+. I have a series of experimental data (u as velocity and y as normal distance+ u free stream+ characteristic length) and I would like to u* and then y+ as well as u+...I was wondering if you or any friend can help me in this respect.
borhanalhoseini is offline   Reply With Quote

Old   February 25, 2014, 17:14
Default
  #97
Senior Member
 
Join Date: Nov 2012
Location: Bavaria
Posts: 145
Rep Power: 13
aylalisa is on a distinguished road
Dear borhanalhoseini,

is there a working version of yPlusUplus (or plusPostRANS) for OF.2.2.2 for LES simulations (NWR) that calculates the fields - as originally intended by florian_krause - to be able to create plots u+ vs. y+?

Quote:
Within the two mentioned lines I calculate fields of y+ and u+. If you think about a fully developed turbulent channel or pipe flow you might want to compare your results to, for instance, available DNS data by means of the near wall velocity profile in wall units or in other words comparing a profile of u+ over y+.
#3

Thank you !

Aylalisa

Last edited by aylalisa; February 26, 2014 at 04:10.
aylalisa is offline   Reply With Quote

Old   April 28, 2014, 06:11
Default
  #98
Senior Member
 
ArathoN
Join Date: Jul 2011
Posts: 137
Rep Power: 15
ArathoN is on a distinguished road
The best way is to simulate the case then use the wallshearstress to calculate tau or the yplus (you'll have yplus and y so you can derive utau). At this point you have all the needed data to calculate uplus.
aylalisa likes this.
ArathoN is offline   Reply With Quote

Old   April 29, 2014, 04:48
Default
  #99
Senior Member
 
Join Date: Nov 2012
Location: Bavaria
Posts: 145
Rep Power: 13
aylalisa is on a distinguished road
Hello ArathoN,

documentation for wallShearStress utility says:
'Calculates and writes the wall shear stress, for the specified times when using RAS turbulence models. ' This description prevented me to use it.

The computation of stress is 'stress tensor * vector', that is what wallShearStress utility does. Why does the documentation insist on RAS turbulence models


Aylalisa
aylalisa is offline   Reply With Quote

Old   April 29, 2014, 07:09
Default
  #100
Senior Member
 
ArathoN
Join Date: Jul 2011
Posts: 137
Rep Power: 15
ArathoN is on a distinguished road
Quote:
Originally Posted by aylalisa View Post
Hello ArathoN,

documentation for wallShearStress utility says:
'Calculates and writes the wall shear stress, for the specified times when using RAS turbulence models. ' This description prevented me to use it.

The computation of stress is 'stress tensor * vector', that is what wallShearStress utility does. Why does the documentation insist on RAS turbulence models


Aylalisa

Never noticed that part, unfortunately my C++ knowledge is almost null so i can't delve into the appropriate code.

If you have the yPlus, y and the cinematic viscosity you can calculate the friction velocity (utau) and from it you can compute the uPlus.

My only problem here would be how to consider the utau, because I'm certain it is variable on the domain so it's better to averag it (in my opinion) but I'm not certain about this.

P.S Another doubt what section would you consider when you want to compare the simulated data with the experimental data for example in a channel?
ArathoN is offline   Reply With Quote

Reply

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

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
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 08:21.