Coordinate-Independent Computations on Differential Equations

This project investigates the computational representation of differentiable manifolds, with the primary goal of solving partial differential equations using multiple coordinate systems on general n- dimensional spaces. In the process, this abstraction is used to perform accurate integrations of ordinary differential equations using multiple coordinate systems. In the case of linear partial differential equations, however, unexpected difficulties arise even with the simplest equations.

A Modern Differential Geometric Approach to Shape from Shading

How the visual system extracts shape information from a single grey-level image can be approached by examining how the information about shape is contained in the image. This technical report considers the characteristic equations derived by Horn as a dynamical system. Certain image critical points generate dynamical system critical points. The stable and unstable manifolds of these critical points correspond to convex and concave solution surfaces, giving more general existence and uniqueness results. A new kind of highly parallel, robust shape from shading algorithm is suggested on neighborhoods of these critical points. The information at bounding contours in the image is also analyzed.

Automatic Qualitative Analysis of Ordinary Differential Equations Using Piecewise Linear Approximations

This paper explores automating the qualitative analysis of physical systems. It describes a program, called PLR, that takes parameterized ordinary differential equations as input and produces a qualitative description of the solutions for all initial values. PLR approximates intractable nonlinear systems with piecewise linear ones, analyzes the approximations, and draws conclusions about the original systems. It chooses approximations that are accurate enough to reproduce the essential properties of their nonlinear prototypes, yet simple enough to be analyzed completely and efficiently. It derives additional properties, such as boundedness or periodicity, by theoretical methods. I demonstrate PLR on several common nonlinear systems and on published examples from mechanical engineering.

Developing the ability for making evaluative judgements.

It is suggested that a more specific emphasis should be placed in undergraduate education on the explicit development of the ability to make evaluative judgements. This higher level cognitive ability is highlighted as the foundation for much sound and successful personal and professional development throughout education, and in lifelong development. Yet it seldom seems to receive the explicit attention which is given in curricula in higher education to analysis and creativity. The paper includes several examples of activities which have been used and judged effective in developing the ability to make evaluative judgements, often at an early stage in undergraduate courses. Keywords: evaluative judgement; personal development; self-assessment; peerassessment; lifelong learning Outline

Public service employment and the public-private wage differential in British regions

Thomas, Dennis, Henley, Andrew, 'Public service employment and the public-private wage differential in British regions', Regional Studies (2001) 35(3) pp.229-240 RAE2008

Singular perturbation of quantum stochastic differential equations with coupling through an oscillator mode

Gough, John; Van Handel, R., (2007) 'Singular perturbation of quantum stochastic differential equations with coupling through an oscillator mode', Journal of Statistical Physics 127(3) pp.575-607 RAE2008

An assessment of the ability of the stay-green phenotype in Lolium species to provide an improved protein supply for ruminants

Kingston-Smith, A. H., Bollard, A. L., Humphreys, M. O,, Theodorou, M. K. (2002). An assessment of the ability of the stay-green phenotype in Lolium species to provide an improved protein supply for ruminants. Annals of Botany, 89(6), 731-740. Sponsorship: BBSRC/MAFF/Milk Development Council/Meat and Livestock Commission/Industry. RAE2008

Uniformly convergent finite element and finite difference methods for singularly perturbed ordinary differential equations

This thesis is concerned with uniformly convergent finite element and finite difference methods for numerically solving singularly perturbed two-point boundary value problems. We examine the following four problems: (i) high order problem of reaction-diffusion type; (ii) high order problem of convection-diffusion type; (iii) second order interior turning point problem; (iv) semilinear reaction-diffusion problem. Firstly, we consider high order problems of reaction-diffusion type and convection-diffusion type. Under suitable hypotheses, the coercivity of the associated bilinear forms is proved and representation results for the solutions of such problems are given. It is shown that, on an equidistant mesh, polynomial schemes cannot achieve a high order of convergence which is uniform in the perturbation parameter. Piecewise polynomial Galerkin finite element methods are then constructed on a Shishkin mesh. High order convergence results, which are uniform in the perturbation parameter, are obtained in various norms. Secondly, we investigate linear second order problems with interior turning points. Piecewise linear Galerkin finite element methods are generated on various piecewise equidistant meshes designed for such problems. These methods are shown to be convergent, uniformly in the singular perturbation parameter, in a weighted energy norm and the usual L2 norm. Finally, we deal with a semilinear reaction-diffusion problem. Asymptotic properties of solutions to this problem are discussed and analysed. Two simple finite difference schemes on Shishkin meshes are applied to the problem. They are proved to be uniformly convergent of second order and fourth order respectively. Existence and uniqueness of a solution to both schemes are investigated. Numerical results for the above methods are presented.

