Recent years have seen a large increase in the popularity of Texas hold ’em poker. It is now the most commonly played variant of the game, both in casinos and through online platforms. In this paper, we present a simulation study for games of Texas hold ’em with between two and 23 players. From these simulations, we estimate the probabilities of each player having been dealt the winning hand. These probabilities are calculated conditional on both partial information (that is, the player only having knowledge of his/her cards) and also on fuller information (that is, the true probabilities of each player winning given knowledge of the cards dealt to each player). Where possible, our estimates are compared to exact analytic results and are shown to have converged to three significant figures.

With these results, we assess the poker strategies described in two recent pieces of popular culture. In comparing the ideas expressed in Taylor Swift’s song, New Romantics, and the betting patterns employed by James Bond in the 2006 film, Casino Royale, we conclude that Ms Swift demonstrates a greater understanding of the true probabilities of winning a game of Texas hold ’em poker.

Penguins are flightless, so they are forced to walk while on land. In particular, they show rather specific behaviours in their homecoming, which are interesting to observe and to describe analytically. We observed that penguins have the tendency to waddle back and forth on the shore to create a sufficiently large group, and then walk home compactly together. The mathematical framework that we introduce describes this phenomenon, by taking into account “natural parameters”, such as the eyesight of the penguins and their cruising speed. The model that we propose favours the formation of conglomerates of penguins that gather together, but, on the other hand, it also allows the possibility of isolated and exposed individuals.

The model that we propose is based on a set of ordinary differential equations. Due to the discontinuous behaviour of the speed of the penguins, the mathematical treatment (to get existence and uniqueness of the solution) is based on a “stop-and-go” procedure. We use this setting to provide rigorous examples in which at least some penguins manage to safely return home (there are also cases in which some penguins remain isolated). To facilitate the intuition of the model, we also present some simple numerical simulations that can be compared with the actual movement of the penguin parade.

We use an integral equation formulation approach to value shout options, which are exotic options giving an investor the ability to “shout” and lock in profits while retaining the right to benefit from potentially favourable movements in the underlying asset price. Mathematically, the valuation is a free boundary problem involving an optimal exercise boundary which marks the region between shouting and not shouting. We also find the behaviour of the optimal exercise boundary for one- and two-shout options close to expiry.

At the ANZIAM conference in Hobart in February 2018, there were several talks on the solution of Laplace problems in multiply connected domains by means of conformal mapping. It appears to be not widely known that such problems can also be solved by the elementary method of series expansions with coefficients determined by least-squares fitting on the boundary. (These are not convergent series; the coefficients depend on the degree of the approximation.) Here we give a tutorial introduction to this method, which converges at an exponential rate if the boundary data are sufficiently well-behaved. The mathematical foundations go back to Runge in 1885 and Walsh in 1929. One of our examples involves an approximate Cantor set with up to 2048 components.

We analyse a scattering problem of electromagnetic waves by a bounded chiral conductive obstacle, which is surrounded by a dielectric, via the quasi-stationary approximation for the Maxwell equations. We prove the reciprocity relations for incident plane and spherical electric waves upon the scatterer. Mixed reciprocity relations have also been proved for a plane wave and a spherical wave. In the case of spherical waves, the point sources are located either inside or outside the scatterer. These relations are used to study the inverse scattering problems.

Some swimming microorganisms are sensitive to light, and this can affect the way in which they negotiate their environment. In particular, photophobic cells are repelled from unfavourable light conditions, and in a quiescent fluid environment this can be observed as elevated cell levels in regions away from these light conditions. This photoresponsive effect is of interest due to its potential technological applications. For example, the use of light to focus and direct cells could be used as a convenient means to separate out the algae used in biofuel production (for example, hydrogen), or exploited within devices for biodetection of environmental contaminants. However, in these types of situations the swimming cells will usually be suspended in a flow with shear. In this environment, it has previously been shown that cells can become hydrodynamically trapped in regions of high fluid shear, and so the extent to which photofocusing can occur under these conditions is not immediately clear. Moreover, in applications where the light must pass through appreciable volumes of the suspension, cells will typically absorb light and so shade each other from the illumination. As such, the intensity at any point in the flow is dependent upon the global cell concentration. Hence, in this study we model the coupled influence of fluid shear and cell photosensitivity on a suspension of swimming microorganisms, and ask under what circumstances a suspension of photophobic cells might be focused into high concentration regions.

Microscale propulsion is integral to numerous biomedical systems, including biofilm formation and human reproduction, where the surrounding fluids comprise suspensions of polymers. These polymers endow the fluid with non-Newtonian rheological properties, such as shear-thinning and viscoelasticity. Thus, the complex dynamics of non-Newtonian fluids present numerous modelling challenges. Here, we demonstrate that neglecting ‘out-of-plane’ effects during swimming through a shear-thinning fluid results in a significant overestimate of fluid viscosity around the undulatory swimmer Caenorhabditis elegans. This miscalculation of viscosity corresponds with an overestimate of the power the swimmer expends, a key biophysical quantity important for understanding the internal mechanics of the swimmer. As experimental flow-tracking techniques improve, accurate experimental estimates of power consumption in similar undulatory systems, such as the planar beating of human sperm through cervical mucus, will be required to probe the interaction between internal power generation, fluid rheology, and the resulting waveform.

Nonsymmetric branching flow through a three-dimensional (3D) vessel is considered at medium-to-high flow rates. The branching is from one mother vessel to two or more daughter vessels downstream, with laminar steady or unsteady conditions assumed. The inherent 3D nonsymmetry is due to the branching shapes themselves, or the differences in the end pressures in the daughter vessels, or the incident velocity profiles in the mother. Computations based on lattice-Boltzmann methodology are described first. A subsequent analysis focuses on small 3D disturbances and increased Reynolds numbers. This reduces the 3D problem to a two-dimensional one at the outer wall in all pressure-driven cases. As well as having broader implications for feeding into a network of vessels, the findings enable predictions of how much swirling motion in the cross-plane is generated in a daughter vessel downstream of a 3D branch junction, and the significant alterations provoked locally in the shear stresses and pressures at the walls. Nonuniform incident wall-shear and unsteady effects are examined. A universal asymptotic form is found for the flux change into each daughter vessel in a 3D branching of arbitrary cross-section with a thin divider.

Heterogeneity in pulmonary microvascular blood flow (perfusion) provides an early indicator of lung disease or disease susceptibility. However, most computational models of the pulmonary vasculature neglect structural heterogeneities, and are thus not accurate predictors of lung function in disease that is not diffuse (spread evenly through the lung). Models that do incorporate structural heterogeneity have either neglected the temporal dynamics of blood flow, or the structure of the smallest blood vessels. Larger than normal oscillations in pulmonary capillary calibre, high oscillatory stress contribute to disease progression. Hence, a model that captures both spatial and temporal heterogeneity in pulmonary perfusion could provide new insights into the early stages of pulmonary vascular disease. Here, we present a model of the pulmonary vasculature, which captures both flow dynamics, and the anatomic structure of the pulmonary blood vessels from the right to left heart including the micro-vasculature. The model is compared to experimental data in normal lungs. We confirm that spatial heterogeneity in pulmonary perfusion is time-dependent, and predict key features of pulmonary hypertensive disease using a simple implementation of increased vascular stiffness.