# Viscous term expansion

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

 March 23, 2011, 08:35 Viscous term expansion #1 Senior Member     Santiago Marquez Damian Join Date: Aug 2009 Location: Santa Fe, Santa Fe, Argentina Posts: 430 Rep Power: 16 Sponsored Links Hi all, I'm studying the implementation of the viscous term in NS equations in the general framework used for example in pisoFoam, nevertheless I can't understand some of the things that are done. Viscous term is: div[nuEff*(grad(U)+grad(U)^T)] due incompressibility and applying div to the product, it becomes laplacian(nuEff, U)+ grad(nuEff)*[(grad(U)+grad(U)^T)] (1) in laminar regime nuEff=nu and in turbulent regime nuEff=nu+nut. Reading the code we have from pisoFoam.C: Code: 00065 fvVectorMatrix UEqn 00066 ( 00067 fvm::ddt(U) 00068 + fvm::div(phi, U) 00069 + turbulence->divDevReff(U) 00070 ); so that the viscous term is contained in turbulence->divDevReff(U) calling. Taking, for example, the kEpsilon.C implementation: Code: 00181 tmp kEpsilon::divDevReff(volVectorField& U) const 00182 { 00183 return 00184 ( 00185 - fvm::laplacian(nuEff(), U) 00186 - fvc::div(nuEff()*dev(fvc::grad(U)().T())) 00187 ); 00188 } the laplacian term is present, but the second term is completely different to the second one in (1). I tried to derive the above formulation without success and my advisor showed me that both expressions [(1) and which is present in divDevRefff] give different results in deed, using a counterexample with a div(U)=0 field. All ideas are welcome. Regards. __________________ Santiago MÁRQUEZ DAMIÁN, Ph.D. Research Scientist Research Center for Computational Methods (CIMEC) - CONICET/UNL Tel: 54-342-4511594 Int. 7032 Colectora Ruta Nac. 168 / Paraje El Pozo (3000) Santa Fe - Argentina. http://www.cimec.org.ar

 December 13, 2011, 21:44 #2 New Member   Join Date: Mar 2009 Location: Sao Jose dos Campos, Brazil Posts: 28 Rep Power: 10 Hello, Santiago. Perhaps you have it already solved, but here goes my try: The laminar part of viscous stress tensor is: The turbulent part of viscous stress tensor is Now, summing both contributions: The last term in OpenFOAM code is not the same, though: - it is computed the deviatoric of and not the gradient itself. This is discussed in another topic: http://www.cfd-online.com/Forums/ope...ivdevreff.html The trace of last term is zero because of the explicit formulation using the velocity field from previous time step, that should be divergence free.

 Thread Tools Display Modes Linear Mode

 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 OffTrackbacks are On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post siw CFX 1 December 6, 2010 07:41 dandalf Main CFD Forum 7 November 17, 2010 10:21 seb62 OpenFOAM Running, Solving & CFD 0 March 19, 2009 04:41 mt Phoenics 10 February 11, 2005 11:46 yann Phoenics 0 June 5, 2003 10:23