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

 June 15, 2014, 19:02 offsetCylinder tutorial, steady state, #1 New Member   Riccardo Join Date: May 2014 Posts: 10 Rep Power: 5 Greetings all, I am trying to run the offsetCylinder "tutorial" with different non-newtonian fluids at flow velocities such that the regimen stays laminar. So, by modifying only the viscosity parameters I am able to run the simulation and to obtain a visual representation of a state that appears to be laminar and steady. However from the log of the simulation only one component of velocity appears to have converged after the simulation. This means that after a certain timestep, only that component of velocity shows 0 iterations. I ran the simulation for a long time so I don't think the simulation duration is of any importance anymore. 1-Is this a normal behavior of the nonNewtonianIcoFoam solver? 2-If yes, without performing new simulations with steady state solvers, how can I understand at what approximate timestep a steady state is acceptably reached (merely for presentation purpouses)? 3-Unfortunately C syntax is still misterious to me so I am not really aware of what I'm doing most of the time. Could someone sum up what is going on when running nonNewtonianIcoFoam, I mean from a mathematical/numerical point of view? Thank you Last edited by tabsquare; June 15, 2014 at 20:02.

June 16, 2014, 11:08
#2
New Member

Riccardo
Join Date: May 2014
Posts: 10
Rep Power: 5
I attach a plot of the residuals. I am not really sure how to interpret the results. Indeed this is a case where both components of velocity converge. But what about pressure?
Attached Images
 residuals.jpg (52.3 KB, 11 views)

 June 18, 2014, 12:27 #3 Senior Member   Joachim Herb Join Date: Sep 2010 Posts: 426 Rep Power: 13 The pressure equations is solved twice (or more often, probably 10 times) per solution of the other equations because you have set (in fvSolutions) either nCorrectors bigger than 1 or nNonOrthogonalCorrectors bigger than 0. tabsquare likes this.

July 15, 2014, 12:45
#4
New Member

Riccardo
Join Date: May 2014
Posts: 10
Rep Power: 5
OK still I have no good knowledge of the math behind this, but I get what you meant to some extent.
I have another question. I attached a plot of velocity residuals as a function of time. Pressure residuals show a similar behavior.
Final residuals of velocity decrease until "tolerance" is reached, then the curve of final residuals merges with the Initial residuals one and both keep decreasing until the end of the simulation(manually stopped) is reached. To what solver is this second decrease due to?

I am using the same files contained in the incompressible/nonNewtonianIcoFoam/offsetCylinder/ tutorial. I only made some modification to transportProperties and controlDict.
Attached Images
 Selection_042.jpg (20.5 KB, 8 views)

 Thread Tools 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 On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post Chander Main CFD Forum 5 December 23, 2013 06:31 Sas CFX 15 July 13, 2010 08:56 Kushagra CFX 1 June 22, 2008 19:06 Lisa Main CFD Forum 11 July 5, 2000 14:37 Jitendra Main CFD Forum 1 June 27, 2000 18:49

All times are GMT -4. The time now is 20:43.