Numerical Simulation of Modified Kortweg-de Vries Equation by Linearized Implicit Schemes

In this paper, we are going to derive four numerical methods for solving the Modified Kortweg-de Vries (MKdV) equation using fourth Pade approximation for space direction and Crank Nicolson in the time direction. Two nonlinear schemes and two linearized schemes are presented. All resulting schemes will be analyzed for accuracy and stability. The exact solution and the conserved quantities are used to highlight the efficiency and the robustness of the proposed schemes. Interaction of two and three solitons will be also conducted. The numerical results show that the interaction behavior is elastic and the conserved quantities are conserved exactly, and this is a good indication of the reliability of the schemes which we derived. A comparison with some existing is presented as well.

Quinpi:Integrating Conservation Laws with CWENO Implicit Methods

Many interesting applications of hyperbolic systems of equations are stiff,and require the time step to satisfy restrictive stability conditions.One way to avoid small time steps is to use implicit time integration.Implicit integration is quite straightforward for first-order schemes.High order schemes instead also need to control spurious oscillations,which requires limiting in space and time also in the linear case.We propose a framework to simplify considerably the application of high order non-oscillatory schemes through the introduction of a low order implicit predictor,which is used both to set up the nonlinear weights of a standard high order space reconstruction,and to achieve limiting in time.In this preliminary work,we concentrate on the case of a third-order scheme,based on diagonally implicit Runge Kutta(DIRK)integration in time and central weighted essentially non-oscillatory(CWENO)reconstruction in space.The numerical tests involve linear and nonlinear scalar conservation laws.

A three dimensional implicit immersed boundary method with application

Most algorithms of the immersed boundary method originated by Peskin are explicit when it comes to the computation of the elastic forces exerted by the immersed boundary to the fluid. A drawback of such an explicit approach is a severe restriction on the time step size for maintaining numerical stability. An implicit immersed boundary method in two dimensions using the lattice Boltzmann approach has been proposed. This paper reports an extension of the method to three dimensions and its application to simulation of a massive flexible sheet interacting with an incompressible viscous flow.

An implicit method using contravariant velocity components and its application to calculations in a harbour-channel area

The key problem in the computation of fluid dynamics using fine boundary-fitted grids is how to improve the numerical stability and decrease the calculating quantity. To solve this problem, implicit schemes should be adopted since explicit schemes may bring about a great increase in computation quantity according to the Courant-FrledrichsLewy condition. Whereas the adoption of implicit schemes is difficult to be realized because of the existence of two partial derivatives of surface elevations with respect to variables of alternative direction coordinates in each momentum equation in non-rectangular coordinates. With an aim to design an implicit scheme in non-reetangular ccordinates in the present paper, new momentum equations with the contravariant components of velocity vector are derived based on the shallow water dynamic equations in generalized curvilinear coordinates. In each equation, the coefficients before the two detivatives of surface elevations have different orders of magnitude, i. e., the derivative with the larger ceefficient rnay play a more important role than that with the smaller one. With this advantage, the ADI scheme can then be easily employed to improve the numerical stability and decrease the calculating quantity. The calculation in a harbour and a channel in Macau nearshore area shows that the implicit model is effective in calculating current fields in small size areas.

Analysis of an Implicit Finite Difference Scheme for Time Fractional Diffusion Equation

Time fractional diffusion equation is usually used to describe the problems involving non-Markovian random walks. This kind of equation is obtained from the standard diffusion equation by replacing the first-order time derivative with a fractional derivative of order α∈(0, 1). In this paper, an implicit finite difference scheme for solving the time fractional diffusion equation with source term is presented and analyzed, where the fractional derivative is described in the Caputo sense. Stability and convergence of this scheme are rigorously established by a Fourier analysis. And using numerical experiments illustrates the accuracy and effectiveness of the scheme mentioned in this paper.

Strong Convergence of an Implicit Iteration Process for a Finite Family of Asymptotically Ф-pseudocontractive Mappings

Strong convergence theorems for approximation of common fixed points of asymptotically Ф-quasi-pseudocontractive mappings and asymptotically C-strictly- pseudocontractive mappings are proved in real Banach spaces by using a new composite implicit iteration scheme with errors. The results presented in this paper extend and improve the main results of Sun, Gu and Osilike published on J. Math. Anal. Appl.

TRANSONIC FLOW CALCULATION OF EULER EQUATIONS BY IMPLICIT ITERATING SCHEME WITH FLUX SPLITTING

Three dimensional Euler equations are solved in the finite volume form with van Leer's flux vector splitting technique. Block matrix is inverted by Gauss-Seidel iteration in two dimensional plane while strongly implicit alternating sweeping is implemented in the direction of the third dimension. Very rapid convergence rate is obtained with CFL number reaching the order of 100. The memory resources can be greatly saved too. It is verified that the reflection boundary condition can not be used with flux vector splitting since it will produce too large numerical dissipation. The computed flow fields agree well with experimental results. Only one or two grid points are there within the shock transition zone.

