# Question about fvmlaplacian operator

 User Name Remember Me Password
 Register Blogs Members List Search Today's Posts Mark Forums Read

 LinkBack Thread Tools Display Modes
 February 1, 2007, 14:40 Hi all, I'm trying to impleme #1 Member   cosimo bianchini Join Date: Mar 2009 Location: Florence, Tuscany, Italy Posts: 88 Rep Power: 10 Hi all, I'm trying to implement an anisotropic turbulence model. What I need to do is to correct turbulence->divRhoR() in such a way: -- read Uv = velocity vector = (U,V,W) -- divRhoR = lam + L + Dev laminarPart lam = fvm::laplacian(mu(), Uv) fvc::div(mu()*dev2(fvc::grad(Uv)().T())) laplacianPart L = mut* [(Uxx + gamma*Uyy + Uzz), (gamma*Vxx + Vyy + Vzz), (Wxx + Wyy +Wzz)] deviatoricPart Dev = 1/3*mut*[ (Uxx + (3*gamma - 2)*Vxy + Wxz), ((3*gamma-2)*Uyx + Vyy + Wyz), (Uzx + Vzy + Wzz)] What I have done is: volTensorField B = mut()*dev2(fvc::grad(U)().T()); B.replace(tensor::XY,(B.component(tensor::XY) * gamma)); B.replace(tensor::YX,(B.component(tensor::YX) * gamma)); Dev = fvc::div(B) Right now I'm stuck with the implementation of the laplacian part because I would like to make it implicit but I don't know what (vector,tensor,..) to put in place of A fvm::laplacian(A,U) Now I get to the questions: 1) How is actually fvm::laplacian Operator working? 2) Is it doing something similar to div(A,grad(U)) ? 3) What kind of product does it use in multiplying A and U? 4) Is div(T) a = (T11x + T12x + T13x, T21y + T22y + T23y, T31z + T32z + T33z) [as written in programmer's guide] or b = (T11x + T12y + T13z, T21x + T22y + T23z, T31x + T32y + T33z) ? Thanks really a lot in advance for your answers Cosimo __________________ Cosimo Bianchini Ergon Research s.r.l. Via Panciatichi, 92 50127 Florence - ITALY Tel: +39 055 0763716 Mob: +39 320 9460153 e-mail: cosimo.bianchini@ergonresearch.it URL: www.ergonresearch.it

 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 liuhuafei OpenFOAM Running, Solving & CFD 6 October 3, 2009 06:58 maka OpenFOAM 1 May 1, 2006 19:22 Maria. Main CFD Forum 5 September 17, 2005 22:10 hemph OpenFOAM 1 September 13, 2005 12:40 lichun Dong FLUENT 7 June 18, 2005 22:04

All times are GMT -4. The time now is 08:28.

 Contact Us - CFD Online - Privacy Statement - Top