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.

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.

A model of workplace safety incorporating worker interactions and simple interventions

Although there was substantial research into the occupational health and safety sector over the past forty years, this generally focused on statistical analyses of data related to costs and/or fatalities and injuries. There is a lack of mathematical modelling of the interactions between workers and the resulting safety dynamics of the workplace. There is also little work investigating the potential impact of different safety intervention programs prior to their implementation. In this article, we present a fundamental, differential equation-based model of workplace safety that treats worker safety habits similarly to an infectious disease in an epidemic model. Analytical results for the model, derived via phase plane and stability analysis, are discussed. The model is coupled with a model of a generic safety strategy aimed at minimising unsafe work habits, to produce an optimal control problem. The optimal control model is solved using the forward-backward sweep numerical scheme implemented in Matlab.

Novel solutions for a model of wound healing angiogenesis

We prove the existence of novel, shock-fronted travelling wave solutions to a model of wound healing angiogenesis studied in Pettet et al (2000 IMA J. Math. App. Med. 17 395–413) assuming two conjectures hold. In the previous work, the authors showed that for certain parameter values, a heteroclinic orbit in the phase plane representing a smooth travelling wave solution exists. However, upon varying one of the parameters, the heteroclinic orbit was destroyed, or rather cut-off, by a wall of singularities in the phase plane. As a result, they concluded that under this parameter regime no travelling wave solutions existed. Using techniques from geometric singular perturbation theory and canard theory, we show that a travelling wave solution actually still exists for this parameter regime. We construct a heteroclinic orbit passing through the wall of singularities via a folded saddle canard point onto a repelling slow manifold. The orbit leaves this manifold via the fast dynamics and lands on the attracting slow manifold, finally connecting to its end state. This new travelling wave is no longer smooth but exhibits a sharp front or shock. Finally, we identify regions in parameter space where we expect that similar solutions exist. Moreover, we discuss the possibility of more exotic solutions.

Footwear modifies coronal plane forefoot and sagittal plane hallux kinematics during stance phase of walking gait

Footwear is designed to reduce injury, and enhance performance. However, the effect footwear has on foot and ankle kinematics currently remains unknown. Acknowledging the need for improved understanding, multi-segment models of the foot-shoe complex need to be established to both describe and quantify the effect footwear has on the foot and ankle during stance phase of gait. The purpose of this study was to quantify how footwear alters the kinematics of the foot inside the shoe during stance phase of walking gait.

Simulation of inplane and out of plane AE source in thin plate

In most materials, short stress waves are generated during the process of plastic deformation, phase transformation, crack formation and crack growth. These phenomena are applied in acoustic emission (AE) for the detection of material defects in a wide spectrum of areas, ranging from nondestructive testing for the detection of materials defects to monitoring of microseismical activity. AE technique is also used for defect source identification and for failure detection. AE waves consist of P waves (primary longitudinal waves), S waves (shear/transverse waves) and Rayleigh (surface) waves as well as reflected and diffracted waves. The propagation of AE waves in various modes has made the determination of source location difficult. In order to use acoustic emission technique for accurate identification of source, an understanding of wave propagation of the AE signals at various locations in a plate structure is essential. Furthermore, an understanding of wave propagation can also assist in sensor location for optimum detection of AE signals along with the characteristics of the source. In real life, as the AE signals radiate from the source it will result in stress waves. Unless the type of stress wave is known, it is very difficult to locate the source when using the classical propagation velocity equations. This paper describes the simulation of AE waves to identify the source location and its characteristics in steel plate as well as the wave modes. The finite element analysis (FEA) is used for the numerical simulation of wave propagation in thin plate. By knowing the type of wave generated, it is possible to apply the appropriate wave equations to determine the location of the source. For a single plate structure, the results show that the simulation algorithm is effective to simulate different stress waves.

The persistence of phase-separation in LiFePO4 with two-dimensional Li+ transport : the Cahn-Hilliard-reaction equation and the role of defects

