Scientific journal
Fundamental research
ISSN 1812-7339
"Перечень" ВАК
ИФ РИНЦ = 1,674

MODELING OF THE ACTIVATION PROCESSES OF NATURAL WATER MIXING CONCRETE IN THE ELECTROMAGNETIC FIELD

Boyarkin D.I. 2 Fomichev V.T. 1 Erofeev V.T. 2 Emelyanov D.V. 2 Matvievskiy A.A. 3
1 Volgograd State Architectural and Construction University
2 Mordovian State University n.a. N.P. Ogarev
3 JSC «MAXMIR»
1993 KB
It is shown that activation the mortar and concrete mixes, and also the components making them is one of effective processing methods. At a certain stage of process of the formation of the cement binder it is possible to influence actively technological, structural and physics and technology properties of the received materials, namely by use for a mixing of the activated water and water solutions. The devices and devices directed on generation of electromagnetic waves have to correspond to a set of criteria that comes to light by means of the models directed on the solution of initial and regional tasks for the three-dimensional wave equation. Implementation of the differential scheme of the numerical solution of a task by means of a splitting method which along with stability possesses property of a minimum of volume of calculations is shown. For the proof of stability of differential schemes the spectral sign of Neumann is used.
activation
natural water
mixing of concrete
electromagnetic field
modeling of processes
solution of the wave equation
splitting method
stability of differential schemes
1. Erofeev V.T., Mitina E.A., Matvievskij A.A., Osipov A.K., Emelyanov D.V., Judin P.V. Kompozicionnye stroitel’nye materialy na aktivirovannoj vode zatvorenija [Composite building materials on an activated water mixing]. Stroitel’nye materialy [Building materials]. 2007. no. 11. рр. 56–57.
2. Erofeev V.T., Mitina E.A., Matvievskij A.A., Emelyanov D.V., Judin P.V. Dolgovechnost’ cementnyh kompozitov na aktivirovannoj vode [Durability of cement composites along a water]. Promyshlennoe i grazhdanskoe stroitel’stvo [Industrial and civil construction]. 2008. no. 7. рр. 51–53.
3. Erofeev V.T., Fomichev V.T., Emelyanov D.V., Rodin A.I., Eremin A.V. Vlijanie aktivirovannoj vody zatvorenija na strukturoobrazovanie cementnyh past [The effect of activated water mixing structure of cement pastes]. Vestnik Volgogradskogo gosudarstvennogo arhitekturno-stroitel’nogo universiteta [Bulletin of Volgograd state architectural and construction University]. Ser.: Construction and architecture. 2013. Vol. 30(49). рр. 179–183.
4. Sedova A.A., Osipov A.K., Emelyanov D.V., Judin P.V., Erofeev V.T. Ustanovlenie pokazatelej fiziko-himicheskih svojstv aktivirovannoj vody dlja sostavlenija matematicheskoj modeli tehnologicheskogo processa [Establishing indicators of physico-chemical properties of activated water for the compilation of mathematical process models]. Vestnik MADI [Bulletin of MADI]. 2011. no. 2. рр. 101–108.
5. Janenko N.N. Metod drobnyh shagov reshenija mnogomernyh zadach matematicheskoj fiziki [Fractional step method for solving multidimensional problems of mathematical physics]. Novosibirsk: Science, 1967.

Активация растворных и бетонных смесей, а также составляющих их компонентов является одним из действенных технологических приемов, позволяющих целенаправленно регулировать свойства изделий на их основе [1, 2].

Анализ теоретических исследований и моделирование процессов воздействия электромагнитного поля на природную воду затворения растворных и бетонных смесей на основе цементного вяжущего позволяют утверждать, что на определенном этапе процесса их структурообразования имеется возможность активно влиять на технологические, структурные и физико-технические свойства получаемого материала [3, 4].

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

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

Классическая постановка начальной задачи для трехмерного волнового уравнения в положительном полупространстве будет иметь вид

boayrkin001.wmf (1)

boayrkin002.wmf (2)

boayrkin003.wmf (3)

Используя отображение boayrkin004.wmf, где boayrkin005.wmf, можем получить постановку этой задачи в обобщенном смысле, в пространстве С.Л. Соболева. Рассмотрим выражение

boayrkin006.wmf

Тогда

boayrkin007.wmf

boayrkin008.wmf

Интегрируя которое и используя формулу Гаусса – Остроградского по произвольной области D, будем иметь

boayrkin009.wmf (4)

Если область D – все положительное полупространство (t > 0) с кусочно-гладкой границей S: boayrkin010.wmf, а boayrkin011.wmf, то получим начальную задачу для трехмерного волнового уравнения в пространстве С.Л. Соболева.

boayrkin012.wmf (5)

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

Пусть γ и γ′ – две непересекающиеся гиперповерхности в положительном полупространстве boayrkin013.wmf.

Определим функцию φ следующим образом:

boayrkin014.wmf

и функция φ гладким образом убывает от 1 до 0 вместе со своими производными boayrkin015.wmf от γ до γ′.

Рассмотрим левый интеграл из уравнения (5):

boayrkin016.wmf

Учитывая свойства функции φ, этот интеграл будет иметь вид

boayrkin017.wmf

Заменяя интегрирование по полосе от γ до γ′ интегрированием по нормали и касательным направлениям к гиперповерхности γ, получим

boayrkin018.wmf

Ширину полосы γγ′ возьмем такой, чтобы можно было предположить, что

boayrkin019.wmf boayrkin020.wmf boayrkin021.wmf boayrkin022.wmf

Так как nx = sx, ny = sy, nz = sz, nt = –st и sxds = dx, syds = dy, szds = dz, stds = dt, то получим

boayrkin023.wmf

Таким образом, принимая во внимание правый интеграл из уравнения (5), получим обобщенную постановку начальной задачи для волнового уравнения на основе интегральных законов сохранения:

boayrkin024.wmf

Реализуем явную разностную схему численного решения начально-краевой задачи для волнового уравнения.

Определим равномерную сетку с шагом h по пространственным переменным и с шагом τ по времени

boayrkin025.wmf

где xi = ih, yj = jh, zk = kh, tn = nτ, i = 0, ±1, ±2, ... , j = 0, ±1, ±2, ... , k= 0, ±1, ±2, ... , n= 0, 1, 2, ...

Для начальной задачи волнового уравнения

boayrkin026.wmf (6)

boayrkin027.wmf (7)

boayrkin028.wmf (8)

построим явную многомерную разностную схему:

boayrkin029.wmf

boayrkin030.wmf (9)

boayrkin031.wmf (10)

boayrkin032.wmf (11)

Из теории разностных схем известно, что невязка между точным решением задачи (6)–(7) и приближенным решением задачи (9)–(11) есть boayrkin033.wmf, а необходимое условие устойчивости разностной схемы (9)–(11) равняется boayrkin034.wmf.

Реализуем метод расщепления численного решения начально-краевой задачи для волнового уравнения.

При решении многомерных задач объем вычислений методом сеток существенно возрастает. Для его снижения Н.Н. Яненко [5] предложил эффективный метод расщепления, который наряду с устойчивостью обладает свойством минимальности объема вычислений. Решение задачи (6)–(7) сводится к решению трех вспомогательных задач. Построим схему метода расщепления для решения задачи (6)–(7). Введем дифференциальный оператор A, такой, что

boayrkin035.wmf

Тогда уравнение (6) можно записать в виде

boayrkin036.wmf

Оператор A является суммой трех операторов

boayrkin037.wmf boayrkin038.wmf boayrkin039.wmf

Рассмотрим три вспомогательные одномерные задачи.

Задача 1:

boayrkin040.wmf (12)

с начальными условиями

boayrkin041.wmf (13)

boayrkin042.wmf (14)

boayrkin043.wmf (15)

Задача 2:

boayrkin044.wmf (16)

с начальными условиями

boayrkin045.wmf (17)

boayrkin046.wmf (18)

Задача 3:

boayrkin047.wmf (19)

с начальными условиями

boayrkin048.wmf (20)

boayrkin049.wmf (21)

Поставленные задачи могут быть решены последовательно: сначала задача 1, потом задача 2, а затем задача 3.

Установим связь между функциями v(x, y, z, tn+1), w(x, y, z, tn+1), f(x, y, z, tn+1) и решением u(x, y, z, tn+1) исходной задачи на (n + 1) слое.

Пользуясь формулой Тейлора, для функции v(x, y, z, tn+1) получим соотношение

boayrkin050.wmf (22)

(23)

где E – единичный оператор.

Для функций w(x, y, z, tn+1) и f(x, y, z, tn+1) можно получить аналогичные выражения:

boayrkin051.wmf (24)

boayrkin052.wmf (25)

Воспользуемся дифференциальным условием (21) для задачи 3:

boayrkin053.wmf

Пользуясь формулой Тейлора, получим соотношение

boayrkin054.wmf

Учитывая условие (20), получим

boayrkin055.wmf

Аналогично

boayrkin056.wmf

boayrkin057.wmf

Рассмотрим выражение

boayrkin058.wmf

С учетом х начальных данных для вспомогательных задачи 1, задачи 2 и задачи 3 получим

boayrkin059.wmf

и так как

boayrkin060.wmf

boayrkin061.wmf

boayrkin062.wmf

boayrkin063.wmf,

получим

boayrkin064.wmf (26)

Следовательно, можем положить, что

boayrkin065.wmf (27)

с соответствующей аппроксимационной оценкой.

Процесс решения исходной задачи с тремя пространственными переменными x, y, z заменен процессом решения трех задач с одной пространственной переменной в силу расщепления дифференциального оператора A на сумму трех операторов A1 + A2 + A3.

