Created validationSuite
First validation test case: backwardFacingStep with simpleFoam, comparison of velocity profiles and pressure coefficient with experimental data. Author: Robert Keser, Merge: Vuko Vukcevic
This commit is contained in:
Vuko Vukcevic
parent
5ba5d4ab07
commit
b6cf22304f
31 changed files with 14823 additions and 0 deletions
|
|
@ -0,0 +1,52 @@
|
|||
/*--------------------------------*- C++ -*----------------------------------*\
|
||||
| ========= | |
|
||||
| \\ / F ield | foam-extend: Open Source CFD |
|
||||
| \\ / O peration | Version: 3.2 |
|
||||
| \\ / A nd | Web: http://www.foam-extend.org |
|
||||
| \\/ M anipulation | |
|
||||
\*---------------------------------------------------------------------------*/
|
||||
FoamFile
|
||||
{
|
||||
version 2.0;
|
||||
format ascii;
|
||||
class volVectorField;
|
||||
object U;
|
||||
}
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
dimensions [0 1 -1 0 0 0 0];
|
||||
|
||||
internalField uniform (44.31525 0 0);
|
||||
|
||||
boundaryField
|
||||
{
|
||||
inlet
|
||||
{
|
||||
type fixedValue;
|
||||
value $internalField;
|
||||
}
|
||||
outlet
|
||||
{
|
||||
type zeroGradient;
|
||||
}
|
||||
upperWall
|
||||
{
|
||||
type fixedValue;
|
||||
value uniform (0 0 0);
|
||||
}
|
||||
lowerWall
|
||||
{
|
||||
type fixedValue;
|
||||
value uniform (0 0 0);
|
||||
}
|
||||
frontAndBack
|
||||
{
|
||||
type empty;
|
||||
}
|
||||
symmetry
|
||||
{
|
||||
type symmetryPlane;
|
||||
}
|
||||
}
|
||||
|
||||
// ************************************************************************* //
|
||||
|
|
@ -0,0 +1,61 @@
|
|||
/*--------------------------------*- C++ -*----------------------------------*\
|
||||
| ========= | |
|
||||
| \\ / F ield | foam-extend: Open Source CFD |
|
||||
| \\ / O peration | Version: 3.2 |
|
||||
| \\ / A nd | Web: http://www.foam-extend.org |
|
||||
| \\/ M anipulation | |
|
||||
\*---------------------------------------------------------------------------*/
|
||||
FoamFile
|
||||
{
|
||||
version 2.0;
|
||||
format ascii;
|
||||
class volScalarField;
|
||||
location "0";
|
||||
object epsilon;
|
||||
}
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
dimensions [0 2 -3 0 0 0 0];
|
||||
|
||||
internalField uniform 0.079598;
|
||||
|
||||
boundaryField
|
||||
{
|
||||
inlet
|
||||
{
|
||||
type fixedValue;
|
||||
value $internalField;
|
||||
}
|
||||
outlet
|
||||
{
|
||||
type zeroGradient;
|
||||
}
|
||||
upperWall
|
||||
{
|
||||
type epsilonWallFunction;
|
||||
refValue uniform 0;
|
||||
value $internalField;
|
||||
Cmu 0.09;
|
||||
kappa 0.41;
|
||||
E 9.8;
|
||||
}
|
||||
lowerWall
|
||||
{
|
||||
type epsilonWallFunction;
|
||||
refValue uniform 0;
|
||||
value $internalField;
|
||||
Cmu 0.09;
|
||||
kappa 0.41;
|
||||
E 9.8;
|
||||
}
|
||||
frontAndBack
|
||||
{
|
||||
type empty;
|
||||
}
|
||||
symmetry
|
||||
{
|
||||
type symmetryPlane;
|
||||
}
|
||||
}
|
||||
|
||||
// ************************************************************************* //
|
||||
|
|
@ -0,0 +1,53 @@
|
|||
/*--------------------------------*- C++ -*----------------------------------*\
|
||||
| ========= | |
|
||||
| \\ / F ield | foam-extend: Open Source CFD |
|
||||
| \\ / O peration | Version: 3.2 |
|
||||
| \\ / A nd | Web: http://www.foam-extend.org |
|
||||
| \\/ M anipulation | |
|
||||
\*---------------------------------------------------------------------------*/
|
||||
FoamFile
|
||||
{
|
||||
version 2.0;
|
||||
format ascii;
|
||||
class volScalarField;
|
||||
location "0";
|
||||
object k;
|
||||
}
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
dimensions [0 2 -2 0 0 0 0];
|
||||
|
||||
internalField uniform 294.57e-03;
|
||||
|
||||
boundaryField
|
||||
{
|
||||
inlet
|
||||
{
|
||||
type fixedValue;
|
||||
value $internalField;
|
||||
}
|
||||
outlet
|
||||
{
|
||||
type zeroGradient;
|
||||
}
|
||||
upperWall
|
||||
{
|
||||
type kqRWallFunction;
|
||||
value $internalField;
|
||||
}
|
||||
lowerWall
|
||||
{
|
||||
type kqRWallFunction;
|
||||
value $internalField;
|
||||
}
|
||||
frontAndBack
|
||||
{
|
||||
type empty;
|
||||
}
|
||||
symmetry
|
||||
{
|
||||
type symmetryPlane;
|
||||
}
|
||||
}
|
||||
|
||||
// ************************************************************************* //
|
||||
|
|
@ -0,0 +1,53 @@
|
|||
/*--------------------------------*- C++ -*----------------------------------*\
|
||||
| ========= | |
|
||||
| \\ / F ield | foam-extend: Open Source CFD |
|
||||
| \\ / O peration | Version: 3.2 |
|
||||
| \\ / A nd | Web: http://www.foam-extend.org |
|
||||
| \\/ M anipulation | |
|
||||
\*---------------------------------------------------------------------------*/
|
||||
FoamFile
|
||||
{
|
||||
version 2.0;
|
||||
format ascii;
|
||||
class volScalarField;
|
||||
location "0";
|
||||
object omega;
|
||||
}
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
dimensions [0 0 -1 0 0 0 0];
|
||||
|
||||
internalField uniform 97.37245;
|
||||
|
||||
boundaryField
|
||||
{
|
||||
inlet
|
||||
{
|
||||
type fixedValue;
|
||||
value $internalField;
|
||||
}
|
||||
outlet
|
||||
{
|
||||
type zeroGradient;
|
||||
}
|
||||
upperWall
|
||||
{
|
||||
type omegaWallFunction;
|
||||
value $internalField;
|
||||
}
|
||||
lowerWall
|
||||
{
|
||||
type omegaWallFunction;
|
||||
value $internalField;
|
||||
}
|
||||
frontAndBack
|
||||
{
|
||||
type empty;
|
||||
}
|
||||
symmetry
|
||||
{
|
||||
type symmetryPlane;
|
||||
}
|
||||
}
|
||||
|
||||
// ************************************************************************* //
|
||||
|
|
@ -0,0 +1,54 @@
|
|||
/*--------------------------------*- C++ -*----------------------------------*\
|
||||
| ========= | |
|
||||
| \\ / F ield | foam-extend: Open Source CFD |
|
||||
| \\ / O peration | Version: 3.2 |
|
||||
| \\ / A nd | Web: http://www.foam-extend.org |
|
||||
| \\/ M anipulation | |
|
||||
\*---------------------------------------------------------------------------*/
|
||||
FoamFile
|
||||
{
|
||||
version 2.0;
|
||||
format ascii;
|
||||
class volScalarField;
|
||||
object p;
|
||||
}
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
dimensions [0 2 -2 0 0 0 0];
|
||||
|
||||
internalField uniform 0;
|
||||
|
||||
boundaryField
|
||||
{
|
||||
inlet
|
||||
{
|
||||
type zeroGradient;
|
||||
}
|
||||
|
||||
outlet
|
||||
{
|
||||
type fixedValue;
|
||||
value uniform 0;
|
||||
}
|
||||
|
||||
upperWall
|
||||
{
|
||||
type zeroGradient;
|
||||
}
|
||||
|
||||
lowerWall
|
||||
{
|
||||
type zeroGradient;
|
||||
}
|
||||
|
||||
frontAndBack
|
||||
{
|
||||
type empty;
|
||||
}
|
||||
symmetry
|
||||
{
|
||||
type symmetryPlane;
|
||||
}
|
||||
}
|
||||
|
||||
// ************************************************************************* //
|
||||
12
validationSuite/incompressible/simpleFoam/backwardFacingStep/Allclean
Executable file
12
validationSuite/incompressible/simpleFoam/backwardFacingStep/Allclean
Executable file
|
|
@ -0,0 +1,12 @@
|
|||
#!/bin/bash
|
||||
|
||||
# Source tutorial run functions
|
||||
. $WM_PROJECT_DIR/bin/tools/CleanFunctions
|
||||
|
||||
cleanCase
|
||||
\rm -rf 0
|
||||
\rm -rf surfaces sets
|
||||
\rm -f validationData/*.raw
|
||||
\rm -f validationData/*.xy
|
||||
\rm -rf validationResults
|
||||
cp -r save 0
|
||||
25
validationSuite/incompressible/simpleFoam/backwardFacingStep/Allrun
Executable file
25
validationSuite/incompressible/simpleFoam/backwardFacingStep/Allrun
Executable file
|
|
@ -0,0 +1,25 @@
|
|||
#!/bin/sh
|
||||
# Source tutorial run functions
|
||||
. $WM_PROJECT_DIR/bin/tools/RunFunctions
|
||||
|
||||
application=simpleFoam
|
||||
|
||||
runApplication blockMesh
|
||||
runApplication potentialFoam
|
||||
runApplication $application
|
||||
runApplication sample -latestTime
|
||||
|
||||
timeStep=$(grep "^endTime" system/controlDict | awk '{ print $2 }' | sed 's/;//g')
|
||||
|
||||
cp -v surfaces/$timeStep/p_lowerWall.raw validationData/
|
||||
cp -v sets/$timeStep/profile_?_U.xy validationData/
|
||||
|
||||
./validationData/gnuplot
|
||||
|
||||
mkdir -v validationResults
|
||||
|
||||
mv -v *.png validationResults/
|
||||
|
||||
cp -v validationData/BFSLines.pdf validationResults/
|
||||
|
||||
gnome-open validationResults/Cp.png
|
||||
|
|
@ -0,0 +1,21 @@
|
|||
/*--------------------------------*- C++ -*----------------------------------*\
|
||||
| ========= | |
|
||||
| \\ / F ield | foam-extend: Open Source CFD |
|
||||
| \\ / O peration | Version: 3.2 |
|
||||
| \\ / A nd | Web: http://www.foam-extend.org |
|
||||
| \\/ M anipulation | |
|
||||
\*---------------------------------------------------------------------------*/
|
||||
FoamFile
|
||||
{
|
||||
version 2.0;
|
||||
format ascii;
|
||||
class dictionary;
|
||||
object RASProperties;
|
||||
}
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
RASModel kOmegaSST;//kEpsilon;
|
||||
|
||||
turbulence on;
|
||||
|
||||
// ************************************************************************* //
|
||||
|
|
@ -0,0 +1,104 @@
|
|||
/*--------------------------------*- C++ -*----------------------------------*\
|
||||
| ========= | |
|
||||
| \\ / F ield | foam-extend: Open Source CFD |
|
||||
| \\ / O peration | Version: 3.2 |
|
||||
| \\ / A nd | Web: http://www.foam-extend.org |
|
||||
| \\/ M anipulation | |
|
||||
\*---------------------------------------------------------------------------*/
|
||||
FoamFile
|
||||
{
|
||||
version 2.0;
|
||||
format ascii;
|
||||
class dictionary;
|
||||
object blockMeshDict;
|
||||
}
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
convertToMeters 1;//0.0127;
|
||||
|
||||
vertices
|
||||
(
|
||||
(-110 9 0)
|
||||
(-110 1 0)
|
||||
(0 1 0)
|
||||
(0 0 0)
|
||||
(50 0 0)
|
||||
(50 1 0)
|
||||
(50 9 0)
|
||||
(0 9 0)
|
||||
|
||||
(-110 9 1)
|
||||
(-110 1 1)
|
||||
(0 1 1)
|
||||
(0 0 1)
|
||||
(50 0 1)
|
||||
(50 1 1)
|
||||
(50 9 1)
|
||||
(0 9 1)
|
||||
|
||||
(-130 9 0)
|
||||
(-130 1 0)
|
||||
|
||||
(-130 9 1)
|
||||
(-130 1 1)
|
||||
);
|
||||
|
||||
blocks
|
||||
(
|
||||
hex (17 1 0 16 19 9 8 18) (15 45 1) simpleGrading (0.85925 7.95 1) //(100 40 1)
|
||||
hex (1 2 7 0 9 10 15 8) (262 45 1) simpleGrading (0.04368 7.95 1) //(550 40 1)
|
||||
hex (2 5 6 7 10 13 14 15) (330 45 1) simpleGrading (5.1341 7.95 1) //(250 40 1)
|
||||
hex (3 4 5 2 11 12 13 10) (330 19 1) simpleGrading (5.1341 1 1) //(250 10 1)
|
||||
);
|
||||
|
||||
edges
|
||||
(
|
||||
);
|
||||
|
||||
patches
|
||||
(
|
||||
patch inlet
|
||||
(
|
||||
(16 17 19 18)
|
||||
)
|
||||
patch outlet
|
||||
(
|
||||
(5 6 14 13)
|
||||
(4 5 13 12)
|
||||
|
||||
)
|
||||
patch symmetry
|
||||
(
|
||||
(0 16 18 8)
|
||||
(1 9 19 17)
|
||||
|
||||
)
|
||||
wall upperWall
|
||||
(
|
||||
(0 8 15 7)
|
||||
(7 15 14 6)
|
||||
)
|
||||
wall lowerWall
|
||||
(
|
||||
(1 9 10 2)
|
||||
(2 10 11 3)
|
||||
(3 11 12 4)
|
||||
)
|
||||
empty frontAndBack
|
||||
(
|
||||
(0 7 2 1)
|
||||
(2 7 6 5)
|
||||
(2 5 4 3)
|
||||
(8 9 10 15)
|
||||
(10 13 14 15)
|
||||
(10 11 12 13)
|
||||
(8 18 19 9)
|
||||
(1 17 16 0)
|
||||
)
|
||||
);
|
||||
|
||||
mergePatchPairs
|
||||
(
|
||||
);
|
||||
|
||||
// ************************************************************************* //
|
||||
|
|
@ -0,0 +1,59 @@
|
|||
/*--------------------------------*- C++ -*----------------------------------*\
|
||||
| ========= | |
|
||||
| \\ / F ield | foam-extend: Open Source CFD |
|
||||
| \\ / O peration | Version: 4.0 |
|
||||
| \\ / A nd | Web: http://www.foam-extend.org |
|
||||
| \\/ M anipulation | |
|
||||
\*---------------------------------------------------------------------------*/
|
||||
FoamFile
|
||||
{
|
||||
version 2.0;
|
||||
format ascii;
|
||||
class polyBoundaryMesh;
|
||||
location "constant/polyMesh";
|
||||
object boundary;
|
||||
}
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
6
|
||||
(
|
||||
inlet
|
||||
{
|
||||
type patch;
|
||||
nFaces 45;
|
||||
startFace 66499;
|
||||
}
|
||||
outlet
|
||||
{
|
||||
type patch;
|
||||
nFaces 64;
|
||||
startFace 66544;
|
||||
}
|
||||
symmetry
|
||||
{
|
||||
type patch;
|
||||
nFaces 30;
|
||||
startFace 66608;
|
||||
}
|
||||
upperWall
|
||||
{
|
||||
type wall;
|
||||
nFaces 592;
|
||||
startFace 66638;
|
||||
}
|
||||
lowerWall
|
||||
{
|
||||
type wall;
|
||||
nFaces 611;
|
||||
startFace 67230;
|
||||
}
|
||||
frontAndBack
|
||||
{
|
||||
type empty;
|
||||
nFaces 67170;
|
||||
startFace 67841;
|
||||
}
|
||||
)
|
||||
|
||||
|
||||
// ************************************************************************* //
|
||||
|
|
@ -0,0 +1,21 @@
|
|||
/*--------------------------------*- C++ -*----------------------------------*\
|
||||
| ========= | |
|
||||
| \\ / F ield | foam-extend: Open Source CFD |
|
||||
| \\ / O peration | Version: 3.2 |
|
||||
| \\ / A nd | Web: http://www.foam-extend.org |
|
||||
| \\/ M anipulation | |
|
||||
\*---------------------------------------------------------------------------*/
|
||||
FoamFile
|
||||
{
|
||||
version 2.0;
|
||||
format ascii;
|
||||
class dictionary;
|
||||
object transportProperties;
|
||||
}
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
transportModel Newtonian;
|
||||
|
||||
nu nu [0 2 -1 0 0 0 0] 1.230979e-03;
|
||||
|
||||
// ************************************************************************* //
|
||||
|
|
@ -0,0 +1,52 @@
|
|||
/*--------------------------------*- C++ -*----------------------------------*\
|
||||
| ========= | |
|
||||
| \\ / F ield | foam-extend: Open Source CFD |
|
||||
| \\ / O peration | Version: 3.2 |
|
||||
| \\ / A nd | Web: http://www.foam-extend.org |
|
||||
| \\/ M anipulation | |
|
||||
\*---------------------------------------------------------------------------*/
|
||||
FoamFile
|
||||
{
|
||||
version 2.0;
|
||||
format ascii;
|
||||
class volVectorField;
|
||||
object U;
|
||||
}
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
dimensions [0 1 -1 0 0 0 0];
|
||||
|
||||
internalField uniform (44.31525 0 0);
|
||||
|
||||
boundaryField
|
||||
{
|
||||
inlet
|
||||
{
|
||||
type fixedValue;
|
||||
value $internalField;
|
||||
}
|
||||
outlet
|
||||
{
|
||||
type zeroGradient;
|
||||
}
|
||||
upperWall
|
||||
{
|
||||
type fixedValue;
|
||||
value uniform (0 0 0);
|
||||
}
|
||||
lowerWall
|
||||
{
|
||||
type fixedValue;
|
||||
value uniform (0 0 0);
|
||||
}
|
||||
frontAndBack
|
||||
{
|
||||
type empty;
|
||||
}
|
||||
symmetry
|
||||
{
|
||||
type symmetryPlane;
|
||||
}
|
||||
}
|
||||
|
||||
// ************************************************************************* //
|
||||
|
|
@ -0,0 +1,61 @@
|
|||
/*--------------------------------*- C++ -*----------------------------------*\
|
||||
| ========= | |
|
||||
| \\ / F ield | foam-extend: Open Source CFD |
|
||||
| \\ / O peration | Version: 3.2 |
|
||||
| \\ / A nd | Web: http://www.foam-extend.org |
|
||||
| \\/ M anipulation | |
|
||||
\*---------------------------------------------------------------------------*/
|
||||
FoamFile
|
||||
{
|
||||
version 2.0;
|
||||
format ascii;
|
||||
class volScalarField;
|
||||
location "0";
|
||||
object epsilon;
|
||||
}
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
dimensions [0 2 -3 0 0 0 0];
|
||||
|
||||
internalField uniform 0.079598;
|
||||
|
||||
boundaryField
|
||||
{
|
||||
inlet
|
||||
{
|
||||
type fixedValue;
|
||||
value $internalField;
|
||||
}
|
||||
outlet
|
||||
{
|
||||
type zeroGradient;
|
||||
}
|
||||
upperWall
|
||||
{
|
||||
type epsilonWallFunction;
|
||||
refValue uniform 0;
|
||||
value $internalField;
|
||||
Cmu 0.09;
|
||||
kappa 0.41;
|
||||
E 9.8;
|
||||
}
|
||||
lowerWall
|
||||
{
|
||||
type epsilonWallFunction;
|
||||
refValue uniform 0;
|
||||
value $internalField;
|
||||
Cmu 0.09;
|
||||
kappa 0.41;
|
||||
E 9.8;
|
||||
}
|
||||
frontAndBack
|
||||
{
|
||||
type empty;
|
||||
}
|
||||
symmetry
|
||||
{
|
||||
type symmetryPlane;
|
||||
}
|
||||
}
|
||||
|
||||
// ************************************************************************* //
|
||||
|
|
@ -0,0 +1,53 @@
|
|||
/*--------------------------------*- C++ -*----------------------------------*\
|
||||
| ========= | |
|
||||
| \\ / F ield | foam-extend: Open Source CFD |
|
||||
| \\ / O peration | Version: 3.2 |
|
||||
| \\ / A nd | Web: http://www.foam-extend.org |
|
||||
| \\/ M anipulation | |
|
||||
\*---------------------------------------------------------------------------*/
|
||||
FoamFile
|
||||
{
|
||||
version 2.0;
|
||||
format ascii;
|
||||
class volScalarField;
|
||||
location "0";
|
||||
object k;
|
||||
}
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
dimensions [0 2 -2 0 0 0 0];
|
||||
|
||||
internalField uniform 294.57e-03;
|
||||
|
||||
boundaryField
|
||||
{
|
||||
inlet
|
||||
{
|
||||
type fixedValue;
|
||||
value $internalField;
|
||||
}
|
||||
outlet
|
||||
{
|
||||
type zeroGradient;
|
||||
}
|
||||
upperWall
|
||||
{
|
||||
type kqRWallFunction;
|
||||
value $internalField;
|
||||
}
|
||||
lowerWall
|
||||
{
|
||||
type kqRWallFunction;
|
||||
value $internalField;
|
||||
}
|
||||
frontAndBack
|
||||
{
|
||||
type empty;
|
||||
}
|
||||
symmetry
|
||||
{
|
||||
type symmetryPlane;
|
||||
}
|
||||
}
|
||||
|
||||
// ************************************************************************* //
|
||||
|
|
@ -0,0 +1,53 @@
|
|||
/*--------------------------------*- C++ -*----------------------------------*\
|
||||
| ========= | |
|
||||
| \\ / F ield | foam-extend: Open Source CFD |
|
||||
| \\ / O peration | Version: 3.2 |
|
||||
| \\ / A nd | Web: http://www.foam-extend.org |
|
||||
| \\/ M anipulation | |
|
||||
\*---------------------------------------------------------------------------*/
|
||||
FoamFile
|
||||
{
|
||||
version 2.0;
|
||||
format ascii;
|
||||
class volScalarField;
|
||||
location "0";
|
||||
object omega;
|
||||
}
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
dimensions [0 0 -1 0 0 0 0];
|
||||
|
||||
internalField uniform 97.37245;
|
||||
|
||||
boundaryField
|
||||
{
|
||||
inlet
|
||||
{
|
||||
type fixedValue;
|
||||
value $internalField;
|
||||
}
|
||||
outlet
|
||||
{
|
||||
type zeroGradient;
|
||||
}
|
||||
upperWall
|
||||
{
|
||||
type omegaWallFunction;
|
||||
value $internalField;
|
||||
}
|
||||
lowerWall
|
||||
{
|
||||
type omegaWallFunction;
|
||||
value $internalField;
|
||||
}
|
||||
frontAndBack
|
||||
{
|
||||
type empty;
|
||||
}
|
||||
symmetry
|
||||
{
|
||||
type symmetryPlane;
|
||||
}
|
||||
}
|
||||
|
||||
// ************************************************************************* //
|
||||
|
|
@ -0,0 +1,54 @@
|
|||
/*--------------------------------*- C++ -*----------------------------------*\
|
||||
| ========= | |
|
||||
| \\ / F ield | foam-extend: Open Source CFD |
|
||||
| \\ / O peration | Version: 3.2 |
|
||||
| \\ / A nd | Web: http://www.foam-extend.org |
|
||||
| \\/ M anipulation | |
|
||||
\*---------------------------------------------------------------------------*/
|
||||
FoamFile
|
||||
{
|
||||
version 2.0;
|
||||
format ascii;
|
||||
class volScalarField;
|
||||
object p;
|
||||
}
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
dimensions [0 2 -2 0 0 0 0];
|
||||
|
||||
internalField uniform 0;
|
||||
|
||||
boundaryField
|
||||
{
|
||||
inlet
|
||||
{
|
||||
type zeroGradient;
|
||||
}
|
||||
|
||||
outlet
|
||||
{
|
||||
type fixedValue;
|
||||
value uniform 0;
|
||||
}
|
||||
|
||||
upperWall
|
||||
{
|
||||
type zeroGradient;
|
||||
}
|
||||
|
||||
lowerWall
|
||||
{
|
||||
type zeroGradient;
|
||||
}
|
||||
|
||||
frontAndBack
|
||||
{
|
||||
type empty;
|
||||
}
|
||||
symmetry
|
||||
{
|
||||
type symmetryPlane;
|
||||
}
|
||||
}
|
||||
|
||||
// ************************************************************************* //
|
||||
|
|
@ -0,0 +1,47 @@
|
|||
/*--------------------------------*- C++ -*----------------------------------*\
|
||||
| ========= | |
|
||||
| \\ / F ield | foam-extend: Open Source CFD |
|
||||
| \\ / O peration | Version: 3.2 |
|
||||
| \\ / A nd | Web: http://www.foam-extend.org |
|
||||
| \\/ M anipulation | |
|
||||
\*---------------------------------------------------------------------------*/
|
||||
FoamFile
|
||||
{
|
||||
version 2.0;
|
||||
format ascii;
|
||||
class dictionary;
|
||||
object controlDict;
|
||||
}
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
application simpleFoam;
|
||||
|
||||
startFrom startTime;
|
||||
|
||||
startTime 0;
|
||||
|
||||
stopAt endTime;
|
||||
|
||||
endTime 1000;
|
||||
|
||||
deltaT 1;
|
||||
|
||||
writeControl timeStep;
|
||||
|
||||
writeInterval 1000;
|
||||
|
||||
purgeWrite 0;
|
||||
|
||||
writeFormat ascii;
|
||||
|
||||
writePrecision 6;
|
||||
|
||||
writeCompression uncompressed;
|
||||
|
||||
timeFormat general;
|
||||
|
||||
timePrecision 6;
|
||||
|
||||
runTimeModifiable yes;
|
||||
|
||||
// ************************************************************************* //
|
||||
|
|
@ -0,0 +1,48 @@
|
|||
/*--------------------------------*- C++ -*----------------------------------*\
|
||||
| ========= | |
|
||||
| \\ / F ield | foam-extend: Open Source CFD |
|
||||
| \\ / O peration | Version: 3.2 |
|
||||
| \\ / A nd | Web: http://www.foam-extend.org |
|
||||
| \\/ M anipulation | |
|
||||
\*---------------------------------------------------------------------------*/
|
||||
FoamFile
|
||||
{
|
||||
version 2.0;
|
||||
format ascii;
|
||||
class dictionary;
|
||||
location "system";
|
||||
object decomposeParDict;
|
||||
}
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
numberOfSubdomains 2;
|
||||
|
||||
method metis;
|
||||
|
||||
simpleCoeffs
|
||||
{
|
||||
n ( 2 1 1 );
|
||||
delta 0.001;
|
||||
}
|
||||
|
||||
hierarchicalCoeffs
|
||||
{
|
||||
n ( 2 1 1 );
|
||||
delta 0.001;
|
||||
order xyz;
|
||||
}
|
||||
|
||||
metisCoeffs
|
||||
{
|
||||
}
|
||||
|
||||
manualCoeffs
|
||||
{
|
||||
dataFile "";
|
||||
}
|
||||
|
||||
distributed no;
|
||||
|
||||
roots ( );
|
||||
|
||||
// ************************************************************************* //
|
||||
|
|
@ -0,0 +1,55 @@
|
|||
/*--------------------------------*- C++ -*----------------------------------*\
|
||||
| ========= | |
|
||||
| \\ / F ield | foam-extend: Open Source CFD |
|
||||
| \\ / O peration | Version: 3.2 |
|
||||
| \\ / A nd | Web: http://www.foam-extend.org |
|
||||
| \\/ M anipulation | |
|
||||
\*---------------------------------------------------------------------------*/
|
||||
FoamFile
|
||||
{
|
||||
version 2.0;
|
||||
format ascii;
|
||||
class dictionary;
|
||||
object fvSchemes;
|
||||
}
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
ddtSchemes
|
||||
{
|
||||
default steadyState;
|
||||
}
|
||||
|
||||
gradSchemes
|
||||
{
|
||||
default Gauss linear;
|
||||
}
|
||||
|
||||
divSchemes
|
||||
{
|
||||
default none;
|
||||
// div(phi,U) Gauss linearUpwind Gauss linear;
|
||||
div(phi,U) Gauss upwind;
|
||||
div(phi,k) Gauss upwind;
|
||||
div(phi,epsilon) Gauss upwind;
|
||||
div(phi,omega) Gauss upwind;
|
||||
div((nuEff*dev(grad(U).T()))) Gauss linear;
|
||||
div((nuEff*dev(T(grad(U))))) Gauss linear;
|
||||
}
|
||||
|
||||
laplacianSchemes
|
||||
{
|
||||
default Gauss linear corrected;
|
||||
}
|
||||
|
||||
interpolationSchemes
|
||||
{
|
||||
default linear;
|
||||
}
|
||||
|
||||
snGradSchemes
|
||||
{
|
||||
default corrected;
|
||||
}
|
||||
|
||||
|
||||
// ************************************************************************* //
|
||||
|
|
@ -0,0 +1,70 @@
|
|||
/*--------------------------------*- C++ -*----------------------------------*\
|
||||
| ========= | |
|
||||
| \\ / F ield | foam-extend: Open Source CFD |
|
||||
| \\ / O peration | Version: 3.2 |
|
||||
| \\ / A nd | Web: http://www.foam-extend.org |
|
||||
| \\/ M anipulation | |
|
||||
\*---------------------------------------------------------------------------*/
|
||||
FoamFile
|
||||
{
|
||||
version 2.0;
|
||||
format ascii;
|
||||
class dictionary;
|
||||
object fvSolution;
|
||||
}
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
solvers
|
||||
{
|
||||
p
|
||||
{
|
||||
solver GAMG;
|
||||
tolerance 1e-7;
|
||||
relTol 0.1;
|
||||
smoother GaussSeidel;
|
||||
nPreSweeps 0;
|
||||
nPostSweeps 2;
|
||||
cacheAgglomeration on;
|
||||
agglomerator faceAreaPair;
|
||||
nCellsInCoarsestLevel 10;
|
||||
mergeLevels 1;
|
||||
}
|
||||
"U|k|epsilon|omega"
|
||||
{
|
||||
solver smoothSolver;
|
||||
smoother GaussSeidel;
|
||||
tolerance 1e-8;
|
||||
relTol 0.1;
|
||||
nSweeps 1;
|
||||
|
||||
}
|
||||
}
|
||||
|
||||
SIMPLE
|
||||
{
|
||||
nNonOrthoCorrectors 0;
|
||||
}
|
||||
|
||||
relaxationFactors
|
||||
{
|
||||
fields
|
||||
{
|
||||
p 0.3;
|
||||
}
|
||||
equations
|
||||
{
|
||||
U 0.7;
|
||||
k 0.7;
|
||||
epsilon 0.7;
|
||||
omega 0.7;
|
||||
}
|
||||
}
|
||||
|
||||
fieldBounds
|
||||
{
|
||||
U 500;
|
||||
p -5e4 5e4;
|
||||
}
|
||||
|
||||
|
||||
// ************************************************************************* //
|
||||
|
|
@ -0,0 +1,249 @@
|
|||
/*--------------------------------*- C++ -*----------------------------------*\
|
||||
| ========= | |
|
||||
| \\ / F ield | foam-extend: Open Source CFD |
|
||||
| \\ / O peration | Version: 4.0 |
|
||||
| \\ / A nd | Web: http://www.foam-extend.org |
|
||||
| \\/ M anipulation | |
|
||||
\*---------------------------------------------------------------------------*/
|
||||
FoamFile
|
||||
{
|
||||
version 2.0;
|
||||
format ascii;
|
||||
class dictionary;
|
||||
location system;
|
||||
object sampleDict;
|
||||
}
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
// Set output format : choice of
|
||||
// xmgr
|
||||
// jplot
|
||||
// gnuplot
|
||||
// raw
|
||||
setFormat raw;
|
||||
|
||||
// Surface output format. Choice of
|
||||
// null : suppress output
|
||||
// foamFile : separate points, faces and values file
|
||||
// dx : DX scalar or vector format
|
||||
// vtk : VTK ascii format
|
||||
// raw : x y z value format for use with e.g. gnuplot 'splot'.
|
||||
//
|
||||
// Note:
|
||||
// other formats such as obj, stl, etc can also be written (by proxy)
|
||||
// but without any values!
|
||||
surfaceFormat raw;
|
||||
|
||||
// interpolationScheme. choice of
|
||||
// cell : use cell-centre value only; constant over cells (default)
|
||||
// cellPoint : use cell-centre and vertex values
|
||||
// cellPointFace : use cell-centre, vertex and face values.
|
||||
// 1] vertex values determined from neighbouring cell-centre values
|
||||
// 2] face values determined using the current face interpolation scheme
|
||||
// for the field (linear, gamma, etc.)
|
||||
interpolationScheme cellPoint;
|
||||
|
||||
// Fields to sample.
|
||||
fields
|
||||
(
|
||||
p
|
||||
U
|
||||
);
|
||||
|
||||
// Set sampling definition: choice of
|
||||
// uniform evenly distributed points on line
|
||||
// face one point per face intersection
|
||||
// midPoint one point per cell, inbetween two face intersections
|
||||
// midPointAndFace combination of face and midPoint
|
||||
//
|
||||
// curve specified points, not nessecary on line, uses
|
||||
// tracking
|
||||
// cloud specified points, uses findCell
|
||||
//
|
||||
// axis: how to write point coordinate. Choice of
|
||||
// - x/y/z: x/y/z coordinate only
|
||||
// - xyz: three columns
|
||||
// (probably does not make sense for anything but raw)
|
||||
// - distance: distance from start of sampling line (if uses line) or
|
||||
// distance from first specified sampling point
|
||||
//
|
||||
// type specific:
|
||||
// uniform, face, midPoint, midPointAndFace : start and end coordinate
|
||||
// uniform: extra number of sampling points
|
||||
// curve, cloud: list of coordinates
|
||||
sets
|
||||
(
|
||||
profile_1
|
||||
{
|
||||
type uniform;
|
||||
axis distance;
|
||||
|
||||
//- cavity. Slightly perturbed so not to align with face or edge.
|
||||
start (-4 1.0 0.5);
|
||||
end (-4 8.2 0.5);
|
||||
nPoints 200;
|
||||
}
|
||||
profile_2
|
||||
{
|
||||
type uniform;
|
||||
axis distance;
|
||||
|
||||
//- cavity. Slightly perturbed so not to align with face or edge.
|
||||
start (1 0.0 0.5);
|
||||
end (1 8.2 0.5);
|
||||
nPoints 200;
|
||||
}
|
||||
profile_3
|
||||
{
|
||||
type uniform;
|
||||
axis distance;
|
||||
|
||||
//- cavity. Slightly perturbed so not to align with face or edge.
|
||||
start (4 0.0 0.5);
|
||||
end (4 8.2 0.5);
|
||||
nPoints 200;
|
||||
}
|
||||
profile_4
|
||||
{
|
||||
type uniform;
|
||||
axis distance;
|
||||
|
||||
//- cavity. Slightly perturbed so not to align with face or edge.
|
||||
start (6 0.0 0.5);
|
||||
end (6 8.2 0.5);
|
||||
nPoints 200;
|
||||
}
|
||||
profile_5
|
||||
{
|
||||
type uniform;
|
||||
axis distance;
|
||||
|
||||
//- cavity. Slightly perturbed so not to align with face or edge.
|
||||
start (10 0.0 0.5);
|
||||
end (10 8.2 0.5);
|
||||
nPoints 200;
|
||||
}
|
||||
//
|
||||
// lineX2
|
||||
// {
|
||||
// type face;
|
||||
// axis x;
|
||||
//
|
||||
// //- flangeHex
|
||||
// //start (0 20 -20);
|
||||
// //end (0 20 10);
|
||||
//
|
||||
// //- nablaCavity
|
||||
// //start (-1 0.05 0.005);
|
||||
// //end ( 1 0.05 0.005);
|
||||
//
|
||||
// //- cavity
|
||||
// start (0.001 0.5101 0.00501);
|
||||
// end (2.01 0.5101 0.00501);
|
||||
// nPoints 10;
|
||||
// }
|
||||
//
|
||||
// somePoints
|
||||
// {
|
||||
// type cloud;
|
||||
// axis xyz;
|
||||
// points ((0.049 0.049 0.00501)(0.051 0.049 0.00501));
|
||||
// }
|
||||
|
||||
);
|
||||
|
||||
// Surface sampling definition: choice of
|
||||
// plane : values on plane defined by point, normal.
|
||||
// patch : values on patch.
|
||||
//
|
||||
// 1] patches are not triangulated by default
|
||||
// 2] planes are always triangulated
|
||||
// 3] iso-surfaces are always triangulated
|
||||
surfaces
|
||||
(
|
||||
// constantPlane
|
||||
// {
|
||||
// type plane; // always triangulated
|
||||
// basePoint (0.0501 0.0501 0.005);
|
||||
// normalVector (0.1 0.1 1);
|
||||
//
|
||||
// //- Optional: restrict to a particular zone
|
||||
// // zoneName zone1;
|
||||
// }
|
||||
//
|
||||
// interpolatedPlane
|
||||
// {
|
||||
// type plane; // always triangulated
|
||||
// // make plane relative to the coordinateSystem (Cartesian)
|
||||
// coordinateSystem
|
||||
// {
|
||||
// origin (0.0501 0.0501 0.005);
|
||||
// }
|
||||
// basePoint (0 0 0);
|
||||
// normalVector (0.1 0.1 1);
|
||||
// interpolate true;
|
||||
// }
|
||||
//
|
||||
// movingWall_constant
|
||||
// {
|
||||
// type patch;
|
||||
// patchName movingWall;
|
||||
// // Optional: whether to leave as faces (=default) or triangulate
|
||||
// // triangulate false;
|
||||
// }
|
||||
|
||||
lowerWall
|
||||
{
|
||||
type patch;
|
||||
patchName lowerWall;
|
||||
interpolate false;
|
||||
// Optional: whether to leave as faces (=default) or triangulate
|
||||
// triangulate false;
|
||||
}
|
||||
|
||||
// interpolatedIso
|
||||
// {
|
||||
// // Iso surface for interpolated values only
|
||||
// type isoSurface; // always triangulated
|
||||
// isoField rho;
|
||||
// isoValue 0.5;
|
||||
// interpolate true;
|
||||
//
|
||||
// //zone ABC; // Optional: zone only
|
||||
// //exposedPatchName fixedWalls; // Optional: zone only
|
||||
//
|
||||
// // regularise false; // Optional: do not simplify
|
||||
// }
|
||||
// constantIso
|
||||
// {
|
||||
// // Iso surface for constant values.
|
||||
// // Triangles guaranteed not to cross cells.
|
||||
// type isoSurfaceCell; // always triangulated
|
||||
// isoField rho;
|
||||
// isoValue 0.5;
|
||||
// interpolate false;
|
||||
// regularise false; // do not simplify
|
||||
// }
|
||||
//
|
||||
// triangleCut
|
||||
// {
|
||||
// // Cutingplane using iso surface
|
||||
// type cuttingPlane;
|
||||
// planeType pointAndNormal;
|
||||
// pointAndNormalDict
|
||||
// {
|
||||
// basePoint (0.4 0 0.4);
|
||||
// normalVector (1 0.2 0.2);
|
||||
// }
|
||||
// interpolate true;
|
||||
//
|
||||
// //zone ABC; // Optional: zone only
|
||||
// //exposedPatchName fixedWalls; // Optional: zone only
|
||||
//
|
||||
// // regularise false; // Optional: do not simplify
|
||||
// }
|
||||
|
||||
);
|
||||
|
||||
|
||||
// *********************************************************************** //
|
||||
File diff suppressed because it is too large
Load diff
|
|
@ -0,0 +1,23 @@
|
|||
# Exp data from Driver and Seegmiller, bottom wall
|
||||
# x means x/H
|
||||
variables="x","cf","error"
|
||||
-3.956 2.88e-3 .202e-3
|
||||
-1.804 2.85e-3 .171e-3
|
||||
-.804 3.11e-3 .311e-3
|
||||
.484 .09e-3 .014e-3
|
||||
1.804 -.11e-3 .019e-3
|
||||
2.804 -.52e-3 .109e-3
|
||||
3.804 -1.02e-3 .163e-3
|
||||
4.804 -.82e-3 .123e-3
|
||||
5.882 -.22e-3 .062e-3
|
||||
7.090 .45e-3 .167e-3
|
||||
7.978 1.02e-3 .122e-3
|
||||
9.196 1.18e-3 .106e-3
|
||||
11.196 1.54e-3 .123e-3
|
||||
13.196 1.73e-3 .121e-3
|
||||
16.196 1.91e-3 .134e-3
|
||||
20.196 2.02e-3 .141e-3
|
||||
24.136 2.04e-3 .143e-3
|
||||
28.130 2.05e-3 .205e-3
|
||||
32.196 2.12e-3 .148e-3
|
||||
35.994 2.30e-3 .184e-3
|
||||
|
|
@ -0,0 +1,59 @@
|
|||
# Exp data from Driver and Seegmiller
|
||||
# x means x/H
|
||||
# cp_orig is the original Cp data, referenced to conditions near x/H=-4
|
||||
# cp has been shifted for convenience so that Cp=0 near x/H=40 (see Eca et al: AIAA-2009-3647, AIAA-2007-4089, AIAA-2005-4728)
|
||||
#variables="x","cp","cp_orig"
|
||||
#zone,t="bottom wall"
|
||||
-8.5 -0.1102 0.0039
|
||||
-6.5 -0.1141 0.0000
|
||||
-4.5 -0.1189 -0.0048
|
||||
-2.5 -0.1372 -0.0231
|
||||
-0.5 -0.1613 -0.0472
|
||||
0.5 -0.1777 -0.0636
|
||||
1.0 -0.1806 -0.0665
|
||||
2.0 -0.1883 -0.0742
|
||||
3.0 -0.1903 -0.0762
|
||||
3.5 -0.1806 -0.0665
|
||||
4.0 -0.1565 -0.0424
|
||||
4.5 -0.1064 0.0077
|
||||
5.0 -0.0659 0.0482
|
||||
5.5 -0.0359 0.0782
|
||||
6.0 -0.0012 0.1129
|
||||
6.5 0.0162 0.1303
|
||||
7.0 0.0248 0.1389
|
||||
8.0 0.0374 0.1515
|
||||
9.0 0.0394 0.1535
|
||||
11.0 0.0336 0.1477
|
||||
13.0 0.0268 0.1409
|
||||
15.0 0.0201 0.1342
|
||||
17.0 0.0162 0.1303
|
||||
19.5 0.0144 0.1285
|
||||
21.5 0.0105 0.1246
|
||||
23.5 0.0086 0.1227
|
||||
25.5 0.0086 0.1227
|
||||
27.5 0.0077 0.1218
|
||||
29.5 0.0077 0.1218
|
||||
31.5 0.0058 0.1199
|
||||
33.5 0.0019 0.1160
|
||||
35.5 0.0048 0.1189
|
||||
37.5 0.0000 0.1141
|
||||
#zone,t="top wall"
|
||||
# -5.0 -0.1054 0.0087
|
||||
# 3.0 -0.0727 0.0414
|
||||
# 5.0 -0.0630 0.0511
|
||||
# 7.0 -0.0543 0.0598
|
||||
# 9.0 -0.0331 0.0810
|
||||
# 11.0 -0.0148 0.0993
|
||||
# 13.0 -0.0023 0.1118
|
||||
# 15.0 -0.0003 0.1138
|
||||
# 17.0 0.0065 0.1206
|
||||
# 19.0 0.0075 0.1216
|
||||
# 21.0 0.0075 0.1216
|
||||
# 23.0 0.0066 0.1207
|
||||
# 25.0 0.0095 0.1236
|
||||
# 27.0 0.0076 0.1217
|
||||
# 29.0 0.0067 0.1208
|
||||
# 31.0 0.0048 0.1189
|
||||
# 33.0 0.0019 0.1160
|
||||
# 35.0 0.0038 0.1179
|
||||
# 37.0 -0.0010 0.1131
|
||||
|
|
@ -0,0 +1,148 @@
|
|||
#!/bin/bash
|
||||
|
||||
# If you want to change the output format change
|
||||
# set term pngcairo dashed size 1024,768 font "Arial,12"
|
||||
# to
|
||||
# set terminal postscript eps color enhanced
|
||||
# and change the format in the output file name
|
||||
|
||||
gnuplot -persist << EOF
|
||||
|
||||
# Legend position
|
||||
# set key left top
|
||||
# Turn off legend
|
||||
# set nokey
|
||||
|
||||
#set term x11 0
|
||||
|
||||
set style line 1 lt 2 lc rgb "black" lw 1
|
||||
set style line 2 lt 2 lc rgb "black" lw 1 dt 2
|
||||
set style line 3 lt 2 lc rgb "black" lw 1 dt 3
|
||||
set style line 4 lt 2 lc rgb "blue" lw 1 dt 1
|
||||
set style line 5 lt 2 lc rgb "blue" lw 1 dt 2
|
||||
set style line 6 lt 2 lc rgb "blue" lw 1 dt 3
|
||||
|
||||
set term pngcairo dashed size 1024,768 font "Arial,12"
|
||||
set output "Cp.png"
|
||||
set key left top
|
||||
set title "Lower Wall"
|
||||
set xlabel "{/:Italic x}/{/:Italic H}"
|
||||
set ylabel "{/:Italic C_p}"
|
||||
set xrange [-5:30]
|
||||
set yrange [-0.25:0.1]
|
||||
#set logscale y 10; set format y '%.0e'
|
||||
|
||||
plot "./validationData/p_lowerWall.raw" using 1:(\$4/(0.5*2165.6)) ls 1 title 'CFD'with lines, "./validationData/cp_exp.dat" using 1:2 ls 1 title 'EXP'
|
||||
|
||||
reset
|
||||
#set term x11 1
|
||||
|
||||
set style line 1 lt 2 lc rgb "black" lw 1
|
||||
set style line 2 lt 2 lc rgb "black" lw 1 dt 2
|
||||
set style line 3 lt 2 lc rgb "black" lw 1 dt 3
|
||||
set style line 4 lt 2 lc rgb "blue" lw 1 dt 1
|
||||
set style line 5 lt 2 lc rgb "blue" lw 1 dt 2
|
||||
set style line 6 lt 2 lc rgb "blue" lw 1 dt 3
|
||||
|
||||
|
||||
set term pngcairo dashed size 1024,768 font "Arial,12"
|
||||
set output "Profile_1.png"
|
||||
set key left top
|
||||
set title "Line 1 ({/:Italic x}/{/:Italic H} = -4)"
|
||||
set xlabel "{/:Italic u}/{/:Italic U_{ref}}"
|
||||
set ylabel "({/:Italic y}-{/:Italic y_0})/{/:Italic H}"
|
||||
set xrange [-0.4:1.2]
|
||||
set yrange [0:3]
|
||||
#set logscale y 10; set format y '%.0e'
|
||||
|
||||
plot "./validationData/profile_1_U.xy" using (\$2/46.536):1 ls 1 title 'CFD'with lines, "./validationData/profile_1_exp.dat" using 3:(\$2-1) ls 1 title 'EXP'
|
||||
|
||||
reset
|
||||
#set term x11 2
|
||||
|
||||
set style line 1 lt 2 lc rgb "black" lw 1
|
||||
set style line 2 lt 2 lc rgb "black" lw 1 dt 2
|
||||
set style line 3 lt 2 lc rgb "black" lw 1 dt 3
|
||||
set style line 4 lt 2 lc rgb "blue" lw 1 dt 1
|
||||
set style line 5 lt 2 lc rgb "blue" lw 1 dt 2
|
||||
set style line 6 lt 2 lc rgb "blue" lw 1 dt 3
|
||||
|
||||
set term pngcairo dashed size 1024,768 font "Arial,12"
|
||||
set output "Profile_2.png"
|
||||
set key left top
|
||||
set title "Line 2 ({/:Italic x}/{/:Italic H} = 1)"
|
||||
set xlabel "{/:Italic u}/{/:Italic U_{ref}}"
|
||||
set ylabel "{/:Italic y}/{/:Italic H}"
|
||||
set xrange [-0.4:1.2]
|
||||
set yrange [0:3]
|
||||
#set logscale y 10; set format y '%.0e'
|
||||
|
||||
plot "./validationData/profile_2_U.xy" using (\$2/46.536):1 ls 1 title 'CFD'with lines, "./validationData/profile_2_exp.dat" using 3:2 ls 1 title 'EXP'
|
||||
|
||||
reset
|
||||
#set term x11 3
|
||||
|
||||
set style line 1 lt 2 lc rgb "black" lw 1
|
||||
set style line 2 lt 2 lc rgb "black" lw 1 dt 2
|
||||
set style line 3 lt 2 lc rgb "black" lw 1 dt 3
|
||||
set style line 4 lt 2 lc rgb "blue" lw 1 dt 1
|
||||
set style line 5 lt 2 lc rgb "blue" lw 1 dt 2
|
||||
set style line 6 lt 2 lc rgb "blue" lw 1 dt 3
|
||||
|
||||
set term pngcairo dashed size 1024,768 font "Arial,12"
|
||||
set output "Profile_3.png"
|
||||
set key left top
|
||||
set title "Line 3 ({/:Italic x}/{/:Italic H} = 4)"
|
||||
set xlabel "{/:Italic u}/{/:Italic U_{ref}}"
|
||||
set ylabel "{/:Italic y}/{/:Italic H}"
|
||||
set xrange [-0.4:1.2]
|
||||
set yrange [0:3]
|
||||
#set logscale y 10; set format y '%.0e'
|
||||
|
||||
plot "./validationData/profile_3_U.xy" using (\$2/46.536):1 ls 1 title 'CFD'with lines, "./validationData/profile_3_exp.dat" using 3:2 ls 1 title 'EXP'
|
||||
|
||||
reset
|
||||
#set term x11 4
|
||||
|
||||
set style line 1 lt 2 lc rgb "black" lw 1
|
||||
set style line 2 lt 2 lc rgb "black" lw 1 dt 2
|
||||
set style line 3 lt 2 lc rgb "black" lw 1 dt 3
|
||||
set style line 4 lt 2 lc rgb "blue" lw 1 dt 1
|
||||
set style line 5 lt 2 lc rgb "blue" lw 1 dt 2
|
||||
set style line 6 lt 2 lc rgb "blue" lw 1 dt 3
|
||||
|
||||
set term pngcairo dashed size 1024,768 font "Arial,12"
|
||||
set output "Profile_4.png"
|
||||
set key left top
|
||||
set title "Line 4 ({/:Italic x}/{/:Italic H} = 6)"
|
||||
set xlabel "{/:Italic u}/{/:Italic U_{ref}}"
|
||||
set ylabel "{/:Italic y}/{/:Italic H}"
|
||||
set xrange [-0.4:1.2]
|
||||
set yrange [0:3]
|
||||
#set logscale y 10; set format y '%.0e'
|
||||
|
||||
plot "./validationData/profile_4_U.xy" using (\$2/46.536):1 ls 1 title 'CFD'with lines, "./validationData/profile_4_exp.dat" using 3:2 ls 1 title 'EXP'
|
||||
|
||||
reset
|
||||
#set term x11 5
|
||||
|
||||
set style line 1 lt 2 lc rgb "black" lw 1
|
||||
set style line 2 lt 2 lc rgb "black" lw 1 dt 2
|
||||
set style line 3 lt 2 lc rgb "black" lw 1 dt 3
|
||||
set style line 4 lt 2 lc rgb "blue" lw 1 dt 1
|
||||
set style line 5 lt 2 lc rgb "blue" lw 1 dt 2
|
||||
set style line 6 lt 2 lc rgb "blue" lw 1 dt 3
|
||||
|
||||
set term pngcairo dashed size 1024,768 font "Arial,12"
|
||||
set output "Profile_5.png"
|
||||
set key left top
|
||||
set title "Line 5 ({/:Italic x}/{/:Italic H} = 10)"
|
||||
set xlabel "{/:Italic u}/{/:Italic U_{ref}}"
|
||||
set ylabel "{/:Italic y}/{/:Italic H}"
|
||||
set xrange [-0.4:1.2]
|
||||
set yrange [0:3]
|
||||
#set logscale y 10; set format y '%.0e'
|
||||
|
||||
plot "./validationData/profile_5_U.xy" using (\$2/46.536):1 ls 1 title 'CFD'with lines, "./validationData/profile_5_exp.dat" using 3:2 ls 1 title 'EXP'
|
||||
|
||||
EOF
|
||||
|
|
@ -0,0 +1,25 @@
|
|||
# Exp data drom Driver and Seegmiller,
|
||||
# y is y/H, u is u/Uref, v is v/Uref,
|
||||
# uu is 1000*u'u'/Uref**2
|
||||
# vv is 1000*v'v'/Uref**2
|
||||
# uv is 1000*u'v'/Uref**2
|
||||
# triple moments are 10000*u'u'u'/Uref**3 (for example)
|
||||
# variables="ID","y","u","v","uu","vv","uv","uuu","uvv","vuu","vvv"
|
||||
# zone,t="exp, x/H=-4"
|
||||
13 1.10 0.657 0.000 5.16 2.70 -1.02 -0.18 -0.14 0.12 0.22
|
||||
14 1.15 0.696 0.000 5.26 2.19 -1.01 -0.29 -0.12 0.11 0.22
|
||||
15 1.20 0.719 0.001 4.67 2.23 -1.07 -0.36 -0.13 0.19 0.16
|
||||
16 1.30 0.760 0.000 3.86 2.32 -1.07 -0.36 -0.17 0.18 0.14
|
||||
17 1.40 0.790 0.001 3.49 2.10 -0.97 -0.17 -0.03 0.02 0.07
|
||||
18 1.50 0.818 -0.001 3.06 2.05 -0.92 -0.15 -0.12 0.11 0.13
|
||||
19 1.70 0.870 0.000 2.58 1.64 -0.71 -0.25 -0.13 0.14 0.13
|
||||
20 2.00 0.926 0.000 1.80 1.18 -0.44 -0.23 -0.10 0.10 0.11
|
||||
21 2.40 0.982 0.003 0.86 0.63 -0.10 -0.12 -0.06 0.05 0.04
|
||||
22 2.80 1.003 0.007 0.42 0.32 0.02 0.02 0.01 0.00 0.00
|
||||
23 3.20 1.003 0.006 0.39 0.30 0.03 0.03 0.02 0.00 0.01
|
||||
24 3.60 1.000 0.006 0.37 0.28 0.04 0.02 0.01 0.00 0.00
|
||||
25 4.00 1.001 0.003 0.37 0.28 0.03 0.02 0.01 -0.01 0.00
|
||||
26 5.00 1.001 0.000 0.37 0.27 0.03 0.02 0.01 -0.01 0.00
|
||||
27 6.00 1.000 -0.004 0.36 0.26 0.04 0.01 0.01 -0.01 0.00
|
||||
28 7.50 1.001 -0.006 0.37 0.28 0.04 0.02 0.01 0.00 0.00
|
||||
29 8.20 0.943 0.002 1.60 1.28 0.46 -0.28 -0.17 -0.13 -0.08
|
||||
|
|
@ -0,0 +1,36 @@
|
|||
# Exp data drom Driver and Seegmiller,
|
||||
# y is y/H, u is u/Uref, v is v/Uref,
|
||||
# uu is 1000*u'u'/Uref**2
|
||||
# vv is 1000*v'v'/Uref**2
|
||||
# uv is 1000*u'v'/Uref**2
|
||||
# triple moments are 10000*u'u'u'/Uref**3 (for example)
|
||||
# variables="ID","y","u","v","uu","vv","uv","uuu","uvv","vuu","vvv"
|
||||
# zone,t="exp, x/H=1"
|
||||
2 0.10 -0.008 -0.005 3.18 2.94 -0.72 0.00 -0.35 0.14 0.19
|
||||
3 0.15 -0.017 -0.004 4.71 3.22 0.46 -0.46 -0.14 0.15 0.06
|
||||
4 0.20 -0.027 -0.001 5.04 3.29 0.21 -0.38 -0.11 0.22 -0.01
|
||||
5 0.30 -0.035 0.004 5.12 3.49 -0.14 -0.31 -0.21 0.31 0.05
|
||||
6 0.40 -0.035 0.014 4.50 3.67 -0.43 -0.33 -0.24 0.33 0.31
|
||||
7 0.50 -0.035 0.022 4.94 2.90 -0.65 -0.48 -0.25 0.32 0.20
|
||||
8 0.60 -0.027 0.028 4.56 3.00 -0.74 -0.61 -0.23 0.28 0.18
|
||||
9 0.70 -0.005 0.025 5.03 3.44 -0.92 -0.39 -0.15 -0.10 -0.04
|
||||
10 0.80 0.081 0.006 9.24 7.38 -3.71 2.63 2.22 -2.73 -2.24
|
||||
11 0.90 0.292 -0.027 13.80 11.77 -7.05 1.18 0.94 -0.65 -1.00
|
||||
12 1.00 0.519 -0.033 13.93 10.41 -5.85 -4.18 -2.73 2.91 2.32
|
||||
13 1.10 0.704 -0.023 6.18 3.79 -1.54 -1.19 -0.49 0.45 0.62
|
||||
14 1.15 0.740 -0.022 4.87 2.84 -1.25 -0.07 -0.16 0.12 0.25
|
||||
15 1.20 0.762 -0.023 4.27 2.51 -1.05 -0.21 -0.11 0.09 0.24
|
||||
16 1.30 0.793 -0.031 3.87 2.26 -0.98 -0.22 -0.11 0.14 0.16
|
||||
17 1.40 0.827 -0.027 3.38 2.05 -0.82 -0.26 -0.09 0.13 0.15
|
||||
18 1.50 0.851 -0.031 3.08 1.98 -0.81 -0.26 -0.11 0.14 0.13
|
||||
19 1.70 0.891 -0.034 2.49 1.57 -0.64 -0.25 -0.11 0.09 0.09
|
||||
20 2.00 0.942 -0.033 1.76 1.21 -0.44 -0.24 -0.09 0.11 0.09
|
||||
21 2.40 0.990 -0.031 0.91 0.67 -0.12 -0.11 -0.04 0.05 0.04
|
||||
22 2.80 1.009 -0.028 0.45 0.36 0.01 0.01 0.01 0.00 0.01
|
||||
23 3.20 1.003 -0.028 0.41 0.33 0.04 0.02 0.02 0.00 0.00
|
||||
24 3.60 1.002 -0.026 0.40 0.33 0.03 0.01 0.01 0.00 0.00
|
||||
25 4.00 1.000 -0.028 0.38 0.30 0.03 0.02 0.01 0.00 0.00
|
||||
26 5.00 0.997 -0.027 0.36 0.29 0.01 0.02 0.01 0.00 0.00
|
||||
27 6.00 0.994 -0.027 0.33 0.27 -0.01 0.01 0.01 -0.01 -0.01
|
||||
28 7.50 0.993 -0.025 0.36 0.30 0.01 0.01 0.01 0.00 0.00
|
||||
29 8.20 0.919 -0.013 1.83 1.56 0.59 -0.30 -0.19 -0.15 -0.14
|
||||
|
|
@ -0,0 +1,36 @@
|
|||
# Exp data drom Driver and Seegmiller,
|
||||
# y is y/H, u is u/Uref, v is v/Uref,
|
||||
# uu is 1000*u'u'/Uref**2
|
||||
# vv is 1000*v'v'/Uref**2
|
||||
# uv is 1000*u'v'/Uref**2
|
||||
# triple moments are 10000*u'u'u'/Uref**3 (for example)
|
||||
# variables="ID","y","u","v","uu","vv","uv","uuu","uvv","vuu","vvv"
|
||||
# zone,t="exp, x/H=4"
|
||||
2 0.10 -0.171 -0.014 9.71 7.67 -2.64 4.26 2.49 -2.60 -2.47
|
||||
3 0.15 -0.154 -0.020 13.20 10.10 -4.45 10.50 4.70 -5.18 -4.94
|
||||
4 0.20 -0.134 -0.027 16.50 12.98 -6.30 13.67 5.61 -6.52 -5.34
|
||||
5 0.30 -0.057 -0.046 25.91 16.16 -9.94 16.29 4.72 -5.66 -4.41
|
||||
6 0.40 0.052 -0.068 32.92 18.95 -12.02 -4.83 -2.74 3.50 2.13
|
||||
7 0.50 0.169 -0.082 30.67 18.32 -11.11 -26.75 -8.98 11.86 6.02
|
||||
8 0.60 0.279 -0.093 23.24 15.90 -8.99 -19.61 -9.14 9.42 6.90
|
||||
9 0.70 0.381 -0.094 18.21 13.03 -6.93 -5.06 -4.23 3.71 4.03
|
||||
10 0.80 0.473 -0.097 17.03 10.98 -5.98 -3.65 -3.38 3.11 3.89
|
||||
11 0.90 0.568 -0.097 14.45 8.30 -4.52 -5.19 -2.98 2.45 3.04
|
||||
12 1.00 0.656 -0.092 11.55 6.33 -3.48 -6.85 -2.78 2.86 2.55
|
||||
13 1.10 0.730 -0.089 7.70 4.20 -2.10 -3.69 -1.47 1.52 1.57
|
||||
14 1.15 0.754 -0.087 6.27 3.52 -1.65 -2.27 -0.94 1.00 1.14
|
||||
15 1.20 0.778 -0.084 4.85 2.83 -1.19 -0.84 -0.42 0.48 0.71
|
||||
16 1.30 0.811 -0.083 3.80 2.19 -0.87 -0.44 -0.19 0.16 0.32
|
||||
17 1.40 0.835 -0.081 3.17 1.80 -0.67 -0.29 -0.11 0.13 0.20
|
||||
18 1.50 0.856 -0.079 2.86 1.52 -0.56 -0.31 -0.08 0.10 0.14
|
||||
19 1.70 0.895 -0.076 2.27 1.25 -0.45 -0.34 -0.10 0.14 0.13
|
||||
20 2.00 0.940 -0.072 1.47 0.78 -0.23 -0.22 -0.07 0.07 0.07
|
||||
21 2.40 0.978 -0.063 0.61 0.39 0.00 -0.06 -0.01 0.03 0.02
|
||||
22 2.80 0.986 -0.056 0.38 0.26 0.06 0.01 0.00 -0.01 0.00
|
||||
23 3.20 0.984 -0.050 0.36 0.23 0.05 0.01 0.00 0.00 0.00
|
||||
24 3.60 0.985 -0.047 0.36 0.23 0.07 0.01 0.00 0.00 0.00
|
||||
25 4.00 0.985 -0.043 0.33 0.21 0.06 0.00 0.00 0.00 0.00
|
||||
26 5.00 0.983 -0.037 0.33 0.22 0.06 0.01 0.00 0.00 0.00
|
||||
27 6.00 0.980 -0.031 0.34 0.23 0.06 0.01 0.00 0.00 0.00
|
||||
28 7.50 0.980 -0.025 0.35 0.25 0.07 0.01 0.00 -0.01 0.00
|
||||
29 8.20 0.892 -0.013 2.23 1.59 0.65 -0.40 -0.22 -0.19 -0.13
|
||||
|
|
@ -0,0 +1,36 @@
|
|||
# Exp data drom Driver and Seegmiller,
|
||||
# y is y/H, u is u/Uref, v is v/Uref,
|
||||
# uu is 1000*u'u'/Uref**2
|
||||
# vv is 1000*v'v'/Uref**2
|
||||
# uv is 1000*u'v'/Uref**2
|
||||
# triple moments are 10000*u'u'u'/Uref**3 (for example)
|
||||
# variables="ID","y","u","v","uu","vv","uv","uuu","uvv","vuu","vvv"
|
||||
# zone,t="exp, x/H=6"
|
||||
2 0.10 0.147 -0.052 17.58 12.81 -4.10 0.03 0.65 -0.11 -2.05
|
||||
3 0.15 0.169 -0.061 19.50 15.02 -5.48 0.04 0.24 0.11 -1.31
|
||||
4 0.20 0.188 -0.067 19.53 16.50 -6.79 -2.21 -1.70 1.31 0.58
|
||||
5 0.30 0.248 -0.084 20.73 16.98 -7.86 -5.45 -4.28 3.36 3.25
|
||||
6 0.40 0.306 -0.091 19.70 16.60 -8.23 -6.80 -5.78 4.97 5.60
|
||||
7 0.50 0.369 -0.097 20.03 15.89 -8.24 -9.08 -7.74 7.14 6.87
|
||||
8 0.60 0.435 -0.099 18.36 13.86 -7.36 -8.93 -6.86 6.50 6.48
|
||||
9 0.70 0.501 -0.095 17.26 12.00 -6.57 -10.43 -7.51 7.32 7.26
|
||||
10 0.80 0.569 -0.094 14.96 9.63 -5.06 -10.37 -5.91 6.20 5.55
|
||||
11 0.90 0.632 -0.091 11.96 7.62 -3.79 -10.19 -4.74 5.11 4.09
|
||||
12 1.00 0.686 -0.089 8.88 5.57 -2.80 -7.74 -3.45 3.55 2.93
|
||||
13 1.10 0.730 -0.086 6.11 3.95 -1.76 -3.81 -1.72 1.93 1.87
|
||||
14 1.15 0.747 -0.086 5.23 3.35 -1.51 -2.68 -1.23 1.39 1.33
|
||||
15 1.20 0.763 -0.085 4.36 2.75 -1.25 -1.54 -0.74 0.84 0.79
|
||||
16 1.30 0.794 -0.083 3.40 2.19 -0.80 -0.47 -0.23 0.25 0.44
|
||||
17 1.40 0.815 -0.082 2.87 1.79 -0.65 -0.33 -0.12 0.20 0.25
|
||||
18 1.50 0.834 -0.081 2.50 1.53 -0.61 -0.29 -0.10 0.14 0.14
|
||||
19 1.70 0.869 -0.077 2.01 1.24 -0.43 -0.27 -0.10 0.13 0.13
|
||||
20 2.00 0.915 -0.074 1.28 0.80 -0.21 -0.22 -0.06 0.08 0.04
|
||||
21 2.40 0.949 -0.066 0.53 0.44 -0.04 -0.01 0.00 0.01 0.02
|
||||
22 2.80 0.954 -0.060 0.42 0.33 0.03 0.02 0.02 0.00 0.00
|
||||
23 3.20 0.956 -0.054 0.41 0.33 0.03 0.02 0.02 0.01 0.00
|
||||
24 3.60 0.959 -0.050 0.39 0.31 0.03 0.03 0.02 0.00 0.00
|
||||
25 4.00 0.961 -0.046 0.38 0.29 0.02 0.03 0.02 0.00 0.01
|
||||
26 5.00 0.963 -0.040 0.38 0.30 0.02 0.02 0.02 0.00 0.00
|
||||
27 6.00 0.965 -0.035 0.38 0.31 0.02 0.02 0.02 0.00 -0.01
|
||||
28 7.50 0.969 -0.027 0.42 0.34 0.03 0.02 0.01 0.00 0.00
|
||||
29 8.20 0.870 -0.018 2.75 1.84 0.76 -0.55 -0.20 -0.20 -0.11
|
||||
|
|
@ -0,0 +1,36 @@
|
|||
# Exp data drom Driver and Seegmiller,
|
||||
# y is y/H, u is u/Uref, v is v/Uref,
|
||||
# uu is 1000*u'u'/Uref**2
|
||||
# vv is 1000*v'v'/Uref**2
|
||||
# uv is 1000*u'v'/Uref**2
|
||||
# triple moments are 10000*u'u'u'/Uref**3 (for example)
|
||||
# variables="ID","y","u","v","uu","vv","uv","uuu","uvv","vuu","vvv"
|
||||
# zone,t="exp, x/H=10"
|
||||
2 0.10 0.418 -0.029 8.89 5.63 -1.43 2.05 0.73 -0.34 -0.46
|
||||
3 0.15 0.433 -0.027 9.97 6.15 -1.84 2.59 0.74 -0.35 -1.10
|
||||
4 0.20 0.447 -0.032 10.51 6.97 -2.35 2.19 0.83 -0.25 -0.99
|
||||
5 0.30 0.471 -0.036 11.57 8.34 -3.43 1.23 0.55 -0.24 -0.37
|
||||
6 0.40 0.496 -0.041 12.24 9.08 -4.19 0.67 -0.37 0.07 0.46
|
||||
7 0.50 0.526 -0.043 12.86 9.16 -4.51 -1.12 -1.29 0.99 1.47
|
||||
8 0.60 0.560 -0.045 13.37 9.05 -4.61 -3.54 -2.59 2.08 2.43
|
||||
9 0.70 0.596 -0.047 13.03 8.44 -4.61 -5.62 -3.18 3.10 3.00
|
||||
10 0.80 0.638 -0.047 12.00 7.75 -4.30 -7.23 -4.00 3.65 3.93
|
||||
11 0.90 0.676 -0.046 10.96 6.44 -3.68 -8.76 -3.81 3.96 3.37
|
||||
12 1.00 0.712 -0.047 9.04 5.35 -2.94 -8.54 -3.51 3.72 2.94
|
||||
13 1.10 0.745 -0.046 7.23 4.32 -2.33 -6.66 -2.57 2.72 2.36
|
||||
14 1.15 0.758 -0.045 6.44 3.92 -2.02 -5.85 -2.28 2.36 2.02
|
||||
15 1.20 0.772 -0.044 5.66 3.52 -1.71 -5.04 -2.00 2.01 1.68
|
||||
16 1.30 0.798 -0.042 4.22 2.80 -1.36 -2.94 -1.28 1.48 1.38
|
||||
17 1.40 0.819 -0.040 3.29 2.21 -0.96 -1.61 -0.81 0.80 0.96
|
||||
18 1.50 0.839 -0.039 2.54 1.80 -0.72 -0.83 -0.53 0.51 0.65
|
||||
19 1.70 0.871 -0.038 1.60 1.03 -0.35 -0.28 -0.09 0.11 0.15
|
||||
20 2.00 0.908 -0.034 0.96 0.66 -0.14 -0.13 -0.03 0.08 0.09
|
||||
21 2.40 0.929 -0.032 0.43 0.35 0.02 0.02 0.02 0.02 0.03
|
||||
22 2.80 0.933 -0.031 0.35 0.27 0.04 0.03 0.02 0.01 0.02
|
||||
23 3.20 0.937 -0.032 0.33 0.25 0.04 0.04 0.03 0.02 0.03
|
||||
24 3.60 0.939 -0.030 0.32 0.24 0.04 0.04 0.03 0.02 0.02
|
||||
25 4.00 0.940 -0.031 0.31 0.22 0.02 0.05 0.02 0.01 0.01
|
||||
26 5.00 0.944 -0.029 0.30 0.22 0.03 0.04 0.02 0.01 0.01
|
||||
27 6.00 0.947 -0.028 0.30 0.23 0.03 0.03 0.02 0.01 0.02
|
||||
28 7.50 0.949 -0.024 0.38 0.32 0.05 0.00 0.01 0.01 0.01
|
||||
29 8.20 0.827 -0.019 3.14 2.18 0.95 -0.49 -0.25 -0.20 -0.16
|
||||
|
|
@ -0,0 +1,18 @@
|
|||
Author: Robert Keser, robert.keser@gmail.com
|
||||
|
||||
Reference experimental data:
|
||||
@article{DRIVER-1985,
|
||||
doi = {10.2514/3.8890},
|
||||
title = {Features of a reattaching turbulent shear layer in divergent
|
||||
channelflow},
|
||||
author = {D. M. Driver and H. L. Seegmiller},
|
||||
publisher = {American Institute of Aeronautics and Astronautics},
|
||||
journal = {AIAA Journal},
|
||||
issnp = {0001-1452},
|
||||
issne = {1533-385X},
|
||||
year = {1985},
|
||||
month = {02},
|
||||
volume = {23},
|
||||
issue = {2},
|
||||
page = {163--171}
|
||||
}
|
||||
Reference in a new issue