Phyein Usul← Numa Koudsie

Research projects

Portfolios.

Projects of months rather than afternoons: a question from research, a model rebuilt and checked step by step, a solver you can run, and results you can reproduce. Each project keeps its own pages; this page introduces them.

01

Bachelor thesis, Leipzig 2022 · remodelled

cloud physics · stochastic collection · numerical methods

Droplet Collection

How cloud droplets become rain

A large drop falls faster than a small one and sweeps it up; that is how a cloud of tiny droplets can turn into rain within minutes. Numa Koudsie’s Bachelor thesis, Collisional Ice Multiplication Processes in Clouds (Leipzig, 2022), modelled this collection and the ice splinters it can produce, with Stokes’ law for the fall speed and the geometric cross-section for the chance of a hit. The remodel replaces each ingredient with physics that holds at the sizes that matter: a drag balance for the fall speed, collision efficiencies from the trajectory of the small droplet in the air flow around the large one, and the full stochastic collection equation for the bookkeeping. Every step is derived on the theory page and tested on the checks page, and the solver runs in your browser.

6 / 14 min
to 50 / 90 % of the mass above 40 µm, extended kernel
3 / 5 min
the same, thesis kernel
2.1 %
largest fall-speed error against measurement, from 0.3 mm up
10−7
mass error against Golovin’s exact solution

1.1The warm cloud

The test case is a warm cloud of \(300\ \mathrm{cm^{-3}}\) droplets of radius \(8\ \mathrm{\mu m}\) with \(1\ \mathrm{cm^{-3}}\) collectors of \(40\ \mathrm{\mu m}\), at \(10\,{}^{\circ}\mathrm{C}\) and \(850\ \mathrm{hPa}\). Growth is measured by the mass fraction in drops of radius \(r \ge 40\ \mathrm{\mu m}\): a diagnostic of growth, not a rainfall prediction. With the thesis kernel the mass piles up near \(500\ \mathrm{\mu m}\) and stops, because the Weber number computed with Stokes speeds switches coalescence off; with the extended kernel the drops keep growing to millimetres. The kernel matters more than the bookkeeping.

1.2What changed, one line each

  1. Fall speed. Stokes’ law gives \(\mathrm{Re}\approx 150\) for a \(200\ \mathrm{\mu m}\) drop, outside its own validity. A force balance with Beard’s drag coefficient \(C_D(\mathrm{Re})\), solved through the Davies number, reproduces the measurements of Gunn and Kinzer within 2.1 % from 0.3 mm up.
  2. Collision efficiency. The small droplet obeys \(\tau\,\dot{\mathbf v} = \mathbf u - \mathbf v - \tau g^{*}\mathbf e_z\) in the flow around the collector; scaled by \(R\) and \(U\) this is the Stokes number \(\mathrm{St} = \tau U/R\). The grazing trajectory gives \(E\); the solver reproduces the critical Stokes numbers \(1/12\) and \(1.21\).
  3. Bookkeeping. The thesis’s two fixed populations become one spectrum obeying the Smoluchowski collection equation; the thesis equations are recovered as its limit of rare collectors.
  4. Numerics. A logarithmic mass grid, pair fluxes computed in place, a deposit that conserves number and mass, and a limiter that is derived rather than guessed; tested against Golovin’s exact solution, with a mass error of \(10^{-7}\) and a spectrum error of 5 % at 8 bins per doubling.
  5. Ice. Splinter fragments join the interacting population. The app’s fragmentation is event-based: a fixed mass fraction \(\lambda\) per ice–ice collision and a fixed fragment mass. The thesis’s code drains mass at a random rate per second instead, so the cold runs test the app’s closure, not the thesis’s model. With fragments in the race the ice number grows from 167 to 1 700–2 100 cm⁻³ in an hour and the mass above 40 µm falls to 0.2–4 %: the closure, not the solver, drives that result.

1.3Four figures from the project

Fall speed against drop size for the two models and the measurements of Gunn and Kinzer, and the error of the drag balance.
Figure 1Left: the two fall-speed models against Gunn and Kinzer. Right: the error of the drag balance, at most 2.07 % at the listed sizes from 0.3 mm up and −7.57 % at 0.1 mm.
Map of the collision efficiency against the radii of the collector and the collected droplet.
Figure 2Collision efficiency for water on water at \(10\,{}^{\circ}\mathrm{C}\), \(850\ \mathrm{hPa}\), on the \(45 \times 40\) grid used by the solver.
Mass spectra of the warm cloud at several times, for the thesis kernel and the new kernel.
Figure 3Evolution of the mass spectrum in the warm cloud, stochastic collection equation, thesis kernel (left) and new kernel (right).
Golovin test: numerical spectra against the exact solution.
Figure 4Golovin test with 8 bins per doubling of mass: numerical (solid) against exact (dashed).

