Towards a Modeling Environment for Earth Systems Simulations

The developments of models in Earth Sciences, e.g. for earthquake prediction and for the simulation of mantel convection, are fare from being finalized. Therefore there is a need for a modelling environment that allows scientist to implement and test new models in an easy but flexible way. After been verified, the models should be easy to apply within its scope, typically by setting input parameters through a GUI or web services. It should be possible to link certain parameters to external data sources, such as databases and other simulation codes. Moreover, as typically large-scale meshes have to be used to achieve appropriate resolutions, the computational efficiency of the underlying numerical methods is important. Conceptional this leads to a software system with three major layers: the application layer, the mathematical layer, and the numerical algorithm layer. The latter is implemented as a C/C++ library to solve a basic, computational intensive linear problem, such as a linear partial differential equation. The mathematical layer allows the model developer to define his model and to implement high level solution algorithms (e.g. Newton-Raphson scheme, Crank-Nicholson scheme) or choose these algorithms form an algorithm library. The kernels of the model are generic, typically linear, solvers provided through the numerical algorithm layer. Finally, to provide an easy-to-use application environment, a web interface is (semi-automatically) built to edit the XML input file for the modelling code. In the talk, we will discuss the advantages and disadvantages of this concept in more details. We will also present the modelling environment escript which is a prototype implementation toward such a software system in Python (see www.python.org). Key components of escript are the Data class and the PDE class. Objects of the Data class allow generating, holding, accessing, and manipulating data, in such a way that the actual, in the particular context best, representation is transparent to the user. They are also the key to establish connections with external data sources. PDE class objects are describing (linear) partial differential equation objects to be solved by a numerical library. The current implementation of escript has been linked to the finite element code Finley to solve general linear partial differential equations. We will give a few simple examples which will illustrate the usage escript. Moreover, we show the usage of escript together with Finley for the modelling of interacting fault systems and for the simulation of mantel convection.

Soliton generation from randomly modulated return-to-zero pulses

We consider return-to-zero (RZ) pulses with random phase modulation propagating in a nonlinear channel (modelled by the integrable nonlinear Schrödinger equation, NLSE). We suggest two different models for the phase fluctuations of the optical field: (i) Gaussian short-correlated fluctuations and (ii) generalized telegraph process. Using the rectangular-shaped pulse form we demonstrate that the presence of phase fluctuations of both types strongly influences the number of solitons generated in the channel. It is also shown that increasing the correlation time for the random phase fluctuations affects the coherent content of a pulse in a non-trivial way. The result obtained has potential consequences for all-optical processing and design of optical decision elements.

Lattice approach to the dynamics of phase-coded soliton trains

We address the collective dynamics of a soliton train propagating in a medium described by the nonlinear Schrödinger equation. Our approach uses the reduction of train dynamics to the discrete complex Toda chain (CTC) model for the evolution of parameters for each train constituent: such a simplification allows one to carry out an approximate analysis of the dynamics of positions and phases of individual interacting pulses. Here, we employ the CTC model to the problem which has relevance to the field of fibre optics communications where each binary digit of transmitted information is encoded via the phase difference between the two adjacent solitons. Our goal is to elucidate different scenarios of the train distortions and the subsequent information garbling caused solely by the intersoliton interactions. First, we examine how the structure of a given phase pattern affects the initial stage of the train dynamics and explain the general mechanisms for the appearance of unstable collective soliton modes. Then we further discuss the nonlinear regime concentrating on the dependence of the Lax scattering matrix on the input phase distribution; this allows one to classify typical features of the train evolution and determine the distance where the soliton escapes from its slot. In both cases, we demonstrate deep mathematical analogies with the classical theory of crystal lattice dynamics.

Numerical modelling of dispersion managed soliton transmission

