CFD Online Logo CFD Online URL
Home > Forums > Main CFD Forum

solving non-newtonian fluid

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

LinkBack Thread Tools Display Modes
Old   April 1, 2003, 20:34
Default solving non-newtonian fluid
Posts: n/a

I want to solve the flow with non-newtonian fluid in a tube using numerical method. The carrau model is used as a constitute equation. mu=viscosity gamma=shear rate

(mu-mu_inf)/(mu_0 - mu_inf) = (1 + (lamuda*gamma)^2)^((n-1)/2)

The governing eq.

0=-dp/dz + 1/r * d(mu*r*dw/dr)/dr

I would like to solve this, but have no idea. For newtonian fluid, it's so simple to analyze this. The constitue eq. and governing eq. are coupled. How can I deal with this problem using numerical approach(tips or algorithm)? If anyone know how to solve this or some open sources availalbe on site or at ftp, please let me know. I would much appreciate your help.

Thanks in advance. Peace in God.
  Reply With Quote

Old   April 2, 2003, 13:37
Default Re: solving non-newtonian fluid
Posts: n/a
Your problem equations are: (1)constitutive equation; (2)momentum conservation equation; (3)mass conservation equation.

Solve these three equations by Newton method since they are non-linear. Your unknowns are pressure, velocity and viscosity. You don't need boundary codnition for your constitutive equation since it's a linear algebraic equation.

From mathematical view, this way should work. I don't know if it's practical or not in your case.
  Reply With Quote

Old   April 2, 2003, 20:01
Default Re: solving non-newtonian fluid
Posts: n/a
I just want to know the solving procedure in numerical approach. Normally, I prefer interative method(SOR), not newton method.
  Reply With Quote

Old   April 3, 2003, 04:48
Default Re: solving non-newtonian fluid
Posts: n/a
Why not just solve it by iteration/quadrature:

(1) Guess w (w_g) and calculate mu

(2) update w (w_*) by integrating the momentum equation (it's a simple quadrature).

(3) obtain a new guess for w from w = (1-q)w_* + q*w_g

(4) repeat from 1 until converged.

In 3 q is a relaxation parameter used to stop the iteration oscillating (which it will do if q=1 ).
  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
Orifice Plate with a fully developed flow - Problems with convergence jonmec OpenFOAM Running, Solving & CFD 3 July 28, 2011 05:24
calculation diverge after continue to run zhajingjing OpenFOAM 0 April 28, 2010 04:35
ForcesCoeffs ronaldo OpenFOAM 4 September 14, 2009 07:11
On the damBreak4phaseFine cases paean OpenFOAM Running, Solving & CFD 0 November 14, 2008 22:14
Could anybody help me see this error and give help liugx212 OpenFOAM Running, Solving & CFD 3 January 4, 2006 19:07

All times are GMT -4. The time now is 14:40.