Diffusion-Wave Type Solutions to the Second-Order Evolutionary Equation with Power Nonlinearities

The paper deals with a nonlinear second-order one-dimensional evolutionary equation related to applications and describes various diffusion, filtration, convection, and other processes. The particular cases of this equation are the well-known porous medium equation and its generalizations. We construct solutions that describe perturbations propagating over a zero background with a finite velocity. Such effects are known to be atypical for parabolic equations and appear as a consequence of the degeneration of the equation at the points where the desired function vanishes. Previously, we have constructed it, but here the case of power nonlinearity is considered. It allows for conducting a more detailed analysis. We prove a new theorem for the existence of solutions of this type in the class of piecewise analytical functions, which generalizes and specifies the earlier statements. We find and study exact solutions having the diffusion wave type, the construction of which is reduced to the second-order Cauchy problem for an ordinary differential equation (ODE) that inherits singularities from the original formulation. Statements that ensure the existence of global continuously differentiable solutions are proved for the Cauchy problems. The properties of the constructed solutions are studied by the methods of the qualitative theory of differential equations. Phase portraits are obtained, and quantitative estimates are determined by constructing and analyzing finite difference schemes. The most significant result is that we have shown that all the special cases for incomplete equations take place for the complete equation, and other configurations of diffusion waves do not arise.

Removing the stability limit of the time–space domain explicit finite difference schemes for acoustic modeling with stability condition-based spatial operators

The time step and grid spacing in explicit finite-difference (FD) modeling are constrained by the Courant-Friedrichs-Lewy (CFL) condition. Recently, it has been found that spatial FD coefficients may be designed through simultaneously minimizing the spatial dispersion error and maximizing the CFL number. This allows one to stably use a larger time step or a smaller grid spacing than usually possible. However, when using such a method, only second-order temporal accuracy is achieved. To address this issue, I propose a method to determine the spatial FD coefficients, which simultaneously satisfy the stability condition of the whole wavenumber range and the time–space domain dispersion relation of a given wavenumber range. Therefore, stable modeling can be performed with high-order spatial and temporal accuracy. The coefficients can adapt to the variation of velocity in heterogeneous models. Additionally, based on the hybrid absorbing boundary condition, I develop a strategy to stably and effectively suppress artificial reflections from the model boundaries for large CFL numbers. Stability analysis, accuracy analysis and numerical modeling demonstrate the accuracy and effectiveness of the proposed method.

A Numerical Model with Finite Difference Schemes for Multi-Species Solute Transport in Porous Media

A precise forecast of contaminant and solute transport has an inevitable role in the management of water resources. In accordance with this purpose, in this paper, a novel one-dimensional numerical model for the transport of a decay chain through homogeneous porous media is proposed. To develop the suggested model, two different schemes of the finite difference method, namely the Lax-Wendroff scheme and Fourth-Order scheme, are used. The verification and validation of the established model are examined by the analytical results of three multi-species solute dispersion problems with three- and four-chain members. The total mean square error, L2- and L∞-norms are applied to assess the results. Although analyses show that both schemes provided reliable results, the numerical results of the Lax-Wendroff scheme are more accurate.

Nonstandard Finite Difference Schemes for an SIR Epidemic Model

This paper aims to present two nonstandard finite difference (NFSD) methods to solve an SIR epidemic model. The proposed methods have important properties such as positivity and boundedness and they also preserve conservation law. Numerical comparisons confirm that the accuracy of our method is better than that of other existing standard methods such as the second-order Runge–Kutta (RK2) method, the Euler method and some ready-made MATLAB codes.

Modelling and numerical synchronization of chaotic system with fractional-order operator

Numerical solution of nonlinear chaotic fractional in space reaction–diffusion system is considered in this paper on a large but finite spatial domain size x ∈ [0, L] for L ≫ 0, x = x(x, y) and t ∈ [0, T]. The classical order chaotic ordinary differential equation is formulated by introducing the second-order spatial fractional derivative with order β ∈ (1, 2]. This second order spatial derivative is modelled by using the definition of the Riesz fractional derivative. The method of approximation combines the Fourier spectral method with the novel exponential time difference schemes. The proposed technique is known to have gained spectral accuracy over finite difference schemes. Applicability and suitability of the suggested methods are tested on Rössler chaotic system of recurring interests in one and two dimensions.