How do we solve the differential equation \(E^2(E^2\psi)=0\)?

We have, in an earlier article, defined the Stokes stream function for an incompressible fluid at the boundaries of *r* = a and *r* = ∞ where:

To develop a Stokes stream function that satisfies eq36 for *a* ≤ *r* ≤ ∞ , let’s consider a solution of the form:

where is a function of *r*.

Substitute eq37 in eq36 and noting that , we have:

Letting and solving eq38, we get *n* = ±1, 2, 4. Therefore, the functiontakes the form:

Substitute eq39 in eq37,

To satisfy the function at the top of the page at *r* → ∞, we let *D* = *0* and *C* = ** u**/2. Eq40 becomes:

Substitute eq41 in eq15 and eq16, we have:

To satisfy the function at the top of the page at *r* = *a*, *u _{r}* =

*0*and

*u*=

_{θ}*0*, i.e. velocity of the fluid on the surface of the object is zero relative to the velocity of the object. This is because the force of attraction between the fluid particles and solid particles of the object is greater than that between the fluid particles. Eq42 and eq43 become:

Solving the simultaneous equations gives:

Substitute *A* and *B* in eq41,

Therefore, we have a Stokes stream function that satisfies the differential equation eq36 for *a* ≤ *r* ≤ ∞ .

Substitute eq44 in eq15 and eq16, we have:

Eq45 and eq46 are the component flow velocities of the viscous incompressible fluid.

##### NExt article: derivation of stokes’ law

##### Previous article: derivation of the differential equation, *E*^{2}(*E*^{2}ψ) = *0*

*E*(

^{2}*E*) =

^{2}ψ*0*