|
[Sponsors] |
curvature correction term, material derivative of a tensor |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
April 1, 2010, 11:08 |
curvature correction term, material derivative of a tensor
|
#1 |
New Member
Volker Tritschler
Join Date: Jan 2010
Posts: 20
Rep Power: 16 |
Hi,
I'm planning to implement a curvature correction term to an explicit algebraic reynolds stress model based on the existing k-e-model of Launder and Sharma in order to sensitize the Launder-Sharma-model to streamline curvature. The curvature correction term is to be implemented according to P.E. Smirnov and F.R. Menter "Sensitization of the SST Turbulence Model to Rotation and Curvature by Applying the Spalart-Shur Correction Term". I ran into problems quite fast. There is a term, which is the material derivative of the shear strain tensor. The implementation is not straight-forward, especially in tensor notation. Either one has somehow to implement a summing over all faces of a control volume or to evaluate the gradient of the shear strain tensor, which leads to a tensor of third order. And right now, I do not see a way to get rid of this problem. Does anyone has already tried to implement a curvature correction or has any experience on that? Or even know how to implement a material derivative of a tensor? I'm glad for any hint. Greets, volker |
|
April 1, 2010, 11:46 |
|
#2 |
Senior Member
|
I have try RC correction as you mentioned in ref paper. you can send me a mail. if interesting.
my email address: waynezw0618@163.com |
|
May 6, 2010, 16:55 |
Gradient of a second order tensor
|
#3 |
New Member
Michael B Martell Jr
Join Date: Feb 2010
Location: Amherst, MA
Posts: 18
Rep Power: 16 |
I am attempting to add a term to a Reynolds stress transport model which involves taking the gradient of the Reynolds stress tensor, which should result in a tensor of rank three.
Even a simple attempt fails: Code:
volTensorField test = fvc::grad(R); Code:
error: no type named ‘type’ in ‘class Foam::typeOfRank<double, 3>’ error: no matching function for call to 'grad(Foam::GeometricField<Foam::Tensor<double>, Foam::fvPatchField, Foam::volMesh>&)’ Any ideas or help would be appreciated! |
|
May 6, 2010, 17:55 |
|
#4 |
New Member
Volker Tritschler
Join Date: Jan 2010
Posts: 20
Rep Power: 16 |
hi michael,
I had exactly the same problem, though I found it kind of surprising because OpenFOAM is supposed to handle tensors of rank three. But the gradient seems to be restricted to tensors of rank one. I hope that there might be a workaround. hope someone can help us. cheers, volker |
|
May 6, 2010, 18:03 |
|
#5 |
New Member
Michael B Martell Jr
Join Date: Feb 2010
Location: Amherst, MA
Posts: 18
Rep Power: 16 |
Hi Volker. I have discussed this problem with some of my colleagues here and we agree that this would be non-trivial to implement. I am going to try, and I'll keep you posted on my progress. Feel free to contact me through http://www.ecs.umass.edu/~mmartell if you'd like.
FYI, I have tried not assigning the gradient to a particular type by simple putting Code:
fvc::grad(R); Cheers, Mike |
|
May 7, 2010, 08:29 |
|
#6 |
New Member
Volker Tritschler
Join Date: Jan 2010
Posts: 20
Rep Power: 16 |
would be great, if you would come up with something new.
keep me posted! |
|
May 27, 2010, 04:50 |
|
#7 |
New Member
Volker Tritschler
Join Date: Jan 2010
Posts: 20
Rep Power: 16 |
hi michael,
did you manage to implement the gradient of a tensor? cheers, volker |
|
June 3, 2010, 09:08 |
|
#8 |
New Member
Michael B Martell Jr
Join Date: Feb 2010
Location: Amherst, MA
Posts: 18
Rep Power: 16 |
Hi Volker,
I managed to do what I needed "by hand" but it is sloppy and not generic. I'll attach the source for it. It has no external dependencies (you can ignore the inclusion of OEC.H) and it doesn't need to be cast in the incompressible or RASModels name space - I simply took it right from my existing code. In my case, I needed something like fvc::grad(kINT) & fvc::grad(RINT) with kINT being a volScalarField and RINT being a volSymmTensorField. Note that this inner product (of grad(kINT), a vector, and grad(RINT), a rank 3 tensor) results in a rank 2 symmetric tensor. I know this isn't exactly what you were looking for, but the structure is there. Let me know if you need any help Mike |
|
Thread Tools | Search this Thread |
Display Modes | |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
mass flow in is not equal to mass flow out | saii | CFX | 12 | March 19, 2018 06:21 |
Curvature Correction & Reattachment Mod. | Felix | CFX | 5 | December 19, 2010 17:04 |
How to use SST Curvature Correction | mike wilson | CFX | 11 | December 19, 2010 17:04 |
ATTENTION! Reliability problems in CFX 5.7 | Joseph | CFX | 14 | April 20, 2010 16:45 |
Two-Phase Buoyant Flow Issue | Miguel Baritto | CFX | 4 | August 31, 2006 13:02 |