
[Sponsors] 
July 22, 2011, 11:15 
Laplacian term in NS equation

#1 
New Member
Diego Villa
Join Date: Mar 2010
Location: Genova Italy
Posts: 28
Rep Power: 8 
Hi All,
I have a little question on the NavierStokes implementation in OF. In the interPhaseChangeFoam solver the NS equation has written as: fvVectorMatrix UEqn ( fvm::ddt(rho, U) + fvm::div(rhoPhi, U)  fvm::Sp(fvc::ddt(rho) + fvc::div(rhoPhi), U)  fvm::laplacian(muEff, U)  (fvc::grad(U) & fvc::grad(muEff)) ); My question is: what does the term " (fvc::grad(U) & fvc::grad(muEff))" mean? In the theory that term shouldn't exist. if it's true the relationship laplacian(muEff, U)=muEff*laplacian(U)+grad(U)*grad(muEff) (see P38 in Programmers Guide) I make some mistake? Diego 

July 25, 2011, 23:26 

#2  
Senior Member
Mirko Vukovic
Join Date: Mar 2009
Posts: 159
Rep Power: 9 
Quote:
Mirko 

July 26, 2011, 05:38 

#3 
New Member
Diego Villa
Join Date: Mar 2010
Location: Genova Italy
Posts: 28
Rep Power: 8 
Hi Mirko
I'm not sure of that. Why this terms should be zero when the solution is reached? Both therms are not zero at the end, and why the inner product should be zero? Could you send me a reference about that? Thanks a lot for the replay Diego 

July 26, 2011, 06:25 

#4 
Senior Member
Pablo
Join Date: Mar 2009
Posts: 102
Rep Power: 9 
Hi Diego,
The answer is here "http://www.cfdonline.com/Forums/openfoam/82640interfoamueqn.html" and here "http://powerlab.fsb.hr/ped/kturbo/OpenFOAM/SummerSchool2009/presentations/MitjaMorgut2009.pdf" 

July 26, 2011, 11:45 

#5 
New Member
Diego Villa
Join Date: Mar 2010
Location: Genova Italy
Posts: 28
Rep Power: 8 
Hi Pablo
I'm agree that if the muEff is not constant the terms "div(muEff*grad(U))" and "muEff*lapalcian(U)" are different, but in the Navie Stokes equation the first term is the right one. But as the OF Programmers Guide tell at P38, the openfoam command laplacian(muEff,U), should be intended as the div(muEff*grad(U)) and not asmuEff*lapalcian(U), so should be not necessary split the terms in two contribute in the equation. If I say something wrong...could you say me where? Diego 

July 26, 2011, 14:36 

#6 
Senior Member
Pablo
Join Date: Mar 2009
Posts: 102
Rep Power: 9 

July 27, 2011, 03:44 

#7 
Senior Member
Anton Kidess
Join Date: May 2009
Location: Delft, Netherlands
Posts: 1,139
Rep Power: 20 
Pablo, you have errors in the bracketing of the equations in the three lines of your eqnarray. The final result is correct though


July 27, 2011, 04:25 

#8 
New Member
Diego Villa
Join Date: Mar 2010
Location: Genova Italy
Posts: 28
Rep Power: 8 
Thank you very much Pablo,
I found the little mistake in the brackets, but now all it is more clear. In every books i found that this therm is always hides in the source terms without explicit its. Thank you again. 

July 27, 2011, 06:56 

#9 
Senior Member
Pablo
Join Date: Mar 2009
Posts: 102
Rep Power: 9 
Sorry, brackets are wrong at the stress term.
Pablo 

July 20, 2012, 14:32 
transpose term

#10 
Senior Member
ehsan
Join Date: Mar 2009
Posts: 106
Rep Power: 9 
Dear All
I like to ask a question: What happen to the term of transpose of (Grad (U))? To what term it is converted in equation given in "http://imageshack.us/photo/myimages/3/formulasns.jpg/" Regards Ehsan 

Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
Calculation of the Governing Equations  Mihail  CFX  7  September 7, 2014 06:27 
Large source term in species equation  MACFD  FLUENT  4  January 4, 2011 15:16 
Is it possible to use divergence in the source term of an scalar transport equation?  jannnesss  CFX  0  January 8, 2010 20:53 
What is the Boussinesq Term in Momentum Equation  CFDtoy  Main CFD Forum  0  August 11, 2008 09:56 
bouyancy term in epsilon equation  Michael  Main CFD Forum  1  June 25, 1999 10:20 