calculatedChaos Gaming

Sastora wrote:Liquid, you are a lie.

MY COVER IS BLOWN!1!

Late night win.

Late night fail

Win.

Nice placement of the math stuff Jesus. Here is some fun math:
% Initialize Variables
fintempnew=zeros(21,61); %create current fin node matrix
fintempnew(1:21,1:61)=300; %set matrix to ambient temp to start
fintempold=zeros(21,61); %create previous fin node matrix
fintempold(1:21,1:61)=300; %set initial temp in previous iteration
errorfin=zeros(21,61); %create fin error matrix
errorfin(1:11,21:61)=1; %set air temp error to acceptable
m = 1;
n = 1;
dx=.001;
dy=.001;
h=60;
k=13.4;
Tinfinite=300;
Tempbase=300;
qinbase=88.3573*(500-Tempbase);
epsilon=.00000025;
error=0;
iterations=0;
errorsum=0;
while errorsum < 1200
qinbase=88.3573*(500-Tempbase);
%Top Left Node
if m==1 && n==1
fintempnew(m,n)=(((qinbase*dx)/k)+fintempnew(m+1,n)+fintempnew(m,n+1))/2;
error=abs(fintempnew(m,n) - fintempold(m,n))/fintempnew(m,n);
if error < epsilon
errorfin(m,n)=1;
end
fintempold(m,n)=fintempnew(m,n);
n=2;
%Top Left Boundary
elseif m==1 && 2 <= n && n <= 20
fintempnew(m,n)=(2*fintempnew(m+1,n) + fintempnew(m,n-1) + fintempnew(m,n+1))/4;
error=abs(fintempnew(m,n) - fintempold(m,n))/fintempnew(m,n);
if error < epsilon
errorfin(m,n)=1;
end
fintempold(m,n)=fintempnew(m,n);
n=n+1;
%Top Right Node
elseif m == 1 && n == 21
fintempnew(m,n)=(fintempnew(m,n-1) + fintempnew(m+1,n) + ((h*dy*Tinfinite)/k))/(2+(h*dy)/k);
error=abs(fintempnew(m,n) - fintempold(m,n))/fintempnew(m,n);
if error < epsilon
errorfin(m,n)=1;
end
fintempold(m,n)=fintempnew(m,n);
n = 1;
m = 2;
%Right Hand Boundaries
elseif 2 <= m && m <= 10 && n == 21 || 12 <= m && m <= 20 && n == 61
fintempnew(m,n) = (2*fintempnew(m,n-1) + fintempnew(m-1,n) + fintempnew(m+1,n) + ((2*h*dy*Tinfinite)/k))/(4 + ((2*h*dy)/k));
error=abs(fintempnew(m,n) - fintempold(m,n))/fintempnew(m,n);
if error < epsilon
errorfin(m,n)=1;
end
fintempold(m,n)=fintempnew(m,n);
n=1;
m=m+1;
%Inside Corner Node
elseif m == 11 && n == 21
fintempnew(m,n) = (k*(fintempnew(m,n-1) + fintempnew(m+1,n) +.5*(fintempnew(m-1,n) + fintempnew(m,n+1)) + ((h*dy*Tinfinite)/k)))/(3*k + h*dy);
error=abs(fintempnew(m,n) - fintempold(m,n))/fintempnew(m,n);
if error < epsilon
errorfin(m,n)=1;
end
fintempold(m,n)=fintempnew(m,n);
n=n+1;
%Top Right Boundary
elseif m == 11 && 22 <= n && n <= 60
fintempnew(m,n) = (fintempnew(m,n-1) + (2*h*dx*Tinfinite)/k + fintempnew(m,n+1) + 2*fintempnew(m+1,n))/(4 + (2*h*dx)/k);
error=abs(fintempnew(m,n) - fintempold(m,n))/fintempnew(m,n);
if error < epsilon
errorfin(m,n)=1;
end
fintempold(m,n)=fintempnew(m,n);
n=n+1;
% Middle Right Node
elseif m == 11 && n == 61
fintempnew(m,n)=(fintempnew(m,n-1) + fintempnew(m+1,n) + (2*h*dx*Tinfinite)/k)/(2 + (2*h*dx)/k);
error=abs(fintempnew(m,n) - fintempold(m,n))/fintempnew(m,n);
if error < epsilon
errorfin(m,n)=1;
end
fintempold(m,n)=fintempnew(m,n);
n=1;
m=m+1;
% Bottom Right Node
elseif m == 21 && n == 61
fintempnew(m,n)=(fintempnew(m,n-1) + fintempnew(m-1,n) + ((h*dy*Tinfinite)/k))/(2+(h*dy)/k);
error=abs(fintempnew(m,n) - fintempold(m,n))/fintempnew(m,n);
if error < epsilon
errorfin(m,n)=1;
end
fintempold(m,n)=fintempnew(m,n);
m=1;
n=1;
clc
errorsum=sum(sum(errorfin));
iterations = iterations +1 %#ok<NOPTS>
% Bottom Boundary
elseif m == 21 && 2<= n && n <= 60
fintempnew(m,n) = (fintempnew(m,n-1) + 2*fintempnew(m-1,n) + fintempnew(m,n+1))/4;
error=abs(fintempnew(m,n) - fintempold(m,n))/fintempnew(m,n);
if error < epsilon
errorfin(m,n)=1;
end
fintempold(m,n)=fintempnew(m,n);
n=n+1;
% Bottom Left Node
elseif m == 21 && n == 1
fintempnew(m,n) = ((qinbase*dy)/k + fintempnew(m-1,n) + fintempnew(m,n+1))/2;
error=abs(fintempnew(m,n) - fintempold(m,n))/fintempnew(m,n);
if error < epsilon
errorfin(m,n)=1;
end
fintempold(m,n)=fintempnew(m,n);
n=n+1;
% Left Boundary
elseif 2 <= m && m <= 20 && n == 1
fintempnew(m,n) = ((2*qinbase*dy)/k + fintempnew(m-1,n) + 2*fintempnew(m,n+1) + fintempnew(m+1,n))/4;
error=abs(fintempnew(m,n) - fintempold(m,n))/fintempnew(m,n);
if error < epsilon
errorfin(m,n)=1;
end
fintempold(m,n)=fintempnew(m,n);
n=n+1;
% Inside Nodes
else
fintempnew(m,n) = (fintempnew(m,n+1) + fintempnew(m,n-1) + fintempnew(m+1,n) + fintempnew(m-1,n))/4;
error=abs(fintempnew(m,n) - fintempold(m,n))/fintempnew(m,n);
if error < epsilon
errorfin(m,n)=1;
end
fintempold(m,n)=fintempnew(m,n);
n=n+1;
end
insideplatetemp=mean(mean(fintempnew(1:21,1)));
Tempbase=insideplatetemp;
end
% Heat Xfer Variables
wattsbaseplate=0;
wattsfinside=0;
wattsfintip=0;
heatxferperfin=0;
totalheatxfer=0;
finbasetemp=0;
insideplatetemp=0;
averagefintemp=0;
finbasetemperature=0;
totalrelativeerror=0;
temperaturerange=0;
totalaverage=0;
% Heat Xfer Calculations
wattsbaseplate=h*.001*.2*(fintempnew(1,21)-Tinfinite)*.5 + sum(h*.001*.2*(fintempnew(2:10,21)-Tinfinite));
wattsfinside=sum(h*.001*.2*(fintempnew(11,21:61)-Tinfinite));
wattsfintip=sum(h*.001*.2*(fintempnew(12:20,61)-Tinfinite))+h*.001*.2*(fintempnew(21,61)-Tinfinite);
heatxferperfin=wattsbaseplate+wattsfinside+wattsfintip;
totalheatxfer=heatxferperfin*20;
finbasetemperature=mean(mean(fintempnew(11:21,21:22)));
insideplatetemp=mean(mean(fintempnew(1:21,1)));
averagefintemp=mean(mean(fintempnew(11:21,21:61)));
temperaturerange=range(range(fintempnew));
Tempbase=finbasetemperature;
% Output
clc
fintempnew %#ok<NOPTS>
fprintf('The calculation took %5.0f iterations.\n',iterations)
disp('Steady State Results:')
fprintf('The total heat transferred is %4.4f watts.\n',totalheatxfer)
fprintf('The average fin base temperature is %3.2f kelvin.\n',finbasetemperature)
fprintf('The inside plate temperature is %3.2f kelvin.\n',insideplatetemp)
fprintf('The average fin temperature is %3.2f kelvin. \n',averagefintemp)
fprintf('The total temperature range from inside surface to fin tip is %3.2f kelvin. \n', temperaturerange)
disp('The figure shows the center of the fin tip')
disp('at the top right corner of each plot.')
subplot(1,2,1)
pcolor(fintempnew')
colorbar
xlabel('mm')
ylabel('mm')
title('Temperature Distribution')
subplot(1,2,2)
contourf(fintempnew')
xlabel('mm')
ylabel('mm')
title('Temperature Ranged')
colorbar

Very late night post.
1+1= Piss off, this is a competition of endurance, not math. 'Cuz I'm cool.

The answer is 3!

The answer is that is a steady state heat transfer program. I didn't have the transient segment redily available, otherwise I would have posted it as well.

steiner949 wrote:The answer is that is a steady state heat transfer program. I didn't have the transient segment redily available, otherwise I would have posted it as well.

Otherwise known as 3!

Steady state means there is now change in state, therefor you are incorrect and ss=0

Laft 5 Daed turns! The safe audio brands Tame Frotress 2. Tame Frotress 2 prepares Laft 5 Daed inside the entertained castle. Tame Frotress 2 untidies Laft 5 Daed beneath an ideal dish.

???? ????

Hey dog, where'd you find your avatar? I know an owl freak i want to show it to.

A little late to the fight, but I'm win'en.

I cannot wait for the summer to start so I can be playing more often. Stupid work prevents me from playing durring the days.

steiner949 wrote:I cannot wait for the summer to start so I can be playing more often. Stupid work prevents me from playing durring the days.

AGREED!

ONLY A MONTH AWAY!

I NEED A SUMMER JOB! NOOO

Be a spambot?

Liquid Death wrote:Be a spambot?

nah, too much competition. I TYPED THIS ON MY PSP!

Hey, I think we should talk about something that will get this thread locked so I can win.

The most controversial topic I know.

When loading toilet paper into a dispenser, should the roll be in front rolling towards you, or towards the back rolling away from you?

I am a front man myself.

Front

How old were you when you lost your virginity?

15

Bragging?

Just using this thread as a more permanent version of the shoutbox for conversing...I mean...

i win
i'm winning
no you aren't, I am

is only mildly entertaining for a page or 2

your wrong now...

Haha, front all the way Jesus. It is retarded how people put it behind because if you don't pay attention you are ripping it more than you need when it is in the back, and no body likes dirty hands.

It's also easier to rip off if it's hanging from the front. I mean not that I am some sort of weakling, but I buy toilet paper with more than one ply, and if you are trying to rip off a piece of the good stuff, trying to do that from the backside you are going to end up with more than you need and that's just a waste of toilet paper.

Agreed, the school and work use the worst possible paper because they are too cheap. They say that if you want 2-ply, you need to fold it over and then you have 2-ply. Lazy SOB's!