Uniformly convergent finite element methods for singularly perturbed parabolic partial differential equations

This thesis is concerned with uniformly convergent finite element methods for numerically solving singularly perturbed parabolic partial differential equations in one space variable. First, we use Petrov-Galerkin finite element methods to generate three schemes for such problems, each of these schemes uses exponentially fitted elements in space. Two of them are lumped and the other is non-lumped. On meshes which are either arbitrary or slightly restricted, we derive global energy norm and L2 norm error bounds, uniformly in the diffusion parameter. Under some reasonable global assumptions together with realistic local assumptions on the solution and its derivatives, we prove that these exponentially fitted schemes are locally uniformly convergent, with order one, in a discrete L∞norm both outside and inside the boundary layer. We next analyse a streamline diffusion scheme on a Shishkin mesh for a model singularly perturbed parabolic partial differential equation. The method with piecewise linear space-time elements is shown, under reasonable assumptions on the solution, to be convergent, independently of the diffusion parameter, with a pointwise accuracy of almost order 5/4 outside layers and almost order 3/4 inside the boundary layer. Numerical results for the above schemes are presented. Finally, we examine a cell vertex finite volume method which is applied to a model time-dependent convection-diffusion problem. Local errors away from all layers are obtained in the l2 seminorm by using techniques from finite element analysis.

Variations in total and differential milk somatic cell counts and plasmin levels and their role in proteolysis and quality of milk and cheese

Increased plasmin and plasminogen levels and elevated somatic cell counts (SCC) and polymorphonuclear leucocyte levels (PMN) were evident in late lactation milk. Compositional changes in these milks were associated with increased SCC. The quality of late lactation milks was related to nutritional status of herds, with milks from herds on a high plane of nutrition having composition and clotting properties similar to, or superior to, early-mid lactation milks. Nutritionally-deficient cows had elevated numbers of polymorphonuclear leucocytes (PMNs) in their milk, elevated plasmin levels and increased overall proteolytic activity. The dominant effect of plasmin on proteolysis in milks of low SCC was established. When present in elevated numbers, somatic cells and PMNs in particular had a more significant influence on the proteolysis of both raw and pasteurised milks than plasmin. PMN protease action on the caseins showed proteolysis products of two specific enzymes, cathepsin B and elastase, which were also shown in high SCC milk. Crude extracts of somatic cells had a high specificity on αs1-casein. Cheeses made from late lactation milks had increased breakdown of αs1-casein, suggestive of the action of somatic cell proteinases, which may be linked to textural defects in cheese. Late lactation cheeses also showed decreased production of small peptides and amino acids, the reason for which is unknown. Plasmin, which is elevated in activity in late lactation milk, accelerated the ripening of Gouda-type cheese, but was not associated with defects of texture or flavour. The retention of somatic cell enzymes in cheese curd was confirmed, and a potential role in production of bitter peptides identified. Cheeses made from milks containing high levels of PMNs had accelerated αs1-casein breakdown relative to cheeses made from low PMN milk of the same total SCC, consistent with the demonstrated action of PMN proteinases. The two types of cheese were determined significantly different by blind triangle testing.

Differential and numerical models of hysteretic systems with stochastic and deterministic inputs

Many deterministic models with hysteresis have been developed in the areas of economics, finance, terrestrial hydrology and biology. These models lack any stochastic element which can often have a strong effect in these areas. In this work stochastically driven closed loop systems with hysteresis type memory are studied. This type of system is presented as a possible stochastic counterpart to deterministic models in the areas of economics, finance, terrestrial hydrology and biology. Some price dynamics models are presented as a motivation for the development of this type of model. Numerical schemes for solving this class of stochastic differential equation are developed in order to examine the prototype models presented. As a means of further testing the developed numerical schemes, numerical examination is made of the behaviour near equilibrium of coupled ordinary differential equations where the time derivative of the Preisach operator is included in one of the equations. A model of two phenotype bacteria is also presented. This model is examined to explore memory effects and related hysteresis effects in the area of biology. The memory effects found in this model are similar to that found in the non-ideal relay. This non-ideal relay type behaviour is used to model a colony of bacteria with multiple switching thresholds. This model contains a Preisach type memory with a variable Preisach weight function. Shown numerically for this multi-threshold model is a pattern formation for the distribution of the phenotypes among the available thresholds.

