Recent content by range.rover

Hi friends, i have developed an code for a non linear heat conduction in 2 dimensions with dirichlet boundary condition by finite difference method in Matlab. my code is running slow to give output. If anybody has any idea of solving this equation or have written any Code for this equation...

i appreciate for your conversation, i have done my homework and solved it numerically as you said taking 3*3 matrix and i happen to program my code in the same way i have calculated, now i get to know why you are facing the the problem,if you just run my program line to line you can check in...

Hey gsal, I appreciate your reply when i discritised this equation "d2t/dx2 + d2t/dx[SUP]2 = -Q " i wrote the linear eauation ((T(i+1,j) +T(i-1,j) +T(i,j+1)+T(i,j-1) -4*T(i,j))*k)/ h[SUP]2 = Q, above in my program h is the distance between the nodes this has lead me to the tridiagonal...

thanks for your reply i am solving in a two dimensional staedy state heat conduction equtation with a heat source at the center in square plate with dirichlet boundary conditions, i mean to say values on boundries of four sides are same and constant. dt2 /dx2 + dt2/dy2 = -Q(x,y) i solved it...

NO it has'nt solved my problem. if you can be kind enough i am posting my code have a look at it%%%%%% STRICTLY ONLY FOR ODD NUMBERS %%%%%%%%% clear all close all clc n = 153;% n grids and has n - 2 interior points per dimension x = linspace(0,1,n);% line in x direction of 0 to 1 and divided...

HI guys,this is my first programming experience , i have developed an MATLAB code for steady state heat conduction equation , on governing equation dt2 /dx2 + dt2/dy2 = -Q(x,y) i have solved this equation with finite difference method, As far as i know if we increase the mesh size it leads...

yes, you are true the whole term is divided by 4k my above equation is like T = (((a+b+c+d)*k)+(B(i,j)*dx^2))/(4*k) a = T((i+1,j) b =T(i-1,j) c =T(i,j+1) d =T(i,j-1) B(i,j) is right hand side where dx is space between the nodes.

Thanx for your reply Chestermiller, B(i,j) is the right hand side matrix, I just applied central defference to the equation and took out the value of point(i,j) Tnew(i,j) = ((T((i+1,j) +T(i-1,j) +T(i,j+1)+T(i,j-1) )*k)+(B(i,j)*〖dx〗^2)))/(4*k)

Hi guys , i am solving this equation by Finite difference method. (dt2/dx2 + dt2/dy2 )= -Q(x,y) i have developed a program on this to calculate the maximum temperature, when i change the mesh size the maximum temperature is also changing, Should the maximum temperature change with mesh...

the error is exponentially decreasing and some times it is raising up at certain values of n,like you can see in the images down below at n = 25.

I have already checked my program and found that my ittirations are going way beond than 100000,it is the reason why it is taking long time. I am trying a lot to vectorize my while loop but i am not able to do because of my lack of knowledge and experience in programming skills.Thats why i...

does negative divergence of gradient tempearature gives to lalace equation...? -div(∇T) = [∂^2T/∂x^2+∂^2T/∂y^2]

Hi guys,I am new for programming. I made my first attempt to solve an heat equation by finite difference method and wrote a code for it in Matlab.I got a solution but i need help from you guys to (vectorize) increase the speed of my program. ((∂^2 T)/(∂x^2 )+(∂^2 T)/(∂y^2 )+(∂^2 T)/(∂z^2...

Hi friends i am trying to solve heat equation for a single time domain,with Finite Difference Method. I wrote a program in MATLAB,it sucessfully worked, but as i increse the number of points, my program is taking lots of time to give out result. Can anybody help me in increasing the speed of my...