Models of collective cell motion for cell populations with different aspect ratio : diffusion, proliferation and travelling waves

Continuum, partial differential equation models are often used to describe the collective motion of cell populations, with various types of motility represented by the choice of diffusion coefficient, and cell proliferation captured by the source terms. Previously, the choice of diffusion coefficient has been largely arbitrary, with the decision to choose a particular linear or nonlinear form generally based on calibration arguments rather than making any physical connection with the underlying individual-level properties of the cell motility mechanism. In this work we provide a new link between individual-level models, which account for important cell properties such as varying cell shape and volume exclusion, and population-level partial differential equation models. We work in an exclusion process framework, considering aligned, elongated cells that may occupy more than one lattice site, in order to represent populations of agents with different sizes. Three different idealizations of the individual-level mechanism are proposed, and these are connected to three different partial differential equations, each with a different diffusion coefficient; one linear, one nonlinear and degenerate and one nonlinear and nondegenerate. We test the ability of these three models to predict the population level response of a cell spreading problem for both proliferative and nonproliferative cases. We also explore the potential of our models to predict long time travelling wave invasion rates and extend our results to two dimensional spreading and invasion. Our results show that each model can accurately predict density data for nonproliferative systems, but that only one does so for proliferative systems. Hence great care must be taken to predict density data for with varying cell shape.

Numerical method for predicting the elastic lateral distortional buckling moment of a mono-symmetric beam with web openings

When used as floor joists, the new mono-symmetric LiteSteel beam (LSB) sections require web openings to provide access for inspections and various services. The LSBs consist of two rectangular hollow flanges connected by a slender web, and are subjected to lateral distortional buckling effects in the intermediate span range. Their member capacity design formulae developed to date are based on their elastic lateral buckling moments, and only limited research has been undertaken to predict the elastic lateral buckling moments of LSBs with web openings. This paper addresses this research gap by reporting the development of web opening modelling techniques based on an equivalent reduced web thickness concept and a numerical method for predicting the elastic buckling moments of LSBs with circular web openings. The proposed numerical method was based on a formulation of the total potential energy of LSBs with circular web openings. The accuracy of the proposed method’s use with the aforementioned modelling techniques was verified through comparison of its results with those of finite strip and finite element analyses of various LSBs.

Elastic lateral buckling of cantilever LiteSteel beams under transverse loading

The LiteSteel Beam (LSB) is a new hollow flange channel section developed by OneSteel Australian Tube Mills using its patented dual electric resistance welding and automated continuous roll-forming technologies. The LSB has a unique geometry consisting of torsionally rigid rectangular hollow flanges and a relatively slender web. Its flexural strength for intermediate spans is governed by lateral distortional buckling characterised by simultaneous lateral deflection, twist and web distortion. Recent research on LSBs has mainly focussed on their lateral distortional buckling behaviour under uniform moment conditions. However, in practice, LSB flexural members are subjected to non-uniform moment distributions and load height effects as they are often under transverse loads applied above or below their shear centre. These loading conditions are known to have significant effects on the lateral buckling strength of beams. Many steel design codes have adopted equivalent uniform moment distribution and load height factors based on data for conventional hot-rolled, doubly symmetric I-beams subject to lateral torsional buckling. The non-uniform moment distribution and load height effects of transverse loading on cantilever LSBs, and the suitability of the current design modification factors to include such effects are not known. This paper presents a numerical study based on finite element analyses of the elastic lateral buckling strength of cantilever LSBs subject to transverse loading, and the results. The applicability of the design modification factors from various steel design codes was reviewed, and suitable recommendations are presented for cantilever LSBs subject to transverse loading.

Cold and heat waves in the United States

