Statistical relationships between variations of the geomagnetic field, auroral electrojet, and geomagnetically induced currents

Using observations from the IMAGE magnetic observatories and the station for recording geomagnetically induced currents (GIC) in the electric transmission line in 2015, we examine relationships between geomagnetic field and GIC variations. The GIC intensity is highly correlated (R>0.7) with the field variability |dB/dt| and closely correlated with variations in the time derivatives of X and Y components. Daily variations in the mean geomagnetic field variability |dB/dt| and GIC intensity have a wide night maximum, associated with the electrojet, and a wide morning maximum, presumably caused by intense Pc5–Pi3 geomagnetic pulsations. We have constructed a regression linear model to estimate GIC from the time derivative of the geomagnetic field and AE index. Statistical distributions of the probability density of the AE index, geomagnetic field derivative, and GIC correspond to the log-normal law. The constructed distributions are used to evaluate the probabilities of extreme values of GIC and |dB/dt|

Investigation of the influence of the nucleon pairing force on the octupole deformation of uranium and thorium isotopes in the mean-field approximation

In the Hartree-Fock-Bogolyubov approximation and the axial symmetry of nuclei with Skyrm forces SkM∗, the properties of U with A = 220 – 232 and Th with A = 218 – 230 isotopes were calculated. In addition, the properties of U and Th near the neutron drip line A = 280 – 292 were calculated. Pairing of nucleons in nuclei is described by pairing forces of zero radius of action of a mixed type with different sets of pairing force constants. In our calculations, we used constrained conditions on the parameters of quadrupole β₂ and octupole β₃ deformations of nuclei. It is shown that for the considered U isotopes, as well as Th isotopes, the deformation of the β₃ nuclei strongly depends on the choice of the parameters of the nucleon pairing force.У наближенні Хартрі-Фока-Боголюбова та аксіальної симетрії ядер із силами Скірма SkM∗ проведено розрахунки властивостей ізотопів U з A = 220 – 232 і Th з A = 218 – 230. Додатково проведені розрахунки властивостей U і Th біля границі нейтронної стабільності A = 280 – 292. Спарювання нуклонів у ядрах описується силами спарювання нульового радіуса дії змішаного типу з різними наборами констант сил спарювання. У розрахунках використовувалися накладені умови на параметри квадрупольної β₂ і октупольної β₃ деформацій ядер. Показано, що для розглянутих ізотопів U, як і ізотопів Th, деформація ядер β₃ сильно залежить від вибору параметрів сили спарювання нуклонів

Role of dislocations in formation of ohmic contacts to heavily doped n-Si

We present experimental results concerning a high density of structural defects (in particular, dislocations) in the near-contact region of heavily doped n-silicon. They appear in the course of firing Au Pd Ti Pd -Si n ohmic contact at 450С for 10 min in a vacuum of ~10 Pa⁻⁴ . These defects lead to appearance of metal shunts that determine the current flow mechanism in these ohmic contacts. The calculated and experimental temperature dependences of contact resistivity, ρс(Т), are in good agreement. It is shown that ρс increases with temperature. This is characteristic of a model of ohmic contacts with a high dislocation density in the near-contact region of semiconducto

On a feature of temperature dependence of contact resistivity for ohmic contacts to n-Si with an n⁺ -n doping step

We present both theoretical and experimental temperature dependences of contact resistivity ρс(Т) for ohmic contacts to the silicon n⁺ -n-structures whose n⁺ -layer was formed using phosphorus diffusion or ion implantation. The ρс(Т) dependence was measured in the 125–375 K temperature range with the transmission line method, with allowance made for conduction in both the n⁺ -layer and n⁺ -n doping step

Тривимірне моделювання часового поля променевим і кінцево-різницевим методами для вирішення завдань сейсмології

