Все выпуски

Для трубчатых пороховых элементов большого удлинения, используемых в артиллерийских метательных зарядах, имеют место условия неравномерного горения. Здесь необходимо параллельно рассматривать процессы горения и движения пороховых газов внутри и вне каналов пороховых трубок. Без этого невозможно адекватно поставить и решить задачи о воспламенении, эрозионном горении и напряженно-деформированном состоянии трубчатых пороховых элементов в процессе выстрела. В работе представлена физико-математическая постановка основной задачи внутренней баллистики артиллерийского выстрела для заряда, состоящего из совокупности пороховых трубок. Горение и движение пучка пороховых трубок по каналу ствола моделируются эквивалентным трубчатым зарядом всестороннего горения. Площади торца и сечения канала такого заряда (эквивалентной трубки) равны сумме площадей торцов и сечений каналов пороховых трубок соответственно. Поверхность горения канала равна сумме внутренних поверхностей трубок в пучке. Внешняя поверхность горения эквивалентной трубки равна сумме внешних поверхностей трубок в пучке. Предполагается, что эквивалентная трубка движется по оси канала ствола. Скорость движения эквивалентного трубчатого заряда и его текущее положение определяются из второго закона Ньютона. Для расчета параметров течения использованы двумерные осесимметричные уравнения газовой динамики, для решения которых строится осесимметричная ортогонализированная разностная сетка, адаптирующаяся к условиям течения. При перемещении и горении трубки разностная сетка перестраивается с учетом изменяющихся областей интегрирования. Для численного решения системы газодинамических уравнений применяется метод контрольного объема. Параметры газа на границах контрольных объемов определяются с использованием автомодельного решения задачи о распаде произвольного разрыва С.К. Годунова. Разработанная методика использована при расчетах внутрибаллистических параметров артиллерийского выстрела. Данный подход рассмотрен впервые и позволяет по-новому подойти к проектированию трубчатых артиллерийских зарядов, поскольку позволяет получить необходимую информацию в виде полей скорости и давления пороховых газов для расчета процесса постепенного воспламенения, нестационарного эрозионного горения, напряженно-деформированного состояния и прочности пороховых элементов при выстреле. Представлены временные зависимости параметров внутрибаллистического процесса и распределения основных параметров течения продуктов горения в различные моменты времени.

There are conditions of uneven combustion for tubular powder elements of large elongation used in artillery propelling charges. Here it is necessary to consider in parallel the processes of combustion and movement of powder gases inside and outside the channels of the powder tubes. Without this, it is impossible to adequately formulate and solve the problems of ignition, erosive combustion and stress-strain state of tubular powder elements in the shot process. The paper presents a physical and mathematical formulation of the main problem of the internal ballistics of an artillery shot for a charge consisting of a set of powder tubes. Combustion and movement of a bundle of powder tubes along the barrel channel is modeled by an equivalent tubular charge of all-round combustion. The end and cross-sectional areas of the channel of such a charge (equivalent tube) are equal to the sum of the areas of the ends and cross-sections of the channels of the powder tubes, respectively. The combustion surface of the channel is equal to the sum of the inner surfaces of the tubes in the bundle. The outer combustion surface of the equivalent tube is equal to the sum of the outer surfaces of the tubes in the bundle. It is assumed that the equivalent tube moves along the axis of the bore. The speed of motion of an equivalent tubular charge and its current position are determined from Newton’s second law. To calculate the flow parameters, we used two-dimensional axisymmetric equations of gas dynamics, for the solution of which an axisymmetric orthogonalized difference mesh is constructed, which adapts to the flow conditions. When the tube moves and burns, the difference grid is rearranged taking into account the changing regions of integration. The control volume method is used for the numerical solution of the system of gas-dynamic equations. The gas parameters at the boundaries of the control volumes are determined using a self-similar solution to the Godunov problem of decay for an arbitrary discontinuity. The developed technique was used to calculate the internal ballistics parameters of an artillery shot. This approach is considered for the first time and allows a new approach to the design of tubular artillery charges, since it allows obtaining the necessary information in the form of fields of velocity and pressure of powder gases for calculating the process of gradual ignition, unsteady erosive combustion, stress-strain state and strength of powder elements during the shot. The time dependences of the parameters of the internal ballistics process and the distribution of the main parameters of the flow of combustion products at different times are presented.

