CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Programming & Development

Question about EEqn.H file

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

Like Tree1Likes
  • 1 Post By alundilong

Reply
 
LinkBack Thread Tools Display Modes
Old   July 7, 2014, 23:56
Default Question about EEqn.H file
  #1
New Member
 
Harshad Lalit
Join Date: May 2013
Posts: 26
Rep Power: 5
harshad88 is on a distinguished road
Dear all,

I am trying to develop a new solver for premixed flame simulations using reactingFoam. I have a question about the EEqn.H file in the default code. It reads:

+ (
he.name() == "e"
? fvc::div
(
fvc::absolute(phi/fvc::interpolate(rho), U),
p,
"div(phiv,p)"
)
: -dpdt
)


What is the meaning of fvc::absolute(phi/fvc::interpolate(rho), U), and the line : -dpdt ???

Your help would be really appreciated. Thanks
harshad88 is offline   Reply With Quote

Old   July 8, 2014, 17:26
Default
  #2
Senior Member
 
Wouter van der Meer
Join Date: May 2009
Location: Elahuizen, Netherlands
Posts: 144
Rep Power: 9
wouter is on a distinguished road
Hello harshad88,

I am no expert.
It is part of a C or C++ construct A?b:c meaning
if A is true than b else c in this case c = "-dpdt"

hope this helps
Wouter
wouter is offline   Reply With Quote

Old   July 9, 2014, 11:10
Default
  #3
New Member
 
Harshad Lalit
Join Date: May 2013
Posts: 26
Rep Power: 5
harshad88 is on a distinguished road
Thanks for your reply Wouter.

Could you also tell me what the meaning of fvc::absolute(phi/fvc::interpolate(rho), U) is? I cant understand what this terms corresponds to in the energy equation.

Thanks
harshad88 is offline   Reply With Quote

Old   July 9, 2014, 17:34
Default
  #4
Senior Member
 
Wouter van der Meer
Join Date: May 2009
Location: Elahuizen, Netherlands
Posts: 144
Rep Power: 9
wouter is on a distinguished road
As i wrote I am no expert but I think this means phiv because of the "div(phiv,p)" but I have no clue what this means.

Hope someone else can take over
Wouter
wouter is offline   Reply With Quote

Old   July 10, 2014, 19:25
Default
  #5
Member
 
yijin Mao
Join Date: May 2010
Location: Columbia, MO
Posts: 48
Rep Power: 8
alundilong is on a distinguished road
Quote:
Originally Posted by harshad88 View Post
Thanks for your reply Wouter.

Could you also tell me what the meaning of fvc::absolute(phi/fvc::interpolate(rho), U) is? I cant understand what this terms corresponds to in the energy equation.

Thanks
fvc::absolute(phi/fvc::interpolate(rho), U) return absolute face flux.
if the mesh is moving, it will add up the contribution of moving mesh,
which mean the face flux is
phi/fvc::interpolate(rho) + fvc::meshPhi(U)
;
if static mesh, it will return phi/fvc::interpolate(rho)
fumiya likes this.
alundilong is offline   Reply With Quote

Old   July 10, 2014, 22:07
Default
  #6
Member
 
Daniel
Join Date: Jun 2014
Posts: 60
Rep Power: 4
Dan1788 is on a distinguished road
Thanks Yijin, That was really helpful
Dan1788 is offline   Reply With Quote

Reply

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
what is swap4foam ?? AB08 OpenFOAM 28 February 2, 2016 02:22
groovyBC in openFOAM-2.0 for parabolic velocity bc ofslcm OpenFOAM Pre-Processing 24 August 23, 2015 22:41
patch error OF v2.2.2 mac hewei OpenFOAM Installation on Windows, Mac and other Unsupported Platforms 4 November 30, 2013 17:55
Annoying issue of automatic "Rescale to Data Range " with paraFoam/paraview 3.12 keepfit OpenFOAM Paraview & paraFoam 60 September 18, 2013 03:23
Version 15 on Mac OS X gschaider OpenFOAM Installation 120 December 2, 2009 11:23


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