Three-dimensional modeling of temporal field allows reconstructing the kinamatics of wave processes observed in the Earth during seismic studies and determine in such a way spatial position into it of the studied objects. Modeling is also an important tool for inspection of correctness of methods for solving the inverse problem during the process of interpretation both seismological and seismic exploration data. Seismic rays which determine the direction of the flow of high-frequency part of seismic wave field energy are very important in seismology.  Tracking the rays and their calculation allows solving different problems of seismology as well as checking the accuracy of results obtained while different methods of processing and interpretation of data observed on the Earth surface are used.A version of finite-differential modeling of temporal field in three-dimensional spherical Earth is based on direct net approximation of the eikonal equation and it is this approach to reconstruction of values of seismic waves arrivals to any point of Earth depths is the most stabile and as a result can guarantee correctness of solving lots of applied problems of seismology. At the same time continuation of temporal field is a part of calculative process of finite-differential migration designed at the Institute of Geophysics named after S.I.Subbotin of NAS of Ukraine and is used while processing seismic exploration data.Development of computer technologies during recent decades brought the appearance of theories, algorithms, and software-based complexes realizing the solving of three-dimensional geophysical in particular seismological problems.  The paper gives theoretical foundations, algorithms, and results of application  of three-dimensional temporal fields modeling developed for both radial and finite-differential methods by practical examples.Трехмерное моделирование временного поля позволяет воспроизвести кинематику волновых процессов, наблюдаемых на Земле при сейсмических исследованиях, и, таким образом, определить пространственное размещение в ней целевых объектов. Также моделирования является важным инструментом для проверки корректности методов решения обратных задач в процессе интерпретации как сейсмологических, так и сейсморазведочных данных. В сейсмологии особое значение имеют сейсмические лучи, которые определяют направление потока энергии высокочастотной части сейсмического волнового поля. Проследить лучей и их расчет позволяют решать широкий круг задач сейсмологии, а также проверять точность полученных результатов при применении различных методов обработки и интерпретации наблюдений на земной поверхности данных. Вариант скинченноризницевого моделирования временного поля в трехмерной сферической Земли основывается на непосредственном сеточной аппроксимации уравнения эйконала, и именно этот подход к определению значений времени прихода сейсмических волн в каждую точку земных глубин является устойчивым, а следовательно, может гарантировать корректность решения многих прикладных задач сейсмологии. В то же время продолжение временного поля является частью вычислительного процесса скинченноризницевои миграции, которая разработана в Институте геофизики им. С. И. Субботина НАН Украины и которую используют при обработке данных сейсморазведки. В течение последних десятилетий с развитием компьютерных технологий связано появление теорий, алгоритмов и программных комплексов, которые позволили реализовать решения трехмерных геофизических, в частности сейсмологических, задач. В статье приведены теоретические основы, алгоритмы и результаты применения разработанных вариантов трехмерного моделирования временных полей как лучевой, так и скинченноризницевим методом на практических примерах.Тривимірне моделювання часового поля дає змогу відтворити кінематику хвильових процесів, які спостерігаються в Землі при сейсмічних дослідженнях, і, таким чином, визначити просторове розміщення в ній цільових об'єктів. Також моделювання є важливим інструментом для перевірки коректності методів розв'язання обернених задач у процесі інтерпретацій як сейсмологічних , так і сейсморозвідувальних даних. У її особливе значення мають сейсмічні промені, які визнача­ють напрямок потоку енергії високочастотної частини сейсмічного хвильового поля. Простеження променів та їх розрахунок дають змогу розв'язувати широке коло завдань сейсмології, а також перевіряти точність отриманих результатів при застосуванні різних методів обробки та інтерпретації спостережених на земній поверхні даних. Варіант скінченно-різницевого моделювання часового поля у тривимірний сферичній Землі ґрунтується  на безпосередній сітковий апроксимації рівняння ейконалу, та саме цей підхід до визначення значень часу приходу сейсмічних хвиль у кожну точку земних глибин є найстійкішим, а отже, може гарантувати коректність розв'язання багатьох прикладних завдань її. Водночас продовження часового поля е частиною обчислювального процесу скінченно-різницевої міграцій, яка розроблена в Інституті геофізики ім. С.І. Субботіна НАН України і яку використовують при обробці даних сейсморозвідки. Протягом останніх десятиліть з розвитком комп'ютерних технологій пов'язана поява теорій, алгоритмів та програмних комплексів, які дали змо­гу реалізувати розв'язання тривимірних геофізичних, зокрема сейсмологічних, за­дач. У статті наведено теоретичні основи, алгоритми та результати застосування роз роблених варіантів тривимірного моделювання часових полів як променевим, так і скінченно-різницевим методом на практичних прикладах

