# Why I get unbounded solution with bounded scheme?

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

 August 24, 2012, 23:54 Why I get unbounded solution when I shouldn't? #1 Member   George Pichurov Join Date: Jul 2010 Posts: 40 Rep Power: 9 Sponsored Links 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.

 August 26, 2012, 21:53 #2 Member   Guiliguili Join Date: Aug 2010 Location: Montréal Posts: 94 Rep Power: 9 Show the geometry. You may have recirculation.

August 27, 2012, 00:28
#3
Member

George Pichurov
Join Date: Jul 2010
Posts: 40
Rep Power: 9
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
 6hu-swf-inlet.jpg (55.7 KB, 11 views) 6hu-swf-inlet2.jpg (44.9 KB, 5 views)

 August 27, 2012, 12:40 #4 Member   Guiliguili Join Date: Aug 2010 Location: Montréal Posts: 94 Rep Power: 9 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?

 August 27, 2012, 17:42 #5 Member   George Pichurov Join Date: Jul 2010 Posts: 40 Rep Power: 9 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.

 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 gueynard a. Main CFD Forum 19 June 27, 2014 21:22 su_junwei OpenFOAM Running, Solving & CFD 3 November 18, 2011 03:35 analyse In India Main CFD Forum 10 March 9, 2007 04:01 liugx212 OpenFOAM Running, Solving & CFD 3 January 4, 2006 19:07 Abhijit Tilak Main CFD Forum 6 February 5, 1999 02:16