В настоящей работе исследован турбулентный поток Тейлора – Куэтта с помощью двухмерного моделирования на базе осредненных уравнений Навье – Стокса (RANS) и нового двухжидкостного подхода к турбулентности при числах Рейнольдса в диапазоне от 1000 до 8000. Исследуется течение, обусловленное вращающимся внутренним и неподвижным внешним цилиндрами. Рассмотрен случай соотношения диаметров цилиндров 1:2. Известно, что возникающее круговое течение характеризуется анизотропной турбулентностью и математическое моделирование таких потоков является сложной задачей. Для описания таких потоков используются либо методы прямого моделирования, которые требуют больших вычислительных затрат, либо достаточно трудоемкие методы рейнольдсовых напряжений или же линейные RANS-модели со специальными поправками на вращение, которые способны описывать анизотропную турбулентность. В работе для сравнения различных подходов к моделированию турбулентности представлены численные результаты линейных RANS-моделей SARC, SST-RC, метода рейнольдсовых напряжений SSG/LRR-RSM-w2012, прямого моделирования турбулентности DNS, а также новой двухжидкостной модели. Показано, что недавно разработанная двухжидкостная модель адекватно описывает рассматриваемый поток. Помимо этого, двухжидкостная модель проста для численной реализации и имеет хорошую сходимость.

In this paper, the turbulent Taylor – Couette flow is investigated using two-dimensional modeling based on the averaged Navier – Stokes (RANS) equations and a new two-fluid approach to turbulence at Reynolds numbers in the range from 1000 to 8000. The flow due to a rotating internal and stationary external cylinders. The case of ratio of cylinder diameters 1:2 is considered. It is known that the emerging circular flow is characterized by anisotropic turbulence and mathematical modeling of such flows is a difficult task. To describe such flows, either direct modeling methods are used, which require large computational costs, or rather laborious Reynolds stress methods, or linear RANS models with special corrections for rotation, which are able to describe anisotropic turbulence. In order to compare different approaches to turbulence modeling, the paper presents the numerical results of linear RANS models SARC, SST-RC, Reynolds stress method SSG/LRR-RSM-w2012, DNS direct turbulence modeling, as well as a new two-fluid model. It is shown that the recently developed twofluid model adequately describes the considered flow. In addition, the two-fluid model is easy to implement numerically and has good convergence.

Разработана двумерная математическая модель для оценки напряжений в сварных соединениях, формируемых при многопроходной сварке многослойных сталей. Основой модели является система уравнений, которая включает вариационное уравнение Лагранжа инкрементальной теории пластичности и вариационное уравнение теплопроводности, выражающее принцип М. Био. Вариационно-разностным методом решается задача теплопроводности для расчета нестационарного температурного поля, а затем на каждом шаге по времени – квазистатическая задача термопластичности. Разностная схема построена на треугольных сетках, что дает некоторое повышение точности при описании положения границ раздела структурных элементов.

A two-dimensional mathematical model was developed for estimating the stresses in welded joints formed during multipass welding of multilayer steels. The basis of the model is the system of equations that includes the Lagrange variational equation of incremental plasticity theory and the variational equation of heat conduction, which expresses the principle of M. Biot. Variational-difference method was used to solve the problems of heat conductivity and calculation of the transient temperature field, and then at each time step – for the quasi-static problem of thermoplasticity. The numerical scheme is based on triangular meshes, which gives a more accuracy in describing the boundaries of structural elements as compared to rectangular grids.

