Viscoelastic Fluid Flows using OpenFOAM The solver viscoelasticFluidFoam

ata · May 1, 2011, 02:52

It seems that the problem of the code has been solved.
Enjoy

anijdon · May 2, 2011, 02:15

hello;
yes, it does well; I extremely thankful for your kindness;
kind regards

anijdon · May 2, 2011, 10:04

hello M.r kamyabi;
excuse me; I did some change in the solver's code and compiled it eventually, but it seems has main problem when running a case; It doesn't do well and gives this error:

nano.zip
I don't know what to do

; would you guide me?
thanks;

ata · May 3, 2011, 03:14

It seems that the problem of the transportProperties file has been solved.
Enjoy

/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 1.7.0 |
| \\ / A nd | Web: www.OpenFOAM.com |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
location "constant";
object transportProperties;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

transportModel Newtonian;

nuF nuF [0 2 -1 0 0 0 0] 9.772e-4;
kF kF [1 1 -3 -1 0 0 0] 0.6047;
rhoF rhoF [1 -3 0 0 0 0 0] 9.9778e2;
muF muF [1 -1 -1 0 0 0 0] 9.7937e-7;
CpF CpF [0 2 -2 -1 0 0 0] 4076;
DTF DT [0 2 -1 0 0 0 0] 1.4868e-7;
Kb Kb [1 2 -2 -1 0 0 0] 1.4868e-7;

kP kP [1 1 -3 -1 0 0 0] 401;
rhoP rhoP [1 -3 0 0 0 0 0] 8933;
dP dP [0 1 0 0 0 0 0] 0.00002;
CpP CpP [0 2 -2 -1 0 0 0] 385;

//*************************************
nu nu [0 2 -1 0 0 0 0] 9.772e-4;
//*************************************

// ************************************************** *********************** //

anijdon · May 3, 2011, 12:17

hello;
thank you for reply;I want it read the nu from that formula which determined in the createFields not from transportProperties; If I change it's name there would be no problem?
thanks

ata · May 4, 2011, 03:50

Hi
I think you should use "IOobject::NO_READ" for it.
Good luck

ata · May 9, 2011, 06:44

Hi dear FOAMers
I have two questions about viscoelasticFluidFoam solver.
1. What does "p.boundaryField().updateCoeffs();" do?
2. I think after the line "p.relax();" the line "p.correctBoundaryConditions();" must be added. Am I true?
Best regards

Ata

fredgauss · May 9, 2011, 07:30

How far is the developmentof the multiphase viscoelastic solver? Would you be releasing is sometime?

Quote:

Originally Posted by jovani

Hello,

What you need is a multiphase solver with viscoelastic treatment. The viscoelasticFluidFoam solver doesn't do this, because your formulation is to closed flow of 1 only fluid. You cannot use interFoam too, because it doesn't support viscoelastic flow, only multiphase newtonian or GNF fluid.

I have a solver to treat multiphase viscoelastic flow, but a can't release it now.

Jovani

ata · May 10, 2011, 02:19

Hi
I have a two phase viscoelastic solver too. I will release it after the end of my Ph.D.
Any response to my question?
Good luck

fredgauss · May 10, 2011, 03:21

Ata,
What methodology are you using for your solver? When will you finish your Ph.D.? We tried some development here but it didn't work at least for high Weissenberg numbers. How does your solver perform?

ata · May 13, 2011, 03:14

Hi
I used VOF method. I used two solvers interFoam and viscoelasticFluidFoam for creating my new solver. I will finish my Ph.D. about 2012. Use low tolerance for tau equation (e.g. 1e-15). May be this solve your problem.
Good luck

Ata

usefkeyghobadi · August 1, 2011, 16:33

hi
any body have the PTT-Linear Model of viscoelasticfluidFoam paper?
PTT-Linear used in above solver is many different to another same of this model
i cant understand " taufirst"
please help me.
thank you

ata · August 2, 2011, 00:27

Hi
taufirst and tausecond and ... are used for multy mode constitutive equations.
Good luck

Ata

usefkeyghobadi · August 2, 2011, 06:04

Hi
Dear Ata
in Constitutive Equation used just tau
foe example LPTT model
(1+(ε_k λ_k)/(ηp_k ) tr(τp_κ ) )τp_κ+λ_k τp_k=2ηp_k D
tau first?

usefkeyghobadi · August 2, 2011, 06:06

Hi
Dear Ata
in Constitutive Equation used just tau
foe example LPTT model
(1+(ε_k λ_k)/(ηp_k ) tr(τp_κ ) )τp_κ+λ_k τp_k=2ηp_k D
tau first?

