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

Viscoelastic Fluid Flows using OpenFOAM The solver viscoelasticFluidFoam

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

Like Tree23Likes

Reply
 
LinkBack Thread Tools Display Modes
Old   November 2, 2013, 11:00
Default
  #321
Senior Member
 
Dongyue Li
Join Date: Jun 2012
Location: Torino, Italy
Posts: 661
Rep Power: 8
sharonyue is on a distinguished road
Thanks Mostafa,

BTW, do you know that if the surface tension is important in this viscoelasticInterFoam?
sharonyue is offline   Reply With Quote

Old   November 2, 2013, 14:46
Default
  #322
Senior Member
 
adambarfi's Avatar
 
Mostafa Mahmoudi
Join Date: Jan 2012
Posts: 314
Rep Power: 6
adambarfi is on a distinguished road
Send a message via Yahoo to adambarfi Send a message via Skype™ to adambarfi
Quote:
Originally Posted by sharonyue View Post
Thanks Mostafa,

BTW, do you know that if the surface tension is important in this viscoelasticInterFoam?
sorry, I didn't hear about viscoelasticInterFoam before!!!! but for multiphase flows the surface tension is necessary.
adambarfi is offline   Reply With Quote

Old   November 4, 2013, 10:03
Default
  #323
Senior Member
 
Dongyue Li
Join Date: Jun 2012
Location: Torino, Italy
Posts: 661
Rep Power: 8
sharonyue is on a distinguished road
Okay, back to the UEqn. I mean the code implanted in 1.6ext is:

Code:
fvc::div(tau_, "div(tau)")
      - fvc::laplacian(etaP(), U)
      + fvm::laplacian(etaP() + etaS(), U)
Why dont we just use:
Code:
fvc::div(tau_, "div(tau)")
      + fvm::laplacian(etaS(), U)
sharonyue is offline   Reply With Quote

Old   November 14, 2013, 04:59
Default viscoelasticProperties Problem
  #324
New Member
 
KimJuHyun
Join Date: Nov 2013
Posts: 2
Rep Power: 0
JuHyun is on a distinguished road
Hello Foamers,

i'm korean student, Want you to help me.

Perhaps it may be a very easy question.

i'm using the viscoelasticFluidFOAM multi-phase flow simulation tries.

but, this file. I did different from the content of the file.

'transportProperties' is the file I was using.

phase1
{
transportModel HerschelBulkley;
nu nu [ 0 2 -1 0 0 0 0 ] 3e-4;
rho rho [1 -3 0 0 0 0 0] 1000;

HerschelBulkleyCoeffs
{
nu0 nu0 [0 2 -1 0 0 0 0] 0.006;
tau0 tau0 [0 2 -2 0 0 0 0] 0.025031;
k k [0 2 -1 0 0 0 0] 3e-4;
n n [0 0 0 0 0 0 0] 1;
}
}
phase2
{
transportModel Newtonian;
nu nu [0 2 -1 0 0 0 0] 1.46e-05;
rho rho [1 -3 0 0 0 0 0] 1.225;
}

refPhase phase2;

sigma sigma [1 0 -2 0 0 0 0] 0.066;


but, the viscoelasticFluidFOAM in transportProperties without using viscoelasticProerties[Giesekus] is used.

I Have used a lot of different ways. but Always viscoelasticProperties in the error appeared.

viscoelasticProperties

