Due to the increasing application of fractional calculus in engineering and biomedical processes, we analyze a new method for the numerical simulation of a large class of coupled systems of fractional-order partial differential equations. In this paper, we study shifted Jacobi polynomials in the case of two variables and develop some new operational matrices of fractional-order integrations as well as fractional-order differentiations. By the use of these operational matrices, we present a new and easy method for solving a generalized class of coupled systems of fractional-order partial differential equations subject to some initial conditions. We convert the system under consideration to a system of easily solvable algebraic equation without discretizing the system, and obtain a highly accurate solution. Also, the proposed method is compared with some other well-known differential transform methods. The proposed method is computer oriented. We use MatLab to perform the necessary calculation. The next two parts will appear soon.

This article is intended to fill in the blank of the numerical schemes with second-order convergence accuracy in time for nonlinear Stokes’ first problem for a heated generalized second grade fluid with fractional derivative. A linearized difference scheme is proposed. The time fractional-order derivative is discretized by second-order shifted and weighted Gr¨unwald-Letnikov difference operator. The convergence accuracy in space is improved by performing the average operator. The presented numerical method is unconditionally stable with the global convergence order of in maximum norm, where τ and h are the step sizes in time and space, respectively. Finally, numerical examples are carried out to verify the theoretical results, showing that our scheme is efficient indeed.

We consider a possibly anisotropic integrodifferential semilinear equation, driven by a non-decreasing nonlinearity. We prove that if the solution grows less than the order of the operator at infinity, then it must be affine (possibly constant).

This paper is concerned with the invisibility cloaking in acoustic wave scattering from a new perspective. We are especially interested in achieving the invisibility cloaking by completely regular and isotropic mediums. It is shown that an interior transmission eigenvalue problem arises in our study, which is the one considered theoretically in Cakoni et al. (Transmission eigenvalues for inhomogeneous media containing obstacles, Inverse Problems and Imaging, 6 (2012), 373–398). Based on such an observation, we propose a cloaking scheme that takes a three-layer structure including a cloaked region, a lossy layer and a cloaking shell. The target medium in the cloaked region can be arbitrary but regular, whereas the mediums in the lossy layer and the cloaking shell are both regular and isotropic. We establish that if a certain non-transparency condition is satisfied, then there exists an infinite set of incident waves such that the cloaking device is nearly invisible under the corresponding wave interrogation. The set of waves is generated from the Herglotz approximation of the associated interior transmission eigenfunctions. We provide both theoretical and numerical justifications.

We consider the unique solvability of an inverse-source problem with integral transmitting condition for a time-fractional mixed type equation in rectangular domain where the unknown source term depends only on the space variable. The solution is based on a series expansion using a bi-orthogonal basis in space, corresponding to a non-self-adjoint boundary value problem. Under certain regularity conditions on the given data, we prove the uniqueness and existence of the solution for the given problem. The influence of the transmitting condition on the solvability of the problem is also demonstrated. Two different transmitting conditions are considered — viz. a full integral form and a special case. In order to simplify the bulky expressions appearing in the proof of our main result, we establish a new property of the recently introduced Mittag-Leffler type function in two variables.

An initial-boundary value problem for a time-fractional diffusion equation is discretized in space, using continuous piecewise-linear finite elements on a domain with a re-entrant corner. Known error bounds for the case of a convex domain break down, because the associated Poisson equation is no longer $H^{2}$-regular. In particular, the method is no longer second-order accurate if quasi-uniform triangulations are used. We prove that a suitable local mesh refinement about the re-entrant corner restores second-order convergence. In this way, we generalize known results for the classical heat equation.

Given a compact Lie group $G$, in this paper we establish $L^{p}$-bounds for pseudo-differential operators in $L^{p}(G)$. The criteria here are given in terms of the concept of matrix symbols defined on the noncommutative analogue of the phase space $G\times \widehat{G}$, where $\widehat{G}$ is the unitary dual of $G$. We obtain two different types of $L^{p}$ bounds: first for finite regularity symbols and second for smooth symbols. The conditions for smooth symbols are formulated using $\mathscr{S}_{\unicode[STIX]{x1D70C},\unicode[STIX]{x1D6FF}}^{m}(G)$ classes which are a suitable extension of the well-known $(\unicode[STIX]{x1D70C},\unicode[STIX]{x1D6FF})$ ones on the Euclidean space. The results herein extend classical $L^{p}$ bounds established by C. Fefferman on $\mathbb{R}^{n}$. While Fefferman’s results have immediate consequences on general manifolds for $\unicode[STIX]{x1D70C}>\max \{\unicode[STIX]{x1D6FF},1-\unicode[STIX]{x1D6FF}\}$, our results do not require the condition $\unicode[STIX]{x1D70C}>1-\unicode[STIX]{x1D6FF}$. Moreover, one of our results also does not require $\unicode[STIX]{x1D70C}>\unicode[STIX]{x1D6FF}$. Examples are given for the case of $\text{SU}(2)\cong \mathbb{S}^{3}$ and vector fields/sub-Laplacian operators when operators in the classes $\mathscr{S}_{0,0}^{m}$ and $\mathscr{S}_{\frac{1}{2},0}^{m}$ naturally appear, and where conditions $\unicode[STIX]{x1D70C}>\unicode[STIX]{x1D6FF}$ and $\unicode[STIX]{x1D70C}>1-\unicode[STIX]{x1D6FF}$ fail, respectively.

