The Sinc‐Galerkin method for solving an inverse parabolic problem with unknown source term

The inverse problem of determining an unknown source term depending on space variable in a parabolic equation is considered. A numerical algorithm is presented for recovering the unknown function and obtaining a solution of the problem. As this inverse problem is ill‐posed, Tikhonov regularization is used for finding a stable solution. For solving the direct problem, a Galerkin method with the Sinc basis functions in both the space and time domains is presented. This approximate solution displays an exponential convergence rate and is valid on the infinite time interval. Finally, some examples are presented to illustrate the ability and efficiency of this numerical method. © 2012 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2013     

Analytical and numerical solutions of a two‐dimensional multi‐term time‐fractional Oldroyd‐B model

In this paper, we consider a two‐dimensional multi‐term time‐fractional Oldroyd‐B equation on a rectangular domain. Its analytical solution is obtained by the method of separation of variables. We employ the finite difference method with a discretization of the Caputo time‐fractional derivative to obtain an implicit difference approximation for the equation. Stability and convergence of the approximation scheme are established in the L_∞ ‐norm. Two examples are given to illustrate the theoretical analysis and analytical solution. The results indicate that the present numerical method is effective for this general two‐dimensional multi‐term time‐fractional Oldroyd‐B model.     

D'Alembert‐type solution of the Cauchy problem for the Boussinesq‐Klein‐Gordon equation

In this paper, we construct a weakly‐nonlinear d'Alembert‐type solution of the Cauchy problem for the Boussinesq‐Klein‐Gordon (BKG) equation. Similarly to our earlier work based on the use of spatial Fourier series, we consider the problem in the class of periodic functions on an interval of finite length (including the case of localized solutions on a large interval), and work with the nonlinear partial differential equation with variable coefficients describing the deviation from the oscillating mean value. Unlike our earlier paper, here we develop a novel multiple‐scales procedure involving fast characteristic variables and two slow time scales and averaging with respect to the spatial variable at a constant value of one or another characteristic variable, which allows us to construct an explicit and compact d'Alembert‐type solution of the nonlinear problem in terms of solutions of two Ostrovsky equations emerging at the leading order and describing the right‐ and left‐propagating waves. Validity of the constructed solution in the case when only the first initial condition for the BKG equation may have nonzero mean value follows from our earlier results, and is illustrated numerically for a number of instructive examples, both for periodic solutions on a finite interval, and localized solutions on a large interval. We also outline an extension of the procedure to the general case, when both initial conditions may have nonzero mean values. Importantly, in all cases, the initial conditions for the leading‐order Ostrovsky equations by construction have zero mean, while initial conditions for the BKG equation may have nonzero mean values.     

Construction of global‐in‐time solutions to Kolmogorov‐Feller pseudodifferential equations with a small parameter using characteristics

Using an idea going back to Madelung, we construct global in time solutions to the transport equation corresponding to the asymptotic solution of the Kolmogorov‐Feller equation describing a system with diffusion, potential and jump terms. To do that we use the construction of a generalized delta‐shock solution of the continuity equation for a discontinuous velocity field. We also discuss corresponding problem of asymptotic solution construction (Maslov tunnel asymptotics).     

Diffusion in poro‐plastic media

A model is developed for the flow of a slightly compressible fluid through a saturated inelastic porous medium. The initial‐boundary‐value problem is a system that consists of the diffusion equation for the fluid coupled to the momentum equation for the porous solid together with a constitutive law which includes a possibly hysteretic relation of elasto‐visco‐plastic type. The variational form of this problem in Hilbert space is a non‐linear evolution equation for which the existence and uniqueness of a global strong solution is proved by means of monotonicity methods. Various degenerate situations are permitted, such as incompressible fluid, negligible porosity, or a quasi‐static momentum equation. The essential sufficient conditions for the well‐posedness of the system consist of an ellipticity condition on the term for diffusion of fluid and either a viscous or a hardening assumption in the constitutive relation for the porous solid. Copyright © 2004 John Wiley & Sons, Ltd.     

