Nusselt number and skin friction for bvp4c

2 次查看（过去 30 天）
Syed Mohiuddin
Syed Mohiuddin 2024-10-29
编辑： Syed Mohiuddin 2024-10-30
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
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 中查找有关 Analog Devices ADALM1000 Support from Data Acquisition Toolbox 的更多信息

标签

产品


版本

R2023b

Community Treasure Hunt

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

Start Hunting!

Translated by