# Velocity variance of potential flow in a channel

Consider incompressible potential flow ($\vec{u} = \nabla \phi$) in a complex (no-penetration) geometry with an inlet and an outlet. At these locations, the potential-field value is set as a boundary condition with some difference $\Delta\phi = \phi_{in} - \phi_{out}$. After solving the Laplace equation for the potential ($\phi$), the total variance of the flow field ($E$) may be defined as the volume ($V$) integral $E = \iiint_V \vec{v}\cdot\vec{v} \mathrm{d}V$, and the signed flux ($\Phi$) as a surface integral of (say) the inlet area $A_{in}$, $\Phi = \iint_{A_{in}} \nabla\phi \cdot \vec{n} \mathrm{d}A$, with $\vec{n}$ the normal to the inlet plane. Given the system parameters $\Delta\phi, E$ and $\Phi$, we can define a single dimensionless number,

$\Pi = \frac{\Delta\phi\Phi}{E}.$

Intuitively, I assumed that this parameter could “summarize” the transport efficiency of a complex geometry by comparing the volumetric flow rate against the kinetic energy stored in the flow field. This should have be relevant for starting and/or pulsating flows. However, I found after some experimentation that this number is always $\Pi = 1$, which I now understand.

Using the product rule for the inner product ($\nabla \cdot s\vec{a} = \nabla s \cdot \vec{a} + s \nabla \cdot \vec{a}$), the fact that we study incompressible flow (i.e. $\nabla \cdot \nabla\phi = 0$) and Gauss’s divergence theorem, we can rewrite the volume integral,

$E = \iiint_V \nabla\phi \cdot \nabla \phi \mathrm{d}V = \iiint_V \nabla \cdot \phi\nabla\phi - \phi \nabla\cdot\nabla\phi \mathrm{d}V =$ $= \iiint_V \nabla \cdot \phi\nabla\phi \mathrm{d}V = \oiint_A \phi \nabla \phi \cdot \vec{n} \mathrm{d}A.$

This surface integral only has contributions at the inlet and outlet, as the $\nabla \phi$ vector is perpendicular to $\vec{n}$ on the channel’s walls. Such that,

$E = \iint_{A_{in}} \phi \nabla \phi \cdot \vec{n} \mathrm{d}A + \iint_{A_{out}} \phi \nabla \phi \cdot \vec{n} \mathrm{d}A$ $= \phi_{in} \iint_{A_{in}} \nabla \phi \cdot \vec{n} \mathrm{d}A + \phi_{out} \iint_{A_{out}} \nabla \phi \cdot \vec{n} \mathrm{d}A$ $= \phi_{in} \Phi - \phi_{out}\Phi = \Delta\phi \Phi,$

which explains the unitary value for the aforementioned dimensionless ratio. Well done vector calculus for providing clarity on this topic!

### But what does “summarize” a flowability of a complex geometry?

Say that for some non-zero value of $\Delta \phi$, we have some total flux $\Phi_1$ and variance $E_1$. We can compute the volume flux per unit of energy of a flow field with a “Normalized” variance $O = \frac{\Phi_1}{\sqrt{E}_1}$. This Omicron ($O$) value may be compared between various geometries to say something about their ability to facilitate flow with low kinetic energy.

# Game of life on a sparse grid

I am currently starting to implement an adaptive Cartesian binary space partitioning (CBSP) mesh for a PDE-solver project. For this purpose I want to combine potentially small sections of Cartesian meshes of $N \times M$ cells, where $N$ and $M$ can be two different powers of two. This entails a dynamic memory management such that only portions of the grid are allocated in memory. The movie below shows a Conway’s Game of Life session on a $128 \times 64$ cells domain on such a grid. The color coding shows the allocated cells (Blue, magenta and green), the cells on which the game is played (magenta and green) and the cells that are “alive” (green). Meaning that no memory is allocated for the cells in the red region.

The memory management is performed using a pointer tree. In this example the data of a field is accessed via 32 pointers, each pointing up to 32 other pointers which could point to arrays of 8 cells. The (seemingly excessive) ratio of the number of pointers in the tree and the array size is tuned for the most convincing movie. The next admissible ratio in the code would be 16/32 as $16^2\times32 = 128\times64$.

# Code a flow solver

I have finished most of the writing for the simple flow solver pages. The purpose of these pages is to educate any body who is interested in the workings of flow-problem solvers. In fact, I actually coded this solver in basilisk before attempting to write a useful one.

# Global temperature distribution

In preparation for a job interview, I made a movie which illustrated my so-called “research vision”. It is nice enough to share, and displayed below.

# Paramo Sheffield

I invested some time in setting up a work bench for future DIY
projects. For this purpose, I bought a small second-hand vice and gave
it some tender love and care. Altough I *tought* I was finished,
the Marktplaats.nl algorithm kept on pushing vice advertisements. When I
came across the bigger brother of my (now) little vice, I knew it was
meant to be…

The restoration entailed:

- Paint stripping
- Full decrease and clean
- Jaw-nut thread retreading
- Repaint
- Adding white highlights
- Soft-jaw upgrade(?) for the big vice.

Side note: The replaceable jaws are a joke. New jaws for the large one, cost more than an entire new vice (albeit made from Chinesium).