Все выпуски

Работа посвящена вопросу об абсолютной непрерывности спектра двумерного обобщенного периодического оператора Шрёдингера $H_g+V=-\nabla g\nabla+V$, где непрерывная положительная функция $g$ и скалярный потенциал $V$ имеют общую решетку периодов $Λ$. Решения уравнения $(H_g+V)\varphi=0$ определяют, в частности, электрическое и магнитное поля для электромагнитных волн, распространяющихся в двумерных фотонных кристаллах. При этом функция $g$ и скалярный потенциал $V$ выражаются через диэлектрическую проницаемость $\varepsilon$ и магнитную проницаемость $\mu$ ($V$ также зависит от частоты электромагнитной волны). Диэлектрическая проницаемость $\varepsilon$ может быть разрывной функцией (и обычно выбирается кусочно-постоянной), поэтому возникает задача об ослаблении известных условий гладкости для функции $g$, обеспечивающих абсолютную непрерывность спектра оператора $H_g+V$. В настоящей работе предполагается, что коэффициенты Фурье функций $g^{\pm\frac12}$ при некотором $q\in[1, \frac43)$ удовлетворяют условию $\sum\left(|N|^\frac12\left|\left(g^{\pm\frac12}\right)_N\right|\right)^q<+\infty$ и скалярный потенциал $V$ имеет нулевую грань относительно оператора $-Δ$ в смысле квадратичных форм. Пусть $K$ - элементарная ячейка решетки $Λ$, $K^*$ - элементарная ячейка обратной решетки $\Lambda^*$. Оператор $H_g+V$ унитарно эквивалентен прямому интегралу операторов $H_g(k)+V$, где $k$ - квазиимпульс из $2\pi K^*$, действующих в $L^2(K)$. Последние операторы можно также рассматривать при комплексных векторах $k+ik'\in \mathbb{C}^2$. В статье используется метод Томаса. Доказательство абсолютной непрерывности спектра оператора $H_g+V$ сводится к доказательству обратимости операторов $H_g(k+ik')+V-\lambda$, $\lambda\in \mathbb{R}$, при определенным образом выбираемых комплексных векторах $k+ik'\in \mathbb{C}^2$ (зависящих от $g$, $V$ и числа $\lambda$) с достаточно большой мнимой частью $k'$.

The paper is concerned with the problem of absolute continuity of the spectrum of the two-dimensional generalized periodic Schrodinger operator $H_g+V=-\nabla g\nabla+V$ where the continuous positive function $g$ and the scalar potential $V$ have a common period lattice $\Lambda$. The solutions of the equation $(H_g+V)\varphi=0$ determine, in particular, the electric field and the magnetic field of electromagnetic waves propagating in two-dimensional photonic crystals. The function $g$ and the scalar potential $V$ are expressed in terms of the electric permittivity $\varepsilon$ and the magnetic permeability $\mu$ ($V$ also depends on the frequency of the electromagnetic wave). The electric permittivity $\varepsilon$ may be a discontinuous function (and usually it is chosen to be piecewise constant) so the problem to relax the known smoothness conditions on the function $g$ that provide absolute continuity of the spectrum of the operator $H_g+V$ arises. In the present paper we assume that the Fourier coefficients of the functions $g^{\pm\frac12}$ for some $q\in[1, \frac43)$ satisfy the condition $\sum\left(|N|^\frac12\left|\left(g^{\pm\frac12}\right)_N\right|\right)^q<+\infty$, and the scalar potential $V$ has relative bound zero with respect to the operator $-\Delta$ in the sense of quadratic forms. Let $K$ be the fundamental domain of the lattice $\Lambda$, and assume that $K^*$ is the fundamental domain of the reciprocal lattice $\Lambda^*$. The operator $H_g+V$ is unitarily equivalent to the direct integral of operators $H_g(k)+V$, with quasimomenta $k\in 2\pi K^*$, acting on the space $L^2(K)$. The last operators can be also considered for complex vectors $k+ik'\in \mathbb{C}^2$. We use the Thomas method. The proof of absolute continuity of the spectrum of the operator $H_g+V$ amounts to showing that the operators $H_g(k+ik')+V-\lambda$, $\lambda\in \mathbb{R}$, are invertible for some appropriately chosen complex vectors $k+ik'\in \mathbb{C}^2$ (depending on $g$, $V$, and the number $\lambda$) with sufficiently large imaginary parts $k'$.

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

The formation process of a spherical shock impulse in gas and its interaction with a protective aqueous foam barrier, accompanied by formation of vortex flows, are numerically investigated. The problem is solved in a two-dimensional axisymmetric formulation using a two-phase model of a gas-liquid mixture, which includes the laws of conservation of mass, momentum and energy of the mixture and an equation for the dynamics of volume content of phases.The numerical implementation is carried out on the basis of the OpenFOAM package using the standard compressibleMultiphaseInterFoam solver, modified in accordance with the conditions of the problem and model representations. The discretization of the system of equations in the chosen solver is carried out by the method of finite volumes using the computational Pimple algorithm. A significant decrease in the intensity of the shock wave in its interaction with the aqueous foam barrier is shown and the causes leading to vortex formation in the gas region are revealed. The reliability of the results obtained is estimated by comparison with solutions of a similar problem by other numerical methods.

