A fourth‐order compact algorithm for nonlinear reaction‐diffusion equations with Neumann boundary conditions

In this article, we discuss a scheme for dealing with Neumann and mixed boundary conditions using a compact stencil. The resulting compact algorithm for solving systems of nonlinear reaction‐diffusion equations is fourth‐order accurate in both the temporal and spatial dimensions. We also prove that the standard second‐order approximation to zero Neumann boundary conditions provides fourth‐order accuracy when the nonlinear reaction term is independent of the spatial variables. Numerical examples, including an application of this algorithm to a mathematical model describing frontal polymerization process, are presented in the article to demonstrate the accuracy and efficiency of the scheme. © 2005 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2005     

A numerical solution for advection‐diffusion equation using dual reciprocity method

This article describes a numerical method based on the boundary integral equation and dual reciprocity method(DRM) for solving the one‐dimensional advection‐diffusion equations. The concept of DRM is used to convert the domain integral to the boundary that leads to an integration free method. The time derivative is approximated by the time‐stepping method. Numerical results are presented for some problems to demonstrate the usefulness and accuracy of the new approach. © 2012 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2013     

A system of boundary integral equations for the transmission problem in acoustics

In this paper, we reduce the classical two-dimensional transmission problem in acoustic scattering to a system of coupled boundary integral equations (BIEs), and consider the weak formulation of the resulting equations. Uniqueness and existence results for the weak solution of corresponding variational equations are established. In contrast to the coupled system in Costabel and Stephan (1985) [4], we need to take into account exceptional frequencies to obtain the unique solvability. Boundary element methods (BEM) based on both the standard and a two-level fast multipole Galerkin schemes are employed to compute the solution of the variational equation. Numerical results are presented to verify the efficiency and accuracy of the numerical methods.     

A class of nonlinear singularly perturbed initial boundary value problems for reaction diffusion equations with boundary perturbation

A class of nonlinear initial boundary value problems for reaction diffusion equations with boundary perturbation is considered. Under suitable conditions and using the theory of differential inequalities the asymptotic solution of the initial boundary value problems is studied.     

Critical exponents for nonlinear diffusion equations with nonlinear boundary sources

This paper is concerned with the large time behavior of solutions to two types of nonlinear diffusion equations with nonlinear boundary sources on the exterior domain of the unit ball. We are interested in the critical global exponent q₀ and the critical Fujita exponent qc for the problems considered, and show that q₀=qc for the multi-dimensional porous medium equation and non-Newtonian filtration equation with nonlinear boundary sources. This is quite different from the known results that q₀<qc for the one-dimensional case.     

A coupled continuous-discontinuous FEM approach for convection diffusion equations

In this article, we introduce a coupled approach of local discontinuous Galerkin and standard finite element method for solving convection diffusion problems. The whole domain is divided into two disjoint subdomains. The discontinuous Galerkin method is adopted in the subdomain where the solution varies rapidly, while the standard finite element method is used in the other subdomain due to its lower computational cost. The stability and a priori error estimate are established. We prove that the coupled method has O((ε1 / 2 + h 1 / 2 )h k ) convergence rate in an associated norm, where ε is the diffusion coefficient, h is the mesh size and k is the degree of polynomial. The numerical results verify our theoretical results. Moreover, 2k-order superconvergence of the numerical traces at the nodes, and the optimal convergence of the errors under L 2 norm are observed numerically on the uniform mesh. The numerical results also indicate that the coupled method has the same convergence order and almost the same errors as the purely LDG method.     

Approximate controllability for fractional diffusion equations by interior control

We consider the approximate controllability by interior control of a partial differential equation with time derivatives of non-integer order. First, we prove the well-posedness of the inhomogeneous problem for the controlled system. Next, we study the dual system and show a weak type of unique continuation property. Finally, we prove the approximate controllability.     

An inverse diffusion problem with nonlocal boundary conditions

This article considers the inverse problem of identification of a time‐dependent thermal diffusivity together with the temperature in an one‐dimensional heat equation with nonlocal boundary and integral overdetermination conditions when a heat exchange takes place across boundary of the material. The well‐posedness of the problem is studied under some regularity, and consistency conditions on the data of the problem together with the nonnegativity condition on the Fourier coefficients of the initial data and source term. The inverse problem is also studied numerically by using the Crank–Nicolson finite difference scheme combined with predictor‐corrector technique. The numerical examples are presented and discussed. © 2015 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 32: 564–590, 2016     

Numerical stability of the BEM for advection‐diffusion problems

