|
[Sponsors] |
October 19, 2016, 12:12 |
Problem in 3D extension of pitzDaily
|
#1 |
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, 12:39 |
|
#2 |
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 |
|
Tags |
k-epsilon model, settings, turbulence |
|
|
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 |