Answered step by step
Verified Expert Solution
Link Copied!

Question

1 Approved Answer

Exercise 4-1: Now, compute the solution to the 2D heat equation on a circular disk in Matlab. Recall that the heat equation is given by:

image text in transcribed

Exercise 4-1: Now, compute the solution to the 2D heat equation on a circular disk in Matlab. Recall that the heat equation is given by: Ut = advu. For this problem, we will solve the heat equation using a finite difference scheme on a Cartesian grid, as in class. We will use a grid of 300 x 300 with the circular disk in the center. The radius of the circle is r =1, a = 1, and the domain is [-1.5, 1.5) in x and [-1.5, 1.5] in y. 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 find points that are outside the disk, because they satisfy x2 + y2 > 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(0) = cos(0). 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. Exercise 4-1: Now, compute the solution to the 2D heat equation on a circular disk in Matlab. Recall that the heat equation is given by: Ut = advu. For this problem, we will solve the heat equation using a finite difference scheme on a Cartesian grid, as in class. We will use a grid of 300 x 300 with the circular disk in the center. The radius of the circle is r =1, a = 1, and the domain is [-1.5, 1.5) in x and [-1.5, 1.5] in y. 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 find points that are outside the disk, because they satisfy x2 + y2 > 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(0) = cos(0). 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

Step by Step Solution

There are 3 Steps involved in it

Step: 1

blur-text-image

Get Instant Access to Expert-Tailored Solutions

See step-by-step solutions with expert insights and AI powered tools for academic success

Step: 2

blur-text-image_2

Step: 3

blur-text-image_3

Ace Your Homework with AI

Get the answers you need in no time with our AI-driven, step-by-step assistance

Get Started

Recommended Textbook for

More Books

Students also viewed these Databases questions

Question

5. Have you stressed the topics relevance to your audience?

Answered: 1 week ago