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 Tree21Likes

Reply
 
LinkBack Thread Tools Display Modes
Old   August 21, 2013, 20:03
Post
  #281
New Member
 
胡长旭
Join Date: Aug 2013
Posts: 26
Rep Power: 3
hcx552362 is on a distinguished road
Quote:
Originally Posted by adambarfi View Post
what do you mean of 'it doesn't work'?
any error?


Hi, Thanks a lot.
I have solved my question. I am so glad with your help.
Thanks again! but now another question is appear.
When I calculating the example of viscoelasticFluidFoam, I run this program,
code:
viscoelasticFluidFoam > log &
the results show that:
FOAM FATAL TO ARROR:
Unkown asymmetric matrix solver BICGSTAB
Valid asymmeyric matrix solvers are:
.........


what is wrong?
Have a nice day!
hcx552362 is offline   Reply With Quote

Old   August 21, 2013, 21:14
Default
  #282
Senior Member
 
adambarfi's Avatar
 
Mostafa Mahmoudi
Join Date: Jan 2012
Posts: 281
Rep Power: 5
adambarfi is on a distinguished road
Send a message via Yahoo to adambarfi Send a message via Skype™ to adambarfi
Quote:
Originally Posted by hcx552362 View Post
Hi, Thanks a lot.
I have solved my question. I am so glad with your help.
Thanks again! but now another question is appear.
When I calculating the example of viscoelasticFluidFoam, I run this program,
code:
viscoelasticFluidFoam > log &
the results show that:
FOAM FATAL TO ARROR:
Unkown asymmetric matrix solver BICGSTAB
Valid asymmeyric matrix solvers are:
.........


what is wrong?
Have a nice day!
hi
I didn't hear before about BICGSTAB matrix solver. go to solutionDict and look for which equation you used BICGSTAB solver. change it to a asymmetric matrix solver.
it had guid you and had awared you about valid solvers in the last error.
adambarfi is offline   Reply With Quote

Old   August 22, 2013, 02:58
Post
  #283
New Member
 
胡长旭
Join Date: Aug 2013
Posts: 26
Rep Power: 3
hcx552362 is on a distinguished road
Quote:
Originally Posted by adambarfi View Post
hi
I didn't hear before about BICGSTAB matrix solver. go to solutionDict and look for which equation you used BICGSTAB solver. change it to a asymmetric matrix solver.
it had guid you and had awared you about valid solvers in the last error.
Hi, I have sloved this question. My program is run already. Thanks a lot

Have a nice day!
hcx552362 is offline   Reply With Quote

Old   August 23, 2013, 02:05
Post the same
  #284
New Member
 
k zhang
Join Date: Jun 2012
Location: China
Posts: 9
Rep Power: 4
skykzhang is on a distinguished road
Quote:
Originally Posted by skykzhang View Post
i am interested in the viscoelasticFluid. But i feel uncerain about the N1, PSD and so on. ( N1=σ11-σ22) .
Like: If the N1 is greater , is that means the elastic characteristic of the poly melt powerfuler? Can someone help me clearly know the question? Or recommand any books or journals about these issues.
Best regards!
Hi, viscoelasticFormers, there is another question,
How can i get the parameters of the Giesekus Model about one melt through rheological experiments, like etaS, etaP, alfar, lambda. I want to fit the parameters through the polymat module.
skykzhang is offline   Reply With Quote

Old   August 23, 2013, 11:13
Default ViscoelasticInterFoam
  #285
New Member
 
Erick Lorenzato
Join Date: Mar 2013
Posts: 8
Rep Power: 3
EricKl. is on a distinguished road
Hello everyone,

I'm working with viscoelastic fluids in my thesis and i need compile viscoelasticInterFoam