В работе рассмотрено приложение методов кинетической теории к задачам гемодинамики. Для моделирования выбраны решеточные уравнения Больцмана. Данные модели описывают дискретизированную по пространственной и временной координате динамику движения частиц на одномерной решетке. Хорошо известно, что в пределе малых длин свободного пробега решеточные уравнения Больцмана описывают уравнения гидродинамики. Если течение достаточно медленное (мало число Маха), то данные уравнения гидродинамики переходят в уравнения Навье – Стокса для сжимаемого газа. Если в получающихся гидродинамических уравнениях переменные, отвечающие плотности и скорости звука, считать площадью поперечного сечения сосуда и скоростью распространения пульсовой волны давления, то выводятся хорошо известные в биомеханике нелинейные уравнения распространения несжимаемой вязкой жидкости (крови) в эластичном сосуде для частного случая постоянной пульсовой скорости.

В общем случае скорость распространения пульсовой волны зависит от площади просвета сосуда. Следует отметить интересную аналогию: уравнение состояния решеточного газа в новых переменных становится законом, связывающим давление и площадь поперечного сечения сосуда. Таким образом, в общем случае требуется модифицировать уравнение состояния для решеточного уравнения Больцмана. Данная процедура хорошо известна в теории неидеального газа и многофазных течений и эквивалентна введению в уравнения виртуальной силы. Получающиеся уравнения могут использоваться для моделирования любых законов, связывающих скорость пульсовой волны и площадь просвета сосуда.

В качестве тестовых задач рассмотрено распространение уединенной нелинейной пульсовой волны в сосуде с упругими свойствами, описываемыми законом Лапласа. Во второй задаче рассмотрено распространение пульсовых волн для бифуркации сосудов. Показано, что результаты расчетов хорошо совпадают с данными из предыдущих исследований.

In the present paper the application of the kinetic methods to the blood flow problems in elastic vessels is studied. The Lattice Boltzmann (LB) kinetic equation is applied. This model describes the discretized in space and time dynamics of particles traveling in a one-dimensional Cartesian lattice. At the limit of the small times between collisions LB models describe hydrodynamic equations which are equivalent to the Navier – Stokes for compressible if the considered flow is slow (small Mach number). If one formally changes in the resulting hydrodynamic equations the variables corresponding to density and sound wave velocity by luminal area and pulse wave velocity then a well-known 1D equations for the blood flow motion in elastic vessels are obtained for a particular case of constant pulse wave speed.

In reality the pulse wave velocity is a function of luminal area. Here an interesting analogy is observed: the equation of state (which defines sound wave velocity) becomes pressure-area relation. Thus, a generalization of the equation of state is needed. This procedure popular in the modeling of non-ideal gas and is performed using an introduction of a virtual force. This allows to model arbitrary pressure-area dependence in the resulting hemodynamic equations.

Two test case problems are considered. In the first problem a propagation of a sole nonlinear pulse wave is studied in the case of the Laplace pressure-area response. In the second problem the pulse wave dynamics is considered for a vessel bifurcation. The results show good precision in comparison with the data from literature.

В данной работе с помощью алгоритма Метрополиса авторами были изучены магнитные системы, в которых из-за конкуренции между прямым гейзенберговским обменом и взаимодействием Дзялошинского–Мория возникают магнитные вихревые структуры — скирмионы.

В статье рассматриваются условия зарождения и стабильного существования магнитных скирмионов в двумерных магнитных пленках в рамках классической модели Гейзенберга. Изучена термическая стабильность скирмионов в магнитной пленке. Были рассмотрены процессы формирования различных состояний в изучаемой системе при варьировании величины внешнего магнитного поля, выделены различные фазы, в которые переходит система спинов Гейзенберга. Было выделено семь фаз: парамагнитная, спиральная, лабиринтная, спираль-скирмионная, скирмионная, скирмион-ферромагнитная и ферромагнитная фазы, подробный анализ конфигураций которых приводится в статье.

