We are an active university Mathematics Department with a strong teaching and research reputation. We offer students the chance to study at undergraduate or postgraduate level on degree programmes leading to: BSc in Mathematics, BSc/BA joint courses in Mathematics or Applied Statistics and a wide range of other subjects. We have an active research group focusing on Computational Applied Mathematics, with research students studying for the degrees of MPhil and PhD, postdoctoral workers and associated collaborators from across the world.

Collections in this community

Recent Submissions

  • Composite Matrices from Group Rings, Composite G-Codes and Constructions of Self-Dual Codes

    Dougherty, Steven; Gildea, Joe; Korban, Adrian; Kaya, Abidin; University of Scranton; University of Chester; Sampoerna University
    In this work, we define composite matrices which are derived from group rings. We extend the idea of G-codes to composite G-codes. We show that these codes are ideals in a group ring, where the ring is a finite commutative Frobenius ring and G is an arbitrary finite group. We prove that the dual of a composite G-code is also a composite G-code. We also define quasi-composite G-codes. Additionally, we study generator matrices, which consist of the identity matrices and the composite matrices. Together with the generator matrices, the well known extension method, the neighbour method and its generalization, we find extremal binary self-dual codes of length 68 with new weight enumerators for the rare parameters $\gamma$ = 7; 8 and 9: In particular, we find 49 new such codes. Moreover, we show that the codes we find are inaccessible from other constructions.
  • High order algorithms for numerical solution of fractional differential equations

    Asl, Mohammad Shahbazi; Javidi, Mohammad; Yan, Yubin; University of Chester; University of Tabriz
    In this paper, two novel high order numerical algorithms are proposed for solving fractional differential equations where the fractional derivative is considered in the Caputo sense. The total domain is discretized into a set of small subdomains and then the unknown functions are approximated using the piecewise Lagrange interpolation polynomial of degree three and degree four. The detailed error analysis is presented, and it is analytically proven that the proposed algorithms are of orders 4 and 5. The stability of the algorithms is rigorously established and the stability region is also achieved. Numerical examples are provided to check the theoretical results and illustrate the efficiency and applicability of the novel algorithms.
  • G-Codes, self-dual G-Codes and reversible G-Codes over the Ring Bj,k

    Dougherty, Steven; Gildea, Joe; Korban, Adrian; Sahinkaya, Serap; Tarsus University; University of Chester
    In this work, we study a new family of rings, Bj,k, whose base field is the finite field Fpr . We study the structure of this family of rings and show that each member of the family is a commutative Frobenius ring. We define a Gray map for the new family of rings, study G-codes, self-dual G-codes, and reversible G-codes over this family. In particular, we show that the projection of a G-code over Bj,k to a code over Bl,m is also a G-code and the image under the Gray map of a self-dual G-code is also a self-dual G-code when the characteristic of the base field is 2. Moreover, we show that the image of a reversible G-code under the Gray map is also a reversible G2j+k-code. The Gray images of these codes are shown to have a rich automorphism group which arises from the algebraic structure of the rings and the groups. Finally, we show that quasi-G codes, which are the images of G-codes under the Gray map, are also Gs-codes for some s.
  • The multi-dimensional Stochastic Stefan Financial Model for a portfolio of assets

    Antonopoulou, Dimitra; Bitsaki, Marina; Karali, Georgia; University of Chester; University of Crete
    The financial model proposed in this work involves the liquidation process of a portfolio of n assets through sell or (and) buy orders placed, in a logarithmic scale, at a (vectorial) price with volatility. We present the rigorous mathematical formulation of this model in a financial setting resulting to an n-dimensional outer parabolic Stefan problem with noise. The moving boundary encloses the areas of zero trading, the so-called solid phase. We will focus on a case of financial interest when one or more markets are considered. In particular, our aim is to estimate for a short time period the areas of zero trading, and their diameter which approximates the minimum of the n spreads of the portfolio assets for orders from the n limit order books of each asset respectively. In dimensions n = 3, and for zero volatility, this problem stands as a mean field model for Ostwald ripening, and has been proposed and analyzed by Niethammer in [25], and in [7] in a more general setting. There in, when the initial moving boundary consists of well separated spheres, a first order approximation system of odes had been rigorously derived for the dynamics of the interfaces and the asymptotic pro le of the solution. In our financial case, we propose a spherical moving boundaries approach where the zero trading area consists of a union of spherical domains centered at portfolios various prices, while each sphere may correspond to a different market; the relevant radii represent the half of the minimum spread. We apply It^o calculus and provide second order formal asymptotics for the stochastic version dynamics, written as a system of stochastic differential equations for the radii evolution in time. A second order approximation seems to disconnect the financial model from the large diffusion assumption for the trading density. Moreover, we solve the approximating systems numerically.
  • Numerical approximation of the Stochastic Cahn-Hilliard Equation near the Sharp Interface Limit

    Antonopoulou, Dimitra; Banas, Lubomir; Nurnberg, Robert; Prohl, Andreas; University of Chester; University of Bielefeld; Imperial College London; University of Tuebingen
    Abstract. We consider the stochastic Cahn-Hilliard equation with additive noise term that scales with the interfacial width parameter ε. We verify strong error estimates for a gradient flow structure-inheriting time-implicit discretization, where ε only enters polynomially; the proof is based on higher-moment estimates for iterates, and a (discrete) spectral estimate for its deterministic counterpart. For γ sufficiently large, convergence in probability of iterates towards the deterministic Hele-Shaw/Mullins-Sekerka problem in the sharp-interface limit ε → 0 is shown. These convergence results are partly generalized to a fully discrete finite element based discretization. We complement the theoretical results by computational studies to provide practical evidence concerning the effect of noise (depending on its ’strength’ γ) on the geometric evolution in the sharp-interface limit. For this purpose we compare the simulations with those from a fully discrete finite element numerical scheme for the (stochastic) Mullins-Sekerka problem. The computational results indicate that the limit for γ ≥ 1 is the deterministic problem, and for γ = 0 we obtain agreement with a (new) stochastic version of the Mullins-Sekerka problem.
  • Entropy-driven cell decision-making predicts "fluid-to-solid" transition in multicellular systems

    Kavallaris, Nikos; Barua, Arnab; Syga, Simon; Mascheroni, Pietro; Meyer-Hermann, Michael; Deutsch, Andreas; Hatzikirou, Haralampos; University of Chester; Helmholtz Centre for Infection Research; Technische Univesität Dresden; Technische Universität Braunschweig; Khalifa University
    Cellular decision making allows cells to assume functionally different phenotypes in response to microenvironmental cues, with or without genetic change. It is an open question, how individual cell decisions influence the dynamics at the tissue level. Here, we study spatio-temporal pattern formation in a population of cells exhibiting phenotypic plasticity, which is a paradigm of cell decision making. We focus on the migration/resting and the migration/proliferation plasticity which underly the epithelial-mesenchymal transition (EMT) and the go or grow dichotomy. We assume that cells change their phenotype in order to minimize their microenvironmental entropy following the LEUP (Least microEnvironmental Uncertainty Principle) hypothesis. In turn, we study the impact of the LEUP-driven migration/resting and migration/proliferation plasticity on the corresponding multicellular spatiotemporal dynamics with a stochastic cell-based mathematical model for the spatio-temporal dynamics of the cell phenotypes. In the case of the go or rest plasticity, a corresponding mean-field approximation allows to identify a bistable switching mechanism between a diffusive (fluid) and an epithelial (solid) tissue phase which depends on the sensitivity of the phenotypes to the environment. For the go or grow plasticity, we show the possibility of Turing pattern formation for the "solid" tissue phase and its relation with the parameters of the LEUP-driven cell decisions.
  • Extending an Established Isomorphism between Group Rings and a Subring of the n × n Matrices

    Dougherty, Steven; Gildea, Joe; Korban, Adrian; University of Scranton; University of Chester
    In this work, we extend an established isomorphism between group rings and a subring of the n × n matrices. This extension allows us to construct more complex matrices over the ring R. We present many interesting examples of complex matrices constructed directly from our extension. We also show that some of the matrices used in the literature before can be obtained by a direct application of our extended isomorphism.
  • Two high-order time discretization schemes for subdiffusion problems with nonsmooth data

    Yan, Yubin; Wang, Yanyong; Yang, Yan; University of Chester; Lvliang University
    Two new high-order time discretization schemes for solving subdiffusion problems with nonsmooth data are developed based on the corrections of the existing time discretization schemes in literature. Without the corrections, the schemes have only a first order of accuracy for both smooth and nonsmooth data. After correcting some starting steps and some weights of the schemes, the optimal convergence orders $O(k^{3- \alpha})$ and $O(k^{4- \alpha})$ with $0< \alpha <1$ can be restored for any fixed time $t$ for both smooth and nonsmooth data, respectively. The error estimates for these two new high-order schemes are proved by using Laplace transform method for both homogeneous and inhomogeneous problem. Numerical examples are given to show that the numerical results are consistent with the theoretical results.
  • Dynamics of shadow system of a singular Gierer-Meinhardt system on an evolving domain

    Kavallaris, Nikos I.; Bareira, Raquel; Madzvamuse, Anotida; University of Chester; Polytechnic Institute of Setubal; University of Lisbon; Sussex University
    The main purpose of the current paper is to contribute towards the comprehension of the dynamics of the shadow system of a singular Gierer-Meinhardt model on an isotropically evolving domain. In the case where the inhibitor's response to the activator's growth is rather weak, then the shadow system of the Gierer-Meinhardt model is reduced to a single though non-local equation whose dynamics is thoroughly investigated throughout the manuscript. The main focus is on the derivation of blow-up results for this non-local equation, which can be interpreted as instability patterns of the shadow system. In particular, a diffusion-driven instability (DDI), or Turing instability, in the neighbourhood of a constant stationary solution, which then is destabilised via diffusion-driven blow-up, is observed. The latter indicates the formation of some unstable patterns, whilst some stability results of global-in-time solutions towards non-constant steady states guarantee the occurrence of some stable patterns. Most of the theoretical results are verified numerically, whilst the numerical approach is also used to exhibit the dynamics of the shadow system when analytical methods fail.
  • DOMestic Energy Systems and Technologies InCubator (DOMESTIC) and indoor air quality of the built environment

    Li, Jinghua; Khalid, Yousaf; Phillips, Gavin J.; University of Chester
    Oral presentation at RMetS Students and Early Career Scientists Conference 2020 on research project DOMESTIC (DOMestic Energy Systems and Technologies InCubator), which aims to build a facility for the demonstration of domestic technologies and design methodologies (i.e. air quality, energy efficiency).
  • New Self-Dual Codes of Length 68 from a 2 × 2 Block Matrix Construction and Group Rings

    Bortos, Maria; Gildea, Joe; Kaya, Abidin; Korban, Adrian; Tylyshchak, Alexander; Uzhgorod National University, University of Chester, Harmony School of Technology, University of Chester, Uzhgorod National University
    Many generator matrices for constructing extremal binary self-dual codes of different lengths have the form G = (In | A); where In is the n x n identity matrix and A is the n x n matrix fully determined by the first row. In this work, we define a generator matrix in which A is a block matrix, where the blocks come from group rings and also, A is not fully determined by the elements appearing in the first row. By applying our construction over F2 +uF2 and by employing the extension method for codes, we were able to construct new extremal binary self-dual codes of length 68. Additionally, by employing a generalised neighbour method to the codes obtained, we were able to con- struct many new binary self-dual [68,34,12]-codes with the rare parameters $\gamma = 7$; $8$ and $9$ in $W_{68,2}$: In particular, we find 92 new binary self-dual [68,34,12]-codes.
  • Self-Dual Codes using Bisymmetric Matrices and Group Rings

    Gildea, Joe; Kaya, Abidin; Korban, Adrian; Tylyshchak, Alexander; University of Chester ; Sampoerna University ; University of Chester: Uzhgorod National University (Elsevier, 2020-08-14)
    In this work, we describe a construction in which we combine together the idea of a bisymmetric matrix and group rings. Applying this construction over the ring F4 + uF4 together with the well known extension and neighbour methods, we construct new self-dual codes of length 68: In particular, we find 41 new codes of length 68 that were not known in the literature before.
  • New Extremal binary self-dual codes of length 68 from generalized neighbors

    Gildea, Joe; Kaya, Abidin; Korban, Adrian; Yildiz, Bahattin; University of Chester; Sampoerna University; Northern Arizona University
    In this work, we use the concept of distance between self-dual codes, which generalizes the concept of a neighbor for self-dual codes. Using the $k$-neighbors, we are able to construct extremal binary self-dual codes of length 68 with new weight enumerators. We construct 143 extremal binary self-dual codes of length 68 with new weight enumerators including 42 codes with $\gamma=8$ in their $W_{68,2}$ and 40 with $\gamma=9$ in their $W_{68,2}$. These examples are the first in the literature for these $\gamma$ values. This completes the theoretical list of possible values for $\gamma$ in $W_{68,2}$.
  • Higher Order Time Stepping Methods for Subdiffusion Problems Based on Weighted and Shifted Grünwald–Letnikov Formulae with Nonsmooth Data

    Yan, yubin; Wang, Yanyong; Yan, Yuyuan; Pani, Amiya K.; University of Chester, Lvliang University, Jimei University, Indian Institute of Technology Bombay (Springer Link, 2020-05-19)
    Two higher order time stepping methods for solving subdiffusion problems are studied in this paper. The Caputo time fractional derivatives are approximated by using the weighted and shifted Gr\"unwald-Letnikov formulae introduced in Tian et al. [Math. Comp. 84 (2015), pp. 2703-2727]. After correcting a few starting steps, the proposed time stepping methods have the optimal convergence orders $O(k^2)$ and $ O(k^3)$, respectively for any fixed time $t$ for both smooth and nonsmooth data. The error estimates are proved by directly bounding the approximation errors of the kernel functions. Moreover, we also present briefly the applicabilities of our time stepping schemes to various other fractional evolution equations. Finally, some numerical examples are given to show that the numerical results are consistent with the proven theoretical results.
  • An analysis of the L1 scheme for stochastic subdiffusion problem driven by integrated space-time white noise

    Yan, Yubin; Yan, Yuyuan; Wu, Xiaolei; University of Chester, Lvliang University, Jimei University (Elsevier, 2020-06-02)
    We consider the strong convergence of the numerical methods for solving stochastic subdiffusion problem driven by an integrated space-time white noise. The time fractional derivative is approximated by using the L1 scheme and the time fractional integral is approximated with the Lubich's first order convolution quadrature formula. We use the Euler method to approximate the noise in time and use the truncated series to approximate the noise in space. The spatial variable is discretized by using the linear finite element method. Applying the idea in Gunzburger \et (Math. Comp. 88(2019), pp. 1715-1741), we express the approximate solutions of the fully discrete scheme by the convolution of the piecewise constant function and the inverse Laplace transform of the resolvent related function. Based on such convolution expressions of the approximate solutions, we obtain the optimal convergence orders of the fully discrete scheme in spatial multi-dimensional cases by using the Laplace transform method and the corresponding resolvent estimates.
  • New binary self-dual codes via a generalization of the four circulant construction

    Gildea, Joe; Kaya, Abidin; Yildiz, Bahattin; University of Chester ; Sampoerna University ; Northern Arizona University (Croatian Mathematical Society, 2020-05-31)
    In this work, we generalize the four circulant construction for self-dual codes. By applying the constructions over the alphabets $\mathbb{F}_2$, $\mathbb{F}_2+u\mathbb{F}_2$, $\mathbb{F}_4+u\mathbb{F}_4$, we were able to obtain extremal binary self-dual codes of lengths 40, 64 including new extremal binary self-dual codes of length 68. More precisely, 43 new extremal binary self-dual codes of length 68, with rare new parameters have been constructed.
  • 2^n Bordered Constructions of Self-Dual codes from Group Rings

    Dougherty, Steven; Gildea, Joe; Kaya, Abidin; University of Scranton; University of Chester; Sampoerna Academy (Elsevier, 2020-08-04)
    Self-dual codes, which are codes that are equal to their orthogonal, are a widely studied family of codes. Various techniques involving circulant matrices and matrices from group rings have been used to construct such codes. Moreover, families of rings have been used, together with a Gray map, to construct binary self-dual codes. In this paper, we introduce a new bordered construction over group rings for self-dual codes by combining many of the previously used techniques. The purpose of this is to construct self-dual codes that were missed using classical construction techniques by constructing self-dual codes with different automorphism groups. We apply the technique to codes over finite commutative Frobenius rings of characteristic 2 and several group rings and use these to construct interesting binary self-dual codes. In particular, we construct some extremal self-dual codes length 64 and 68, constructing 30 new extremal self-dual codes of length 68.
  • Finite-time blow-up of a non-local stochastic parabolic problem

    Kavallaris, Nikos I.; Yan, Yubin; University of Chester (Elsevier, 2020-04-13)
    The main aim of the current work is the study of the conditions under which (finite-time) blow-up of a non-local stochastic parabolic problem occurs. We first establish the existence and uniqueness of the local-in-time weak solution for such problem. The first part of the manuscript deals with the investigation of the conditions which guarantee the occurrence of noise-induced blow-up. In the second part we first prove the $C^{1}$-spatial regularity of the solution. Then, based on this regularity result, and using a strong positivity result we derive, for first in the literature of SPDEs, a Hopf's type boundary value point lemma. The preceding results together with Kaplan's eigenfunction method are then employed to provide a (non-local) drift term induced blow-up result. In the last part of the paper, we present a method which provides an upper bound of the probability of (non-local) drift term induced blow-up.
  • New Extremal Self-Dual Binary Codes of Length 68 via Composite Construction, F2 + uF2 Lifts, Extensions and Neighbors

    Dougherty, Steven; Gildea, Joe; Korban, Adrian; Kaya, Abidin; University of Scranton; University of Chester; University of Chester; Sampoerna Academy; (Inderscience, 2020-02-29)
    We describe a composite construction from group rings where the groups have orders 16 and 8. This construction is then applied to find the extremal binary self-dual codes with parameters [32, 16, 8] or [32, 16, 6]. We also extend this composite construction by expanding the search field which enables us to find more extremal binary self-dual codes with the above parameters and with different orders of automorphism groups. These codes are then lifted to F2 + uF2, to obtain extremal binary images of codes of length 64. Finally, we use the extension method and neighbor construction to obtain new extremal binary self-dual codes of length 68. As a result, we obtain 28 new codes of length 68 which were not known in the literature before.
  • The diffusion-driven instability and complexity for a single-handed discrete Fisher equation

    Yan, Yubin; Zhang, Guang; Zhang, Ruixuan; University of Chester; Tianjin University of Commerce (Elsevier, 2019-12-19)
    For a reaction diffusion system, it is well known that the diffusion coefficient of the inhibitor must be bigger than that of the activator when the Turing instability is considered. However, the diffusion-driven instability/Turing instability for a single-handed discrete Fisher equation with the Neumann boundary conditions may occur and a series of 2-periodic patterns have been observed. Motivated by these pattern formations, the existence of 2-periodic solutions is established. Naturally, the periodic double and the chaos phenomenon should be considered. To this end, a simplest two elements system will be further discussed, the flip bifurcation theorem will be obtained by computing the center manifold, and the bifurcation diagrams will be simulated by using the shooting method. It proves that the Turing instability and the complexity of dynamical behaviors can be completely driven by the diffusion term. Additionally, those effective methods of numerical simulations are valid for experiments of other patterns, thus, are also beneficial for some application scientists.

View more