The document contains MATLAB code for analyzing thermal and flow characteristics in a laminar flow scenario, focusing on the calculation of various parameters such as thermal boundary layer thickness, Nusselt number, and wall shear stress. It includes symbolic variable definitions, numerical methods for solving differential equations, and visualizations of results through plots. Additionally, it discusses the independence of the Nusselt number from initial conditions and provides references for further reading.
The document contains MATLAB code for analyzing thermal and flow characteristics in a laminar flow scenario, focusing on the calculation of various parameters such as thermal boundary layer thickness, Nusselt number, and wall shear stress. It includes symbolic variable definitions, numerical methods for solving differential equations, and visualizations of results through plots. Additionally, it discusses the independence of the Nusselt number from initial conditions and provides references for further reading.
The document contains MATLAB code for analyzing thermal and flow characteristics in a laminar flow scenario, focusing on the calculation of various parameters such as thermal boundary layer thickness, Nusselt number, and wall shear stress. It includes symbolic variable definitions, numerical methods for solving differential equations, and visualizations of results through plots. Additionally, it discusses the independence of the Nusselt number from initial conditions and provides references for further reading.
The document contains MATLAB code for analyzing thermal and flow characteristics in a laminar flow scenario, focusing on the calculation of various parameters such as thermal boundary layer thickness, Nusselt number, and wall shear stress. It includes symbolic variable definitions, numerical methods for solving differential equations, and visualizations of results through plots. Additionally, it discusses the independence of the Nusselt number from initial conditions and provides references for further reading.
% Constants U = 0.25; % m/s r0 = 0.1; % meter mu = 855E-6; % Dynamic viscosity of water nu = 855E-6/(1000/1.003); % Kinematic viscosity of water k = 0.613; % W/(m.K), thermal conductivity of water Pr = 5.83; % Prandtl number of water
for i = 2:num_steps % Get previous values Z_prev = Z_values(i-1); x_prev = x_values(i-1); lamda_prev = lamda_values(i-1); H_prev = H_values(i-1);
% Update x x_current = x_values(i); a = 2 * U * sin(x_prev / r0); % Uinf b = 2 * U * cos(x_prev / r0) / r0; % dUinf/dx
% Solve the differential equation using ode23 for the current step odefun = @(x, Z) dZdx(x, Z, H_prev); [x_sol, Z_sol] = ode23(odefun, [x_prev, x_current], Z_prev); Z_values(i) = Z_sol(end);
% Define the equation involving lamda_f
eq = ((37/315 - lamda_f/945 - lamda_f^2/9072)^2 * (lamda_f / b) == Z_values(i)); lamda_sol = vpasolve(eq, lamda_f, lamda_prev); % Use previous lambda as initial guess
% Filter roots within the range -12 to 12
lamda_filtered = double(lamda_sol(lamda_sol >= -12 & lamda_sol <= 12 & imag(lamda_sol) == 0)); if isempty(lamda_filtered) disp(['No valid lambda found at iteration ', num2str(i)]); break; % Exit the loop if no valid solution is found end
% Update lamda_values lamda_values(i) = lamda_filtered(1); % Use the first valid solution
I think the jumps are due to cos term in the differential of the dUinf/dx term because jump at teta=90 degree does not appear on lambda graph.
Figure 3 – Shear stress vs. Theta
Figure 4 – Boundary layer thickness vs. Theta
Figure 5 – Thermal boundary layer thickness vs. Theta Figure 6 – Convection coefficient vs. Theta
Figure 7 – Nusselt number vs. Theta
Q2- Φ(η) graph for random Φ(0) condition:
Nu is independent from the value of the Φ at the center. All Φ(0) values give Nu=9.8696.
Figure 8 – Φ(η) vs. Theta
Q2 Code:
function find_nu % Initial guess for Nu Nu_guess = 6;
% Use the fzero function to find the root of the boundary condition function options = optimset('Display','iter'); Nu_corrected = fzero(@boundary_condition, Nu_guess, options);
fprintf('The true value of Nu is: %.10f\n', Nu_corrected);
% Solve the system with the found nu
[x_vals, y_vals] = solve_odes(Nu_corrected);
% Plot the solution
figure; plot(x_vals, y_vals(1,:), 'DisplayName', 'Φ(x)'); xlabel('η'); ylabel('Φ(η)'); title(sprintf('Solution of the Energy Equation (For Φ=T-Ts) with Nusselt Number=%.4f', Nu_corrected)); legend; grid on;
% Verify boundary conditions
y_at_1 = y_vals(1,end); y_prime_at_0 = y_vals(2,1); fprintf('y(1) = %.10f, y''(0) = %.10f\n', y_at_1, y_prime_at_0); end function res = boundary_condition(Nu) % Boundary condition to be satisfied at x = 1 [~, y_vals] = solve_odes(Nu); y_at_1 = y_vals(1,end); res = y_at_1; % We want this to be 0 end function [x_vals, y_vals] = solve_odes(Nu) % Define the system of ODEs function dydx = odes(x, y) dydx = [y(2); -Nu/4 * y(1)]; end % Initial conditions y0 = [1; 0]; % y(0) = 1 and y'(0) = 0
% Solve the system using ode45
[x_vals, y_vals] = ode45(@odes, [0, 1], y0);
% Transpose to match the MATLAB output format
y_vals = y_vals'; end
REFERENCES
Al Ahmar, R., & Majdalani, J. (2022). On the Kármán momentum-integral approach and the Pohlhausen paradox: Extension to a cylinder in Crossflow with a potential Farfield Motion. Physics of Fluids, 34(6). https://doi.org/10.1063/5.0096780
Mecili, M., & Mezaache, E. H. (2013). Slug flow-heat transfer in parallel plate microchannel including slip effects and axial conduction. Energy Procedia, 36, 268–277. https://doi.org/10.1016/j.egypro.2013.07.031