Построены две фазовые диаграммы: на первой показано поведение системы при постоянном $D$ в зависимости от величин внешнего магнитного поля и температуры: $(T, B)$, на второй — изменение кон- фигураций системы при постоянной температуре $T$ в зависимости от величины взаимодействия Дзялошинского–Мории и внешнего магнитного поля: $(D, B)$.

Полученные в ходе численных экспериментов данные будут использованы в дальнейших исследованиях при определении модельных параметров системы для формирования стабильного скирмионного состояния и разработки методов контроля скирмионов в магнитной пленке.

Magnetic systems, in which due to competition between the direct Heisenberg exchange and the Dzyaloshinskii –Moriya interaction, magnetic vortex structures — skyrmions appear, were studied using the Metropolis algorithm.

The conditions for the nucleation and stable existence of magnetic skyrmions in two-dimensional magnetic films in the frame of the classical Heisenberg model were considered in the article. A thermal stability of skyrmions in a magnetic film was studied. The processes of the formation of various states in the system at different values of external magnetic fields were considered, various phases into which the Heisenberg spin system passes were recognized. The authors identified seven phases: paramagnetic, spiral, labyrinth, spiralskyrmion, skyrmion, skyrmion-ferromagnetic and ferromagnetic phases, a detailed analysis of the configurations is given in the article.

Two phase diagrams were plotted: the first diagram shows the behavior of the system at a constant $D$ depending on the values of the external magnetic field and temperature $(T, B)$, the second one shows the change of the system configurations at a constant temperature $T$ depending on the magnitude of the Dzyaloshinskii – Moriya interaction and external magnetic field: $(D, B)$.

The data from these numerical experiments will be used in further studies to determine the model parameters of the system for the formation of a stable skyrmion state and to develop methods for controlling skyrmions in a magnetic film.

Белки главного комплекса гистосовместимости (ГКГС) играют ключевую роль в работе адаптивной иммунной системы, и определение связывающихся с ними пептидов — важный шаг в разработке вакцин и понимании механизмов аутоиммунных заболеваний. На сегодняшний день существует ряд методов для предсказания связывания определенной аллели ГКГС с пептидом. Одним из лучших таких методов является NetMHCpan-4.0, основанный на ансамбле искусственных нейронных сетей. В данной работе представлена методология качественного улучшения архитектуры нейронной сети, лежащей в основе NetMHCpan-4.0. Предлагаемый метод использует технику построения ансамбля и добавляет в качестве входных данных оценку модели Поттса, взятой из статистической механики и являющейся обобщением модели Изинга. В общем случае модельо тражает взаимодействие спинов в кристаллической решетке. Применительно к задаче белок-пептидного взаимодействия вместо спинов используются типы аминокислот, находящихся в кармане связывания. В предлагаемом методе модель Поттса используется для более всестороннего представления физической природы взаимодействия полипептидных цепей, входящих в состав комплекса. Для оценки взаимодействия комплекса «ГКГС + пептид» нами используется двумерная модель Поттса с 20 состояниями (соответствующими основным аминокислотам). Решая обратную задачу с использованием данных об экспериментально подтвержденных взаимодействующих парах, мы получаем значения параметров модели Поттса, которые затем применяем для оценки новой пары «ГКГС + пептид», и дополняем этим значением входные данные нейронной сети. Такой подход, в сочетании с техникой построения ансамбля, позволяет улучшитьт очность предсказания, по метрике положительной прогностической значимости (PPV), по сравнению с базовой моделью.