Extreme cold and heat waves, characterised by a number of cold or hot days in succession, place a strain on people’s cardiovascular and respiratory systems. The increase in deaths due to these waves may be greater than that predicted by extreme temperatures alone. We examined cold and heat waves in 99 US cities for 14 years (1987–2000) and investigated how the risk of death depended on the temperature threshold used to define a wave, and a wave’s timing, duration and intensity. We defined cold and heat waves using temperatures above and below cold and heat thresholds for two or more days. We tried five cold thresholds using the first to fifth percentiles of temperature, and five heat thresholds using the ninety-fifth to ninety-ninth percentiles. The extra wave effects were estimated using a two-stage model to ensure that their effects were estimated after removing the general effects of temperature. The increases in deaths associated with cold waves were generally small and not statistically significant, and there was even evidence of a decreased risk during the coldest waves. Heat waves generally increased the risk of death, particularly for the hottest heat threshold. Cold waves of a colder intensity or longer duration were not more dangerous. Cold waves earlier in the cool season were more dangerous, as were heat waves earlier in the warm season. In general there was no increased risk of death during cold waves above the known increased risk associated with cold temperatures. Cold or heat waves earlier in the cool or warm season may be more dangerous because of a build up in the susceptible pool or a lack of preparedness for cold or hot temperatures.

Tactically-driven nonmonotone travelling waves

Models of cell invasion incorporating directed cell movement up a gradient of an external substance and carrying capacity-limited proliferation give rise to travelling wave solutions. Travelling wave profiles with various shapes, including smooth monotonically decreasing, shock-fronted monotonically decreasing and shock-fronted nonmonotone shapes, have been reported previously in the literature. The existence of tacticallydriven shock-fronted nonmonotone travelling wave solutions is analysed for the first time. We develop a necessary condition for nonmonotone shock-fronted solutions. This condition shows that some of the previously reported shock-fronted nonmonotone solutions are genuine while others are a consequence of numerical error. Our results demonstrate that, for certain conditions, travelling wave solutions can be either smooth and monotone, smooth and nonmonotone or discontinuous and nonmonotone. These different shapes correspond to different invasion speeds. A necessary and sufficient condition for the travelling wave with minimum wave speed to be nonmonotone is presented. Several common forms of the tactic sensitivity function have the potential to satisfy the newly developed condition for nonmonotone shock-fronted solutions developed in this work.

Beat phenomena in metal nanowires, and their implications for resonance-based elastic property measurements

The elastic properties of 1D nanostructures such as nanowires are often measured experimentally through actuation of the nanowire at its resonance frequency, and then relating the resonance frequency to the elastic stiffness using elementary beam theory. In the present work, we utilize large scale molecular dynamics simulations to report a novel beat phenomenon in [110]oriented Ag nanowires. The beat phenomenon is found to arise from the asymmetry of the lattice spacing in the orthogonal elementary directions of the [110] nanowire, i.e. the [-110] and [001] directions, which results in two different principal moments of inertia. Because of this, actuations imposed along any other direction are found to decompose into two orthogonal vibrational components based on the actuation angle relative to these two elementary directions, with this phenomenon being generalizable to <110> FCC nanowires of different materials (Cu, Au, Ni, Pd and Pt). The beat phenomenon is explained using a discrete moment of inertia model based on the hard sphere assumption, the model is utilized to show that surface effects enhance the beat phenomenon, while the effect is reduced with increasing nanowires cross-sectional size or aspect ratio. Most importantly, due to the existence of the beat phenomena, we demonstrate that in resonance experiments only a single frequency component is expected to be observed, particularly when the damping ratio is relatively large or very small. Furthermore, for a large range of actuation angles, the lower frequency is more likely to be detected than the higher one, which implies that experimental predictions of Young’s modulus obtained from resonance may in fact be under predictions. The present study therefore has significant implications for experimental interpretations of Young’s modulus as obtained via resonance testing.

Locally weighted learning model predictivek control for elastic joint robots

This paper proposes an efficient and online learning control system that uses the successful Model Predictive Control (MPC) method in a model based locally weighted learning framework. The new approach named Locally Weighted Learning Model Predictive Control (LWL-MPC) has been proposed as a solution to learn to control complex and nonlinear Elastic Joint Robots (EJR). Elastic Joint Robots are generally difficult to learn to control due to their elastic properties preventing standard model learning techniques from being used, such as learning computed torque control. This paper demonstrates the capability of LWL-MPC to perform online and incremental learning while controlling the joint positions of a real three Degree of Freedom (DoF) EJR. An experiment on a real EJR is presented and LWL-MPC is shown to successfully learn to control the system to follow two different figure of eight trajectories.

Travelling waves for a velocity–jump model of cell migration and proliferation

