CFD Online URL
[Sponsors]
Home > Forums > OpenFOAM

viscoelastic fluid abbreviations

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

Like Tree1Likes
  • 1 Post By JonW

Reply
 
LinkBack Thread Tools Display Modes
Old   October 4, 2012, 03:43
Post viscoelastic fluid abbreviations
  #1
Senior Member
 
Tushar C.
Join Date: Jul 2009
Location: India
Posts: 281
Blog Entries: 1
Rep Power: 8
Tushar@cfd is on a distinguished road
Hello Everybody,

I am stuck with the terms:

"// Twice the rate of deformation tensor
volSymmTensorField twoD = twoSymm(L);

// Stress transport equation
fvSymmTensorMatrix tauEqn
(
fvm::ddt(tau_)
+ fvm::div(phi(), tau_)
==
etaP_/lambda_*twoD
+ twoSymm(C)
- fvm::Sp(1/lambda_, tau_)
);

- I do know the term, it is: rate of deformation of stress and is expressed as: 1/2 * (grad(U) + grad(U).T).

-Does any one know from where, I will get the OpenFOAM coding of this term "twoSymm(L)" in of OF directory?

-I want to see the code definition of it and from which file it is compiled in the OF.

----
Thanks in Advance.
Tushar@cfd is offline   Reply With Quote

Old   October 4, 2012, 08:21
Default general viscoelastic material model - OpenFOAM-1.6-ext
  #2
Member
 
Jon Elvar Wallevik
Join Date: Nov 2010
Location: Reykjavik, ICELAND
Posts: 52
Rep Power: 6
JonW is on a distinguished road
dear Tushar Chourushi, it appears that you want to do a general viscoleastic simulation. This is a very difficult subject, both theoretically as well as numerically and there are lot of instability issues surrounding it. My recommendation is use rather OpenFOAM-1.6-ext. There is a general viscoelastic solver there.
It is more difficult to install the OF-ext version relative to the standard OF version, so my second recommendation is use ubuntu 12.04. I just recently installed it on one of mine Ubuntu server 12.04.1 (64bit) machines without problems. My third recommendation is, use lot of time to prepare, both for the installation of OF-ext as well in learning on the viscoelastic solver. Just jumping into this and expecting things to go you way is most certainly a killer. Be patient in your difficult work ahead.
good luck,
Jon
JonW is offline   Reply With Quote

Old   October 4, 2012, 08:33
Default Dear Jon
  #3
Senior Member
 
Tushar C.
Join Date: Jul 2009
Location: India
Posts: 281
Blog Entries: 1
Rep Power: 8
Tushar@cfd is on a distinguished road
Thanks for your suggestion, I have already installed the viscoelasticFoam on my machine. I wanted to see how the equations are written in openFOAM as compared to the constitutive equation for different models?
That too.. I was having problem with the term which I have mentioned in my earlier thread, as I was not able to co-relate it with the theory.

Anyways, Thanks for your help.

---
Regards,
Tushar
Tushar@cfd is offline   Reply With Quote

Old   October 4, 2012, 09:27
Default about twoSymm(L)
  #4
Member
 
Jon Elvar Wallevik
Join Date: Nov 2010
Location: Reykjavik, ICELAND
Posts: 52
Rep Power: 6
JonW is on a distinguished road
hi again,
I dont think there is any twoSymm(L) used in the code. But if you look at the

~/OpenFOAM/OpenFOAM-2.1.x/src/transportModels/incompressible/viscosityModels/viscosityModel/viscosityModel.C
there is a component,...
sqrt(2.0)*mag(symm(fvc::grad(U_)));
..., which calculates the shear rate.
thus the rate of deformation tensor 1/2 * (grad(U) + grad(U).T) is "symm(fvc::grad(U_))"
The mag(symm(fvc::grad(U_))); would be its double dot, square root.

Any you can see a bit by searching for "shear rate"
J.
JonW is offline   Reply With Quote

Old   January 5, 2013, 00:41
Smile
  #5
Senior Member
 
Tushar C.
Join Date: Jul 2009
Location: India
Posts: 281
Blog Entries: 1
Rep Power: 8
Tushar@cfd is on a distinguished road
Thanks Jon,

Actually I wasn't talking about the available code with OpenFOAM, I was looking at the extended source code of OpenFOAM. There I found that variable, anyways thanks for your help.
I got the answer while exploring the forum.
Tushar@cfd is offline   Reply With Quote

Old   February 25, 2014, 03:30
Default shear rate
  #6
Member
 
s.rasoul_varedi
Join Date: Feb 2010
Posts: 82
Rep Power: 5
desert_1250 is an unknown quantity at this point
Send a message via Yahoo to desert_1250
Hi foamers
How can i extract shear rate from my results using viscoelasticFluidFoam?
tnx
desert_1250 is offline   Reply With Quote