Тривимірне скінченно-різницеве моделювання хвильового поля з урахуванням сферичності Землі

Three dimensional modeling of the wave field makes possible to obtain a concept on distribution of seismic waves in geological medium with different complexity of its structure. The location of the source of fluctuations may be defined in this case at any depth. Therefore, three dimensional modeling of the wave field may be used while solving the problems of seismology related to the nature of strains in the Earth’s interior, while studying the dangerous phenomena of endogenous origin as well as in the case of studies of the deep structure of the Earth, the boundaries of its division and its tectonic structures.All existing methods of numerical simulation of the wave field are based on the solving the equation of elastodynamics with the help of different methods defined by interfacial and initial conditions and are subdivided into axial-symmetrical and pseudo-spectral. The first ones are based on approximation of structural model as a rotary-symmetrical along the vertical axis of the source of fluctuations and solving of the equation of elastodynamics in sylindric or spheric coordinates and in the second ones while solving the equation of elastodynamics Fourier transformation is used.At the Institute of Geophysics named after S.I. Subbotin NAS of Ukraine a method of finite-differential simulation of the wave field has been elaborated both in two-dimensional and three-dimensional variants. The choice of finite-differential method for proceeding of the wave field is explained by its high stability and accuracy as well as with possibility of its application to the wide class of heterogenous models of the medium. Finite-differential modeling is based on the solving of the scalar equation with application of the spacio-temporal net. While solving the problems of seismology in the case of three-dimensional modeling of the wave field we are to take sphericity of the Earth into account. For this purpose in the differential wave equation a transition from Cartesian to spherical coordinate system has been provided. The theory and algorithm of the elaborated three-dimensional modeling of the wave field with finite differential method taking into account sphericity of the Earth has been considered in the paper with demonstration on the model example.Трехмерное моделирование волнового поля дает возможность получить представление о распространении сейсмических волн в геологической среде разной степени сложности строения. При этом расположение источника колебаний может быть задано на любой глубине. Таким образом, трехмерное моделирование волнового поля можно применять при решении задач сейсмологии, связанных с природой напряжений внутри Земли, при исследованиях опасных явлений ендогенного происхождения, а также при изучении глибинного строения Земли, границ ее раздела и тектонических структур. Все известные методы численного моделирования волнового поля основаны на решении уравнения эластодинамики разными методами с заданными граничными и начальными условиями и подразделяются на осесимметричные и псевдоспектральные. Первые основаны на аппроксимации структурной модели как вращательно-симметричной вдоль вертикальной оси источника колебаний и решении уравнения эластодинамики в цилиндрических или сферических координатах, а вторые - на решении уравнения эластодинамики с использованием преобразования Фурье. В Институте геофизики им. С.И. Субботина НАН Украины разработан метод конечно-разностного моделирования волнового поля как в двумерном, так и трехмерном варианте. Выбор конечно-разностного метода для продолжения волнового поля объясняется его высокой устойчивостью и точностью, а также возможностью использовать при его применении широкий класс неоднородных моделей среды. Конечно-разностное моделирование основано на решении скалярного волнового уравнения с применением пространственно-временной сетки. При решении задач сейсмологии в случае трехмерного моделирования волнового поля необходимо учитывать сферичность Земли. Для этого в дифференциальном волновом уравнении предусмотрен переход с декартовой на сферическую систему координат. В статье рассмотрены теория и алгоритм разработанного трехмерного моделирования волнового поля конечно-разностным методом с учетом сферичности Земли. Эффективность метода показана на модельном примере.Тривимірне моделювання хвильового поля дає змогу отримати уявлення про поширення сейсмічних хвиль у геологічному середовищі різного ступеня складності будови. При цьому розташування джерела коливань може бути задано на будь-якій глибині. Таким чином, тривимірне моделювання хвильового поля можна застосовувати при розв’язанні задач сейсмології, пов’язаних з природою напружень всередині Землі, під час досліджень небезпечних явищ ендогенного походження, а також вивчення глибинної будови Землі, меж її поділу і тектонічних структур. Всі існуючі методи чисельного моделювання хвильового поля засновані на розв’язанні рівняння еластодинаміки різними методами із заданими граничними і початковими умовами. Їх поділяють на осесиметричні і псевдоспектральні. Перше засноване на апроксимації структурної моделі як обертально-симетричної уздовж вертикальної осі джерела коливань і розв’язанні рівняння еластодінамікі в циліндричних або сферичних координатах, а друге — на розв’язанні рівняння еластодинаміки з використанням перетворення Фур’є. В Інституті геофізики ім. С.І. Суботіна НАН України розроблено метод скінченно-різницевого моделювання хвильового поля як у двовимірному, так і тривимірному варіанті. Вибір скінченно-різницевого методу для продовження хвильового поля пояснюється його високою стійкістю і точністю, а також можливістю використовувати при його застосуванні широкий клас неоднорідних моделей середовища. Скінченно-різницеве моделювання ґрунтується на розв’язанні скалярного хвильового рівняння із застосуванням просторово-часової сітки. При розв’язанні задач сейсмології в разі тривимірного моделювання хвильового поля необхідно враховувати сферичність Землі. Для цього в диференціальному хвильовому рівнянні передбачено перехід з декартової на сферичну систему координат. У статті розглянуто теорію і алгоритм розробленого тривимірного моделювання хвильового поля скінченно-різницевим методом з урахуванням сферичності Землі, ефективність якого показана на модельному прикладі

