CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > ANSYS > CFX

Porous domain energy calculation and specific heat capacity temperature dependent

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

Reply
 
LinkBack Thread Tools Display Modes
Old   November 28, 2017, 13:27
Default Porous domain energy calculation and specific heat capacity temperature dependent
  #1
Member
 
VB
Join Date: Jul 2016
Posts: 35
Rep Power: 3
Vishnu_bharathi is on a distinguished road
Hello people,

I have a simple model where porous domain is inbetween two fluid domains. Hot fluid is given as Input. I am using Specific heat capacities for fluid and solid material in porous domain to be temperature dependent. Also the density of fluid is temperature dependent. Energy calculation formula used 'Integral over (mass_flow*C_p*DeltaT )dt. Is there a way to calculate C_p for each cell?(in ansys cfx pre or post).

The opearating temperature range is about 800C, so I do not want to take single average C_p value as it varies quite significantly between the temp range and compromises the Energy calculation.

In function Calculator I tried Function:volumeInt , Location: porous, Variable:Specific Heat Capacity at Constant Pressure. What volumeInt(SHC)@Domain returns actually?

Ultimately I need to calculate Energy in porous region/solid, where it is calculated using C_p that is temperature dependent. Any suggestion?
Vishnu_bharathi is offline   Reply With Quote

Old   November 28, 2017, 17:29
Default
  #2
Super Moderator
 
Glenn Horrocks
Join Date: Mar 2009
Location: Sydney, Australia
Posts: 14,321
Rep Power: 110
ghorrocks is a jewel in the roughghorrocks is a jewel in the roughghorrocks is a jewel in the roughghorrocks is a jewel in the rough
In the material model make Cp/specific heat a function of temperature. It is simple as that. But be aware that if you are specifying density as well you need to make it thermodynamically consistent across the temperature range.
ghorrocks is offline   Reply With Quote

Old   November 28, 2017, 17:45
Default
  #3
Member
 
VB
Join Date: Jul 2016
Posts: 35
Rep Power: 3
Vishnu_bharathi is on a distinguished road
Quote:
Originally Posted by ghorrocks View Post
In the material model make Cp/specific heat a function of temperature. It is simple as that. But be aware that if you are specifying density as well you need to make it thermodynamically consistent across the temperature range.
Hi. What do you mean by make specific heat, a function of temperature in material model? I have specified Cp as temperature dependent in expression -CFX Pre for porous solid domain but I am sure you are not specifing this, are you?

I have completed my simulation run using C_p of Porous domain as temperature dependent(C_p is temperature dependent for both solid and fluid materials in porous). I need the formula 'Integral over (mass_flow*C_p*DeltaT )dt for Energy calculation in porous domain where C_p is not constant. To put in other way, Is there a way to access C_p during each iteration/each node so that energy of total porous domain is calculated at each time step or at the end of simulation, either one is fine?
Vishnu_bharathi is offline   Reply With Quote

Old   November 28, 2017, 17:54
Default
  #4
Senior Member
 
Join Date: Jun 2009
Posts: 732
Rep Power: 15
Opaque is on a distinguished road
Not clear what you are trying to compute?

When working with varying specific heat capacities, you are better off using energy quantities directly than simplified formulas unless you want to keep the thermodynamics consistency yourself.

For example, the integral over volume of static enthalpy, i.e. volumeInt ( density * static enthalpy)@Volume of Interest
Opaque is offline   Reply With Quote

Old   November 28, 2017, 18:06
Default
  #5
Member
 
VB
Join Date: Jul 2016
Posts: 35
Rep Power: 3
Vishnu_bharathi is on a distinguished road
Quote:
Originally Posted by Opaque View Post
Not clear what you are trying to compute?

When working with varying specific heat capacities, you are better off using energy quantities directly than simplified formulas unless you want to keep the thermodynamics consistency yourself.

For example, the integral over volume of static enthalpy, i.e. volumeInt ( density * static enthalpy)@Volume of Interest
Hi. I need to calculate Energy contained in porous domain after my simulation or during the run, either one is fine. "you are better off using energy quantities directly than simplified formulas" - I would love to keep it simple . The formula that I mentioned is used to compute Energy in each node and in turn the whole domain. (eg., in Open foam or similar software where we have access to each node/cell in a domain, unlike Ansys). I believe Ansys is even simpler and accurate in Energy calculations over interested domain. What this returns? "volumeInt ( density * static enthalpy)@Volume of Interest"
Vishnu_bharathi is offline   Reply With Quote

