## Inclined flow without constant height assumption. Setting up the problem, but not actually solving it (reworking vorticity equations)

Posted by peeterjoot on March 20, 2012

## Pressure and vorticity equations with the non-linear term retained. (Reworked slightly)

In section 40-2 [1], the identity

is used to put the vorticity equation into a form with one additional portion expressed as a gradient. This is a superior way to handle the inertial term because the curl of that gradient is then killed off.

I’ve expressed the curl as a wedge product, and not a cross product (either works since they related by a constant duality transformation). With the wedge product the identity 4.37 has different signs. That is

Here I’ve used the Hestenes overdot notation [2] to mark the operational range of the gradient (i.e. indicating that the gradient acts only on one of the terms initially). That gives us

Navier-Stokes (for incompressible flows) now takes the form

where in our problem we have killed the time dependence and have

The divergence of unfortunately isn’t zero. For exposition purposes, let’s write this out explicitly as a function of the vorticity components

Expanding out that divergence we have

or

Let’s also, for exposition, expand out the curl of this remaining non-linear term in coordinates. Being a bit smarter this time, we can avoid expressing in terms of and and leave it as a bivector explicitly. We have

This is

The Navier-Stokes equations are now recast in terms of vorticity as

\begin{subequations}

\end{subequations}

Having restated things with the term moved to the RHS, let’s also now write out 4.44a and 4.44b in coordinate form (we want this for the 2D case). This is

\begin{subequations}

\end{subequations}

For our problem where we have only and components, and any operations are zero, we find

and

So we have

\begin{subequations}

\end{subequations}

Here I’ve used again, so that the pair of differential operators on the LHS of the respective equations above are

Note that since could equal zero (as in our problem) we will likely have additional work to ensure that any solution that we find to this set of equations is also still a solution to our original first order Navier-Stokes equation.

# Now what?

The strategy that I’d thought to attempt to tackle this problem, when I’d left it like 4.36 was something along the following lines

\begin{itemize}

\item First ignore the non-linear terms. Find solutions for the homogeneous vorticity and pressure Laplacian equations that satisfy our boundary value conditions, and use that to find a first solution for .

\item Use this to solve for and from the vorticity.

\end{itemize}

However, after reworking it using the identity found in Feynman’s dry water chapter, I think it’s best not to try to solve it yet, and study some more first. I have a feeling that there are likely more such techniques that have been developed that will be useful to know before I try to plow my way through things.

Regardless, it’s interesting to see just how tricky the equations of motion become when one doesn’t make unrealistic assumptions. I have a feeling that to actually attempt this specific problem, I may very well need a computer and numerical techniques.

# References

[1] R.P. Feynman, R.B. Leighton, and M.L. Sands. *Feynman lectures on physics.[Lectures on physics]*, chapter The flow of dry water. Addison-Wesley Publishing Company. Reading, Massachusetts, 1963.

[2] D. Hestenes. *New Foundations for Classical Mechanics*. Kluwer Academic Publishers, 1999.

## Leave a Reply