This thesis presents the results of numerical modelling of the propagation of dispersion managed solitons. The theory of optical pulse propagation in single mode optical fibre is introduced specifically looking at the use of optical solitons for fibre communications. The numerical technique used to solve the nonlinear Schrödinger equation is also introduced. The recent developments in the use of dispersion managed solitons are reviewed before the numerical results are presented. The work in this thesis covers two main areas; (i) the use of a saturable absorber to control the propagation of dispersion managed solutions and (ii) the upgrade of the installed standard fibre network to higher data rates through the use of solitons and dispersion management. Saturable absorbe can be used to suppress the build up of noise and dispersive radiation in soliton transmission lines. The use of saturable absorbers in conjunction with dispersion management has been investigated both as a single pulse and for the transmission of a 10Gbit/s data pattern. It is found that this system supports a new regime of stable soliton pulses with significantly increased powers. The upgrade of the installed standard fibre network to higher data rates through the use of fibre amplifiers and dispersion management is of increasing interest. In this thesis the propagation of data at both 10Gbit/s and 40Gbit/s is studied. Propagation over transoceanic distances is shown to be possible for 10Gbit/s transmission and for more than 2000km at 40Gbit/s. The contribution of dispersion managed solitons in the future of optical communications is discussed in the thesis conclusions.

Exploitation of fibre characteristics in temperature sensing and radio-over-fibre applications

The development of sensing devices is one of the instrumentation fields that has grown rapidly in the last decade. Corresponding to the swift advance in the development of microelectronic sensors, optical fibre sensors are widely investigated because of their advantageous properties over the electronics sensors such as their wavelength multiplexing capability and high sensitivity to temperature, pressure, strain, vibration and acoustic emission. Moreover, optical fibre sensors are more attractive than the electronics sensors as they can perform distributed sensing, in terms of covering a reasonably large area using a single piece of fibre. Apart from being a responsive element in the sensing field, optical fibre possesses good assets in generating, distributing, processing and transmitting signals in the future broadband information network. These assets include wide bandwidth, high capacity and low loss that grant mobility and flexibility for wireless access systems. Among these core technologies, the fibre optic signal processing and transmission of optical and radio frequency signals have been the subjects of study in this thesis. Based on the intrinsic properties of single-mode optical fibre, this thesis aims to exploit the fibre characteristics such as thermal sensitivity, birefringence, dispersion and nonlinearity, in the applications of temperature sensing and radio-over-fibre systems. By exploiting the fibre thermal sensitivity, a fully distributed temperature sensing system consisting of an apodised chirped fibre Bragg grating has been implemented. The proposed system has proven to be efficient in characterising grating and providing the information of temperature variation, location and width of the heat source applied in the area under test.To exploit the fibre birefringence, a fibre delay line filter using a single high-birefringence optical fibre structure has been presented. The proposed filter can be reconfigured and programmed by adjusting the input azimuth of launched light, as well as the strength and direction of the applied coupling, to meet the requirements of signal processing for different purposes in microwave photonic and optical filtering applications. To exploit the fibre dispersion and nonlinearity, experimental investigations have been carried out to study their joint effect in high power double-sideband and single-sideband modulated links with the presence of fibre loss. The experimental results have been theoretically verified based on the in-house implementation of the split-step Fourier method applied to the generalised nonlinear Schrödinger equation. Further simulation study on the inter-modulation distortion in two-tone signal transmission has also been presented so as to show the effect of nonlinearity of one channel on the other. In addition to the experimental work, numerical simulations have also been carried out in all the proposed systems, to ensure that all the aspects concerned are comprehensively investigated.

Numerically optimized implementations, mitigation and simulations of polarisation mode dispersion in real-time transmission systems

