Majumdar's correction for underrelaxation
These days, I try to implement the approach proposed by S. Majumdar (Numerical heat transfer,vol13,pp125, 1988) to solve the underrelaxation factor problem with nonstaggered grids. But it does not work well. I really need your help to analyze my problem-- How to implement the term '(1-relaxation factor)*(V*[i+1/2]-f-V*[i]-V*[i+1])'. When I add it to my code, the results still change with the relaxation factor. Could you give me a favor for this?
|
Re: Majumdar's correction for underrelaxation
fluxMaj = (1 - urf) * ( flux - density ( UE * Ax + VE * Ay + WE * Az) )
flux = ( terms u would have calculated even otherwise) + fluxMaj ; // fluxMaj is addtional // UE VE WE are vel components at face E and Ax Ay Az are area component of the face // flux is flux thorugh the face hope this help |
All times are GMT -4. The time now is 16:47. |