All issues

Efficient method for numerical solving of the steady transport equation in x-y-z-geometry has been suggested. The equation is being solved on hexagonal mesh, reflecting real structure of the reactor active zone cross-section. Method of characteristics is used, that inherits all the outcomes from the two-dimensional r-z-geometry calculation. Two variants of the method of characteristics have been applied for solving the transport equation in a cell: method of short characteristics and its conservative modification. It has been confirmed that in three-dimensional geometry conservative method has advantage over pure characteristic and it produces highly accurate solution, especially for quasi-diffusion tensor components.

The results of the experimental and theoretical researches of kinetics of erythrosine penetration into collagenase molecules have represented in this paper. The case with introduction of the compound (fluorescent probe) which has dimers to enzyme solution as an injection has been considered. It was shown that monomers and dimers can penetrate into enzyme molecules with formation complexes monomer — enzyme, dimer- enzyme. Moreover, transformation of probe fluorescence spectra is at each time moment. Spectrum maximum shift, and its form change. At a time, the immobilized dye dimers greatly impact to formation of end fluorescence spectrum. Well correlation between experimental and theoretical results confirms reality of the obtained data.

Functional modeling of blood clotting and fibrin-polymer mesh formation is of a significant value for medical and biophysics applications. Despite the fact of some discrepancies present in simplified functional models their results are of the great interest for the experimental science as a handy tool of the analysis for research planning, data processing and verification. Under conditions of the good correspondence to the experiment functional models can be used as an element of the medical treatment methods and biophysical technologies. The aim of the paper in hand is a modeling of a point system of the fibrin-polymer formation as a multistage polymerization process with a sol-gel transition at the final stage. Complex-value Rosenbroke method of second order (CROS) used for computational experiments. The results of computational experiments are presented and discussed. It was shown that in the physiological range of the model coefficients there is a lag period of approximately 20 seconds between initiation of the reaction and fibrin gel appearance which fits well experimental observations of fibrin polymerization dynamics. The possibility of a number of the consequent $(n = 1–3)$ sol-gel transitions demonstrated as well. Such a specific behavior is a consequence of multistage nature of fibrin polymerization process. At the final stage the solution of fibrin oligomers of length 10 can reach a semidilute state, leading to an extremely fast gel formation controlled by oligomers’ rotational diffusion. Otherwise, if the semidilute state is not reached the gel formation is controlled by significantly slower process of translational diffusion. Such a duality in the sol-gel transition led authors to necessity of introduction of a switch-function in an equation for fibrin-polymer formation kinetics. Consequent polymerization events can correspond to experimental systems where fibrin mesh formed gets withdrawn from the volume by some physical process like precipitation. The sensitivity analysis of presented system shows that dependence on the first stage polymerization reaction constant is non-trivial.

The study tested correctness of three interatomic potentials available in the scientific literature in reproducing a vacancy diffusion in concentrated Fe–Cr alloys by molecular dynamic simulations. It was necessary for further detailed study of vacancy diffusion mechanism in these alloys with Cr content 5–25 at.% at temperatures in the range of 600–1000 K. The analysis of the potentials was performed on alloys models with Cr content 10, 20, 50 at.%. The consideration of the model with chromium content 50 at.% was necessary for further study of diffusion processes in chromium-rich precipitates in these alloys. The formation energies and the atomic mobilities of iron and chromium atoms were calculated and analyzed in the alloys via an artificially created vacancy for all used potentials. A time dependence of mean squared displacement of atoms was chosen as а main characteristic for the analysis of atomic mobilities. The simulation of vacancy formation energies didn’t show qualitative differences between the investigated potentials. The study of atomic mobilities showed a poor reproduction of vacancy diffusion in the simulated alloys by the concentration-dependent model (CDM), which strongly underestimated the mobility of chromium atoms via vacancy in the investigated range of temperature and chromium content. Also it was established, that the two-band model (2BM) of potentials in its original and modified version doesn’t have such drawbacks. This allows one to use these potentials in simulations of vacancy diffusion mechanism in Fe–Cr alloys. Both potentials show a significant dependence of the ratio of chromium and iron atomic mobilities on temperature and Cr content in simulated alloys. The quantitative data of the diffusion coefficients of atoms obtained by these potentials also differ significantly.