We deal with the problem of determining an unknown part of the boundary of an electrical conductor that is inaccessible for external observation and where a corrosion process is going on. We obtain estimates of the size of this damaged region from above and below.

In Angulo-Ardoy et al. [Anal. PDE, 9(3) (2016), 575–596], we found some necessary conditions for a Riemannian manifold to admit a local limiting Carleman weight (LCW), based on the Cotton–York tensor in dimension 3 and the Weyl tensor in dimension 4. In this paper, we find further necessary conditions for the existence of local LCWs that are often sufficient. For a manifold of dimension 3 or 4, we classify the possible Cotton–York, or Weyl tensors, and provide a mechanism to find out whether the manifold admits local LCW for each type of tensor. In particular, we show that a product of two surfaces admits an LCW if and only if at least one of the two surfaces is of revolution. This provides an example of a manifold satisfying the eigenflag condition of Angulo-Ardoy et al. [Anal. PDE, 9(3) (2016), 575–596] but not admitting LCW.

In this paper we consider a system of reaction–diffusion–advection equations with a free boundary, which arises in a competition ecological model in heterogeneous environment. The evolution of the free-boundary problem is discussed, which is an extension of the results of Du and Lin (Discrete Contin. Dynam. Syst. B19 (2014), 3105–3132). Precisely, when u is an inferior competitor, we prove that (u, v) → (0, V) as t→∞. When u is a superior competitor, we prove that a spreading–vanishing dichotomy holds, namely, as t→∞, either h(t)→∞ and (u, v) → (U, 0), or limt→∞h(t) < ∞ and (u, v) → (0, V). Moreover, in a weak competition case, we prove that two competing species coexist in the long run, while in a strong competition case, two species spatially segregate as the competition rates become large. Furthermore, when spreading occurs, we obtain some rough estimates of the asymptotic spreading speed.

The computational work and storage of numerically solving the time fractional PDEs are generally huge for the traditional direct methods since they require total memory and work, where NT and NS represent the total number of time steps and grid points in space, respectively. To overcome this difficulty, we present an efficient algorithm for the evaluation of the Caputo fractional derivative of order α∈(0,1). The algorithm is based on an efficient sum-of-exponentials (SOE) approximation for the kernel t–1–α on the interval [Δt, T] with a uniform absolute error ε. We give the theoretical analysis to show that the number of exponentials Nexp needed is of order for T≫1 or for TH1 for fixed accuracy ε. The resulting algorithm requires only storage and work when numerically solving the time fractional PDEs. Furthermore, we also give the stability and error analysis of the new scheme, and present several numerical examples to demonstrate the performance of our scheme.

In this paper, we present an efficient rescaling scheme for computing the long-time dynamics of expanding interfaces. The idea is to design an adaptive time-space mapping such that in the new time scale, the interfaces evolves logarithmically fast at early growth stage and exponentially fast at later times. The new spatial scale guarantees the conservation of the area/volume enclosed by the interface. Compared with the original rescaling method in [J. Comput. Phys. 225(1) (2007) 554–567], this adaptive scheme dramatically improves the slow evolution at early times when the size of the interface is small. Our results show that the original three-week computation in [J. Comput. Phys. 225(1) (2007) 554–567] can be reproduced in about one day using the adaptive scheme. We then present the largest and most complicated Hele-Shaw simulation up to date.

In this paper, we consider the numerical method that is proposed and analyzed in [J. Cao and C. Xu, J. Comput. Phys., 238 (2013), pp. 154–168] for the fractional ordinary differential equations. It is based on the so-called block-by-block approach, which is a common method for the integral equations. We extend the technique to solve the nonlinear system of fractional ordinary differential equations (FODEs) and present a general technique to construct high order schemes for the numerical solution of the nonlinear coupled system of fractional ordinary differential equations (FODEs). By using the present method, we are able to construct a high order schema for nonlinear system of FODEs of the order α, α>0. The stability and convergence of the schema is rigorously established. Under the smoothness assumption f,g ∈ C 4[0,T], we prove that the numerical solution converges to the exact solution with order 3+α for 0<α≤1 and order 4 for α>1. Some numerical examples are provided to confirm the theoretical claims.

