Use quiver () function to plot vector field (xy plane) of O2 flux in lab5.

2 次查看（过去 30 天）
Joy Odigbo
Joy Odigbo 2021-11-9
What I have currently:
rt= 50; % in micrometers
mu= 3.75*10.^8; %in mPa
delta_p= 1000; %delta p mmHg/um
l= 200; %length of the vessel (0-200 um)
v= ((rt.^2)/4.*u).*delta_p.*(1-(r/rt).^2); %velocity of vessel
x=linspace(0,200,20); % x points for the meshgrid
y=linspace(-rt,rt,-200); %y point for the meshgrid
%graph selection
[X, Y]=meshgrid (x,y);
vrx=blood(Y,rt,mu,delta_p);
U=Vrx;
V=U;
V(:,:)=0;
%vrx=v;
figure (1)
quiver(X,Y,u,v);
%part b graphing
%velocity of blood function
function [Vrx]=blood(r,rt,mu,delta_p)
Vrx=((rt^2)./(4*mu))*delta_p*(1-r/rt.^2);
end
What I had for Lab 5:
function[j]=lab05(x0,y0)
syms x y;
%muscle cell radius (R) - um - 1*10^-6m
R=25*(10^(-6));
%diffusion coeff(D) - cm^2*s^-1 - (1*10^-2m)^2*s^-1
D=(10^(-5))*(10^(-4));
%oxygen concentration outside cell (co2)0 - mol O2*mL^-1
co2=2.22*(10^(-7));
%oxygen consumption rate (sigma)- mole*mL^-1*s^-1
sigma=1.65*10^-8;
%oxygen concentration equation(phi)
phi=((sigma*(R^2))/4*D)*((((x^2)+(y^2))/(R^2))-1)+co2;
rtest= (x0^2)+(y0^2);
if rtest < (R)^2
disp('Will Work')
else
rtest > (R)^2; disp('Wont Work')
end
gradphi=gradient(phi,[x,y]);
%diffusing oxygen flux equation(j)
j=-D*gradphi;
jr0=subs(j,{x,y},{x0,y0});
divergence(j,[x,y])-phi;
disp('Flux of diffusing oxygen:') end

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

回答（0 个）

类别

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

标签

Community Treasure Hunt

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

Start Hunting!

Translated by