Numerical Solutions of Space-Fractional Advection–Diffusion–Reaction Equations

Background: solute transport in highly heterogeneous media and even neutron diffusion in nuclear environments are among the numerous applications of fractional differential equations (FDEs), being demonstrated by field experiments that solute concentration profiles exhibit anomalous non-Fickian growth rates and so-called “heavy tails”. Methods: a nonlinear-coupled 3D fractional hydro-mechanical model accounting for anomalous diffusion (FD) and advection–dispersion (FAD) for solute flux is described, accounting for a Riesz derivative treated through the Grünwald–Letnikow definition. Results: a long-tailed solute contaminant distribution is displayed due to the variation of flow velocity in both time and distance. Conclusions: a finite difference approximation is proposed to solve the problem in 1D domains, and subsequently, two scenarios are considered for numerical computations.

A problem of non – linear deformation of five–layer conical shells with allowance for discrete ribs

A problem of non – linear deformation of multiplayer conical shells with allowance for discrete ribs under non – stationary loading is considered. The system of non – linear differential equations is based on the Timoshenko type theory of rods and shells. The Reissner’s variational principle is used for deductions of the motion equations. An efficient numerical method with using Richardson type finite difference approximation for solution of problems on nonstationary behaviour of multiplayer shells of revolution with allowance distcrete ribs which permit to realize solution of the investigated wave problems with the use of personal computers. As a numerical example, the problem of dynamic deformation of a five-layer conical shell with rigidly clamped ends under the action of an internal distributed load was considered.

A high order finite-difference scheme for time-domain modelling of time-varying seismoelectric waves

Simulation of the seismoelectric effect serves as a useful tool to capture the observed seismoelectric conversion phenomenon in porous media, thus offering promising potential in underground exploration activities to detect pore fluids such as water, oil and gas. The static electromagnetic (EM) approximation is among the most widely used methods for numerical simulation of the seismoelectric responses. However, the static approximation ignores the accompanying electric field generated by the shear wave, resulting in considerable errors when compared to analytical results, particularly under high salinity conditions. To mitigate this problem, we propose a spatial high-order finite-difference time-domain (FDTD) method based on Maxwell's full equations of time-varying EM fields to simulate the seismoelectric response in 2D mode. To improve the computational efficiency influenced by the velocity differences between seismic and electromagnetic waves, different time steps are set according to the stability conditions, and the seismic feedback values of EM time nodes are obtained by linear approximation within the seismic unit time step. To improve the simulation accuracy of the seismoelectric response with the time-varying EM calculation method, finite-difference coefficients are obtained by solving the spatial high-order difference approximation based on Taylor expansion. The proposed method yields consistent simulation results compared to those obtained from the analytical method under different salinity conditions, thus indicating its validity for simulating seismoelectric responses in porous media. We further apply our method to both layered and anomalous body models and extend our algorithm to 3D. Results show that the time-varying EM calculation method could effectively capture the reflection and transmission phenomena of the seismic and EM wavefields at the interfaces of contrasting media. This may allow for the identification of abnormal locations, thus highlighting the capability of seismoelectric response simulation to detect subsurface properties.

Accelerated nonstandard finite difference method for singularly perturbed Burger-Huxley equations

Objective The main purpose of this paper is to present an accelerated nonstandard finite difference method for solving the singularly perturbed Burger-Huxley equation in order to produce more accurate solutions. Results The quasilinearization technique is used to linearize the nonlinear term. A nonstandard methodology of Mickens procedure is used in the spatial direction and also within the first order temporal direction that construct the first-order finite difference approximation to solve the considered problem numerically. To accelerate the rate of convergence from first to second-order, the Richardson extrapolation technique is applied. Numerical experiments were conducted to support the theoretical results.

Mathematical Analysis and Optimal Control of Giving up the Smoking Model

