Mathematical modeling of two-phase compressible fluid filtration based on modified adaptive method of minimum amendments

The work objective is to build and investigate the modified adaptive method of minimum amendments (MAMMA) which is destined for the numerical simulation of the two-phase compressible fluid filtration in porous media. This approach allows overcoming the known use limitations of other methods of the finite-difference equations solution, such as: crucial differential pressures acting on the oil-and-water bearing formation; and the compressibility of the medium at the considerable gas content in the oil phase. An approximation method - an explicit one for defining the function of water saturation, and an implicit one for the pressure function computation - is selected as the research basis. When setting the initial boundary value problem and its sampling, the process of the two-phase compressible fluid filtration in the space-dimensional domain with the lateral area bounded below by the subface of stratum, and above - by the bed top, is considered. A two-layer iterative method of the variational type - a modified method of minimal amendments adapted for solving finite-difference equations of the two-phase compressible fluid with a non-selfadjoint operator under the most general assumptions on the properties of the grid-problem operator is built. It is shown that a MAMMA has the asymptotic convergence rate characteristic of the “classical” alternate triangular method that does not use the Chebyshev acceleration technique and can be applied to the problems with a self-adjoint operator. Numerical experiments have confirmed the high efficiency of MAMMA. It is established that to achieve the specified accuracy, the number of iterations at the MAMMA reduces to 3-20 times as compared to the method of Seidel and the overrelaxation method

Sufficient conditions for convergence of positive solutions to linearized two-dimensional sediment transport problem

Introduction. The sediment transportation is one of the major processes that define the magnitude and back surface changing rate for water bodies. The most used prognostic studies in this field are based on the mathematical models that allow reducing, and in some cases, eliminating the expensive and often environmentally burdensome experiments. Spatially one-dimensional models are mainly used to predict changes in bottom topography. For actual coastal systems with irregular coastal configuration, the flow vector is generally not orthogonal to the tangent line for the coastline at each of its points. It also may not coincide with the wind stress vector. Therefore, to solve lots of practically important problems associated with the prediction of the dynamics of the back surface of water basins, it is necessary to use spatially two-dimensional models of sediment transportation and effective numerical methods of their implementation. Materials and Methods. The authors (A.I. Sukhinov, A.E. Chistyakov, E.A. Protsenko, and V.V. Sidoryakina) have earlier proposed a spatially two-dimensional model of sediment transport that satisfies the basic conservation laws (of material balance and momentum) which is a quasilinear equation of parabolic type. The linear difference schemes are constructed and studied; model and some practically important tasks are solved. However, the theoretical study on the proximity of solutions for the original nonlinear initial-boundary value problem and the linearized continuous task - on which basis a discrete model (difference scheme) was built - remained in the shadow. The study of the linearized problem correctness and the determination of sufficient conditions for positivity of solutions are of special interest because only positive solutions to this sediment transport problem have physical value within the framework of the considered models. Research Results. The investigated nonlinear two-dimensional model of sediment transport in the coastal zone of shallow reservoirs takes account of the following physically significant factors and parameters: soil porosity; critical value of the tangent stress at which load transport starts; turbulent interaction; dynamically varying of the bottom geometry; wind currents; and bottom friction. Linearization is carried out on the time grid - nonlinear coefficients of the parabolic equation are taken for the previous time grid step. Next, a chain of tasks connected by the initial data - final solutions of the linearized mixed Cauchy problems chain on a uniform time grid is constructed, and thus, the linearization for the initial 2D nonlinear model is carried out in large. Earlier, the authors have proved the existence and uniqueness of the solution to a linear tasks chain. Prior estimate of the proximity of the linearized problem chain solution to the initial non-linear task solution has been also obtained. The conditions of its solution positivity and their convergence to the nonlinear sediment transport problem are defined in the norm of the Hilbert space L1 with the rate O(τ) where τ is a time step. Discussion and Conclusions. The obtained research results of the spatially two-dimensional nonlinear sediment transport model can be used for predicting the nonlinear hydrodynamic processes, improving their accuracy and reliability due to the availability of new accounting functionality of physically important factors, including the specification of the boundary conditions

Development and correctness analysis of the mathematical model of transport and suspension sedimentation depending on bottom relief variation

