Commit dd3fb622 authored by Andrew Heather's avatar Andrew Heather
Browse files

ENH: Tutorials - added new channel tutorial to show turbulence DFSEM inlet behaviour

parent 85490129
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: plus |
| \\ / A nd | Web: www.OpenFOAM.com |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volVectorField;
location "1";
object U;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [ 0 1 -1 0 0 0 0 ];
// Re = 395, L = pi, utau = 1, nu = pi/395 = 7.9534e-3
// Ub/utau = 17.55
internalField uniform ( 17.55 0 0 );
boundaryField
{
bottomWall
{
type fixedValue;
value uniform ( 0 0 0 );
}
topWall
{
type fixedValue;
value uniform ( 0 0 0 );
}
sides_half0
{
type cyclic;
}
sides_half1
{
type cyclic;
}
inlet
{
type turbulentDFSEMInlet;
delta 2;
nCellPerEddy 1;
mapMethod nearestCell;
value $internalField;
}
outlet
{
type inletOutlet;
inletValue uniform (0 0 0);
value uniform (0 0 0);
}
}
// ************************************************************************* //
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: plus |
| \\ / A nd | Web: www.OpenFOAM.com |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
location "1";
object k;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [ 0 2 -2 0 0 0 0 ];
internalField uniform 1e-5;
boundaryField
{
bottomWall
{
type fixedValue;
value uniform 0;
}
topWall
{
type fixedValue;
value uniform 0;
}
sides_half0
{
type cyclic;
}
sides_half1
{
type cyclic;
}
inlet
{
type fixedValue;
value $internalField;
}
outlet
{
type fixedValue;
value $internalField;
}
}
// ************************************************************************* //
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: plus |
| \\ / A nd | Web: www.OpenFOAM.com |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
location "1";
object nut;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [ 0 2 -1 0 0 0 0 ];
internalField uniform 0;
boundaryField
{
bottomWall
{
type zeroGradient;
}
topWall
{
type zeroGradient;
}
sides_half0
{
type cyclic;
}
sides_half1
{
type cyclic;
}
inlet
{
type calculated;
value uniform 1e-8;
}
outlet
{
type calculated;
value uniform 1e-8;
}
}
// ************************************************************************* //
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: plus |
| \\ / A nd | Web: www.OpenFOAM.com |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
location "1";
object p;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [ 0 2 -2 0 0 0 0 ];
internalField uniform 0;
boundaryField
{
bottomWall
{
type zeroGradient;
}
topWall
{
type zeroGradient;
}
sides_half0
{
type cyclic;
}
sides_half1
{
type cyclic;
}
inlet
{
type zeroGradient;
}
outlet
{
type fixedValue;
value uniform 0;
}
}
// ************************************************************************* //
#!/bin/sh
cd ${0%/*} || exit 1 # run from this directory
# Source tutorial run functions
. $WM_PROJECT_DIR/bin/tools/RunFunctions
# Get application directory
application=`getApplication`
runApplication blockMesh
#- Run serial
#runApplication $application
#- Run parallel
runApplication decomposePar -cellDist
runParallel $application
#runApplication reconstructPar -latestTime
# ----------------------------------------------------------------- end-of-file
Channel test case for ReTau=395,based on the reference:
Poletto, R., Craft, T., and Revell, A.,
"A New Divergence Free Synthetic Eddy Method for the
Reproduction of Inlet Flow Conditions for LES",
Flow Turbulence Combust (2013) 91:519-539
Inlet patch Reynolds stress, velocity and turbulence length scale data has been
extracted from DNS data of:
Moser, Kim & Mansour
"DNS of Turbulent Channel Flow up to Re_tau=590",
Physics of Fluids (1999) vol 11, 943-945.
Data available from (last checked 28 June 2016)
http://turbulence.ices.utexas.edu/MKM_1999.html
/*--------------------------------*- C++ -*----------------------------------*| ========= | |
| \ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \ / O peration | Version: plus |
| \ / A nd | Web: www.OpenFOAM.com |
| \/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class scalarAverageField;
object values;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
// LMean
0
(
0
1.76693e-06
1.06447e-07
2.85582e-08
1.73393e-08
3.56243e-08
3.79499e-08
1.50479e-07
1.55614e-07
1.93408e-07
4.50118e-08
2.6295e-10
6.18128e-08
4.4709e-08
5.98657e-07
2.25711e-06
4.20946e-06
6.292e-06
7.5414e-06
9.02814e-06
1.04737e-05
1.25187e-05
1.37559e-05
1.40338e-05
1.24681e-05
1.08272e-05
9.10144e-06
7.81661e-06
6.08627e-06
4.68729e-06
3.01241e-06
1.93855e-06
1.03945e-06
2.49601e-07
2.33919e-07
3.47823e-08
3.05604e-06
8.25413e-06
1.04748e-05
1.36651e-05
2.35252e-05
3.87371e-05
5.21582e-05
6.71737e-05
6.89399e-05
4.57512e-05
1.94891e-05
9.7778e-06
1.47858e-05
3.40409e-05
6.06823e-05
8.59513e-05
8.29571e-05
5.42278e-05
3.46015e-05
2.20432e-05
1.73878e-05
1.05749e-05
4.67814e-06
8.66793e-06
1.79361e-05
1.22464e-05
3.93103e-06
1.00778e-06
1.81983e-06
2.70815e-05
9.56468e-05
0.000139304
0.00012493
9.50459e-05
4.78307e-05
8.71962e-06
1.24275e-05
4.91383e-05
8.77264e-05
0.000114449
0.00014567
0.000201758
0.000251863
0.000272594
0.000213599
0.000145126
0.000115603
0.000122779
0.000116865
8.36843e-05
4.06009e-05
3.49149e-05
5.22978e-05
5.07525e-05
2.5995e-05
2.6683e-06
9.29144e-06
1.14821e-05
9.41939e-06
1.01946e-05
6.64024e-06
1.1913e-06
3.25066e-06
8.47834e-06
1.42023e-06
4.30742e-05
0.000106228
0.000115468
9.07632e-05
7.04511e-05
9.21776e-05
0.000104486
0.000111678
0.000109852
0.000111867
9.10747e-05
6.82206e-05
4.09085e-05
2.3961e-05
2.0467e-06
7.74345e-06
1.66716e-05
1.54967e-05
2.95089e-05
4.82299e-05
6.99781e-05
7.16947e-05
7.33475e-05
7.40551e-05
9.45846e-05
0.000107202
0.000120068
0.000122517
0.000120068
0.000107202
9.45846e-05
7.40551e-05
7.33475e-05
7.16947e-05
6.99781e-05
4.82299e-05
2.95089e-05
1.54967e-05
1.66716e-05
7.74345e-06
2.0467e-06
2.3961e-05
4.09085e-05
6.82206e-05
9.10747e-05
0.000111867
0.000109852
0.000111678
0.000104486
9.21776e-05
7.04511e-05
9.07632e-05
0.000115468
0.000106228
4.30742e-05
1.42023e-06
8.47834e-06
3.25066e-06
1.1913e-06
6.64024e-06
1.01946e-05
9.41939e-06
1.14821e-05
9.29144e-06
2.6683e-06
2.5995e-05
5.07525e-05
5.22978e-05
3.49149e-05
4.06009e-05
8.36843e-05
0.000116865
0.000122779
0.000115603
0.000145126
0.000213599
0.000272594
0.000251863
0.000201758
0.00014567
0.000114449
8.77264e-05
4.91383e-05
1.24275e-05
8.71962e-06
4.78307e-05
9.50459e-05
0.00012493
0.000139304
9.56468e-05
2.70815e-05
1.81983e-06
1.00778e-06
3.93103e-06
1.22464e-05
1.79361e-05
8.66793e-06
4.67814e-06
1.05749e-05
1.73878e-05
2.20432e-05
3.46015e-05
5.42278e-05
8.29571e-05
8.59513e-05
6.06823e-05
3.40409e-05
1.47858e-05
9.7778e-06
1.94891e-05
4.57512e-05
6.89399e-05
6.71737e-05
5.21582e-05
3.87371e-05
2.35252e-05
1.36651e-05
1.04748e-05
8.25413e-06
3.05604e-06
3.47823e-08
2.33919e-07
2.49601e-07
1.03945e-06
1.93855e-06
3.01241e-06
4.68729e-06
6.08627e-06
7.81661e-06
9.10144e-06
1.08272e-05
1.24681e-05
1.40338e-05
1.37559e-05
1.25187e-05
1.04737e-05
9.02814e-06
7.5414e-06
6.292e-06
4.20946e-06
2.25711e-06
5.98657e-07
4.4709e-08
6.18128e-08
2.6295e-10
4.50118e-08
1.93408e-07
1.55614e-07
1.50479e-07
3.79499e-08
3.56243e-08
1.73393e-08
2.85582e-08
1.06447e-07
1.76693e-06
0
)
/*--------------------------------*- C++ -*----------------------------------*| ========= | |
| \ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \ / O peration | Version: plus |
| \ / A nd | Web: www.OpenFOAM.com |
| \/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class symmTensorAverageField;
object values;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
// UU, UV, UZ, VV, VW, WW
// Average
(0 0 0 0 0 0)
(
(5.3376e-28 1.1206e-30 -6.1282e-30 4.9018e-29 6.2793e-32 5.9575e-28)
(1.3684e-04 -2.4414e-08 5.9464e-08 9.4734e-11 -2.8525e-11 5.2943e-05)
(2.1828e-03 -1.5728e-06 8.3228e-07 2.3603e-08 -1.7141e-09 8.2715e-04)
(1.1020e-02 -1.8113e-05 3.3616e-06 5.7878e-07 -1.9051e-08 4.0392e-03)
(3.4718e-02 -1.0324e-04 1.0732e-05 5.4400e-06 -1.0411e-07 1.2147e-02)
(8.4454e-02 -4.0027e-04 2.8297e-05 3.0014e-05 -3.8272e-07 2.7852e-02)
(1.7438e-01 -1.2143e-03 6.6865e-05 1.1756e-04 -1.0969e-06 5.3564e-02)
(3.2123e-01 -3.1009e-03 1.4506e-04 3.6204e-04 -2.6288e-06 9.0947e-02)
(5.4336e-01 -6.9486e-03 2.9097e-04 9.3203e-04 -5.4596e-06 1.4062e-01)
(8.5851e-01 -1.4015e-02 5.4723e-04 2.0879e-03 -1.0091e-05 2.0209e-01)
(1.2802e+00 -2.5858e-02 9.7373e-04 4.1852e-03 -1.7119e-05 2.7384e-01)
(1.8132e+00 -4.4135e-02 1.6529e-03 7.6610e-03 -2.7764e-05 3.5369e-01)
(2.4490e+00 -7.0276e-02 2.6847e-03 1.3005e-02 -4.4587e-05 4.3919e-01)
(3.1644e+00 -1.0513e-01 4.1721e-03 2.0724e-02 -7.1044e-05 5.2804e-01)
(3.9229e+00 -1.4869e-01 6.1889e-03 3.1300e-02 -1.0893e-04 6.1834e-01)
(4.6811e+00 -1.9999e-01 8.7562e-03 4.5157e-02 -1.5460e-04 7.0871e-01)
(5.3953e+00 -2.5724e-01 1.1834e-02 6.2636e-02 -1.9552e-04 7.9819e-01)
(6.0286e+00 -3.1808e-01 1.5387e-02 8.3971e-02 -2.1029e-04 8.8610e-01)
(6.5558e+00 -3.8001e-01 1.9449e-02 1.0928e-01 -1.7262e-04 9.7178e-01)
(6.9642e+00 -4.4073e-01 2.4133e-02 1.3854e-01 -6.0490e-05 1.0545e+00)
(7.2524e+00 -4.9837e-01 2.9501e-02 1.7162e-01 1.3268e-04 1.1333e+00)
(7.4274e+00 -5.5162e-01 3.5418e-02 2.0824e-01 3.9086e-04 1.2075e+00)
(7.5016e+00 -5.9969e-01 4.1511e-02 2.4803e-01 6.7818e-04 1.2762e+00)
(7.4897e+00 -6.4224e-01 4.7330e-02 2.9050e-01 9.5511e-04 1.3387e+00)
(7.4069e+00 -6.7924e-01 5.2512e-02 3.3510e-01 1.1989e-03 1.3948e+00)
(7.2679e+00 -7.1090e-01 5.6926e-02 3.8121e-01 1.4016e-03 1.4444e+00)
(7.0865e+00 -7.3759e-01 6.0659e-02 4.2818e-01 1.5532e-03 1.4878e+00)
(6.8750e+00 -7.5980e-01 6.3835e-02 4.7535e-01 1.6360e-03 1.5251e+00)
(6.6441e+00 -7.7810e-01 6.6418e-02 5.2211e-01 1.6465e-03 1.5568e+00)
(6.4025e+00 -7.9302e-01 6.8207e-02 5.6790e-01 1.6026e-03 1.5833e+00)
(6.1573e+00 -8.0503e-01 6.9001e-02 6.1229e-01 1.5073e-03 1.6047e+00)
(5.9141e+00 -8.1452e-01 6.8896e-02 6.5493e-01 1.3268e-03 1.6216e+00)
(5.6767e+00 -8.2184e-01 6.8285e-02 6.9554e-01 1.0303e-03 1.6344e+00)
(5.4477e+00 -8.2733e-01 6.7481e-02 7.3388e-01 6.4119e-04 1.6439e+00)
(5.2291e+00 -8.3133e-01 6.6450e-02 7.6971e-01 2.3642e-04 1.6512e+00)
(5.0223e+00 -8.3412e-01 6.4937e-02 8.0285e-01 -1.3612e-04 1.6567e+00)
(4.8279e+00 -8.3590e-01 6.2875e-02 8.3316e-01 -4.7569e-04 1.6606e+00)
(4.6462e+00 -8.3683e-01 6.0710e-02 8.6061e-01 -7.8412e-04 1.6623e+00)
(4.4768e+00 -8.3700e-01 5.9182e-02 8.8520e-01 -1.0539e-03 1.6618e+00)
(4.3194e+00 -8.3649e-01 5.8527e-02 9.0700e-01 -1.2655e-03 1.6594e+00)
(4.1730e+00 -8.3523e-01 5.8437e-02 9.2607e-01 -1.4126e-03 1.6553e+00)
(4.0370e+00 -8.3321e-01 5.8512e-02 9.4247e-01 -1.4930e-03 1.6500e+00)
(3.9106e+00 -8.3044e-01 5.8376e-02 9.5626e-01 -1.5118e-03 1.6439e+00)
(3.7930e+00 -8.2687e-01 5.8084e-02 9.6754e-01 -1.5196e-03 1.6373e+00)
(3.6834e+00 -8.2252e-01 5.8125e-02 9.7651e-01 -1.6228e-03 1.6299e+00)
(3.5821e+00 -8.1762e-01 5.8537e-02 9.8341e-01 -1.8295e-03 1.6209e+00)
(3.4896e+00 -8.1249e-01 5.8809e-02 9.8850e-01 -2.0145e-03 1.6105e+00)
(3.4053e+00 -8.0726e-01 5.8494e-02 9.9197e-01 -2.0826e-03 1.5991e+00)
(3.3278e+00 -8.0180e-01 5.7876e-02 9.9396e-01 -2.1101e-03 1.5872e+00)
(3.2557e+00 -7.9592e-01 5.7034e-02 9.9455e-01 -2.1628e-03 1.5754e+00)
(3.1886e+00 -7.8961e-01 5.5746e-02 9.9379e-01 -2.2441e-03 1.5637e+00)
(3.1269e+00 -7.8296e-01 5.4259e-02 9.9183e-01 -2.3220e-03 1.5510e+00)
(3.0702e+00 -7.7613e-01 5.2957e-02 9.8887e-01 -2.3519e-03 1.5367e+00)