In this study, we are going to explore mathematically the dynamics of giving up smoking behavior. For this purpose, we will perform a mathematical analysis of a smoking model and suggest some conditions to control this serious burden on public health. The model under consideration describes the interaction between the potential smokers P , the occasional smokers L , the chain smokers S , the temporarily quit smokers Q T , and the permanently quit smokers Q P . Existence, positivity, and boundedness of the proposed problem solutions are proved. Local stability of the equilibria is established by using Routh–Hurwitz conditions. Moreover, the global stability of the same equilibria is fulfilled through using suitable Lyapunov functionals. In order to study the optimal control of our problem, we will take into account a two controls’ strategy. The first control will represent the government prohibition of smoking in public areas which reduces the contact between nonsmokers and smokers, while the second will symbolize the educational campaigns and the increase of cigarette cost which prevents occasional smokers from becoming chain smokers. The existence of the optimal control pair is discussed, and by using Pontryagin minimum principle, these two optimal controls are characterized. The optimality system is derived and solved numerically using the forward and backward difference approximation. Finally, numerical simulations are performed in order to check the equilibria stability, confirm the theoretical findings, and show the role of optimal strategy in controlling the smoking severity.

Mathematical Model for Metal Transfer Study in Additive Manufacturing with Electron Beam Oscillation

A computer model has been developed to investigate the processes of heat and mass transfer under the influence of concentrated energy sources on materials with specified thermophysical characteristics, including temperature-dependent ones. The model is based on the application of the volume of fluid (VOF) method and finite-difference approximation of the Navier–Stokes differential equations formulated for a viscous incompressible medium. The “predictor-corrector” method has been used for the coordinated determination of the pressure field which corresponds to the continuity condition and the velocity field. The modeling technique of the free liquid surface and boundary conditions has been described. The method of calculating surface tension forces and vapor recoil pressure has been presented. The algorithm structure is given, the individual modules of which are currently implemented in the Microsoft Visual Studio environment. The model can be applied for studying the metal transfer during the deposition processes, including the processes with electron beam spatial oscillation. The model was validated by comparing the results of computational experiments and images obtained by a high-speed camera.

Interpolating Stabilized Element Free Galerkin Method for Neutral Delay Fractional Damped Diffusion-Wave Equation

A numerical solution for neutral delay fractional order partial differential equations involving the Caputo fractional derivative is constructed. In line with this goal, the drift term and the time Caputo fractional derivative are discretized by a finite difference approximation. The energy method is used to investigate the rate of convergence and unconditional stability of the temporal discretization. The interpolation of moving Kriging technique is then used to approximate the space derivative, yielding a meshless numerical formulation. We conclude with some numerical experiments that validate the theoretical findings.

Mathematical modeling and analysis of operation of cylindric pyrolysis reactor with radial heating

The vector of development of solid-fuel energy is currently directed towards expanding the range of renewable fuels used. Along with the direct combustion of fuel, the processes of controlled thermal transformation of the raw biomass in an oxygen-free surrounding to obtain a new fuel based on it (liquid, solid, gaseous) are widely spread. A significant part of research in this sphere is related to the study of the formal kinetics of such processes, at the same time, the hardware design of the process is no less important, but less studied. Thus, development of mathematical models of pyrolysis equipment operation is relevant. A decisive difference approximation of these processes in the framework of an axisymmetric formulation of the problem is chosen as a mathematical basis for modeling physical and chemical transformations and transfer processes in the radial direction of a cylindrical pyrolysis reactor. The material constants of the processes are borrowed from the well-known literature references The authors studied the modes of reactor operation not covered by a full-scale experiment, using the previously proposed and verified one-dimensional mathematical model of a cylindrical pyrolysis reactor. The issues of the influence of the dimensionless kinetic function of the process (reaction model) on the thermal transformation of the material in the apparatus are considered. The significant influence of the chosen reaction model on the kinetic nature of the process is pointed out. The mutual influence of drying and pyrolysis the presence of which is due to the energy effects of the processes is considered. A significant spatial heterogeneity of the process is defined and the possibility of the existence of a non-trivial effect of advanced heating of the internal zones of the apparatus in comparison with the peripheral ones is specified. The paper shows that a computational experiment can help to detect non-trivial effects and identify the variability of the process implementation even within the framework of a single design and technological solution of the pyrolysis process. According to the authors, the results of the obtained numerical experiments indicate that mathematical modeling can be the basis of making technological solution. However, further research is also needed to determine reliably the material constants of the process.