|
[Sponsors] | |||||
trouble running simpleFoam with imported polyMesh |
![]() |
|
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
|
|
|
#1 |
|
New Member
N/A
Join Date: Nov 2021
Posts: 3
Rep Power: 6 ![]() |
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)
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
Last edited by noble6ix; November 9, 2022 at 01:24. Reason: *attaching boundary file and 0/ files |
|
|
|
|
|
|
|
|
#2 |
|
New Member
N/A
Join Date: Nov 2021
Posts: 3
Rep Power: 6 ![]() |
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
|
|
|
|
|
|
![]() |
| Thread Tools | Search this Thread |
| Display Modes | |
|
|
Similar Threads
|
||||
| Thread | Thread Starter | Forum | Replies | Last Post |
| Floating point exception while running simpleFoam | Shif | OpenFOAM Running, Solving & CFD | 0 | March 24, 2022 15:40 |
| using METIS functions in fortran | dokeun | Main CFD Forum | 7 | January 29, 2013 05:06 |
| [Gmsh] Import problem | ARC | OpenFOAM Meshing & Mesh Conversion | 0 | February 27, 2010 11:56 |
| Statically Compiling OpenFOAM Issues | herzfeldd | OpenFOAM Installation | 21 | January 6, 2009 10:38 |
| Error running simpleFoam in parallel | skabilan | OpenFOAM Running, Solving & CFD | 2 | August 29, 2008 10:42 |