Homotopy Analysis Solution to Radial Diffusivity Equation of Slightly Compressible Fluid

The salient significance of the solution of radial diffusivity equation to well testing analysis done in oil and gas industry cannot be over-emphasized. Varieties of solutions have been proposed to the radial diffusivity equation, of which the Van Everdingen-Hurst constant terminal rate solution is the most widely accepted and the others are approximate solution having their respective limitations. The main objective of this project, being its first application to oil and gas industry, is to use a new mathematical technique, the homotopy analysis method (HAM) to solve the radial diffusivity equation for slightly compressible fluid. In Using HAM, the Boltzmann transformation method was used to transform the radial PDE to ODE, then a homotopy series was then constructed for the new equation with the linear boundary condition from the original radial diffusivity equation of slightly compressible fluid and the final equation then solved using computation software Maple. The result gotten reveals that the homotopy analysis method gives good results compared to the Van Everdingen and Hurst Solution (Exact solution) and thus proves to be very effective, simple, and accurate when compared to other form of solutions. Hence from the results gotten, Homotopy Analysis Method can therefore be applied in solving other non-linear equations in the petroleum engineering field since it is simple and accurate.

Diffusion Equations of the Electric Charges and Magnetic Flux

Innovative definitions of the electric and magnetic diffusivities through conducting mediums and innovative diffusion equations of the electric charges and magnetic flux are verified in this article. Such innovations depend on the analogy of the governing laws of diffusion of the thermal, electrical, and magnetic energies and newly defined natures of the electric charges and magnetic flux as energy, or as electromagnetic waves, that have electric and magnetic potentials. The introduced diffusion equations of the electric charges and magnetic flux involve Laplacian operator and the introduced diffusivities. Both equations are applied to determine the electric and magnetic fields in conductors as the heat diffusion equation which is applied to determine the thermal field in steady and unsteady heat diffusion conditions. The use of electric networks for experimental modeling of thermal networks represents sufficient proof of similarity of the diffusion equations of both fields. By analysis of the diffusion phenomena of the three considered modes of energy transfer;the rates of flow of these energies are found to be directly proportional to the gradient of their volumetric concentration, or density, and the proportionality constants in such relations are the diffusivity of each energy. Such analysis leads also to find proportionality relations between the potentials of such energies and their volumetric concentrations. Validity of the introduced diffusion equations is verified by correspondence their solutions to the measurement results of the electric and magnetic fields in microwave ovens.

ANALYSIS AND DISCRETIZATION FOR AN OPTIMAL CONTROL PROBLEM OF A VARIABLE-COEFFICIENT RIESZ-FRACTIONAL DIFFUSION EQUATION WITH POINTWISE CONTROL CONSTRAINTS

We present a mathematical and numerical study for a pointwise optimal control problem governed by a variable-coefficient Riesz-fractional diffusion equation.Due to the impact of the variable diffusivity coefficient,existing regularity results for their constantcoefficient counterparts do not apply,while the bilinear forms of the state(adjoint)equation may lose the coercivity that is critical in error estimates of the finite element method.We reformulate the state equation as an equivalent constant-coefficient fractional diffusion equation with the addition of a variable-coefficient low-order fractional advection term.First order optimality conditions are accordingly derived and the smoothing properties of the solutions are analyzed by,e.g.,interpolation estimates.The weak coercivity of the resulting bilinear forms are proven via the Garding inequality,based on which we prove the optimal-order convergence estimates of the finite element method for the(adjoint)state variable and the control variable.Numerical experiments substantiate the theoretical predictions.

A LOCAL DISCONTINUOUS GALERKIN METHOD FOR TIME-FRACTIONAL DIFFUSION EQUATIONS

In this paper,a local discontinuous Galerkin(LDG)scheme for the time-fractional diffusion equation is proposed and analyzed.The Caputo time-fractional derivative(of orderα,with 0<α<1)is approximated by a finite difference method with an accuracy of order3-α,and the space discretization is based on the LDG method.For the finite difference method,we summarize and supplement some previous work by others,and apply it to the analysis of the convergence and stability of the proposed scheme.The optimal error estimate is obtained in the L2norm,indicating that the scheme has temporal(3-α)th-order accuracy and spatial(k+1)th-order accuracy,where k denotes the highest degree of a piecewise polynomial in discontinuous finite element space.The numerical results are also provided to verify the accuracy and efficiency of the considered scheme.

Finite Difference Schemes for Time-Space Fractional Diffusion Equations in One-and Two-Dimensions

