Environmental Engineering Reference
In-Depth Information
Q = 1;
% source/ sink rate [L^3/T]
Qx0 = 0.2;
% baseflow [L^2/T]
i = linspace(1,3,50);
[x,y,z] = meshgrid (i,i,i);
r = sqrt((x-2).^2+(y-2).^2+(z-2).^2);
xslice = [1.3;1.7;2.4];
yslice = [3];
zslice = [1.1;1.9];
phi = -Qx0*x-Q/4/pi/r;
slice (x,y,z,phi,xslice,yslice,zslice); hold on;
[u,v,w] = gradient (-phi);
h = streamline (x,y,z,u,v,w,ones(1,50),i,i);
set (h,'Color','y')
or
coneplot (x,y,z,u,v,w,)
The M-file is included in the accompanying software under the name
'3D_flow.m'
.
Figure
14.7
shows the flow field produced by the M-file after some manual
changes using the MATLAB
Figure editor. The color range for the color map was
restricted and lines on slices were omitted.
®
References
Guyon E, Hulin JP, Petit L (1997) Hydrodynamik. Vieweg, Braunschweig, p 466, in German
Prandtl L, Tietjens OG (1934) Fundamentals of hydro- and aeromechanics. Dover, New York,
p 270