Cell invasion, characterised by moving fronts of cells, is an essential aspect of development, repair and disease. Typically, mathematical models of cell invasion are based on the Fisher–Kolmogorov equation. These traditional parabolic models can not be used to represent experimental measurements of individual cell velocities within the invading population since they imply that information propagates with infinite speed. To overcome this limitation we study combined cell motility and proliferation based on a velocity–jump process where information propagates with finite speed. The model treats the total population of cells as two interacting subpopulations: a subpopulation of left–moving cells, $L(x,t)$, and a subpopulation of right–moving cells, $R(x,t)$. This leads to a system of hyperbolic partial differential equations that includes a turning rate, $\Lambda \ge 0$, describing the rate at which individuals in the population change direction of movement. We present exact travelling wave solutions of the system of partial differential equations for the special case where $\Lambda = 0$ and in the limit that $\Lambda \to \infty$. For intermediate turning rates, $0 < \Lambda < \infty$, we analyse the travelling waves using the phase plane and we demonstrate a transition from smooth monotone travelling waves to smooth nonmonotone travelling waves as $\Lambda$ decreases through a critical value $\Lambda_{crit}$. We conclude by providing a qualitative comparison between the travelling wave solutions of our model and experimental observations of cell invasion. This comparison indicates that the small $\Lambda$ limit produces results that are consistent with experimental observations.

Thermal-elastic stresses and the criticality of the continental crust

Heating or cooling can lead to high stresses in rocks due to the different thermal-elastic properties of minerals. In the upper 4 km of the crust, such internal stresses might cause fracturing. Yet it is unclear if thermal elasticity contributes significantly to critical stresses and failure deeper in Earth's continental crust, where ductile creep causes stress relaxation. We combined a heating experiment conducted in a Synchrotron microtomograph (Advanced Photon Source, USA) with numerical simulations to calculate the grain-scale stress field in granite generated by slow burial. We find that deviatoric stresses >100 MPa can be stored during burial, with relaxation times from 100's to 1000's ka, even in the ductile crust. Hence, grain-scale thermal-elastic stresses may serve as nuclei for instabilities, thus rendering the continental crust close to criticality.

Comparison between analytical and 3D finite element solutions for borehole stresses in anisotropic elastic rock

We present a rigorous validation of the analytical Amadei solution for the stress concentration around an arbitrarily orientated borehole in general anisotropic elastic media. First, we revisit the theoretical framework of the Amadei solution and present analytical insights that show that the solution does indeed contain all special cases of symmetry, contrary to previous understanding, provided that the reduced strain coefficients b11 and b55 are not equal. It is shown from theoretical considerations and published experimental data that the b11 and b55 are not equal for realistic rocks. Second, we develop a 3D finite element elastic model within a hybrid analytical–numerical workflow that circumvents the need to rebuild and remesh the model for every borehole and material orientation. Third, we show that the borehole stresses computed from the numerical model and the analytical solution match almost perfectly for different borehole orientations (vertical, deviated and horizontal) and for several cases involving isotropic, transverse isotropic and orthorhombic symmetries. It is concluded that the analytical Amadei solution is valid with no restriction on the borehole orientation or the symmetry of the elastic anisotropy.

Lattice-free models of cell invasion : discrete simulations and travelling waves

Invasion waves of cells play an important role in development, disease and repair. Standard discrete models of such processes typically involve simulating cell motility, cell proliferation and cell-to-cell crowding effects in a lattice-based framework. The continuum-limit description is often given by a reaction–diffusion equation that is related to the Fisher–Kolmogorov equation. One of the limitations of a standard lattice-based approach is that real cells move and proliferate in continuous space and are not restricted to a predefined lattice structure. We present a lattice-free model of cell motility and proliferation, with cell-to-cell crowding effects, and we use the model to replicate invasion wave-type behaviour. The continuum-limit description of the discrete model is a reaction–diffusion equation with a proliferation term that is different from lattice-based models. Comparing lattice based and lattice-free simulations indicates that both models lead to invasion fronts that are similar at the leading edge, where the cell density is low. Conversely, the two models make different predictions in the high density region of the domain, well behind the leading edge. We analyse the continuum-limit description of the lattice based and lattice-free models to show that both give rise to invasion wave type solutions that move with the same speed but have very different shapes. We explore the significance of these differences by calibrating the parameters in the standard Fisher–Kolmogorov equation using data from the lattice-free model. We conclude that estimating parameters using this kind of standard procedure can produce misleading results.

