Solving Flow Velocity in 2D for Flow Past a Circle

[riverboat]

Homework Statement

hello everybody,

I am trying to obtain the flow velocity in two dimensions u(x,y) for the case of a flow past a circle. The equations to solve are:

\vec{\nabla} p = \mu \vec{\nabla}^2 u

\vec{\nabla} \cdot \vec{u} = 0

I am very blocked at the moment and I don't become any idea at which can I start to calculate u(x,y). I want to solve it, but could some one give some sugerence to start to solve alone.

Homework Equations

The Attempt at a Solution

 

[CFDFEAGURU]

Do you have access to the book "Transport Phenomena" by Bird, Stewart, and Lightfoot?

 

[CFDFEAGURU]

First list you assumptions/postulates. This will allow you to stary simplfying the equations.

 

[kyiydnlm]

you can write down the stream function first since the flow passes a known circle; then, it won't be difficult.

 

[jalalmalo]

this reminds me of the heatring problem that Fourier studied

 

[riverboat]

My attempt following the case of 3 dimensions that I found in the book Fluid Mechanics by Kundu.

First of all, I have to obtain the stream function, for that we now

\vec{u} = \nabla \times \phi

in cartesians coordinates:
u_x = \frac{\partial \phi}{\partial y}
u_y = -\frac{\partial \phi}{\partial x}

The Navier - Stokes equation can be written as a function of the vorticity:
\vec{\nabla} p = \mu \vec{\nabla} \times (\vec{\nabla} \times \vec{u}) = \mu \vec{\nabla} \times \vec{w}

with \vec{w} the vorticity:

\vec{w} = -\frac{\partial^2 \phi}{\partial x^2} - \frac{\partial^2 \phi}{\partial y^2} = \vec{\nabla}^2 \phi

now I insert this result in the momentum equation and

\vec{\nabla} \times \vec{w} = \frac{\partial}{\partial y} \left( \frac{\partial^2 \phi}{\partial x^2}+ \frac{\partial^2 \phi}{\partial y^2}\right) e_x - \frac{\partial}{\partial x} \left( \frac{\partial^2 \phi}{\partial x^2}+ \frac{\partial^2 \phi}{\partial y^2}\right) e_y = \left( \frac{\partial^3 \phi}{\partial x^2 \partial y}+ \frac{\partial^3 \phi}{\partial y^3}\right) e_x - \left( \frac{\partial^3 \phi}{\partial x^3}+ \frac{\partial^3 \phi}{\partial y^2 \partial x}\right) e_y

\vec{\nabla} p = \mu \vec{\nabla} \times \vec{w} = \mu \left[ \left( \frac{\partial^3 \phi}{\partial x^2 \partial y}+ \frac{\partial^3 \phi}{\partial y^3}\right) e_x - \left( \frac{\partial^3 \phi}{\partial x^3}+ \frac{\partial^3 \phi}{\partial y^2 \partial x}\right) e_y \right]

Equating terms for e_x and e_y:
\frac{\partial p}{\partial x} = \mu \left( \frac{\partial^3 \phi}{\partial x^2 \partial y}+ \frac{\partial^3 \phi}{\partial y^3}\right)

\frac{\partial p}{\partial y} = \mu \left( \frac{\partial^3 \phi}{\partial x^3}+ \frac{\partial^3 \phi}{\partial y^2 \partial x}\right)

How can I resolve this two equations?

The next step is to fix the boundary conditions. In the circle surface it should be no slip \frac{\partial \phi}{\partial x} = \frac{\partial \phi}{\partial y} =0

and far of the circle is the boundary condition like this u_x = u_0 x and u_y = u_0 y

 

[gabbagabbahey]

You haven't stated whether you are trying to model a compressible or an incompressible flow yet (or any other assumptions you are making about the fluid flow!). For the latter case, what is \mathbf{\nabla}\times\textbf{u}?

 

[kyiydnlm]

hi, rb

from the second, third equations and the last two equations in your last line, exact stream function can be found. when you get the stream equation, you can put it into rest equations.

 

[riverboat]

Sorry, I forgot to said that I am considering a incompresible and viscous flow with very small Reynolds numbers.

\mathbf{\nabla}\times\textbf{u} is the curl for u, that you can become from the properties of the laplacian operator:
\mathbf{\nabla}^2 \mathbf{u} = \mathbf{\nabla} ( \mathbf{\nabla} \mathbf{u})- \mathbf{\nabla} \times ( \mathbf{\nabla} \times \mathbf{u})


Second, I am not sure if my boundary conditions far from the circle are correct or not.

 

[CFDFEAGURU]

This problem is laid out in great detail in the book "Transport Phenomena" by Bird, Stewart, and Lightfoot.

 

[minger]

I've cited this paper many of times. Taken from NASA Report No. 496, "General Theory of Aerodynamic Instability and the Mechanism of Flutter" by Theodorsen. (Available on the NASA technical report server).

\\Start quoted
Let us temporarily represent the wing by a circle. The potential of a source \epsilon at the origin is given by:
<br /> \phi = \frac{\epsilon}{4\pi}\log(x^2 + y^2)<br />
For a source \epsilon at (x_1,y_1) on the circle:
<br /> \phi = \frac{\epsilon}{4\pi}\log[(x-x_1)^2 + (y-y_1)^2)]<br />
Putting a double source 2\epsilon at (x_1,y_1) and a double negative source -2\epsilon at (x_1,-y_1) we obtain for the flow around a circle:
<br /> \phi = \frac{\epsilon}{4\pi}\log\frac{(x-x_1)^2 + (y-y_1)^2)}{x-x_1)^2 + (y+y_1)^2)}<br />
The function \phi on the circle gives directly the surface potential of a straight line pq, the projection of the circle on the horizontal diameter. In this case, y=\sqrt{1-x^2} and \phi is a function of x only.
\\end quote

The paper goes on to talk about pitching and plunging an actual airfoil shape, and the math involved is...fun. Either way, the link to the paper is (helpful for the figure)
http://ntrs.nasa.gov/archive/nasa/casi.ntrs.nasa.gov/19800006788_1980006788.pdf

Good luck!

 

[germana2006]

This problem is laid out in great detail in the book "Transport Phenomena" by Bird, Stewart, and Lightfoot.

In this book, they threat the problem in spherical coordinates (r, \vartheta, \phi) or cartesian coordinates (x,y,z) and the problem proposed here is to resolve in two dimensions only, polar coordinates (only r, \vartheta) or cartesian coordinates (x,y).

 

[CFDFEAGURU]

In the Transport Phenomena book this problem is outlined as an example of how to solve the flow of an inviscid fluid (one with negligible viscosity) through the use of a velocity potential. The problem is not solved in spherical cooridinates. (Why whould you use spherical coordinates for a cylindrical problem?)

Anyways, I am not going to post what the book has done but if you want to view it yourself then pick up a copy of the 2nd edition printed in 2007 and read section 4.3.

Thanks
Matt