 |
 |
 |
|
|
|
[Sponsors] | ||||
October 5, 2016, 05:19
|
LES case diverges after 5 flow-through-time
|
#1 |
|
New Member
James Zh
Join Date: Mar 2016
Posts: 3
Rep Power: 10  |
Hi Foamers,
This is my first post  I'm doing Ma7 LES on rhoCentralFoam with kEqn model. The geometry is quite simple. The mapped inflow boundary condition is used as in the attached geometry. The back and forth patches are cyclic, while the top one is supersonicfreestream. details of physical and numerical boundary conditions: Code:
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: v3.0+ |
| \\ / A nd | Web: www.OpenFOAM.com |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class polyBoundaryMesh;
location "constant/polyMesh";
object boundary;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
6
(
INLET
{
type mappedPatch;
inGroups 1(mappedPatch);
nFaces 5355;
startFace 5836184;
sampleMode nearestCell;
sampleRegion region0;
samplePatch none;
offsetMode uniform;
offset (0.06 0 0);
}
SYM
{
type cyclic;
inGroups 1(cyclic);
nFaces 31195;
startFace 5841539;
matchTolerance 0.0006;
transform unknown;
neighbourPatch SYM2;
}
SYM2
{
type cyclic;
inGroups 1(cyclic);
nFaces 31195;
startFace 5872734;
matchTolerance 0.0006;
transform unknown;
neighbourPatch SYM;
}
WALL_D
{
type wall;
inGroups 1(wall);
nFaces 23121;
startFace 5903929;
}
OUTLET
{
type patch;
nFaces 5355;
startFace 5927050;
}
FAROUT
{
type patch;
nFaces 23121;
startFace 5932405;
}
)
// ************************************************************************* //
Code:
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: v3.0+ |
| \\ / A nd | Web: www.OpenFOAM.com |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
location "0";
object alphat;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [1 -1 -1 0 0 0 0];
internalField uniform 0;
boundaryField
{
INLET
{
type calculated;
value uniform 0;
}
FAROUT
{
type calculated;
value uniform 0;
}
OUTLET
{
type calculated;
value uniform 0;
}
SYM
{
type cyclic;
}
SYM2
{
type cyclic;
}
WALL_D
{
type zeroGradient;
}
}
// ************************************************************************* //
Code:
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: v3.0+ |
| \\ / A nd | Web: www.OpenFOAM.com |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
location "0";
object k;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 2 -2 0 0 0 0];
internalField uniform 4925.8;
boundaryField
{
INLET
{
type mapped;
value uniform 4925.8;
interpolationScheme cell;
setAverage true;
average 4925.8;
}
FAROUT
{
type zeroGradient;
}
OUTLET
{
type zeroGradient;
}
SYM
{
type cyclic;
}
SYM2
{
type cyclic;
}
WALL_D
{
type zeroGradient;
}
}
// ************************************************************************* //
Code:
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: v3.0+ |
| \\ / A nd | Web: www.OpenFOAM.com |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
location "0";
object nut;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 2 -1 0 0 0 0];
internalField uniform 0;
boundaryField
{
INLET
{
type zeroGradient;
}
FAROUT
{
type zeroGradient;
}
OUTLET
{
type zeroGradient;
}
SYM
{
type cyclic;
}
SYM2
{
type cyclic;
}
WALL_D
{
type zeroGradient;
}
}
// ************************************************************************* //
Code:
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: v3.0+ |
| \\ / A nd | Web: www.OpenFOAM.com |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
object p;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [1 -1 -2 0 0 0 0];
internalField uniform 1341.28589;
boundaryField
{
INLET
{
type fixedValue;
value uniform 1341.28589;
}
FAROUT
{
type zeroGradient;
}
OUTLET
{
type waveTransmissive;
field p;
phi phi;
rho rho;
psi thermo:psi;
gamma 1.4;
fieldInf 1341.28589;
lInf 10;
value uniform 1341.28589;
}
SYM
{
type cyclic;
}
SYM2
{
type cyclic;
}
WALL_D
{
type zeroGradient;
}
}
// ************************************************************************* //
Code:
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: v3.0+ |
| \\ / A nd | Web: www.OpenFOAM.com |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
object T;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 0 0 1 0 0 0];
internalField uniform 62.9;
boundaryField
{
INLET
{
type fixedValue;
value uniform 62.9;
}
FAROUT
{
type zeroGradient;
}
OUTLET
{
type zeroGradient;
}
SYM
{
type cyclic;
}
SYM2
{
type cyclic;
}
WALL_D
{
type zeroGradient;
}
}
// ************************************************************************* //
Code:
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: v3.0+ |
| \\ / A nd | Web: www.OpenFOAM.com |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volVectorField;
object U;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
dimensions [0 1 -1 0 0 0 0];
internalField uniform (1146.1 0 0);
boundaryField
{
INLET
{
type mapped;
value uniform (1146.1 0 0);
interpolationScheme cell;
setAverage false;
average (1146.1 0 0);
}
FAROUT
{
type supersonicFreestream;
pInf 1341.28589;
TInf 62.9;
UInf (1146.1 0 0);
gamma 1.4;
value uniform (1146.1 0 0);
}
OUTLET
{
type pressureInletOutletVelocity;
inletValue uniform (0 0 0);
value uniform (0 0 0);
}
SYM
{
type cyclic;
}
SYM2
{
type cyclic;
}
WALL_D
{
type fixedValue;
value uniform (0 0 0);
}
}
// ************************************************************************* //
Code:
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: v3.0+ |
| \\ / A nd | Web: www.OpenFOAM.com |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
location "system";
object fvSchemes;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
fluxScheme Kurganov;
ddtSchemes
{
default backward;
}
gradSchemes
{
default Gauss linear;
}
divSchemes
{
default none;
div(tauMC) Gauss linear;
div(phi,k) Gauss Gamma 0.2;
}
laplacianSchemes
{
default Gauss linear corrected;
}
interpolationSchemes
{
default none;
reconstruct(rho) vanLeer;
reconstruct(U) vanLeerV;
reconstruct(T) vanLeer;
interpolate(rhoU) linear;
interpolate(muEff) linear;
interpolate(tauMC) linear;
interpolate(rho) linear;
}
snGradSchemes
{
default corrected;
}
// ************************************************************************* //
Code:
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: v3.0+ |
| \\ / A nd | Web: www.OpenFOAM.com |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
location "system";
object fvSolution;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
solvers
{
"(rho|rhoU|rhoE)"
{
solver diagonal;
}
U
{
solver smoothSolver;
smoother GaussSeidel;
nSweeps 2;
tolerance 1e-07;
relTol 0;
}
e
{
$U;
tolerance 1e-7;
relTol 0;
}
k
{
$U;
tolerance 1e-7;
relTol 0;
}
}
// ************************************************************************* //
I monitored the mass flow rate of all the patches and it seemed OK with the inlet, outlet and the freestream one. The mass flow rate of cyclic patch kept increasing up to around 0.01 of the oulet's, is it common? Maybe because the turbulent boundary layer is still developing, i am not sure. The case has been running for about two days(around 5 flow-through-time), finally blow up...   Delta t is set to 1e-8 and the maxCo is less than 0.2 when it diverged. I have stuck here for a week... Any help will be appreciated. Regards James Zh Last edited by mcshidada; October 5, 2016 at 08:01. |
|
|
|
|
|
October 5, 2016, 07:36
|
|
#2 |
|
New Member
James Zh
Join Date: Mar 2016
Posts: 3
Rep Power: 10 |
Any comment is welcomed
|
|
|
|
|
|
|
July 17, 2018, 22:50
|
|
#3 | |
|
New Member
Ontario
Join Date: Jul 2018
Posts: 3
Rep Power: 8 |
Quote:
|
||
|
|
|
||
 |
| Tags |
| high-speed, les, rhocentralfoam |
|
|
Similar Threads
|
||||
| Thread | Thread Starter | Forum | Replies | Last Post |
| Extrusion with OpenFoam problem No. Iterations 0 | Lord Kelvin | OpenFOAM Running, Solving & CFD | 8 | March 28, 2016 12:08 |
| Moving mesh | Niklas Wikstrom (Wikstrom) | OpenFOAM Running, Solving & CFD | 122 | June 15, 2014 07:20 |
| Same SimpleFOAM Case converges with openFOAM 2.1 but diverges with openFOAM 2.0.1 | alsdia | OpenFOAM Running, Solving & CFD | 3 | October 22, 2012 12:25 |
| Orifice Plate with a fully developed flow - Problems with convergence | jonmec | OpenFOAM Running, Solving & CFD | 3 | July 28, 2011 06:24 |
| Could anybody help me see this error and give help | liugx212 | OpenFOAM Running, Solving & CFD | 3 | January 4, 2006 19:07 |
Linear Mode