The aim of this thesis is to present numerical investigations of the polarisation mode dispersion (PMD) effect. Outstanding issues on the side of the numerical implementations of PMD are resolved and the proposed methods are further optimized for computational efficiency and physical accuracy. Methods for the mitigation of the PMD effect are taken into account and simulations of transmission system with added PMD are presented. The basic outline of the work focusing on PMD can be divided as follows. At first the widely-used coarse-step method for simulating the PMD phenomenon as well as a method derived from the Manakov-PMD equation are implemented and investigated separately through the distribution of a state of polarisation on the Poincaré sphere, and the evolution of the dispersion of a signal. Next these two methods are statistically examined and compared to well-known analytical models of the probability distribution function (PDF) and the autocorrelation function (ACF) of the PMD phenomenon. Important optimisations are achieved, for each of the aforementioned implementations in the computational level. In addition the ACF of the coarse-step method is considered separately, based on the result which indicates that the numerically produced ACF, exaggerates the value of the correlation between different frequencies. Moreover the mitigation of the PMD phenomenon is considered, in the form of numerically implementing Low-PMD spun fibres. Finally, all the above are combined in simulations that demonstrate the impact of the PMD on the quality factor (Q=factor) of different transmission systems. For this a numerical solver based on the coupled nonlinear Schrödinger equation is created which is otherwise tested against the most important transmission impairments in the early chapters of this thesis.

Controlling soliton refraction in optical lattices

We show in the framework of the 1D nonlinear Schrödinger equation that the value of the refraction angle of a fundamental soliton beam passing through an optical lattice can be controlled by adjusting either the shape of an individual waveguide or the relative positions of the waveguides. In the case of the shallow refractive index modulation, we develop a general approach for the calculation of the refraction angle change. The shape of a single waveguide crucially affects the refraction direction due to the appearance of a structural form factor in the expression for the density of emitted waves. For a lattice of scatterers, wave-soliton interference inside the lattice leads to the appearance of an additional geometric form factor. As a result, the soliton refraction is more pronounced for the disordered lattices than for the periodic ones.

Impact of third-order fibre dispersion on the evolution of parabolic optical pulses

We develop a perturbation analysis that describes the effect of third-order dispersion on the similariton pulse solution of the nonlinear Schrodinger equation in a fibre gain medium. The theoretical model predicts with sufficient accuracy the pulse structural changes induced, which are observed through direct numerical simulations.

Parabolic optical pulses under the action of the third-order dispersion

Recent developments in nonlinear optics reveal an interesting class of pulses with a parabolic intensity profile in the energy-containing core and a linear frequency chirp that can propagate in a fiber with normal group-velocity dispersion. Parabolic pulses propagate in a stable selfsimilar manner, holding certain relations (scaling) between pulse power, width, and chirp parameter. In the additional presence of linear amplification, they enjoy the remarkable property of representing a common asymptotic state (or attractor) for arbitrary initial conditions. Analytically, self-similar (SS) parabolic pulses can be found as asymptotic, approximate solutions of the nonlinear Schr¨odinger equation (NLSE) with gain in the semi-classical (largeamplitude/small-dispersion) limit. By analogy with the well-known stable dynamics of solitary waves - solitons, these SS parabolic pulses have come to be known as similaritons. In practical fiber systems, inherent third-order dispersion (TOD) in the fiber always introduces a certain degree of asymmetry in the structure of the propagating pulse, eventually leading to pulse break-up. To date, there is no analytic theory of parabolic pulses under the action of TOD. Here, we develop aWKB perturbation analysis that describes the effect of weak TOD on the parabolic pulse solution of the NLSE in a fiber gain medium. The induced perturbation in phase and amplitude can be found to any order. The theoretical model predicts with sufficient accuracy the pulse structural changes induced by TOD, which are observed through direct numerical NLSE simulations.

Self-similar parabolic optical solitary waves

We study solutions of the nonlinear Schrödinger equation (NLSE) with gain, describing optical pulse propagation in an amplifying medium. We construct a semiclassical self-similar solution with a parabolic temporal variation that corresponds to the energy-containing core of the asymptotically propagating pulse in the amplifying medium. We match the self-similar core through Painlevé functions to the solution of the linearized equation that corresponds to the low-amplitude tails of the pulse. The analytic solution accurately reproduces the numerically calculated solution of the NLSE.

A perturbative analysis of dispersion-managed solitons

