# successful time step for 2D problem failing for 3D case of same problem.

 User Name Remember Me Password
 Register Blogs Members List Search Today's Posts Mark Forums Read

 March 9, 2015, 08:02 successful time step for 2D problem failing for 3D case of same problem. #1 New Member   Manoj Dhiman Join Date: Dec 2014 Posts: 4 Rep Power: 11 Hello, I am trying to simulate 3D convective diffusive unsteady problem using SIMPLER technique (Finite Volume Method). I got successful 2D case of the problem. For 3D I am keeping the size of space step same as in 2D, with one extra direction. The problem is in choosing the time step. When I keep it same as in 2D the simulation diverges and in order to converge it I have to reduce the time step by order of 4. Is this thing common or I am committing some other problem. Unluckily if its common than it will take ages for this simulation. BTW I am running this code with gfortran compilers.

March 9, 2015, 11:17
#2
Senior Member

Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,815
Rep Power: 73
Quote:
 Originally Posted by mandman Hello, I am trying to simulate 3D convective diffusive unsteady problem using SIMPLER technique (Finite Volume Method). I got successful 2D case of the problem. For 3D I am keeping the size of space step same as in 2D, with one extra direction. The problem is in choosing the time step. When I keep it same as in 2D the simulation diverges and in order to converge it I have to reduce the time step by order of 4. Is this thing common or I am committing some other problem. Unluckily if its common than it will take ages for this simulation. BTW I am running this code with gfortran compilers.

yes, the numerical stability involves the relation of dimensionality therefore the constraint is more strict when going from 1D, 2D, 3D.

 March 9, 2015, 21:32 #3 Senior Member   adrin Join Date: Mar 2009 Posts: 115 Rep Power: 17 Though it is true that the stability constraint for 3D is stricter than in 2D, it is difficult to believe that the time step size would have to drop by 4 _orders_. By 4 _orders_ I understand you're saying 10^4 (10,000 times). Otherwise, if the drop is _factor_ 4 (4 times) it is somewhat more reasonable, though it depends (very much) on the problem and what physics are introduced by the addition of the third dimension (and 4x increase doesn't really translate to "ages" of simulation time in added cost) adrin

 March 10, 2015, 08:31 #4 Senior Member   Join Date: Jul 2009 Posts: 357 Rep Power: 19 If you are saying that the time step has to be 10000 times smaller rather than 4 times smaller, I would look for an error in the coding. A factor of 10000 is not reasonable for a jump to 3D, in my experience.

 March 10, 2015, 14:18 #5 New Member   Manoj Dhiman Join Date: Dec 2014 Posts: 4 Rep Power: 11 by O(4) I mean 10000.

 March 10, 2015, 14:40 #6 Senior Member   adrin Join Date: Mar 2009 Posts: 115 Rep Power: 17 Then you most likely (definitely?) have a bug in the code! What are the 2- and 3-D problems you were solving and what are the corresponding CFL values? adrin

 March 10, 2015, 19:01 #7 New Member   Manoj Dhiman Join Date: Dec 2014 Posts: 4 Rep Power: 11 Ya!! you guys are right. I found THE "bug", it was resting in expression calculating the difference between pressure and pressure correction. I was struggling with this problem for last 2 months, and I was so much confident about my code that I thought that its common for 3D problems. But you guys conviction made me doubtful. I even learned OpenMP and parallelised this code, but there’s no way I could have compensated for O(4) time increase. Thanks all I am working on double diffusive convection problem, let me know if you want to know more. You can email me on manoj_d@students.iitmandi.ac.in

 Thread Tools Search this Thread Search this Thread: Advanced Search Display Modes Linear Mode

 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 OffTrackbacks are Off Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post skuznet OpenFOAM Running, Solving & CFD 6 November 15, 2017 12:12 Artex85 OpenFOAM Running, Solving & CFD 11 February 16, 2017 13:40 setzberg OpenFOAM Running, Solving & CFD 6 October 15, 2015 13:29 msrinath80 OpenFOAM Running, Solving & CFD 18 March 3, 2015 05:36 Forrest_Lei OpenFOAM 3 July 19, 2011 06:00

All times are GMT -4. The time now is 16:37.

 Contact Us - CFD Online - Privacy Statement - Top