# specific heat=f(temperature)-how to do???

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

 June 12, 2003, 09:06 specific heat=f(temperature)-how to do??? #1 Micha 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 gas-mixture 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 1-41 or star-commands-help: 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 m-hack 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 4xF 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 Micha Guest   Posts: n/a Hi 4xF, I got a polynom for cp/R fitted with MS-Excel 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 4xF 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 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 Duan J Q CD-adapco 2 November 20, 2007 07:59 grzesiek FLUENT 2 February 16, 2006 12:18 Andrea CFX 1 May 30, 2004 01:55 gino FLUENT 3 November 23, 2001 00:16 Steve Amphlett Main CFD Forum 0 August 25, 2000 09:51

All times are GMT -4. The time now is 23:09.