Browsing Mathematics by Publication Date
Now showing items 4160 of 195

Malliavin Calculus for the stochastic Cahn Hilliard/AllenCahn equation with unbounded noise diffusionThe stochastic partial di erential equation analyzed in this work, is motivated by a simplified mesoscopic physical model for phase separation. It describes pattern formation due to adsorption and desorption mechanisms involved in surface processes, in the presence of a stochastic driving force. This equation is a combination of CahnHilliard and AllenCahn type operators with a multiplicative, white, spacetime noise of unbounded di usion. We apply Malliavin calculus, in order to investigate the existence of a density for the stochastic solution u. In dimension one, according to the regularity result in [5], u admits continuous paths a.s. Using this property, and inspired by a method proposed in [8], we construct a modi ed approximating sequence for u, which properly treats the new second order AllenCahn operator. Under a localization argument, we prove that the Malliavin derivative of u exists locally, and that the law of u is absolutely continuous, establishing thus that a density exists.

A posteriori error estimates for fully discrete schemes for the time dependent Stokes problemThis work is devoted to a posteriori error analysis of fully discrete finite element approximations to the time dependent Stokes system. The space discretization is based on popular stable spaces, including Crouzeix–Raviart and Taylor–Hood finite element methods. Implicit Euler is applied for the time discretization. The finite element spaces are allowed to change with time steps and the projection steps include alternatives that is hoped to cope with possible numerical artifices and the loss of the discrete incompressibility of the schemes. The final estimates are of optimal order in L∞(L2) for the velocity error.

Fourier spectral methods for stochastic space fractional partial differential equations driven by special additive noisesFourier spectral methods for solving stochastic space fractional partial differential equations driven by special additive noises in onedimensional 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 spacetime 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.

The sharp interface limit for the stochastic CahnHilliard EquationWe study the two and three dimensional stochastic CahnHilliard equation in the sharp interface limit, where the positive parameter \eps tends to zero, which measures the width of transition layers generated during phase separation. We also couple the noise strength to this parameter. Using formal asymptotic expansions, we identify the limit. In the right scaling we indicate that the solutions of stochastic CahnHilliard converge to a solution of a HeleShaw problem with stochastic forcing. In the case when the noise is sufficiently small, we rigorously prove that the limit is a deterministic HeleShaw problem. Finally, we discuss which estimates are necessary in order to extend the rigorous result to larger noise strength.

Constructions for SelfDual Codes Induced from Group RingsIn this work, we establish a strong connection between group rings and selfdual codes. We prove that a group ring element corresponds to a selfdual code if and only if it is a unitary unit. We also show that the doublecirculant and fourcirculant constructions come from cyclic and dihedral groups, respectively. Using groups of order 8 and 16 we find many new construction methods, in addition to the wellknown methods, for selfdual codes. We establish the relevance of these new constructions by finding many extremal binary selfdual codes using them, which we list in several tables. In particular, we construct 10 new extremal binary selfdual codes of length 68.

An analysis of the modified L1 scheme for timefractional partial differential equations with nonsmooth dataWe introduce a modified L1 scheme for solving time fractional partial differential equations and obtain error estimates for smooth and nonsmooth initial data in both homogeneous and inhomogeneous cases. Jin \et (2016, An analysis of the L1 scheme for the subdiffusion equation with nonsmooth data, IMA J. of Numer. Anal., 36, 197221) established an $O(k)$ convergence rate for the L1 scheme for smooth and nonsmooth initial data for the homogeneous problem, where $k$ denotes the time step size. We show that the modified L1 scheme has convergence rate $O(k^{2\alpha}), 0< \alpha <1$ for smooth and nonsmooth initial data in both homogeneous and inhomogeneous cases. Numerical examples are given to show that the numerical results are consistent with the theoretical results.

A novel highorder algorithm for the numerical estimation of fractional differential equationsThis paper uses polynomial interpolation to design a novel highorder algorithm for the numerical estimation of fractional differential equations. The RiemannLiouville fractional derivative is expressed by using the Hadamard finitepart integral and the piecewise cubic interpolation polynomial is utilized to approximate the integral. The detailed error analysis is presented and it is established that the convergence order of the algorithm is O(h4−a). Asymptotic expansion of the error for the presented algorithm is also investigated. Some numerical examples are provided and compared with the exact solution to show that the numerical results are in well agreement with the theoretical ones and also to illustrate the accuracy and efficiency of the proposed algorithm.

