heat diffusion

logistic_guy

Full Member
Joined
Apr 17, 2024
Messages
588
here is the question

Derive the solution of the two dimensional steady-state heat diffusion equation using the method of separation of variables with the boundary conditions in figure 1.16. Use this solution to find the temperature at the midpoint T(1,0.5)\displaystyle T(1,0.5) with the boundary temperatures in figure 1.17. Show that the temperature at the midpoint is 94.5C\displaystyle 94.5^{\circ}C by considering the first five nonzero terms of the infinite series that must be evaluated.

Hint: Use the transformation θTT1T2T1\displaystyle \theta \equiv \frac{T - T_1}{T_2 - T_1} and solve 2θx2+2θy2=0.\displaystyle \frac{\partial^2 \theta}{\partial x^2} + \frac{\partial^2 \theta}{\partial y^2} = 0.

figure 1.16

figure16.png
figure 1.17

figure17.png

my attemb

this is the solution i must get:eek:

solution.png

i think they want me to use the first part

θ(x,y)=2πn=1(1)n+1+1n\displaystyle \theta(x,y) = \frac{2}{\pi}\sum_{n=1}^{\infty}\frac{(-1)^{n+1} + 1}{n}

which look like infinite series of calculus
 
Last edited:
Read this thread, it will help you a lot to derive the solution.

thank

it look different and it don't explain any steps
 
Open the link in post #3. Example 2 is exactly the same as the partial differential equation in your OP. Just follow the steps. If you got stuck anywhere, tell us.
i open link, example 2 is very different
 
i open link, example 2 is very different
🤦‍♂️

It is the same, it just uses different variables. I will show you.

This is your partial differential equation:

2θx2+2θy2=0\displaystyle \frac{\partial^2 \theta}{\partial x^2} + \frac{\partial^2 \theta}{\partial y^2} = 0

θ(0,y)=0\displaystyle \theta(0,y) = 0
θ(L,y)=0\displaystyle \theta(L,y) = 0
θ(x,0)=0\displaystyle \theta(x,0) = 0
θ(x,W)=1\displaystyle \theta(x,W) = 1

This is the partial differential equation in Example 2:

2u3x2+2u3y2=0\displaystyle \frac{\partial^2 u_3}{\partial x^2} + \frac{\partial^2 u_3}{\partial y^2} = 0

u3(0,y)=0\displaystyle u_3(0,y) = 0
u3(L,y)=0\displaystyle u_3(L,y) = 0
u3(x,0)=0\displaystyle u_3(x,0) = 0
u3(x,H)=f2(x)\displaystyle u_3(x,H) = f_2(x)

Compare between them and you will notice only this difference:

u3=θ\displaystyle u_3 = \theta
H=W\displaystyle H = W
f2(x)=1\displaystyle f_2(x) = 1

Read the solution of Example 2, and anywhere you see u3,H, or f2(x),\displaystyle u_3, H, \ \text{or} \ f_2(x), just replace it with θ,W, or 1,\displaystyle \theta, W, \ \text{or} \ 1, respectively.

Note: I took the boundary conditions of your partial differential equation from figure 1.16.
 
🤦‍♂️

It is the same, it just uses different variables. I will show you.

This is your partial differential equation:

2θx2+2θy2=0\displaystyle \frac{\partial^2 \theta}{\partial x^2} + \frac{\partial^2 \theta}{\partial y^2} = 0

θ(0,y)=0\displaystyle \theta(0,y) = 0
θ(L,y)=0\displaystyle \theta(L,y) = 0
θ(x,0)=0\displaystyle \theta(x,0) = 0
θ(x,W)=1\displaystyle \theta(x,W) = 1

This is the partial differential equation in Example 2:

2u3x2+2u3y2=0\displaystyle \frac{\partial^2 u_3}{\partial x^2} + \frac{\partial^2 u_3}{\partial y^2} = 0

u3(0,y)=0\displaystyle u_3(0,y) = 0
u3(L,y)=0\displaystyle u_3(L,y) = 0
u3(x,0)=0\displaystyle u_3(x,0) = 0
u3(x,H)=f2(x)\displaystyle u_3(x,H) = f_2(x)

Compare between them and you will notice only this difference:

u3=θ\displaystyle u_3 = \theta
H=W\displaystyle H = W
f2(x)=1\displaystyle f_2(x) = 1

Read the solution of Example 2, and anywhere you see u3,H, or f2(x),\displaystyle u_3, H, \ \text{or} \ f_2(x), just replace it with θ,W, or 1,\displaystyle \theta, W, \ \text{or} \ 1, respectively.

Note: I took the boundary conditions of your partial differential equation from figure 1.16.
give me some time to read this
 
i'm back

you don't tell me about this variable

u3(x,y)=h(x)φ(y)\displaystyle u_3(x,y) = h(x)\varphi(y)

what is h(x)\displaystyle h(x)?
what is φ(x)\displaystyle \varphi(x)?
 
what is φ(x)\displaystyle \varphi(x)φ(x)?
There is no Φ(x) - it is Φ(y)

This is separable variable method of solving PDE.

If you are not familiar with it - do a google search with "separable variable method of solving PDE" as key words.

Please let us know what you find
 
There is no Φ(x) - it is Φ(y)

This is separable variable method of solving PDE.

If you are not familiar with it - do a google search with "separable variable method of solving PDE" as key words.

Please let us know what you find
i don't understand

i'm using this site

 
i'm back

you don't tell me about this variable

u3(x,y)=h(x)φ(y)\displaystyle u_3(x,y) = h(x)\varphi(y)

what is h(x)\displaystyle h(x)?
what is φ(x)\displaystyle \varphi(x)?
This is the main step in any PDE when solved by the method of separation of variables. You can choose any two variables you want, not necessarily h(x) and φ(y)h(x) \ \text{and} \ \varphi(y).

