Now, compute the solution to the 2D heat equation on a circular disk in Matlab. Recall that the heat equation is given by: ut = alpha^2 * grenadian^2 * u: For this problem, we will solve the heat equ
Recall that the heat equation is given by:
ut = alpha^2 * grenadian^2 * u:
For this problem, we will solve the heat equation using a fnite difference scheme on a Cartesian
grid, as in class. We will use a grid of 300 300 with the circular disk in the center. The radius of
the circle is r = 1, alpha = 1, and the domain is [-1:5; 1:5] in x and [-1:5; 1:5] in y. Numerical Analysis Homework Help.
You can impose the boundary conditions by enforcing the temperature at points that are outside
of the disk at the beginning of each new time-step. It should be easy to nd points that are outside
the disk, because they satisfy x^2 + y^2 > 1.
Simulate the unsteady heat equation for the following boundary conditions:
• The left half of the boundary of the disk is fixed at a temperature of u = 1 and the right
half of the boundary is fixed at u = 2. Try simulating this with zero initial conditions first.
Next, try initial conditions inside the disk where the top half is u = 1 and the bottom half
is u = 1.
• The temperature at the boundary of the disk is fixed at u(theta) = cos(theta).
Include your Matlab code and show some plots of your solutions to the heat equation. Plot the
temperature distribution for each case 1) early on in the diffusion process, 2) near steady state,
and 3) somewhere in the middle.
More details please see the picture.