A higher order numerical method for time fractional partial differential equations with nonsmooth dataGao et al. (2014) introduced a numerical scheme to approximate the Caputo fractional derivative with the convergence rate $O(k^{3\alpha}), 0< \alpha <1$ by directly approximating the integerorder derivative with some finite difference quotients in the definition of the Caputo fractional derivative, see also Lv and Xu (2016), where $k$ is the time step size. Under the assumption that the solution of the time fractional partial differential equation is sufficiently smooth, Lv and Xu (2016) proved by using energy method that the corresponding numerical method for solving time fractional partial differential equation has the convergence rate $O(k^{3\alpha}), 0< \alpha <1$ uniformly with respect to the time variable $t$. However, in general the solution of the time fractional partial differential equation has low regularity and in this case the numerical method fails to have the convergence rate $O(k^{3 \alpha}), 0 < \alpha <1$ uniformly with respect to the time variable $t$. In this paper, we first obtain a similar approximation scheme to the RiemannLiouville fractional derivative with the convergence rate $O(k^{3 \alpha}), 0 < \alpha <1$ as in Gao \et \cite{gaosunzha} (2014) by approximating the Hadamard finitepart integral with the piecewise quadratic interpolation polynomials. Based on this scheme, we introduce a time discretization scheme to approximate the time fractional partial differential equation and show by using Laplace transform methods that the time discretization scheme has the convergence rate $O(k^{3 \alpha}), 0 < \alpha <1$ for any fixed $t_{n}>0$ for smooth and nonsmooth data in both homogeneous and inhomogeneous cases. Numerical examples are given to show that the theoretical results are consistent with the numerical results.

NonLocal Partial Differential Equations for Engineering and Biology: Mathematical Modeling and AnalysisThis book presents new developments in nonlocal mathematical modeling and mathematical analysis on the behavior of solutions with novel technical tools. Theoretical backgrounds in mechanics, thermodynamics, game theory, and theoretical biology are examined in details. It starts off with a review and summary of the basic ideas of mathematical modeling frequently used in the sciences and engineering. The authors then employ a number of models in bioscience and material science to demonstrate applications, and provide recent advanced studies, both on deterministic nonlocal partial differential equations and on some of their stochastic counterparts used in engineering. Mathematical models applied in engineering, chemistry, and biology are subject to conservation laws. For instance, decrease or increase in thermodynamic quantities and nonlocal partial differential equations, associated with the conserved physical quantities as parameters. These present novel mathematical objects are engaged with rich mathematical structures, in accordance with the interactions between species or individuals, selforganization, pattern formation, hysteresis. These models are based on various laws of physics, such as mechanics of continuum, electromagnetic theory, and thermodynamics. This is why many areas of mathematics, calculus of variation, dynamical systems, integrable systems, blowup analysis, and energy methods are indispensable in understanding and analyzing these phenomena. This book aims for researchers and upper grades students in mathematics, engineering, physics, economics, and biology.

Group Rings, GCodes and Constructions of SelfDual and Formally SelfDual CodesWe describe Gcodes, 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 Gcode is also a Gcode. We give constructions of selfdual and formally selfdual 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 selfdual codes obtained by the construction. We show that several of the standard constructions of selfdual 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 selfdual 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 quasiG codes and give a construction of these codes.

An approach to construct higher order time discretisation schemes for time fractional partial differential equations with nonsmooth dataIn this paper, we shall review an approach by which we can seek higher order time discretisation schemes for solving time fractional partial differential equations with nonsmooth data. The low regularity of the solutions of time fractional partial differential equations implies standard time discretisation schemes only yield first order accuracy. To obtain higher order time discretisation schemes when the solutions of time fractional partial differential equations have low regularities, one may correct the starting steps of the standard time discretisation schemes to capture the singularities of the solutions. We will consider these corrections of some higher order time discretisation schemes obtained by using Lubich's fractional multistep methods, L1 scheme and its modification, discontinuous Galerkin methods, etc. Numerical examples are given to show that the theoretical results are consistent with the numerical results.