The proteins of the Major Histocompatibility Complex (MHC) play a key role in the functioning of the adaptive immune system, and the identification of peptides that bind to them is an important step in the development of vaccines and understanding the mechanisms of autoimmune diseases. Today, there are a number of methods for predicting the binding of a particular MHC allele to a peptide. One of the best such methods is NetMHCpan-4.0, which is based on an ensemble of artificial neural networks. This paper presents a methodology for qualitatively improving the underlying neural network underlying NetMHCpan-4.0. The proposed method uses the ensemble construction technique and adds as input an estimate of the Potts model taken from static mechanics, which is a generalization of the Ising model. In the general case, the model reflects the interaction of spins in the crystal lattice. Within the framework of the proposed method, the model is used to better represent the physical nature of the interaction of proteins included in the complex. To assess the interaction of the MHC + peptide complex, we use a two-dimensional Potts model with 20 states (corresponding to basic amino acids). Solving the inverse problem using data on experimentally confirmed interacting pairs, we obtain the values of the parameters of the Potts model, which we then use to evaluate a new pair of MHC + peptide, and supplement this value with the input data of the neural network. This approach, combined with the ensemble construction technique, allows for improved prediction accuracy, in terms of the positive predictive value (PPV) metric, compared to the baseline model.

Предлагается простая модель на основе уравнения кинетического типа для описания распространения вируса в пространстве посредством миграции носителей вируса из выделенного центра. Рассматриваются страны, для которых применима одномерная модель: Россия, Италия, Чили. Одномерный подход возможен из-за географического расположения этих стран и их протяженности в направлениях от центров заражения (Москвы, Ломбардии и Сантьяго соответственно). Определяется изменение плотности зараженных во времени и пространстве. Применяется двухпараметрическая модель. Первый параметр — величина средней скорости распространения, соответствующий переносу инфицированных транспортными средствами. Второй параметр — частота уменьшения количества инфицированных элементов по мере продвижения по территории страны, что связано с прибытием пассажиров в места назначения, а также с карантинными мерами, препятствующими их перемещению по стране. Параметры модели определяются по фактически известным данным. Строится аналитическое решение, для получения серии расчетов применяются также простые численные методы. В модели рассматривается пространственное распространение заболевания, при этом заражения на местах не учитываются. Поэтому вычисленные значения на начальном этапе хорошо соответствуют экспериментальным данным, а затем плотность заболевших начинает быстрее возрастать из-за заражений на местах. Тем не менее модельные расчеты позволяют делать некоторые предсказания. Помимо скорости заражения, возможна аналогичная «скорость выздоровления». По моменту времени достижения охвата большей части населения страны при движении фронта выздоровления делается вывод о начале глобального выздоровления, что соответствует реальным данным.

A simple model based on a kinetic-type equation is proposed to describe the spread of a virus in space through the migration of virus carriers from a certain center. The consideration is carried out on the example of three countries for which such a one-dimensional model is applicable: Russia, Italy and Chile. The geographical location of these countries and their elongation in the direction from the centers of infection (Moscow, Milan and Lombardia in general, as well as Santiago, respectively) makes it possible to use such an approximation. The aim is to determine the dynamic density of the infected in time and space. The model is two-parameter. The first parameter is the value of the average spreading rate associated with the transfer of infected moving by transport vehicles. The second parameter is the frequency of the decrease of the infected as they move through the country, which is associated with the passengers reaching their destination, as well as with quarantine measures. The parameters are determined from the actual known data for the first days of the spatial spread of the epidemic. An analytical solution is being built; simple numerical methods are also used to obtain a series of calculations. The geographical spread of the disease is a factor taken into account in the model, the second important factor is that contact infection in the field is not taken into account. Therefore, the comparison of the calculated values with the actual data in the initial period of infection coincides with the real data, then these data become higher than the model data. Those no less model calculations allow us to make some predictions. In addition to the speed of infection, a similar “speed of recovery” is possible. When such a speed is found for the majority of the country's population, a conclusion is made about the beginning of a global recovery, which coincides with real data.

