
[Sponsors] 
June 12, 2003, 09:06 
specific heat=f(temperature)how to do???

#1 
Guest
Posts: n/a

Hello,
I would like to model the specific heat as a function of temperature. Therefore I calculated the specific heat of my gasmixture at different points with an other programm, calculated the quotient Cp/R according to the manual and also fitted a polynom to these values (4th degree). So I got 5 polynomial coefficients. I tried to define the polynom with command popcp (without the n+1 coefficient), but Unfortunately PROSTAR calculates the wrong values for the specific heat...! Probably the reason is, that I have to define another coefficient for enthalpy (>methology equation 141 or starcommandshelp: popcp). Why do I need this one? Where can I get it from? What do I have to do to define the polynom for cp=f(t) correctly? Thanks for every advice! Micha 

June 12, 2003, 10:21 
Re: specific heat=f(temperature)how to do???

#2 
Guest
Posts: n/a

Try a Fortan code in the CONVTE.f subroutine. You define both enthalpy and cp there.
Greetings, Marc 

June 13, 2003, 06:29 
Re: specific heat=f(temperature)how to do???

#3 
Guest
Posts: n/a

Hmmm.. Quite tricky this one... The polynomial definition is set accordingly to the definition of H/RT used by the Chemkin database. The relation to be used corresponds to the following formula:
H/RT = Cp/R = C1 + C2 T /2 + C3 T^2 /3 + C4 T^3 /4 + C5 T^4 /5 + C6/T So, you have to calculate your coefficients so that they match the definition above. 

June 13, 2003, 07:48 
Re: specific heat=f(temperature)how to do???

#4 
Guest
Posts: n/a

Hi 4xF,
I got a polynom for cp/R fitted with MSExcel to something like Cp/R=a1+a2*T+a3*T^2+a4*T^3+a5*T^4 with cp...specific heat capacity and R...specific gas constant. So according to you, I can't use this simple polynomial fitting, but I have to do another one which includes also the term: a6/T . Then I have to calculate my values for PROSTAR like this: C1=a1 C2=a2*2 C3=a3*3 C4=a4*4 C5=a5*5 C6=a6 Is that correct? I just want to make sure I got it...! Thanks. Micha 

June 13, 2003, 10:07 
Re: specific heat=f(temperature)how to do???

#5 
Guest
Posts: n/a

This is correct, when you do it like that. You can set a6 to 0. if you want to. a6 should be used to calculate the reference enthalpy at T =298.15. For a pure specie, this value should give you a value of H/RT that is nearly 0.


Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
specific heat and mean specific heat  Duan J Q  CDadapco  2  November 20, 2007 07:59 
specific heat  grzesiek  FLUENT  2  February 16, 2006 12:18 
Specific heat  Andrea  CFX  1  May 30, 2004 01:55 
specific heat in UDF  gino  FLUENT  3  November 23, 2001 00:16 
Integral specific heat  Steve Amphlett  Main CFD Forum  0  August 25, 2000 09:51 