Old   November 28, 2017, 18:12
Default
  #6
Senior Member
 
Erik
Join Date: Feb 2011
Location: Earth (Land portion)
Posts: 718
Rep Power: 14
evcelica is on a distinguished road
Quote:
Originally Posted by ghorrocks View Post
But be aware that if you are specifying density as well you need to make it thermodynamically consistent across the temperature range.
Pay attention to this, you can't specify both density and specific heat as variable unless you are very sure you have no thermodynamic inconsistencies.

Ways to get around this include:
Use the Boussinesq approximation for buoyant flow,
Use the ideal gas model,
Use a real gas model
Use an RGP table.

You can get some weird stuff happening with the energy equations if you specify both density and Specific heat and functions of temperature.
evcelica is offline   Reply With Quote

Old   November 28, 2017, 18:25
Default
  #7
Member
 
VB
Join Date: Jul 2016
Posts: 35
Rep Power: 3
Vishnu_bharathi is on a distinguished road
Quote:
Originally Posted by evcelica View Post
Pay attention to this, you can't specify both density and specific heat as variable unless you are very sure you have no thermodynamic inconsistencies.

Ways to get around this include:
Use the Boussinesq approximation for buoyant flow,
Use the ideal gas model,
Use a real gas model
Use an RGP table.

You can get some weird stuff happening with the energy equations if you specify both density and Specific heat and functions of temperature.
Hi Evcelica, Thank you for the suggestion. I have defined density and specific heat in such a way it does'nt behave weird in the temp range. I believed so. now the simulation has been completed. How do you calculate Energy in each domains now? I need this Energy value to cross check and few other stuffs.
Vishnu_bharathi is offline   Reply With Quote

Old   November 29, 2017, 05:29
Default VolumeInt , Energy
  #8
Member
 
VB
Join Date: Jul 2016
Posts: 35
Rep Power: 3
Vishnu_bharathi is on a distinguished road
Quote:
Originally Posted by Vishnu_bharathi View Post
Hi. I need to calculate Energy contained in porous domain after my simulation or during the run, either one is fine. "you are better off using energy quantities directly than simplified formulas" - I would love to keep it simple . The formula that I mentioned is used to compute Energy in each node and in turn the whole domain. (eg., in Open foam or similar software where we have access to each node/cell in a domain, unlike Ansys). I believe Ansys is even simpler and accurate in Energy calculations over interested domain. What this returns? "volumeInt ( density * static enthalpy)@Volume of Interest"
Update: In order to verify this expression before using in Porous domain: I used volumeInt(density*static enthalpy)@my domain and compared with Energy calculated by m(mass)*C_p*dt. Both have difference. Why do I have difference here? Note: Here, the domain has constanct specific heat capacity and density, its a solid domain.
Vishnu_bharathi is offline   Reply With Quote

Old   November 29, 2017, 11:15
Default
  #9
Senior Member
 
Join Date: Jun 2009
Posts: 732
Rep Power: 15
Opaque is on a distinguished road
What is the reference temperature of the solid material definition?

Static Enthalpy = Integral from Tref to T of (Cp dT).

It is only

Static Enthalpy = Cp * (T - T_ref)

uniform Cp

Simplicity is plagued by the assumptions along the way. Assumptions are usually forgotten or implied and violated.
Opaque is offline   Reply With Quote

Old   November 29, 2017, 14:24
Default
  #10
Member
 
VB
Join Date: Jul 2016
Posts: 35
Rep Power: 3
Vishnu_bharathi is on a distinguished road
I formula volumeInt(density*Static Enthalpy)@domain worked. I am able to arrive the Energy in Joules. I have some other problem will post as new thread.
Vishnu_bharathi is offline   Reply With Quote

Reply

Tags
energy function, integral computing, porous domain, specific heats, temperature dependent

Thread Tools
Display Modes

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 Off
Trackbacks are On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
JANAF Coefficient of Air for Specific Heat Calculation. rhoCentralFoam cyril2304 OpenFOAM Running, Solving & CFD 1 October 8, 2017 07:44
Calculation of the Governing Equations Mihail CFX 7 September 7, 2014 06:27
poor calculation with temparture-dependant specific heat ryanmoser FLUENT 3 January 3, 2013 02:32
CFX doesn't continue calculation... mactech001 CFX 6 November 15, 2009 22:25
air bubble is disappear increasing time using vof xujjun CFX 9 June 9, 2009 07:59


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