CFD Online Discussion Forums (http://www.cfd-online.com/Forums/)
-   FLUENT (http://www.cfd-online.com/Forums/fluent/)
-   -   Residuals querry (http://www.cfd-online.com/Forums/fluent/74098-residuals-querry.html)

 Mohsin March 24, 2010 09:56

Residuals querry

Hello

My residuals are set to 10^-4 for continuity, x,y,z, and turbulent dissipation rate. while doing iterations the residuals goes down just near 10^-4 after 1000 iterations and again moves upwards towards 10^-2 and then doesn't come back.

Is it ok If i stop after 1000 iterations or I should stop the iterations when I m getting constant results (like straight line and unchanged residuals).

Pelase guide me

 -mAx- March 25, 2010 03:01

No you cannot stop iterating after 1000 iterations.
If it drops towards, then there is a problem. You have to find why...

 Mohsin March 25, 2010 03:27

You mean to say that I have to find out why the residuals move from 10^-5 to 10^-3 and goes up and up until 10^-2.

How to find it...any suggestions?

For Velocity, Turbulence What parameters are neeeded to be seen?
Modifying Geometry or mesh? what can be done?

Thanks

 -mAx- March 25, 2010 04:12

grid/check
*check if the mass continuity is conserved (massflow in = massflow out)
...
It is typical debugging, and it's a part of CFD.

 Mohsin March 25, 2010 04:31

Flux report indicate the mass flow rate of all the boundaries= -6.5*10^-7 kg/s. So the continuity seems to conserve.

What else I can see?

 -mAx- March 25, 2010 04:55

-what is the max skewness of your mesh?
-what kind of BC did you used
-is your massflow or pressure at inlet in agreement with your expectations?
etc...

 Mohsin March 25, 2010 05:13

Max skewness = 0.89
At inlet = velocity Inlet
At outlets = Pressure Outlet
Others are wall and fluid.

Ya they are in accordance with my expectations.

Really exhasted:)

 -mAx- March 25, 2010 06:07

*did you compared the results for iterations = 10000, and then for iterations where the residuals drop?
*compressible/uncompressible?
*what kinds of solver/schema/turbulence model

 Mohsin March 26, 2010 03:35

Sorry fro late reply as i was comparing the results of iterations.

At 2500 iterations I got better results as compared with doing 10,000 iterations (at 10,000 iteration the residual plot was constant and was not fluctuating).

I think we should look at the point of minimum in residual plot. After seeing at the residual plots of 10,000 iterations I initialized the calculations and gave the number of iterations until the minumum value (2500 iterations) which i saw from 10,000 iteration plot and got the better result as compared to doing 10,000 iterations.

I m using the following:

Pressure based solver
Incompressible
K-E 2 equation Realizable model
DPM model

 -mAx- March 26, 2010 04:28

ah ok you are running DPM model.
I have 0-experience with it, so I won't be helpful from this side.
Sorry.
I would check tutorial about DPM, if one exists...

 Mohsin March 26, 2010 05:00

http://my.fit.edu/itresources/manual...pdf/ug/pdf.htm

This is Fluent's User guide Chapter 22 is about DPM modeling.

But i was asking generally about the residuals behaviour. once it starts to move up can we stop it there or we should wait for it until it becomes constant.

In my case it became constant until 7000 iterations. But i was getting better results at 2500 iterations. Does it allways supposed to go towards convergence(means the residual plot contnually goes on decreasing), If not then whether something is wrong or its normal.

 -mAx- March 26, 2010 05:07

No the residuals should decrease and then stay constant. But they are not suppose to drop towards.
Monitor the pressure at your inlet, and check its convergence

 Mohsin March 26, 2010 05:13

it means that it should always decrease then remain constant. In my case it first decreased to 10^-5 until 2500 iterations then started to increase until 10^-3 and then remained constant until 10,000 iterations.

So it means that there is some problem somwhere that originates at 2500th iteration which casues my solution to diverge at that point.

I did not specify pressure at the inlet (velcity inlet) . its 0. I just specified the operating pressure.

my Conditions

Max skewness = 0.89
At inlet = velocity Inlet
At outlets = Pressure Outlet
Others are wall and fluid.

??

 -mAx- March 26, 2010 05:44

if you set a velocity at your pressure, the solver will compute the pressure inlet.
Monitor it, to see if the computed pressure converges or not

 Mohsin March 26, 2010 10:01

Welocity is different as pressure. I have just specified the operating pressure of 1 atmosphere. My inlet is velocity inlet whose velocity is 6m/s. However i used pressure outlet at the outlet and there i specified a guage pressure of 0 Pa. As in my case there is no reaction so Inlet and outlet pressures are same. So I dont think i need to specify it.

What do u thinK?

I bothered you so much. Thanks for ur help.

 -mAx- March 26, 2010 10:28

without talking about DPM, if you have a massflow, then you have a drop pressure.
As your outlet-pressure is fixed (your BC), then the drop pressure will be dp= pressure_inlet - pressure_outlet --> that means fluent will adjust the pressure at your inlet.
Maybe with DPM it is different...

 Mohsin March 26, 2010 10:34

I dont know the pressure at outlet. I just know the pressure at inlet which i included in operating pressure. 1 atm + 7 atm-g. I couldnt specify the pressure at inlet coz of velocity inlet boundary condition and on the outlet i dont know the pressure so i just let it be 0 guage at the outlet (the outlet pressure is with relative to operating pressure so as it is 0 it means that outlet pressure at the outlet is 1 atm +7 atm.g.

Thats what i know about the pressure. If u can give me some additional information i will be thankful.

 Mohsin March 26, 2010 10:36

Also please tell me. I am using parallal processing on windows 7 with 1TB hard disk and 4 GB ram Intel core quad 4 cores. So i m using fluent 3d -t4 in CMD but still it is slow. At first it was fast but suddenly it became very slow. I restarted my computer but still the same result. What may have caused this to happen?

 -mAx- March 26, 2010 14:30

Quote:
 Originally Posted by Mohsin (Post 251840) I dont know the pressure at outlet. I just know the pressure at inlet which i included in operating pressure. 1 atm + 7 atm-g. I couldnt specify the pressure at inlet coz of velocity inlet boundary condition and on the outlet i dont know the pressure so i just let it be 0 guage at the outlet (the outlet pressure is with relative to operating pressure so as it is 0 it means that outlet pressure at the outlet is 1 atm +7 atm.g.
If you don't know the pressure at your outlet, then you set at it at 0. You give the massflow at outlet, and fluent computes the pressure at inlet (eg: your drop pressure). Finally your "real" pressure at outlet will be <the pressure inlet you know (1 atm + 7 atm-g)> - <the pressure inlet computed from fluent>

Quote:
 Originally Posted by Mohsin (Post 251841) Also please tell me. I am using parallal processing on windows 7 with 1TB hard disk and 4 GB ram Intel core quad 4 cores. So i m using fluent 3d -t4 in CMD but still it is slow. At first it was fast but suddenly it became very slow. I restarted my computer but still the same result. What may have caused this to happen?
Do a test with another case (a basic one), maybe it comes from your model which has some convergence issue

 Mohsin March 26, 2010 14:41

I changed my opearting system to 32 bit windows 7. It is still slow so it was not the operating systmes problem. I m using 0.4 million cells with parallel processing of 4 cores fluent 3d -t4 but still it is slow. earlier it was fast but i donty know...Its the basic case which i m doing but its ok.

Also The mass flow rate i cannot give at the outlet coz this is what i have to find. Basically i have to know the concentration of particles at the outlet so that they are uniformly distributed among the 3 outlets of the cylinder. I will give non uniform flow at the inlet and check whether the particles can be uniform at the outlet. So i cant give mass flow rate.

My result with a simple case (when the flow at the inlet is also uniform) goes down on the residual plot until 10^-5 and then goes up (as u said it shouldn't go up) but it goes up and then becomes constant. Ssshhh:)

All times are GMT -4. The time now is 11:55.