This paper studies a system of semi-linear fractional diffusion equations which arise in competitive predator-prey models by replacing the second-order derivatives in the spatial variables with fractional derivatives of order less than two. Moving finite element methods are proposed to solve the system of fractional diffusion equations and the convergence rates of the methods are proved. Numerical examples are carried out to confirm the theoretical findings. Some applications in anomalous diffusive Lotka-Volterra and Michaelis-Menten-Holling predator-prey models are studied.

We extend the full wavefield modeling with forward scattering theory and Volterra Renormalization to a vertically varying two-parameter (velocity and density) acoustic medium. The forward scattering series, derived by applying Born-Neumann iterative procedure to the Lippmann-Schwinger equation (LSE), is a well known tool for modeling and imaging. However, it has limited convergence properties depending on the strength of contrast between the actual and reference medium or the angle of incidence of a plane wave component. Here, we introduce the Volterra renormalization technique to the LSE. The renormalized LSE and related Neumann series are absolutely convergent for any strength of perturbation and any incidence angle. The renormalized LSE can further be separated into two sub-Volterra type integral equations, which are then solved noniteratively. We apply the approach to velocity-only, density-only, and both velocity and density perturbations. We demonstrate that this Volterra Renormalization modeling is a promising and efficient method. In addition, it can also provide insight for developing a scattering theory-based direct inversion method.

This paper deals with numerical solution to the multi-term time fractional diffusion equation in a finite domain. An implicit finite difference scheme is established based on Caputo's definition to the fractional derivatives, and the upper and lower bounds to the spectral radius of the coefficient matrix of the difference scheme are estimated, with which the unconditional stability and convergence are proved. The numerical results demonstrate the effectiveness of the theoretical analysis, and the method and technique can also be applied to other kinds of time/space fractional diffusion equations.

It is known that large time-stepping method are useful for simulating phase field models. In this work, an adaptive time-stepping strategy is proposed based on numerical energy stability and equi-distribution principle. The main idea is to use the energy variation as an indicator to update the time step, so that the resulting algorithm is free of user-defined parameters, which is different from several existing approaches. Some numerical experiments are presented to illustrate the effectiveness of the algorithms.

Consider the time-harmonic acoustic scattering from an extended elastic body surrounded by a finite number of point-like obstacles in a fluid. We assume point source waves are emitted from arrayed transducers and the signals of scattered near-field data are recorded by receivers not far away from the scatterers (compared to the incident wavelength). The forward scattering can be modeled as an interaction problem between acoustic and elastic waves together with a multiple scattering problem between the extend solid and point scatterers. We prove a necessary and sufficient condition that can be used simultaneously to recover the shape of the extended elastic solid and to locate the positions of point scatterers. The essential ingredient in our analysis is the outgoing-to-incoming (OtI) operator applied to the resulting near-field response matrix (or operator). In the first part, we justify the MUSIC algorithm for locating point scatterers from near-field measurements. In the second part, we apply the factorization method, the continuous analogue of MUSIC, to the two-scale scattering problem for determining both extended and point scatterers. Numerical examples in 2D are demonstrated to show the validity and accuracy of our inversion algorithms.

The purpose of this paper is to numerically realize the inverse scattering scheme proposed in [19] of reconstructing complex elastic objects by a single far-field measurement. The unknown elastic scatterers might consist of both rigid bodies and traction-free cavities with components of multiscale sizes presented simultaneously. We conduct extensive numerical experiments to show the effectiveness and efficiency of the imaging scheme proposed in [19]. Moreover, we develop a two-stage technique, which can significantly speed up the reconstruction to yield a fast imaging scheme.

The first-order cross correlation and corresponding applications in the passive imaging are deeply studied by Garnier and Papanicolaou in their pioneer works. In this paper, the results of the first-order cross correlation are generalized to the second-order cross correlation. The second-order cross correlation is proven to be a statistically stable quantity, with respective to the random ambient noise sources. Specially, with proper time scales, the stochastic fluctuation for the second-order cross correlation converges much faster than the first-order one. Indeed, the convergent rate is of order , with 0 < α < 1. Besides, by using the stationary phase method in both homogeneous and scattering medium, similar behaviors of the singular components for the second-order cross correlation are obtained. Finally, two imaging methods are proposed to search for a target point reflector: One method is based on the imaging function, and has a better signal-to-noise rate; Another method is based on the geometric property, and can improve the bad range resolution of the imaging results.