Vol 28, No 4 (2020)

This article is devoted to some aspects of using the renovation mechanism (different types of renovation are considered, definitions and brief overview are also given) with one or several thresholds as the mathematical models of active queue management mechanisms. The attention is paid to the queuing systems in which a threshold mechanism with renovation is implemented. This mechanism allows to adjust the number of packets in the system by dropping (resetting) them from the queue depending on the ratio of a certain control parameter with specified thresholds at the moment of the end of service on the device (server) (in contrast to standard RED-like algorithms, when a possible drop of a packet occurs at the time of arrivals of next packets in the system). The models with one, two and three thresholds with different types of renovation are under consideration. It is worth noting that the thresholds determine not only from which place in the buffer the packets are dropped, but also to which the reset of packets occurs. For some of the models certain analytical and numerical results are obtained (the references are given), some of them are only under investigation, so only the mathematical model and current results may be considered. Some results of comparing classic RED algorithm with renovation mechanism are presented.

We implement several explicit Runge-Kutta schemes that preserve quadratic invariants of autonomous dynamical systems in Sage. In this paper, we want to present our package ex.sage and the results of our numerical experiments. In the package, the functions rrk_solve, idt_solve and project_1 are constructed for the case when only one given quadratic invariant will be exactly preserved. The function phi_solve_1 allows us to preserve two specified quadratic invariants simultaneously. To solve the equations with respect to parameters determined by the conservation law we use the elimination technique based on Gröbner basis implemented in Sage. An elliptic oscillator is used as a test example of the presented package. This dynamical system has two quadratic invariants. Numerical results of the comparing of standard explicit Runge-Kutta method RK(4,4) with rrk_solve are presented. In addition, for the functions rrk_solve and idt_solve, that preserve only one given invariant, we investigated the change of the second quadratic invariant of the elliptic oscillator. In conclusion, the drawbacks of using these schemes are discussed.

This paper investigates the waveguide propagation of polarized electromagnetic radiation in a thin-film integral optical waveguide. To describe this propagation, the adiabatic approximation of solutions of Maxwell’s equations is used. The construction of a reduced model for adiabatic waveguide modes that retains all the properties of the corresponding approximate solutions of the Maxwell system of equations was carried out by the author in a previous publication in DCM & ACS, 2020, No 3. In this work, for a special case when the geometry of the waveguide and the electromagnetic field are invariant in the transverse direction. In this case, there are separate nontrivial TEand TM-polarized solutions of this reduced model. The paper describes the parametrically dependent on longitudinal coordinates solutions of problems for eigenvalues and eigenfunctions - adiabatic waveguide TE and TM polarizations. In this work, we present a statement of the problem of finding solutions to the model of adiabatic waveguide modes that describe the stationary propagation of electromagnetic radiation. The paper presents solutions for the single-mode propagation of TE and TM polarized adiabatic waveguide waves.

In this paper, we report a numerical simulation of laser ablation of a material by ultrashort laser pulses. The thermal mechanism of laser ablation is described in terms of a one-dimensional nonstationary heat conduction equation in a coordinate system associated with a moving evaporation front. The laser action is taken into account through the functions of the source in the thermal conductivity equation that determine the coordinate and time dependence of the laser source. For a given dose of irradiation of the sample, the profiles of the sample temperature at different times, the dynamics of the displacement of the sample boundary due to evaporation, the velocity of this boundary, and the temperature of the sample at the moving boundary are obtained. The dependence of the maximum temperature on the sample surface and the thickness of the ablation layer on the radiation dose of the incident laser pulse is obtained. Numerical calculations were performed using the finite difference method. The obtained results agree with the results of other works obtained by their authors.