Finite Temperature Variation - Heat Transfer

1 次查看(过去 30 天)
I am writing a code that displays temperature variation. There is a triangle section cut off of a square(nxn) that is exposed to convection heat transfer. The part that is cut off is a 90 degree triangle between i and n-.6*i. I need help defining this section for every row because the slope is not in one row. The cut out part is after a certain meter distance.
  3 个评论
Kaelyn
Kaelyn 2013-10-13
Yes, the triangle is at the top right corner of the square(nxn). I am having trouble identifying the cut off if the triangle on each row. The area is 1meter by 1meter and the missing area cut off is 1/2*.4^2. Giving a 45 degree angle.
Image Analyst
Image Analyst 2013-10-13
编辑:Image Analyst 2013-10-13
If the triangle is formed by two sides and one side is "i" elements long, and the other side is "n-0.6*i" elements long, how is that a 45 degree angle for the hypoteneuse? To get 45 degrees, both sides would have to be the same length.

请先登录,再进行评论。

回答(3 个)

Image Analyst
Image Analyst 2013-10-13
编辑:Image Analyst 2013-10-13
Do you mean like this:
clc;
n = 100
% Create sample random data.
myArray = randi(9, [n, n], 'uint8') + 100;
subplot(2, 2, 1);
imshow(myArray);
title('Original Image', 'FontSize', 20);
% Chop off i by n-6*i triangle
i = 4
xvertices = [0, n-6*i, 0];
yvertices = [0,0, i];
mask = poly2mask(xvertices, yvertices, n, n);
subplot(2,2,2);
imshow(mask);
title('Mask Image', 'FontSize', 20);
out = myArray .* uint8(~mask);
subplot(2,2,3);
imshow(out);
title('Output Image', 'FontSize', 20);
% Enlarge figure to full screen.
set(gcf, 'Units', 'Normalized', 'OuterPosition', [0 0 1 1]);
Note: requires Image Processing Toolbox because of poly2mask.
  2 个评论
Image Analyst
Image Analyst 2013-10-13
Kaelyn's comment moved here, since it's not an answer to her original question, but a comment on mine.
Yes but n and i will change according to the size of the mesh. The cut off of the triangle will be identified and used to calculate other temperatures. So I have been trying to use the mod function.
Image Analyst
Image Analyst 2013-10-13
编辑:Image Analyst 2013-10-13
n and i are variables in my code. You can change them. You can do everything with a for loop, instead of poly2mask(), if you want.

请先登录,再进行评论。


Youssef  Khmou
Youssef Khmou 2013-10-13
编辑:Youssef Khmou 2013-10-13
Kaelyn,
OK the problem now is clear , here is fast way
N=500; % more resolution better 2D heat conduction resolution .
H=ones(N);
M=flipud(triu(H)); % TOP RIGHT TRIANGLE
surface(M);
shading interp
  6 个评论
Youssef  Khmou
Youssef Khmou 2013-10-13
编辑:Youssef Khmou 2013-10-13
N=100;
M=rot90(triu(ones(N),-60));
surface(M), shading interp
Youssef  Khmou
Youssef Khmou 2013-10-13
is the problem solved by this last instruction?

请先登录,再进行评论。


Youssef  Khmou
Youssef Khmou 2013-10-13
N=1000; % 1 meter sampled with 1000 Hz
p=0.6*N; % your 0.6 starting point .
M=zeros(N); % initial matrix
A=(N-p)/N; % the Coefficient for constructing the triangle
for n=1:N
M(n,1:A*n)=1;
end

类别

Help CenterFile Exchange 中查找有关 Geometry and Mesh 的更多信息

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by