Existence and zero‐electron‐mass limit of strong solutions to the stationary compressible Navier‐Stokes‐Poisson equation with large external force

In this study, we consider a viscous compressible model of plasma and semiconductors, which is expressed as a compressible Navier‐Stokes‐Poisson equation. We prove that there exists a strong solution to the boundary value problem of the steady compressible Navier‐Stokes‐Poisson equation with large external forces in bounded domain, provided that the ratio of the electron/ions mass is appropriately small. Moreover, the zero‐electron‐mass limit of the strong solutions is rigorously verified. The main idea in the proof is to split the original equation into 4 parts, a system of stationary incompressible Navier‐Stokes equations with large forces, a system of stationary compressible Navier‐Stokes equations with small forces, coupled with 2 Poisson equations. Based on the known results about linear incompressible Navier‐Stokes equation, linear compressible Navier‐Stokes, linear transport, and Poisson equations, we try to establish uniform in the ratio of the electron/ions mass a priori estimates. Further, using Schauder fixed point theorem, we can show the existence of a strong solution to the boundary value problem of the steady compressible Navier‐Stokes‐Poisson equation with large external forces. At the same time, from the uniform a priori estimates, we present the zero‐electron‐mass limit of the strong solutions, which converge to the solutions of the corresponding incompressible Navier‐Stokes‐Poisson equations.     

A non‐standard finite difference scheme for an advection‐diffusion‐reaction equation

In this note, a non‐standard finite difference (NSFD) scheme is proposed for an advection‐diffusion‐reaction equation with nonlinear reaction term. We first study the diffusion‐free case of this equation, that is, an advection‐reaction equation. Two exact finite difference schemes are constructed for the advection‐reaction equation by the method of characteristics. As these exact schemes are complicated and are not convenient to use, an NSFD scheme is derived from the exact scheme. Then, the NSFD scheme for the advection‐reaction equation is combined with a finite difference space‐approximation of the diffusion term to provide a NSFD scheme for the advection‐diffusion‐reaction equation. This new scheme could preserve the fixed points, the positivity, and the boundedness of the solution of the original equation. Numerical experiments verify the validity of our analytical results. Copyright © 2014 JohnWiley & Sons, Ltd.     

Determination of a time‐ and space‐dependent heat flux relaxation function by means of a restricted Dirichlet‐to‐Neumann operator

We consider an inverse problem to recover a space‐ and time‐dependent relaxation function of heat flux in a three‐dimensional body on the basis of the restriction of the Dirichlet‐to‐Neumann operator of the related equation of heat flow onto a set of Dirichlet data of the form of a product of a fixed time‐dependent coefficient and a free space‐dependent function. Uniqueness of the solution of this inverse problem is proved. Copyright © 2004 John Wiley & Sons, Ltd.     

Determination of an unknown function in a parabolic inverse problem by Sinc‐collocation method

In this article, an inverse problem of determining an unknown time‐dependent source term of a parabolic equation is considered. We change the inverse problem to a Volterra integral equation of convolution‐type. By using Sinc‐collocation method, the resulting integral equation is replaced by a system of linear algebraic equations. The convergence analysis is included, and it is shown that the error in the approximate solution is bounded in the infinity norm by the condition number and the norm of the inverse of the coefficient matrix multiplied by a factor that decays exponentially with the size of the system. Some examples are given to demonstrate the computational efficiency of the method. © 2010 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 27: 1584–1598, 2010     

Stability of non‐constant steady‐state solutions for non‐isentropic Euler–Maxwell system with a temperature damping term

