Numerical Analysis and Computation Commons^™

Mathematical Modelling Of Hybrid Photonic Structures For Holographic Sensors, Jack Lyons

Doctoral

This thesis outlines a mathematical framework for modelling the formation of holographic gratings in hybrid photopolymer based nanocomposites with the aim of optimising their holographic recording properties for optical sensing applications. Thus, the second aim of the work is to model the change in optical properties of the grating in response to exposure to a target analyte. This work has been a collaborative research project between the School of Mathematics & Statistics at Technological University Dublin and the Centre for Industrial and Engineering Optics that have done extensive experimental work with holographic gratings recorded in photopolymer materials.

In recent years, …

Mathematical Modelling Of Disease Outbreak, Favour Christian, Matthew Molloy

SURE Journal: Science Undergraduate Research Experience Journal

Establishing a model framework for more research necessitates a thorough understanding of the causes, distribution, prevalence, and evolution of infectious illnesses. The main mathematical concept used in this modelling simulation is ordinary differential equations (ODEs). The purpose of this study was to investigate the significance of the many criteria linked to a zombie virus spread. The zombie framework provides an accessible and relatively simple representation of the nature of infectious disease spread, allowing for tractable assumptions and the development of more complex situations.

The models are designed around a zombie outbreak in which the zombie virus is spread through a …

Modelling Spherical Aberration Detection In An Analog Holographic Wavefront Sensor, Emma Branigan, Suzanne Martin, Matthew Sheehan, Kevin Murphy

Conference Papers

The analog holographic wavefront sensor (AHWFS) is a simple and robust solution to wavefront sensing in turbulent environments. Here, the ability of a photopolymer based AHWFS to detect refractively generated spherical aberration is modelled and verified.

Five-Wave Resonances In Deep Water Gravity Waves: Integrability, Numerical Simulations And Experiments, Dan Lucas, Marc Perlin, Dian-Yong Liu, Shane Walsh, Rossen Ivanov, Miguel D. Bustamante

Articles

In this work we consider the problem of finding the simplest arrangement of resonant deep water gravity waves in one-dimensional propagation, from three perspectives: Theoretical, numerical and experimental. Theoretically this requires using a normal-form Hamiltonian that focuses on 5-wave resonances. The simplest arrangement is based on a triad of wave vectors K₁ + K₂ = K₃ (satisfying specific ratios) along with their negatives, corresponding to a scenario of encountering wave packets, amenable to experiments and numerical simulations. The normal-form equations for these encountering waves in resonance are shown to be non-integrable, but they admit an integrable reduction …

Runge–Kutta–Gegenbauer Explicit Methods For Advection-Diffusion Problems, Stephen O'Sullivan

Articles

In this paper, Runge-Kutta-Gegenbauer (RKG) stability polynomials of arbitrarily high order of accuracy are introduced in closed form. The stability domain of RKG polynomials extends in the the real direction with the square of polynomial degree, and in the imaginary direction as an increasing function of Gegenbauer parameter. Consequently, the polynomials are naturally suited to the construction of high order stabilized Runge-Kutta (SRK) explicit methods for systems of PDEs of mixed hyperbolic-parabolic type.

We present SRK methods composed of L ordered forward Euler stages, with complex-valued stepsizes derived from the roots of RKG stability polynomials of degree $L$. Internal stability …

Surface Waves Over Currents And Uneven Bottom, Alan Compelli, Rossen Ivanov, Calin I. Martin, Michail D. Todorov

Articles

The propagation of surface water waves interacting with a current and an uneven bottom is studied. Such a situation is typical for ocean waves where the winds generate currents in the top layer of the ocean. The role of the bottom topography is taken into account since it also influences the local wave and current patterns. Specific scaling of the variables is selected which leads to approximations of Boussinesq and KdV types. The arising KdV equation with variable coefficients, dependent on the bottom topography, is studied numerically when the initial condition is in the form of the one soliton solution …

Factorized Runge-Kutta-Chebyshev Methods, Stephen O'Sullivan

Conference papers

The second-order extended stability Factorized Runge-Kutta-Chebyshev (FRKC2) class of explicit schemes for the integration of large systems of PDEs with diffusive terms is presented. FRKC2 schemes are straightforward to implement through ordered sequences of forward Euler steps with complex stepsizes, and easily parallelised for large scale problems on distributed architectures.

Preserving 7 digits for accuracy at 16 digit precision, the schemes are theoretically capable of maintaining internal stability at acceleration factors in excess of 6000 with respect to standard explicit Runge-Kutta methods. The stability domains have approximately the same extents as those of RKC schemes, and are a third longer …

Congruent Visual Speech Enhances Entrainment To Continuous Auditory Speech In Noise-Free Conditions, Michael Crosse, John S. Butler, Edmumd Lalor

Articles

Congruent audiovisual speech enhances our ability to comprehend a speaker, even in noise-free conditions. When incongruent auditory and visual information is presented concurrently, it can hinder a listener’s perception and even cause him or her to perceive information that was not presented in either modality. Efforts to investigate the neural basis of these effects have often focused on the special case of discrete audiovisual syllables that are spatially and temporally congruent, with less work done on the case of natural, continuous speech. Recent electrophysiological studies have demonstrated that cortical response measures to continuous auditory speech can be easily obtained using …

A Class Of High-Order Runge-Kutta-Chebyshev Stability Polynomials, Stephen O'Sullivan

Articles

The analytic form of a new class of factorized Runge-Kutta-Chebyshev (FRKC) stability polynomials of arbitrary order N is presented. Roots of FRKC stability polynomials of degree L = MN are used to construct explicit schemes comprising L forward Euler stages with internal stability ensured through a sequencing algorithm which limits the internal amplification factors to ~ L². The associated stability domain scales as M² along the real axis. Marginally stable real-valued points on the interior of the stability domain are removed via a prescribed damping procedure. By construction, FRKC schemes meet all linear order conditions; for nonlinear …

