Friday, November 30, 2018

The earth rings

https://www.nationalgeographic.com/science/2018/11/strange-earthquake-waves-rippled-around-world-earth-geology/

More Scale the Mesh 3


After: transformPoints -scale '(0.001 0.001 0.001)'    (Scale by mm)

/*--------------------------------*- 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       vectorField;
    location    "constant/polyMesh";
    object      points;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //


882
(
(0 0 0)
(0.02 0 0)
(0 0 0.2)
(0.02 0 0.2)
(0 0.2 0)
(0.02 0.2 0)
(0 0.2 0.2)
(0.02 0.2 0.2)
(0 0 0.01)
(0 0 0.02)
(0 0 0.03)
(0 0 0.04)
(0 0 0.05)
(0 0 0.06)
(0 0 0.07)
(0 0 0.08)
(0 0 0.09)
(0 0 0.1)
(0 0 0.11)
(0 0 0.12)
(0 0 0.13)
(0 0 0.14)
(0 0 0.15)
(0 0 0.16)
(0 0 0.17)
(0 0 0.18)
(0 0 0.19)


-----------------------

Before:


/*--------------------------------*- 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       vectorField;
    location    "constant/polyMesh";
    object      points;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //


882
(
(0 0 0)
(20 0 0)
(0 0 200)
(20 0 200)
(0 200 0)
(20 200 0)
(0 200 200)
(20 200 200)
(0 0 10)
(0 0 20)
(0 0 30)
(0 0 40)
(0 0 50)
(0 0 60)
(0 0 70)
(0 0 80)
(0 0 90)
(0 0 100)
(0 0 110)
(0 0 120)
(0 0 130)
(0 0 140)
(0 0 150)
(0 0 160)
(0 0 170)
(0 0 180)
(0 0 190)

-------------------------

transformPoints -scale '(0.01 0.01 0.01)'      scale by cm

Thursday, November 29, 2018

More Scale the Mesh 2

/*--------------------------------*- 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       dictionary;
    object      blockMeshDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

convertToMeters 0.1;

vertices
(
    (0 0 0)
    (1 0 0)
    (1 1 0)
    (0 1 0)
    (0 0 0.1)
    (1 0 0.1)
    (1 1 0.1)
    (0 1 0.1)
);
-----------------------------------------------------------------------------------

/*--------------------------------*- 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       vectorField;
    location    "constant/polyMesh";
    object      points;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //


882
(
(0 0 0)
(0.005 0 0)
(0.01 0 0)
(0.015 0 0)
(0.02 0 0)
(0.025 0 0)
(0.03 0 0)
(0.035 0 0)
(0.04 0 0)
(0.045 0 0)
(0.05 0 0)
(0.055 0 0)
(0.06 0 0)
(0.065 0 0)
(0.07 0 0)
(0.075 0 0)
(0.08 0 0)
(0.085 0 0)
(0.09 0 0)
(0.095 0 0)
(0.1 0 0)
----------------------------------------------------------------------------------------------------

The step size is 0.1 / 20 or the "convertToMeters 0.1" divided by the number of patches "20".

What number of patches would one use for an irregular shape? Would it be based on the number of patches in the shortest face?

----------------------------------------------------------------------------------------------------

/*--------------------------------*- 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       polyBoundaryMesh;
    location    "constant/polyMesh";
    object      boundary;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

3
(
    movingWall
    {
        type            wall;
        inGroups        1(wall);
        nFaces          20;
        startFace       760;
    }
    fixedWalls
    {
        type            wall;
        inGroups        1(wall);
        nFaces          60;
        startFace       780;
    }
    frontAndBack
    {
        type            empty;
        inGroups        1(empty);
        nFaces          800;
        startFace       840;
    }
)

More Scale the MEsh

incompressible/icoFoam/Mesh_1
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

Processing tag:164
Starting reading units at line 3.
l:1
units:"  SI: Meter (newton)"
unitType:2
Unit factors:
    Length scale       : 1
    Force scale        : 1
    Temperature scale  : 1
    Temperature offset : 273.15


Processing tag:2420
Skipping tag 2420 on line 9
Skipping section at line 9.

Processing tag:2411
Starting reading points at line 20.
Read 882 points.

Processing tag:2412
Starting reading cells at line 1787.
First occurrence of element type 11 for cell 1 at line 1788
First occurrence of element type 44 for cell 85 at line 2280
First occurrence of element type 115 for cell 1045 at line 4040
Read 400 cells and 880 boundary faces.

Processing tag:2467
Starting reading patches at line 4842.
For group 1 named sourceFaceleft trying to read 400 patch face indices.
For group 2 named destinationFaceleft trying to read 400 patch face indices.
For group 3 named channelWallsleft trying to read 40 patch face indices.
For group 4 named leftRightWallleft trying to read 20 patch face indices.
For group 5 named inlet trying to read 20 patch face indices.

Sorting boundary faces according to group (patch)
0: sourceFaceleft is patch
1: destinationFaceleft is patch
2: channelWallsleft is patch
3: leftRightWallleft is patch
4: inlet is patch

Constructing mesh with non-default patches of size:
    sourceFaceleft    400
    destinationFaceleft    400
    channelWallsleft    40
    leftRightWallleft    20
    inlet    20

End
---------------------------------------------
/*--------------------------------*- 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       vectorField;
    location    "constant/polyMesh";
    object      points;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //


(0 0 0)
(0 0 10)
(0 0 20)
(0 0 30)
(0 0 40)
(0 0 50)
(0 0 60)
(0 0 70)
(0 0 80)
(0 0 90)
(0 0 100)
(0 0 110)
(0 0 120)
(0 0 130)
(0 0 140)
(0 0 150)
(0 0 160)
(0 0 170)
(0 0 180)
(0 0 190)
(0 0 200)



Saturday, November 24, 2018

Scale the Mesh

Scale the mesh:




By looking at:

https://www.cfd-online.com/Forums/openfoam-meshing/198268-how-scale-stl-when-importing-salome-any-cad-tool-openfoam.html

try:

transformPoints -scale '(0.001 0.001 0.001)'

https://cfd.direct/openfoam/user-guide/v6-basic-file-format/

https://openfoamwiki.net/index.php/TransformPoints4


What is the default mesh size???

http://www.salome-platform.org/forum/forum_10/549857363

-------------------------------------------------------------------------------------------------------
  -1
   164
         1  SI: Meter (newton)         2
    1.0000000000000000E+0    1.0000000000000000E+0    1.0000000000000000E+0
    2.7314999999999998E+2
    -1
    -1
  2420
         1
SMESH_Mesh
         1         0         0
Global Cartesian Coordinate System
    1.0000000000000000E+0    0.0000000000000000E+0    0.0000000000000000E+0
    0.0000000000000000E+0    1.0000000000000000E+0    0.0000000000000000E+0
    0.0000000000000000E+0    0.0000000000000000E+0    1.0000000000000000E+0
    0.0000000000000000E+0    0.0000000000000000E+0    0.0000000000000000E+0
    -1
    -1
  2411
         1         1         1        11
   2.0000000000000000E+01   4.5500000000000000E+02   0.0000000000000000E+00
         2         1         1        11
   0.0000000000000000E+00   4.5500000000000000E+02   0.0000000000000000E+00
         3         1         1        11
   2.0000000000000000E+01   5.2000000000000000E+02   0.0000000000000000E+00
         4         1         1        11
   0.0000000000000000E+00   5.2000000000000000E+02   0.0000000000000000E+00
         5         1         1        11
   2.0000000000000000E+01   1.6600000000000000E+02   5.0000000000000000E+02

-----------------------------------------------------------------------------------------------------
/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  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   : ideasUnvToFoam InletMesh.unv
Date   : Nov 26 2018
Time   : 18:32:30
Host   : "caelinux"
PID    : 22307
Case   : /home/caelinux/OpenFOAM/caelinux-4.1/run/tutorials/incompressible/icoFoam/test-InletMesh
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

Processing tag:164
Starting reading units at line 3.
l:1
units:"  SI: Meter (newton)"
unitType:2
Unit factors:
    Length scale       : 1
    Force scale        : 1
    Temperature scale  : 1
    Temperature offset : 273.15


-------------------------------------------------------------------------------
/*--------------------------------*- 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       vectorField;
    location    "constant/polyMesh";
    object      points;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //


882
(
(0 0 0)
(0.005 0 0)
(0.01 0 0)
(0.015 0 0)
(0.02 0 0)
(0.025 0 0)
(0.03 0 0)
(0.035 0 0)
(0.04 0 0)
(0.045 0 0)
(0.05 0 0)
(0.055 0 0)
(0.06 0 0)
(0.065 0 0)
(0.07 0 0)
(0.075 0 0)
(0.08 0 0)
(0.085 0 0)
(0.09 0 0)
(0.095 0 0)
(0.1 0 0)

---------------------------------

See: https://cfd.direct/openfoam/user-guide/v6-cavity/   ... length is 0.1 meters

To do: Nov 24th

Merge:

https://github.com/bshambaugh/openfoam-experiments/tree/master/InletMesh

with
https://github.com/bshambaugh/openfoam-experiments/tree/master/NGpitzDaily

This gives:

https://github.com/bshambaugh/openfoam-experiments/tree/master/BNGpitzDaily





with the error:

smoothSolver:  Solving for Ux, Initial residual = 8.87405e-07, Final residual = 8.87405e-07, No Iterations 0
smoothSolver:  Solving for Uy, Initial residual = 0.168394, Final residual = 4.90847e-07, No Iterations 11
smoothSolver:  Solving for Uz, Initial residual = 0.762311, Final residual = 8.70185e-07, No Iterations 13
#0  Foam::error::printStack(Foam::Ostream&) at ??:?
#1  Foam::sigFpe::sigHandler(int) at ??:?
#2  ? in "/lib/x86_64-linux-gnu/libc.so.6"
#3  double Foam::sumProd<double>(Foam::UList<double> const&, Foam::UList<double> const&) at ??:?
#4  Foam::PCG::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
#5  Foam::GAMGSolver::solveCoarsestLevel(Foam::Field<double>&, Foam::Field<double> const&) const at ??:?
#6  Foam::GAMGSolver::Vcycle(Foam::PtrList<Foam::lduMatrix::smoother> const&, Foam::Field<double>&, Foam::Field<double> const&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::PtrList<Foam::Field<double> >&, Foam::PtrList<Foam::Field<double> >&, unsigned char) const at ??:?
#7  Foam::GAMGSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
#8  Foam::fvMatrix<double>::solveSegregated(Foam::dictionary const&) at ??:?
#9  Foam::fvMatrix<double>::solve(Foam::dictionary const&) at ??:?
#10  Foam::fvMatrix<double>::solve() at ??:?
#11  ? at ??:?
#12  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#13  ? at ??:?
Floating point exception (core dumped)


Wednesday, November 21, 2018

Tuesday, November 20, 2018

freeStreamPressure

https://www.simscale.com/docs/content/simulation/model/boundaryConditionTypes/OF_pressure.html

Monday, November 19, 2018

Method of creation

1) Create folders:

InletMesh

FineMesh

CourseMesh


2) mergeMeshs -overwrite  FineMesh InletMesh
    mergeMeshs -overwrite  CourseMesh FineMesh


edit: mergeMesh to mergeMeshs on Dec. 1st, 2018

Error with RAS Model

 For squareCylinder_kEpsilon/   :

Error with RAS Model.

https://stackoverflow.com/questions/48421407/openfoam-foam-fatal-error-unknown-turbulencemodel-type-rasmodel

the correction gives a new error.

https://www.cfd-online.com/Forums/openfoam-solving/165495-error-div-nueff-dev2-t-grad-u.html

--------

For CFDwOpenFoam/LaminarVortexShedding :

https://bugs.openfoam.org/view.php?id=2347


0002347: Error while calculating forces - could not find rho

Saturday, November 17, 2018

stl or cad tool --- scaling | Turbulence Modelling

https://www.cfd-online.com/Forums/openfoam-meshing/198268-how-scale-stl-when-importing-salome-any-cad-tool-openfoam.html


Cad tool --- scaling


---------------------

https://www.youtube.com/playlist?list=PLOMAD1NLiFnw4sfxSkiCZwj20CzHjOuGh

Turbulence Modelling


https://github.com/meinarsve/CFDwOpenFoam


Also see Martin Einarsve on YouTube.

Tuesday, November 13, 2018

ideasToFoam length scale

ideasUnvToFoam LeftBottom.unv
/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  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   : ideasUnvToFoam LeftBottom.unv
Date   : Nov 14 2018
Time   : 06:07:50
Host   : "caelinux"
PID    : 16488
Case   : /home/caelinux/OpenFOAM/caelinux-4.1/run/tutorials/incompressible/icoFoam/LeftBottom
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

Processing tag:164
Starting reading units at line 3.
l:1
units:"  SI: Meter (newton)"
unitType:2
Unit factors:
    Length scale       : 1
    Force scale        : 1
    Temperature scale  : 1
    Temperature offset : 273.15


Processing tag:2420
Skipping tag 2420 on line 9
Skipping section at line 9.

Processing tag:2411
Starting reading points at line 20.
Read 882 points.

Processing tag:2412
Starting reading cells at line 1787.
First occurrence of element type 11 for cell 1 at line 1788
First occurrence of element type 44 for cell 165 at line 2280
First occurrence of element type 115 for cell 1045 at line 4040
Read 400 cells and 880 boundary faces.

Processing tag:2467
Starting reading patches at line 4842.
For group 1 named symmetryWallLeftBottom trying to read 20 patch face indices.
For group 2 named freeBoundryLeftBottom trying to read 20 patch face indices.
For group 3 named leftSide trying to read 20 patch face indices.
For group 4 named bottomToTop trying to read 20 patch face indices.
For group 5 named backFaceLeftBottom trying to read 400 patch face indices.
For group 6 named frontFaceLeftBottom trying to read 400 patch face indices.

Sorting boundary faces according to group (patch)
0: symmetryWallLeftBottom is patch
1: freeBoundryLeftBottom is patch
2: leftSide is patch
3: bottomToTop is patch
4: backFaceLeftBottom is patch
5: frontFaceLeftBottom is patch

Constructing mesh with non-default patches of size:
    symmetryWallLeftBottom    20
    freeBoundryLeftBottom    20
    leftSide    20
    bottomToTop    20
    backFaceLeftBottom    400
    frontFaceLeftBottom    400

End

error message for continuity

2elbow-LT-for-holly$ cat log.icoFoam
/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  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   : icoFoam
Date   : Nov 13 2018
Time   : 23:29:30
Host   : "caelinux"
PID    : 15899
Case   : /home/caelinux/OpenFOAM/caelinux-4.1/run/tutorials/incompressible/icoFoam/2elbow-LT-for-holly
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 mesh for time = 0


PISO: Operating solver in PISO mode

Reading transportProperties

Reading field p

AMI: Creating addressing and weights between 40 source faces and 5 target faces
AMI: Patch source sum(weights) min/max/average = 1, 1, 1
AMI: Patch target sum(weights) min/max/average = 1, 1, 1
Reading field U

Reading/calculating face flux field phi


Starting time loop

Time = 0.05

Courant Number mean: 0.00886767 max: 2.5
smoothSolver:  Solving for Ux, Initial residual = 1, Final residual = 1.26206e-11, No Iterations 1
smoothSolver:  Solving for Uy, Initial residual = 0, Final residual = 0, No Iterations 0
smoothSolver:  Solving for Uz, Initial residual = 1, Final residual = 1.26206e-11, No Iterations 1


--> FOAM FATAL ERROR:
Continuity error cannot be removed by adjusting the outflow.
Please check the velocity boundary conditions and/or run potentialFoam to initialise the outflow.
Total flux              : 2.63698e+08
Specified mass inflow   : 1.3e+07
Specified mass outflow  : 0
Adjustable mass outflow : 0


    From function bool Foam::adjustPhi(Foam::surfaceScalarField&, const volVectorField&, Foam::volScalarField&)
    in file cfdTools/general/adjustPhi/adjustPhi.C at line 107.

FOAM exiting

I need to tweak for the symmetryPlane boundary condition

/2elbow-LT-for-holly$ cat log.icoFoam
/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  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   : icoFoam
Date   : Nov 13 2018
Time   : 21:33:59
Host   : "caelinux"
PID    : 15652
Case   : /home/caelinux/OpenFOAM/caelinux-4.1/run/tutorials/incompressible/icoFoam/2elbow-LT-for-holly
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 mesh for time = 0


PISO: Operating solver in PISO mode

Reading transportProperties

Reading field p

AMI: Creating addressing and weights between 40 source faces and 5 target faces
AMI: Patch source sum(weights) min/max/average = 1, 1, 1
AMI: Patch target sum(weights) min/max/average = 1, 1, 1
Reading field U

Reading/calculating face flux field phi


Starting time loop

Time = 0.05

Courant Number mean: 0.00886767 max: 2.5
smoothSolver:  Solving for Ux, Initial residual = 1, Final residual = 1.26206e-11, No Iterations 1
smoothSolver:  Solving for Uy, Initial residual = 0, Final residual = 0, No Iterations 0
smoothSolver:  Solving for Uz, Initial residual = 1, Final residual = 1.26206e-11, No Iterations 1


--> FOAM FATAL ERROR:
Continuity error cannot be removed by adjusting the outflow.
Please check the velocity boundary conditions and/or run potentialFoam to initialise the outflow.
Total flux              : 2.63698e+08
Specified mass inflow   : 1.3e+07
Specified mass outflow  : 0
Adjustable mass outflow : 0


    From function bool Foam::adjustPhi(Foam::surfaceScalarField&, const volVectorField&, Foam::volScalarField&)
    in file cfdTools/general/adjustPhi/adjustPhi.C at line 107.

FOAM exiting

------------------------------------------------------------------
fvSolution:
----------------------------------------------------------------


/2elbow-LT-for-holly/system$ cat fvSolution
/*--------------------------------*- 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       dictionary;
    location    "system";
    object      fvSolution;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

solvers
{
    p
    {
        solver          PCG;
        preconditioner  DIC;
        tolerance       1e-06;
        relTol          0.05;
    }

    pFinal
    {
        $p;
        relTol          0;
    }

    U
    {
        solver          smoothSolver;
        smoother        symGaussSeidel;
        tolerance       1e-05;
        relTol          0;
    }
}

PISO
{
    nCorrectors     2;
    nNonOrthogonalCorrectors 2;
    pRefCell  0;
    pRefValue 0;
}


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


Sunday, November 11, 2018

search all tutorials for the cyclicAMI boundary condition to figure out how to use it...

/tutorials$ grep -rnw '.' -e 'cyclicAMI'
./multiphase/multiphaseInterDyMFoam/laminar/mixerVesselAMI2D/system/blockMeshDict.m4:791:        type            cyclicAMI;
./multiphase/multiphaseInterDyMFoam/laminar/mixerVesselAMI2D/system/blockMeshDict.m4:818:        type            cyclicAMI;
./multiphase/multiphaseInterDyMFoam/laminar/mixerVesselAMI2D/0.orig/p_rgh:45:        type            cyclicAMI;
./multiphase/multiphaseInterDyMFoam/laminar/mixerVesselAMI2D/0.orig/p_rgh:51:        type            cyclicAMI;
./multiphase/multiphaseInterDyMFoam/laminar/mixerVesselAMI2D/0.orig/alpha.water:45:        type            cyclicAMI;
./multiphase/multiphaseInterDyMFoam/laminar/mixerVesselAMI2D/0.orig/alpha.water:51:        type            cyclicAMI;
./multiphase/multiphaseInterDyMFoam/laminar/mixerVesselAMI2D/0.orig/alpha.oil:45:        type            cyclicAMI;
./multiphase/multiphaseInterDyMFoam/laminar/mixerVesselAMI2D/0.orig/alpha.oil:51:        type            cyclicAMI;
./multiphase/multiphaseInterDyMFoam/laminar/mixerVesselAMI2D/0.orig/alpha.mercury:45:        type            cyclicAMI;
./multiphase/multiphaseInterDyMFoam/laminar/mixerVesselAMI2D/0.orig/alpha.mercury:51:        type            cyclicAMI;
./multiphase/multiphaseInterDyMFoam/laminar/mixerVesselAMI2D/0.orig/alpha.air:45:        type            cyclicAMI;
./multiphase/multiphaseInterDyMFoam/laminar/mixerVesselAMI2D/0.orig/alpha.air:51:        type            cyclicAMI;
./multiphase/multiphaseInterDyMFoam/laminar/mixerVesselAMI2D/0.orig/U:47:        type            cyclicAMI;
./multiphase/multiphaseInterDyMFoam/laminar/mixerVesselAMI2D/0.orig/U:53:        type            cyclicAMI;
./multiphase/interPhaseChangeDyMFoam/propeller/system/createPatchDict:30:            type            cyclicAMI;
./multiphase/interPhaseChangeDyMFoam/propeller/system/createPatchDict:44:            type            cyclicAMI;
./multiphase/interDyMFoam/ras/mixerVesselAMI/system/createBafflesDict:38:                type            cyclicAMI;
./multiphase/interDyMFoam/ras/mixerVesselAMI/system/createBafflesDict:48:                type            cyclicAMI;
./multiphase/interDyMFoam/ras/mixerVesselAMI/0.orig/p_rgh:36:        type            cyclicAMI;
./multiphase/interDyMFoam/ras/mixerVesselAMI/0.orig/nut:44:        type            cyclicAMI;
./multiphase/interDyMFoam/ras/mixerVesselAMI/0.orig/k:44:        type            cyclicAMI;
./multiphase/interDyMFoam/ras/mixerVesselAMI/0.orig/epsilon:44:        type            cyclicAMI;
./multiphase/interDyMFoam/ras/mixerVesselAMI/0.orig/alpha.water:43:        type            cyclicAMI;
./multiphase/interDyMFoam/ras/mixerVesselAMI/0.orig/U:57:        type            cyclicAMI;
./mesh/foamyHexMesh/mixerVessel/system/createBafflesDict:80:                type            cyclicAMI;
./mesh/foamyHexMesh/mixerVessel/system/createBafflesDict:92:                type            cyclicAMI;
./mesh/foamyHexMesh/mixerVessel/0.orig/p_rgh:38:        type            cyclicAMI;
./mesh/foamyHexMesh/mixerVessel/0.orig/nut:44:        type            cyclicAMI;
./mesh/foamyHexMesh/mixerVessel/0.orig/k:44:        type            cyclicAMI;
./mesh/foamyHexMesh/mixerVessel/0.orig/epsilon:44:        type            cyclicAMI;
./mesh/foamyHexMesh/mixerVessel/0.orig/alpha.phase1:43:        type            cyclicAMI;
./mesh/foamyHexMesh/mixerVessel/0.orig/U:57:        type            cyclicAMI;
./incompressible/simpleFoam/pipeCyclic/system/blockMeshDict:80:        type cyclicAMI;
./incompressible/simpleFoam/pipeCyclic/system/blockMeshDict:94:        type cyclicAMI;
./incompressible/pimpleDyMFoam/propeller/system/createPatchDict:30:            type            cyclicAMI;
./incompressible/pimpleDyMFoam/propeller/system/createPatchDict:44:            type            cyclicAMI;
./incompressible/pimpleDyMFoam/oscillatingInletACMI2D/system/createBafflesDict:26:    // NOTE: cyclicAMI patches MUST BE defined PRIOR to their associted
./incompressible/pimpleDyMFoam/mixerVesselAMI2D/system/blockMeshDict.m4:791:        type            cyclicAMI;
./incompressible/pimpleDyMFoam/mixerVesselAMI2D/system/blockMeshDict.m4:818:        type            cyclicAMI;
./incompressible/pimpleDyMFoam/mixerVesselAMI2D/0/p:45:        type            cyclicAMI;
./incompressible/pimpleDyMFoam/mixerVesselAMI2D/0/p:51:        type            cyclicAMI;
./incompressible/pimpleDyMFoam/mixerVesselAMI2D/0/nut:48:        type            cyclicAMI;
./incompressible/pimpleDyMFoam/mixerVesselAMI2D/0/nut:54:        type            cyclicAMI;
./incompressible/pimpleDyMFoam/mixerVesselAMI2D/0/k:48:        type            cyclicAMI;
./incompressible/pimpleDyMFoam/mixerVesselAMI2D/0/k:54:        type            cyclicAMI;
./incompressible/pimpleDyMFoam/mixerVesselAMI2D/0/epsilon:48:        type            cyclicAMI;
./incompressible/pimpleDyMFoam/mixerVesselAMI2D/0/epsilon:54:        type            cyclicAMI;
./incompressible/pimpleDyMFoam/mixerVesselAMI2D/0/U:47:        type            cyclicAMI;
./incompressible/pimpleDyMFoam/mixerVesselAMI2D/0/U:53:        type            cyclicAMI;

Friday, November 9, 2018

example of more (for two faces: preliminary)

    object      boundary;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

10
(
    leftRightWallsright
    {
        type            cyclicAMI;
        neighbourPatch  leftRightWallleft;
        transform       noOrdering;
        nFaces          40;
        startFace       3880;
    }
    outlet
    {
        type            patch;
        nFaces          40;
        startFace       3920;
    }
    channelWallsright
    {
        type            patch;
        nFaces          80;
        startFace       3960;
    }
    destinationFaceright
    {
        type            patch;
        nFaces          1600;
        startFace       4040;
    }
    sourceFaceright
    {
        type            patch;
        nFaces          1600;
        startFace       5640;
    }
    sourceFaceleft
    {
        type            patch;
        nFaces          400;
        startFace       7240;
    }
    destinationFaceleft
    {
        type            patch;
        nFaces          400;
        startFace       7640;
    }
    channelWallsleft
    {
        type            patch;
        nFaces          40;
        startFace       8040;
    }
    leftRightWallleft
    {
        type            cyclicAMI;
        neighbourPatch  leftRightWallsright;
        transform       noOrdering;
        nFaces          20;
        startFace       8080;
    }
    inlet
    {
        type            patch;
        nFaces          20;
        startFace       8100;
    }
)

Thursday, November 8, 2018

start with one square, then work onto more

cat boundary
/*--------------------------------*- 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       polyBoundaryMesh;
    location    "constant/polyMesh";
    object      boundary;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

5
(
    sourceFaceleft
    {
        type            patch;
        nFaces          400;
        startFace       760;
    }
    destinationFaceleft
    {
        type            patch;
        nFaces          400;
        startFace       1160;
    }
    channelWallsleft
    {
        type            patch;
        nFaces          40;
        startFace       1560;
    }
    leftRightWallleft
    {
        type            patch;
        nFaces          20;
        startFace       1600;
    }
    inlet
    {
        type            patch;
        nFaces          20;
        startFace       1620;
    }
)

Tuesday, November 6, 2018

foam

cat log.icoFoam /*---------------------------------------------------------------------------*\ | ========= | | | \\ / 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 : icoFoam Date : Nov 07 2018 Time : 00:02:38 Host : "caelinux" PID : 16608 Case : /home/caelinux/OpenFOAM/caelinux-4.1/run/tutorials/incompressible/icoFoam/elbow-twoFace-right 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 mesh for time = 0 PISO: Operating solver in PISO mode Reading transportProperties Reading field p AMI: Creating addressing and weights between 40 source faces and 20 target faces AMI: Patch source sum(weights) min/max/average = 1, 1, 1 AMI: Patch target sum(weights) min/max/average = 1, 1, 1 Reading field U --> FOAM FATAL IO ERROR: "ill defined primitiveEntry starting at keyword 'value' on line 50 and ending at line 89" file: /home/caelinux/OpenFOAM/caelinux-4.1/run/tutorials/incompressible/icoFoam/elbow-twoFace-right/0/U at line 89. From function void Foam::primitiveEntry::readEntry(const Foam::dictionary&, Foam::Istream&) in file db/dictionary/primitiveEntry/primitiveEntryIO.C at line 189. FOAM exiting

possible solution

 caelinux@caelinux:~/OpenFOAM/caelinux-4.1/run/tutorials/incompressible/icoFoam/elbow-twoFace-right$ cp system/changeDictionaryDict 0/dictionaryReplacement
caelinux@caelinux:~/OpenFOAM/caelinux-4.1/run/tutorials/incompressible/icoFoam/elbow-twoFace-right$ cd 0
caelinux@caelinux:~/OpenFOAM/caelinux-4.1/run/tutorials/incompressible/icoFoam/elbow-twoFace-right/0$ ls
caelinux@caelinux:~/OpenFOAM/caelinux-4.1/run/tutorials/incompressible/icoFoam/elbow-twoFace-right$ changeDictionary

0$ cat dictionaryReplacement
/*--------------------------------*- 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       dictionary;
    location    "0";
    object      dictionaryReplacement;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dictionaryReplacement
{
    boundary
    {
        leftRightWallleft
        {
            type            cyclicAMI;
            neighbourPatch  leftRightWallsright;
            transform       noOrdering;
        }
        leftRightWallsright
        {
            type            cyclicAMI;
            neighbourPatch  leftRightWallleft;
            transform       noOrdering;
        }
    }
}

boundary
{
    leftRightWallleft
    {
        type            cyclicAMI;
        neighbourPatch  leftRightWallsright;
        transform       noOrdering;
    }
    leftRightWallsright
    {
        type            cyclicAMI;
        neighbourPatch  leftRightWallleft;
        transform       noOrdering;
    }
}

Monday, November 5, 2018

mergeMeshes experiment

/incompressible/icoFoam$ mergeMeshes elbow-twoFace-right/ elbow-twoFace-left/
/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  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   : mergeMeshes elbow-twoFace-right/ elbow-twoFace-left/
Date   : Nov 06 2018
Time   : 04:40:06
Host   : "caelinux"
PID    : 9319
Case   : /home/caelinux/OpenFOAM/caelinux-4.1/run/tutorials/incompressible/icoFoam
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

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Master:      "elbow-twoFace-right"  region region0
mesh to add: "elbow-twoFace-left"  region region0

Create Times
Reading master mesh for time = 0.1
Create mesh

Reading mesh to add for time = 0
Create mesh

Writing combined mesh to 0.15
patch names:
10
(
leftRightWallsright
outlet
channelWallsright
destinationFaceright
sourceFaceright
sourceFaceleft
destinationFaceleft
channelWallsleft
leftRightWallleft
inlet
)

patch dicts:
10
(

{
    type            patch;
    nFaces          40;
    startFace       4640;
}


{
    type            patch;
    nFaces          40;
    startFace       4680;
}


{
    type            patch;
    nFaces          80;
    startFace       4720;
}


{
    type            patch;
    nFaces          1600;
    startFace       4800;
}


{
    type            patch;
    nFaces          1600;
    startFace       6400;
}


{
    type            patch;
    nFaces          800;
    startFace       8000;
}


{
    type            patch;
    nFaces          800;
    startFace       8800;
}


{
    type            patch;
    nFaces          80;
    startFace       9600;
}


{
    type            patch;
    nFaces          40;
    startFace       9680;
}


{
    type            patch;
    nFaces          40;
    startFace       9720;
}

)

point zone names: 0()
face zone names: 0()
cell zone names: 0()
End