• Edge-based nonlinear diffusion for finite element approximations of convection–diffusion equations and its relation to algebraic flux-correction schemes

      Barrenechea, Gabriel; Burman, Erik; Karakatsani, Fotini; University of Strathclyde; University College London; University of Chester (Springer, 2016-05-07)
      For the case of approximation of convection–diffusion equations using piecewise affine continuous finite elements a new edge-based nonlinear diffusion operator is proposed that makes the scheme satisfy a discrete maximum principle. The diffusion operator is shown to be Lipschitz continuous and linearity preserving. Using these properties we provide a full stability and error analysis, which, in the diffusion dominated regime, shows existence, uniqueness and optimal convergence. Then the algebraic flux correction method is recalled and we show that the present method can be interpreted as an algebraic flux correction method for a particular definition of the flux limiters. The performance of the method is illustrated on some numerical test cases in two space dimensions.
    • 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.
    • Error estimates of high-order numerical methods for solving time fractional partial differential equations

      Li, Zhiqiang; Yan, Yubin; Luliang University; Shanghai University; University of Chester (De Gruyter, 2018-07-12)
      Error estimates of some high-order numerical methods for solving time fractional partial differential equations are studied in this paper. We first provide the detailed error estimate of a high-order numerical method proposed recently by Li et al. \cite{liwudin} for solving time fractional partial differential equation. We prove that this method has the convergence order $O(\tau^{3- \alpha})$ for all $\alpha \in (0, 1)$ when the first and second derivatives of the solution are vanish at $t=0$, where $\tau$ is the time step size and $\alpha$ is the fractional order in the Caputo sense. We then introduce a new time discretization method for solving time fractional partial differential equations, which has no requirements for the initial values as imposed in Li et al. \cite{liwudin}. We show that this new method also has the convergence order $O(\tau^{3- \alpha})$ for all $\alpha \in (0, 1)$. The proofs of the error estimates are based on the energy method developed recently by Lv and Xu \cite{lvxu}. We also consider the space discretization by using the finite element method. Error estimates with convergence order $O(\tau^{3- \alpha} + h^2)$ are proved in the fully discrete case, where $h$ is the space step size. Numerical examples in both one- and two-dimensional cases are given to show that the numerical results are consistent with the theoretical results.
    • Error estimates of a high order numerical method for solving linear fractional differential equations

      Li, Zhiqiang; Yan, Yubin; Ford, Neville J.; Luliang University; University of Chester (Elsevier, IMACS, 2016-04-29)
      In this paper, we first introduce an alternative proof of the error estimates of the numerical methods for solving linear fractional differential equations proposed in Diethelm [6] where a first-degree compound quadrature formula was used to approximate the Hadamard finite-part integral and the convergence order of the proposed numerical method is O(∆t 2−α ), 0 < α < 1, where α is the order of the fractional derivative and ∆t is the step size. We then use the similar idea to prove the error estimates of a high order numerical method for solving linear fractional differential equations proposed in Yan et al. [37], where a second-degree compound quadrature formula was used to approximate the Hadamard finite-part integral and we show that the convergence order of the numerical method is O(∆t 3−α ), 0 < α < 1. The numerical examples are given to show that the numerical results are consistent with the theoretical results.
    • Existence and regularity of solution for a Stochastic CahnHilliard / Allen-Cahn equation with unbounded noise diffusion

      Antonopoulou, Dimitra; Karali, Georgia D.; Millet, Annie; University of Chester (Elsevier, 2015-10-24)
      The Cahn-Hilliard/Allen-Cahn equation with noise is a simplified mean field model of stochastic microscopic dynamics associated with adsorption and desorption-spin flip mechanisms in the context of surface processes. For such an equation we consider a multiplicative space-time white noise with diffusion coefficient of linear growth. Applying technics from semigroup theory, we prove local existence and uniqueness in dimensions d = 1,2,3. Moreover, when the diffusion coefficient satisfies a sub-linear growth condition of order α bounded by 1 3, which is the inverse of the polynomial order of the nonlinearity used, we prove for d = 1 global existence of solution. Path regularity of stochastic solution, depending on that of the initial condition, is obtained a.s. up to the explosion time. The path regularity is identical to that proved for the stochastic Cahn-Hilliard equation in the case of bounded noise diffusion. Our results are also valid for the stochastic Cahn-Hilliard equation with unbounded noise diffusion, for which previous results were established only in the framework of a bounded diffusion coefficient. As expected from the theory of parabolic operators in the sense of Petrovsk˘ıı, the bi-Laplacian operator seems to be dominant in the combined model.
    • Existence of time periodic solutions for a class of non-resonant discrete wave equations

      Zhang, Guang; Feng, Wenying; Yan, Yubin; University of Chester (Springer, 2015-04-17)
      In this paper, a class of discrete wave equations with Dirichlet boundary conditions are obtained by using the center-difference method. For any positive integers m and T, when the existence of time mT-periodic solutions is considered, a strongly indefinite discrete system needs to be established. By using a variant generalized weak linking theorem, a non-resonant superlinear (or superquadratic) result is obtained and the Ambrosetti-Rabinowitz condition is improved. Such a method cannot be used for the corresponding continuous wave equations or the continuous Hamiltonian systems; however, it is valid for some general discrete Hamiltonian systems.
    • Existence theory for a class of evolutionary equations with time-lag, studied via integral equation formulations

      Baker, Christopher T. H.; Lumb, Patricia M.; University of Chester (University of Chester, 2006)
      In discussions of certain neutral delay differential equations in Hale’s form, the relationship of the original problem with an integrated form (an integral equation) proves to be helpful in considering existence and uniqueness of a solution and sensitivity to initial data. Although the theory is generally based on the assumption that a solution is continuous, natural solutions of neutral delay differential equations of the type considered may be discontinuous. This difficulty is resolved by relating the discontinuous solution to its restrictions on appropriate (half-open) subintervals where they are continuous and can be regarded as solutions of related integral equations. Existence and unicity theories then follow. Furthermore, it is seen that the discontinuous solutions can be regarded as solutions in the sense of Caratheodory (where this concept is adapted from the theory of ordinary differential equations, recast as integral equations).
    • Explosive solutions of a stochastic non-local reaction–diffusion equation arising in shear band formation

      Kavallaris, Nikos I.; University of Chester (Wiley, 2015-07-07)
      In this paper, we consider a non-local stochastic parabolic equation which actually serves as a mathematical model describing the adiabatic shear-banding formation phenomena in strained metals. We first present the derivation of the mathematical model. Then we investigate under which circumstances a finite-time explosion for this non-local SPDE, corresponding to shear-banding formation, occurs. For that purpose some results related to the maximum principle for this non-local SPDE are derived and afterwards the Kaplan's eigenfunction method is employed.
    • Exponential stability in p-th mean of solutions, and of convergent Euler-type solutions, of stochastic delay differential equations

      Baker, Christopher T. H.; Buckwar, Evelyn; Univesity College Chester ; Humboldt-Universität zu Berlin (Elsevier, 2005-12-15)
      This article carries out an analysis which proceeds as follows: showing that an inequality of Halanay type (derivable via comparison theory) can be employed to derive conditions for p-th mean stability of a solution; producing a discrete analogue of the Halanay-type theory, that permits the development of a p-th mean stability analysis of analogous stochastic difference equations. The application of the theoretical results is illustrated by deriving mean-square stability conditions for solutions and numerical solutions of a constant-coefficient linear test equation.
    • 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.
    • Finite Difference Method for Two-Sided Space-Fractional Partial Differential Equations

      Pal, Kamal; Liu, Fang; Yan, Yubin; Roberts, Graham; University of Chester (Springer International Publishing, 2015-06-17)
      Finite difference methods for solving two-sided space-fractional partial differential equations are studied. The space-fractional derivatives are the left-handed and right-handed Riemann-Liouville fractional derivatives which are expressed by using Hadamard finite-part integrals. The Hadamard finite-part integrals are approximated by using piecewise quadratic interpolation polynomials and a numerical approximation scheme of the space-fractional derivative with convergence order O(Δx^(3−α )),10 , where Δt,Δx denote the time and space step sizes, respectively. Numerical examples are presented and compared with the exact analytical solution for its order of convergence.
    • A finite element method for time fractional partial differential equations

      Ford, Neville J.; Xiao, Jingyu; Yan, Yubin; University of Chester ; Harbin Institute of Technology ; University of Chester (2011-09-01)
      This article considers the finite element method for time fractional differential equations.
    • 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.
    • Fixed point theroms and their application - discrete Volterra applications

      Baker, Christopher T. H.; Song, Yihong; University of Chester (University of Chester, 2006)
      The existence of solutions of nonlinear discrete Volterra equations is established. We define discrete Volterra operators on normed spaces of infinite sequences of finite-dimensional vectors, and present some of their basic properties (continuity, boundedness, and representation). The treatment relies upon the use of coordinate functions, and the existence results are obtained using fixed point theorems for discrete Volterra operators on infinite-dimensional spaces based on fixed point theorems of Schauder, Rothe, and Altman, and Banach’s contraction mapping theorem, for finite-dimensional spaces.
    • Fourier spectral methods for some linear stochastic space-fractional partial differential equations

      Liu, Yanmei; Khan, Monzorul; Yan, Yubin; LuLiang University; University of Chester (MDPI, 2016-07-01)
      Fourier spectral methods for solving some linear stochastic space-fractional partial differential equations perturbed by space-time white noises in one-dimensional case are introduced and analyzed. The space-fractional derivative is defined by using the eigenvalues and eigenfunctions of Laplacian subject to some boundary conditions. We approximate the space-time white noise by using piecewise constant functions and obtain the approximated stochastic space-fractional partial differential equations. The approximated stochastic space-fractional partial differential equations are then solved by using Fourier spectral methods. Error estimates in $L^{2}$- norm are obtained. Numerical examples are given.
    • Fourier spectral methods for stochastic space fractional partial differential equations driven by special additive noises

      Liu, Fang; Yan, Yubin; Khan, Monzorul; Lvliang University, University of Chester (EudoxusPress, 2018-02-28)
      Fourier spectral methods for solving stochastic space fractional partial differential equations driven by special additive noises in one-dimensional case are introduced and analyzed. The space fractional derivative is defined by using the eigenvalues and eigenfunctions of Laplacian subject to some boundary conditions. The space-time noise is approximated by the piecewise constant functions in the time direction and by some appropriate approximations in the space direction. The approximated stochastic space fractional partial differential equations are then solved by using Fourier spectral methods. For the linear problem, we obtain the precise error estimates in the $L_{2}$ norm and find the relations between the error bounds and the fractional powers. For the nonlinear problem, we introduce the numerical algorithms and MATLAB codes based on the FFT transforms. Our numerical algorithms can be adapted easily to solve other stochastic space fractional partial differential equations with multiplicative noises. Numerical examples for the semilinear stochastic space fractional partial differential equations are given.
    • Fractional boundary value problems: Analysis and numerical methods

      Ford, Neville J.; Morgado, Maria L.; University of Chester ; University of Tras-os-Montes e Alto Douro (Springer, 2011-07-28)
      This journal article discusses nonlinear boundary value problems.
    • Fractional pennes' bioheat equation: Theoretical and numerical studies

      Ferras, Luis L.; Ford, Neville J.; Morgado, Maria L.; Rebelo, Magda S.; Nobrega, Joao M.; University of Minho & University of Chester, University of Chester, UTAD, UNL Lisboa, University of Minho (de Gruyter, 2015-08-04)
      In this work we provide a new mathematical model for the Pennes’ bioheat equation, assuming a fractional time derivative of single order. Alternative versions of the bioheat equation are studied and discussed, to take into account the temperature-dependent variability in the tissue perfusion, and both finite and infinite speed of heat propagation. The proposed bio heat model is solved numerically using an implicit finite difference scheme that we prove to be convergent and stable. The numerical method proposed can be applied to general reaction diffusion equations, with a variable diffusion coefficient. The results obtained with the single order fractional model, are compared with the original models that use classical derivatives.
    • G-codes over Formal Power Series Rings and Finite Chain Rings

      Dougherty, Steven; Gildea, Joe; Korban, Adrian; University of Scranton; University of Chester (2020-02-29)
      In this work, we define $G$-codes over the infinite ring $R_\infty$ as ideals in the group ring $R_\infty G$. We show that the dual of a $G$-code is again a $G$-code in this setting. We study the projections and lifts of $G$-codes over the finite chain rings and over the formal power series rings respectively. We extend known results of constructing $\gamma$-adic codes over $R_\infty$ to $\gamma$-adic $G$-codes over the same ring. We also study $G$-codes over principal ideal rings.
    • 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.