A boundary element method (BEM) approach has been developed to solve the time‐dependent 1D advection‐diffusion equation. The 1D solution is part of a 3D numerical scheme for solving advection‐diffusion (AD) problems in fractured porous media. The full 3D scheme includes a 3D solution for the porous matrix, which is coupled with a 2D solution for fractures and a 1D solution for fracture intersections. As the hydraulic conductivity of the fracture intersections is usually higher than the hydraulic conductivity of the fractures and by at least one order of magnitude higher than the hydraulic conductivity of the porous matrix, the fastest flow and solute transport occurs in the fracture intersections. Therefore it is important to have an accurate and stable 1D solution of the transient AD problems. This article presents two different 1D BEM formulations for solution of the AD problems. The particular advantage of these formulations is that they provide one of the most straightforward and simplest ways to couple multiple intersecting 2D Boundary Element problems discretized with linear discontinuous elements. Both formulations are tested and compared for accuracy, stability, and consistency. The analysis helps to select the more suitable formulations according to the properties of the problem under consideration. © 2004 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2004     

The collocation method for the solution of boundary integral equations

This article is devoted to investigating the approximate solutions to a class of boundary integral equations over a closed, bounded and smooth surface found via the collocation method. This article provides sufficient conditions for the convergence of this method in the space of continuous functions.     

Note on exact solutions for a class of nonlinear diffusion equations

An effective characterization is given for a class of generalized nonlinear diffusion equations with power law dependent terms. Further, a new auxiliary equation ansatz is derived. Consequently, new exact traveling wave trigonometric function, solitary-like and Weierstrass elliptic solutions to a subclass are obtained by means of an auxiliary equation method and a generalized Riccati equation expansion method.     

A second order finite difference-spectral method for space fractional diffusion equations

A high order finite difference-spectral method is derived for solving space fractional diffusion equations,by combining the second order finite difference method in time and the spectral Galerkin method in space.The stability and error estimates of the temporal semidiscrete scheme are rigorously discussed,and the convergence order of the proposed method is proved to be O(τ2+Nα-m)in L2-norm,whereτ,N,αand m are the time step size,polynomial degree,fractional derivative index and regularity of the exact solution,respectively.Numerical experiments are carried out to demonstrate the theoretical analysis.     

Amplitude equation for the stochastic reaction‐diffusion equations with random Neumann boundary conditions

In this paper, we consider a quite general class of reaction‐diffusion equations with cubic nonlinearities and with random Neumann boundary conditions. We derive rigorously amplitude equations, using the natural separation of time‐scales near a change of stability and investigate whether additive degenerate noise and random boundary conditions can lead to stabilization of the solution of the stochastic partial differential equation or not. The nonlinear heat equation (Ginzburg–Landau equation) is used to illustrate our result. Copyright © 2015 John Wiley & Sons, Ltd.     

反应扩散方程的奇摄动

§ 1　 IntroductionThe nonlinearsingularly perturbed problem is a very attractive subjectof study in theinternational academic circles[1 ] .During the past decade many approximate methods havebeen developed and refined,including the method of average,boundary layer method,matched asymptotic expansions,and multiple scales.Recently,many scholars,for example,Bohé[2 ] ,Butuzov and Smurov[3] ,O Malley[4] ,Butuzov,Nefedov and Schneider[5] ,Kelley[6]and so on did a great deal of work.Mo consider…     

Combined boundary integral equations for acoustic scattering‐resonance problems

In this paper, boundary integral formulations for a time‐harmonic acoustic scattering‐resonance problem are analyzed. The eigenvalues of eigenvalue problems resulting from boundary integral formulations for scattering‐resonance problems split in general into two parts. One part consists of scattering‐resonances, and the other one corresponds to eigenvalues of some Laplacian eigenvalue problem for the interior of the scatterer. The proposed combined boundary integral formulations enable a better separation of the unwanted spectrum from the scattering‐resonances, which allows in practical computations a reliable and simple identification of the scattering‐resonances in particular for non‐convex domains. The convergence of conforming Galerkin boundary element approximations for the combined boundary integral formulations of the resonance problem is shown in canonical trace spaces. Numerical experiments confirm the theoretical results. Copyright © 2016 John Wiley & Sons, Ltd.     

A nonlocal nonlinear diffusion equation with blowing up boundary conditions

We deal with boundary value problems (prescribing Dirichlet or Neumann boundary conditions) for a nonlocal nonlinear diffusion operator which is analogous to the porous medium equation. First, we prove existence, uniqueness and the validity of a comparison principle for these problems. Next, we impose boundary data that blow up in finite time and study the behavior of the solutions.