Can you help me to fix my Matlab code, please? I get this error :… Can you help me to fix my Matlab code, please? I get this error : quench_1D_cylinder_v1File: quench_1D_cylinder_v1.m Line: 24 Column: 46Invalid expression. When calling a function or indexing a variable, use parentheses. Otherwise, checkfor mismatched delimiters. This is my Matlab code Based on the solution of the screenshot I have attached to this question: clcclearvarsclose all% % This code calculates the cooling rate of a cylinder considering% convective heat losses. It compares the prediction against % benchmark data. The material properties are assumed to be temperature% independent and does not account for radiation.%%% load benchmark datadata = load(‘Results_data.txt’);data_t = data(:,1);data_T = data(:,2:end);labelvec{1} = ‘a) inner radius’;labelvec{2} = ‘b) mid point’;labelvec{3} = ‘c) outer radius’;figure(1)for i = 1: 3 subplot(1,3,i) plot(data_t,data_T(:,i),’k’,’linewidth’,2 hold on title(char(labelvec(i)))end%% Material parametersrho = 7939.8; % Density (kg/m^3)Cp = 472.2; % Specific Heat (J/kg/K)k0 = 14.6; % Thermal conductivity (W/m/K)%% Geometryinner_radius = 2.5; % mouter_radius = inner_radius + 0.225; % m%% Thermocouplestcloc(1) = inner_radius + 0.003;tcloc(2) = inner_radius + (outer_radius -inner_radius)/2;tcloc(3) = outer_radius – 0.003;%% Boundary conditionsT0 = 900 + 273.15; % Initial temperature (K)Tinfinity = 20 + 273.15; % Environment temperature (K)htc_inner = 100; % left heat transfer coefficient W/m^2/Khtc_outer = 300; % right heat transfer coefficient W/m^2/K%% Simulation controln = 20; % Spatial discretizationdt = 1; % time step (s)dr = (outer_radius-inner_radius)/(n-1); % grid spacing (m)r = inner_radius:dr:outer_radius; % vector of grid points (m)%% initializetime = 0; % initialise time (s)T = ones(n,1)*T0; % initialise temperature (K)T_new = zeros(n,1); % initialise future temperature (K)%% Model preparationBiot_inner = htc_inner*dr./k0; % Biot number for inner radiusBiot_outer = htc_outer*dr./k0; % Biot number for outer radiusalpha = k0./rho./Cp; % thermal diffusivity (m^2/s)Fo = alpha*dt/(dr*dr); % Fourier number%% Record temperatureit_out = 1;time_out(it_out) = time;TempTC_out(it_out,1:length(tcloc)) = interp1(r,T,tcloc)-273.15;%% Model calculationwhile max(T) > 40+273.15 %% Update temperature field T = T_new; time = time + dt; %% Solve temperature field % first point An = 1 – 2*Fo – 2*Fo*Biot_inner; Anp1 = 2*Fo; bc = 2*Fo*Biot_inner; T_new(1) = T(1)*An + T(2)*Anp1 + Tinfinity*bc; % last point An = 1 – 2*Fo – 2*Fo*Biot_outer; Anm1 = 2*Fo; bc = 2*Fo*Biot_outer; T_new(n) = T(n)*An + T(n-1)*Anm1 + Tinfinity*bc; % mid points for i = 2: n An = 1-2*Fo; Anp1 = Fo*(1 + 1/2); Anm1 = Fo*(1 – 1/2); T_new(i) = An*T(i) + Anp1*T(i+1) + Anm1*T(i-1); end %% thermocouple data it_out = it_out + 1; time_out(it_out) = time; TempTC_out(it_out,1:length(tcloc)) = interp1(r,T,tcloc)-273.15; endfor i = 1: 3 figure(1) subplot(1,3,i) plot(time_out,TempTC_out(:,i),’–r’) hold offend 1D heat transfer with convection and radiation In this problem we model a 1Dproblem where heat is lost from the inner and outer diameter ofa cylinder duringquenching, considering heat loss through convection. We compare… Show more… Show more Boundary regions We are modelling a case where we can ignore heat fluxesacting in the y and z directions. Consider the heat flowing through the volumeillustrated in Figure 1. – Ax Consider the heat flowing through the fol… Show more… Show more Engineering & Technology Mechanical Engineering MATERIALS 5335
