
[Sponsors] 
June 30, 2011, 13:22 
one question about representation of equations?

#1 
Member

Hi, I am new in openFoam, and I try to understand this equation from XiFoam solver:
{ fvScalarMatrix hEqn ( fvm::ddt(rho, h) + mvConvection>fvmDiv(phi, h)  fvm::laplacian(turbulence>alphaEff(), h) == DpDt ); hEqn.relax(); hEqn.solve(); thermo.correct(); } This is the enthalpy equation, I want to understand what is the meaning of each term, I understand that fvm::ddt(rho, h) is the representation of the time derivative, fvmDiv(phy,h) represent the convective term, but in this stage I am some confused about what is the meaning of the mvConvection>here, the next term is the diffusive term (fvm::laplacian(turbulence>alphaEff(), h)) here I am confused with the meaning of turbulence>alphaEff(). I will be a lot thankful if someone can explain me this. thanks for your help Best Regards
__________________
Felipe G 

June 30, 2011, 17:51 

#2 
Senior Member
David Gaden
Join Date: Apr 2009
Location: Winnipeg, Canada
Posts: 436
Rep Power: 14 
You have to remember, this is still C++ code. They've overloaded the operators and created namespaces in an elegant way so it looks like you are typing in an equation, but it is still just C++, and sometimes the syntax gets in the way a little.
Code:
+ mvConvection>fvmDiv(phi, h) Code:
+ fvm::div(phi, h) Another reason for this may be because there are a lot of variables using a common convectionScheme... so rather than clutter up your fvSchemes file with entries for each variable, they define a common one, and only look that one up. Code:
turbulence>alphaEff() Hope that helps. 

June 30, 2011, 18:52 

#3 
Member

Thanks for the reply David, was very helpful for me. Another question, the last term , DpDt is like an assignment ?
Best regards
__________________
Felipe G 

July 1, 2011, 18:57 

#4  
Senior Member
Daniel P. Combest
Join Date: Mar 2009
Location: St. Louis, USA
Posts: 603
Rep Power: 21 
Quote:
Dan 

July 1, 2011, 23:49 

#5 
Member

thanks for your reply Dan
__________________
Felipe G 

July 2, 2011, 12:02 

#6 
Senior Member
David Gaden
Join Date: Apr 2009
Location: Winnipeg, Canada
Posts: 436
Rep Power: 14 
Thanks Dan. Your article looks like a good read... relevant to my stuff, too.
Felipe, the last term DpDt isn't an assignment. It's the "equal to" operator. It behaves as an equals sign would in the equation you are creating. 

August 18, 2011, 07:21 

#7  
Member
A. Bernath
Join Date: Jun 2011
Location: Karlsruhe, Germany
Posts: 39
Rep Power: 7 
Quote:
Quote:
Thx, Alex EDIT: Just figured it out: DpDt is the total differentiation Dp = dp/dt + dp/dx Dx + dp/dy Dy + dp/dz Dz divided by Dt. With Dx/Dt = u, Dy/Dt = v, ... the above equation follows. My Question now is: Where is the implementation of this operator? I had no success by searching in the User Guide and neither was there a found in the Programming Guide. I also found a table with most of the operators (e.g. ddt(p)) but no word about DpDt. EDITē: Found the implementation in createFields.H. The "find" command under linux seems to be very helpful in such cases =) Last edited by derkermit; August 18, 2011 at 11:22. 

Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
coupled source term in 3 transport equations  mhassani  OpenFOAM Running, Solving & CFD  0  February 26, 2011 05:14 
2 Different Transonic Small Disturbance Equations?  Frank  Main CFD Forum  1  February 23, 2006 00:44 
Solver Equations  Carola  CFX  9  August 12, 2003 08:27 
Pressureless EulerPoisson Equations  Mikael Tushensson  Main CFD Forum  0  April 12, 2001 20:18 
Pictures and Equations?  andy  Main CFD Forum  3  December 11, 1998 10:51 