# assemble it so that it will run

3 vues (au cours des 30 derniers jours)
MINATI PATRA le 11 Jan 2024
Commenté : MINATI PATRA le 11 Jan 2024
% this function estimates the correct value of a for big Pr cases
function f=big_Pr(a)
x0=[0 0 a(1) 1 a(2)];
[t,Y]=ode15s(@governing_equation,[0 7],x0);
% check the BC's at infinity approach zero
f(1)=Y(end,2);
f(2)=Y(end,4)^2+Y(end,5)^2;
end
% this function defines the similarity equations as defined by
% Bejan using correct scaling and Ostrach using incorrect scaling
function xdot=governing_equation(t,x)
global Pr
% Based on Bejan's scaling similarity equation
% The state-space derivation is shown in the document
xdot(1)=x(2);
xdot(2)=x(3);
xdot(3)=-(1/Pr)*(0.5*x(2)^2-(3*x(1)*x(3))/4)+x(4);
xdot(4)=x(5);
xdot(5)=3*x(1)*x(5)/4;
xdot=xdot';
% Based on Ostrach's scaling similarity equation
% % xdot(1)=x(2);
% % xdot(2)=x(3);
% % xdot(3)=-3*x(1)*x(3)+2*x(2)^2-x(4);
% % xdot(4)=x(5);
% % xdot(5)=-3*Pr*x(1)*x(5);
% % xdot=xdot';
end
global Pr
Pr = 1
options=optimset('TolFun',1e-6,'TolX',1e-6);
a=fsolve(@big_Pr,[0.5 -0.6],options);
a
x0=[0 0 a(1) 1 a(2)];
[t,Y]=ode15s(@governing_equation,[0 4],x0);
figure(1);
plot(t,-Y(:,2),'r')
hold on
figure(2);
plot(t,Y(:,4),'r')
hold on
Pr=100
options=optimset('TolFun',1e-6,'TolX',1e-6);
a=fsolve(@big_Pr,[0.25 -0.7],options);
a
x0=[0 0 a(1) 1 a(2)];
[t,Y]=ode15s(@governing_equation,[0 4],x0);
figure(1);
plot(t,-Y(:,2),'b')
figure(2);
plot(t,Y(:,4),'b')
Pr=1000
options=optimset('TolFun',1e-6,'TolX',1e-6);
a=fsolve(@big_Pr,[0.1 -0.65],options);
a
x0=[0 0 a(1) 1 a(2)];
[t,Y]=ode15s(@governing_equation,[0 4],x0);
figure(1);
plot(t,-Y(:,2),'k')
figure(2);
plot(t,Y(:,4),'k')
Pr=10
options=optimset('TolFun',1e-6,'TolX',1e-6);
a=fsolve(@big_Pr,[0.5 -0.52],options);
a
x0=[0 0 a(1) 1 a(2)];
[t,Y]=ode15s(@governing_equation,[0 4],x0);
figure(1);
plot(t,-Y(:,2),'g')
figure(2);
plot(t,Y(:,4),'g')
Pr=0.01
options=optimset('TolFun',1e-6,'TolX',1e-6);
a=fsolve(@small_Pr,[0.15 -0.17],options);
a
x0=[0 0 a(1) 1 a(2)];
[t,Y]=ode15s(@governing_equation,[0 4],x0);
figure(1);
plot(t,-Y(:,2),'y')
legend('Pr = 1','Pr = 100','Pr = 1000','Pr = 10','Pr = 0.01')
title('Vertical Velocity Profiles')
xlabel('η = (x/y) Ra_y^(^0^.^2^5^)')
ylabel('G = v(y/α) Ra_y^(^0^.^2^5^)')
figure(2);
plot(t,Y(:,4),'y')
legend('Pr = 1','Pr = 100','Pr = 1000','Pr = 10','Pr = 0.01')
title('Temperature Profiles')
xlabel('η = (x/y) Ra_y^(^0^.^2^5^)')
ylabel('ϴ')
function f=small_Pr(a)
x0=[0 0 a(1) 1 a(2)];
[t,Y]=ode15s(@governing_equation,[0 30],x0);
f(1)=Y(end,2);
f(2)=Y(end,4)^2+Y(end,5)^2;
end
%% I want to run this code, please assemble it so that it will run in 2023b version
##### 0 commentairesAfficher -2 commentaires plus anciensMasquer -2 commentaires plus anciens

Connectez-vous pour commenter.

### Réponse acceptée

