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

tr() of a volSymmTensorField ??

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

Reply
 
LinkBack Thread Tools Display Modes
Old   January 10, 2013, 19:19
Default tr() of a volSymmTensorField ??
  #1
Member
 
Jack
Join Date: Aug 2012
Posts: 42
Rep Power: 4
JackW is on a distinguished road
Hi foamers,

I have modified the vorticity.C file to create an application to calculate the strain rate tensor:


volSymmTensorField strainRate
(
IOobject
(
"strainRate",
runTime.timeName(),
mesh,
IOobject::NO_READ
),
symm(fvc::grad(U))
);

This compiled fine and ran. However as a sanity check I calculated:
max(tr(strainRate)).value()
min(tr(strainRate)).value()

...hoping it would be close to zero, (i am running an incompressible case)
and it wasn't, for the pisoFoam cavity case ~11, ~-39 were recorded.

So I ran,
max(fvc::div(U)).value()
min(fvc::div(U)).value()

To expect exactly the same value, however this gave ~7 and ~7.

I think I have misinterpreted some behaviour somewhere.

Any help would be greatly appreciated!

Thanks

Jack
JackW 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
RSTM laplacian volSymmTensorField Thomas Baumann OpenFOAM 0 April 11, 2012 05:35
EigenVectors and volSymmTensorField shane5ul OpenFOAM 0 November 25, 2009 12:23
Transformation of volSymmTensorField to volTensorField svens OpenFOAM 2 November 19, 2009 06:47
Problem with volSymmTensorField components cedric_duprat OpenFOAM Running, Solving & CFD 4 August 20, 2009 05:35
msha solve Eqn programming in C using FOAM msha OpenFOAM 2 December 30, 2007 03:39


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