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

FvmlaplaciansurfaceTensorField vf

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

Like Tree1Likes
  • 1 Post By hjasak

LinkBack Thread Tools Display Modes
Old   August 30, 2005, 07:17
Default Hi all, I expected laplacia
Join Date: Mar 2009
Posts: 93
Rep Power: 10
sampaio is on a distinguished road
Hi all,

I expected laplacian scheme when using a tensor Gamma to reproduce, say, an anisotropic material heat transfer.

But when I checked gaussLaplacianScheme.C, I realized that it actually takes the trace of the tensor, instead of the tensor itself...

What is the reason for that? (I suspect it is because otherwise we could not multiply Sf by vf and take advantage of snGrad. Is that right?)

Any suggestions on how to implement a fvm scheme like div ( gammaTensor * grad(T) ) or even div (gammaTensor & grad(U))?

Thanks a lot.
sampaio is offline   Reply With Quote

Old   August 30, 2005, 08:37
Default As you correctly say, if you w
Senior Member
Hrvoje Jasak
Join Date: Mar 2009
Location: London, England
Posts: 1,810
Rep Power: 25
hjasak will become famous soon enough
As you correctly say, if you want to have tensorial viscosity and make use of the snGrad, you can only make it implicit in the trace. What you do is to decompose the tensor into the diagonal bit and the rest and for the diagonal bit the (S_f . grad) maps through to

(phi_N - phi_P)/distance

For the rest, you cannot do the trick because the tensorial viscosity dots the gradient before the face area vector (so there is a rotation). Thus, for the full tensorial viscosity, you will have an implicit and explicit part; you have to be pretty careful because of possible stability/boundedness problems. For an example of how to add the explicit part, have a look at how non-orthogonal correction is handled - pretty easy.

If your variable is a vector, we can recognise that tensorial viscosity actually represents a rotation and couples the components of the vector, making it interesting for the block solver, but that's further down the line.

BTWm can you tell me what kind of tensorial viscosity you are looking at - I suspect it is at least a symmetric tensor. If you align the grid with the eigenvalues (like you can for porous media), you may get much better behaviour because the explicit correction vanishes.


brucechen likes this.
Hrvoje Jasak
Providing commercial FOAM/OpenFOAM and CFD Consulting:
hjasak is offline   Reply With Quote

Old   August 31, 2005, 08:10
Default Thanks Prof. Jasak, Yes I am
Join Date: Mar 2009
Posts: 93
Rep Power: 10
sampaio is on a distinguished road
Thanks Prof. Jasak,
Yes I am using a symmetric tensor.
sampaio 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

All times are GMT -4. The time now is 00:58.