Asymptotic approximation of eigenvalues of vector equations

A vectorial extension of the Keller-Rubinow method of computing asymptotic approximations of eigenvalues in bounded domains is presented. The method is applied to the problem of a multimode step-profile cylindrical optical fibre, including the effects of polarisation. A comparison of the asymptotic results with the exact eigenvalues is made when these are available, and the agreement is shown to be good

A hierarchy of models for superconducting thin films

A hierarchy of models for type-II superconducting thin films is presented. Through appropriate asymptotic limits this hierarchy passes from the mesoscopic Ginzburg--Landau model to the London model with isolated vortices as $\delta$-function singularities to vortex-density models and finally to macroscopic critical-state models. At each stage it is found that a key nondimensional parameter is $\Lambda = \lambda^2/d L$, where $\lambda$ is the penetration depth of the magnetic field, a material parameter, and d and L are a typical thickness and lateral dimension of the film,respectively. The models simplify greatly if this parameter is large or small

The use of a formal sensitivity analysis on epidemic models with immune protection from maternally acquired antibodies

This paper considers the outcome of a formal sensitivity analysis on a series of epidemic model structures developed to study the population level effects of maternal antibodies. The analysis is used to compare the potential influence of maternally acquired immunity on various age and time domain observations of infection and serology, with and without seasonality. The results of the analysis indicate that time series observations are largely insensitive to variations in the average duration of this protection, and that age related empirical data are likely to be most appropriate for estimating these characteristics

Axisymmetric buckling of a spherical shell embedded in an elastic medium under uniaxial stress at infinity

The problem of a thin spherical linearly-elastic shell, perfectly bonded to an infinite linearly-elastic medium is considered. A constant axisymmetric stress field is applied at infinity in the matrix, and the displacement and stress fields in the shell and matrix are evaluated by means of harmonic potential functions. In order to examine the stability of this solution, the buckling problem of a shell which experiences this deformation is considered. Using Koiter's nonlinear shallow shell theory, restricting buckling patterns to those which are axisymmetric, and using the Rayleigh–Ritz method by expanding the buckling patterns in an infinite series of Legendre functions, an eigenvalue problem for the coefficients in the infinite series is determined. This system is truncated and solved numerically in order to analyse the behaviour of the shell as it undergoes buckling, and to identify the critical buckling stress in two cases — namely where the shell is hollow and the stress at infinity is either uniaxial or radial

Development and Verification of a Flight Stack for a High-Altitude Glider in Ada/SPARK 2014

SPARK 2014 is a modern programming language and a new state-of-the-art tool set for development and verification of high-integrity software. In this paper, we explore the capabilities and limitations of its latest version in the context of building a flight stack for a high-altitude unmanned glider. Towards that, we deliberately applied static analysis early and continuously during implementation, to give verification the possibility to steer the software design. In this process we have identified several limitations and pitfalls of software design and verification in SPARK, for which we give workarounds and protective actions to avoid them. Finally, we give design recommendations that have proven effective for verification, and summarize our experiences with this new language

Cardiac Electromechanics: The effect of contraction model on the mathematical problem and accuracy of the numerical scheme

Models of cardiac electromechanics usually contain a contraction model determining the active tension induced at the cellular level, and the equations of nonlinear elasticity to determine tissue deformation in response to this active tension. All contraction models are dependent on cardiac electro-physiology, but can also be dependent on\ud the stretch and stretch-rate in the fibre direction. This fundamentally affects the mathematical problem being solved, through classification of the governing PDEs, which affects numerical schemes that can be used to solve the governing equations. We categorise contraction models into three types, and for each consider questions such as classification and the most appropriate choice from two numerical methods (the explicit and implicit schemes). In terms of mathematical classification, we consider the question of strong ellipticity of the total strain energy (important for precluding ‘unnatural’ material behaviour) for stretch-rate-independent contraction models; whereas for stretch-rate-dependent contraction models we introduce a corresponding third-order problem and explain how certain choices of boundary condition could lead to constraints on allowable initial condition. In terms of suitable numerical methods, we show that an explicit approach (where the contraction model is integrated in the timestep prior to the bulk deformation being computed) is: (i) appropriate for stretch-independent contraction models; (ii) only conditionally-stable, with the stability criterion independent of timestep, for contractions models which just depend on stretch (but not stretch-rate), and (iii) inappropriate for stretch-rate-dependent models