In this paper,finite difference schemes for solving time-space fractional diffusion equations in one dimension and two dimensions are proposed.The temporal derivative is in the Caputo-Hadamard sense for both cases.The spatial derivative for the one-dimensional equation is of Riesz definition and the two-dimensional spatial derivative is given by the fractional Laplacian.The schemes are proved to be unconditionally stable and convergent.The numerical results are in line with the theoretical analysis.

Exponential Time Differencing Method for a Reaction-Diffusion System with Free Boundary

For reaction-diffusion equations in irregular domains with moving boundaries,the numerical stability constraints from the reaction and diffusion terms often require very restricted time step sizes,while complex geometries may lead to difficulties in the accuracy when discretizing the high-order derivatives on grid points near the boundary.It is very challenging to design numerical methods that can efficiently and accurately handle both difficulties.Applying an implicit scheme may be able to remove the stability constraints on the time step,however,it usually requires solving a large global system of nonlinear equations for each time step,and the computational cost could be significant.Integration factor(IF)or exponential time differencing(ETD)methods are one of the popular methods for temporal partial differential equations(PDEs)among many other methods.In our paper,we couple ETD methods with an embedded boundary method to solve a system of reaction-diffusion equations with complex geometries.In particular,we rewrite all ETD schemes into a linear combination of specificФ-functions and apply one state-of-the-art algorithm to compute the matrix-vector multiplications,which offers significant computational advantages with adaptive Krylov subspaces.In addition,we extend this method by incorporating the level set method to solve the free boundary problem.The accuracy,stability,and efficiency of the developed method are demonstrated by numerical examples.

带分数阶Robin边界条件的时间-空间分数阶扩散方程的有限差分方法

In this paper, an efficient numerical method is proposed to solve the Caputo-Riesz fractional diffusion equation with fractional Robin boundary conditions. We approximate the Riesz space fractional derivatives using the fractional central difference scheme with second-order accurate. A priori estimation of the solution of the numerical scheme is given, and the stability and convergence of the numerical scheme are analyzed.Finally, a numerical example is used to verify the accuracy and efficiency of the numerical method.

A CLASS OF NONLINEAR SINGULARLY PERTURBED PROBLEMS FOR REACTION DIFFUSION EQUATIONS

A class of nonlinear singularly perturbed problems for reaction diffusion equations are considered. Under suitable conditions, by using the theory of differential inequalities, the asymptotic behavior of solutions for the initial boundary value problems are studied, reduced problems of which possess two intersecting solutions.

THE LARGE TIME BEHAVIOR OF SPECTRAL APPROXIMATION FOR A CLASS OF PSEUDOPARABOLIC VISCOUS DIFFUSION EQUATION

The asymptotic behavior of the solutions to a class of pseudoparabolic viscous diffusion equation with periodic initial condition is studied by using the spectral method. The semidiscrete Fourier approximate solution of the problem is constructed and the error estimation between spectral approximate solution and exact solution on large time is also obtained. The existence of the approximate attractor AN and the upper semicontinuity d（AN,A） → 0 are proved.

A meshless method based on moving Kriging interpolation for a two-dimensional time-fractional diffusion equation

Fractional diffusion equations have been the focus of modeling problems in hydrology, biology, viscoelasticity, physics, engineering, and other areas of applications. In this paper, a meshfree method based on the moving Kriging inter- polation is developed for a two-dimensional time-fractional diffusion equation. The shape function and its derivatives are obtained by the moving Kriging interpolation technique. For possessing the Kronecker delta property, this technique is very efficient in imposing the essential boundary conditions. The governing time-fractional diffusion equations are transformed into a standard weak formulation by the Galerkin method. It is then discretized into a meshfree system of time-dependent equations, which are solved by the standard central difference method. Numerical examples illustrating the applicability and effectiveness of the proposed method are presented and discussed in detail.

THE QUASI-BOUNDARY VALUE METHOD FOR IDENTIFYING THE INITIAL VALUE OF THE SPACE-TIME FRACTIONAL DIFFUSION EQUATION

In this article,we consider to solve the inverse initial value problem for an inhomogeneous space-time fractional diffusion equation.This problem is ill-posed and the quasi-boundary value method is proposed to deal with this inverse problem and obtain the series expression of the regularized solution for the inverse initial value problem.We prove the error estimates between the regularization solution and the exact solution by using an a priori regularization parameter and an a posteriori regularization parameter choice rule.Some numerical results in one-dimensional case and two-dimensional case show that our method is efficient and stable.

IDENTIFYING AN UNKNOWN SOURCE IN SPACE-FRACTIONAL DIFFUSION EQUATION