Целью работы является исследование монотонной конечно-разностной схемы второго порядка точности, созданной на основе обобщения одномерной Z-схемы. Исследование проведено для модельных уравнений переноса несжимаемой среды. В работе описано двумерное обобщение Z-схемы с нелинейной коррекцией, использующей вместо потоков косые разности, содержащие значения из разных временных слоев. Численно проверена монотонность полученной нелинейной схемы для функций-ограничителей двух видов, как для гладких решений, так и для негладких, и получены численные оценки порядка точности построенной схемы. Построенная схема является абсолютно устойчивой, но теряет свойство монотонности при превышении шага Куранта. Отличительной особенностью предложенной конечно-разностной схемы является минимальность ее шаблона.

Построенная численная схема предназначена для моделей плазменных неустойчивостей различных масштабов в низкоширотной ионосферной плазме Земли. Одна из реальных задач, при решении которых возникают подобные уравнения, — это численное моделирование сильно нестационарных среднемасштабных процессов в земной ионосфере в условиях возникновения неустойчивости Рэлея – Тейлора и плазменных структур с меньшими масштабами, механизмами генерации которых являются неустойчивости других типов, что приводит к явлению F-рассеяния. Вследствие того, что процессы переноса в ионосферной плазме контролируются магнитным полем, в поперечном к магнитному полю направле- нии предполагается выполнение условия несжимаемости плазмы.

The aim of the work is to study a monotone finite-difference scheme of the second order of accuracy, created on the basis of a generalization of the one-dimensional Z-scheme. The study was carried out for model equations of the transfer of an incompressible medium. The paper describes a two-dimensional generalization of the Z-scheme with nonlinear correction, using instead of streams oblique differences containing values from different time layers. The monotonicity of the obtained nonlinear scheme is verified numerically for the limit functions of two types, both for smooth solutions and for nonsmooth solutions, and numerical estimates of the order of accuracy of the constructed scheme are obtained.

The constructed scheme is absolutely stable, but it loses the property of monotony when the Courant step is exceeded. A distinctive feature of the proposed finite-difference scheme is the minimality of its template. The constructed numerical scheme is intended for models of plasma instabilities of various scales in the low-latitude ionospheric plasma of the Earth. One of the real problems in the solution of which such equations arise is the numerical simulation of highly nonstationary medium-scale processes in the earth’s ionosphere under conditions of the appearance of the Rayleigh – Taylor instability and plasma structures with smaller scales, the generation mechanisms of which are instabilities of other types, which leads to the phenomenon F-scattering. Due to the fact that the transfer processes in the ionospheric plasma are controlled by the magnetic field, it is assumed that the plasma incompressibility condition is fulfilled in the direction transverse to the magnetic field.

В работе исследуется влияние быстрого локального выделения тепла вблизи обтекаемой сверхзвуковым потоком газа (воздуха) поверхности на область отрыва, возникающую при быстром его повороте. Данная поверхность состоит из двух плоскостей, образующих при пересечении тупой угол, так что при обтекании этой поверхности сверхзвуковой поток газа поворачивается на положительный угол, что формирует косой скачок уплотнения, взаимодействующий с пограничным слоем и вызывающий отрыв потока. Быстрый локальный нагрев газа над обтекаемой поверхностью моделирует протяженный искровой разряд субмикросекундной длительности, пересекающий поток. Газ, нагретый в зоне разряда, взаимодействует с областью отрыва. Течение можно считать плоским, поэтому численное моделирование проводится в двумерной постановке. Численное моделирование проведено для ламинарного режима течения с использованием солвера sonicFoam пакета программ OpenFOAM.

В работе описан способ построения двумерной расчетной сетки с использованием шестигранных ячеек. Выполнено исследование сеточной сходимости. Приводится методика задания начальных профилей параметров течения на входе в расчетную область, позволяющая сократить время счета при уменьшении количества расчетных ячеек. Описан способ нестационарного моделирования процесса быстрого локального нагрева газа, заключающегося в наложении дополнительных полей повышенных значений давления и температуры, вычисленных из величины энергии, вложенной в набегающий сверхзвуковой поток газа, на соответствующие поля величин, предварительно полученные в стационарном случае. Параметры энерговклада в поток, соответствующие параметрам процесса инициирования электрического разряда, а также параметры набегающего потока близки к экспериментальным величинам.