1.4The project’s own pages

A static copy of the Droplet Collection app. Its saved-runs list is left out: every run is computed fresh in your browser, and nothing is stored.

02

A corrected study by Sol · independently audited · 2026

rotating mechanics · energy harvesting · nonlinear dynamics

Rotor Pendula

Can a swinging weight on a turbine blade power a sensor?

Hang a 100 g weight on a 20 cm rod from a point on a turning blade. Each time the blade goes round, gravity pulls the weight first one way and then the other, so it swings once per turn, like a child on a swing pushed once per pass; a generator at the pivot brakes that swing and turns it into electricity. The study builds everything from the pendulum’s Lagrangian: the equations of motion, the energy books that say the rotor pays for every milliwatt, the driven-oscillator resonance at \(r_p = L\) with the best generator setting \(\kappa = |R-1|\), real turbine sizes from the tip-speed ratio, the double pendulum, chaos, and a sensor budget from two datasheets. An independent audit rebuilt the model from scratch and reproduced every headline number; the solver runs in your browser. Sol’s whole report is included, word for word, with its nine tables rebuilt from the data and all 181 files of the delivery to download.

57.1 mW
best harvest, single pendulum, 100 g on 20 cm
37.6 µW
what a blade sensor node needs
0.64–0.74
double ÷ single, same mass and reach
7.2 × 10−5
largest share of the turbine’s output, paid by the rotor

2.1The answer in four lines

  1. There is enough power. All 18 studied mountings give 9.55–57.1 mW of mechanical power; after an assumed 50 % conversion that is at least 127 times the 37.6 µW a low-duty sensor node draws.
  2. The rotor pays for it. At fixed rotor speed the generator’s energy comes from the rotor’s drive; the pendulum cannot raise the turbine’s net output, and costs it at most \(7.2\times10^{-5}\) of it.
  3. One rod beats two. With the same mass and the same reach, a double pendulum gives 0.64–0.74 of the single one’s power at every mounting: half its mass sits at half the reach.
  4. Chaos does not help. 2 of 54 scanned double pendula are chaotic; at matched forcing the regular orbit harvests about 0.379 W, the chaotic one about 0.24–0.25 W. The catch is force, not power: a fast rotor’s pivot sits in 18–215 g, beyond a ±8 g accelerometer.

2.2What the audit corrected

An independent model, written from scratch, re-ran all 36 mapping cases and all 54 chaos cases; every headline number reproduced, to about one part in \(10^{11}\). Eight things in the report needed correcting, and none changes a conclusion. The main ones: the turbine-size formula is \(56\ \mathrm{m/s}/\Omega\), not mm/s; the centrifugal factor \(\Gamma\) rises with rotor speed, not falls; the chaotic power’s error bar is about \(\pm0.01\ \mathrm{W}\), not \(\pm0.003\ \mathrm{W}\); one mounting puts the bob 66 mm past the blade tip; and one figure is blurred by single-precision storage. The project shows each correction where it occurs, with its evidence.

2.3Four figures from the study

The planar model: a pivot on a circle around the hub, an upper rod and a lower rod, and the generators at the two joints.
Figure 1The planar model: the pivot rides a circle round the hub; the generators feel relative joint speeds, not absolute rod speeds.
Mean harvested power at the 18 mountings against rotor speed, single and double pendulum, small and large rotors.
Figure 2Mean harvested power at the 18 mountings, single and double pendulum; lines only join the sampled points.
Left: power against damping for the matched group; right: all 54 scanned cases.
Figure 3Left: the matched group, power against damping. Right: all 54 scanned cases; the vertical line alone is not the chaos test.
Largest radial specific force at the pivot for each mounting against the accelerometer range.
Figure 4The largest radial specific force at the pivot, \(g(\Gamma+1)\), against the accelerometer’s \(\pm8\,g\) range.

2.4The project’s own pages

A static copy of the Rotor Pendula app. Its saved-runs list is left out: every run is computed fresh in your browser, and nothing is stored.