Efficient Difference Schemes for the Caputo-Tempered Fractional Diffusion Equations Based on Polynomial Interpolation

The tempered fractional calculus has been successfully applied for depicting the time evolution of a system describing non-Markovian diffusion particles.The related governing equations are a series of partial differential equations with tempered fractional derivatives.Using the polynomial interpolation technique,in this paper,we present three efficient numerical formulas,namely the tempered L1 formula,the tempered L1-2 formula,and the tempered L2-1_(σ)formula,to approximate the Caputo-tempered fractional derivative of orderα∈(0,1).The truncation error of the tempered L1 formula is of order 2-α,and the tempered L1-2 formula and L2-1_(σ)formula are of order 3-α.As an application,we construct implicit schemes and implicit ADI schemes for one-dimensional and two-dimensional time-tempered fractional diffusion equations,respectively.Furthermore,the unconditional stability and convergence of two developed difference schemes with tempered L1 and L2-1_(σ)formulas are proved by the Fourier analysis method.Finally,we provide several numerical examples to demonstrate the correctness and effectiveness of the theoretical analysis.

Efficient Finite Difference Methods for the Numerical Analysis of One-Dimensional Heat Equation

In this paper, we investigate and analyze one-dimensional heat equation with appropriate initial and boundary condition using finite difference method. Finite difference method is a well-known numerical technique for obtaining the approximate solutions of an initial boundary value problem. We develop Forward Time Centered Space (FTCS) and Crank-Nicolson (CN) finite difference schemes for one-dimensional heat equation using the Taylor series. Later, we use these schemes to solve our governing equation. The stability criterion is discussed, and the stability conditions for both schemes are verified. We exhibit the results and then compare the results between the exact and approximate solutions. Finally, we estimate error between the exact and approximate solutions for a specific numerical problem to present the convergence of the numerical schemes, and demonstrate the resulting error in graphical representation.

Efficient Simulation of Wave Propagation with Implicit Finite Difference Schemes

Finite difference method is an important methodology in the approximation of waves.In this paper,we will study two implicit finite difference schemes for the simulation of waves.They are the weighted alternating direction implicit(ADI)scheme and the locally one-dimensional(LOD)scheme.The approximation errors,stability conditions,and dispersion relations for both schemes are investigated.Our analysis shows that the LOD implicit scheme has less dispersion error than that of the ADI scheme.Moreover,the unconditional stability for both schemes with arbitrary spatial accuracy is established for the first time.In order to improve computational efficiency,numerical algorithms based on message passing interface(MPI)are implemented.Numerical examples of wave propagation in a three-layer model and a standard complex model are presented.Our analysis and comparisons show that both ADI and LOD schemes are able to efficiently and accurately simulate wave propagation in complex media.

ON THE CONVERGENCE OF IMPLICIT DIFFERENCE SCHEMES FOR HYPERBOLIC CONSERVATION LAWS

Discusses the treatment of implicit difference approximations to hyperbolic conservation laws with non-convex flux. Information on weighted implicit difference schemes; Convergence of difference scheme; Establishment of the convergence of the approximate solution toward the entropy solution.

NUMERICAL SIMULATION ANALYSIS FOR MIGRATION-ACCUMULATION OF OIL AND WATER

Numerical simulation of oil migration and accumulation is to describe the history of oil migration and accumulation in basin evolution. It is of great value in the exploration oil resources and their rational evaluation. This paper, puts forward the mathematical model and the modified method of alternating direction implicit interactive scheme. For the famous hydraulic experiment of secondary migration_accumulation (cut plane and plane problem), it has been done the numerical simulation test, and both the computational and experimental results are identical.

High-order implicit discontinuous Galerkin schemes for unsteady compressible Navier–Stokes equations

Efficient solution techniques for high-order temporal and spatial discontinuous Galerkin（DG） discretizations of the unsteady Navier–Stokes equations are developed. A fourth-order implicit Runge–Kutta（IRK） scheme is applied for the time integration and a multigrid preconditioned GMRES solver is extended to solve the nonlinear system arising from each IRK stage. Several modifications to the implicit solver have been considered to achieve the efficiency enhancement and meantime to reduce the memory requirement. A variety of time-accurate viscous flow simulations are performed to assess the resulting high-order implicit DG methods. The designed order of accuracy for temporal discretization scheme is validate and the present implicit solver shows the superior performance by allowing quite large time step to be used in solving time-implicit systems. Numerical results are in good agreement with the published data and demonstrate the potential advantages of the high-order scheme in gaining both the high accuracy and the high efficiency.

High-order numerical methods of fractional-order Stokes' first problem for heated generalized second grade fluid