A Framework for Adaptation of Australian Households to Heat Waves

Climate change is leading to an increased frequency and severity of heat waves. Spells of several consecutive days of unusually high temperatures have led to increased mortality rates for the more vulnerable in the community. The problem is compounded by the escalating energy costs and increasing peak electrical demand as people become more reliant on air conditioning. Domestic air conditioning is the primary determinant of peak power demand which has been a major driver of higher electricity costs. This report presents the findings of multidisciplinary research which develops a national framework to evaluate the potential impacts of heat waves. It presents a technical, social and economic approach to adapt Australian residential buildings to ameliorate the impact of heat waves in the community and reduce the risk of its adverse outcomes. Through the development of a methodology for estimating the impact of global warming on key weather parameters in 2030 and 2050, it is possible to re-evaluate the size and anticipated energy consumption of air conditioners in future years for various climate zones in Australia. Over the coming decades it is likely that mainland Australia will require more cooling than heating. While in some parts the total electricity usage for heating and cooling may remain unchanged, there is an overall significant increase in peak electricity demand, likely to further drive electricity prices. Through monitoring groups of households in South Australia, New South Wales and Queensland, the impact of heat waves on both thermal comfort sensation and energy consumption for air conditioning has been evaluated. The results show that households are likely to be able to tolerate slightly increased temperature levels indoors during periods of high outside temperatures. The research identified that household electricity costs are likely to rise above what is currently projected due to the impact of climate change. Through a number of regulatory changes to both household design and air conditioners, this impact can be minimised. A number of proposed retrofit and design measures are provided, which can readily reduce electricity usage for cooling at minimal cost to the household. Using a number of social research instruments, it is evident that households are willing to change behaviour rather than to spend money. Those on lower income and elderly individuals are the least able to afford the use of air conditioning and should be a priority for interventions and assistance. Increasing community awareness of cost effective strategies to manage comfort and health during heat waves is a high priority recommended action. Overall, the research showed that a combined approach including behaviour change, dwelling modification and improved air conditioner selection can readily adapt Australian households to the impact of heat waves, reducing the risk of heat related deaths and household energy costs.

Elastic buckling analysis of ideal thin-walled structures under combined loading using a finite strip method

The details of an application of the finite strip method to the elastic buckling analysis of thin-walled structures with various boundary conditions and subjected to single or combined loadings of longitudinal compression, transverse compression, bending and shear are presented. The presence of shear loading is accounted for by modifying the displacement functions which are commonly used in cases when shear is absent. A program based on the finite strip method was used to obtain the elastic buckling stress, buckling plot and buckling mode of thin-walled structures and some of these results are presented.

Waves for change : the role of the South Tarawa-based women's interests program in the decolonisation process of the Gilbert Islands

This thesis has created a space for women in the history of the decolonisation of the Gilbert Islands. It traces the historical development of the national women's interests program in the Republic of Kiribati (formerly of the Gilbert and Ellice Islands Colony (GEIC)) as it was implemented through a network of women's clubs during the 1960s and 1970s. This thesis has provided the first history and interpretation of the Indigenous women's interests movement as it impacted the Gilbert Islands. It offers a narrative of the movement in terms of three overlapping waves of women leaders, based on an analysis of fieldwork, archival research and interviews conducted on South Tarawa, Kiribati.

Higher-order non-linear analysis of steel structures. Part I : elastic second-order formulation

This paper presents a higher-order beam-column formulation that can capture the geometrically non-linear behaviour of steel framed structures which contain a multiplicity of slender members. Despite advances in computational frame software, analyses of large frames can still be problematic from a numerical standpoint and so the intent of the paper is to fulfil a need for versatile, reliable and efficient non-linear analysis of general steel framed structures with very many members. Following a comprehensive review of numerical frame analysis techniques, a fourth-order element is derived and implemented in an updated Lagrangian formulation, and it is able to predict flexural buckling, snap-through buckling and large displacement post-buckling behaviour of typical structures whose responses have been reported by independent researchers. The solutions are shown to be efficacious in terms of a balance of accuracy and computational expediency. The higher-order element forms a basis for augmenting the geometrically non-linear approach with material non-linearity through the refined plastic hinge methodology described in the companion paper.