rheology
{

type multiMode;

models
(
first
{
type Giesekus;
rho rho [1 -3 0 0 0 0 0] 803.87097;
etaS etaS [1 -1 -1 0 0 0 0] 0.002;
etaP etaP [1 -1 -1 0 0 0 0] 0.04;
lambda lambda [0 0 1 0 0 0 0] 0.6855;
alpha alpha [0 0 0 0 0 0 0] 0.5;
}

second
{
type Giesekus;
rho rho [1 -3 0 0 0 0 0] 803.87097;
etaS etaS [1 -1 -1 0 0 0 0] 0.002;
etaP etaP [1 -1 -1 0 0 0 0] 0.2324;
lambda lambda [0 0 1 0 0 0 0] 0.1396;
alpha alpha [0 0 0 0 0 0 0] 0.2;
}

third
{
type Giesekus;
rho rho [1 -3 0 0 0 0 0] 803.87097;
etaS etaS [1 -1 -1 0 0 0 0] 0.002;
etaP etaP [1 -1 -1 0 0 0 0] 0.5664;
lambda lambda [0 0 1 0 0 0 0] 0.0389;
alpha alpha [0 0 0 0 0 0 0] 0.3;
}

fourth
{
type Giesekus;
rho rho [1 -3 0 0 0 0 0] 803.87097;
etaS etaS [1 -1 -1 0 0 0 0] 0.002;
etaP etaP [1 -1 -1 0 0 0 0] 0.5850;
lambda lambda [0 0 1 0 0 0 0] 0.0059;
alpha alpha [0 0 0 0 0 0 0] 0.2;
}

);
}

If you change this file, do let me know how, please

JuHyun.
JuHyun is offline   Reply With Quote

Old   November 14, 2013, 05:02
Default
  #325
Senior Member
 
Dongyue Li
Join Date: Jun 2012
Location: Torino, Italy
Posts: 661
Rep Power: 8
sharonyue is on a distinguished road
Quote:
Originally Posted by JuHyun View Post
Hello Foamers,

i'm korean student, Want you to help me.

Perhaps it may be a very easy question.

i'm using the viscoelasticFluidFOAM multi-phase flow simulation tries.

but, this file. I did different from the content of the file.

'transportProperties' is the file I was using.

phase1
{
transportModel HerschelBulkley;
nu nu [ 0 2 -1 0 0 0 0 ] 3e-4;
rho rho [1 -3 0 0 0 0 0] 1000;

HerschelBulkleyCoeffs
{
nu0 nu0 [0 2 -1 0 0 0 0] 0.006;
tau0 tau0 [0 2 -2 0 0 0 0] 0.025031;
k k [0 2 -1 0 0 0 0] 3e-4;
n n [0 0 0 0 0 0 0] 1;
}
}
phase2
{
transportModel Newtonian;
nu nu [0 2 -1 0 0 0 0] 1.46e-05;
rho rho [1 -3 0 0 0 0 0] 1.225;
}

refPhase phase2;

sigma sigma [1 0 -2 0 0 0 0] 0.066;


but, the viscoelasticFluidFOAM in transportProperties without using viscoelasticProerties[Giesekus] is used.

I Have used a lot of different ways. but Always viscoelasticProperties in the error appeared.

viscoelasticProperties

rheology
{

type multiMode;

models
(
first
{
type Giesekus;
rho rho [1 -3 0 0 0 0 0] 803.87097;
etaS etaS [1 -1 -1 0 0 0 0] 0.002;
etaP etaP [1 -1 -1 0 0 0 0] 0.04;
lambda lambda [0 0 1 0 0 0 0] 0.6855;
alpha alpha [0 0 0 0 0 0 0] 0.5;
}

second
{
type Giesekus;
rho rho [1 -3 0 0 0 0 0] 803.87097;
etaS etaS [1 -1 -1 0 0 0 0] 0.002;
etaP etaP [1 -1 -1 0 0 0 0] 0.2324;
lambda lambda [0 0 1 0 0 0 0] 0.1396;
alpha alpha [0 0 0 0 0 0 0] 0.2;
}

third
{
type Giesekus;
rho rho [1 -3 0 0 0 0 0] 803.87097;
etaS etaS [1 -1 -1 0 0 0 0] 0.002;
etaP etaP [1 -1 -1 0 0 0 0] 0.5664;
lambda lambda [0 0 1 0 0 0 0] 0.0389;
alpha alpha [0 0 0 0 0 0 0] 0.3;
}

fourth
{
type Giesekus;
rho rho [1 -3 0 0 0 0 0] 803.87097;
etaS etaS [1 -1 -1 0 0 0 0] 0.002;
etaP etaP [1 -1 -1 0 0 0 0] 0.5850;
lambda lambda [0 0 1 0 0 0 0] 0.0059;
alpha alpha [0 0 0 0 0 0 0] 0.2;
}

);
}