desert_1250 · August 2, 2011, 06:47

Dear yosof,
tua is the summation of tau at each mode of viscoelastic fluid because of its relaxation spectra

best regard
___________
Rasoul

usefkeyghobadi · August 2, 2011, 14:49

Quote:

Originally Posted by desert_1250

Dear yosof,
tua is the summation of tau at each mode of viscoelastic fluid becouse of its relaxation spectra

best regard
___________
Rasoul

Vielen Dank Herr Varedi!
bitte antwortest nschrischten(SMS!)

desert_1250 · August 3, 2011, 00:56

Quote:

Originally Posted by usefkeyghobadi

Vielen Dank Herr Varedi!
bitte antwortest nschrischten(SMS!)

Will fragen.
Viel Glück
Das Recht auf Asyl

___________
Rasoul

mmasoudian · August 21, 2011, 19:12

Dear Fomers
I am trying to compile VicoelasticFluidFoam in OpenFoam2.0.0 I read the instruction which proposed in forum but It was not compiled.
could you please let me know is this solver work with this version of Openfoam (2.0.0) or i have to downgrate it to OF1.5.
If your answer is yes could you tell me how I can compile it.
Regards

ata · August 22, 2011, 04:46

Hi
It is possible. Tell me what is your problem in detaile?
Good luck

Ata

May 13, 2011, 03:14	Two phase viscoelastic fluid flow	#151
ata Senior Member ata kamyabi Join Date: Aug 2009 Location: Kerman Posts: 323 Rep Power: 17	Hi I used VOF method. I used two solvers interFoam and viscoelasticFluidFoam for creating my new solver. I will finish my Ph.D. about 2012. Use low tolerance for tau equation (e.g. 1e-15). May be this solve your problem. Good luck Ata

August 2, 2011, 00:27	Viscoelastic Fluid Flows using OpenFOAM The solver viscoelasticFluidFoam	#153
ata Senior Member ata kamyabi Join Date: Aug 2009 Location: Kerman Posts: 323 Rep Power: 17	Hi taufirst and tausecond and ... are used for multy mode constitutive equations. Good luck Ata

August 2, 2011, 06:47		#156
desert_1250 Member s.rasoul_varedi Join Date: Feb 2010 Posts: 82 Rep Power: 15	Dear yosof, tua is the summation of tau at each mode of viscoelastic fluid because of its relaxation spectra best regard ___________ Rasoul Last edited by desert_1250; August 3, 2011 at 00:56.

August 21, 2011, 19:12	VicoelasticFluidFoam and OF2.0.0	#159
mmasoudian New Member mmam Join Date: Aug 2011 Location: Porto Posts: 3 Rep Power: 14	Dear Fomers I am trying to compile VicoelasticFluidFoam in OpenFoam2.0.0 I read the instruction which proposed in forum but It was not compiled. could you please let me know is this solver work with this version of Openfoam (2.0.0) or i have to downgrate it to OF1.5. If your answer is yes could you tell me how I can compile it. Regards

August 22, 2011, 04:46	VicoelasticFluidFoam in OpenFoam2.0.0	#160
ata Senior Member ata kamyabi Join Date: Aug 2009 Location: Kerman Posts: 323 Rep Power: 17	Hi It is possible. Tell me what is your problem in detaile? Good luck Ata

May 2, 2011, 02:15		#142
anijdon Member Maryam Mousazadeh Join Date: Oct 2010 Posts: 47 Rep Power: 15	hello; yes, it does well; I extremely thankful for your kindness; kind regards

May 2, 2011, 10:04		#143
anijdon Member Maryam Mousazadeh Join Date: Oct 2010 Posts: 47 Rep Power: 15	hello M.r kamyabi; excuse me; I did some change in the solver's code and compiled it eventually, but it seems has main problem when running a case; It doesn't do well and gives this error: nano.zip I don't know what to do ; would you guide me? thanks;

