Introduction. This paper addresses a unified spatially inhomogeneous, non-stationary model of interaction between genetically modified crop resources (corn) and the corn borer pest, which is also present on a relatively small section of non-modified corn. The model assumes that insect pests influence both types of crops and are capable of independent movement (taxis) towards the gradient of plant resources. It also considers diffusion processes in the dynamics of all components of the unified model, biomass growth, genetic characteristics of both types of plant resources, processes of crop consumption, phenomena of growth and degradation, diffusion, and mutation of pests. The model allows for predictive calculations aimed at reducing crop losses and increasing the resistance of transgenic crops to pests by slowing down the natural mutation rate of the pest.

   Materials and Methods. The mathematical model is an extension of Kostitsin’s model and is formulated as an initial-boundary value problem for a nonlinear system of convection-diffusion equations. These equations describe the spatiotemporal dynamics of biomass density changes in two types of crops — transgenic and non-modified — as well as the specific populations (densities) of three genotypes of pests (the corn borer) resulting from mutations. The authors linearized the convection-diffusion equations by applying a time-lag method on the time grid, with nonlinear terms from each
equation taken from the previous time layer. The terms describing taxis are presented in a symmetric form, ensuring the skew-symmetry of the corresponding continuous operator and, in the case of spatial grid approximation, the finite-difference operator.

   Results. A stable monotonic finite-difference scheme is developed, approximating the original problem with second-order accuracy on a uniform 2D spatial grid. Numerical solutions of model problems are provided, qualitatively corresponding to observed processes. Solutions are obtained for various ratios of modified and non-modified sections of the field.

   Discussion and Conclusion. The obtained results regarding pest behavior, depending on the type of taxis, could significantly extend the time for pests to acquire Bt resistance. The concentration dynamics of pests moving in the direction of the food gradient differs markedly from the concentration of pests moving towards a mate for reproduction.

   Introduction. The paper considers the solution of boundary value problems on an interval for linear ordinary differential equations, in which the coefficients and the right-hand side are continuous functions. The conditions for the orthogonality of the residual equation to the coordinate functions are supplemented by a system of linearly independent boundary conditions. The number of coordinate functions m must exceed the order n of the differential equation.

   Materials and Methods. To numerically solve the boundary value problem, a system of linearly independent coordinate functions is proposed on a symmetric interval [−1,1], where each function has a unit Chebyshev’s norm. A modified Petrov-Galerkin method is applied, incorporating linearly independent boundary conditions from the original problem into the system of linear algebraic equations. An integral quadrature formula with twelfth-order error is used to compute the scalar product of two functions.

   Results. A criterion for the existence and uniqueness of a solution to the boundary value problem is obtained, provided that n linearly independent solutions of the homogeneous differential equation are known. Formulas are derived for the matrix coefficients and the coefficients of the right-hand side in the system of linear algebraic equations for the vector expansion of the solution in terms of the coordinate function system. These formulas are obtained for second- and third-order linear differential equations. The modified Bubnov-Galerkin method is formulated for differential equations of arbitrary order.

   Discussion and Conclusions. The derived formulas for the generalized Bubnov-Galerkin method may be useful for solving boundary value problems involving linear ordinary differential equations. Three boundary value problems with second- and third-order differential equations are numerically solved, with the uniform norm of the residual not exceeding 10^–11.

   Introduction. The study is devoted to the numerical investigation of laser radiation’s effect on an oncoming two-phase flow of nanoparticles and multicomponent hydrocarbon gases. Under such exposure, the hydrogen content in the products increases, and methane is bound into more complex hydrocarbons on the surface of catalytic nanoparticles and in the gas phase. The hot walls of the tube serve as the primary source of heat for the reactive two-phase medium containing catalytic nanoparticles.

   Materials and Methods. The main method used is mathematical modelling, which includes the numerical solution of a system of equations for a viscous gas-dust two-phase medium, taking into account chemical reactions and laser radiation. The model accounts for the two-phase gas-dust medium’s multicomponent and multi-temperature nature, ordinary differential equations (ODEs) for the temperature of catalytic nanoparticles, ODEs of chemical kinetics, endothermic effects of radical chain reactions, diffusion of light methyl radicals CH₃ and hydrogen atoms H, which initiate methane
conversion, as well as absorption of laser radiation by ethylene and particles.

   Results. The distributions of parameters characterizing laminar subsonic flows of the gas-dust medium in an axisymmetric tube with chemical reactions have been obtained. It is shown that the absorption of laser radiation by ethylene in the oncoming flow leads to a sharp increase in methane conversion and a predominance of aromatic compounds in the product output.

   Discussion and Conclusion. Numerical modelling of the dynamics of reactive two-phase media is of interest for the development of theoretical foundations for the processing of methane into valuable products. The results obtained confirm the need for joint use of mathematical modelling and laboratory experiments in the development of new resource-saving and economically viable technologies for natural gas processing.

   Introduction. This paper addresses an initial-boundary value problem for the transport of multifractional suspensions applied to coastal marine systems. This problem describes the processes of transport, deposition of suspension particles, and the transitions between its various fractions. To obtain monotonic finite difference schemes for diffusion-convection problems of suspensions, it is advisable to use schemes that satisfy the maximum principle. When constructing a finite difference scheme that adheres to the maximum principle, it is desirable to achieve second-order spatial accuracy for both
interior and boundary points of the domain under study.

   Materials and Methods. This problem presents certain difficulties when considering the boundaries of the geometric domain, where boundary conditions of the second and third kinds are applied. In these cases, to maintain second-order approximation accuracy, an “extended” grid is introduced (a grid supplemented with fictitious nodes). The guideline
is the approximation of the given boundary conditions using the central difference formula, with the exclusion of the concentration function at the fictitious node from the resulting expressions.

   Results. Second-order accurate finite difference schemes for the diffusion-convection problem of multifractional suspensions in coastal marine systems are constructed.

   Discussion and Conclusion. The proposed schemes are not absolutely stable, and a detailed analysis of stability and convergence, particularly concerning the grid step ratio, remains an important problem that the author plans to address in the future.