Jovani sent me the files. They were: Case (the turorial); Solver (Tutorial - Make folder with; viscoelasticTwoPhaseModels - folder with make, and viscoelasticLaws viscoelasticModel).
I tried to install them as the tutorial viscoelasticFluidFoam (http://lists.cfd-online.com/pipermai...er/029862.html).
So I can not. Is generating a lot of mistakes.

Anyone have any tutorial or way to install it?


I put the folders as viscoelasticFluidFoam, but when I do the command
Wmake libso && && cd -

Generates this error

root @ admin-desktop: ~ / OpenFOAM/OpenFOAM-1.6-ext/src/transportModels/viscoelastic viscoelasticTwoPhaseModels # cd / && wmake libso && cd -
wmakeLnInclude: linking include files to. / lnInclude
make: *** No rule to make target `viscoelasticLaws / DCPP / DCPP.dep ', needed by` Make/linux64GccDPOpt/dependencies'. Stop.

Thanks in advance!
Erick
EricKl. is offline   Reply With Quote

Old   August 26, 2013, 08:31
Default
  #286
Member
 
Sita Drost
Join Date: Mar 2009
Location: Arnhem, The Netherlands
Posts: 46
Rep Power: 7
sita is on a distinguished road
Quote:
Originally Posted by Raymond.Leoi View Post
There are four modes for the default PTT-Exponential. If only one mode is concerned, how to disable the rest higher modes?
Dear Raymond,

In constant/viscoelasticProperties, you just need to remove the "type multiMode" part and only define one mode. It's all rather straightforward, and I think you can even check the viscoelasticProperties-files in the tutorials; there are both single- and multi-mode tutorials.
Furthermore, you need to substitute the taufirst, tausecond, etc. files under 0/ with a single tau file.

Good luck,
Sita
sita is offline   Reply With Quote

Old   August 26, 2013, 08:53
Default
  #287
Member
 
Sita Drost
Join Date: Mar 2009
Location: Arnhem, The Netherlands
Posts: 46
Rep Power: 7
sita is on a distinguished road
Quote:
Originally Posted by skykzhang View Post
Hi, viscoelasticFormers, there is another question,
How can i get the parameters of the Giesekus Model about one melt through rheological experiments, like etaS, etaP, alfar, lambda. I want to fit the parameters through the polymat module.
Hi,

I think you should start with a basic rheology textbook. Some good ones are: Dynamics of Polymeric Liquids (Part I and II), by Bird, Armstrong and Hassager, Rheology: Principles, Measurements and Applications, by Macosko, and An Introduction to Rheology, by Barnes, Hutton and Walters. Also very instructive: Constitutive Equations for Polymer Melts and Solutions, by Larson.


Good luck,
Sita
sita is offline   Reply With Quote

Old   August 26, 2013, 09:18
Default
  #288
Member
 
Join Date: Feb 2012
Posts: 55
Rep Power: 4
Raymond.Leoi is on a distinguished road
Quote:
Originally Posted by sita View Post
Dear Raymond,

In constant/viscoelasticProperties, you just need to remove the "type multiMode" part and only define one mode. It's all rather straightforward, and I think you can even check the viscoelasticProperties-files in the tutorials; there are both single- and multi-mode tutorials.
Furthermore, you need to substitute the taufirst, tausecond, etc. files under 0/ with a single tau file.

Good luck,
Sita
Got it, Sita. Cheers.
Raymond.Leoi is offline   Reply With Quote

Old   August 26, 2013, 10:47
Default Weird inlet flow
  #289
Member
 
Join Date: Feb 2012
Posts: 55
Rep Power: 4
Raymond.Leoi is on a distinguished road
I ran a case w/ cylinder using Oldroyd-B model and specified the fully-developed velocity profile for the inlet and zeroGradient for outlet. The velocity magnitude and streamline was shown in the attache snapshot. It's quite weird to see two circulation zones near the inlet. Any idea for this problem?
Attached Images
File Type: png Screen Shot 2013-08-26 at 15.39.07.png (80.0 KB, 25 views)
Raymond.Leoi is offline   Reply With Quote

Old   August 26, 2013, 21:50
Post
  #290
New Member
 
salame ama
Join Date: Dec 2012
Posts: 28
Rep Power: 3
salame is on a distinguished road
Quote:
Originally Posted by sita View Post
Hi,

I think you should start with a basic rheology textbook. Some good ones are: Dynamics of Polymeric Liquids (Part I and II), by Bird, Armstrong and Hassager, Rheology: Principles, Measurements and Applications, by Macosko, and An Introduction to Rheology, by Barnes, Hutton and Walters. Also very instructive: Constitutive Equations for Polymer Melts and Solutions, by Larson.


Good luck,
Sita

Hi, sita:
Thanks!
These days, i read some books you recommanded. If take the single-modle, I think just the alfar should be fitted. is that right?

Great regards!
salame
salame is offline   Reply With Quote

Old   August 27, 2013, 20:08
Post
  #291
New Member
 
胡长旭
Join Date: Aug 2013
Posts: 26
Rep Power: 3
hcx552362 is on a distinguished road
Hello, everyone!
Did anyone who used the S-MDCPP model in the OpenFoam. When I using this Model, the result is not good! maybe it is not converage!
What is wrong ?
Thanks!
hcx552362 is offline   Reply With Quote

Old   August 27, 2013, 20:13
Default
  #292
New Member
 
胡长旭
Join Date: Aug 2013
Posts: 26
Rep Power: 3
hcx552362 is on a distinguished road
Hi ,
Did you use the S-MDCPP model in the OpenFoam. When I using this model , the result is poor!
Thanks
hcx552362 is offline   Reply With Quote

Old   August 28, 2013, 07:08
Default
  #293
Member
 
Join Date: Feb 2012
Posts: 55
Rep Power: 4
Raymond.Leoi is on a distinguished road
Quote:
Originally Posted by hcx552362 View Post
Hello, everyone!
Did anyone who used the S-MDCPP model in the OpenFoam. When I using this Model, the result is not good! maybe it is not converage!
What is wrong ?
Thanks!
No clue about your problem. It would be helpful to give more details.
Raymond.Leoi is offline   Reply With Quote

Old   August 28, 2013, 20:28
Default
  #294
New Member
 
胡长旭
Join Date: Aug 2013
Posts: 26
Rep Power: 3
hcx552362 is on a distinguished road
Hi, Thanks for youe reply!
When I employing the S-MDCPP model to do this case, and I transfer the result file to Tecplot360, the result is serious. For example, U P and Tau. I think the pattern is not steady!
Thanks a lot!
hcx552362 is offline   Reply With Quote

Old   August 28, 2013, 23:46
Default
  #295
Senior Member
 
adambarfi's Avatar
 
Mostafa Mahmoudi
Join Date: Jan 2012
Posts: 281
Rep Power: 5
adambarfi is on a distinguished road
Send a message via Yahoo to adambarfi Send a message via Skype™ to adambarfi
Quote:
Originally Posted by hcx552362 View Post
Hi, Thanks for youe reply!
When I employing the S-MDCPP model to do this case, and I transfer the result file to Tecplot360, the result is serious. For example, U P and Tau. I think the pattern is not steady!
Thanks a lot!
how about the residuals and convergence?
adambarfi is offline   Reply With Quote

Old   August 29, 2013, 02:12
Default
  #296
New Member
 
胡长旭
Join Date: Aug 2013
Posts: 26
Rep Power: 3
hcx552362 is on a distinguished road
Hi
I used the default file. for p is 1*10-7, and for tau and U is 1*10-6. The controlDict file is also default.
Thanks
hcx552362 is offline   Reply With Quote

Old   August 29, 2013, 02:13
Post
  #297
New Member
 
胡长旭
Join Date: Aug 2013
Posts: 26
Rep Power: 3
hcx552362 is on a distinguished road
Quote:
Originally Posted by adambarfi View Post
how about the residuals and convergence?
Hi
I used the default file. for p is 1*10-7, and for tau and U is 1*10-6. The controlDict file is also default.
Thanks
hcx552362 is offline   Reply With Quote

Old   August 29, 2013, 05:31
Default
  #298
Senior Member
 
adambarfi's Avatar
 
Mostafa Mahmoudi
Join Date: Jan 2012
Posts: 281
Rep Power: 5
adambarfi is on a distinguished road
Send a message via Yahoo to adambarfi Send a message via Skype™ to adambarfi
Quote:
Originally Posted by hcx552362 View Post
Hi
I used the default file. for p is 1*10-7, and for tau and U is 1*10-6. The controlDict file is also default.
Thanks
I mean your residuals during your solution.
could you attach the log file?
adambarfi is offline   Reply With Quote

Old   August 29, 2013, 19:58
Default
  #299
New Member
 
胡长旭
Join Date: Aug 2013
Posts: 26
Rep Power: 3
hcx552362 is on a distinguished road
Quote:
Originally Posted by adambarfi View Post
I mean your residuals during your solution.
could you attach the log file?
Hi, The log file contains a lot of something, Maybe the constitutive equation and parameter has some problems!
Thanks
hcx552362 is offline   Reply With Quote

Old   August 31, 2013, 20:58
Post About the S-MDCPP model
  #300
New Member
 
胡长旭
Join Date: Aug 2013
Posts: 26
Rep Power: 3
hcx552362 is on a distinguished road
Hi, everyone
When I using the S-MDCPP model by the default to solve my question, The results of U, P and tau is serious. the pattern is riots. What is wrong with this problem? Maybe you can run it, and check the result. Can anyone help me with this question?
Thanks a lot
best wishes
Have a nice day
hcx552362 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 11:02
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 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 05:09


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