Рассмотрено волновое уравнение с двумя пространственными и одной временной независимыми переменными и эффектом наследственности вида $$\frac{\partial^2 u}{\partial t^2}=a^2\left(\frac{\partial^2 u}{\partial x^2} + \frac{\partial^2 u}{\partial y^2}\right) + f\big(x,y,t,u(x,y,t),u_t(x,y,\cdot)\big),\\u_t(x,y,\cdot)=\left\{u(x,y,t+\xi),-\tau \leqslant \xi\leqslant 0\right\}. $$На основе идеи разделения текущего состояния и функции-предыстории сконструировано семейство сеточных методов для численного решения этого уравнения. По текущему состоянию строится полный аналог известного для уравнения без запаздывания метода с факторизацией, а влияние предыстории учитывается с помощью интерполяционных конструкций. Исследован порядок локальной погрешности алгоритма. Получена теорема о сходимости и порядке сходимости методов с помощью вложения в общую разностную схему систем с последействием. Приводятся результаты расчетов тестового примера с переменным запаздыванием.

The paper presents the consideration of the wave equation with two space variables and one time variable and with heredity effect  $$\frac{\partial^2 u}{\partial t^2}=a^2\left(\frac{\partial^2 u}{\partial x^2} + \frac{\partial^2 u}{\partial y^2}\right) + f\big(x,y,t,u(x,y,t),u_t(x,y,\cdot)\big),\\u_t(x,y,\cdot)=\left\{u(x,y,t+\xi),-\tau \leqslant \xi\leqslant 0\right\}. $$A family of grid methods is constructed for the numerical solution of this equation; the methods are based on the idea of separating the current state and the history function. A complete analog of the factorization method which is known for an equation without delay is constructed according to the current state. Influence of prehistory is taken into consideration by interpolation constructions. The local error order of the algorithm is investigated. A theorem on the convergence and on the order of convergence of methods is obtained by means of embedding into a general difference scheme with aftereffect. The results of calculating a test example with variable delay are presented.

Определяется параметрическое семейство конечномерных пространств специальных квадратичных сплайнов лагранжевого типа. В каждом пространстве в качестве решения начально-граничной задачи для простейшего волнового уравнения предлагается оптимальный сплайн, дающий наименьшую невязку, представляющую собой квадрат нормы в пространстве L₂. Для коэффициентов этого сплайна и для его невязки получены точные формулы. Формула для коэффициентов сплайна представляет собой линейную форму от конечных разностей дискретно заданных начальных и граничных условий исходной задачи. Формула для невязки J представляет собой положительно определенную квадратичную форму от этих же величин. Коэффициенты обеих форм вычислимы через многочлены Чебышева 2-го рода. Явный вид формулы для невязки позволяет при заданной точности вычислений ε > 0 решить неравенство J < ε² и получить априори достаточное количество узлов разностной схемы.

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

We define the parameter family of finite-dimensional spaces of special quadratic splines of Lagrange’s type. In each space, the optimal spline which gives the smallest residual being a square of the norm in the space L₂, is proposed as a solution to the initial-boundary problem for the simplest wave equation. The exact formulas for the coefficients of the spline and its residual are obtained. The formula for the coefficients of this spline is a linear form of finite differences of the discretely given initial and boundary conditions of the original problem. The formula for the residual J is a positive definite quadratic form of these quantities. The coefficients of both forms are computable via Chebyshev’s polynomials of the second kind. The explicit form of the formula for the residual allows to solve the inequality J < ε² for a given computing accuracy ε > 0 and to receive a priori sufficient number of nodes of a difference scheme.

The investigations were carried out for one time layer, which has two sublayers. We obtained difference formulas of the initial condition for the partial derivative with respect to time. They allow to create a difference scheme for the new layer, which in turn allows to continue the iterative computational process in time as far as desired.

С помощью упрощенного метода возмущений исследуется влияние взаимодействия между пузырьками на распространение волн в однородном слабосжимаемом вязкоупругом пузырьковом потоке. С использованием подхода сохранения кинетической энергии выводится уравнение динамики пузырьков. Динамика пузырьков и уравнения смеси в сочетании с уравнением состояния газа позволяют исследовать явление распространения ударной волны в смеси. Выведено двумерное уравнение Кортевега-де Фриза-Бюргера в терминах профиля давления. Установлено, что при использовании рассматриваемых нами параметров взаимодействие между пузырьками не оказывает влияния.

The effect of bubble-bubble interaction on wave propagation in homogeneous weakly compressible viscoelastic bubbly flow is investigated using the reductive perturbation method. The bubble dynamics equation is derived using the kinetic energy conservation approach. The bubble dynamics and mixture equations are coupled with the equation of state for gas to investigate the shock wave propagation phenomenon in the mixture. A two-dimensional Korteweg-de VriesBurger (KdVB) equation in terms of a pressure profile is derived. It is found that the bubble-bubble interaction has no effect when using the parameters under our consideration.

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

The solution of a boundary value problem for a simple wave equation defined on a rectangle can be represented as a sum of two terms. They are solutions of two boundary value problems: in the first case, the boundary functions are constant, while in the second the initial functions have a special form. Such decomposition allows to apply two-dimensional splines for the numerical solution of both problems. The first problem was studied previously, and an economical algorithm of its numerical solution was developed.
To solve the second problem we define a finite-dimensional space of splines of Lagrangian type, and recommend an optimal spline giving the smallest residual as a solution. We obtain exact formulas for the coefficients of this spline and its residual. The formula for the coefficients of this spline is a linear form of initial finite differences defined on the boundary.
The formula for the residual is a sum of two simple terms and two positive definite quadratic forms of new finite differences defined on the boundary. Elements of matrices of forms are expressed through Chebyshev polynomials, both matrices are invertible and have the property that their inverses matrices are of tridiagonal form. This feature allows us to obtain upper and lower bounds for the spectrum of matrices, and to show that the residual tends to zero when the numerical problem dimension increases. This fact ensures the correctness of the proposed algorithm of numerical solution of the second problem which has linear computational complexity.