====== Differences ====== This shows you the differences between two versions of the page.
Next revision | Previous revision | ||
gibson:teaching:spring-2018:math445:lecture:cylinderflow [2018/04/02 11:38] gibson created |
gibson:teaching:spring-2018:math445:lecture:cylinderflow [2018/04/10 06:52] (current) gibson |
||
---|---|---|---|
Line 8: | Line 8: | ||
p &= 2 \frac{a^2}{r^2} \cos 2 \theta - \frac{a^4}{r^4} | p &= 2 \frac{a^2}{r^2} \cos 2 \theta - \frac{a^4}{r^4} | ||
\end{eqnarray*} | \end{eqnarray*} | ||
+ | where $v_0$ is the speed of the incoming fluid, $a$ is the radius of the cylinder centered at the origin, and $r,\theta$ are the polar coordinates of point $x,y$. | ||
+ | These functions are valid only for points outside the cylinder $(r>a)$. | ||
The following Matlab code plots the velocity as a quiver plot and the pressure with contours. | The following Matlab code plots the velocity as a quiver plot and the pressure with contours. | ||
- | <code matlab> | + | <file matlab plotcylinderflow.m> |
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% | %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% | ||
% Inviscid, irrotational flow around a cylinder | % Inviscid, irrotational flow around a cylinder | ||
Line 85: | Line 87: | ||
title('inviscid 2D flow around a cylinder') | title('inviscid 2D flow around a cylinder') | ||
- | </code> | + | </file> |
{{ :gibson:teaching:spring-2016:math445:lecture:cylinderflow.png?direct&500 |}} | {{ :gibson:teaching:spring-2016:math445:lecture:cylinderflow.png?direct&500 |}} |