Hassaan le 11 Jan 2024
clear
clc
close all
% Define global variable
global Pr
% Set the Prandtl number and solve for different cases
Pr = 1;
options = optimset('TolFun', 1e-6, 'TolX', 1e-6);
a = fsolve(@big_Pr, [0.5 -0.6], options);
Equation solved. fsolve completed because the vector of function values is near zero as measured by the value of the function tolerance, and the problem appears regular as measured by the gradient.
x0 = [0 0 a(1) 1 a(2)];
[t, Y] = ode15s(@governing_equation, [0 4], x0);
plot_results(t, Y, 'r', 'Pr = 1');
Pr = 100;
a = fsolve(@big_Pr, [0.25 -0.7], options);
Equation solved. fsolve completed because the vector of function values is near zero as measured by the value of the function tolerance, and the problem appears regular as measured by the gradient.
x0 = [0 0 a(1) 1 a(2)];
[t, Y] = ode15s(@governing_equation, [0 4], x0);
plot_results(t, Y, 'b', 'Pr = 100');
Pr = 1000;
a = fsolve(@big_Pr, [0.1 -0.65], options);
Equation solved. fsolve completed because the vector of function values is near zero as measured by the value of the function tolerance, and the problem appears regular as measured by the gradient.
x0 = [0 0 a(1) 1 a(2)];
[t, Y] = ode15s(@governing_equation, [0 4], x0);
plot_results(t, Y, 'k', 'Pr = 1000');
Pr = 10;
a = fsolve(@big_Pr, [0.5 -0.52], options);
Equation solved, inaccuracy possible. The vector of function values is near zero, as measured by the value of the function tolerance. However, the last step was ineffective.
x0 = [0 0 a(1) 1 a(2)];
[t, Y] = ode15s(@governing_equation, [0 4], x0);
plot_results(t, Y, 'g', 'Pr = 10');
Pr = 0.01;
a = fsolve(@small_Pr, [0.15 -0.17], options);
No solution found. fsolve stopped because the relative size of the current step is less than the value of the step size tolerance squared, but the vector of function values is not near zero as measured by the value of the function tolerance.
x0 = [0 0 a(1) 1 a(2)];
[t, Y] = ode15s(@governing_equation, [0 4], x0);
plot_results(t, Y, 'y', 'Pr = 0.01');
% Functions
function f = big_Pr(a)
x0 = [0 0 a(1) 1 a(2)];
[t, Y] = ode15s(@governing_equation, [0 7], x0);
f(1) = Y(end, 2);
f(2) = Y(end, 4)^2 + Y(end, 5)^2;
end
function f = small_Pr(a)
x0 = [0 0 a(1) 1 a(2)];
[t, Y] = ode15s(@governing_equation, [0 30], x0);
f(1) = Y(end, 2);
f(2) = Y(end, 4)^2 + Y(end, 5)^2;
end
function xdot = governing_equation(t, x)
global Pr
xdot(1) = x(2);
xdot(2) = x(3);
xdot(3) = -(1/Pr)*(0.5*x(2)^2 - (3*x(1)*x(3))/4) + x(4);
xdot(4) = x(5);
xdot(5) = 3*x(1)*x(5)/4;
xdot = xdot';
end
function plot_results(t, Y, color, label)
figure(1);
plot(t, -Y(:, 2), color);
hold on;
title('Vertical Velocity Profiles');
xlabel('η = (x/y) Ra_y^(^0^.^2^5^)');
ylabel('G = v(y/α) Ra_y^(^0^.^2^5^)');
legend(label);
figure(2);
plot(t, Y(:, 4), color);
hold on;
title('Temperature Profiles');
xlabel('η = (x/y) Ra_y^(^0^.^2^5^)');
ylabel('ϴ');
legend(label);
end
• For Cases Where a Solution is Found: The solutions seem fine, but you should still verify them against expected results or additional criteria to ensure they are accurate and meaningful.
• For the Case with Inaccuracy Possible: You might want to adjust the solver options like increasing the maximum number of iterations (MaxIter) or changing the function tolerance (TolFun). This could help the solver take more steps to find a better solution.
• For the No Solution Found Case: This is a bit trickier. You might need to provide a better initial guess or adjust solver settings. In some cases, the problem formulation itself might need to be revisited, especially if it's a particularly challenging nonlinear problem.
Remember, solving nonlinear equations, especially those from complex systems, can be quite sensitive to initial conditions and solver settings. Fine-tuning these parameters is often a key part of the process.
---------------------------------------------------------------------------------------------------------------------------------------------------
If you find the solution helpful and it resolves your issue, it would be greatly appreciated if you could accept the answer. Also, leaving an upvote and a comment are also wonderful ways to provide feedback.
Professional Interests
• Technical Services and Consulting
• Embedded Systems | Firmware Developement | Simulations
• Electrical and Electronics Engineering
Feel free to contact me.
##### 1 commentaireAfficher -1 commentaires plus anciensMasquer -1 commentaires plus anciens
MINATI PATRA le 11 Jan 2024
Thanks!

Connectez-vous pour commenter.

### Catégories

En savoir plus sur Engines & Motors dans Help Center et File Exchange

### Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by