In this paper, we identify a space-dependent source for a fractional diffusion equation. This problem is ill-posed, i.e., the solution （if it exists） does not depend continuously on the data. The generalized Tikhonov regularization method is proposed to solve this problem. An a priori error estimate between the exact solution and its regularized approximation is obtained. Moreover, an a posteriori parameter choice rule is proposed and a stable error estimate is also obtained, Numerical examples are presented to illustrate the validity and effectiveness of this method.

SINGULAR SOLUTIONS FOR A CONVECTION DIFFUSION EQUATION WITH ABSORPTION

In this paper we discuss the existence and nonexistence of singular solutions for a porous medium equations with convection and absorption terms.

Superconvergence Study of the Direct Discontinuous Galerkin Method and Its Variations for Diffusion Equations

In this paper,we apply the Fourier analysis technique to investigate superconvergence properties of the direct disontinuous Galerkin(DDG)method(Liu and Yan in SIAM J Numer Anal 47(1):475-698,2009),the DDG method with the interface correction(DDGIC)(Liu and Yan in Commun Comput Phys 8(3):541-564,2010),the symmetric DDG method(Vidden and Yan in Comput Math 31(6):638-662,2013),and the nonsymmetric DDG method(Yan in J Sci Comput 54(2):663-683,2013).We also include the study of the interior penalty DG(IPDG)method,due to its close relation to DDG methods.Error estimates are carried out for both P2 and P3 polynomial approximations.By investigating the quantitative errors at the Lobatto points,we show that the DDGIC and symmetric DDG methods are superior,in the sense of obtaining(k+2)th superconvergence orders for both P2 and P3 approximations.Superconvergence order of(k+2)is also observed for the IPDG method with P3 polynomial approximations.The errors are sensitive to the choice of the numerical flux coefficient for even degree P2 approximations,but are not for odd degree P3 approxi-mations.Numerical experiments are carried out at the same time and the numerical errors match well with the analytically estimated errors.

SOME RESULTS ON A DEGENERATE AND SINGULAR DIFFUSION EQUATION

In this article, a degenerate and singular diffusion problem is studied. The existence of solutions is established by parabolic regularization. Some properties of solutions, for instance, asymptotic behavior, are also discussed.

Maximum modulus principle estimates for one dimensional fractional diffusion equation

We present scheme I for solving one-dimensional fractional diffusion equation with variable coefficients based on the maximum modulus principle and two Grunwald approxima- tions. Scheme II is obtained by using classic Crank-Nicolson approximations in order to improve the time convergence. Schemes are proved to be unconditionally stable and second-order accuracy in spatial grid size for the problem with order of fractional derivative belonging to [（√17- 1）/2, 2] using the maximum modulus principle. A numerical example is given to confirm the theoretical analysis result.

New variable separation solutions for the generalized nonlinear diffusion equations

The functionally generalized variable separation of the generalized nonlinear diffusion equations ut = A（u, Ux）Uxx ＋ B（u, ux） is studied by using the conditional Lie-Blicklund symmetry method. The variant forms of the considered equations, which admit the corresponding conditional Lie--Biicklund symmetries, are characterized. To construct functionally gener- alized separable solutions, several concrete examples defined on the exponential and trigonometric invariant subspaces are provided.

THREE HIGH-ORDER SPLITTING SCHEMES FOR 3D TRANSPORT EQUATION

Two high-order splitting schemes based on the idea of the operators splitting method are given. The three-dimensional advection-diffusion equation was split into several one-dimensional equations that were solved by these two schemes, only three computational grid points were needed in each direction but the accuracy reaches the spatial fourth-order. The third scheme proposed is based on the classical ADI scheme and the accuracy of the advection term of it can reach the spatial fourth-order. Finally, two typical numerical experiments show that the solutions of these three schemes compare well with that given by the analytical solution when the Peclet number is not bigger than 5.

Mixed time discontinuous space-time finite element method for convection diffusion equations

A mixed time discontinuous space-time finite element scheme for secondorder convection diffusion problems is constructed and analyzed. Order of the equation is lowered by the mixed finite element method. The low order equation is discretized with a space-time finite element method, continuous in space but discontinuous in time. Stability, existence, uniqueness and convergence of the approximate solutions are proved. Numerical results are presented to illustrate efficiency of the proposed method.

ASYMPTOTIC PROPERTIES OF FUNCTIONAL DIFFERENTIAL EQUATIONS IN BANACH SPACES

The paper is devoted to the asymptotic properties of functional differential equations in Banach spaces.The criteria of the invariant and attracting sets are obtained.Particularly, the sufficient condition of asymptotic stability of the equilibrium point is given as the system has an equilibrium point.Several examples are also worked out to demonstrate the validity of the results.