CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Main CFD Forum

mesh size & urf SIMPLE FVM MATLAB

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

Reply
 
LinkBack Thread Tools Display Modes
Old   August 14, 2012, 09:52
Default mesh size & urf SIMPLE FVM MATLAB
  #1
New Member
 
Michael Moor
Join Date: May 2012
Location: Ireland
Posts: 29
Rep Power: 4
michaelmoor.aero 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')
end
if MassResidual <= 10^-3
alphap = 0.1
disp('alphap is now 0.3')
end

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

I am also solving the equations using Jacobi:
kk=1;
ip=0;
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);
%end
%end
kk=kk+1 ;
mx = max(max(max(abs(pcheck-pdash))));
ip=ip+1;
if (kk>50&&mx<10^-6), break, end
end

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

I'd appreciate any assistance!

Best Regards,
Michael
0565113@gmail.com
michaelmoor.aero is offline   Reply With Quote

Old   August 14, 2012, 11:12
Default
  #2
New Member
 
Michael Moor
Join Date: May 2012
Location: Ireland
Posts: 29
Rep Power: 4
michaelmoor.aero 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:
kk=1;
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);
end
end
kk=kk+1; % Counter
end

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);
michaelmoor.aero is offline   Reply With Quote

Reply

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

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
Layers:problem with curvature giulio.topazio OpenFOAM Native Meshers: snappyHexMesh and Others 10 August 22, 2012 10:03
Proof reading matlab code for SIMPLE and FVM michaelmoor.aero Main CFD Forum 0 July 3, 2012 17:06
unstructure mesh size control near to wall using gambit lalit kumar FLUENT 4 July 6, 2010 16:13
Mesh size for particulate flow simulations Shahri Main CFD Forum 0 March 24, 2009 19:40
Guides for calculation of mesh size for combustion Luk CFX 4 March 17, 2008 04:34


All times are GMT -4. The time now is 03:30.