If you change this file, do let me know how, please

JuHyun.
Hello,

I think it cannot simulate two-phase cases. You have to code your own solver using VOF of two-fluid model.
sharonyue is offline   Reply With Quote

Old   November 14, 2013, 05:15
Default
  #326
New Member
 
KimJuHyun
Join Date: Nov 2013
Posts: 2
Rep Power: 0
JuHyun is on a distinguished road
Quote:
Originally Posted by sharonyue View Post
Hello,

I think it cannot simulate two-phase cases. You have to code your own solver using VOF of two-fluid model.
I understand, I'll try it.

Thank you.
JuHyun is offline   Reply With Quote

Old   November 20, 2013, 16:24
Default
  #327
Member
 
Join Date: Aug 2013
Posts: 60
Rep Power: 3
sur4j is on a distinguished road
Hello,

I would like to simulate a polymer manufacturing process known as curing of rubber, this is a process where rubber is heated resulting in viscosity changes as the molecules react and form bonds (cross-link). When the rubber is initially heated its viscosity decreases and at a certain time starts to flow, some time after this the viscosity starts to increase (cross-linking occurs) up till the point it is solid rubber.

Can I use this solver to simulate this process or will I need a alternative method of simulation? I would appreciate any help at all with this. Thank you.
sur4j is offline   Reply With Quote

Old   December 4, 2013, 11:06
Default questions??
  #328
New Member
 
Heliana Cardenas
Join Date: Jul 2013
Posts: 19
Rep Power: 3
heliana60 is on a distinguished road
Quote:
Originally Posted by amin144 View Post
Hi dear FOAMers

I have basic questions, maybe very amateur questions.
1)What is tau in models like Giesekus? It's total stress or extra stress?

2)If I want calculate wall stress, can I use the utility "wallShearStress" or this give me only Newtonian part of stress?

3)what is best way to calculate eta and shear rate and plot diagram of them according to each other?


Best regards,
Amin
Hi Amin,

I am new to openFoam and using the viscoelasticFluidFoam solver. Could you calculate the total stress as Jovani told you to? or you found maybe another way? I was confused by the "tau" this solver uses because aparently it doesn't take into account the solvent stress. is that true?

Could you also calculate the shear rate?

heliana
heliana60 is offline   Reply With Quote

Old   January 30, 2014, 12:37
Default
  #329
Member
 
Nickolas P
Join Date: Oct 2010
Location: Greece
Posts: 30
Rep Power: 6
NickolasPl is on a distinguished road
Hello everyone,

I have a question regarding the viscoelasticFluidFoam. I understand that in the viscoelastic models (e.g. LPTT) the user must prescribe the BCs for "tau"". My question is if the tau boundary condition corresponds to the total stress tensor (as derived from the EVSS method) or the stress contribution from the polymer.
Then if we have to set the "tau" boundary condition on a free surface (e.g. the free surfaces for planar or annular extrudate swelling) what would reflect more realistic the true extrudate swelling concerning the free surface stresses...?

This.... type fixedValue;
value uniform (0 0 0 0 0 0);

????

..or this... type zeroGradient;

????

Thanking you in advance.

Nickolas
NickolasPl is offline   Reply With Quote

Old   February 9, 2014, 14:02
Unhappy Strain tensor , shear rate ?!
  #330
Member
 
ali jafari
Join Date: Sep 2012
Posts: 50
Rep Power: 4
ali jafari is on a distinguished road
hi every body

is there any utility that it can calculate strain tensor or velocity gradient for viscoelasticFluidFoam solver?
ali jafari is offline   Reply With Quote

Old   February 9, 2014, 14:42
Default
  #331
Member
 
s.rasoul_varedi
Join Date: Feb 2010
Posts: 82
Rep Power: 6
desert_1250 is an unknown quantity at this point
Send a message via Yahoo to desert_1250
Quote:
Originally Posted by ali jafari View Post
hi every body