This work is concerned with the periodic problem for compressible non‐isentropic Euler–Maxwell systems with a temperature damping term arising in plasmas. For this problem, we prove the global in time existence of a smooth solution around a given non‐constant steady state with the help of an induction argument on the order of the mixed time‐space derivatives of solutions in energy estimates. Moreover, we also show the convergence of the solution to this steady state as the time goes to the infinity. This phenomenon on the charge transport shows the essential relation of the systems with the non‐isentropic Euler–Maxwell and the isentropic Euler–Maxwell systems. Copyright © 2015 John Wiley & Sons, Ltd.     

Similar constructing method for solving nonhomogeneous mixed boundary value problem of n‐interval composite ode and its application

This paper studies a simple method—Similar Constructing Method (SCM)—for constructing the exact solutions of the nonhomogeneous mixed boundary value problem for sets of n‐interval composite second‐order ordinary differential equation (ODE) with variable coefficient. Then this paper proves the correctness of the solution obtained by SCM. After that, this paper has done simulation experiment. This section uses the SCM to solve the nonhomogeneous boundary value problem of three‐interval composite Bessel equation. Solutions are presented in graphical form for various parameter values, and the influence of parameters on the solution is analyzed. The example shows that using SCM to solve the class of nonhomogeneous mixed boundary value problems of n‐interval composite second‐order linear ODE is easy, convenient, and effective.     

Inverse problem for a space‐time fractional diffusion equation: Application of fractional Sturm‐Liouville operator

An inverse problem of determining a time‐dependent source term from the total energy measurement of the system (the over‐specified condition) for a space‐time fractional diffusion equation is considered. The space‐time fractional diffusion equation is obtained from classical diffusion equation by replacing time derivative with fractional‐order time derivative and Sturm‐Liouville operator by fractional‐order Sturm‐Liouville operator. The existence and uniqueness results are proved by using eigenfunction expansion method. Several special cases are discussed, and particular examples are provided.     

Spatial behaviour of solutions of the dual‐phase‐lag heat equation

In this paper we study the spatial behaviour of solutions of some problems for the dual‐phase‐lag heat equation on a semi‐infinite cylinder. The theory of dual‐phase‐lag heat conduction leads to a hyperbolic partial differential equation with a third derivative with respect to time. First, we investigate the spatial evolution of solutions of an initial boundary‐value problem with zero boundary conditions on the lateral surface of the cylinder. Under a boundedness restriction on the initial data, an energy estimate is obtained. An upper bound for the amplitude term in this estimate in terms of the initial and boundary data is also established. For the case of zero initial conditions, a more explicit estimate is obtained which shows that solutions decay exponentially along certain spatial‐time lines. A class of non‐standard problems is also considered for which the temperature and its first two time derivatives at a fixed time T are assumed proportional to their initial values. Copyright © 2004 John Wiley & Sons, Ltd.     

Accuracy of multivalue methods during change of time‐step size: Adams‐Moulton

Multivalue methods are a class of time‐stepping procedures for numerical solution of differential equations that progress to a new time level using the approximate solution for the function of interest and its derivatives at a single time level. The methods differ from multistep procedures, which make use of solutions to the differential equation at multiple time levels to advance to the new time level. Multistep methods are difficult to employ when a change in time‐step is desired, because the standard formulas (e.g., Adams‐Moulton or Gear) must be modified to accommodate the change. Multivalue methods seem to possess the desirable feature that the time‐step may be changed arbitrarily as one proceeds, since the solution proceeds from a single time level. However, in practice, changes in the time‐step introduce lower order errors or alter the coefficient in the truncation error term. Here, the multivalue Adams‐Moulton method is presented based on a general interpolation procedure. Modifications required to retain the high‐order accuracy of these methods during a change in time‐step are developed. Additionally, a formula for the unknown initial derivatives is presented. Finally, two examples are provided to illustrate the potential merit of the modification to the standard multivalue methods. © 2000 John Wiley & Sons, Inc. Numer Methods Partials Differential Eq 16: 312–326, 2000     

The quasi‐reversibility method for a final value problem of the time‐fractional diffusion equation with inhomogeneous source

