|
[Sponsors] |
December 2, 2014, 23:07 |
transportModel problem in interFoam
|
#1 |
New Member
Yongxiao Wang
Join Date: Nov 2014
Posts: 16
Rep Power: 11 |
Hi,
I want to simulate a non-Newtonian fluid flow by using interFoam solver. I modified the transportProperties as follows: Code:
water { transportModel poweLaw; nu nu [ 0 2 -1 0 0 0 0 ] 0.001; rho rho [ 1 -3 0 0 0 0 0 ] 1000; powerLawCoeffs { k k [0 2 -1 0 0 0 0] 100; n n [0 0 0 0 0 0 0] 0.1; nuMin nuMin [0 2 -1 0 0 0 0] 1; nuMax nuMax [0 2 -1 0 0 0 0] 10000; } } But the result of alpha.* is keeping the initial state. In other words, the fluid does not flow. And other results (U/P) are also not correct. I modify the water transportModel to Newtonian, and run it again. Everything is fine. Maybe this problem is infantile, but as a beginner I still don't know where is wrong. Any help is greatly appreciated. Thanks! Yongxiao |
|
December 3, 2014, 04:00 |
|
#2 | |
Senior Member
|
Quote:
Do the following changes in your case file: Code:
{ transportModel powerLaw; nu nu [ 0 2 -1 0 0 0 0 ] 0.001; rho rho [ 1 -3 0 0 0 0 0 ] 1000; powerLawCoeffs { k k [0 2 -1 0 0 0 0] 100; n n [0 0 0 0 0 0 0] 0.1; nuMin nuMin [0 2 -1 0 0 0 0] 0.0001; nuMax nuMax [0 2 -1 0 0 0 0] 10000; } } - Best Luck! |
||
December 3, 2014, 06:45 |
|
#3 |
New Member
Yongxiao Wang
Join Date: Nov 2014
Posts: 16
Rep Power: 11 |
Hi Tushar,
Thank you for your reply. But it's still not working properly. Everything is like before. |
|
December 3, 2014, 07:48 |
|
#4 |
Senior Member
|
Hello Yanci,
Are you comparing your results with some standards (papers)? Also, check the values of the viscosity for both the Newtonian and Non-Newtonian fluids. What is the viscosity value for the Newtonian fluid when it runs fine? - Best Regards! |
|
December 3, 2014, 08:57 |
|
#5 |
New Member
Yongxiao Wang
Join Date: Nov 2014
Posts: 16
Rep Power: 11 |
Hi Tushar,
Yes,you are right. I must say I didn't compare my results with any standard. In the fact, this is just a test for my own transportModel. Because I encountered a similar problem, when I implement my viscosity model into the interFoam. So I don't care too much about the specific value in this test. But now I seem to know where is wrong after your reminder. I set a large value for flow index (k) according to the previous. And it makes the viscosity values become too large. In a short period of time,its shape does not change obviously only under the influence of gravity. So the result of alpha.* is keeping the initial state. What do you think, am I right? Thank you for your reply. Best Regards! Yongxiao |
|
December 3, 2014, 23:22 |
|
#6 | |
Senior Member
|
Quote:
Sorry for the late reply, I was offline. You himself found the answer. You are correct with n=0.1, you are simulating the case of shear-thinning fluid. The large value of flow consistency index makes your fluid very thick in the very short period of time (or, for initial fluid viscosity is very high). This could be the reason for it's inability to capture the physics. Anyways Best Luck for your future work. - Best Regards! |
||
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
Flow around Cylinder with interFoam (Flow Recovery Problem) | jimbean | OpenFOAM Running, Solving & CFD | 0 | February 28, 2014 10:22 |
Initialisation problem with interFoam | belkadi | OpenFOAM Programming & Development | 0 | December 2, 2013 05:52 |
Problem with deltaT in interFoam | fedarduino | OpenFOAM Running, Solving & CFD | 0 | November 11, 2013 11:23 |
Interfoam - Problem with mesh quality ? | danvica | OpenFOAM Running, Solving & CFD | 4 | April 9, 2012 13:58 |
Pressure problem in Interfoam | danvica | OpenFOAM Running, Solving & CFD | 12 | March 14, 2012 02:56 |