If I was the one who solved that PDE in that website, I would choose these two variables:

u3(x,y)=X(x)Y(y)u_3(x,y) = X(x)Y(y)

I advise you to do the same, so start with this:

θ(x,y)=X(x)Y(y)\theta(x,y) = X(x)Y(y)
 
This is the main step in any PDE when solved by the method of separation of variables. You can choose any two variables you want, not necessarily h(x) and φ(y)h(x) \ \text{and} \ \varphi(y).

If I was the one who solved that PDE in that website, I would choose these two variables:

u3(x,y)=X(x)Y(y)u_3(x,y) = X(x)Y(y)

I advise you to do the same, so start with this:

θ(x,y)=X(x)Y(y)\theta(x,y) = X(x)Y(y)
thank

i reach this step

θn(x,y)=Bnsinh(nπyL)sin(nπxL)\displaystyle \theta_n(x,y) = B_n \sinh(\frac{n\pi y}{L})\sin(\frac{n\pi x}{L})

the next step i don't understand:(
 
thank

i reach this step

θn(x,y)=Bnsinh(nπyL)sin(nπxL)\displaystyle \theta_n(x,y) = B_n \sinh(\frac{n\pi y}{L})\sin(\frac{n\pi x}{L})

the next step i don't understand:(
This step:

θn(x,y)=Bnsinh(nπyL)sin(nπxL),     n=1,2,3,\displaystyle \theta_n(x,y) = B_n \sinh\left(\frac{n\pi y}{L}\right)\sin\left(\frac{n\pi x}{L}\right), \ \ \ \ \ n = 1, 2, 3, \cdots

tells you that you have infinite eigenfunctions. Each one of them is a solution.

And this step:

θ(x,y)=n=1Bnsinh(nπyL)sin(nπxL)\displaystyle \theta(x,y) = \sum_{n = 1}^{\infty} B_n \sinh\left(\frac{n\pi y}{L}\right)\sin\left(\frac{n\pi x}{L}\right)

tells you that the summation of all eigenfunctions is also a solution which is called the Superposition Principle.

The next step that comes after the Superposition Principle is very important. It will help you to find the constant BnB_n by using the nonhomogeneous boundary condition, θ(x,W)=1\theta(x,W) = 1.

Continue.
 
Last edited:
Bn=2Lsinh(nπWH)0Lsin(nπxL)dx   n=1,2,3,...\displaystyle B_n = \frac{2}{L\sinh(\frac{n\pi W}{H})}\int_{0}^{L}\sin(\frac{n\pi x}{L}) dx \ \ \ n = 1,2,3,...

it say this is Bn\displaystyle B_n but i don't understand why:(
 
Bn=2Lsinh(nπWH)0Lsin(nπxL)dx   n=1,2,3,...\displaystyle B_n = \frac{2}{L\sinh(\frac{n\pi W}{H})}\int_{0}^{L}\sin(\frac{n\pi x}{L}) dx \ \ \ n = 1,2,3,...

it say this is Bn\displaystyle B_n but i don't understand why:(
For now, you don't need to know why BnB_n is equal to that integral. Let your goal be just to arrive to the same solution in post #1. Later, we would fill the gaps that you did not understand about the solution.

Your solution still does not look like the given one. And your task now is to solve this integral:

0Lsin(nπxL) dx= ?\displaystyle \int_{0}^{L}\sin\left(\frac{n\pi x}{L}\right) \ dx = \ ?
 
Last edited:
For now, you don't need to know why BnB_n is equal to that integral. Let your goal be just to arrive to the same solution in post #1. Later, we would fill the gaps that you did not understand about the solution.

Your solution still does not look like the given one. And your task now is to solve this integral:

0Lsin(nπxL) dx= ?\displaystyle \int_{0}^{L}\sin\left(\frac{n\pi x}{L}\right) \ dx = \ ?
this easy

0Lsin(nπxL)dx=cos(nπxL)+c\displaystyle \int_{0}^{L}\sin(\frac{n\pi x}{L}) dx = \cos(\frac{n\pi x}{L}) + c

i'm not sure if i should change the constant c\displaystyle c to L\displaystyle L. i think it make more sense to use the limit of integration, so the answer is cos(nπxL)+L\displaystyle \cos(\frac{n\pi x}{L}) + L

Bn=2Lsinh(nπWH)cos(nπxL)+L   n=1,2,3,...\displaystyle B_n = \frac{2}{L\sinh(\frac{n\pi W}{H})}\cos(\frac{n\pi x}{L}) + L \ \ \ n = 1,2,3,...

i don't know how to simplify it more:(
 
this easy

0Lsin(nπxL)dx=cos(nπxL)+c\displaystyle \int_{0}^{L}\sin(\frac{n\pi x}{L}) dx = \cos(\frac{n\pi x}{L}) + c

i'm not sure if i should change the constant c\displaystyle c to L\displaystyle L. i think it make more sense to use the limit of integration, so the answer is cos(nπxL)+L\displaystyle \cos(\frac{n\pi x}{L}) + L

Bn=2Lsinh(nπWH)cos(nπxL)+L   n=1,2,3,...\displaystyle B_n = \frac{2}{L\sinh(\frac{n\pi W}{H})}\cos(\frac{n\pi x}{L}) + L \ \ \ n = 1,2,3,...

i don't know how to simplify it more:(
no comment

🤐

-----------------------

0Lsin(nπxL) dx=L(1cosnπ)nπ\displaystyle \int_{0}^{L}\sin\left(\frac{n\pi x}{L}\right) \ dx = \frac{L(1 - \cos n\pi)}{n\pi}

Now think about it. What happens if nn is odd and what happens if nn is even?
 
Top