A multispecies chemotaxis system: Lyapunov functionals, duality, critical massWe introduce a multispecies chemotaxis type system admitting an arbitrarily large number of population species, all of which are attracted vs.\ repelled by a single chemical substance. The production vs.\ destruction rates of the chemotactic substance by the species is described by a probability measure. For such a model we investigate the variational structures, in particular we prove the existence of Lyapunov functionals, we establish duality properties as well as a logarithmic HardyLittlewoodSobolev type inequality for the associated free energy. The latter inequality provides the optimal critical value for the conserved total population mass.

A note on finite difference methods for nonlinear fractional differential equations with nonuniform meshesWe consider finite difference methods for solving nonlinear fractional differential equations in the Caputo fractional derivative sense with nonuniform meshes. Under the assumption that the Caputo derivative of the solution of the fractional differential equation is suitably smooth, Li et al. \lq \lq Finite difference methods with nonuniform meshes for nonlinear fractional differential equations\rq\rq, Journal of Computational Physics, 316(2016), 614631, obtained the error estimates of finite difference methods with nonuniform meshes. However the Caputo derivative of the solution of the fractional differential equation in general has a weak singularity near the initial time. In this paper, we obtain the error estimates of finite difference methods with nonuniform meshes when the Caputo fractional derivative of the solution of the fractional differential equation has lower smoothness. The convergence result shows clearly how the regularity of the Caputo fractional derivative of the solution affect the order of convergence of the finite difference methods. Numerical results are presented that confirm the sharpness of the error analysis.

Detailed error analysis for a fractional adams method with graded meshesWe consider a fractional Adams method for solving the nonlinear fractional differential equation $\, ^{C}_{0}D^{\alpha}_{t} y(t) = f(t, y(t)), \, \alpha >0$, equipped with the initial conditions $y^{(k)} (0) = y_{0}^{(k)}, k=0, 1, \dots, \lceil \alpha \rceil 1$. Here $\alpha$ may be an arbitrary positive number and $ \lceil \alpha \rceil$ denotes the smallest integer no less than $\alpha$ and the differential operator is the Caputo derivative. Under the assumption $\, ^{C}_{0}D^{\alpha}_{t} y \in C^{2}[0, T]$, Diethelm et al. \cite[Theorem 3.2]{dieforfre} introduced a fractional Adams method with the uniform meshes $t_{n}= T (n/N), n=0, 1, 2, \dots, N$ and proved that this method has the optimal convergence order uniformly in $t_{n}$, that is $O(N^{2})$ if $\alpha > 1$ and $O(N^{1\alpha})$ if $\alpha \leq 1$. They also showed that if $\, ^{C}_{0}D^{\alpha}_{t} y(t) \notin C^{2}[0, T]$, the optimal convergence order of this method cannot be obtained with the uniform meshes. However, it is well known that for $y \in C^{m} [0, T]$ for some $m \in \mathbb{N}$ and $ 0 < \alpha 1$, we show that the optimal convergence order of this method can be recovered uniformly in $t_{n}$ even if $\, ^{C}_{0}D^{\alpha}_{t} y$ behaves as $t^{\sigma}, 0< \sigma <1$. Numerical examples are given to show that the numerical results are consistent with the theoretical results.

