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

Problem in 3D extension of pitzDaily

Register Blogs Community New Posts Updated Threads Search

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   October 19, 2016, 12:12
Default Problem in 3D extension of pitzDaily
  #1
New Member
 
Fujiang Yu
Join Date: Jan 2014
Posts: 12
Rep Power: 12
yfjok22 is on a distinguished road
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 is offline   Reply With Quote

Old   October 19, 2016, 12:39
Default
  #2
New Member
 
Fujiang Yu
Join Date: Jan 2014
Posts: 12
Rep Power: 12
yfjok22 is on a distinguished road
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
yfjok22 is offline   Reply With Quote

Reply

Tags
k-epsilon model, settings, turbulence


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
area does not match neighbour by ... % -- possible face ordering problem St.Pacholak OpenFOAM 11 September 4, 2024 05:28
UDF compiling problem Wouter Fluent UDF and Scheme Programming 6 June 6, 2012 05:43
Problem in implementing cht tilek CFX 3 May 8, 2011 09:39
natural convection problem for a CHT problem Se-Hee CFX 2 June 10, 2007 07:29
Adiabatic and Rotating wall (Convection problem) ParodDav CFX 5 April 29, 2007 20:13


All times are GMT -4. The time now is 14:55.