No description available.

Reaction-diffusion equations have been used to study various phenomena across different fields. These equations can be posed on the whole real line, or on a subinterval, depending on the situation being studied. For finite intervals, we also impose diverse boundary conditions on the system. In the present thesis, we solely focus on the bistable reaction-diffusion equation while working on a bounded interval of the form $[0,L]$ ($L>0$). Furthermore, we consider both mixed and no-flux boundary conditions, where we extend the former to Dirichlet boundary conditions once our analysis of that system is complete. We first use phase-plane analysis to set up our initial investigation of both systems. This gives us an integral describing the transit time of orbits within the phase-plane. This allows us to determine the bifurcation diagram of both systems. We then transform the integral to ease numerical calculations. Finally, we determine the stability of the steady states of each system.

In the present study, effect of satisfying stress boundary conditions, in addition to displacement boundary conditions, in the axisymmetric vibration analysis of circular and annular plates is investigated. A new axisymmetric finite element, which is based on a combination of the conventional displacement-type variational principle and the Reissner¡¦s principle, is proposed. With this formulation, stresses, like displacements, are primary variables, and both displacement and stress boundary conditions can be easily and exactly imposed. Axisymmetric vibration frequencies of some typical circular and annular plates are then obtained with the present approach and are compared with those by the displacement-type axisymmetric finite element. It is found that the conventional finite element, though not satisfying stress boundary conditions, can still obtain sufficiently accurate vibration frequencies of circular and annular plates.

none

This work presents a method to quantitatively calculate the error induced through application of approximate boundary conditions in quantum charge transport simulations based on the Wigner transport equation (WTE). Except for the special case of homogeneous material, there exists no methodology for the calculation of exact boundary conditions. Consequently, boundary conditions are customarily approximated by equilibrium or near-equilibrium distributions known to be correct in the classical limit. This practice can, however, exert deleterious impact on the accuracy of numerical calculations and can even lead to unphysical results. The Yoder group has recently developed a series expansion for exact boundary conditions which, when truncated, can be used to calculate boundary conditions of successively greater accuracy through consideration of successively higher order terms, the computational penalty for which is however not to be underestimated. This thesis focuses on the calculation and analysis of the second order term of the series expansion. A method is demonstrated to calculate the term for any general device structure in one spatial dimension. In addition, numerical analysis is undertaken to directly compare the first and second order terms. Finally a method to incorporate the first order term into simulation is formulated.

The aim of this work was to develop an envelope function method to calculate the electronic states and optical properties of complex quantum well heterostructures, and to demonstrate its effectiveness by application to some device structures of topical interest. In particular, structures have been considered which might form the basis of intensity modulators and polarization insensitive amplifier devices for light at a wavelength of 1.55 µm. The modulator structures considered all have the general form of two coupled quantum wells of different widths as the active region. The application of an electric field in the growth direction is intended to result in a shift in the energy and spatial localisation of the confined states and produce an increase in the absorption coefficient at longer wavelengths than the zero field absorption edge. The effectiveness of certain structures is examined in terms of field induced absorption increase at 1.55 µm. A system which shows a significant increase in absorption coefficient at this wavelength on application of a practical electric field has been identified as a possible candidate for an intensity modulator. In the case of the amplifier, the active region of the most promising structure considered consists of a stepped well which comprises two layers, one with tensile and one with compressive strain. It is known that the presence of the two oppositely strained layers can result in the TE and TM gain peaks appearing at similar photon energies. Our calculations show that a suitable choice of strain and layer widths can result in a small or zero difference between the TE and TM gains at 1.55 µm, which can be important for the polarization insensitive operation of devices in optical communications applications. In order to predict the optical properties of quantum well devices it is necessary to calculate the electron and hole states for a range of in-plane wavevectors. The calculations developed and carried out in this work are based on a multi-layer (eight band) k.p model including strain effects. The interfacial boundary conditions which result from approximations to Burt's exact envelope function theory are included in the model. The effect of an electric field is modelled by including a potential energy term in each layer Hamiltonian which is equal to the average energy shift across the layer in question due to the presence of the field. The model has been developed with flexibility in mind and has applications beyond the specific devices considered in this thesis.

