CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Running, Solving & CFD

trouble running simpleFoam with imported polyMesh

Register Blogs Community New Posts Updated Threads Search

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   November 5, 2022, 22:23
Default trouble running simpleFoam with imported polyMesh
  #1
New Member
 
N/A
Join Date: Nov 2021
Posts: 3
Rep Power: 4
noble6ix is on a distinguished road
using the pitzDaily tutorial, tried importing my own mesh and running simpleFoam to test it but received an error after computing Uz:
Code:
/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2206                                  |
|   \\  /    A nd           | Website:  www.openfoam.com                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
Build  : _76d719d1e6-20220624 OPENFOAM=2206 version=v2206
Arch   : "LSB;label=32;scalar=64"
Exec   : simpleFoam
Date   : Oct 28 2022
Time   : 14:07:40
Host   : tinkercliffs1
PID    : 149462
I/O    : uncollated
Case   : /home/jag21791/OpenFOAM-v2206/run/pitzDaily
nProcs : 1
trapFpe: Floating point exception trapping enabled (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster (fileModificationSkew 5, maxFileModificationPolls 20)
allowSystemOperations : Allowing user-supplied system call operations

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

Create mesh for time = 0


SIMPLE: convergence criteria
    field p	 tolerance 0.01
    field U	 tolerance 0.001
    field "(k|epsilon|omega|f|v2)"	 tolerance 0.001

Reading field p

Reading field U

Reading/calculating face flux field phi

Selecting incompressible transport model Newtonian
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;
}

No MRF models present

No finite volume options present

Starting time loop

streamLine streamlines:
    Employing velocity field U
    automatic track length specified through number of sub cycles : 5

Time = 1

smoothSolver:  Solving for Ux, Initial residual = 1, Final residual = 4.16282e+07, No Iterations 1000
smoothSolver:  Solving for Uy, Initial residual = 1, Final residual = 9.73246e+06, No Iterations 1000
smoothSolver:  Solving for Uz, Initial residual = 0.999858, Final residual = 0.0456548, No Iterations 1
#0  Foam::error::printStack(Foam::Ostream&) at ??:?
#1  Foam::sigFpe::sigHandler(int) at ??:?
#2  ? in /lib64/libpthread.so.0
#3  Foam::GAMGSolver::scale(Foam::Field<double>&, Foam::Field<double>&, Foam::lduMatrix const&, Foam::FieldField<Foam::Field, double> const&, Foam::UPtrList<Foam::lduInterfaceField const> const&, Foam::Field<double> const&, unsigned char) const at ??:?
#4  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 ??:?
#5  Foam::GAMGSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
#6  Foam::fvMatrix<double>::solveSegregated(Foam::dictionary const&) at ??:?
#7  Foam::fvMatrix<double>::solveSegregatedOrCoupled(Foam::dictionary const&) at ??:?
#8  Foam::fvMesh::solve(Foam::fvMatrix<double>&, Foam::dictionary const&) const at ??:?
#9  ? at ??:?
#10  __libc_start_main in /lib64/libc.so.6
#11  ? at ??:?
Floating point exception (core dumped)
I ran checkMesh before this and received an Ok. Here is the output of checkMesh:
Code:
/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2206                                  |
|   \\  /    A nd           | Website:  www.openfoam.com                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
Build  : _76d719d1e6-20220624 OPENFOAM=2206 version=v2206
Arch   : "LSB;label=32;scalar=64"
Exec   : checkMesh
Date   : Oct 28 2022
Time   : 14:04:44
Host   : tinkercliffs1
PID    : 90255
I/O    : uncollated
Case   : /home/jag21791/OpenFOAM-v2206/run/pitzDaily
nProcs : 1
trapFpe: Floating point exception trapping enabled (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster (fileModificationSkew 5, maxFileModificationPolls 20)
allowSystemOperations : Allowing user-supplied system call operations

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

Create mesh for time = 0

Time = 0

Mesh stats 
    points:           30762
    internal points:  0
    faces:            60680
    internal faces:   29920
    cells:            15100
    faces per cell:   6
    boundary patches: 6
    point zones:      0
    face zones:       0
    cell zones:       0

Overall number of cells of each type:
    hexahedra:     15100
    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                  
    centerline          120      244      ok (non-closed singly connected)  
    farField            170      342      ok (non-closed singly connected)  
    inlet               70       142      ok (non-closed singly connected)  
    outlet              110      222      ok (non-closed singly connected)  
    sideWedge           90       182      ok (non-closed singly connected)  
    symmetry            30200    30762    ok (non-closed singly connected)  

Checking faceZone topology for multiply connected surfaces...
    No faceZones found.

Checking basic cellZone addressing...
    No cellZones found.

Checking geometry...
    Overall domain bounding box (-2 0 0) (5 8.15296 0.05)
    Mesh has 3 geometric (non-empty/wedge) directions (1 1 1)
    Mesh has 3 solution (non-empty) directions (1 1 1)
    Boundary openness (-3.84497e-18 -2.46502e-18 4.94153e-16) OK.
    Max cell openness = 2.90805e-16 OK.
    Max aspect ratio = 50.4182 OK.
    Minimum face area = 2.04966e-05. Maximum face area = 0.0455154.  Face area magnitudes OK.
    Min volume = 1.02483e-06. Max volume = 0.00164667.  Total volume = 2.82968.  Cell volumes OK.
    Mesh non-orthogonality Max: 13.1711 average: 0.817124
    Non-orthogonality check OK.
    Face pyramids OK.
    Max skewness = 0.309378 OK.
    Coupled point location match (average 0) OK.

Mesh OK.

End
What could be causing this error? Thanks in advance
Attached Files
File Type: zip initialConditions.zip (10.2 KB, 0 views)
File Type: zip boundaryFile.zip (948 Bytes, 0 views)

Last edited by noble6ix; November 9, 2022 at 00:24. Reason: *attaching boundary file and 0/ files
noble6ix is offline   Reply With Quote

Old   November 9, 2022, 00:06
Default
  #2
New Member
 
N/A
Join Date: Nov 2021
Posts: 3
Rep Power: 4
noble6ix is on a distinguished road
Update: since the residuals were so high on the velocity, I modified my initial conditions to change the wall functions to non-wall functions, mainly of zeroGradient type. Greatly reduced residuals but produces the same error
noble6ix is offline   Reply With Quote

Reply


Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Floating point exception while running simpleFoam Shif OpenFOAM Running, Solving & CFD 0 March 24, 2022 14:40
using METIS functions in fortran dokeun Main CFD Forum 7 January 29, 2013 04:06
[Gmsh] Import problem ARC OpenFOAM Meshing & Mesh Conversion 0 February 27, 2010 10:56
Statically Compiling OpenFOAM Issues herzfeldd OpenFOAM Installation 21 January 6, 2009 09:38
Error running simpleFoam in parallel skabilan OpenFOAM Running, Solving & CFD 2 August 29, 2008 09:42


All times are GMT -4. The time now is 19:59.