We examine the solution of the two-dimensional Cahn-Hilliard-reaction (CHR) equation in the xy plane as a model of Li+ intercalation into LiFePO4 material. We validate our numerical solution against the solution of the depth-averaged equation, which has been used to model intercalation in the limit of highly orthotropic diffusivity and gradient penalty tensors. We then examine the phase-change behaviour in the full CHR system as these parameters become more isotropic, and find that as the Li+ diffusivity is increased in the x direction, phase separation persists at high currents, even in small crystals with averaged coherency strain included. The resulting voltage curves decrease monotonically, which has previously been considered a hallmark of crystals that fill homogeneously.

A phase-locked-loop design for the smooth operation of a hybrid microgrid

A microgrid contains both distributed generators (DGs) and loads and can be viewed by a controllable load by utilities. The DGs can be either inertial synchronous generators or non-inertial converter interfaced. Moreover, some of them can come online or go offline in plug and play fashion. The combination of these various types of operation makes the microgrid control a challenging task, especially when the microgrid operates in an autonomous mode. In this paper, a new phase locked loop (PLL) algorithm is proposed for smooth synchronization of plug and play DGs. A frequency droop for power sharing is used and a pseudo inertia has been introduced to non-inertial DGs in order to match their response with inertial DGs. The proposed strategy is validated through PSCAD simulation studies.

Integrated decoupling and matching network incorporated in the ground plane of a compact monopole array

Compact arrays enable various applications such as antenna beam-forming and multi-input, multi-output (MIMO) schemes on limited-size platforms. The reduced element spacing in compact arrays introduces high levels of mutual coupling which can affect the performance of the adaptive array. This coupling causes a mismatch at the input ports, which disturbs the performance of the individual elements in the array and affects the implementation of beam steering. In this article, a reactive decoupling network for a 3-element monopole array is used to establish port isolation while simultaneously matching input impedance at each port to the system impendence. The integrated decoupling and matching network is incorporated in the ground plane of the monopole array, providing further development scope for beamforming using phase shifters and power splitters in double-layered circuits.

Condensed phase combustion travelling waves with sequential exothermic or endothermic reactions

The one-dimensional propagation of a combustion wave through a premixed solid fuel for two-stage kinetics is studied. We re-examine the analysis of a single reaction travelling-wave and extend it to the case of two-stage reactions. We derive an expression for the travelling wave speed in the limit of large activation energy for both reactions. The analysis shows that when both reactions are exothermic, the wave structure is similar to the single reaction case. However, when the second reaction is endothermic, the wave structure can be significantly different from single reaction case. In particular, as might be expected, a travelling wave does not necessarily exist in this case. We establish conditions in the limiting large activation energy limit for the non-existence, and for monotonicity of the temperature profile in the travelling wave.

The Behavior of Hydroxyl Units of Synthetic Goethite and its Dehydroxylated Product Hematite

The behavior of the hydroxyl units of synthetic goethite and its dehydroxylated product hematite was characterized using a combination of Fourier transform infrared (FTIR) spectroscopy and X-ray diffraction (XRD) during the thermal transformation over a temperature range of 180-270 degrees C. Hematite was detected at temperatures above 200 degrees C by XRD while goethite was not observed above 230 degrees C. Five intense OH vibrations at 3212-3194, 1687-1674, 1643-1640, 888-884 and 800-798 cm(-1), and a H2O vibration at 3450-3445 cm(-1) were observed for goethite. The intensity of hydroxyl stretching and bending vibrations decreased with the extent of dehydroxylation of goethite. Infrared absorption bands clearly show the phase transformation between goethite and hematite: in particular. the migration of excess hydroxyl units from goethite to hematite. Two bands at 536-533 and 454-452 cm(-1) are the low wavenumber vibrations of Fe-O in the hematite structure. Band component analysis data of FTIR spectra support the fact that the hydroxyl units mainly affect the a plane in goethite and the equivalent c plane in hematite.