Papers › Penguin huddling: a continuum model
Penguin huddling: a continuum model
Samuel J. Harris, N. R. McDonald
The archive published only this paper's code-link row. Authors, date and abstract are from arXiv's metadata (CC0), read from the Kaggle arXiv metadata snapshot of 2026-09-12 where its title matched the archive's; the title is the archive's.
Penguins huddling in a cold wind are represented by a two-dimensional, continuum model. The huddle boundary evolves due to heat loss to the huddle exterior and through the reorganisation of penguins as they seek to regulate their heat production within the huddle. These two heat transfer mechanisms, along with area, or penguin number, conservation, gives a free boundary problem whose dynamics depend on both the dynamics interior and exterior to the huddle. Assuming the huddle shape evolves slowly compared to the advective timescale of the exterior wind, the interior temperature is governed by a Poisson equation and the exterior temperature by the steady advection-diffusion equation. The exterior, advective wind velocity is the gradient of a harmonic, scalar field. The conformal invariance of the exterior governing equations is used to convert the system to a Polubarinova-Galin type equation, with forcing depending on both the interior and exterior temperature gradients at the huddle boundary. The interior Poisson equation is not conformally invariant, so the interior temperature gradient is found numerically using a combined adaptive Antoulas-Anderson and least squares algorithm. The results show that, irrespective of the starting shape, penguin huddles evolve into an egg-like steady shape. This shape is dependent on the wind strength, parameterised by the P\'eclet number Pe, and a parameter \b{eta} which effectively measures the strength of the interior self-generation of heat by the penguins. The numerical method developed is applicable to a further five free boundary problems.
Code
Repository list and official/mentioned flags are the archive's, frozen 2025-07-28. Reachability, where shown, is from one Syntology probe window (2026-09-16 to 2026-09-18); repositories not probed show nothing. GitHub stars are not tracked.
Code Syntology ran Syntology
Not run by Syntology. Nothing on this page verifies that the listed code works.
Results from the paper archive 2025-07-28
No leaderboard rows for this paper in the archive.
Report a problem or propose a change · a person checks every report against the paper or source before anything changes; decisions are listed on /corrections