Search

In 2020, after New Zealand emerged from its short-lived lockdown and students returned to university campuses, mathematicians Jeanette McLeod and Phil Wilson from the University of Canterbury taught a final-year mathematics course using craft. A part of most lectures included a hands-on craft activity to explore a key concept and to foster curiosity. This article describes the design and delivery of the course and reflects on the overwhelmingly positive student feedback. We suggest that the creation of, and interaction with, mathematical craft objects played a significant role in shaping these responses by generating feelings of engagement, reducing stress, and allowing tangible access to mathematical ideas. [ABSTRACT FROM AUTHOR]

• Basal slip affects ice thickness but can be hard to distinguish in surface data. • Modelled ice surfaces are fit to observations using an optimal control approach through an augmented Lagrangian algorithm. • Through diffusion coefficient of a linear elliptic partial differential equation we can infer ice thickness and basal slip. A key parameter in ice flow modelling is basal slipping at the ice-bed interface as it can have a large effect on the resultant ice thickness. Unfortunately, its contribution to surface observations can be hard to distinguish from that of bed undulations. Therefore, inferring the ice thickness from surface measurements is an interesting and non-trivial inverse problem. This paper presents a method for recovering dually the ice thickness, in turn the bed topography, and the basal slip using only surface elevation and speed measurements. The unidirectional shallow ice approximation is first implemented to model steady state ice flow for given bedrock and basal slip profiles. This surface is then taken as synthetic observed data. To invert the data, an optimisation approach is taken, employing a basic Uzawa type algorithm to minimise the augmented Lagrangian objective functional. This optimisation recovers the diffusion coefficient which gives the best fit to observations. Combining this recovered diffusion with observed surface velocity, a simple Newton's method is used to recover both the ice thickness and basal slip. The method was successful in each test case and this implies that it should be possible to recover both of these parameters in two-dimensional cases also. [ABSTRACT FROM AUTHOR]

Invasive tree species pose a grave global threat to biodiversity, particularly impacting vulnerable grassland ecosystems. This paper presents a novel approach to simulate the invasion of ecosystems by these aggressive species. Our innovative coupled ODE-PDE model offers explicit invasion simulations, a significant leap beyond traditional integro-difference matrix models. With robust theoretical properties, including solution positivity, and drawing upon PDE theory, our model offers deeper insights into invasion behavior. Our model reveals that invasions undergo a prolonged quiescent phase before rapidly transitioning to a consistent rate of advancement after a critical point. This study not only sheds light on invasion dynamics but also charts a course for future exploration and potential model extensions. [ABSTRACT FROM AUTHOR]

Models of invasive species spread often assume that landscapes are spatially homogeneous; thus simplifying analysis but potentially reducing accuracy. We extend a recently developed partial differential equation model for invasive conifer spread to account for spatial heterogeneity in parameter values and introduce a method to obtain key outputs (e.g. spread rates) from computational simulations. Simulations produce patterns of spatial spread which appear qualitatively similar to observed patterns in grassland ecosystems invaded by exotic conifers, validating our spatially explicit strategy. We find that incorporating spatial variation in different parameters does not significantly affect the evolution of invasions (which are characterised by a long quiescent period followed by rapid evolution towards to a constant rate of invasion) but that distributional assumptions can have a significant impact on the spread rate of invasions. Our work demonstrates that spatial variation in site-suitability or other parameters can have a significant impact on invasions and must be considered when designing models of invasive species spread. [ABSTRACT FROM AUTHOR]

