CFD Online Logo CFD Online URL
Home > Forums > Software User Forums > OpenFOAM

Problems with turbulent channel flow with pimpleFoam

Register Blogs Members List Search Today's Posts Mark Forums Read

LinkBack Thread Tools Search this Thread Display Modes
Old   April 8, 2022, 12:21
Default Problems with turbulent channel flow with pimpleFoam
New Member
Join Date: Mar 2022
Location: Hamburg
Posts: 2
Rep Power: 0
HeleShawman is on a distinguished road
Hi everyone!

I am trying to solve a 3D turbulent flow in a channel similar to this tutorial:

With the only difference that I am adding a small block downstream where I add an inlet that adds flow from the bottom upwards (to mix with the main stream)-see attached figure (mesh.png).

I use exactly the same controlDict, fvScheme and fvSolution files as in the abovementioned case (I attach the blockMeshDict file, which based in that case analogous file).

It seems a really simple problem very similar to one already solved in the OpenFOAM documentation, however, when I run the case using pimpleFoam: the following error appears:

PIMPLE: iteration 1
Turbulent DFSEM patch: inlet seeded 601 eddies with total volume 1.57943
smoothSolver: Solving for Ux, Initial residual = 1, Final residual = 0.0641366, No Iterations 1
smoothSolver: Solving for Uy, Initial residual = 1, Final residual = 0.0789141, No Iterations 1
smoothSolver: Solving for Uz, Initial residual = 1, Final residual = 0.0165905, No Iterations 1
GAMG: Solving for p, Initial residual = 1, Final residual = 8.55324e-07, No Iterations 99
time step continuity errors : sum local = 3.02402e-09, global = -8.48373e-10, cumulative = -8.48373e-10
PIMPLE: iteration 2
#0 Foam::error:rintStack(Foam::Ostream&) at ??:?
#1 Foam::sigFpe::sigHandler(int) at ??:?
#2 ? in /lib/x86_64-linux-gnu/
#3 Foam::smoothSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
#4 Foam::fvMatrix<Foam::Vector<double> >::solveSegregated(Foam::dictionary const&) at ??:?
#5 Foam::fvMatrix<Foam::Vector<double> >::solveSegregatedOrCoupled(Foam::dictionary const&) at ??:?
#6 Foam::fvMesh::solve(Foam::fvMatrix<Foam::Vector<do uble> >&, Foam::dictionary const&) const at ??:?
#7 ? in ~/OpenFOAM/OpenFOAM-v1912/platforms/linux64GccDPInt32Opt/bin/pimpleFoam
#8 __libc_start_main in /lib/x86_64-linux-gnu/
#9 ? in ~/OpenFOAM/OpenFOAM-v1912/platforms/linux64GccDPInt32Opt/bin/pimpleFoam
Floating point exception (core dumped)

Do anyone know what is happening and how could it be solved?

Sorry for my ignorance as it seems to be a very basic question!

Kind regards,
Attached Files
File Type: c blockMeshDict.c (2.9 KB, 1 views)
HeleShawman is offline   Reply With Quote


channel flow, pimplefoam, turbulent flow

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

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
Open Channel Flow - problems with boundaries conditions fabsl FLUENT 0 December 7, 2021 14:54
Enabling Open Channel Flow Sub-Model in Mixture model cod213 FLUENT 0 January 10, 2017 13:40
Problems with validation of LES for turbulent annular flow fgarita OpenFOAM Verification & Validation 0 October 12, 2016 05:37
problem with DNS turbulent channel flow huangxianbei OpenFOAM Running, Solving & CFD 4 March 4, 2014 23:39
Turbulent channel flow by LES Paul Main CFD Forum 14 July 6, 2009 04:02

All times are GMT -4. The time now is 13:24.