is there any utility that it can calculate strain tensor or velocity gradient for viscoelasticFluidFoam solver?
Hi
use these keywords at terminal window:
foamCalc gradU // velocity gradient contour
foamCalc stressDifferences // normal stress differences N1 & N2
foamCalc wallGradU // velocity gradient at wall boundary
foamCalc components tau // all components of stress tensor
foamCalc components U // velocity components

Best
adambarfi and salame like this.
desert_1250 is offline   Reply With Quote

Old   March 7, 2014, 08:44
Default
  #332
New Member
 
Heliana Cardenas
Join Date: Jul 2013
Posts: 19
Rep Power: 3
heliana60 is on a distinguished road
Hey Jovani!

I hope you are fine. I have some problems with the relaxation factors that are used on your solver. They always affect the steady state (it takes longer to get to steady state than it should be), how do you choose them? When I use 0 for U i get a more or less nice result for the evolution of velocity on time, but for pressure or tau I am kinda lost because they are always unstable. Do you know how to choose them?

sorry this is all still kinda new for me

thanks
heliana60 is offline   Reply With Quote

Old   March 13, 2014, 21:10
Default
  #333
Member
 
Jovani L. Favero
Join Date: Mar 2009
Location: Rio de Janeiro, RJ, Brazil
Posts: 41
Rep Power: 9
jovani is on a distinguished road
Send a message via Skype™ to jovani
Quote:
Originally Posted by heliana60 View Post
Hey Jovani!

I hope you are fine. I have some problems with the relaxation factors that are used on your solver. They always affect the steady state (it takes longer to get to steady state than it should be), how do you choose them? When I use 0 for U i get a more or less nice result for the evolution of velocity on time, but for pressure or tau I am kinda lost because they are always unstable. Do you know how to choose them?

sorry this is all still kinda new for me

thanks
Hello Heliana,

I've never heard about a "magic" rule to choose an optimal relaxation factor, it depends on each case and can be changed along the simulation time. On the other hand, "0" makes not sense, take a look on section -- 4.5.2 Solution under-relaxation, http://www.openfoam.org/docs/user/fvSolution.php -- reproduced bellow. Small relaxation factor values increase the time to reach the steady-state solution but can improve numeric robustness.

Best regards,
Jovani

"4.5.2 Solution under-relaxation
A second sub-dictionary of fvSolution that is often used in OpenFOAM is relaxationFactors which controls under-relaxation, a technique used for improving stability of a computation, particularly in solving steady-state problems. Under-relaxation works by limiting the amount which a variable changes from one iteration to the next, either by modifying the solution matrix and source prior to solving for a field or by modifying the field directly. An under-relaxation factor α, 0 < α ≤ 1 specifies the amount of under-relaxation, ranging from none at all for α = 1 and increasing in strength as α → 0. The limiting case where α = 0 represents a solution which does not change at all with successive iterations. An optimum choice of α is one that is small enough to ensure stable computation but large enough to move the iterative process forward quickly; values of α as high as 0.9 can ensure stability in some cases and anything much below, say, 0.2 are prohibitively restrictive in slowing the iterative process."
jovani is offline   Reply With Quote

Old   March 14, 2014, 09:10
Default relaxation factors - viscoelasticFluidFoam - transients
  #334
New Member
 
Heliana Cardenas
Join Date: Jul 2013
Posts: 19
Rep Power: 3
heliana60 is on a distinguished road
Hi Jovani,

Yes, I read that too. For me your solver works really good regarding the results I get in SS, but the transients are what I would like to understand more and they are very much affected by the relaxation factors I use.

When I don't use them, and go to really small time steps, I still don't get the transient I expect in comparison to analytical results I have. Actually I get "wrong" results and my transients look like a the ones for Newtonian. I am trying to implement a nonlinear Maxwell model in a pressure driven flow.

Any comment is very much welcome,

Cheers,

Heliana
heliana60 is offline   Reply With Quote

Old   March 14, 2014, 12:26
Default viscoelasticinterfoam help wanted
  #335
New Member
 
abdikerim kurbanaliev
Join Date: Jun 2010
Location: Bloomington, IN, USA
Posts: 12
Rep Power: 7
kerim is on a distinguished road
Hi everyone!
now i have installed the openfoam 2.3 succeed, of course i download viscoelasticinterFoam files. but i do not kown how to add it in the openfoam. I think you can help me. Thanks a lot .
Have a nice day
kerim is offline   Reply With Quote