The high-order implicit finite difference schemes for solving the fractional- order Stokes＇ first problem for a heated generalized second grade fluid with the Dirichlet boundary condition and the initial condition are given. The stability, solvability, and convergence of the numerical scheme are discussed via the Fourier analysis and the matrix analysis methods. An improved implicit scheme is also obtained. Finally, two numerical examples are given to demonstrate the effectiveness of the mentioned schemes

THE HIGH ACCURACY EXPLICIT DIFFERENCE SCHEME FOR SOLVING PARABOLIC EQUATIONS 3-DIMENSION

In this paper, an explicit three_level symmetrical differencing scheme with parameters for solving parabolic partial differential equation of three_dimension will be considered. The stability condition and local truncation error for the scheme are r<1/2 and O( Δ t 2+ Δ x 4+ Δ y 4+ Δ z 4) ,respectively.

Efficient high-order immersed interface methods for heat equations with interfaces

An efficient high-order immersed interface method （IIM） is proposed to solve two-dimensional （2D） heat problems with fixed interfaces on Cartesian grids, which has the fourth-order accuracy in the maximum norm in both time and space directions. The space variable is discretized by a high-order compact （HOC） difference scheme with correction terms added at the irregular points. The time derivative is integrated by a Crank-Nicolson and alternative direction implicit （ADI） scheme. In this case, the time accuracy is just second-order. The Richardson extrapolation method is used to improve the time accuracy to fourth-order. The numerical results confirm the convergence order and the efficiency of the method.

Numerical Methods for Solving Logarithmic Nonlinear Schrödinger’s Equation

In this study, we will construct numerical techniques for tackling the logarithmic Schr&#246;dinger’s nonlinear equation utilizing the explicit scheme and the Crank-Nicolson scheme of the finite difference method. These schemes will be subjected to accuracy and stability tests before being used. Efficacy and robustness of the techniques under consideration will be demonstrated using an exact solution, one-Gausson, as well as conserved quantities. Interaction of two-soliton will be conducted. The numerical findings revealed, the interplay behavior is flexible.

Effects of the Nature of Boundaries Conditions and Their Truncation Errors on the Distribution of Minority Carriers in Silicon Solar Cell

In this work, the effects of boundaries conditions and truncation errors in the distribution of minority carriers in the semiconductor are studied. It is a one-dimensional digital study of a polycrystalline silicon solar cell under polychromatic illumination in a dynamic state. Starting from the Boltzmann equation of semiconductors, the author establishes the general equation of particle transport. Assumptions made on the latter allow it to give the equation of distribution of minority carriers in a general way in its case to be studied. This dimensioned distribution equation reveals the parameters of influences on the distribution of carriers. It obtains a partial derivative equation for the carrier distribution function. The boundary conditions are then discretized to order one and then to order two. By considering boundary conditions and the nature of the carriers, the author numerically resolves the discretized general equation by assessing the influence of the nature of the boundary conditions and truncation errors and the influence of the discretization step on the density of the charge carriers by setting certain parameters and varying others. The work ends with a conclusion and logical follow-up to this work.

Analysis of Gray Scott’s Model Numerically

In this paper, a two-dimensional nonlinear coupled Gray Scott system is simulated with a finite difference scheme and a finite volume technique. Pre and post-processing lead to a new solution called GSmFoam by understanding geometry settings and mesh information. The concentration profile changes over time, as does the intensity of the contour patterns. The OpenFoam solver gives you the confidence to compare the pattern result with efficient numerical algorithms on the Gray Scott model.

Alternating Cell Direction Implicit Method using Approximate Factorization on Hybrid Grids

In this study,a novel fast-implicit iteration scheme called the alternating cell direction implicit(ACDI)method is combined with the approximate factorization scheme.This application aims to offer a mathematically well-defined version of the ACDI method and to increase the accuracy of the iteration scheme used for the nu-merical solutions of partial differential equations.The ACDI method is a fast-implicit method that can be used for unstructured grids.The use of fast implicit iteration meth-ods with unstructured grids is not common in the literature.The new ACDI method has been applied to the unsteady diffusion equation to determine its convergence and time-dependent solution ability and character.The numerical tests are conducted for different grid types,such as structured,unstructured quadrilateral,and hybrid polygonal grids.Second,the ACDI was applied to the unsteady advection-diffusion equation to understand the time-dependent and progression capabilities of the pre-sented method.Third,a full potential equation solution is created to understand the complexflow solving ability of the presented method.The results of the numerical study are compared with other fast implicit methods,such as the point Gauss–Seidel(PGS)and line Gauss–Seidel(LGS)methods and the fourth-order Runge-Kutta(RK4)method,which is an explicit scheme,and the Laasonen method,which is a fully im-plicit scheme.The study increased the abilities of the ACDI method.Due to the new ACDI method,the approximate factorization method,which is used only in structural grids that are known to be advantageous,can be applied to any mesh structure.