При анализе данных численного моделирования получено, что быстрый локальный нагрев приводит к возникновению газодинамического возмущения (квазицилиндрической ударной волны и нестационарного завихренного течения), которое при взаимодействии с областью отрыва приводит к смещению точки отрыва вниз по потоку. В работе рассмотрен вопрос о влиянии энергии, затраченной на локальный нагрев газа, и положения места нагрева относительно точки отрыва на величину максимального ее смещения.

The influence of the process of initiating a rapid local heat release near surface streamlined by supersonic gas (air) flow on the separation region that occurs during a fast turn of the flow was investigated. This surface consists of two planes that form obtuse angle when crossing, so that when flowing around the formed surface, the supersonic gas flow turns by a positive angle, which forms an oblique shock wave that interacts with the boundary layer and causes flow separation. Rapid local heating of the gas above the streamlined surface simulates long spark discharge of submicrosecond duration that crosses the flow. The gas heated in the discharge zone interacts with the separation region. The flow can be considered two-dimensional, so the numerical simulation is carried out in a two-dimensional formulation. Numerical simulation was carried out for laminar regime of flow using the sonicFoam solver of the OpenFOAM software package.

The paper describes a method for constructing a two-dimensional computational grid using hexagonal cells. A study of grid convergence has been carried out. A technique is given for setting the initial profiles of the flow parameters at the entrance to the computational domain, which makes it possible to reduce the computation time by reducing the number of computational cells. A method for non-stationary simulation of the process of rapid local heating of a gas is described, which consists in superimposing additional fields of increased pressure and temperature values calculated from the amount of energy deposited in oncoming supersonic gas flow on the corresponding fields of values obtained in the stationary case. The parameters of the energy input into the flow corresponding to the parameters of the electric discharge process, as well as the parameters of the oncoming flow, are close to the experimental values.

During analyzing numerical simulation data it was found that the initiation of rapid local heating leads to the appearance of a gas-dynamic perturbation (a quasi-cylindrical shock wave and an unsteady swirling flow), which, when interacting with the separation region, leads to a displacement of the separation point downstream. The paper considers the question of the influence of the energy spent on local heating of the gas, and of the position on the streamlined surface of the place of heating relative to the separation point, on the value of its maximum displacement.

Проведенное исследование основано на сочетании трехмерной гидродинамической модели глобального климата, включая модель океана с реальными глубинами и конфигурацией континентов, модель эволюции морского льда и энерго-, влагобалансовую модель атмосферы. Концентрация аэрозоля от 2010 г. до 2100 г. рассчитывается как управляющий параметр для стабилизации среднегодовой температуры воздуха у поверхности земли. На основе расчетов предполагается, что выбросы серы от 2010 г. до 2100 г. изменяются линейно для первого сценария и квадратично — для второго роста СО₂. Граница Парето исследована и визуализирована для двух параметров — среднеквадратичного отклонения атмосферной температуры для зимнего и летнего сезонов.

The study is based on a three-dimensional hydrodynamic global climate coupled model, including ocean model with real depths and continents configuration, sea ice evolution model and energy and moisture balance atmosphere model. Aerosol concentration from the year 2010 to 2100 is calculated as a controlling parameter to stabilize mean year surface air temperature. It is shown that by this way it is impossible to achieve the space and seasonal uniform approximation to the existing climate, although it is possible significantly reduce the greenhouse warming effect. Climate will be colder at 0.1–0.2 degrees in the low and mid-latitudes and at high latitudes it will be warmer at 0.2–1.2 degrees. The Pareto frontier is investigated and visualized for two parameters — atmospheric temperature mean square deviation for the winter and summer seasons. The Pareto optimal amount of sulfur emissions would be between 23.5 and 26.5 TgS/year.