CFD Online Discussion Forums

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   OpenFOAM Running, Solving & CFD (https://www.cfd-online.com/Forums/openfoam-solving/)
-   -   Symmetry boundary condition in a heating pipe (https://www.cfd-online.com/Forums/openfoam-solving/182052-symmetry-boundary-condition-heating-pipe.html)

albertokarkour December 29, 2016 11:13

Symmetry boundary condition in a heating pipe
 
Hello! Me and a friend are solving an example from Victor's Pozzobon tutorials, example 8 to be specific and we are having troubles setting the 'symmetry' boundary condition, when we set it in the Temperature file like this:

Code:

FoamFile
{
    version    2.0;
    format      ascii;
    class      volScalarField;
    object      p;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [0 0 0 1 0 0 0];

internalField  uniform 20;

boundaryField
{
    inlet
    {
        type            fixedValue;
        value          uniform 20;
    }

    outlet
    {
        type            zeroGradient;
    }

    wall
    {
        type            convectiveHeatFlux;
    }

    axis
    {
        type            symmetry;
    }
   
    back
    {
        type            wedge;
    }

    front
    {
        type            wedge;
    }

}

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

We keep getting this message after running the solver:

Code:

/*---------------------------------------------------------------------------*\
| =========                |                                                |
| \\      /  F ield        | OpenFOAM: The Open Source CFD Toolbox          |
|  \\    /  O peration    | Version:  3.0.x                                |
|  \\  /    A nd          | Web:      www.OpenFOAM.org                      |
|    \\/    M anipulation  |                                                |
|*---------------------------------------------------------------------------*|
|* OpenFOAM for Windows 16.06 (v1)                                          *|
|* Built by CFD Support, www.cfdsupport.com (based on Symscape).            *|
\*---------------------------------------------------------------------------*/
Build  : 3.0.x-ac3f6c67e02f
Exec  : C:\OpenFOAM\ALBERTO-3.0.x\platforms\cygwin64mingw-w64DPInt32Opt\bin\simpleThermFoam.exe
Date  : Dec 29 2016
Time  : 15:52:54
Host  : "ALBERTOKARKOUR"
PID    : 1660
Case  : C:/OpenFOAM/ALBERTO-3.0.x/run/simpleThermFoam
nProcs : 1
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Allowing user-supplied system call operations

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

Create mesh for time = 0


SIMPLE: convergence criteria
    field p      tolerance 1e-006
    field U      tolerance 0.01
    field T      tolerance 1e-006

Reading field p

--> FOAM Warning :
    From function polyBoundaryMesh::groupPatchIDs() const
    in file meshes/polyMesh/polyBoundaryMesh/polyBoundaryMesh.C at line 461
    Patch wall specifies a group wall which is also a patch name. This might give problems later on.
Reading field U

Reading/calculating face flux field phi

Selecting incompressible transport model Newtonian
Selecting turbulence model type laminar
Reading field T



--> FOAM FATAL IO ERROR:

    patch type 'patch' not constraint type 'symmetry'
    for patch axis of field T in file "C:/OpenFOAM/ALBERTO-3.0.x/run/simpleThermFoam/0/T"

file: C:/OpenFOAM/ALBERTO-3.0.x/run/simpleThermFoam/0/T.boundaryField.axis from line 41 to line 41.

    From function symmetryFvPatchField<Type>::symmetryFvPatchField
(
    const fvPatch& p,
    const Field<Type>& field,
    const dictionary& dict
)

    in file fields/fvPatchFields/constraint/symmetry/symmetryFvPatchField.C at line 99.

FOAM exiting

What can we do to solve the error?

Many Thanks!

arvindpj December 29, 2016 12:51

Based on your 0/T file, I am seeing that you are trying to run a 2D axisymmetric simulation.

Place the createPatchDict file in "system" directory and run "createPatch".
This will remove the patches with no faces in them namely the pathch: axis


Below are the contents of createPatchDict

Code:


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

// Do a synchronisation of coupled points after creation of any patches.
// Note: this does not work with points that are on multiple coupled patches
//      with transformations (i.e. cyclics).
pointSync false;

// Patches to create. An empty patch list just removes patches with zero
// faces from $FOAM_CASE/constant/polyMesh/boundary.
patches
(
);

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

Cheers
:)


All times are GMT -4. The time now is 07:29.