[Xaxte]

Hello,



We tried to generate a mesh around an imported geometry (stl.file) with snappyHexMesh. It says that it generated a mesh that was ok, and it looks ok in paraview. But when we run the rhoPimpleFoam solver we get the following error message:

/*---------------------------------------------------------------------------*\
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     | Website:  https://openfoam.org
    \\  /    A nd           | Version:  7
     \\/     M anipulation  |
\*---------------------------------------------------------------------------*/
Build  : 7-6ef561967074
Exec   : rhoPimpleFoam
Date   : Nov 16 2019
Time   : 18:33:36
Host   : "kjersti-VirtualBox"
PID    : 3426
I/O    : uncollated
Case   : /home/kjersti/OpenFOAM/kjersti-7/run/prism
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster (fileModificationSkew 10)
allowSystemOperations : Allowing user-supplied system call operations

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

Create mesh for time = 0


PIMPLE: No convergence criteria found


PIMPLE: No corrector convergence criteria found
        Calculations will do 2 corrections


PIMPLE: Operating solver in transient mode with 2 outer correctors


Reading thermophysical properties

Selecting thermodynamics package 
{
    type            hePsiThermo;
    mixture         pureMixture;
    transport       const;
    thermo          hConst;
    equationOfState perfectGas;
    specie          specie;
    energy          sensibleInternalEnergy;
}

Reading field U

Reading/calculating face flux field phi

Creating turbulence model

Selecting turbulence model type RAS
Selecting RAS turbulence model kEpsilon
RAS
{
    RASModel        kEpsilon;
    turbulence      on;
    printCoeffs     on;
    Cmu             0.09;
    C1              1.44;
    C2              1.92;
    C3              0;
    sigmak          1;
    sigmaEps        1.3;
}

Creating field dpdt

Creating field kinetic energy K

No MRF models present

No finite volume options present
Courant Number mean: 0.000372342 max: 0.0298496

Starting time loop

Courant Number mean: 0.000372342 max: 0.0298496
Time = 5e-09

PIMPLE: Iteration 1
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
smoothSolver:  Solving for Ux, Initial residual = 1, Final residual = 1.87371e-07, No Iterations 2
smoothSolver:  Solving for Uy, Initial residual = 1, Final residual = 7.07998e-08, No Iterations 2
smoothSolver:  Solving for Uz, Initial residual = 1, Final residual = 4.96554e-07, No Iterations 1
smoothSolver:  Solving for e, Initial residual = 1, Final residual = 2.6378e-07, No Iterations 2
smoothSolver:  Solving for p, Initial residual = 1, Final residual = 7.11171e-07, No Iterations 1
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 8.63045e-12, global = -2.41024e-12, cumulative = -2.41024e-12
PIMPLE: Iteration 2
smoothSolver:  Solving for Ux, Initial residual = 2.91336e-05, Final residual = 1.51895e-08, No Iterations 1
smoothSolver:  Solving for Uy, Initial residual = 0.00102276, Final residual = 6.36141e-07, No Iterations 1
smoothSolver:  Solving for Uz, Initial residual = 1.53354e-05, Final residual = 5.23964e-09, No Iterations 1
smoothSolver:  Solving for e, Initial residual = 0.00015063, Final residual = 1.26379e-07, No Iterations 1
smoothSolver:  Solving for p, Initial residual = 2.38202e-05, Final residual = 1.62645e-10, No Iterations 1
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 1.61197e-13, global = 1.98876e-14, cumulative = -2.39035e-12


--> FOAM FATAL ERROR: 
Cannot find face 165412 in any of the patches 1(rocket_scaled)
It seems your patches are not consistent with the mesh : internalFaces:158590  total number of faces:191060

    From function Foam::label Foam::polyBoundaryMesh::whichPatch(Foam::label) const
    in file meshes/polyMesh/polyBoundaryMesh/polyBoundaryMesh.C at line 747.

FOAM aborting

#0  Foam::error::printStack(Foam::Ostream&) at ??:?
#1  Foam::error::abort() at ??:?
#2  Foam::polyBoundaryMesh::whichPatch(int) const at ??:?
#3  void Foam::fvMatrix<double>::setValuesFromList<Foam::UList>(Foam::UList<int> const&, Foam::UList<double> const&) at ??:?
#4  Foam::epsilonWallFunctionFvPatchScalarField::manipulateMatrix(Foam::fvMatrix<double>&) at ??:?
#5  Foam::fvMatrix<double>::boundaryManipulate(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>::Boundary&) at ??:?
#6  Foam::RASModels::kEpsilon<Foam::EddyDiffusivity<Foam::ThermalDiffusivity<Foam::CompressibleTurbulenceModel<Foam::fluidThermo> > > >::correct() at ??:?
#7  ? in "/opt/openfoam7/platforms/linux64GccDPInt32Opt/bin/rhoPimpleFoam"
#8  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#9  ? in "/opt/openfoam7/platforms/linux64GccDPInt32Opt/bin/rhoPimpleFoam"
Aborted (core dumped)

The constant/polyMesh/boundary file is this

/*--------------------------------*- C++ -*----------------------------------*\
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     | Website:  https://openfoam.org
    \\  /    A nd           | Version:  7
     \\/     M anipulation  |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       polyBoundaryMesh;
    location    "1.5e-08/polyMesh";
    object      boundary;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

1
(
    
    rocket_scaled
    {
        type            wall;
        inGroups        1(wall);
        nFaces          20066;
        startFace       170994;
    }
)

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

Have any of you experienced something similar or know what could be wrong? Attached is a picture of the mesh. If you need any more files, just tell me!


Thank you so much for your help, it would be highly appreciated!