Some time stepping methods for fractional diffusion problems with nonsmooth dataWe consider error estimates for some time stepping methods for solving fractional diffusion problems with nonsmooth data in both homogeneous and inhomogeneous cases. McLean and Mustapha \cite{mclmus} (Timestepping error bounds for fractional diffusion problems with nonsmooth initial data, Journal of Computational Physics, 293(2015), 201217) established an $O(k)$ convergence rate for the piecewise constant discontinuous Galerkin method with nonsmooth initial data for the homogeneous problem when the linear operator $A$ is assumed to be selfadjoint, positive semidefinite and densely defined in a suitable Hilbert space, where $k$ denotes the time step size. In this paper, we approximate the RiemannLiouville fractional derivative by Diethelm's method (or $L1$ scheme) and obtain the same time discretisation scheme as in McLean and Mustapha \cite{mclmus}. We first prove that this scheme has also convergence rate $O(k)$ with nonsmooth initial data for the homogeneous problem when $A$ is a closed, densely defined linear operator satisfying some certain resolvent estimates. We then introduce a new time discretization scheme for the homogeneous problem based on the convolution quadrature and prove that the convergence rate of this new scheme is $O(k^{1+ \alpha}), 0<\alpha <1 $ with the nonsmooth initial data. Using this new time discretization scheme for the homogeneous problem, we define a time stepping method for the inhomogeneous problem and prove that the convergence rate of this method is $O(k^{1+ \alpha}), 0<\alpha <1 $ with the nonsmooth data. Numerical examples are given to show that the numerical results are consistent with the theoretical results.

On a degenerate nonlocal parabolic problem describing infinite dimensional replicator dynamicsWe establish the existence of locally positive weak solutions to the homogeneous Dirichlet problem for \[ u_t = u \Delta u + u \int_\Omega \nabla u^2 \] in bounded domains $\Om\sub\R^n$ which arises in game theory. We prove that solutions converge to $0$ if the initial mass is small, whereas they undergo blowup in finite time if the initial mass is large. In particular, it is shown that in this case the blowup set coincides with $\overline{\Omega}$, i.e. the finitetime blowup is global.

On the dynamics of a nonlocal parabolic equation arising from the GiererMeinhardt systemThe purpose of the current paper is to contribute to the comprehension of the dynamics of the shadow system of an activatorinhibitor system known as a GiererMeinhardt model. Shadow systems are intended to work as an intermediate step between single equations and reactiondiffusion systems. In the case where the inhibitor's response to the activator's growth is rather weak, then the shadow system of the GiererMeinhardt model is reduced to a single though nonlocal equation whose dynamics will be investigated. We mainly focus on the derivation of blowup results for this nonlocal equation which can be seen as instability patterns of the shadow system. In particular, a {\it diffusion driven instability (DDI)}, or {\it Turing instability}, in the neighbourhood of a constant stationary solution, which it is destabilised via diffusiondriven blowup, is obtained. The latter actually indicates the formation of some unstable patterns, whilst some stability results of globalintime solutions towards nonconstant steady states guarantee the occurrence of some stable patterns.

Discontinuous Galerkin time stepping method for solving linear space fractional partial differential equationsIn this paper, we consider the discontinuous Galerkin time stepping method for solving the linear space fractional partial differential equations. The space fractional derivatives are defined by using Riesz fractional derivative. The space variable is discretized by means of a Galerkin finite element method and the time variable is discretized by the discontinuous Galerkin method. The approximate solution will be sought as a piecewise polynomial function in $t$ of degree at most $q1, q \geq 1$, which is not necessarily continuous at the nodes of the defining partition. The error estimates in the fully discrete case are obtained and the numerical examples are given.

Blending loworder stabilised finite element methods: a positivity preserving local projection method for the convectiondiffusion equationIn this work we propose a nonlinear blending of two loworder stabilisation mechanisms for the convection–diffusion equation. The motivation for this approach is to preserve monotonicity without sacrificing accuracy for smooth solutions. The approach is to blend a firstorder artificial diffusion method, which will be active only in the vicinity of layers and extrema, with an optimal order local projection stabilisation method that will be active on the smooth regions of the solution. We prove existence of discrete solutions, as well as convergence, under appropriate assumptions on the nonlinear terms, and on the exact solution. Numerical examples show that the discrete solution produced by this method remains within the bounds given by the continuous maximum principle, while the layers are not smeared significantly.

HighOrder Numerical Methods for Solving Time Fractional Partial Differential EquationsIn this paper we introduce a new numerical method for solving time fractional partial differential equation. The time discretization is based on Diethelm’s method where the Hadamard finitepart integral is approximated by using the piecewise quadratic interpolation polynomials. The space discretization is based on the standard finite element method. The error estimates with the convergence order O(τ^(3−α) +h^2 ),0