Generate a computer code to simulate the solidification

Question

Problem:

Generate a computer code to simulate the following solidification situation during a casting process:

The material is a well-known polymer known as PEEK (polyetheretherketone). The properties of this material are as follows:

k = 0.251 W/mK,   ρ = 1262 kg/m³   C_ps = C_pi = 1339 J/kgK

h_s = C_psT_s = 807,417 J/kg,   h_l = C_plT_l = 827,502 J/kg

Assume that Tinitial = 670 K everywhere (the pouring temperature) and that at time = 0 the wall temperatures on walls AB, BC, AD, and DCare set to Twall = 300 K (the mold temperature). Study this solidification problem numerically, and hand in the following:

• A copy of your source code

• Plots of Tcenter versus time and Tp versus time

• The total solidification time and final point of solidification within the product.

• Contour plots of solidus and liquidus interfaces for a selected time.

• A plot showing the temporal progression of solidus interface fronts during the casting process.