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

Missing gravity in new Teqn in compressibleInterFoam, OpenFOAM v2.2.0

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

Like Tree2Likes
  • 1 Post By dduque
  • 1 Post By dduque

Reply
 
LinkBack Thread Tools Display Modes
Old   March 7, 2013, 09:29
Default Missing gravity in new Teqn in compressibleInterFoam, OpenFOAM v2.2.0
  #1
New Member
 
Daniel Duque
Join Date: Jan 2011
Location: ETSIN, Madrid
Posts: 11
Rep Power: 6
dduque is on a distinguished road
Hi everybody.

We were greatly pleased by the modifications in compressibleInterFoam introduced in the new OpenFOAM v2.2.0 release. We were precisely working in extending the EOS to include non-isoT simulations, and the task is now basically done for us.

However, I notice the equation for temperatura in the new Teqn.H reads:

Code:
       fvm::ddt(rho, T)
      + fvm::div(rhoPhi, T)
      - fvm::laplacian(kByCv, T)
      + p*fvc::div(phi)*(alpha1/Cv1 + alpha2/Cv2)
Don't you have a term

rho gu / c_v

on the right hand side of this eq?

This is the power production due to gravity (well, rho gu is the power, divided by c_v gives a temperature variation).

Thanks. Best,

Daniel
dduque is offline   Reply With Quote

Old   March 11, 2013, 10:42
Default
  #2
New Member
 
Daniel Duque
Join Date: Jan 2011
Location: ETSIN, Madrid
Posts: 11
Rep Power: 6
dduque is on a distinguished road
No answers? Well, I'll try to implement it myself, while listening to this tune: Gravity's Gone.
JR22 likes this.
dduque is offline   Reply With Quote

Old   March 12, 2013, 07:59
Default
  #3
New Member
 
Daniel Duque
Join Date: Jan 2011
Location: ETSIN, Madrid
Posts: 11
Rep Power: 6
dduque is on a distinguished road
OK, I have just modified Teqn.H to solve this equation, including the
gravity power term, and it seems to be working.

Code:
    solve
    (
        fvm::ddt(rho, T)
      + fvm::div(rhoPhi, T)
      - fvm::laplacian(kByCv, T)
      + p*fvc::div(phi)*(alpha1/Cv1 + alpha2/Cv2)
        - rho*(U & g) *(alpha1/Cv1 + alpha2/Cv2)
    );
Best.
dduque is offline   Reply With Quote

Old   March 12, 2013, 17:26
Default
  #4
Senior Member
 
Markus Rehm
Join Date: Mar 2009
Location: Erlangen (Germany)
Posts: 176
Rep Power: 8
markusrehm is on a distinguished road
Hi Daniel,

consider posting a bug report to the Mantis:

http://www.openfoam.org/mantisbt

I am not sure if the developers follow the forum.

Markus
markusrehm is online now   Reply With Quote

Old   March 20, 2013, 14:50
Default
  #5
New Member
 
Daniel Duque
Join Date: Jan 2011
Location: ETSIN, Madrid
Posts: 11
Rep Power: 6
dduque is on a distinguished road
Actually, I was wrong: there is NO gravity term in the energy equation. At least, in this energy equation. There is one in the conservation version of the energy equation (for E=rho e + rho u^2/2), but dissapears in this version. Sorry!
alundilong likes this.
dduque is offline   Reply With Quote

Old   June 20, 2015, 23:37
Default
  #6
Member
 
yijin Mao
Join Date: May 2010
Location: Columbia, MO
Posts: 45
Rep Power: 7
alundilong is on a distinguished road
Quote:
Originally Posted by dduque View Post
Actually, I was wrong: there is NO gravity term in the energy equation. At least, in this energy equation. There is one in the conservation version of the energy equation (for E=rho e + rho u^2/2), but dissapears in this version. Sorry!
Yes. However, I did not see any gravity term, rho g&U, in the "total energy form".
alundilong is offline   Reply With Quote

Reply

Tags
compressibleinterfoam, gravity, openfoam v2.2.0, temperature, teqn

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
OpenFOAM 1.6.x, 1.7.0 and 1.7.x are not fully prepared to work with gcc-4.5.x wyldckat OpenFOAM Bugs 18 October 21, 2010 05:51
Cross-compiling OpenFOAM 1.7.0 on Linux for Windows 32 and 64bits with Mingw-w64 wyldckat OpenFOAM Announcements from Other Sources 3 September 8, 2010 06:25
OpenFoam gravity bug? Whyman OpenFOAM 7 June 21, 2010 04:10
Critical errors during OpenFoam installation in OpenSuse 11.0 amscosta OpenFOAM 5 May 1, 2009 14:06
Adventure of fisrst openfoam installation on Ubuntu 710 jussi OpenFOAM Installation 0 April 24, 2008 14:25


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