CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Running, Solving & CFD

Question about fvmlaplacian operator

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

Reply
 
LinkBack Thread Tools Display Modes
Old   February 1, 2007, 14:40
Default Hi all, I'm trying to impleme
  #1
Member
 
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)
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

Energy Engineering Department "S. Stecco"
University of Florence
Via di S.Marta, 3
50139 Florence - ITALY
Tel: +39 055 4796575
Fax: +39 055 4796342
Mob: +39 320 9460153
e-mail: cosimo.bianchini@htc.de.unifi.it
URL: www.htc.de.unifi.it
cosimobianchini is offline   Reply With Quote

Reply

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 10:02.