Очевидно, что метод расщепления не увеличивает погрешность по пространственным переменным x, y, z. Покажем, что и по времени t аппроксимационная оценка имеет тот же порядок, что и при прямом разностном методе.

Для каждой задачи (1, 2, 3) построим явные разностные схемы на той же сетке, что и для исходной многомерной задачи:

boayrkin066.wmf

где xi = ih, yj = jh, zk = kh, tn = nτ, i = 0, ±1, ±2, ... , j = 0, ±1, ±2, ... , k= 0, ±1, ±2, ... , n= 0, 1, 2, ... .

Рассмотрим задачу 1.

Заменим производную vtt(xi, yj, zk, tn) её разностной аппроксимацией. Пользуясь формулой Тейлора, получим следующие соотношения:

boayrkin067.wmf (28)

boayrkin068.wmf (29)

Сложим (28) и (29) и выразим vtt(xi, yj, zk, tn).

boayrkin069.wmf (30)

boayrkin070.wmf (31)

Аналогично получим аппроксимацию для vxx(xi, yj, zk, tn).

boayrkin071.wmf (32)

boayrkin072.wmf boayrkin073.wmf boayrkin074.wmf (33)

Из (31) и (32) получим

boayrkin075.wmf (34)

Оценим внутреннюю невязку boayrkin076.wmf.

boayrkin077.wmf (35)

Пусть boayrkin078.wmf, тогда получим

boayrkin079.wmf (36)

Таким образом получаем, что разностная схема (34) аппроксимирует исходную задачу 1 со вторым порядком аппроксимации по τ.

Начальное условие boayrkin080.wmf задачи 1 можно представить в следующем виде:

boayrkin081.wmf (37)

Тогда разностная схема для задачи 1 будет иметь следующий вид:

boayrkin082.wmf (38)

Аналогично для вспомогательных задач 2 и 3 соответственно получим

boayrkin083.wmf (39)

boayrkin084.wmf (40)

Разностные схемы (38), (39) и (40) аппроксимируют соответствующие им задачи со вторым порядком аппроксимации по h и τ.

Для доказательства устойчивости каждой из разностных схем воспользуемся спектральным признаком Неймана.

Рассмотрим разностную схему (38) для задачи 1:

boayrkin085.wmf

Определим начальное условие в виде трехмерной гармоники, зависящей от трех вещественных параметров α, β, γ:

boayrkin086.wmf (41)

Тогда решение задачи 1 при начальном условии (41) имеет вид

boayrkin087.wmf (42)

Функция boayrkin088.wmf является собственной функцией разностного оператора

boayrkin089.wmf (43)

λ(α, β, γ) – соответствующее собственное число разностного оператора.

Подставив это выражение в разностную схему, получим

boayrkin090.wmf (44)

Поделим полученное равенство на boayrkin091.wmf, получим

boayrkin092.wmf (45)

Обозначим boayrkin093.wmf и заметим, что

boayrkin094.wmf (46)

получим

boayrkin095.wmf (47)

boayrkin096.wmf (48)

Для выполнения условия устойчивости необходимо, чтобы спектр разностного оператора λ(α, β, γ) лежал в единичном круге, т.е. boayrkin097.wmf.

Произведение корней этого уравнения по теореме Виета равно единице. Если дискриминант

boayrkin098.wmf (49)

квадратного уравнения отрицателен, то корни λ1(α, β, γ) и λ2 (α, β, γ) комплексно-сопряженные и равны единице по модулю.

В случае r < 1 дискриминант остается отрицательным при всех α. В этом случае спектр заполняет часть единичной окружности.

В случае r = 1 спектр заполняет всю окружность.

При r > 1, по мере увеличения α от 0 до π, корни λ1(α, β, γ) и λ2(α, β, γ) двигаются из точки λ = 1 по единичной окружности: один по часовой, а другой против часовой стрелки соответственно и сходятся в точке λ = –1. Затем один из корней перемещается по вещественной оси из точки λ = –1 влево, а другой вправо, т.к. они вещественны и λ1∙λ2 = 1.

Условие устойчивости выполнено при r ≤ 1.

Проведя аналогичные подстановки и преобразования для разностных схем (45) и (46), соответствующих вспомогательным задачам 2 и 3, получим

boayrkin099.wmf – для разностной схемы (45),

boayrkin100.wmf – для разностной схемы (46).

По рассуждениям, аналогичным для (41)–(49), получаем, что условие устойчивости для разностных схем (39) и (40) выполняется также при r ≤ 1.

Рецензенты:

Камбург В.Г., д.т.н., профессор кафедры «Информационно-вычислительные системы», Пензенский государственный университет архитектуры и строительства, г. Пенза;

Монастырев П.В., д.т.н., профессор, директор института архитектуры, строительства и транспорта, Тамбовский государственный технический университет, г. Тамбов.

Работа поступила в редакцию 06.03.2015.