Old   February 25, 2014, 20:38
Default shear rate
  #7
Member
 
Jon Elvar Wallevik
Join Date: Nov 2010
Location: Reykjavik, ICELAND
Posts: 52
Rep Power: 6
JonW is on a distinguished road
Quote:
Originally Posted by desert_1250 View Post
Hi foamers
How can i extract shear rate from my results using viscoelasticFluidFoam?
tnx
Hi there, I am not sure what you are looking for here. But if you want to write the shear rate down for each time step, then:

(1) in "createFields.H" add the following:
------------------------------------------------
volScalarField shearRate
(
IOobject
(
"shearRate",
runTime.timeName(),
mesh,
IOobject::NO_READ,
IOobject::AUTO_WRITE
),
mesh,
dimensionedScalar("shearRate", dimensionSet(0,0,-1,0,0,0,0), scalar(0.0))
);
------------------------------------------------



(2) in viscoelasticFluidFoam.C

before "runTime.write();" add the following
------------------------------------------------
dimensionedScalar sqrtOfTwo("sqrtOfTwo", dimless, 1.4142135623730950488);
shearRate = sqrtOfTwo*mag(symm(fvc::grad(U)));
------------------------------------------------

The runTime.write(); thingi will write the shear rate automatically.


One note, now I am assuming that you accept that the shear rate is "sqrtOfTwo*mag(symm(fvc::grad(U)));"
It is so for incompressible fluid that comply to the Generalized Newtonian Model (GNM), e.g. Bingham viscoplastic model, power law fluid or Newtonian fluid. But viscoelasticity is not GNM. If you are in trouble, then you can just call this
"Magnitude of the rate-of-strain tensor" like I did in
Wallevik, J.E., Minimizing end-effects in the coaxial cylinders viscometer: Viscoplastic flow inside the ConTec BML Viscometer 3. J. Non-Newtonian Fluid Mech. 155 (2008) 116-123.

Hope this helps
cheers
Jon
desert_1250 likes this.
JonW is offline   Reply With Quote

Old   February 26, 2014, 01:59
Default
  #8
Member
 
s.rasoul_varedi
Join Date: Feb 2010
Posts: 82
Rep Power: 5
desert_1250 is an unknown quantity at this point
Send a message via Yahoo to desert_1250
Quote:
Originally Posted by JonW View Post
Hi there, I am not sure what you are looking for here. But if you want to write the shear rate down for each time step, then:

(1) in "createFields.H" add the following:
------------------------------------------------
volScalarField shearRate
(
IOobject
(
"shearRate",
runTime.timeName(),
mesh,
IOobject::NO_READ,
IOobject::AUTO_WRITE
),
mesh,
dimensionedScalar("shearRate", dimensionSet(0,0,-1,0,0,0,0), scalar(0.0))
);
------------------------------------------------



(2) in viscoelasticFluidFoam.C

before "runTime.write();" add the following
------------------------------------------------
dimensionedScalar sqrtOfTwo("sqrtOfTwo", dimless, 1.4142135623730950488);
shearRate = sqrtOfTwo*mag(symm(fvc::grad(U)));
------------------------------------------------

The runTime.write(); thingi will write the shear rate automatically.


One note, now I am assuming that you accept that the shear rate is "sqrtOfTwo*mag(symm(fvc::grad(U)));"
It is so for incompressible fluid that comply to the Generalized Newtonian Model (GNM), e.g. Bingham viscoplastic model, power law fluid or Newtonian fluid. But viscoelasticity is not GNM. If you are in trouble, then you can just call this
"Magnitude of the rate-of-strain tensor" like I did in
Wallevik, J.E., Minimizing end-effects in the coaxial cylinders viscometer: Viscoplastic flow inside the ConTec BML Viscometer 3. J. Non-Newtonian Fluid Mech. 155 (2008) 116-123.

Hope this helps
cheers
Jon

Dear Jan,
I exerted and compiled the code succesefully, now i have shear rate at every time writing.

thankful for your quick and valuable reply.
Kind Regards/ Rasoul
desert_1250 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
Fluid pairs fjalil Main CFD Forum 0 June 10, 2009 14:47
Fluid Structure Interaction Apollo FLUENT 8 July 7, 2004 18:54
My Revised "Time Vs Energy" Article For Review Abhi Main CFD Forum 2 July 9, 2002 10:08
Terrible Mistake In Fluid Dynamics History Abhi Main CFD Forum 12 July 8, 2002 10:11
fluid hot volume in fluid cold volume zahid FLUENT 4 June 1, 2002 10:11


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