CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Running, Solving & CFD

heat capacity in TEqn

Register Blogs Community New Posts Updated Threads Search

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   December 1, 2022, 06:11
Default heat capacity in TEqn
  #1
New Member
 
Join Date: Sep 2022
Posts: 18
Rep Power: 4
newGuyAtCFD is on a distinguished road
Hi,

Can anyone provide some insight on this:

VoF solvers weigh some properties as:

\rho = \rho_1\alpha_1 + \rho_2\alpha_2
\mu = \mu_1\alpha_1 + \mu_2\alpha_2

Whislt the heat capacity appears to be different.

\left(\frac{\alpha_1}{Cv_1} + \frac{\alpha_2}{Cv_2}\right) \left( \nabla \cdot  \left(\vec{v} p \right) + \frac{\partial \left(\rho K \right)}{\partial t}  +  \nabla \cdot \left(\rho \vec{v} K \right) - K \left[ \frac{\partial  \rho}{\partial t} + \nabla \cdot \left(\rho \vec{v}\right)\right] \right)

Why is C_v not being weighted as :
C_v = C_{v_1} \alpha_1 + C_{v_2} \alpha_2

And then:
\frac{1}{C_v}\left(  \nabla \cdot  \left(\vec{v} p \right) + \frac{\partial \left(\rho K  \right)}{\partial t}  +  \nabla \cdot \left(\rho \vec{v} K \right) - K  \left[ \frac{\partial  \rho}{\partial t} + \nabla \cdot \left(\rho  \vec{v}\right)\right] \right)


In the laplacian term \alpha_{eff} appears to be treated as :
\alpha_{eff} = \alpha_{eff_1} \alpha_1 + \alpha_{eff_2} \alpha_2


What is the rationale here?

I need to weigh the inverse of C_v since that is how the formulation is derived?

For instance, If I removed density from inside the parenthesis and diveded every term in the equation with \rho, I would need to use the same treatment as C_v?
newGuyAtCFD is offline   Reply With Quote

Reply


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


Similar Threads
Thread Thread Starter Forum Replies Last Post
Continuity Equation for multicomponent simulation lordluan CFX 15 May 19, 2020 19:36
heat capacity function Red Ember CFX 1 February 9, 2016 06:36
Heat capacity of species atoof OpenFOAM Programming & Development 4 September 15, 2014 13:31
Error finding variable "THERMX" sunilpatil CFX 8 April 26, 2013 08:00
Enforce bounds error with heat loss boundary condition at solid walls Chander CFX 2 May 1, 2012 21:11


All times are GMT -4. The time now is 01:00.