We study soliton solutions of the path-averaged propagation equation governing the transmission of dispersion-managed (DM) optical pulses in the (practical) limit when residual dispersion and nonlinearity only slightly affect the pulse dynamics over one compensation period. In the case of small dispersion map strengths, the averaged pulse dynamics is governed by a perturbed form of the nonlinear Schrödinger equation; applying a perturbation theory – elsewhere developed – based on inverse scattering theory, we derive an analytic expression for the envelope of the DM soliton. This expression correctly predicts the power enhancement arising from the dispersion management. Theoretical results are verified by direct numerical simulations.

Impact of third-order dispersion on the evolution of parabolic pulses

We present a perturbation analysis that describes the effect of third-order dispersion on the similariton pulse solution of the nonlinear Schrödinger equation in a fibre gain medium. The theoretical model predicts with sufficient accuracy the pulse structural changes induced, which are observed through direct numerical simulations.

Parabolic optical pulses under the action of the third-order dispersion

Recent developments in nonlinear optics reveal an interesting class of pulses with a parabolic intensity profile in the energy-containing core and a linear frequency chirp that can propagate in a fiber with normal group-velocity dispersion. Parabolic pulses propagate in a stable selfsimilar manner, holding certain relations (scaling) between pulse power, width, and chirp parameter. In the additional presence of linear amplification, they enjoy the remarkable property of representing a common asymptotic state (or attractor) for arbitrary initial conditions. Analytically, self-similar (SS) parabolic pulses can be found as asymptotic, approximate solutions of the nonlinear Schr¨odinger equation (NLSE) with gain in the semi-classical (largeamplitude/small-dispersion) limit. By analogy with the well-known stable dynamics of solitary waves - solitons, these SS parabolic pulses have come to be known as similaritons. In practical fiber systems, inherent third-order dispersion (TOD) in the fiber always introduces a certain degree of asymmetry in the structure of the propagating pulse, eventually leading to pulse break-up. To date, there is no analytic theory of parabolic pulses under the action of TOD. Here, we develop aWKB perturbation analysis that describes the effect of weak TOD on the parabolic pulse solution of the NLSE in a fiber gain medium. The induced perturbation in phase and amplitude can be found to any order. The theoretical model predicts with sufficient accuracy the pulse structural changes induced by TOD, which are observed through direct numerical NLSE simulations.

Lattice approach to the dynamics of phase-coded soliton trains

We address the collective dynamics of a soliton train propagating in a medium described by the nonlinear Schrödinger equation. Our approach uses the reduction of train dynamics to the discrete complex Toda chain (CTC) model for the evolution of parameters for each train constituent: such a simplification allows one to carry out an approximate analysis of the dynamics of positions and phases of individual interacting pulses. Here, we employ the CTC model to the problem which has relevance to the field of fibre optics communications where each binary digit of transmitted information is encoded via the phase difference between the two adjacent solitons. Our goal is to elucidate different scenarios of the train distortions and the subsequent information garbling caused solely by the intersoliton interactions. First, we examine how the structure of a given phase pattern affects the initial stage of the train dynamics and explain the general mechanisms for the appearance of unstable collective soliton modes. Then we further discuss the nonlinear regime concentrating on the dependence of the Lax scattering matrix on the input phase distribution; this allows one to classify typical features of the train evolution and determine the distance where the soliton escapes from its slot. In both cases, we demonstrate deep mathematical analogies with the classical theory of crystal lattice dynamics.

A perturbative analysis of dispersion-managed solitons

We study soliton solutions of the path-averaged propagation equation governing the transmission of dispersion-managed (DM) optical pulses in the (practical) limit when residual dispersion and nonlinearity only slightly affect the pulse dynamics over one compensation period. In the case of small dispersion map strengths, the averaged pulse dynamics is governed by a perturbed form of the nonlinear Schrödinger equation; applying a perturbation theory – elsewhere developed – based on inverse scattering theory, we derive an analytic expression for the envelope of the DM soliton. This expression correctly predicts the power enhancement arising from the dispersion management. Theoretical results are verified by direct numerical simulations.