Introduction. The paper is devoted to the study on the three-dimensional model of transport and suspension sedimentation in the coastal area due to changes in the bottom relief. The model considers the following processes: advective transfer caused by the aquatic medium motion, micro-turbulent diffusion, and gravity sedimentation of suspended particles, as well as the bottom geometry variation caused by the particle settling or bottom sediment rising. The work objective was to conduct an analytical study of the correctness of the initial-boundary value problem corresponding to the constructed model.Materials and Methods. The change in the bottom relief aids in solution to the initial-boundary value problem for a parabolic equation with the lowest derivatives in a domain whose geometry depends on the desired function of the solution, which in general leads to a nonlinear formulation of the problem. The model is linearized on the time grid due to the “freezing” of the bottom relief within a single step in time and the subsequent recalculation of the bottom surface function on the basis of the changed function of the suspension concentration, as well as a possible change in the velocity vector of the aquatic medium.Research Results. For the linearized problem, a quadratic functional is constructed, and the uniqueness of the solution to the corresponding initial boundary value problem is proved within the limits of an unspecified time step. On the basis of the quadratic functional transformation, we obtain a prior estimate of the solution norm in the functional space L2 as a function of the integral time estimates of the right side, and the initial condition. Thus, the stability of the solution to the initial problem from the change of the initial and boundary conditions, the right-hand side function, is established.Discussion and Conclusions. The model can be of value for predicting the spread of contaminants and changes in the bottom topography, both under an anthropogenic impact and due to the natural processes in the coastal area

Small Innovative Enterprises at Universities: Barriers and Opportunities for Development

The given article is of research character and aims at identifying the problems and development mechanism prospects of commercialization of university intellectual activity results through the establishment of small innovative enterprises. The research is based on the analysis of small innovative enterprises creation and functioning mechanism carried out on the basis of empirical data of Don State Technical University. The applied methods of systemic and situational analysis increase the objectivity of the results obtained. It is demonstrated that the main problems of small innovative enterprises functioning have organizational and financial character. The prospects of small innovative enterprises are seen in university innovation ecosystems development allowing creating the necessary environment for their functioning and, thus, improving their attractiveness. The results obtained can be useful to the heads of the universities and public authority representatives when updating or elaborating development programs and strategies.Данная статья носит исследовательский характер и направлена на выявление проблем и перспектив развития механизма коммерциализации результатов интеллектуальной деятельности университетов на основе создания ими малых инновационных предприятий. В основе исследования лежит анализ механизма создания и функционирования малых инновационных предприятий, выполненный на базе эмпирических данных Донского государственного технического университета. Применяемые методы системного и ситуационного анализа повышают объективность полученных результатов. Показано, что основные проблемы функционирования малых инновационных предприятий имеют организационный и финансовый характер. Перспективы малых инновационных предприятий видятся в развитии инновационных экосистем университетов, позволяющих сформировать необходимую среду для их функционирования и повысить их привлекательность. Полученные результаты исследования могут быть полезны руководителям университетов и представителям органов государственной власти при актуализации или разработке ими программ и стратегий развития

Supercomputer Simulation of Oil Spills in the Azov Sea

We present the research on microbiological destruction of oil pollution in shallow water. In order to conduct the research, we use a multiprocessor computer system with distributed memory. The research takes into account the oil fractional composition as well as hydrodynamic and chemical-biological features of water. In order to simulate the dynamics of hydrocarbon microbiological degradation in the Azov Sea, we propose the complex of interrelated precision models. For model discretization, we use the space splitting schemes taking into account the partial filling cells of computational domain. Therefore, the computational accuracy significantly increases, while the computational time decreases. On supercomputer, we implement an experimental software for predictive modelling the ecological situation under oil and other pollution conditioned by natural and industrial challenges in shallow water.Исследование микробиологической деструкции нефтяных загрязнений в мелководном водоеме с учетом фракционного состава нефти, а также гидродинамических и химико-биологических свойств воды проводилось на многопроцессорной вычислительной системе с распределенной памятью. Предложен комплекс взаимосвязанных прецизионных моделей для моделирования динамики микробиологической деградации углеводородов в Азовском море. Для дискретизации моделей использовались схемы расщепления по пространству с учетом частичной заполненности ячеек расчетной области, благодаря чему точность вычислений была значительно повышена, а время выполнения расчетов сократилось. На супер-ЭВМ разработано экспериментальное программное обеспечение, предназначенное для прогнозного моделирования экологической ситуации в мелководном водоеме при нефтяном и других загрязнениях при природных и индустриальных вызовах

Practical Aspects of Implementation of the Parallel Algorithm for Solving Problem of Ctenophore Population Interaction in the Azov Sea

