All issues

This paper analyzes some issues in developing numerical methods for solving problems with a Boltzmann-type linear kinetic transport equation. Existing applications of this type of equation are listed. The focus is on the problem of radiative transfer in a flat layer, which are important for experimental research practice. Key definitions and traditional limitations applied to radiative transfer problems are presented. Some features of formulating radiative transfer problems for flat layers of irregular heterogeneous composite materials that are partially transparent to electromagnetic radiation are considered. The main approaches to the numerical and numerical-analytical solution of the linear kinetic transport equation are outlined.

Some variants of the simplest grid numerical methods for solving of nonstationary kinetic problems of transport a flat layer of a medium with strong attenuation are considered. Problems with one- and two-step variants of these iterative methods are analyzed, for some of them the causes of instability and convergence absence in some of them are investigated and established. It is shown that in the explicit conservative one-step method for a layer of a homogeneous absorbing, but neither radiating nor scattering, medium, unstable modes always exist in the spectrum of harmonic solutions. These modes arise in the region of radiation propagating almost parallel to the layer boundaries, and their instability increases with increasing attenuation effects and is caused by the presence of a small coefficient before the spatial derivative in the transport equation. To limit the undesirable influence of this component, various variants of splitting the equation into two and three fractional steps are considered.

It is shown that the most preferable options are those with explicitly organized fractional steps, for which a proof of their stability and convergence, that based on the Lax’s equivalence theorem is presented. It is demonstrated that the correct building of the fractional step sequence in explicit schemes for numerical solving of the nonstationary linear kinetic transport problems can provide additional stabilization, with the scattering integral plays an important role in stabilizing them. So, when solving kinetic transport problems in media with high scattering albedo, the explicit grid method of settling with splitting the iterations into three fractional steps, that were based on physical processes proved to be the simplest and most effective. The method is implemented as Matlab code, which performs quality control during the generation of the numerical solution process. The most significant modeling results are presented, confirming that the three-step method imposes relatively moderate requirements on resources and numerical integration accuracy, and ensures conditional convergence of iterations. Its mathematical correctness is confirmed by the behavior of the equation residuals and direct control of the convergence of numerical solutions. Its physical correctness is confirmed by ensuring, for ergodic systems, the property of convergence to an invariant steady state independent of the initial conditions. Some discovered and possible limitations of the method are listed.

The work will be useful to specialists in the field of mathematical modeling, numerical methods, kinetic theory, combined heat and mass transfer, dealing with issues of interpretation of experimental data, graduate students and senior students specializing in the indicated areas.

A spatially distributed mathematical model is proposed for the evolution of five myocardial tissue states — healthy myocardium, inflammatory infiltrate, necrosis, replacement fibrosis, and interstitial fibrosis — in the form of a multicomponent parabolic reaction – diffusion system. At each spatial point the tissue state is described by a vector of Dirichlet concentration parameters, which allows one to encode both the expected tissue composition and the uncertainty of classification associated with transition zones of damage. For a truncated system that coincides with the original one inside an invariant region, local Lipschitz continuity and quasi-positivity of the extended reaction operator are established. Global existence of a weak solution, non-negativity, an explicit a priori upper bound, and a quantitative exponential lower bound are proved. These results imply invariance of the admissible region and hence existence of a global weak solution to the original model. Uniqueness is established in a strengthened class of solutions characterized by additional gradient regularity. For an auxiliary reversible transition matrix, a formal entropy balance is derived for smooth positive solutions; for the clinically motivated irreversible transition matrix, a modal spectral stability condition is formulated for the linearized problem. On the numerical side, a positivity condition is derived for the reaction substep of a splitting scheme relative to an adaptive upper trajectory bound, and numerical experiments are conducted for a reproducible calibrated computational variant. In the one-dimensional scenario, by day $60$ the replacement-fibrosis fraction at the lesion center increases from $0.091$ to $0.855$, while the total concentration parameter $\alpha_0$ rises from $11.0$ to $15.9$; in the two-dimensional scenario, on day $14$ the mean value of $\overline p_4^{}$ is $0.431$ in the lesion core versus $0.097$ in the healthy region. These results confirm the mathematical well-posedness of the proposed model and its applicability for the quantitative description of post-infarction scar formation and heterogeneous transition zones.

Modern power transportation systems are the complex engineering systems. Such systems include both point facilities (power producers, consumers, transformer substations, etc.) and the distributed elements (f.e. power lines). Such structures are presented in the form of the graphs with different types of nodes under creating the mathematical models. It is necessary to solve the system of partial differential equations of the hyperbolic type to study the dynamic effects in such systems.

An approach similar to one already applied in modeling similar problems earlier used in the work. New variant of the splitting method was used proposed by the authors. Unlike most known works, the splitting is not carried out according to physical processes (energy transport without dissipation, separately dissipative processes). We used splitting to the transport equations with the drain and the exchange between Reimann’s invariants. This splitting makes possible to construct the hybrid schemes for Riemann invariants with a high order of approximation and minimal dissipation error. An example of constructing such a hybrid differential scheme is described for a single-phase power line. The difference scheme proposed is based on the analysis of the properties of the schemes in the space of insufficient coefficients.

Examples of the model problem numerical solutions using the proposed splitting and the difference scheme are given. The results of the numerical calculations shows that the difference scheme allows to reproduce the arising regions of large gradients. It is shown that the difference schemes also allow detecting resonances in such the systems.