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

mesh size & urf SIMPLE FVM MATLAB

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

LinkBack Thread Tools Search this Thread Display Modes
Old   August 14, 2012, 08:52
Default mesh size & urf SIMPLE FVM MATLAB
Michael Moor
Join Date: May 2012
Location: Ireland
Posts: 30
Rep Power: 13 is on a distinguished road
Hello everyone,
I have developed a laminar 2D solver for Poiseulle flow using FVM and SIMPLE, and the results match nicely to the Hagen-Poiseulle velocity profile, but if i increase the mesh size, the solution blows up... My results are coming out nicely, but i cannot seem to be able to increase the mesh size beyond 36*36... I have to decrease the pressure under-relaxation to 0.01 to get it to iterate, but then it just oscillates at a continuity residual of around e-2... do you know the reason for this and why it happens? or how to fix it?

I am trying to maybe implement a method such that:
if MassResidual <= 10^-2
alphap = 0.05
disp('alphap is now 0.1')
if MassResidual <= 10^-3
alphap = 0.1
disp('alphap is now 0.3')

in an effort to try and march the process along? is this a correct approach?

I am also solving the equations using Jacobi:
while kk<=200;

pcheck = pdash;
%for i=2:nXp-1;
%for j=2:nYp-1;
%pdash(i,j)=(bdash(i,j)+saE(i,j)*pdash(i+1,j)+saW( i,j)*pdash(i-1,j)+saN(i,j)*pdash(i,j+1)+saS(i,j)*pdash(i,j-1))/saP(i,j);
kk=kk+1 ;
mx = max(max(max(abs(pcheck-pdash))));
if (kk>50&&mx<10^-6), break, end

and for residuals:
for i=2:nXp-1;
for j=2:nYp-1;
MassResidual=max([MassResidual abs(bdash(i,j))]);

I'd appreciate any assistance!

Best Regards,
Michael is offline   Reply With Quote

Old   August 14, 2012, 10:12
Michael Moor
Join Date: May 2012
Location: Ireland
Posts: 30
Rep Power: 13 is on a distinguished road
Apologies, I did not add all of the relevant information...

This is the momentum equation that i solve for, including under-relaxation:
while kk<=100;
%Solve iteratively using the Jacobi Method to yield u, which will now
%be the most up to date value, i.e. u*
for i=3:nX-1;
for j=2:nYp-1;
u(i,j)=((biJ(i,j)+uaE(i,j)*u(i+1,j)+uaW(i,j)*u(i-1,j)+uaN(i,j)*u(i,j+1)+uaS(i,j)*u(i,j-1) + (1-alphau)*uaP(i,j)*u(i,j)/alphau))*alphau/uaP(i,j);
kk=kk+1; % Counter

and i correct the pressure with under-relaxation:
%% Correct the Pressure, including under-relaxation

pnew(2:nXp,2:nYp-1) = pstar(2:nXp,2:nYp-1) + alphap*pdash(2:nXp,2:nYp-1); is offline   Reply With Quote


finite volume method, mesh sizing, residuals, simple, under relaxation factors

Thread Tools Search this Thread
Search this Thread:

Advanced Search
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 Off
Pingbacks are On
Refbacks are On

Similar Threads
Thread Thread Starter Forum Replies Last Post
[snappyHexMesh] Layers:problem with curvature giulio.topazio OpenFOAM Meshing & Mesh Conversion 10 August 22, 2012 09:03
Proof reading matlab code for SIMPLE and FVM Main CFD Forum 0 July 3, 2012 16:06
unstructure mesh size control near to wall using gambit lalit kumar FLUENT 4 July 6, 2010 15:13
Mesh size for particulate flow simulations Shahri Main CFD Forum 0 March 24, 2009 17:40
Guides for calculation of mesh size for combustion Luk CFX 4 March 17, 2008 02:34

All times are GMT -4. The time now is 00:09.