Old   March 14, 2014, 15:18
Default High stress values and instability of Viscoelastic Solver with VOF
  #336
Member
 
Ovie Doro
Join Date: Jul 2009
Posts: 99
Rep Power: 7
ovie is on a distinguished road
Hi all:

I am having issues with the stability of the coupled VOF - ViscoelasticFluidFoam solver. The problem typically starts with very high stress values at sharp corners of the geometry. This inturn makes the MULES solver unstable and the simulation crashes. Has anyone had similar problems? I know someone mentioned this high stress problem in another thread but I dont know if the problem has been solved yet.

I have also not been able to reproduce the results in the paper published by Jovani et al for the free surface viscoelastic solver. Has anyone succeeded with this?

Thanks.
ovie is offline   Reply With Quote

Old   March 14, 2014, 20:25
Default Looking for "viscoelasticinterFoam"
  #337
New Member
 
CFDGuy
Join Date: Nov 2013
Posts: 5
Rep Power: 3
mehdinekooee is on a distinguished road
Hi Foamers,

I am looking for the solver "viscoelasticinterFoam". What I can find online is "viscoelasticFluidFoam". Does anyone know how or where could I find viscoelasticinterFoam ?

Thanks
mehdinekooee is offline   Reply With Quote

Old   March 14, 2014, 20:31
Default How sensitive is the viscoelastic-VOF solver to geometry?
  #338
Member
 
Ovie Doro
Join Date: Jul 2009
Posts: 99
Rep Power: 7
ovie is on a distinguished road
Hi all:

I have also encountered the strange behaviour where my viscoelastic-VOF solver shows high sensibility to slight changes in geometry. Has anyone observed a similar problem? And did you manage to resolve it?
ovie is offline   Reply With Quote

Old   March 14, 2014, 20:41
Default
  #339
Member
 
Ovie Doro
Join Date: Jul 2009
Posts: 99
Rep Power: 7
ovie is on a distinguished road
Quote:
Originally Posted by mehdinekooee View Post
Hi Foamers,

I am looking for the solver "viscoelasticinterFoam". What I can find online is "viscoelasticFluidFoam". Does anyone know how or where could I find viscoelasticinterFoam ?

Thanks
Hi:

I dont think the viscoelasticInterFoam solver has been made available yet. So you might want to consider building one yourself by combining the viscoelasticFluidFoam with VOF from interFoam. Its not that difficult. I have manage to build one but its having serious stability issues with high viscoelastic stresses (see my posts on top).

Regards.
ovie is offline   Reply With Quote

Old   March 18, 2014, 06:43
Default foam-extend-3.0 preconditioner problem
  #340
Member
 
Sita Drost
Join Date: Mar 2009
Location: Arnhem, The Netherlands
Posts: 51
Rep Power: 8
sita is on a distinguished road
Dear all,

Recently, I upgraded to foam-extend-3.0. To test if viscoelasticFluidFoam works all right under this version, I tried to run the Oldroyd-B tutorial case. It gives an error regarding the Cholesky preconditioner ("unknown asymmetric matrix preconditioner..."). Was this preconditioner removed in foam-extend-3.0? And if so, which preconditioner is most suitable to use instead?

Thanks in advance for your help,

Regards,

Sita
sita 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
VOF simulation of a viscoelastic fluid sinah OpenFOAM Running, Solving & CFD 10 November 25, 2010 12:02
FREE SURFACE VISCOELASTIC FLOWS Valdemir G. Ferreira Main CFD Forum 6 December 18, 2009 07:14
Viscoelastic flow modeling in OpenFOAM vulda OpenFOAM Running, Solving & CFD 1 March 17, 2008 08:32
Polyflow & OpenFoam on Viscoelastic flow modeling Sumeshen Main CFD Forum 0 March 14, 2008 09:29
Viscoelastic fluid codes joel davison Main CFD Forum 0 November 6, 2001 06:09


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