It is often argued that the indispensability of continuum models comes from their empirical adequacy despite their decoupling from the microscopic details of the modelled physical system. There is thus a commonly held misconception that temperature varying across a region of space or time can always be accurately represented as a continuous function. We discuss three inter-related cases of temperature modelling — in phase transitions, thermal boundary resistance and slip flows — and show that the continuum view is fallacious on the ground that the microscopic details of a physical system are not necessarily decoupled from continuum models. We show how temperature discontinuities are present in both data (experiments and simulations) and phenomena (theory and models) and how discontinuum models of temperature variation may have greater empirical adequacy and explanatory power. The conclusions of our paper are: a) continuum idealisations are not indispensable to modelling physical phenomena and both continuous and discontinuous representations of phenomena work depending on the context; b) temperature is not necessarily a continuously defined function in our best scientific representations of the world; and c) that its continuity, where applicable, is a contingent matter. We also raise a question as to whether discontinuous representations should be considered truly de-idealised descriptions of physical phenomena. [ABSTRACT FROM AUTHOR]

In this work, the flow conditions are numerically investigated inside uniform and non-uniform street canyons well within the atmospheric boundary layer. The numerical simulations use the steady-RANS method with the near-wall modelling approach to simulate wall roughness at the boundary. With the aim of investigating both flow structure in broad terms and pedestrian comfort in the street canyon between parallel buildings, we test different canyon configurations with varied street width, building width and building height. Turbulent conditions are broadly expected to hold within the physically realistic range of Reynolds number of order 10 6 considered here, where we take the building height to be a characteristic length scale and the free stream velocity as the characteristic velocity. In addition to discussing the features of the canyon and wake flow, we investigate the effects of canyon geometry on pedestrian comfort by using the Extended Land Beaufort Scale for this purpose. We present and compare pedestrian comfort 'maps' for each of our geometries. [ABSTRACT FROM AUTHOR]

Glacier ice flow is shaped and defined by several properties, including the bedrock elevation profile and the basal slip distribution. The effect of these two basal properties can be present in similar ways in the surface. For bedrock recovery, this makes distinguishing between them an interesting and complex problem. The results of this paper show that in some synthetic test cases it is indeed possible to distinguish and recover both bedrock elevation and basal slip given free surface elevation and free surface velocity. The unidirectional shallow ice approximation is used to compute steady-state surface data for a number of synthetic cases with different bedrock profiles and basal slip distributions. A simple inversion method based on Newton's method is applied to the known surface data to return the bedrock profile and basal slip distribution. In each synthetic test case, the inversion was successful in recovering both the bedrock elevation profile and the basal slip distribution variables. These results imply that there are a unique bedrock profile and basal slip which give rise to a unique combination of free surface velocity and free surface elevation. [ABSTRACT FROM AUTHOR]

Cortical spreading depression (CSD) is a slow-moving ionic and metabolic disturbance that propagates in cortical brain tissue. In addition to massive cellular depolarizations, CSD also involves significant changes in perfusion and metabolism—aspects of CSD that had not been modeled and are important to traumatic brain injury, subarachnoid hemorrhage, stroke, and migraine. In this study, we develop a mathematical model for CSD where we focus on modeling the features essential to understanding the implications of neurovascular coupling during CSD. In our model, the sodium-potassium–ATPase, mainly responsible for ionic homeostasis and active during CSD, operates at a rate that is dependent on the supply of oxygen. The supply of oxygen is determined by modeling blood flow through a lumped vascular tree with an effective local vessel radius that is controlled by the extracellular potassium concentration. We show that during CSD, the metabolic demands of the cortex exceed the physiological limits placed on oxygen delivery, regardless of vascular constriction or dilation. However, vasoconstriction and vasodilation play important roles in the propagation of CSD and its recovery. Our model replicates the qualitative and quantitative behavior of CSD—vasoconstriction, oxygen depletion, extracellular potassium elevation, prolonged depolarization—found in experimental studies. We predict faster, longer duration CSD in vivo than in vitro due to the contribution of the vasculature. Our results also help explain some of the variability of CSD between species and even within the same animal. These results have clinical and translational implications, as they allow for more precise in vitro, in vivo, and in silico exploration of a phenomenon broadly relevant to neurological disease. [ABSTRACT FROM AUTHOR]