We use cookies to distinguish you from other users and to provide you with a better experience on our websites. Close this message to accept cookies or find out how to manage your cookie settings.
To save content items to your account,
please confirm that you agree to abide by our usage policies.
If this is the first time you use this feature, you will be asked to authorise Cambridge Core to connect with your account.
Find out more about saving content to .
To save content items to your Kindle, first ensure no-reply@cambridge.org
is added to your Approved Personal Document E-mail List under your Personal Document Settings
on the Manage Your Content and Devices page of your Amazon account. Then enter the ‘name’ part
of your Kindle email address below.
Find out more about saving to your Kindle.
Note you can select to save to either the @free.kindle.com or @kindle.com variations.
‘@free.kindle.com’ emails are free but can only be saved to your device when it is connected to wi-fi.
‘@kindle.com’ emails can be delivered even when you are not connected to wi-fi, but note that service fees apply.
Thin spray-on liners (TSLs) have been found to be effective for structurally supporting the walls of mining tunnels and thus reducing the occurrence of rock bursts, an effect primarily due to the penetration of cracks by the liner. Surface tension effects are thus important. However, TSLs are also used to simply stabilize rock surfaces, for example, to prevent rock fall, and in this context crack penetration is desirable but not necessary, and the tensile and shearing strength and adhesive properties of the liner determine its effectiveness. We examine the effectiveness of nonpenetrating TSLs in a global lined tunnel and in a local rock support context. In the tunnel context, we examine the effect of the liner on the stress distribution in a tunnel subjected to a geological or mining event. We show that the liner has little effect on stresses in the surrounding rock and that tensile stresses in the rock surface are transmitted across the liner, so that failure is likely to be due to liner rupture or detachment from the surface. In the local rock support context, loose rock movements are shown to be better achieved using a liner with small Young’s modulus, but high rupture strength.
The effects of apparatus-induced dispersion on nonuniform, density-dependent flow in a cylindrical soil column were investigated using a finite-element model. To validate the model, the results with an analytical solution and laboratory column test data were analysed. The model simulations confirmed that flow nonuniformities induced by the apparatus are dissipated within the column when the distance to the apparatus outlet exceeds $3R/2$, where R represents the radius of the cylindrical column. Furthermore, the simulations revealed that convergent flow in the vicinity of the outlet introduces additional hydrodynamic dispersion in the soil column apparatus. However, this effect is minimal in the region where the column height exceeds $3R/2$. Additionally, it is found that an increase in the solution density gradient during the solute breakthrough period led to a decrease in flow velocity, which stabilized the flow and ultimately reduced dispersive mixing. Overall, this study provides insights into the behaviour of apparatus-induced dispersion in nonuniform, density-dependent flow within a cylindrical soil column, shedding light on the dynamics and mitigation of flow nonuniformities and dispersive mixing phenomena.
In this paper, the pricing of equity warrants under a class of fractional Brownian motion models is investigated numerically. By establishing a new nonlinear partial differential equation (PDE) system governing the price in terms of the observable stock price, we solve the pricing system effectively by a robust implicit-explicit numerical method. This is fundamentally different from the documented methods, which first solve the price with respect to the firm value analytically, by assuming that the volatility of the firm is constant, and then compute the price with respect to the stock price and estimate the firm volatility numerically. It is shown that the proposed method is stable in the maximum-norm sense. Furthermore, a sharp theoretical error estimate for the current method is provided, which is also verified numerically. Numerical examples suggest that the current method is efficient and can produce results that are, overall, closer to real market prices than other existing approaches. A great advantage of the current method is that it can be extended easily to price equity warrants under other complicated models.
We consider a generalization of the well-known nonlinear Nicholson blowflies model with stochastic perturbations. Stability in probability of the positive equilibrium of the considered equation is studied. Two types of stability conditions: delay-dependent and delay-independent conditions are obtained, using the method of Lyapunov functionals and the method of linear matrix inequalities. The obtained results are illustrated by numerical simulations by means of some examples. The results are new, and complement the existing ones.
Three typical elastic problems, including beam bending, truss extension and compression, and two-rings collision are simulated with smoothed particle hydrodynamics (SPH) using Lagrangian and Eulerian algorithms. A contact-force model for elastic collisions and equation of state for pressure arising in colliding elastic bodies are also analytically derived. Numerical validations, on using the corresponding theoretical models, are carried out for the beam bending, truss extension and compression simulations. Numerical instabilities caused by largely deformed particle configurations in finite/large elastic deformations are analysed. The numerical experiments show that the algorithms handle small deformations well, but only the Lagrangian algorithm can handle large elastic deformations. The numerical results obtained from the Lagrangian algorithm also show a good agreement with the theoretical values.
We provide an analytic solution of the Rössler equations based on the asymptotic limit $c\to \infty $ and we show in this limit that the solution takes the form of multiple pulses, similar to “burst” firing of neurons. We are able to derive an approximate Poincaré map for the solutions, which compares reasonably with a numerically derived map.
Mathematical modelling has been used to support the response to the COVID-19 pandemic in countries around the world including Australia and New Zealand. Both these countries have followed similar pandemic response strategies, using a combination of strict border measures and community interventions to minimize infection rates until high vaccine coverage was achieved. This required a different set of modelling tools to those used in countries that experienced much higher levels of prevalence throughout the pandemic.
In this article, we provide an overview of some of the mathematical modelling and data analytics work that has helped to inform the policy response to the pandemic in Australia and New Zealand. This is a reflection on our experiences working at the modelling–policy interface and the impact this has had on the pandemic response. We outline the various types of model outputs, from short-term forecasts to longer-term scenario models, that have been used in different contexts. We discuss issues relating to communication between mathematical modellers and stakeholders such as health officials and policymakers. We conclude with some future challenges and opportunities in this area.
Six patents were secured by E. H. Lanier from 1930 to 1933 for aeroplane designs that were intended to be exceptionally stable. A feature of five of these was a flow-induced “vacuum chamber” which was thought to provide superior stability and increased lift compared to typical wing designs. Initially, this chamber was in the fuselage, but later designs placed it in the wing by replacing a section of the upper skin of the wing with a series of angled slats. We report upon an investigation of the Lanier wing design using inviscid aerodynamic theory and viscous numerical simulations. This took place at the 2005 Australia–New Zealand Mathematics-in-Industry Study Group. The evidence from this investigation does not support the claims but, rather, suggests that any improvement in lift and/or stability seen in the few prototypes that were built was, most probably, due to thicker airfoils than were typical at the time.
Double-masking may be used to reduce the transmission of a virus. If additionally the masks are compressible, with different permeabilities and behaviour under compression, then it may be possible to design a mask that allows for easy breathing under normal breathing conditions, but is relatively impermeable under coughing or sneezing conditions. Such a mask could be both comfortable to wear and effective. We obtain analytical solutions for the steady-state flow-through behaviour of such a double mask under flow-out conditions. The results show that the reduction in permeability required to produce a relatively impermeable mask under high flux expulsion (sneezing) conditions could be achieved using either a single filter compressible mask or two filters with different poroelastic parameters. The parameters can be more easily adjusted using a double mask. For both single- and double-mask cases, there is an abrupt cut off, whereby through-flux levels reduce from a maximum value to zero as pressure drop levels increase beyond a critical value. Additionally, in the double-mask case, there exists a second steady-state solution for particular parameter ranges. This second solution is unlikely to occur under normal circumstances.
Exact solutions are constructed for a class of nonlinear hyperbolic reaction-diffusion equations in two-space dimensions. Reduction of variables and subsequent solutions follow from a special nonclassical symmetry that uncovers a conditionally integrable system, equivalent to the linear Helmholtz equation. The hyperbolicity is commonly associated with a speed limit due to a delay, $\tau $, between gradients and fluxes. With lethal boundary conditions on a circular domain wherein a species population exhibits logistic growth of Fisher–KPP type with equal time lag, the critical domain size for avoidance of extinction does not depend on $\tau $. A diminishing exact solution within a circular domain is also constructed, when the reaction represents a weak Allee effect of Huxley type. For a combustion reaction of Arrhenius type, the only known exact solution that is finite but unbounded is extended to allow for a positive $\tau $.
To compute the maximum speed threshold for helicopters, we model the lift produced by the rotor blades. Using this model, we derive limits for each method traditionally used to alleviate dissymmetry of lift. Additionally, we find the minimum rotor angular velocity required to produce a prescribed lift at a given forward velocity. We derive conditions on the coefficient of lift for helicopter airfoils that maintain altitude. Further considerations are also made with regard to the properties of the air and its effect on helicopter dynamics.
Edited by
Alik Ismail-Zadeh, Karlsruhe Institute of Technology, Germany,Fabio Castelli, Università degli Studi, Florence,Dylan Jones, University of Toronto,Sabrina Sanchez, Max Planck Institute for Solar System Research, Germany
Abstract: In this chapter, I present an overview of waveform tomography, in the context of imaging of the Earth‘s whole mantle at the global scale. In this context, waveform tomography is defined utilising entire wide-band filtered records of the seismic wavefield, generated by natural earthquakes and observed at broadband receivers located at teleseismic distances. This is in contrast to imaging methodologies that first extract secondary observables, such as, most commonly, travel times of the most prominent energy arrivals (i.e. seismic phases), that can be easily identified and isolated in the records. Waveform tomography is a non-linear process that requires the ability to compute the predicted wavefield in a given three-dimensional Earth model and compare it to the observed wavefield. One of its main challenges, is the computational cost involved. I first review the history of methodological developments, specifically focusing on the global, whole mantle Earth imaging problem. I then discuss and contrast the two recent methodologies that have led to the development of the first three-dimensional elastic global shear velocity models that have been published to-date using numerical integration of the wave equation, specifically, using the spectral element method. I discuss how the forward problem is addressed, the data selection approaches, definitions of the misfit function, and computation of kernels for the inverse step of the imaging procedure, as well as the choice of the optimisation method. I also discuss model parametrisation, and, in particular, the important topic of how the strongly heterogeneous crust is modelled. In the final parts of this chapter, I discuss efforts towards resolving the difficult problem of model evaluation and present my views on promising directions and remaining challenges in this rapidly evolving field, aiming at further improving resolution of deep mantle elastic structure with the goal of informing our understanding of the dynamics of our planet.
Edited by
Alik Ismail-Zadeh, Karlsruhe Institute of Technology, Germany,Fabio Castelli, Università degli Studi, Florence,Dylan Jones, University of Toronto,Sabrina Sanchez, Max Planck Institute for Solar System Research, Germany
Abstract: Hydrological sciences cover a wide variety of water-driven processes at the Earth’s surface, above, and below it. Data assimilation techniques in hydrology have developed over the years along many quite independent paths, following not only different data availabilities but also a plethora of problem-specific model structures. Most hydrologic problems that are addressed through data assimilation, however, share some distinct peculiarities: scarce or indirect observation of most important state variables (soil moisture, river discharge, groundwater level, to name a few), incomplete or conceptual modelling, extreme spatial heterogeneity, and uncertainty of controlling physical parameters. On the other side, adoption of simplified and scale-specific models allows for substantial problem reduction that partially compensates these difficulties, opening the path to the assimilation of very indirect observations (e.g. from satellite remote sensing) and efficient model inversion for parameter estimation. This chapter illustrates the peculiarities of data assimilation for state estimation and model inversion in hydrology, with reference to a number of representative applications. Sequential ensemble filters and variational methods are recognised to be the most common choices in hydrologic data assimilation, and the motivations for these choices are also discussed, with several examples.
Edited by
Alik Ismail-Zadeh, Karlsruhe Institute of Technology, Germany,Fabio Castelli, Università degli Studi, Florence,Dylan Jones, University of Toronto,Sabrina Sanchez, Max Planck Institute for Solar System Research, Germany
Abstract: Although many geophysical observations and models are available for the Tibetan Plateau (TP) and its surroundings regions, our knowledge and understanding of the uplift and deformation in the TP caused by the India–Asia collision is still incomplete. Due to the environmental complexity, the gravity method is indispensable to investigate the evolution of the TP. This study concentrates on the Moho depth and crustal density structure in the TP from gravity inversion of Bouguer anomalies. The results show Moho deeper than 60 km in the regions of the TP, suggesting thickening of the crust. Two sinking Moho belts in the southern and northern plateau regions and the linearly increasing Moho depth from the Indian Plate (IP) to the Indus-Yalu suture can be used to infer the crustal fold that has resulted from the India–Asia collision. On the other hand, the density structures show the lower density is commonly found in the crust and the underlying lithospheric mantle beneath the TP, contrasting with the high density in the surrounding blocks. Notably, the high density of the IP is observed underneath the Himalayas, suggesting that the Indian lithosphere extends northward, at least reaching the Indus-Yalu suture. Corresponding to the sinking Moho belts, the crustal densities in these regions present relatively low, which may be evidence for the absence of the eclogites in the lower crust beneath the Himalayas and the Lhasa terrane. In contrast, the relatively high densities underneath the Bangong-Nujiang suture are potentially contributions to the interpretation of the eclogitised lower crust.
Edited by
Alik Ismail-Zadeh, Karlsruhe Institute of Technology, Germany,Fabio Castelli, Università degli Studi, Florence,Dylan Jones, University of Toronto,Sabrina Sanchez, Max Planck Institute for Solar System Research, Germany
Abstract: This chapter provides a broad introduction to Bayesian data assimilation that will be useful to practitioners in interpreting algorithms and results, and for theoretical studies developing novel schemes with an understanding of the rich history of geophysical data assimilation and its current directions. The simple case of data assimilation in a ‘perfect’ model is primarily discussed for pedagogical purposes. Some mathematical results are derived at a high level in order to illustrate key ideas about different estimators. However, the focus of this chapter is on the intuition behind these methods, where more formal and detailed treatments of the data assimilation problem can be found in the various references. In surveying a variety of widely used data assimilation schemes, the key message of this chapter is how the Bayesian analysis provides a consistent framework for the estimation problem and how this allows one to formulate its solution in a variety of ways to exploit the operational challenges in the geosciences.
Edited by
Alik Ismail-Zadeh, Karlsruhe Institute of Technology, Germany,Fabio Castelli, Università degli Studi, Florence,Dylan Jones, University of Toronto,Sabrina Sanchez, Max Planck Institute for Solar System Research, Germany
Abstract: During the past two decades, there have been significant efforts to better quantify emissions of environmentally important trace gases along with their trends. In particular, there has been a clear need for robust estimates of emissions on policy-relevant scales of trace gases that impact air quality and climate. This need has driven the expansion of the observing network to better monitor the changing composition of the atmosphere. This chapter will discuss the use of various data assimilation and inverse modelling approaches to quantify these emissions, with a focus on the use of satellite observations. It will discuss the inverse problem of retrieving the atmospheric trace gas information from the satellite measurements, and the subsequent use of these satellite data for quantifying sources and sinks of the trace gases.
Edited by
Alik Ismail-Zadeh, Karlsruhe Institute of Technology, Germany,Fabio Castelli, Università degli Studi, Florence,Dylan Jones, University of Toronto,Sabrina Sanchez, Max Planck Institute for Solar System Research, Germany
Abstract: The primary observables of the Global Positioning System (GPS) ground tracking sites for geodynamics are the Earth’s surface motions, and their geophysical interpretation is based on the numerical models of various tectonic processes. The key issues for geophysical interpretation of the GPS observations are adequate mechanical models of brittle and ductile rock behaviour used to predict surface motions related to various tectonic processes, and the corresponding inversion techniques which allow separation of the processes, and evaluation of their parameters. For large-scale heterogeneous processes, the inversion of the GPS observations requires regularisation because it implies evaluation of some complicated distributed underground motions from their discrete manifestation at the surface. One of the fastest growing applications of the satellite geodetic observations is investigation of the seismotectonic deformation associated with great earthquakes worldwide at all stages of the seismic cycle – inter-seismic, co-seismic, post-seismic. The inversion techniques based on dislocation models in elastic or viscoelastic medium is one of the approaches that may be widely used for GPS-based studies of various seismotectonic deformations.