Problem in 3D extension of pitzDaily

yfjok22 · October 19, 2016, 11:12

Dear all,

For the testing of new model, a 3D incompressible turbulent case would be needed. To be simple, I'm planning to extend the formal case of simpleFoam, pitzDaily to 3D. However, some problems have been met

Actually, I just extend the geometry of pitzDaily in z direction from 1.0 to 10.0, and also discrete it from 1 to 10. Except that, no further modification of mesh, initial value and boundary condition. However, when I run this case with simpleFoam, it always get diverged, the output would be like this

Time = 0.5

smoothSolver: Solving for Ux, Initial residual = 0.999999, Final residual = 0.0801313, No Iterations 7
smoothSolver: Solving for Uy, Initial residual = 0.999998, Final residual = 0.080132, No Iterations 7
GAMG: Solving for p, Initial residual = 1, Final residual = 0.0977397, No Iterations 3
time step continuity errors : sum local = 3.01736e+113, global = 4.62262e+112, cumulative = 4.62262e+112
#0 Foam::error:rintStack(Foam::Ostream&) at ??:?
#1 Foam::sigFpe::sigHandler(int) at ??:?
#2 ? in "/lib/x86_64-linux-gnu/libc.so.6"
#3 Foam::multiply(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) at ??:?
#4 Foam::tmp<Foam:imensionedField<double, Foam::volMesh> > Foam:perator*<Foam::volMesh>(Foam::tmp<Foam:im ensionedField<double, Foam::volMesh> > const&, Foam:imensionedField<double, Foam::volMesh> const&) at ??:?
#5 Foam::RASModels::kEpsilon<Foam::IncompressibleTurb ulenceModel<Foam::transportModel> >::correct() at ??:?
#6 ? at ??:?
#7 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#8 ? at ??:?
Floating point exception (core dumped)

Based on the output, I think this should be some problem from the turbulent model setting.
So guys, would you give me some suggestion, after extending this pitzDaily case to 3D, what kind of modification would be need for the k-epsilon model, or the problem would be something else?
Thanks a lot and look forward to replies^^

yfjok22 · October 19, 2016, 11:39

I have figured out, it's because I didn't change faces in z direction from empty to wall
But again, I got a problem
If I didn't make any local mesh densification of the mesh near wall, but just make the simpgrading equals 1 in this direction, will it bring some problem to the turbulence model calculating like k-epsilon

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
area does not match neighbour by ... % -- possible face ordering problem	St.Pacholak	OpenFOAM	10	February 7, 2024 21:50
UDF compiling problem	Wouter	Fluent UDF and Scheme Programming	6	June 6, 2012 04:43
Problem in implementing cht	tilek	CFX	3	May 8, 2011 08:39
natural convection problem for a CHT problem	Se-Hee	CFX	2	June 10, 2007 06:29
Adiabatic and Rotating wall (Convection problem)	ParodDav	CFX	5	April 29, 2007 19:13

October 19, 2016, 11:12	Problem in 3D extension of pitzDaily	#1
yfjok22 New Member Fujiang Yu Join Date: Jan 2014 Posts: 12 Rep Power: 12	Dear all, For the testing of new model, a 3D incompressible turbulent case would be needed. To be simple, I'm planning to extend the formal case of simpleFoam, pitzDaily to 3D. However, some problems have been met Actually, I just extend the geometry of pitzDaily in z direction from 1.0 to 10.0, and also discrete it from 1 to 10. Except that, no further modification of mesh, initial value and boundary condition. However, when I run this case with simpleFoam, it always get diverged, the output would be like this Time = 0.5 smoothSolver: Solving for Ux, Initial residual = 0.999999, Final residual = 0.0801313, No Iterations 7 smoothSolver: Solving for Uy, Initial residual = 0.999998, Final residual = 0.080132, No Iterations 7 GAMG: Solving for p, Initial residual = 1, Final residual = 0.0977397, No Iterations 3 time step continuity errors : sum local = 3.01736e+113, global = 4.62262e+112, cumulative = 4.62262e+112 #0 Foam::error:rintStack(Foam::Ostream&) at ??:? #1 Foam::sigFpe::sigHandler(int) at ??:? #2 ? in "/lib/x86_64-linux-gnu/libc.so.6" #3 Foam::multiply(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) at ??:? #4 Foam::tmp<Foam:imensionedField<double, Foam::volMesh> > Foam:perator<Foam::volMesh>(Foam::tmp<Foam:im ensionedField<double, Foam::volMesh> > const&, Foam:imensionedField<double, Foam::volMesh> const&) at ??:? #5 Foam::RASModels::kEpsilon<Foam::IncompressibleTurb ulenceModel<Foam::transportModel> >::correct() at ??:? #6 ? at ??:? #7 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6" #8 ? at ??:? Floating point exception (core dumped) Based on the output, I think this should be some problem from the turbulent model setting*. So guys, would you give me some suggestion, after extending this pitzDaily case to 3D, what kind of modification would be need for the k-epsilon model, or the problem would be something else? Thanks a lot and look forward to replies^^

October 19, 2016, 11:39		#2
yfjok22 New Member Fujiang Yu Join Date: Jan 2014 Posts: 12 Rep Power: 12	I have figured out, it's because I didn't change faces in z direction from empty to wall But again, I got a problem If I didn't make any local mesh densification of the mesh near wall, but just make the simpgrading equals 1 in this direction, will it bring some problem to the turbulence model calculating like k-epsilon