Profiling side-channel attacks on cryptographic algorithms

Traditionally, attacks on cryptographic algorithms looked for mathematical weaknesses in the underlying structure of a cipher. Side-channel attacks, however, look to extract secret key information based on the leakage from the device on which the cipher is implemented, be it smart-card, microprocessor, dedicated hardware or personal computer. Attacks based on the power consumption, electromagnetic emanations and execution time have all been practically demonstrated on a range of devices to reveal partial secret-key information from which the full key can be reconstructed. The focus of this thesis is power analysis, more specifically a class of attacks known as profiling attacks. These attacks assume a potential attacker has access to, or can control, an identical device to that which is under attack, which allows him to profile the power consumption of operations or data flow during encryption. This assumes a stronger adversary than traditional non-profiling attacks such as differential or correlation power analysis, however the ability to model a device allows templates to be used post-profiling to extract key information from many different target devices using the power consumption of very few encryptions. This allows an adversary to overcome protocols intended to prevent secret key recovery by restricting the number of available traces. In this thesis a detailed investigation of template attacks is conducted, along with how the selection of various attack parameters practically affect the efficiency of the secret key recovery, as well as examining the underlying assumption of profiling attacks in that the power consumption of one device can be used to extract secret keys from another. Trace only attacks, where the corresponding plaintext or ciphertext data is unavailable, are then investigated against both symmetric and asymmetric algorithms with the goal of key recovery from a single trace. This allows an adversary to bypass many of the currently proposed countermeasures, particularly in the asymmetric domain. An investigation into machine-learning methods for side-channel analysis as an alternative to template or stochastic methods is also conducted, with support vector machines, logistic regression and neural networks investigated from a side-channel viewpoint. Both binary and multi-class classification attack scenarios are examined in order to explore the relative strengths of each algorithm. Finally these machine-learning based alternatives are empirically compared with template attacks, with their respective merits examined with regards to attack efficiency.

Evaluation of rhizobacterial colonisation and the ability to induce Globodera pallida hatch

Three bacterial isolates, SB13 (Acinetobacter sp.), SB14 (Arthrobacter sp.) and SB15 (Bacillus sp.), were previously isolated from the rhizosphere of sugar beet (Beta vulgaris ssp. vulgaris) plants and shown to increase hatch of potato cyst nematodes in vitro. In this study, the three isolates were assayed for rhizosphere competence. Each isolate was applied to seeds at each of four concentrations (105-108 CFU ml−1) and the inoculated seeds were planted in plastic microcosms containing coarse sand. All three isolates were shown to colonise the rhizosphere, although to differing degrees, with the higher inoculation densities providing significantly better colonisation. The isolates increased sugar beet root and shoot dry weight. Isolates SB14 and SB15 were analysed for their ability to induce in vivo hatch of Globodera pallida in non-sterile soil planted with sugar beet. After 4 and 6 weeks, both isolates had induced significantly greater percentage hatch compared to controls.

Optimal assessment of the ability of children with recurrent respiratory tract infections to produce anti-polysaccharide antibodies.

Specific anti-polysaccharide antibody deficiency (SPAD) is an immune disorder. Diagnostic criteria have not yet been defined clearly. One hundred and seventy-six children evaluated for recurrent respiratory tract infections were analysed retrospectively. For each subject, specific anti-pneumococcal antibodies had been measured with two enzyme-linked immunosorbent assays (ELISAs), one overall assay (OA) using the 23-valent pneumococcal polysaccharide vaccine (23-PPSV) as detecting antigen and the other purified pneumococcal polysaccharide serotypes (serotype-specific assay, SSA) (serotypes 14, 19F and 23F). Antibody levels were measured before (n = 176) and after (n = 93) immunization with the 23-PPSV. Before immunization, low titres were found for 138 of 176 patients (78%) with OA, compared to 20 of 176 patients (11%) with the SSA. We found a significant correlation between OA and SSA results. After immunization, 88% (71 of 81) of the patients considered as responders in the OA test were also responders in the SSA; 93% (71 of 76) of the patients classified as responders according to the SSA were also responders in the OA. SPAD was diagnosed in 8% (seven of 93) of patients on the basis of the absence of response in both tests. Thus, we propose to use OA as a screening test for SPAD before 23-PPSV immunization. After immunization, SSA should be used only in case of a low response in OA. Only the absence of or a very low antibody response detected by both tests should be used as a diagnostic criterion for SPAD.

Differential inducibility of HLA class I and II antigens by r-IFN gamma in type III bare lymphocyte syndrome.

Case Reports