A correct model of lipid molecule (distearoylphosphatidylcholine, DSPC) and lipid membrane in water was constructed. Model lipid membrane is stable and has a reliable energy distribution among degrees of freedom. Also after equilibration model system has spatial parameters very similar to those of real DSPC membrane in liquid-crystalline phase. This model was used for studying of lipid membrane permeability to oxygen and water molecules and sodium ion. We obtained the values for transmembrane mobility and diffusion coefficients profiles, which we used for effective permeability coefficients calculation. We found lipid membranes to have significant diffusional resistance to penetration not only by charged particles, such as ions, but also by nonpolar molecules, such as oxygen molecule. We propose theoretical approach for calculation of particle flow across a membrane, as well as methods for estimation of distribution coefficients between bilayer and water phase.

The possibility to simplify the integration of Langevin equation using the variation of correlation between augmentation was researched. The analytical expression for a set of numerical schemes is presented. It’s shown that asymptotic limits for squared velocity depend on step size. The region of convergence and the convergence orders were estimated. It turned out that the incorrect correlation between increments decrease the accuracy down to the level of first-order methods for schemes based on precise solution.

The paper presents the methodology of «line-by-line» calculations of the Earth and atmosphere thermal radiation. Intensity of radiation is computed by numerical integration of the radiative transfer kinetic equation and the system of the angular momentum equations using quasi-diffusion method. Data from HITRAN molecular spectroscopic database [Rothman et al., 2009] are used to calculate the atmosphere optical parameters.

The paper considers models of platelet thrombus formation in blood plasma flow in a cylindrical vessel, based on the “advection–diffusion” equation and the Fokker–Planck equation. The comparison of the results of calculations based on these models is given. Considered models show qualitatively similar behavior at the initial stage of thrombus formation. А detailed investigation of large clots requires models’ improvement.

Theory of anomalous diffusion, which describe a vast number of transport processes with power law mean squared displacement, is actively advancing in recent years. Diffusion of liquids in porous media, carrier transport in amorphous semiconductors and molecular transport in viscous environments are widely known examples of anomalous deceleration of transport processes compared to the standard model.

Direct Monte Carlo simulation is a convenient tool for studying such processes. An efficient stochastic simulation algorithm is developed in the present paper. It is based on simple renewal process with interarrival times that have power law asymptotics. Analytical derivations show a deep connection between this class of random process and equations with fractional derivatives. The algorithm is further generalized by coupling it with chemical reaction simulation. It makes stochastic approach especially useful, because the exact form of integrodifferential evolution equations for reaction — subdiffusion systems is still a matter of debates.

Proposed algorithm relies on non-markovian random processes, hence one should carefully account for qualitatively new effects. The main question is how molecules leave the system during chemical reactions. An exact scheme which tracks all possible molecule combinations for every reaction channel is computationally infeasible because of the huge number of such combinations. It necessitates application of some simple heuristic procedures. Choosing one of these heuristics greatly affects obtained results, as illustrated by a series of numerical experiments.

This paper presents a methodology for modeling centrifugal pumps using the example of the NM 1250 260 main oil centrifugal pump. We use FlowVision CFD software as the numerical modeling instrument. Bench tests and numerical modeling use water as a working fluid. The geometrical model of the pump is fully three-dimensional and includes the pump housing to account for leakages. In order to reduce the required computational resources, the methodology specifies leakages using flow rate rather than directly modeling them. Surface roughness influences flow through the wall function model. The wall function model uses an equivalent sand roughness, and a formula for converting real roughness into equivalent sand roughness is applied in this work. FlowVision uses the sliding mesh method for simulation of the rotation of the impeller. This approach takes into account the nonstationary interaction between the rotor and diffuser of the pump, allowing for accurate resolution of recirculation vortices that occur at low flow rates.

The developed methodology has achieved high consistency between numerical simulations results and experiments at all pump operating conditions. The deviation in efficiency at nominal conditions is 0.42%, and in head is 1.9%. The deviation of calculated characteristics from experimental ones increases as the flow rate increases and reaches a maximum at the far-right point of the characteristic curve (up to 4.8% in head). This phenomenon occurs due to a slight mismatch between the geometric model of the impeller used in the calculation and the real pump model from the experiment. However, the average arithmetic relative deviation between numerical modeling and experiment for pump efficiency at 6 points is 0.39%, with an experimental efficiency measurement error of 0.72%. This meets the accuracy requirements for calculations. In the future, this methodology can be used for a series of optimization and strength calculations, as modeling does not require significant computational resources and takes into account the non-stationary nature of flow in the pump.