2. For the 2D analysis, due to the double symmetry of the problem, consider a quarter of the model as shown on the right. Assume the symmetric lines as insulated boundary and use the grid provided in the Scheme ($\Delta x = 1$ m), except the grid next to the symmetric lines, $\Delta x = 0.5$ m.
a. For the steady state, show the discretization of the quarter model - please show the numbering system used to obtain the scheme. Separate the internal nodes from the boundary nodes.
b. In solving the temperature distribution, use Matlab coding that would allow mesh refinement. But remember that you must convert the system of linear equations into matrix formation.
c. For the transient analysis, consider the explicit and implicit schemes and experiment with the $\Delta t$ that would provide convergence results. Again, use Matlab coding to solve the temperature distribution. Assume that initially, the temperatures were all 25°C.
3. Perform the transient analysis of the quarter model (as you did in Problem 2) using ANSYS and COMSOL and using fine mesh. Report the following:
a. Verification of the results (particularly support reactions) using energy balance
b. Temperature distribution along the symmetry line (parallel to X) at steady state - compare this to results in Problems 1 and 2
c. Temperature vs time for the corner node of the outer concrete wall and corner of the coolant (interacting with the concrete).
d. Plane view of the temperature distribution at steady state - compare this to result in Problem 2
4. Extend (or extrude) the 2D model in the third direction by 5 meter and study the effects of end conditions by providing different boundary conditions at these ends, namely insulation on one end and convective boundary condition on the other end. Use the convective condition imposed in Problem 2 for this 3D model. And report the following:
a. Planar plot of temperature 2 meters from the insulation end
b. Planar plot of temperature at the convective surface
c. One plane plot of either ZY or ZX section - showing distribution from the insulation to the convective surface