|
[Sponsors] |
trouble running simpleFoam with imported polyMesh |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
November 5, 2022, 23:23 |
trouble running simpleFoam with imported polyMesh
|
#1 |
New Member
N/A
Join Date: Nov 2021
Posts: 3
Rep Power: 5 |
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 |
|
November 9, 2022, 01:06 |
|
#2 |
New Member
N/A
Join Date: Nov 2021
Posts: 3
Rep Power: 5 |
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
|
|
|
|
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 |