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.