Question

In this project, you will perform transient analysis on a L 3 cm steel cube. The steel has the following properties: k 50 WImK, p 7800 kg/m3, cp 500 J/kgK. The cube is initially at a uniform temperature of 800 deg C, then is dropped into an oil bath at To 40 deg C. I. Create a program that will solve the three-dimensional temperature distribution of the cube for dimensionless times of 0.01, 0.1, and 1, and Bi 0.1, 1, and 10 (a total of 9 configurations NOTE: Use the equations below for Tau and Bi. Use a grid size of 1 a. For each configuration, present a temperature contour plot of the temperature on one face of the cube, and at the center symmetry surface of the cube (18 figures total b. In a table (or series of tables), present the center temperature, the center edge temperature, and the corner temperature for each configuration c. Compare the data from your program to theoretically determined values using either the Heisler charts (or one term approximation method) or lumped system analysis. Explain any discrepancies described by the results. 2. For the 3 cases where Bi 10, calculate the amount of heat lost by the cube in your program. Hint: calculate the heat lost in each element using specific heat. Compare these values to those from the Heisler charts 3. For the 3 Biot number cases, determine the amount of time (in seconds) until the system reaches steady state. Compare these values to the expected time to steady state due to theory and comment on any differences

Need help. Stuck and not sure how to finish this program. I'm using MatLab R2016a and would appreciate any help. Thanks.

Answer 1

#include<iostream>
#include<stdlib.h>

using namespace std;

int main()

for(x=1; x<n+1; x++)

c=(n*4)-4;
for(int a1=1; a1<n+1; a1++)

for(int i=1; i<n+1; i++)