May 3, 2011, 03:14		#144
ata Senior Member ata kamyabi Join Date: Aug 2009 Location: Kerman Posts: 323 Rep Power: 17	It seems that the problem of the transportProperties file has been solved. Enjoy /--------------------------------- C++ -----------------------------------\ \| ========= \| \| \| \\ / F ield \| OpenFOAM: The Open Source CFD Toolbox \| \| \\ / O peration \| Version: 1.7.0 \| \| \\ / A nd \| Web: www.OpenFOAM.com \| \| \\/ M anipulation \| \| \---------------------------------------------------------------------------/ FoamFile { version 2.0; format ascii; class dictionary; location "constant"; object transportProperties; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // transportModel Newtonian; nuF nuF [0 2 -1 0 0 0 0] 9.772e-4; kF kF [1 1 -3 -1 0 0 0] 0.6047; rhoF rhoF [1 -3 0 0 0 0 0] 9.9778e2; muF muF [1 -1 -1 0 0 0 0] 9.7937e-7; CpF CpF [0 2 -2 -1 0 0 0] 4076; DTF DT [0 2 -1 0 0 0 0] 1.4868e-7; Kb Kb [1 2 -2 -1 0 0 0] 1.4868e-7; kP kP [1 1 -3 -1 0 0 0] 401; rhoP rhoP [1 -3 0 0 0 0 0] 8933; dP dP [0 1 0 0 0 0 0] 0.00002; CpP CpP [0 2 -2 -1 0 0 0] 385; //*********************************** nu nu [0 2 -1 0 0 0 0] 9.772e-4; //********************************* // ********************************************** ********************* //

May 3, 2011, 12:17		#145
anijdon Member Maryam Mousazadeh Join Date: Oct 2010 Posts: 47 Rep Power: 15	hello; thank you for reply;I want it read the nu from that formula which determined in the createFields not from transportProperties; If I change it's name there would be no problem? thanks

May 4, 2011, 03:50		#146
ata Senior Member ata kamyabi Join Date: Aug 2009 Location: Kerman Posts: 323 Rep Power: 17	Hi I think you should use "IOobject::NO_READ" for it. Good luck

May 9, 2011, 06:44		#147
ata Senior Member ata kamyabi Join Date: Aug 2009 Location: Kerman Posts: 323 Rep Power: 17	Hi dear FOAMers I have two questions about viscoelasticFluidFoam solver. 1. What does "p.boundaryField().updateCoeffs();" do? 2. I think after the line "p.relax();" the line "p.correctBoundaryConditions();" must be added. Am I true? Best regards Ata

May 10, 2011, 02:19		#149
ata Senior Member ata kamyabi Join Date: Aug 2009 Location: Kerman Posts: 323 Rep Power: 17	Hi I have a two phase viscoelastic solver too. I will release it after the end of my Ph.D. Any response to my question? Good luck

May 10, 2011, 03:21		#150
fredgauss New Member Ra Join Date: Apr 2010 Posts: 2 Rep Power: 0	Ata, What methodology are you using for your solver? When will you finish your Ph.D.? We tried some development here but it didn't work at least for high Weissenberg numbers. How does your solver perform?

August 1, 2011, 16:33		#152
usefkeyghobadi New Member Youse Keyghobadi Join Date: Jul 2011 Location: Shahrood-Iran Posts: 23 Rep Power: 14	hi any body have the PTT-Linear Model of viscoelasticfluidFoam paper? PTT-Linear used in above solver is many different to another same of this model i cant understand " taufirst" please help me. thank you

August 2, 2011, 06:04		#154
usefkeyghobadi New Member Youse Keyghobadi Join Date: Jul 2011 Location: Shahrood-Iran Posts: 23 Rep Power: 14	Hi Dear Ata in Constitutive Equation used just tau foe example LPTT model (1+(ε_k λ_k)/(ηp_k ) tr(τp_κ ) )τp_κ+λ_k τp_k=2ηp_k D tau first?

August 2, 2011, 06:06		#155
usefkeyghobadi New Member Youse Keyghobadi Join Date: Jul 2011 Location: Shahrood-Iran Posts: 23 Rep Power: 14	Hi Dear Ata in Constitutive Equation used just tau foe example LPTT model (1+(ε_k λ_k)/(ηp_k ) tr(τp_κ ) )τp_κ+λ_k τp_k=2ηp_k D tau first?

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
VOF simulation of a viscoelastic fluid	sinah	OpenFOAM Running, Solving & CFD	11	December 25, 2017 03:00
FREE SURFACE VISCOELASTIC FLOWS	Valdemir G. Ferreira	Main CFD Forum	6	December 18, 2009 06:14
Viscoelastic flow modeling in OpenFOAM	vulda	OpenFOAM Running, Solving & CFD	1	March 17, 2008 07:32
Polyflow & OpenFoam on Viscoelastic flow modeling	Sumeshen	Main CFD Forum	0	March 14, 2008 08:29
Viscoelastic fluid codes	joel davison	Main CFD Forum	0	November 6, 2001 05:09