CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > FLUENT

Why I get unbounded solution with bounded scheme?

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

Reply
 
LinkBack Thread Tools Display Modes
Old   August 24, 2012, 23:54
Question Why I get unbounded solution when I shouldn't?
  #1
Member
 
George Pichurov
Join Date: Jul 2010
Posts: 39
Rep Power: 7
jorkolino is on a distinguished road
I simulate a scalar field with inlet value equal to 1 and no sources/sinks in the domain (there is one more inlet with zero scalar value). Naturally, I'd expect to have solution within [0;1], however I get some negative values, the same order of magnitude as the positives, in some inlet-adjacent cells Most curious, the solution converges to a high accuracy (1e-15), which means it should be correct WTF?!! I use first order upwind for the scalar, nothing special to note, apart from that inlet diffusion is disabled. But even enabling it only diminished the problem, did not eliminate it.

I've noticed that the scalar diffusion coefficient plays role. Setting higher diff. coeff. leads to higher oscillation, i.e. more and higher magnitude negative values. Yet again, solution converges perfectly. I;ve tried relaxing, but as you guess, it helps convergence, which is no problem by itself, and does not change the final converged solution.

Last edited by jorkolino; August 25, 2012 at 15:53.
jorkolino is offline   Reply With Quote

Old   August 26, 2012, 21:53
Default
  #2
Member
 
Touré
Join Date: Aug 2010
Location: Montréal
Posts: 85
Rep Power: 7
Touré is on a distinguished road
Show the geometry. You may have recirculation.
Touré is offline   Reply With Quote

Old   August 27, 2012, 00:28
Default
  #3
Member
 
George Pichurov
Join Date: Jul 2010
Posts: 39
Rep Power: 7
jorkolino is on a distinguished road
Here is the inlet colored by scalar cell values after convergence. The inlet supplies fixed flow velocity (0.0113m/s) with fixed value of scalar (=1). I don't get how circulation can change the values to negatives as no parts of the domain have sinks. Nothing can decrease value below zero, physically. The only possible twist is I've disabled inlet diffusion. So no diffusive flux from the inlet is expected, only convective. But still that should only possibly increase the scalar value in the nearby cells.
Attached Images
File Type: jpg 6hu-swf-inlet.jpg (55.7 KB, 11 views)
File Type: jpg 6hu-swf-inlet2.jpg (44.9 KB, 5 views)
jorkolino is offline   Reply With Quote

Old   August 27, 2012, 12:40
Default
  #4
Member
 
Touré
Join Date: Aug 2010
Location: Montréal
Posts: 85
Rep Power: 7
Touré is on a distinguished road
What is your scalar function? I work with the level set method and you need to initialize the scalar level set function everywhere in the domain (volume) including all boundaries at the initial time step, but not only on the inlet boundaries. At the initial, the scalar function varies between 0 and 1, but when it evolves in transient (unsteady), it may have negative values depending of the evolution of the flow or for example if it's an outlet. The maximum value of you scalar is 0.0981 instead of 1. Is the value of the scalar of the inlet really fixed to 1?
Touré is offline   Reply With Quote

Old   August 27, 2012, 17:42
Default
  #5
Member
 
George Pichurov
Join Date: Jul 2010
Posts: 39
Rep Power: 7
jorkolino is on a distinguished road
I solve steady state. The value of the boundary is 1, as is the boundary condition, but you won't see it on the pictures, becuse they show cell center values after convergence. I have tried different strategies to initialize, with 0, 1 and in between for all cells in the domain. With the same result. It simply converges to what you see. Negative values at some cells adjacent to inlet. And it converges well, reeeeally well (res 1e-15). I continue to think that the disabled inlet diffusion (which you control via the UDS panel) is the culprit, but I can not see why and how it will create negatives of that order.
jorkolino is offline   Reply With Quote

Reply

Thread Tools
Display Modes

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 On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
grid dependancy gueynard a. Main CFD Forum 19 June 27, 2014 21:22
Bounded scheme gives unbounded solution su_junwei OpenFOAM Running, Solving & CFD 3 November 18, 2011 03:35
Doubt on Implicit Methods analyse In India Main CFD Forum 10 March 9, 2007 04:01
Could anybody help me see this error and give help liugx212 OpenFOAM Running, Solving & CFD 3 January 4, 2006 19:07
Wall functions Abhijit Tilak Main CFD Forum 6 February 5, 1999 02:16


All times are GMT -4. The time now is 02:28.