• 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.
    • 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.
    • 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; Yan, Yubin; University of Chester
      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.
    • Galerkin methods for a Schroedinger-type equation with a dynamical boundary condition in two dimensions

      Antonopoulou, Dimitra; University of Chester (EDP Sciences / SMAI, 2015-06-30)
      In this paper, we consider a two-dimensional Schodinger-type equation with a dynamical boundary condition. This model describes the long-range sound propagation in naval environments of variable rigid bottom topography. Our choice for a regular enough finite element approximation is motivated by the dynamical condition and therefore, consists of a cubic splines implicit Galerkin method in space. Furthermore, we apply a Crank-Nicolson time stepping for the evolutionary variable. We prove existence and stability of the semidiscrete and fully discrete solution.
    • A genetic-algorithm approach to simulating human immunodeficiency virus evolution reveals the strong impact of multiply infected cells and recombination

      Bocharov, Gennady; Ford, Neville J.; Edwards, John T.; Breinig, Tanja; Wain-Hobson, Simon; Meyerhans, Andreas; Institute of Numerical Mathematics, Russian Academy of Sciences ; University of Chester ; University of Chester ; University of the Saarland ; Unité de Rétrovirologie Moléculaire, Institut Pasteur ; University of the Saarland (Society for General Microbiology / High Wire Press, 2005-11-01)
      It has been previously shown that the majority of human immunodeficiency virus type 1 (HIV-1)-infected splenocytes can harbour multiple, divergent proviruses with a copy number ranging from one to eight. This implies that, besides point mutations, recombination should be considered as an important mechanism in the evolution of HIV within an infected host. To explore in detail the possible contributions of multi-infection and recombination to HIV evolution, the effects of major microscopic parameters of HIV replication (i.e. the point-mutation rate, the crossover number, the recombination rate and the provirus copy number) on macroscopic characteristics (such as the Hamming distance and the abundance of n-point mutants) have been simulated in silico. Simulations predict that multiple provirus copies per infected cell and recombination act in synergy to speed up the development of sequence diversity. Point mutations can be fixed for some time without fitness selection. The time needed for the selection of multiple mutations with increased fitness is highly variable, supporting the view that stochastic processes may contribute substantially to the kinetics of HIV variation in vivo.
    • Group Rings, G-Codes and Constructions of Self-Dual and Formally Self-Dual Codes

      Dougherty, Steven; Gildea, Joe; Taylor, Rhian; Tylyshchak, Alexander; University of Scranton; University of Chester; Uzhgorod State University (Springer, 2017-11-15)
      We describe G-codes, which are codes that 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 G-code is also a G-code. We give constructions of self-dual and formally self-dual codes in this setting and we improve the existing construction given in [13] by showing that one of the conditions given in the theorem is unnecessary and, moreover, it restricts the number of self-dual codes obtained by the construction. We show that several of the standard constructions of self-dual codes are found within our general framework. We prove that our constructed codes must have an automorphism group that contains G as a subgroup. We also prove that a common construction technique for producing self-dual codes cannot produce the putative [72, 36, 16] Type II code. Additionally, we show precisely which groups can be used to construct the extremal Type II codes over length 24 and 48. We define quasi-G codes and give a construction of these codes.