CFD Online Logo CFD Online URL
Home > Forums > OpenFOAM Running, Solving & CFD

Question about fvmlaplacian operator

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

LinkBack Thread Tools Display Modes
Old   February 1, 2007, 14:40
Default Hi all, I'm trying to impleme
cosimo bianchini
Join Date: Mar 2009
Location: Florence, Tuscany, Italy
Posts: 88
Rep Power: 8
cosimobianchini is on a distinguished road
Send a message via Skype™ to cosimobianchini
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)

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


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 Bianchini

Ergon Research s.r.l.
Via Panciatichi, 92
50127 Florence - ITALY
Tel: +39 055 0763716
Mob: +39 320 9460153
cosimobianchini is offline   Reply With Quote


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
Question about the fvmatrix and Laplacian operator liuhuafei OpenFOAM Running, Solving & CFD 6 October 3, 2009 06:58
RunTimevalue and runTimetimeName with %3d%3d operator maka OpenFOAM 1 May 1, 2006 19:22
Operator Splitting. Maria. Main CFD Forum 5 September 17, 2005 22:10
Operator precedence hemph OpenFOAM 1 September 13, 2005 12:40
a math operator in UDF lichun Dong FLUENT 7 June 18, 2005 22:04

All times are GMT -4. The time now is 18:35.