Doctor of Philosophy / We consider the problem of minimising the eigenvalues of the Laplacian with Robin boundary conditions $\frac{\partial u}{\partial \nu} + \alpha u = 0$ and generalised Wentzell boundary conditions $\Delta u + \beta \frac{\partial u}{\partial \nu} + \gamma u = 0$ with respect to the domain $\Omega \subset \mathbb R^N$ on which the problem is defined. For the Robin problem, when $\alpha > 0$ we extend the Faber-Krahn inequality of Daners [Math. Ann. 335 (2006), 767--785], which states that the ball minimises the first eigenvalue, to prove that the minimiser is unique amongst domains of class $C^2$. The method of proof uses a functional of the level sets to estimate the first eigenvalue from below, together with a rearrangement of the ball's eigenfunction onto the domain $\Omega$ and the usual isoperimetric inequality. We then prove that the second eigenvalue attains its minimum only on the disjoint union of two equal balls, and set the proof up so it works for the Robin $p$-Laplacian. For the higher eigenvalues, we show that it is in general impossible for a minimiser to exist independently of $\alpha > 0$. When $\alpha < 0$, we prove that every eigenvalue behaves like $-\alpha^2$ as $\alpha \to -\infty$, provided only that $\Omega$ is bounded with $C^1$ boundary. This generalises a result of Lou and Zhu [Pacific J. Math. 214 (2004), 323--334] for the first eigenvalue. For the Wentzell problem, we (re-)prove general operator properties, including for the less-studied case $\beta < 0$, where the problem is ill-posed in some sense. In particular, we give a new proof of the compactness of the resolvent and the structure of the spectrum, at least if $\partial \Omega$ is smooth. We prove Faber-Krahn-type inequalities in the general case $\beta, \gamma \neq 0$, based on the Robin counterpart, and for the ``best'' case $\beta, \gamma > 0$ establish a type of equivalence property between the Wentzell and Robin minimisers for all eigenvalues. This yields a minimiser of the second Wentzell eigenvalue. We also prove a Cheeger-type inequality for the first eigenvalue in this case.

Doctor of Philosophy / We consider the problem of minimising the eigenvalues of the Laplacian with Robin boundary conditions $\frac{\partial u}{\partial \nu} + \alpha u = 0$ and generalised Wentzell boundary conditions $\Delta u + \beta \frac{\partial u}{\partial \nu} + \gamma u = 0$ with respect to the domain $\Omega \subset \mathbb R^N$ on which the problem is defined. For the Robin problem, when $\alpha > 0$ we extend the Faber-Krahn inequality of Daners [Math. Ann. 335 (2006), 767--785], which states that the ball minimises the first eigenvalue, to prove that the minimiser is unique amongst domains of class $C^2$. The method of proof uses a functional of the level sets to estimate the first eigenvalue from below, together with a rearrangement of the ball's eigenfunction onto the domain $\Omega$ and the usual isoperimetric inequality. We then prove that the second eigenvalue attains its minimum only on the disjoint union of two equal balls, and set the proof up so it works for the Robin $p$-Laplacian. For the higher eigenvalues, we show that it is in general impossible for a minimiser to exist independently of $\alpha > 0$. When $\alpha < 0$, we prove that every eigenvalue behaves like $-\alpha^2$ as $\alpha \to -\infty$, provided only that $\Omega$ is bounded with $C^1$ boundary. This generalises a result of Lou and Zhu [Pacific J. Math. 214 (2004), 323--334] for the first eigenvalue. For the Wentzell problem, we (re-)prove general operator properties, including for the less-studied case $\beta < 0$, where the problem is ill-posed in some sense. In particular, we give a new proof of the compactness of the resolvent and the structure of the spectrum, at least if $\partial \Omega$ is smooth. We prove Faber-Krahn-type inequalities in the general case $\beta, \gamma \neq 0$, based on the Robin counterpart, and for the ``best'' case $\beta, \gamma > 0$ establish a type of equivalence property between the Wentzell and Robin minimisers for all eigenvalues. This yields a minimiser of the second Wentzell eigenvalue. We also prove a Cheeger-type inequality for the first eigenvalue in this case.

