SIMPLE for variable viscosity, poor convergence ?
Dear all, I got a serious poor convergence problem while dealing with highly viscous fluid flow with viscosity depending on temperature. It seldom diverge actually but oscillate after the nomalized residual of velocities(or pressure) reduced to aroud 0.1. Collocated FVM and SIMPLE is adopted. A simple 3d tube geometry with cold wall temp and hot fluid is tested. The isothermal version of this code shows very good convergence. The nonlinearity caused by the vaiable viscosity may reduce the coupling ability of SIMPLE algorithm, I guess. Any comments or suggestions would be highly appreciated. Thanks.

Re: SIMPLE for variable viscosity, poor convergenc
Hi,
You really should'n have a problem with SIMPLE and variable viscosity  after all, our turbulence models introduce a variable effective viscosity AND are highly nonlinear (and everything works fine). I suspect you may have a problem with high viscosity: say, your diffusion Co number is very high and you're doing transient/underrelaxation based on the convection flux. If your viscosity is high (even locally), it will change the nature of the momentum equation from convection to diffusiondominated. Take a look at the Re number (is it ~1?) and the amount of work you're doing in the momentum predictor (it should be almost negligible compared to the rpessure solver). Hrv 
Re: SIMPLE for variable viscosity, poor convergenc
Thanks for your response. Numerically, the diffusoindominated problems should be more stable than the convectiondominated problems. In convectiondominated problem, the convective fluxe will contribute negtive values to the coefficient and reduce the diagonaldominated charateristic of the resulting matrix. That's why the upwinding technique is needed. The diffusiondominated problems will give an alway diagonaldominated matrix, and hence the high viscosity should stablize the computation. I use the "deferred correction" technique(1st upwind during iteration and 2nd CDS after converged) to retain the diagonal domination, and indeed very good convergence is observed during inner iteration for each velocity components (only 12 iterations are needed to achieve tolerence of 10e3 by BICGSTAB). However, the variable viscosity perturbs the resulting coefficient matrix considerablely during each SIMPLE outer iteration, and hence reduces the convergence. I ever tried to underrelax the viscosity but it only slightly improve the convergence.

Re: SIMPLE for variable viscosity, poor convergenc
Well, in that case very little can go wrong. Underrelaxation with standard relaxation parameters is tuned for incompressible turbulent flows; maybe you need to retune the relaxation parameters (any oscillatory behaviour in residuals?). The only other possibility is that your model is not well posed in the numerical sense or that there is more than one solution and you're switching between valid answers.
Hrv 
Re: SIMPLE for variable viscosity, poor convergenc
(1). If you plot the viscosity from cell to cell, what is the ratio of viscosity between cells. What is the range of this ratio? (2). What is the Reynods number of the flow? If it is turbulent, what turbulence model are you using?

Re: SIMPLE for variable viscosity, poor convergenc
If you are using the arthmetic interpolation for viscosity to obtain the nodal viscosity, you should try the harmonic averagebased interpolation. This make sense since viscosity is a diffusion coefficient ( same as thermal conductivity in heat conduction).. See Patanker for the Harmonic averaging.

Re: SIMPLE for variable viscosity, poor convergenc
hi John, the viscosity ratio is about 10 in the vincity of wall and 2 elsewhere. Acutally, my system is within creeping flow range. Full NavierStokes are solved.

Re: SIMPLE for variable viscosity, poor convergenc
hi Saud. Yes, I ever thinked that the harmonic interpolation should be better according to Patankar, and tried it, but no obvious improvement was obtained.

Re: SIMPLE for variable viscosity, poor convergenc
Hi Dr Jasak, can you tell me the range of the viscosity ratio within your system. will it diverge if the vsicosity ration is too large ? How you approximate the viscosity at the cell face center?

Re: SIMPLE for variable viscosity, poor convergenc
Hi Anthony,
you said that the isothermal case converged well. Have you tested your temp or scalar transport solver independantly of the variable visc part of your momentum solver? Say a given conservative velocity field or with constant viscosity? There are also some special issues with termperaturevelocity coupling see the work by Galpin and Raithby....sorry I do not have my reference list handy but I think in Numerical Heat Transfer about 1986...don't quote me on that! Regards..................Duane 
Re: SIMPLE for variable viscosity, poor convergenc
(1). There are possible areas: one is the algorithm itself (the form of the equation, the solution algorithm), the other is the boundary condition. (2). You can vary the temperature field (thus the viscosity variation) to check the sentivity of the solution. You can also refine the mesh in the area with high viscosity gradients.

Re: SIMPLE for variable viscosity, poor convergenc
Based on your answers to Mr. Chien's questions, I think you have similar problems with mine before. I have been doing NonNewtonian fluid flow simulations and they have the same poor convergence problem. My viscotiy varies enormously across yaxis for a channel(example). If I use a lower power index in PowerLaw fluid model, the convergence deteriorates seriously. It is well known phenomena and there are a few recommendations.
(1) decrease a relaxation parameter (2) if you use transient code, use smaller time step. (3) mix two nonlinear solvers such as Picard and NewtonRaphson. The former has bigger radius convergence. So use the Picard intiallly to get your intermediate solutions close engough to final solutions and switch to the NewtonRaphson, which has much faster convergence rate(2nd order). It is based on my experience(mainly FEM code). I could not find any more fancy or sophisticated ones that make sense to everyone. Jaret 
Re: SIMPLE for variable viscosity, poor convergenc
You may want to modify SIMPLE to SIMPLET. This modification intended to carry out the temperature effect more effectively when correcting the pressure

Re: SIMPLE for variable viscosity, poor convergenc
Hi, Duane. I ever thinked that It maybe the thermal boundary layer that reduces the convergence. So I tried a constant viscosity case, that is, velocity field is independent on temperature, and the convergence was very good just like that in the isothermal case. Both converged velocity and temperature fields can be obtained. Consquently, I think that it should be the variable viscosity that deteriorates the convergence. The variable viscosity problem should be very common in the field of crystal growth, soldification, casting and polymer processing. Is there anyone try solve these problem with FVM+SIMPLE?

Re: SIMPLE for variable viscosity, poor convergenc
hi Jaret. Actually, I am fighting the nonisothermal, nonNewtonian flow by FVM for my PhD program (Permanent Head Damage :( ). Several months ago, the suceess in the isothermal code really makes me filled with the confidence that I will soon finish my job. After the numerous fights with nonlinearity in the coutless dark nights,however, at present, I thinks I am still far way from the land of the Grail.

Literature for SIMPLET?
Hi, Saud
Can you suggest the literature for SIMPLET. Zeng 
All times are GMT 4. The time now is 09:00. 