This paper is devoted to discuss a multidimensional backward heat conduction problem for time‐fractional diffusion equation with inhomogeneous source. This problem is ill‐posed. We use quasi‐reversibility regularization method to solve this inverse problem. Moreover, the convergence estimates between regularization solution and the exact solution are obtained under the a priori and the a posteriori choice rules. Finally, the numerical examples for one‐dimensional and two‐dimensional cases are presented to show that our method is feasible and effective.     

Stability estimates for the solution to the problem of determining the kernel of a viscoelastic equation

For an integrodifferential equation corresponding to a two-dimensional viscoelastic problem, we study the problem of defining the spatial part of the kernel involved in the integral term of the equation. The support of the sought function is assumed to belong to a compact domain Ω. As information for solving this inverse problem, the traces of the solution to the direct Cauchy problem and its normal derivative are given for some finite time interval on the boundary of Ω. An important feature in the statement of the problem is the fact that the solution of the direct problem corresponds to the zero initial data and a force impulse in time localized on a fixed straight line disjoint with Ω. The main result of the article consists in obtaining a Lipschitz estimate for the conditional stability of the solution to the inverse problem under consideration.     

The Cauchy problem of a fluid‐particle interaction model with external forces

In this paper, we consider the Cauchy problem of a fluid‐particle interaction model with external forces. We first construct the asymptotic profile of the system. The global existence and uniqueness theorem for the solution near the profile is given. Finally, optimal decay rate of the solution to the background profile is obtained by combining the decay rate analysis of a linearized equation with energy estimates for the nonlinear terms. The main method used in this paper is the energy method combining with the macro‐micro decomposition.     

Ritz–Galerkin method for solving a parabolic equation with non‐local and time‐dependent boundary conditions

The paper is devoted to the investigation of a parabolic partial differential equation with non‐local and time‐dependent boundary conditions arising from ductal carcinoma in situ model. Approximation solution of the present problem is implemented by the Ritz–Galerkin method, which is a first attempt at tackling parabolic equation with such non‐classical boundary conditions. In the process of dealing with the difficulty caused by integral term in non‐local boundary condition, we use a trick of introducing the transition function G(x,t) to convert non‐local boundary to another non‐classical boundary, which can be handled with the Ritz–Galerkin method. Illustrative examples are included to demonstrate the validity and applicability of the technique in this paper. Copyright © 2015 John Wiley & Sons, Ltd.     

A priori error estimates of an extrapolated space‐time discontinuous galerkin method for nonlinear convection‐diffusion problems

We deal with the numerical solution of a scalar nonstationary nonlinear convection‐diffusion equation. We employ a combination of the discontinuous Galerkin finite element (DGFE) method for the space as well as time discretization. The linear diffusive and penalty terms are treated implicitly whereas the nonlinear convective term is treated by a special higher order explicit extrapolation from the previous time step, which leads to the necessity to solve only a linear algebraic problem at each time step. We analyse this scheme and derive a priori asymptotic error estimates in the L^∞(L²) –norm and the L²(H¹) –seminorm with respect to the mesh size h and time step τ. Finally, we present an efficient solution strategy and numerical examples verifying the theoretical results. © 2010 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 27: 1456–1482, 2010     

Numerical treatment of two‐parameter singularly perturbed parabolic convection diffusion problems with non‐smooth data

In the present work, we consider a parabolic convection‐diffusion‐reaction problem where the diffusion and convection terms are multiplied by two small parameters, respectively. In addition, we assume that the convection coefficient and the source term of the partial differential equation have a jump discontinuity. The presence of perturbation parameters leads to the boundary and interior layers phenomena whose appropriate numerical approximation is the main goal of this paper. We have developed a uniform numerical method, which converges almost linearly in space and time on a piecewise uniform space adaptive Shishkin‐type mesh and uniform mesh in time. Error tables based on several examples show the convergence of the numerical solutions. In addition, several numerical simulations are presented to show the effectiveness of resolving layer behavior and their locations.