total sensible heat transfer rate to determine convergence

  • 44 Views
  • Last Post 25 May 2020
Weiqiang Liu posted this 20 May 2020

Hi all,

I am doing a propane combustion simulation with fluent. Actually I am trying to reproduce literature results. My maximum temperature is 100K lower than literature results. 

I am wondering is this because my case is not converged yet. all the residuals in my case are flat. The monitor variable kept constant as well. Then I checked mass balance of the combustor, it's balanced. Finally I checked the total sensible heat transfer rate. The detailed data set is as follows:

inlet: -58.4W

outlet:-200W

external_wall:-182W

heat of reaction source:461.55 W

net:9.74 W

I know that the net heat transfer rate of the whole system should be 0 if the case is fully converged. However, I can not make it smaller with further iteration.

I am wondering is this an acceptable value to determine convergence or are there any ways to decrease this value further in order to get converged case.

Best regards

Weiqiang

Order By: Standard | Newest | Votes
rahkumar posted this 22 May 2020

Hello, 

The solution not converging, might be the issue. You can reduce the residuals further and run the simulation and check flux balance again. Let us know what you get. Thanks. 

Regards, Rahul

Weiqiang Liu posted this 22 May 2020

Hi, rahkumar,

I can not reduce the residuals by further iterations under the current energy under relaxation factor which is 0.9

If I increase the under relaxation factor of energy to 0.95, the whole system would become unsteady. The maximum temperature just vibrated between 1600K and 2100K back and forth in just tens of iterations.

However, when it vibrated to 2100K which is very close to literature results, the total sensilble heat flux would be very close to zero, which indicates convergence I think. 

I am wondering how can I make this system steady and approach to converged results without oscillation. 

Thanks very much for your concern!

Weiqiang

rahkumar posted this 22 May 2020

Hello, 

What is the solution method you using?

Regards, Rahul

rwoolhou posted this 22 May 2020

If the solution is unstable have you checked the mesh resolution and whether the flame is stable?

Weiqiang Liu posted this 23 May 2020

Hi,

I think I just followed the solution method settings in the paper. Velocity and pressure is coupled with SIMPLE. All terms are discretized with second upwind method.

I don't have access to my computer now. I will confirm the solution method further.

Best regards

weiqiang

Weiqiang Liu posted this 23 May 2020

Hi ,

I have exactly the same mesh with literature. Because the author listed the node number and density in every edge. The author claimed the flame is stable. However, the author also simulated unsteady results and finally the flame would be stable.

Best regards

Weiqiang

Weiqiang Liu posted this 23 May 2020

Hi rahkumar,

Can I ask you a question? how small the net value of sensible heat transfer rate should be at which I can determine that this case is converged according to your experience?

Best regards

Weiqiang

Weiqiang Liu posted this 25 May 2020

hi rahkumar,

I put the screen shot of solution methods panel below. 

Could you please give me some suggestions on how to improve solution methods in order to get converged results.

Best regards

weiqiang

abenhadj posted this 25 May 2020

1% of the most relevant input at your boundaries. Best to have it as small as possible.

Best regards, Amine

Weiqiang Liu posted this 25 May 2020

hi rahkumar,

I set the URF of energy to be greater than 0.93. Then the already-steady residuals starts oscillating like the picture shown below

Weiqiang Liu posted this 25 May 2020

hi abenhadj,

Do you mean the net heat flux should be less than 1% of the relevant input? For example, the heat flux of inlet is 58W, therefore , the net heat flux should be 0.58W? 

Best regards

Weiqiang

abenhadj posted this 25 May 2020

Do not alter the implicit energy URF. Keep it default. 

Switch to Coupled Solver and check if the behavior is enhanced, if you are concerned about convergence and speed-up of steady-state cases.

 

And yes regarding the net heat flux and even lower.

Best regards, Amine

Weiqiang Liu posted this 25 May 2020

Hi Amine,

Thanks very much for your answer! Do you mean that energy URF does not really effect the steady of the whole system? 

Best regards

Weiqiang

abenhadj posted this 25 May 2020

You should not reduce the URF of the energy equation as it has a negative impact on the diffusive load (energy equation is diffusive). You will require a lot of iterations to get back to the full load. 

 

Best regards, Amine

Weiqiang Liu posted this 25 May 2020

Hi Amine,

Thanks very much! By the way, is the default URF of energy 1 in fluent?

Best regards

Weiqiang

abenhadj posted this 25 May 2020

For segregated solver: Yes.

Best regards, Amine

Weiqiang Liu posted this 25 May 2020

Thanks for your patience

Best regards

Weiqiang

Close