# negative densities

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

 December 9, 2003, 06:54 negative densities #1 CMB Guest   Posts: n/a Hi I dont know why I am getting this 'negative densities' error everytime I try running my model. Densities in all materials are set as ideal gases, den=p/(rt(sum(mm/Mm))). I have not got a big pressure drops or big temperature changes along my model. I got inlets and outlets at atmospheric pressures. I dont understand, becuase when I do the calculations by myself this should not be happening. Has there anyone who have come across something like this? any suggestions will be appreciated. CM

 December 9, 2003, 10:46 Re: negative densities #2 Frederik Arbeiter Guest   Posts: n/a I get this error in the following setup inlet fixed velocity outlet = outlet rho = f(T,p) high pressure drop I can cure this problem with the following procedure 1. rho=const, run until convergence 2. rho=f(T,p) 3. diminish relaxation factors for p and rho.

 December 9, 2003, 10:48 Re: negative densities #3 Frederik Arbeiter Guest   Posts: n/a I get this error in the following setup inlet fixed velocity outlet = outlet rho = f(T,p) high pressure drop I can cure this problem with the following procedure 1. rho=const, run until convergence 2. rho=f(T,p) 3. diminish relaxation factors for p and rho. In my case down to 0.005 for p and 0.001 for rho -> very slow convergence If there is a better approach I'd be very interested.

 December 9, 2003, 11:30 Re: negative densities #4 CMB Guest   Posts: n/a Thank you, that is quite useful. Indeed decreasing density's and pressure relaxation factor makes it easier to converge. But I got a question, first I have to admit that i have very little idea about the control volume method. Is'nt this decrease in the relaxation factor very high?, I mean does this still valid, as an approximation? in other words is it unrealistic to use this values. I might as well be very wrong. We have noted that convergence is easily obtained by using the AMG method to resolve the pressure equation, you can activate it by 'solu,scalar,AMG' and we have been working all this time with 0.05 pressure r.f. and 0.5 for density, but we might give it a try with lower values and see if it works. Many thanks CM

 January 4, 2004, 03:41 Re: negative densities #5 azmir Guest   Posts: n/a Hello! Without changing your physical BCs, you could try pseudo-transient or single-transient with very small time steps. I found out that this is a way to resolve my neg. density issue for compressible flow case. A big jump in grid sizes (e.g. from small grid size to very big grid size at coupled interphase) might also cause the error.

 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 john Siemens 0 March 19, 2008 07:32 Riyaz Ahamed Siemens 2 August 14, 2007 01:47 Chris Siemens 4 June 25, 2007 08:39 shyam Siemens 2 May 8, 2005 06:17 Mihir Siemens 4 December 9, 2002 22:04

All times are GMT -4. The time now is 17:45.