0.6 Solution of the partial differential equations (Page 8/13)

Page 8 / 13

\begin{matrix} \frac{\partial^{2} ϕ}{\partial x^{2}} + \frac{\partial^{2} ϕ}{\partial y^{2}} = 0 \\ \frac{\partial^{2} ψ}{\partial x^{2}} + \frac{\partial^{2} ψ}{\partial y^{2}} = 0 \end{matrix}

The Cauchy-Riemann conditions also imply that the gradient of the velocity potential and the stream function are orthogonal.

(\nabla ϕ) • (\nabla ψ) = \frac{\partial ϕ}{\partial x} \frac{\partial ψ}{\partial x} + \frac{\partial ϕ}{\partial y} \frac{\partial ψ}{\partial y} = 0

If the gradients are orthogonal then the equipotential lines and the streamlines are also orthogonal, with the exception of stagnation points where the velocity is zero.

Since the derivative

\frac{d w}{d z} = \begin{matrix} lim \\ |δ z \to 0| \end{matrix} \frac{δ w}{δ z}

is independent of the direction of the differential $δ z$ in the $(x, y)$ -plane, we may imagine the limit to be taken with $δ z$ remaining parallel to the x-axis $(δ z = δ x)$ giving

\frac{d w}{d z} = \frac{\partial ϕ}{\partial x} + i \frac{\partial ψ}{\partial x} = v_{x} - i v_{y}

Now choosing z to be parallel to the y-axis $(δ z = i δ y)$ ,

\frac{d w}{d z} = \frac{1}{i} \frac{\partial ϕ}{\partial y} + \frac{\partial ψ}{\partial y} = - i v_{y} + v_{x} = v_{x} - i v_{y}

These equations are a restatement that an analytical function has a derivative defined in the complex plane. Moreover, we see that the real part of $w^{'} (z)$ is equal to $v_{x}$ and the imaginary part of $w^{'} (z)$ is equal to $- v_{y}$ . If v is written for the magnitude of $v$ and for the angle between the direction of $v$ and the $x$ -axis, the expression for $d w / d z$ becomes

\begin{matrix} \frac{d w}{d z} = v_{x} - i v_{y} = v e^{- i θ} \\ or \\ v_{x} = real [\frac{d w}{d z}] \\ v_{y} = - imaginary [\frac{d w}{d z}] \end{matrix}

Flow Fields. The simplest flow field that we can imagine is just a constant translation, $w = (U - i V) z$ where $U$ and $V$ are real constants. The components of the velocity vector can be determined from the differential.

\begin{matrix} w (z) = (U - i V) z = (U - i V) (x + i y) = U x + V y + i (- V x + U y) = ϕ + i ψ \\ \frac{d w}{d z} = (U - i V) = v_{x} - i v_{y} \\ v_{x} = U, v_{y} = V \\ ϕ = U x + V y, ψ = - V x + U y \end{matrix}

Another simple function that is analytical with the exception at the origin is

\begin{matrix} w (z) & = & A z^{n} = A r^{n} e^{i n θ} \\ = & A r^{n} cos n θ + i A r^{n} sin n θ \\ = & ϕ + i ψ \\ thus \\ ϕ & = & A r^{n} cos n θ \\ ψ & = & A r^{n} sin n θ \\ \frac{d w}{d z} & = & n A z^{n - 1} = v_{x} - i v_{y} \end{matrix}

where $A$ and $n$ are real constants. The boundary condition at stationary solid surfaces for irrotational flow is that the normal component of velocity is zero or the surface coincides with a streamline. The expression above for the stream function is zero for all $r$ when $θ = 0$ and when $θ = π / n$ . Thus these equations describe the flow between these boundaries are illustrated below.

Fig. 6.5.1 (Batchelor, 1967) Irrotational flow in the region between two straight zero-flux boundaries intersecting at an angle $π / n$ .

Earlier we discussed the Green's function solution of a line source in two dimensions. The same solution can be found in the complex domain. A function that is analytical everywhere except the singularity at $z_{o}$ is the function for a line source of strength $m$ .

\begin{matrix} w (z) = \frac{m}{2 π} ln (z - z_{o}), line source \\ \frac{dw}{dz} = \frac{m}{2 π} \frac{1}{(z - z_{o})} = v_{x} - i v_{y} \end{matrix}

This results can be generalized to multiple line sources or sinks by superposition of solutions. A special case is that of a source and sink of the same magnitude.

\begin{matrix} w (z) = \sum_{i} \frac{m_{i}}{2 π} ln (z - z_{i}), multiple line sources \\ w (z) = \frac{m}{2 π} ln (\frac{z - z_{o}}{z + z_{o}}), source - sink pair \end{matrix}

The above flow fields can be viewed with the MATLAB code corner.m, linesource.m, and multiple.m in the complex subdirectory.

Assignment 7.5

Line Source Solution For

z_{o}

at the origin, derive expressions for the flow potential, stream function, components of velocity, and magnitude of velocity for the solution to the line source in terms of

r

and

θ

. Plot the flow potentials and stream functions. Compute and plot the flow potentials and stream function for the superposition of multiple line sources corresponding to the zero flux boundary conditions at

y = + 1

and

- 1

of the earlier assignment.

<< Chapter < Page Page > Chapter >>

Read also:

Get Jobilize Job Search Mobile App in your pocket Now!

100% Free Mobile Applications
Receive real-time job alerts and never miss the right job again

Source: OpenStax, Transport phenomena. OpenStax CNX. May 24, 2010 Download for free at http://cnx.org/content/col11205/1.1

Google Play and the Google Play logo are trademarks of Google Inc.

Notification Switch

Would you like to follow the 'Transport phenomena' conversation and receive update notifications?

Ask

	Cultural Evolution Contemporary World By Richley Crapo Start Assignment
	8 Java Persistence API By JavaChamp Team Start Quiz
©flickr: Abraham	Multicellular Organisms Test By Monty Hartfield Start Test
	NCE Ch 04 Social and Cultural Foundations By Anh Dao Start Quiz
	7 AP 07 Axial Skeleton MCQ By OpenStax Start Quiz
	3 Microbiology Final 3 By Madison Christian Start Test
	NCE Ch 10 Professional Orientation By Anh Dao Start Quiz
	Business Statistics By David Bourgeois Start Quiz
	9 Sociology 09 Social Stratification in the US MCQ By OpenStax Start Quiz
©flickr: Bertram	Chemistry Ch 1 2 Test 1 By Madison Christian Start Quiz