Integrability, Recursion Operators And Soliton Interactions, Boyka Aneva, Georgi Grahovski, Rossen Ivanov, Dimitar Mladenov

Book chapter/book

This volume contains selected papers based on the talks,presentedat the Conference Integrability, Recursion Operators and Soliton Interactions, held in Sofia, Bulgaria (29-31 August 2012) at the Institute for Nuclear Research and Nuclear Energy of the Bulgarian Academy of Sciences. Included are also invited papers presenting new research developments in the thematic area. The Conference was dedicated to the 65-th birthday of our esteemed colleague and friend Vladimir Gerdjikov. The event brought together more than 30 scientists, from 6 European countries to celebrate Vladimir's scientific achievements. All participants enjoyed a variety of excellent talks in a friendly and stimulating atmosphere. …

A Linearised Singularly Perturbed Convection-Diffusion Problem With An Interior Layer, Eugene O'Riordan, Jason Quinn

Articles

A linear time dependent singularly perturbed convection-diffusion problem is examined. The convective coefficient contains an interior layer (with a hyperbolic tangent profile), which in turn induces an interior layer in the solution. A numerical method consisting of a monotone finite difference operator and a piecewise-uniform Shishkin mesh is constructed and analysed. Neglecting logarithmic factors, first order parameter uniform convergence is established.

A Numerical Method For A Nonlinear Singularly Perturbed Interior Layer Problem Using An Approximate Layer Location, Jason Quinn

Articles

A class of nonlinear singularly perturbed interior layer problems is examined in this paper. Solutions exhibit an interior layer at an a priori unknown location. A numerical method is presented that uses a piecewise uniform mesh refined around approximations to the first two terms of the asymptotic expansion of the interior layer location. The first term in the expansion is used exactly in the construction of the approximation which restricts the range of problem data considered. The method is shown to converge point-wise to the true solution with a first order convergence rate (overlooking a logarithmic factor) for sufficiently small …

Singular Solutions Of Coss-Coupled Epdiff Equations: Waltzing Peakons And Compacton Pairs, Colin Cotter, Darryl Holm, Rossen Ivanov, James Percival

Conference papers

We introduce EPDiff equations as Euler-Poincare´ equations related to Lagrangian provided by a metric, invariant under the Lie Group Diff(Rⁿ). Then we proceed with a particular form of EPDiff equations, a cross coupled two-component system of Camassa-Holm type. The system has a new type of peakon solutions, 'waltzing' peakons and compacton pairs.

A Generalized Nonlinear Model For The Evolution Of Low Frequency Freak Waves, Jonathan Blackledge

Articles

This paper presents a generalized model for simulating wavefields associated with the sea surface. This includes the case when `freak waves' may occur through an effect compounded in the nonlinear (cubic) Schrodinger equation. After providing brief introductions to linear sea wave models, `freak waves' and the linear and nonlinear Schrodinger equations, we present a unified model that provides for a piecewise continuous transition from a linear to a nonlinear state. This is based on introducing a fractional time derivative to develop a fractional nonlinear partial differential equation with a stochastic source function. In order to explore the characteristics of this …

Tight Lower Bound For The Sparse Travelling Salesman Problem, Fredrick Mtenzi

Conference papers

The Sparse Travelling Salesman Problem (Sparse TSP) which is a variant of the classical Travelling Salesman Problem (TSP) is the problem of finding the shortest route of the salesman when visiting cities in a region making sure that each city is visited at least once and returning home at the end. In the Sparse TSP, the distance between cities may not obey the triangle inequality; this makes the use of algorithms and formulations designed for the TSP to require modifications in order to produce near-optimal results. A lower bound for optmisation problems gives us the quality guarantee of the near-optimal …

A Covert Encryption Method For Applications In Electronic Data Interchange, Jonathan Blackledge, Dmitry Dubovitskiy

Articles

A principal weakness of all encryption systems is that the output data can be ‘seen’ to be encrypted. In other words, encrypted data provides a ‘flag’ on the potential value of the information that has been encrypted. In this paper, we provide a new approach to ‘hiding’ encrypted data in a digital image.

In conventional (symmetric) encryption, the plaintext is usually represented as a binary stream and encrypted using an XOR type operation with a binary cipher. The algorithm used is ideally designed to: (i) generate a maximum entropy cipher so that there is no bias with regard to any …

Diffusion And Fractional Diffusion Based Models For Multiple Light Scattering And Image Analysis, Jonathan Blackledge

Articles

This paper considers a fractional light diffusion model as an approach to characterizing the case when intermediate scattering processes are present, i.e. the scattering regime is neither strong nor weak. In order to introduce the basis for this approach, we revisit the elements of formal scattering theory and the classical diffusion problem in terms of solutions to the inhomogeneous wave and diffusion equations respectively. We then address the significance of these equations in terms of a random walk model for multiple scattering. This leads to the proposition of a fractional diffusion equation for modelling intermediate strength scattering that is based …

Simulation Of Dynamic Electrochemical Processes, John Cassidy

Articles

This work is designed to introduce electrochemists in a tutorial manner to the basics of modeling of electrochemical systems based primarily on diffusion equations. There is an introduction to analytical and numerical methods with examples taken from typical electrochemical experiments. The Laplace transform is used to derive the Cottrell equation and chronopotentiometry. The response of an electrode to a Gaussian concentration profile is detailed. Laplace’s equation is solved for a simple cell to determine the potential distribution. Discrete methods are employed to calculate the current time behavior following a potential step using the explicit finite difference method. Cyclic voltammetry is …