View Issue Details
| ID | Project | Category | View Status | Date Submitted | Last Update |
|---|---|---|---|---|---|
| 0003119 | OpenFOAM | Bug | public | 2018-12-04 16:26 | 2018-12-05 00:19 |
| Reporter | weyburne | Assigned To | wyldckat | ||
| Priority | normal | Severity | major | Reproducibility | always |
| Status | closed | Resolution | no change required | ||
| Platform | blueCFD-Core-2017 | OS | win10 | ||
| Summary | 0003119: paraFoam reporting an incorrect pressure and velocity value at the multi-grading boundary | ||||
| Description | The problem occurs with paraFoam (paraView 5.4.1, Win 10) when requesting pressure and velocity slices through a mult-graded region. In the slices, the one data point at the the boundary between the multi-grading areas is slightly incorrect. It maybe an interpolation error, not sure. | ||||
| Steps To Reproduce | The problem can be seen using a simple 2-D channel, one block geometry simulation. BlockMesh is used with the multi-grading option used in both the flow direction, x-plane, and normal to the channel walls, y-plane. Air flow at RT is used with the channel length and the air velocity chosen to keep the flow laminar (Rex=5e5 at exit). SimpleFoam with the laminar option is run until P and U residuals are less than 2x10^7. paraFoam is then run and center line slices in both flow direction and wall normal are obtained. Attached files include all relevant OpenFoam files as well as two jpg graphs of the pressure illustrating the problem, In simple pressure or velocity plots, the problem is difficult to see unless the boundary region is blown up. It is very noticeable if one takes the derivative of the pressure or velocity. | ||||
| Tags | No tags attached. | ||||
|
|
U (1,372 bytes)
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 1.4.1-dev |
| \\ / A nd | Web: http://www.openfoam.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
root "/home/mjaworsk/OpenFOAM/mjaworsk-1.4.1-dev/run/conjugateTray2D";
case "lithium-5mm";
instance "0";
local "";
class volVectorField;
object U;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 1 -1 0 0 0 0];
internalField uniform (0.9375 0 0);
boundaryField
{
inlet
{
type fixedValue;
value uniform (0.9375 0 0);
}
outlet
{
type zeroGradient;
}
plate
{
type fixedValue;
value uniform (0 0 0);
}
frontAndBack
{
type empty;
}
}
// ************************************************************************* //
fvSolution (2,137 bytes)
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 5 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
location "system";
object fvSolution;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
solvers
{
p
{
solver GAMG;
tolerance 1e-010;
relTol 0.1;
smoother GaussSeidel;
nPreSweeps 1;
nPostSweeps 3;
nFinestSweeps 3;
scaleCorrection true;
directSolveCoarsest off;
catchAgglomeration on;
nCellsInCoarstLevel 30;
agglomerator faceAreaPair;
mergeLevels 1;
}
pFinal
{
$p;
relTol 0;
}
"(U|k|epsilon|omega|f|v2)"
{
solver smoothSolver;
smoother symGaussSeidel;
tolerance 1e-08;
relTol 0.1;
}
}
SIMPLE
{
nNonOrthogonalCorrectors 0;
consistent yes;
residualControl
{
p 1e-10;
U 1e-8;
"(k|epsilon|omega|f|v2)" 1e-8;
}
}
relaxationFactors
{
equations
{
P 0.8
U 0.95; // 0.9 is more stable but 0.95 more convergent
".*" 0.8; // 0.9 is more stable but 0.95 more convergent
}
}
PISO
{
nCorrectors 2;
nNonOrthogonalCorrectors 0;
pRefCell 0;
pRefValue 0;
}
potentialFlow
{
nNonOrthogonalCorrectors 10;
}
// ************************************************************************* //
turbulenceProperties (1,093 bytes)
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 5 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
location "constant";
object turbulenceProperties;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
simulationType laminar;
RAS
{
// Tested with kEpsilon, realizableKE, kOmega, kOmegaSST, v2f,
// ShihQuadraticKE, LienCubicKE.
RASModel kEpsilon;
turbulence off;
printCoeffs off;
}
// ************************************************************************* //
transportProperties (1,035 bytes)
/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 1.3 |
| \\ / A nd | Web: http://www.openfoam.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
root "";
case "";
instance "";
local "";
class dictionary;
object transportProperties;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
transportModel Newtonian;
nu nu [0 2 -1 0 0 0 0] 1.5e-5;
// ************************************************************************* //
fvSchemes (1,555 bytes)
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 5 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
location "system";
object fvSchemes;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
ddtSchemes
{
default steadyState;
}
gradSchemes
{
default Gauss linear;
}
divSchemes
{
default none;
div(phi,U) bounded Gauss linearUpwind grad(U);
div(phi,k) bounded Gauss limitedLinear 1;
div(phi,epsilon) bounded Gauss limitedLinear 1;
div(phi,omega) bounded Gauss limitedLinear 1;
div(phi,v2) bounded Gauss limitedLinear 1;
div((nuEff*dev2(T(grad(U))))) Gauss linear;
div(nonlinearStress) Gauss linear;
}
laplacianSchemes
{
default Gauss linear corrected;
}
interpolationSchemes
{
default linear;
}
snGradSchemes
{
default corrected;
}
wallDist
{
method meshWave;
}
// ************************************************************************* //
controlDict (1,359 bytes)
/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 1.3 |
| \\ / A nd | Web: http://www.openfoam.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
root "";
case "";
instance "";
local "";
class dictionary;
object controlDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
application simpleFoam;
startFrom startTime;
startTime 0;
stopAt endTime;
endTime 12.0;
deltaT 0.002;
writeControl runTime;
writeInterval 3.0;
purgeWrite 0;
writeFormat ascii;
writePrecision 15;
writeCompression uncompressed;
timeFormat general;
timePrecision 6;
runTimeModifiable yes;
functions
{
#includeFunc residuals
}
// ************************************************************************* //
blockMeshDict (1,746 bytes)
/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 1.3 |
| \\ / A nd | Web: http://www.openfoam.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
root "";
case "";
instance "";
local "";
class dictionary;
object blockMeshDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
convertToMeters 0.001;
vertices
(
// Block 1
(0 0 0) //0
(8000 0 0) //1
(8000 1000 0) //2
(0 1000 0) //3
(0 0 1) //4
(8000 0 1) //5
(8000 1000 1) //6
(0 1000 1) //7
);
blocks
(
// Block 1
hex (0 1 2 3 4 5 6 7) (1000 600 1)
simpleGrading
(
((0.05 0.3 10) (0.9 0.4 1) (0.05 0.3 0.1))
((0.05 0.3 10) (0.9 0.4 1) (0.05 0.3 0.1))
1
)
);
edges
(
);
patches
(
wall inlet
(
(0 4 7 3)
)
wall outlet
(
(1 5 6 2)
)
wall plate
(
(7 6 2 3)
(0 1 5 4)
)
empty frontAndBack
(
(0 1 2 3)
(4 5 6 7)
)
);
// ************************************************************************* //
p (1,351 bytes)
/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 1.3 |
| \\ / A nd | Web: http://www.openfoam.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
root "/home/mjaworsk/OpenFOAM/mjaworsk-1.4.1-dev/run/conjugateTray2D";
case "lithium-5mm";
instance "0";
local "";
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;
}
plate
{
type zeroGradient;
}
frontAndBack
{
type empty;
}
}
// ************************************************************************* //
|
|
|
The blueCFD-Core project has its own issue tracker: https://github.com/blueCFD/Core/issues Please do not use the OpenFOAM Foundation's issue tracker to report problems that occur with independent projects, unless of course the same issue was reproduced with a version that is provided by the OpenFOAM Foundation. |
| Date Modified | Username | Field | Change |
|---|---|---|---|
| 2018-12-04 16:26 | weyburne | New Issue | |
| 2018-12-04 16:26 | weyburne | File Added: U | |
| 2018-12-04 16:26 | weyburne | File Added: fvSolution | |
| 2018-12-04 16:26 | weyburne | File Added: turbulenceProperties | |
| 2018-12-04 16:26 | weyburne | File Added: transportProperties | |
| 2018-12-04 16:26 | weyburne | File Added: fvSchemes | |
| 2018-12-04 16:26 | weyburne | File Added: controlDict | |
| 2018-12-04 16:26 | weyburne | File Added: blockMeshDict | |
| 2018-12-04 16:26 | weyburne | File Added: p | |
| 2018-12-04 16:26 | weyburne | File Added: GraphPy4.jpg | |
| 2018-12-05 00:19 | wyldckat | Note Added: 0010204 | |
| 2018-12-05 00:19 | wyldckat | Assigned To | => wyldckat |
| 2018-12-05 00:19 | wyldckat | Status | new => closed |
| 2018-12-05 00:19 | wyldckat | Resolution | open => no change required |
