airfoil.m
. Mark where a
cylinder potential is set in a complex airfoil.m
is astonishingly simple for the complexity of the
flow and graphics that it produces.
Then in Matlab, select parameters that produce the flow around a
slightly cambered Joukowski airfoil of roughly 10% thickness ratio
at 15 degrees angle of attack. List your parameters. Plot out the
picture of the streamlines and isobars. Next set the circulation to
zero by setting variables Gamma
and auto
both to zero
and replot to show the effect of not satisfying the Kutta condition.
Do the streamlines still come off smoothly from the trailing edge?
What happens to the pressure?
So, use the Matlab program streamlines.m
to first plot the streamlines for potential flow around a cylinder.
Then create similar programs to plot the streamlines for both
potential and Stokes flow around a sphere. You will need to put in
the appropriate streamfunctions at the mesh points; potential flow
around a sphere is in section 19.8, and Stokes flow in 21.8. You
also need to change the values of the streamfunction that are
plotted as streamlines; your streamlines should be equally spaced
far upstream. The simplest way to do so is write the axysymmetric
streamfunction far upstream, and then figure out what streamfunction
values over there correspond to values equal to
. Comment on the
differences in streamline spacing between potential and Stokes flow
near the surface. Also comment on the differences in streamline
curvature one or two radii away from the sphere.