[random_ran]

Hi there:

I have two questions regarding with the relation between boundary/intial
conditions and the global matrix.

Question1:What is the size of Matrix A (see attached file, from Moukalled 2015)?
Question2: How many unknown number are there?

To give you more information with this question, I want use the
benchmark example: flow over a circular cylinder.

-version: OpenFOAM v4.1.
-application: pisoFoam
-geometry: see attached picture (Yeon, Seong Mo 2014 PHD Thesis)
-checkMeshResut: see attached file
-boundary/initial condition: see attached file
-turbulent model: LES,kEqn see attached file

I also think time step should be included in the discussion. Here, I
just want to make it sample. Let's assume we just finish one time
step, 0.001s.



-P.S. please let me known if you need further information.

Best regards,

Ran

Attach file:

matrix.png

coord.png

geo.png


/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 4.1 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
Build : 4.1
Exec : checkMesh
Date : Jun 28 2017
Time : 17:53:43
Host :
PID : 15901
Case :
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Allowing user-supplied system call operations

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create polyMesh for time = 0

Time = 0

Mesh stats
points: 5100480
faces: 15045220
internal faces: 14790380
cells: 4972600
faces per cell: 6
boundary patches: 5
point zones: 0
face zones: 0
cell zones: 0

Overall number of cells of each type:
hexahedra: 4972600
prisms: 0
wedges: 0
pyramids: 0
tet wedges: 0
tetrahedra: 0
polyhedra: 0

Checking topology...
Boundary definition OK.
Cell to face addressing OK.
Point usage OK.
Upper triangular ordering OK.
Face vertices OK.
Number of regions: 1 (OK).

Checking patch topology for multiply connected surfaces...
Patch Faces Points Surface topology
INLET 10810 11088 ok (non-closed singly connected)
OUTLET 10810 11088 ok (non-closed singly connected)
CYLINDER 21620 22080 ok (non-closed singly connected)
FRONT_CYC 105800 106260 ok (non-closed singly connected)
BACK_CYC 105800 106260 ok (non-closed singly connected)

Checking geometry...
Overall domain bounding box (-20 -20 0) (20 20 3.1396)
Mesh has 3 geometric (non-empty/wedge) directions (1 1 1)
Mesh has 3 solution (non-empty) directions (1 1 1)
Boundary openness (-4.50402e-17 2.18592e-18 9.74372e-16) OK.
Max cell openness = 3.11523e-16 OK.
Max aspect ratio = 75.1683 OK.
Minimum face area = 6.12404e-06. Maximum face area = 0.145038. Face area magnitudes OK.
Min volume = 4.09086e-07. Max volume = 0.00968854. Total volume = 3942.75. Cell volumes OK.
Mesh non-orthogonality Max: 0.106719 average: 0.00668153
Non-orthogonality check OK.
Face pyramids OK.
Max skewness = 0.0172544 OK.
Coupled point location match (average 8.88178e-16) OK.

Mesh OK.

End



\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\ \\\\\\\\\\\\\\\\\\\\\\\\\\
k file
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 4.1 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
location "0.001";
object k;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions [0 2 -2 0 0 0 0];

internalField uniform 0;

boundaryField
{
INLET
{
type fixedValue;
value uniform 2e-05;
}
OUTLET
{
type inletOutlet;
inletValue uniform 0;
value uniform 0;
}
CYLINDER
{
type fixedValue;
value uniform 0;
}
FRONT_CYC
{
type cyclic;
}
BACK_CYC
{
type cyclic;
}
}


// ************************************************** *********************** //


U file

/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 4.1 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volVectorField;
location "0.001";
object U;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions [0 1 -1 0 0 0 0];

internalField uniform (0 0 0);

boundaryField
{
INLET
{
type turbulentInlet;
fluctuationScale (0.02 0.01 0.01);
referenceField uniform (1.5 0 0);
alpha 0.1;
value uniform (1.5 0 0);
}
OUTLET
{
type zeroGradient;
}
CYLINDER
{
type noSlip;
}
FRONT_CYC
{
type cyclic;
}
BACK_CYC
{
type cyclic;
}
}


// ************************************************** *********************** //



nuTilda file
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 4.1 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
location "0.001";
object nuTilda;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions [0 2 -1 0 0 0 0];

internalField uniform 0;

boundaryField
{
INLET
{
type fixedValue;
value uniform 0;
}
OUTLET
{
type inletOutlet;
inletValue uniform 0;
value uniform 0;
}
CYLINDER
{
type fixedValue;
value uniform 0;
}
FRONT_CYC
{
type cyclic;
}
BACK_CYC
{
type cyclic;
}
}


// ************************************************** *********************** //


nut file:

/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 4.1 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
location "0.001";
object nut;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions [0 2 -1 0 0 0 0];

internalField uniform 0;

boundaryField
{
INLET
{
type zeroGradient;
}
OUTLET
{
type zeroGradient;
}
CYLINDER
{
type zeroGradient;
}
FRONT_CYC
{
type cyclic;
}
BACK_CYC
{
type cyclic;
}
}


// ************************************************** *********************** //

p file
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 4.1 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
location "0.001";
object p;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions [0 2 -2 0 0 0 0];

internalField uniform 0;

boundaryField
{
INLET
{
type zeroGradient;
}
OUTLET
{
type fixedValue;
value uniform 0;
}
CYLINDER
{
type zeroGradient;
}
FRONT_CYC
{
type cyclic;
}
BACK_CYC
{
type cyclic;
}
}


// ************************************************** *********************** //

s file

/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 4.1 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
location "0.001";
object s;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions [0 0 0 0 0 0 0];

internalField uniform 0;

boundaryField
{
INLET
{
type fixedValue;
value uniform 1;
}
OUTLET
{
type inletOutlet;
inletValue uniform 0;
value uniform 0;
}
CYLINDER
{
type zeroGradient;
}
FRONT_CYC
{
type cyclic;
}
BACK_CYC
{
type cyclic;
}
}


// ************************************************** *********************** //

turbulence file:

/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 4.x |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
location "constant";
object turbulenceProperties;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

simulationType LES;

LES
{
LESModel kEqn;

turbulence on;

printCoeffs on;

delta cubeRootVol;

cubeRootVolCoeffs
{
deltaCoeff 1;
}

PrandtlCoeffs
{
delta cubeRootVol;
cubeRootVolCoeffs
{
deltaCoeff 1;
}

smoothCoeffs
{
delta cubeRootVol;
cubeRootVolCoeffs
{
deltaCoeff 1;
}

maxDeltaRatio 1.1;
}

Cdelta 0.158;
}

vanDriestCoeffs
{
delta cubeRootVol;
cubeRootVolCoeffs
{
deltaCoeff 1;
}

smoothCoeffs
{
delta cubeRootVol;
cubeRootVolCoeffs
{
deltaCoeff 1;
}

maxDeltaRatio 1.1;
}

Aplus 26;
Cdelta 0.158;
}

smoothCoeffs
{
delta cubeRootVol;
cubeRootVolCoeffs
{
deltaCoeff 1;
}

maxDeltaRatio 1.1;
}
}


// ************************************************** *********************** //