Finite-difference wave migration of initial seismograms of a common point of explosion in the time domain

The final procedure of the standard production processing and interpretation of seismic data is the post-stack migration by the central deep point method (CDP). However, in practice, a more informative image of the deep section is given by the pre-stack migration according to initial seismograms. The actual task of processing and interpretation of the observed wave field to produce the image of the deep structure without the procedure of the CDP stack is presented. This solution is implemented using a new version of the finite-difference pre-stack migration of the shot record in time domain. The advantages and possibilities of applications of this method are shown in model and real examples

МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ И АНАЛИЗ УСТОЙЧИВОСТИ НИЗКОЧАСТОТНЫХ ПРОЦЕССОВ В МАРШЕВОМ ЖРД С ДОЖИГАНИЕМ ГЕНЕРАТОРНОГО ГАЗА

Possible reasons for the loss of stability of low-frequency (up to 50 Hz) dynamic processes in staged liquid-propellant rocket engines (LRE) are considered. It is shown that the theoretical determination of the LRE stability and clarify the nature of its low-frequency oscillations is necessary to investigate the midtilinked dynamic system «feed lines - LRE». An approach to the solution of this problem is presented in linear formulation. The approach is based on the calculation of the spectrum of the matrix, which is an operator of the linear dynamic system «feed lines - LRE» and the decomposition of this system. The mathematical modeling of the low- frequency dynamics of the liquid rocket propulsion system inclusive the perspective staged LRE in an oxidizer rich combustion cycle. The investigated liquid rocket engine is stable. It is shown that taking into account cavitation phenomena in the LRE pumps led to a significant decrease in natural frequencies and decrements offluid oscillations in the LRE feed lines. The low-frequency dynamic processes in the flow regulator are characterized by large decrements of oscillations. The dynamical system under investigation is sensitive to the value of the gas time delays in the gas generator.Рассмотрены возможные причины потери устойчивости низкочастотных (до 50Гц) динамических процессов в жидкостных ракетных двигателях (ЖРД) с дожиганием генераторного газа. Показано, что для теоретического определения устойчивости ЖРД и выяснения природы его низкочастотных колебаний необходимо исследовать многосвязную динамическую систему «питающие магистрали — ЖРД». Изложен подход к решению этой задачи в линейной постановке, основанный на расчете спектра матрицы, представляющей собой оператор линейной динамической системы «питающие магистрали — ЖРД» и декомпозиции этой системы. Выполнено математическое моделирование низкочастотной динамики жидкостной ракетной двигательной установки (ЖРДУ), включающей перспективный маршевый ЖРД с дожиганием окислительного генераторного газа. Установлено, что исследуемая ЖРДУ устойчива. Показано, что учет кавитационных явлений в насосах ЖРД привел к существенному уменьшению собственных частот и декрементов колебаний жидкости в линиях питания окислителем и горючим. Определено, что низкочастотные динамические процессы в регуляторе расхода характеризуются относительно большими декрементами колебаний. Установлено, что исследуемая динамическая система чувствительна к значению времени пребывания газа в газогенераторе