Oseen's approximation

In 1910 Carl Wilhelm Oseen proposed Oseen's Approximation to treat problems in which a flow field involves a small disturbance of a constant mean flow, as in a stream of liquid. His work was based on the experiments of G.G. Stokes, who had studied a sphere of radius "$a$" falling in a fluid of viscosity ($\mu \,$ ). Oseen developed a correction term, which included inertial factors, for the velocity used in Stokes' calculations, to solve the problem. His approximation lead to an improvement to Stokes' calculations.

Importance

The method and formulation for analysis of flow at a very low Reynolds number is important. The slow motion of small particles in a fluid is common in bio-engineering. Oseen's drag formulation can be used in connection with flow of fluids under various special conditions, such as: containing particles, sedimentation of particles, centrifugation or ultracentrifugation of suspensions, colloids, and blood through isolation of tumors and antigens. The fluid does not even have to be a liquid, and the particles do not need to be solid. It can be used in a number of applications, such as smog formation and atomization of liquids.

Bio-engineering application

Blood flow in small vessels, such as capillaries, is characterized by small Reynolds and Womersley numbers. A vessel of diameter of 10 µm with a flow of 1 mm/second, viscosity of 0.02 poise for blood, density of 1 g/$cm^{3}$ and a heart rate of 2 Hz, will have a Reynolds number of 0.005 and a Womersley number of 0.0126. At these small Reynolds and Womersley numbers, the viscous effects of the fluid become predominant. Understanding the movement of these particles is essential for drug delivery and studying metastasis movements of cancers.

Calculations

Oseen considered the sphere to be stationary and the fluid to be flowing with a velocity ($U$) at an infinite distance from the sphere. Inertial terms were neglected in Stokes’ calculations. It is a limiting solution when the Reynolds number tends to zero. When the Reynolds number is small and finite, such as 0.1, correction for the inertial term is needed. Oseen substituted the following velocity values into the Navier-Stokes equations.

$v_{1}$ = $U$ + $v_{1}^{'}$, $v_{2}$ = $v_{2}'$, $v_{3}$ = $v_{3}'$

Inserting these into the Navier-Stokes equations and neglecting the quadratic terms in the primed quantities leads to the derivation of Oseen’s approximation:

$U{\partial v_{1}' \over \partial x_{1}}=-{1 \over \rho }{\partial p \over \partial x_{1}}+\nu \triangledown ^{2}v_{i}'\left({i=1,2,3}\right)$

When Stokes' solution was solved on the basis of Oseen's approximation, it showed that the resultant hydrodynamic force (drag) is given by

$F=6\pi \,\mu \,aU\left(1+{3 \over 8}N_{R}\right)$,

where:

• $N_{R}$ is the Reynolds number based on radius of the sphere
• $F$ is the hydrodynamic force
• $U$ is the flow
• $a$ is the radius of the sphere
• $\mu \,$ is the fluid viscosity

The force from Oseen's equation differs from that of Stokes by a factor of

$1+\left({3 \over 8}\right)N_{R}$.

Error in Stokes' solution

The Navier Stokes equations read:

$\triangledown v'~=0$

$v\triangledown v'~=-\triangledown p+\nu \triangledown ^{2}v'$.

but when the velocity field is:

$v_{y}=U\cos \theta \left({1+{a^{3} \over 2r^{3}}-{3a \over 2r}}\right)$

$v_{z}=-U\sin \theta \left({1-{a^{3} \over 4r^{3}}-{3a \over 4r}}\right)$

In the far field ${r \over a}$ >> 1, the viscous stress is dominated by the last term. That is:

$\triangledown ^{2}v'=O\left({a^{3} \over r^{3}}\right)$

The inertia term is dominated by the term:

$U{\partial v' \over \partial z_{1}}\sim O\left({a^{2} \over r^{2}}\right)$

The error is then given by the ratio:

$U{{\partial v' \over \partial z_{1}} \over {\nu \triangledown ^{2}v'}}=O\left({r \over a}\right)$

This becomes unbounded for ${r \over a}$ >> 1, therefore the inertia cannot be ignored in the far field. By taking the curl, Stokes equation gives $\triangledown ^{2}\zeta \,=0$ Since the body is a source of vorticity,$\zeta \,$ would become unbounded logarithmically for large ${r \over a}$. This is certainly unphysical and is known as Stokes’ paradox.

Modifications to Oseen's Approximation

One may question, however, whether the correction term was chosen by chance, because in a frame of reference moving with the sphere, the fluid near the sphere is almost at rest, and in that region inertial force is negligible and Stokes' equation is well justified. Far away from the sphere, the flow velocity approaches U and Oseen's approximation is more accurate. But Oseen's equation was obtained applying the equation for the entire flow field. This question was answered by Proudman and Pearson in 1957, who solved the Navier-Stokes equations and gave an improved Stokes' solution in the neighborhood of the sphere and an improved Oseen’s solution at infinity, and matched the two solutions in a supposed common region of their validity. They obtained:

$F=6\pi \,\mu \,aU\left(1+{3 \over 8}N_{R}+{9 \over 40}*N_{R}^{2}lnN_{R}+0(N_{R}^{2})\right)$

References

Sources

• Fung, Yuan-cheng. Biomechanics: Circulation. 2nd ed. New York, NY: Springer-Verlag, 1997.
• C.C. Mei, "Advanced Environmental Fluid Mechanics", Web.Mit.edu, 4 April 2011, webpage: MIT-Oseen (PDF).

External links

• "Advanced Environmental Fluid Mechanics", at Web.Mit.edu.