# FvmlaplaciansurfaceTensorField vf

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

 August 30, 2005, 07:17 Hi all, I expected laplacia #1 Member   diablo80@web.de Join Date: Mar 2009 Posts: 93 Rep Power: 9 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.

 August 30, 2005, 08:37 As you correctly say, if you w #2 Senior Member   Hrvoje Jasak Join Date: Mar 2009 Location: London, England Posts: 1,783 Rep Power: 22 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. Enjoy, Hrv brucechen likes this. __________________ Hrvoje Jasak Providing commercial FOAM/OpenFOAM and CFD Consulting: http://wikki.co.uk

 August 31, 2005, 08:10 Thanks Prof. Jasak, Yes I am #3 Member   diablo80@web.de Join Date: Mar 2009 Posts: 93 Rep Power: 9 Thanks Prof. Jasak, Yes I am using a symmetric tensor. cheers, luiz

 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

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