Nusselt number and skin friction for bvp4c

37 次查看（过去 30 天）
Syed Mohiuddin
Syed Mohiuddin 2024-10-29，7:01
编辑： Syed Mohiuddin 2024-10-30，9:33
Hi,
I need to include the nusselt number and skin friction in the program.
How should I include these in program. Please help me in this regard
Program is given below
close all
clc
phi = 0.01;
Betaf = 207;
Betas = 17;
Beta = 0.5;
kof = 0.613;
kos = 400;
m = -1;
bk= 0.5;
bv = 0.5;
GRT = 10;
GRC = 10;
Alpha = 0.5;
Rhof = 997.1;
Rhos = 8933;
P = 0.5;
n = 0.3;
Br = 0.01;
A1 = (1-phi).^2.5;
A2 = 1/(1 + 1/Beta);
A3 = (1-phi)+phi.*((Rhos.*Betas)./(Rhof.*Betaf));
G1 = A2.*A3;
G3 = A1.*A2;
G4 = (kos + 2*kof - 2*phi.*(kof - kos))/(kos + 2*kof + phi.*(kof - kos));
G5 = (1 + 1/Beta)./A1;
G6 = G5./G4;
dydx=@(x,y)[y(4);
y(5);
y(6);
bv.*y(4).*y(5)-(1+bv.*y(2)).*(GRT.*G1.*y(2)+GRC.*G1.*y(3)-G3.*P);
bk.*y(5).^2-Br.*G6.*y(4).^2-(bk-bv).*y(2).*Br.*G6.*y(4).^2+bk.*bv.*Br.*G6.*y(2).^2*y(4).^2;
Alpha.*y(3)];
BC = @(ya,yb)[ya(1);yb(1);ya(2)-(1+m);yb(2)-1.0;ya(3)-(1+n);yb(3)-1.0];
yinit = [0.01;0.01;0.01;0.01;0.01;0.01];
solint = bvpinit(linspace(-1,1,11),yinit);
S1 = bvp4c(dydx,BC,solint);

请先登录，再回答此问题。

采纳的回答

Torsten
Torsten 2024-10-29，9:30
Nu1 = -(1+bk*(1+m))*S1.y(5,1)
Nu2 = -(1+bk)*S1.y(5,end)
Tau1 = (1-bv*(1+m))*S1.y(4,1)
Tau2 = (1-bv)*S1.y(4,end)

更多回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 Simscape 的更多信息

标签

产品


版本

R2023b

Community Treasure Hunt

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

Start Hunting!

Translated by