
[Sponsors] 
October 28, 2008, 17:22 
Energy balance in a pistoncylinder system

#1 
Guest
Posts: n/a

Hi, I am simulating the behaviour of oscillating flow and heat transfer in a pistoncylinder system using CFX11.0. The 2D axisymmetric model includes a moving boundary and isothermal boundary conditions are applied. My problem is that the energy balance of a whole cycle can't be reached while temperature and pressure seem to get to a steady state. In a simple pistoncylinder system, the net PV work during a cycle should be equal to the neat heat in and out of the system as the system gets to a steady state. In my cases, the net heat is nearly 2~9 times larger than the net PV work. It is noted that the net PV work is calculated by integrating PdV over a cycle and the net heat is calculated by integrating the heat flux over whole boundary area using CFXpost function and time period. I had simulated the MIT single space experiments using CFX. The results of pressure agree well with the experimental data, but the net loss is a few times larger than the net PV work over a cycle.
Any help will be highly appreciated. Best regards Taotao 

October 29, 2008, 08:17 
Re: Energy balance in a pistoncylinder system

#2 
Guest
Posts: n/a

How many cycles you have done in your calculation? If energy balance is not reached, net heat flux number is bogus. Is mass flow balanced? Temperature reaches stable condition while energy is not balanced. Things do not add up.


October 29, 2008, 09:40 
Re: Energy balance in a pistoncylinder system

#3 
Guest
Posts: n/a

1. Are you solving for "Total Energy" or "Thermal Energy"? 2. Have you converged within each timestep? 3. How many coefficient loops is the solver using? Does it reach the limit of coefficient loops you have provided? 4. What are you calculating in Post?


October 29, 2008, 16:05 
Re: Energy balance in a pistoncylinder system

#4 
Guest
Posts: n/a

Hi Longbow
Thank you for your response. I had calculated 120 cycles. I don't know how to check mass flow balance because this is a closed system. Regards Taotao 

October 29, 2008, 16:14 
Re: Energy balance in a pistoncylinder system

#5 
Guest
Posts: n/a

Hi Cyclone,
Thank you for your response. 1. I am solving for Total energy. 2. The solution is converged under RMS 1e6 3. The number of coefficient loops is 10 and it converged less than the limit. 4. I am using the function of AreaInv(heat flux)@boundary and then integrate it over the whole cycle to calculate the net heat loss. Regards Tonny 

October 29, 2008, 21:26 
Re: Energy balance in a pistoncylinder system

#6 
Guest
Posts: n/a

Hi,
Are you using a moving mesh with second order time differencing? There has been some problems with heat conservation with fully compressible fluids in moving mesh simulations with second order timestepping in previous versions of CFX, not sure if the current version still has the problem. Glenn Horrocks 

October 29, 2008, 21:53 
Re: Energy balance in a pistoncylinder system

#7 
Guest
Posts: n/a

Hi Glen,
Thank you for your response. Yes, I am using a moving mesh with second order time difference. I am looking at heat transfer in a pistoncylinder system by using CFX11.0. Any suggestion for dealing with this problem? Regards Tonny 

October 30, 2008, 20:13 
Re: Energy balance in a pistoncylinder system

#8 
Guest
Posts: n/a

Hi,
Not recommended. Use first order time differencing. At least try it and see if it gives you better conservation. Glenn Horrocks 

October 30, 2008, 21:32 
Re: Energy balance in a pistoncylinder system

#9 
Guest
Posts: n/a

Hi Glenn,
Thank your very much for your help. I am going to try it using first order time difference. Why is CFX not recommended for this kind of problem? Is it mainly due to moving mesh? Which commercial package such as Fluent, Star CD is more suitable for this kind of problem? Tonny 

October 31, 2008, 01:23 
Re: Energy balance in a pistoncylinder system

#10 
Guest
Posts: n/a

Hi,
No, the "not recommended" refers to the combination of second order time stepping, moving mesh and compressible flow. The easy fix is to go to first order time stepping. CFX can model this type of flow just fine, as long as you are careful and do the normal checks (which you have to do with any software anyway). Glenn Horrocks 

October 31, 2008, 08:26 
Re: Energy balance in a pistoncylinder system

#11 
Guest
Posts: n/a

You don't have inlet and outlet in your model? In other words, you have a completely closed cylinder with piston compressing and decompressing fluid inside? Is piston isothermal boundary too?


October 31, 2008, 09:32 
Re: Energy balance in a pistoncylinder system

#12 
Guest
Posts: n/a

I forgot about that limitation. This has been fixed for version 12, so second order time discretization can be used.


October 31, 2008, 21:14 
Re: Energy balance in a pistoncylinder system

#13 
Guest
Posts: n/a

Hi Longbow,
No, I don't have inlet and outlet boundary. The piston is isothermal boundary too. Thanks Tonny 

November 4, 2008, 16:59 
Re: Energy balance in a pistoncylinder system

#14 
Guest
Posts: n/a

Hi Glenn,
I calculated it using firstorder time differencing scheme which gave better results than secondorder time differencing scheme, but the results doesn't seem to be satisfying while you look at heat transfer in the system. Tonny 

November 4, 2008, 17:03 
Re: Energy balance in a pistoncylinder system

#15 
Guest
Posts: n/a

Hi Cyclone,
I tried version 12.0 (beta) with second order time difference scheme. The results for energy conservation is similar to those from version 11.0 and the convergence is even worse. Tonny 

November 4, 2008, 18:00 
Re: Energy balance in a pistoncylinder system

#16 
Guest
Posts: n/a

Hi,
It was found the area of side wall is constant while I extracted the area of side wall from CFXpost using the internal function of area()@, but the volume of domain varies with time.It doesn't make sense. Have you noticed that problem? Many thanks Tonny 

November 4, 2008, 19:18 
Re: Energy balance in a pistoncylinder system

#17 
Guest
Posts: n/a

Hi,
What do you mean "results don't seem to be satisfying"? You mean you have a heat imbalance? Have you done a convergence sensitivity check? Have you done a timestep sensitivity check? I bet you have too big timesteps and/or incomplete convergence. Glenn Horrocks 

November 10, 2008, 16:55 
Re: Energy balance in a pistoncylinder system

#18 
Guest
Posts: n/a

Hi Glenn,
Sorry for slow response. I have been doing timestep sensitivity check recently using first order time differencing. I found the net heat flow of one cycle strongly depend on the timestep. Up to now, I still have a heat imbalance of 30%. In addition, it make convergence more difficult using first order time differencing. It was found the area of side wall is constant while I extracted the area of side wall from CFXpost using the internal function of area()@, but the volume of domain varies with time.It doesn't make sense. Have you noticed that problem? Thank you very much for your help Taotao 

Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
help to simulate 4stroke piston cylinder  san  Siemens  2  September 29, 2006 13:37 
HELP needed: Energy balance in gaseous combustion  James Willie  FLUENT  11  August 10, 2005 04:58 
Energy balance  Julie Polyakh  Siemens  9  July 7, 2004 04:00 
Energy Balance  Cholo  FLUENT  0  April 9, 2004 13:53 
Why FVM for highRe flows?  Zhong Lei  Main CFD Forum  23  May 14, 1999 13:22 