CFD Online Logo CFD Online URL
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Running, Solving & CFD

Difficulties with viscoelasticFluidFoam solver

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

LinkBack Thread Tools Display Modes
Old   November 24, 2009, 18:43
Default Difficulties with viscoelasticFluidFoam solver
Senior Member
Antonio Martins
Join Date: Mar 2009
Location: Porto, Porto, Portugal
Posts: 112
Rep Power: 10
titio is on a distinguished road
Send a message via MSN to titio Send a message via Skype™ to titio
Hi foamers,

I am trying to use the viscoelasticFluidFoam to model the flow of complex fluids in microfluidic systems. To test the code, and to get an idea of its behaviour and performance, I am running the case of the flow between infinite plates, where several analytical solutions for different constitutive equation models are available.

So far, I got the following results of my comparisons.

- For low values of the Deborah number (De), simulated results agree well with the theory.
- When I increase the value of the De number, it starts to deviate, increasing the deviations with De.
- When I run the solver with a null value for the solvent viscosity, I only can get convergence for very small values of De.
-Depending on the value of the density, the solver converges or diverges.

I used the linear version of the PTT model in my simulations, based on the example given in the tutorial, changing of course the mesh.

After talking with some people with knowledge in the area, they refered that there may be a problem in the wall boundary conditions. In particular, for viscoelastic fluids it is not correct to impose zero gradient boundary conditions for the extra tension. I am reading the literature now concerning this question and it looks like one needs to implement specific boundary conditions in the walls for these fluids.

Does anyone has any idea what kind of boundary conditions I should use, and how to do it.

Thanks in advance for any assistance,


titio is offline   Reply With Quote

Old   November 25, 2009, 05:19
Default Difficulties with viscoelasticFluidFoam solver
Senior Member
ata's Avatar
ata kamyabi
Join Date: Aug 2009
Location: Kerman
Posts: 323
Rep Power: 11
ata is on a distinguished road
Hi titio
Thatís true. I already thought about this subject. In steady state you can use no-slip boundary condition in stationary walls and set div(tau)=-grad(p) that must be modified during the solving. I must tell you this is my opinion and I am not sure if it works or even correct.
Please inform me if it works.
Best regards

ata is offline   Reply With Quote

Old   March 26, 2012, 10:27
New Member
Michael Stiehm
Join Date: Sep 2010
Posts: 13
Rep Power: 9
miael is on a distinguished road
Hello everybody,

I just want o ask, if there is something new regarding theboundary conditions for viscoelastic simulations. I want to simulate a pipe-flow and I don't know, if my inletBC ( tau --> fixedValue uniform (0 0 0 0 0 0); ) and outletBC (tau--> zeroGradient are correct. For the walls I also use zeroGradient.


miael is offline   Reply With Quote


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
Working directory via command line Luiz CFX 4 March 6, 2011 21:02
why the solver reject it? Anyone with experience? bearcat CFX 6 April 28, 2008 14:08
compressible two phase flow in CFX4.4 youngan CFX 0 July 1, 2003 23:32
CFX 5.5 Roued CFX 1 October 2, 2001 16:49
Setting a B.C using UserFortran in 4.3 tokai CFX 10 July 17, 2001 16:25

All times are GMT -4. The time now is 03:59.