Wave radiation is an important process in many geophysical flows. In particular, it is by wave radiation that flows may adjust to a state for which the dynamics is slow. Such a state is described as “balanced”, meaning there is an approximate balance between the Coriolis force and horizontal pressure gradients, and between buoyancy and vertical pressure gradients. In this thesis, wave radiation processes relevant to these enormously complex flows are studied through the use of some highly simplified models, and a parallel aim is to develop accurate numerical techniques for doing so. This thesis is divided into three main parts. 1. We consider accurate numerical boundary conditions for various equations which support wave radiation to infinity. Particular attention is given to discretely non-reflecting boundary conditions, which are derived directly from a discretised scheme. Such a boundary condition is studied in the case of the 1-d Klein-Gordon equation. The limitations concerning the practical implementation of this scheme are explored and some possible improvements are suggested. A stability analysis is developed which yields a simple stability criterion that is useful when tuning the boundary condition. The practical use of higher-order boundary conditions for the 2-d shallow water equations is also explored; the accuracy of such a method is assessed when combined with a particular interior scheme, and an analysis based on matrix pseudospectra reveals something of the stability of such a method. 2. Large-scale atmospheric and oceanic flows are examples of systems with a wide timescale separation, determined by a small parameter. In addition they both undergo constant random forcing. The five component Lorenz-Krishnamurthy system is a system with a timescale separation controlled by a small parameter, and we employ it as a model of the forced ocean by further adding a random forcing of the slow variables, and introduce wave radiation to infinity by the addition of a dispersive PDE. The dynamics are reduced by deriving balance relations, and numerical experiments are used to assess the effects of energy radiation by fast waves. 3. We study quasimodes, which demonstrate the existence of associated Landau poles of a system. In this thesis, we consider a simple model of wave radiation that exhibits quasimodes, that allows us to derive some explicit analytical results, as opposed to physically realistic geophysical fluid systems for which such results are often unavailable, necessitating recourse to numerical techniques. The growth rates obtained for this system, which is an extension of one considered by Lamb, are confirmed using numerical experiments.

This thesis studies the quantum reflection factor of the sinh-Gordon model under boundary conditions consistent with integrability. First, we review the affine Toda field theory in Chapter One. In particular, the classical and quantum integrability of the theory are reviewed on the whole line and on the half-line as well, that is, in the presence of a boundary. We next consider the sinh-Gordon model which is restricted to a half-line by boundary conditions maintaining integrability in Chapter Two. A perturbative calculation of the reflection factor is given to one loop order in the bulk coupling and to first order in the difference of the two parameters introduced at the boundary. The result provides a further verification of Ghoshal's formula. The calculation is consistent with a conjecture for the general dependence of the reflection factor on the boundary parameters and the bulk coupling. In Chapter Three, quantum corrections to the classical reflection factor of the sinh-Gordon model are studied up to second order in the difference of boundary data and to one loop order in the bulk coupling. Chapter Four deals with the quantum reflection factor for the sinh-Gordon model with general boundary conditions. The model is studied under boundary conditions which are compatible with integrability and in the framework of the conventional perturbation theory generalised to the affine Toda field theory. It is found that the general form of a subset of the related quantum corrections are hypergeometric functions. Finally, we sum up this thesis in Chapter Five along with some conclusions and suggestions for further future studies.