Discrete and Continuous Models and Applied Computational Science
Editor-in-Chief: Yuriy P. Rybakov, Doctor of Science (Physics and Mathematics), Professor, Honored Scientist of Russia
ISSN: 2658-4670 (Print). ISSN: 2658-7149 (Online)
Founded in 1993. Publication frequency: quarterly.
Peer-Review: double blind. Publication language: English.
APC: no article processing charge. Open Access: Open Access
, DOAJ SEAL ![]()
PUBLISHER: Peoples’ Friendship University of Russia named after Patrice Lumumba (RUDN University)
See the Journal History to get information on previous journal titles.
Indexation: White List, Russian Index of Science Citation, Scopus (Q3 SJR), VINITI RAS, DOAJ, Google Scholar, Ulrich's Periodicals Directory, WorldCat, Cyberleninka, Dimensions, ResearchBib, Lens, Research4Life, JournalTOCs
Discrete and Continuous Models and Applied Computational Science was created in 2019 by renaming RUDN Journal of Mathematics, Information Sciences and Physics. RUDN Journal of Mathematics, Information Sciences and Physics was created in 2006 by combining the series "Physics", "Mathematics", "Applied Mathematics and Computer Science", "Applied Mathematics and Computer Mathematics".
Discussed issues affecting modern problems of physics, mathematical modeling, computer science. The widely discussed issues Teletraffic theory, queuing systems design, software and databases design and development.
Discussed problems in physics related to quantum theory, nuclear physics and elementary particle physics, astrophysics, statistical physics, the theory of gravity, plasma physics and the interaction of electromagnetic fields with matter, radio physics and electronics, nonlinear optics.
Journal has a high qualitative and quantitative indicators. The Editorial Board consists of well-known scientists of world renown, whose works are highly valued and are cited in the scientific community. Articles are indexed in the Russian and foreign databases. Each paper is reviewed by at least two reviewers, the composition of which includes PhDs, are well known in their circles. Author's part of the magazine includes both young scientists, graduate students and talented students, who publish their works, and famous giants of world science.
Subject areas:
- Mathematics
- Modeling and Simulation
- Mathematical Physics
- Computer Science
- Computer Science (miscellaneous)
Current Issue
Vol 34, No 1 (2026)
- Year: 2026
- Articles: 11
- URL: https://journals.rudn.ru/miph/issue/view/2100
- DOI: https://doi.org/10.22363/2658-4670-2026-34-1
Full Issue
Editorial
Physical dimensional quantities typesetting
Abstract
The siunitx package is designed for typographically correct and consistent typesetting of physical quantities (numbers with units of measurement) in LaTeX documents. It automates formatting according to the rules of the International System of Units (SI), eliminating the need to manually manage spaces, fonts, and separators.
5-11
Computer Science
Usage of polynomial representation of numbers for approximate homomorphic encryption
Abstract
\emph {Introduction} In the modern world of computers and networks the idea of expanding of personal computer resources with the help of cloud storages and computation looks more and more lucrative. However, usage of these resources may endanger data being processed. In last twenty years several algorithms of homomorphic encryption were developed allowing solving of this problem among other applications. However such algorithms are usually constructed as public key systems for long term storage and processing of data. In this article two algorithms of homomorphic encryption optimized for single data processing are proposed. \emph {Purpose} The target of research is development of data coding system which allows safe data processing in public clouds. \emph {Results} Two homomorphic coding systems had been developed, first is based on representation of numbers in the form of polynomials, second based on further representation of polynomials in the form of sets of values. Developed systems allow approximate calculations of coded data without decryption allowing processing of real numbers. System has high level of protection and provides high precision of calculations, comparable with standard personal computer calculation precision. Structure of coded data allows parallel computing. Proposed system allows safe data processing in public networks. Question of finding of optimal parameters for the system stands open both for high precision calculation of limited sets of operations and repeatedly good precision for big sets of operations.
12-23
The waiting time extremal index in GI/G/1 system
Abstract
In this paper the conditions to compare the extremal index of the stationary waiting time in the $M/G/1$ and $GI/M/1$ systems are obtained. These conditions include exponential asymptotic behaviour of waiting time tail and the order in failure rates for the interarrival intervals and for the service times in the systems to be compared. For $M/G/1$ system the obtained result is extended to the mixed service times with ordered components. If, in a $GI/G/1$ system, the service time is determined by a finite mixture whose dominant component of the equilibrium distribution belongs to the class of subexponential distributions then the tail of the limiting distribution of the stationary waiting time is equivalent to the tail of this distribution up to a constant obtained explicitly. Furthermore, the limiting distribution of the maximum of the stationary waiting time belongs to the maximum domain of attraction of the distribution of extreme values of the same type as the maximum of the random variables defined by the dominant component.
24-39
Modeling and Simulation
Derivative-free iterations in $R^n$ with point-wise operations for solving systems of nonlinear equations
Abstract
In this paper, we develop a new family of high-order derivative-free iterative methods for solving systems of nonlinear equations. Specifically, we propose four two-step derivative-free schemes with convergence orders four and five, together with twelve three-step derivative-free schemes achieving convergence orders six, seven, and eight. The main specific of these iterations is that they include a vector or even a scalar iteration parameter instead of the matrix parameter inherent to other existing iterative methods. This structural simplification significantly reduces computational cost, storage requirements, and matrix operations, thereby improving overall computational efficiency. A convergence analysis is presented, establishing the theoretical order of convergence of the proposed methods. The efficiency indices of the proposed schemes are derived and compared with those of several well-known derivative-free iterative methods. The numerical experiments on standard academic problems confirm the theoretical results and demonstrate that the proposed methods are competitive and, in many cases, superior in terms of efficiency and robustness.
40-54
Simulation of the evacuation of passengers and crew from aircraft during a fire on the ground
Abstract
Background. Currently, incidents, including fires on board of aircrafts during takeoff and landing, are becoming more frequent. To address this issue, we introduce new models of fire propagation dynamics and the evacuation process for aircraft passengers, accounting for their physical interactions, along with an integrated model combining such processes as the spread of fire, smoke, and temperature. Nowadays aviation incidents involving onboard fires occur regularly, often resulting in traumas among passengers, as well as material damage. Purpose. The main purpose of this study is to create integrated models that enable analysis of aircraft evacuation under various fire hazard scenarios. Much attention is given to using these models to analyze the process of leaving the aircraft, taking into account various scenarios of the spread of damaging fire factors, which will allow us to develop an optimal sequence of actions for each particular situation. Method. It uses mathematical apparatus of the multi-dimensional cellular automata to describe fire spread, dividing the aircraft into cubic cells with 4 states: burning, burned, consisting of combustible, and non-combustible materials. Calculation of the probabilities of combustion is based on the influence of the neighboring cells, while evacuation models incorporate multi-agent approaches considering passengers' movements, physical contacts, and hazardous factor distributions. The model was created, and graphs were obtained using Python 3.12. Results. The results indicate that the integrated model accurately simulates fire dynamics and evacuation interactions, allowing us to analyze different scenarios to make scenario-based predictions of optimal post-accident exit routes. The model was implemented for two scenarios: a fire in the left engine of the Embraer E-190 and Airbus A320-100 aircraft. Conclusions. Based on the findings, it can be concluded that this approach facilitates decision support systems for enhancing safety during ground-based aircraft fires, providing the model for analyzing and minimizing risks in sudden emergencies.
55-69
Dual quaternion representation of geometrical motion in 3D space
Abstract
\emph {Background} In a previous article we discussed the use of dual quaternions for modeling points, lines and planes and solving standard geometric problems. This article is a logical continuation and reveals the use of dual quaternions to describe isometries of three-dimensional space. \emph {Purpose} The derivation of all necessary formulas for the screw motion of points, straight lines and planes, as well as reflection relative to the plane. Refinement of notation and formalism. \emph {Method} The algebra of dual numbers, quaternions and dual quaternions is used, as well as elements of the theory of screws and sliding vectors. \emph {Results} Formulas for rotation, translation, reflection, helical motion, and mirror rotation are obtained and systematized. \emph {Conclusions} Dual quaternions can serve as a full-fledged tool for describing helical motion in space. Due to the possibility of expressing dual quaternion operations in terms of standard vector and scalar products, the formulas obtained allow for effective software implementation.
70-97
On a finite-difference scheme defining a birational non-quadratic map between time layers
Abstract
The article considers reversible difference schemes for dynamical systems based on the system doubling method proposed by V.N. Abrashin and S.N. Sytova. The method duplicates the original variables, leading to an extended system whose finite-difference approximation defines a birational map between time layers. The preservation of algebraic integrals in such schemes is investigated. It is proved that if the original system admits a homogeneous quadratic first integral, the corresponding bilinear form is exactly preserved by the discrete scheme. This property is demonstrated on the Jacobi oscillator, where the geometric mean of the duplicated variables ensures exact conservation of the quadratic integral. A more detailed analysis is performed on the non-trivial Vanhaecke system, an integrable Hamiltonian system with two degrees of freedom and higher-degree polynomial integrals. Numerical experiments carried out in the computer algebra system Sage using the package fdm.sage confirm that the two copies oscillate synchronously around the exact values of the first integrals, and averaging reduces the oscillation amplitude. For separable Hamiltonian systems, the scheme is shown to be symplectic. The results obtained allow recommending the doubling method for constructing stable and structure-preserving numerical integrators for a wide class of dynamical systems with polynomial right-hand sides, including high-dimensional systems.
98-112
Physics and Astronomy
Interaction of relativistic electrons with intense electromagnetic fields: ponderomotive\,effect, acceleration,\,refraction,\,reflection, dependence\,on\,initial\,conditions
Abstract
The rigorous theory and characterization of charged-particle dynamics in high-intensity electromagnetic fields are fundamental for the development of advanced plasma-based applications. Accurate analytical models must bridge the gap between smoothed trajectories and exact particle motion to predetermine injection and energy gain. The main objective of this review is to establish a rigorous framework for the averaged relativistic motion of electrons, focusing on the strict derivation of ponderomotive forces and the impact of fast-oscillating periodic additions on dynamical variables. By making use of the Krylov--Bogoliubov--Mitropolsky averaging method to obtain the equations of motion, the study analyzes relativistic effects in laser beams and waveguides. These theoretical predictions are substantiated through numerical validation, including test-particle simulations and three-dimensional particle-in-cell simulations (PIC) of relativistic self-trapping regimes such as “laser bullet” and “bubble” structures. The review details the independence of the results on the formulation framework, the strict dependence on wave polarization, and the non-strict potential character of the relativistic ponderomotive force. The analysis demonstrates that periodic fast-oscillating additions are essential for a complete description, accurately setting initial conditions in averaged equations and enabling precise predictions of electron reflection and refraction. Simulations confirm that these fast-oscillating corrections determine electron injection and beam charge in realistic laser–plasma acceleration scenarios. The present review clearly shows that the dual framework of test-particle and PIC models is vital for probing the limits of averaged motion theory. The findings are of direct practical relevance for the optimization of radiation sources and guide the development of future theories incorporating non-adiabatic and field topology dependent effects.
113-124
Mathematical models of low-pressure discharge in a magnetic field supported by UHF electromagnetic field
Abstract
Electron cyclotron resonance (ECR) discharges are an effective way to generate plasma at low working gas pressure. The aim of this work is to develop a mathematical model of the ECR discharge implemented at the RAPIRA facility (RUDN University), which is used for a wide range of scientific research. The evolution of plasma particles is described within the framework of the hydrodynamic approximation (a two-dimensional model with cylindrical symmetry). A three-dimensional model of cold plasma is used to calculate the spatial distribution of the electromagnetic field. Calculations have shown that in the operating mode of the facility (gas pressures from $4\cdot 10^{-4}$ to $10^{-2}$ Torr, magnetic field up to 2500 G), the electron temperature is equalized along the magnetic field lines, and at the same time, the magnetic field ensures a decrease in energy losses to the side walls of the facility. The spatial distributions of the electron density and temperature and the electromagnetic field in the plasma are calculated. The implemented model can serve as a basis for developing a more advanced set of software codes that take into account the non-Maxwellian nature of the electron velocity distribution function, caused by the non-adiabatic nature of their heating in a non-uniform magnetic field.
125-138
Solution of the one-dimensional Schrödinger equation for a heterostructure with a triangular potential function by the power series method
Abstract
In the work by the power series method the one-dimensional Schrödinger equation is solved with a triangular potential function which is applied in various modern heterostructures, in particular for GaAs and the others. By varying available parameters it is possible to obtain the desired precision of the numerical solution of the Schrödinger equation with any type of potential function for modern heterostructures. For the original Schrödinger equation are obtained wave functions in the form Airy functions and the analytical formula for the energy levels through the zeros of the Airy function. The values energy levels from this analytical formula agree with its results obtained by direct power series method with precision up to $10^{{-4}}$ percents, that is, up to 5 decimal signs. However, it is more rational and easier to use the Schrödinger equation solution, because the numerical calculations zeros of Airy function present separate complex and complicated numerical problem. But in order to achieve high numerical accuracy, it is necessary to set the Digits flag to several dozen significant digits and increasing the number of power series, that leads to an increasing in the time spent on the computer.
139-144
Letters
145-149