The paper covers the development and researching mathematical model of interaction processes between plankton and ctenophore populations based on the modern information technologies and computational methods, which leads to increase of the accuracy of predictive modeling of the ecology situation in shallow water in summer. The model takes into account the following: the transport of water environment; microturbulent diffusion; nonlinear interaction of plankton and ctenophore populations; biogenic, temperature and oxygen regimes; influence of salinity. The computational accuracy is significantly increased, and computational time is decreased at using the calculation method based on partially filled cells for discretization of model. The practical significance is the software implementation of the proposed model, the limits and prospects of its practical use are defined. Experimental software was developed based on multiprocessor computer system, which is intended for mathematical modeling of possible progress scenarios in shallow waters ecosystems on the example of the Azov Sea in summer. We used decomposition methods of grid domains in parallel implementation for computationally laborious convection-diffusion problems, taking into account the architecture and parameters of multiprocessor computer system.The paper covers the development and researching mathematical model of interaction processes between plankton and ctenophore populations based on the modern information technologies and computational methods, which leads to increase of the accuracy of predictive modeling of the ecology situation in shallow water in summer. The model takes into account the following: the transport of water environment; microturbulent diffusion; nonlinear interaction of plankton and ctenophore populations; biogenic, temperature and oxygen regimes; influence of salinity. The computational accuracy is significantly increased, and computational time is decreased at using the calculation method based on partially filled cells for discretization of model. The practical significance is the software implementation of the proposed model, the limits and prospects of its practical use are defined. Experimental software was developed based on multiprocessor computer system, which is intended for mathematical modeling of possible progress scenarios in shallow waters ecosystems on the example of the Azov Sea in summer. We used decomposition methods of grid domains in parallel implementation for computationally laborious convection-diffusion problems, taking into account the architecture and parameters of multiprocessor computer system

Modeling and numerical analysis of the effect of dissociation/recombination of water molecules on the transport of salt ions in diffusion layer

Introduction. The paper presents a theoretical study on binary salt ion transport considering the water dissociation/recombination reaction. The work objectives are as follows: to build a mathematical model; to develop an algorithm for the numerical solution to the boundary value problem corresponding to the mathematical model; to work out the similarity theory including the transition to a dimensionless form using characteristic quantities; to determine a physical meaning of trivial similarity criteria; to find nontrivial similarity criteria; to build and analyze the volt-ampere characteristic (VAC).Materials and Methods. The theoretical study and numerical analysis of the transport of binary salt ions consider the dissociation/recombination reaction of water. In this case, the heat transfer equation and the mathematical model of electrodiffusion of four types of ions simultaneously (two salt ions, as well as ????+ and ????????−ions) in the diffusion layer of electromembrane systems with a perfectly selective membrane are used. For the first-order differential equations, a singularly perturbed boundary-value problem is set. In the equation for the electric field, the right side is independent of the intensity. In the numerical solution to the digitized system of equations by the Newton-Kantorovich method, this causes the stability of the method. In this regard, the boundary-value problem is reduced for numerical solution: a transition to a system of the second-order equations is provided, and the missing boundary conditions for the electric field strength are calculated.Research Results. A new mathematical model, a numerical algorithm to solve a boundary value problem, and software are developed. A numerical analysis is carried out, and fundamental laws of the transport of salt ions are determined considering the dissociation/recombination reaction of water molecules, temperature effects, and Joule heating. The VAC is built and analyzed.Discussion and Conclusions. The transport of binary salt ions through a diffusion layer near a cation exchange membrane is considered. A mathematical model of this process is proposed. It takes into account the temperature effects due to dissociation/recombination reactions of water molecules and Joule heating in a solution. The basic laws of the transport of salt ions are established considering the dissociation/recombination reaction of water molecules and temperature effects. The temperature effects of the dissociation/recombination reaction and the Joule heating in the electroneutrality region (ENR) are almost imperceptible (with the exception of the recombination region, RR). The Joule heating in the space-charge region (SCR) is by two orders of magnitude larger than the cooling effect of the water dissociation reaction. Upon recombination, approximately the same heat is released in the RR as during Joule heating in the expanded SCR. However, due to the small size of the RR, the effect of this heat is imperceptible. Therefore, we can assume that there is only one heat source at the interface in the SCR, which, due to its noticeable size, causes a significant increase in temperature in the entire diffusion layer. It follows that the emergence and development of gravitational convection is possible. General conclusions, following from the results obtained, open up the possibility of intensifying the process of transport of salt ions in the electrodialysis machines