Development of mathematical models and research of strongly nonequilibrium phenomena taking into account space-time nonlocality

Full Text

Введение. Изучение различного рода локально-неравновесных систем и процессов переноса (энергии, массы, импульса или их аналогов) основано на различных моделях: локально-неравновесной термодинамики (Sybren R. de Groot, Peter Mazur [1], Istv´an Gyarmati [2], William A. Day [3], Nikola Petrov, Jordan G. Brankov [4], David Jou, Jos´e Casas-V´azquez, Georgy Lebon [5, 6], С. Л. Соболев [7-10]); кинетических (James C. Maxwell [11], В. А. Бубнов [12], Borivoje B. Mikic [13], D. Urushev, M. Borisov, A. Vavrek [14]); феноменологических (Morton E. Gurtin, A. C. Pipkin [15], Jace W. Nunziato [16], Helmut D. Weymann [17], Yehuda Taitel [18], Е. И. Леванов, Е. Н. Сотский [19]). В некоторых из перечисленных направлений рассматривается понятие «тепловой памяти», среды с внутренними степенями свободы, среды скоростного типа. Существуют методы, основанные на предположении, что удельная энтропия является одновременно функцией равновесных экстенсивных параметров и их потоков - «скоростных переменных», описывающих инерционные свойства систем при установлении локального равновесия. Уравнения локальнонеравновесного переноса получают также из уравнения Больцмана путем использования теории случайных блужданий, а также молекулярно-кинетическими методами. В различных версиях перечисленных теорий решается вопрос о том, что понимать под температурой локально-неравновесных систем. В классической термодинамике температура определяется только для равновесных или локально-равновесных систем. В различных версиях локально-неравновесной термодинамики вводятся понятия локально-неравновесных энтропии и внутренней энергии, и локально-неравновесная температура определяется как производная от внутренней энергии по энтропии (David Jou, Jos´e Casas-V´azquez, Georgy Lebon [5, 6]). В средах скоростного типа, когда состояние системы зависит от скорости изменения температуры, абсолютная температура заменяется на «термодинамическую» (Nikola Petrov, Jordan G. Brankov [4]) Обе эти температуры совпадают, если абсолютная температура неизменна во времени. В случае, когда температура изменяется во времени, термодинамическая температура представляется в виде разложения в ряд Тейлора по степеням абсолютной температуры. Так как величина абсолютной температуры пропорциональна скорости ее изменения, величина отклонения термодинамической температуры от абсолютной характеризует степень локальной неравновесности системы. В молекулярной динамике вводится понятие «кинетической» температуры как меры средней кинетической энергии молекул, определяемой по локально-неравновесной функции распределения (С. Л. Черешнев, А. П. Генич, С. В. Куликов, Г. Б. Манелис [20], Rosemary A. MacDonald, Donald H. Tsai [21]). Необходимо отметить, что из достаточно большого числа версий локально-неравновесных процессов некоторые из них приводят к противоречивым результатам. Так, устраняя проблемы параболических уравнений, связанные с бесконечной скоростью распространения потенциалов исследуемых полей, появляются новые проблемы, связанные со скачками искомой функции, с отрицательными ее значениями в обратной волне, а также с бесконечными величинами искомой функции в обратной волне. В настоящее время наиболее перспективным, по мнению автора данной статьи, является направление исследований, связанное с использованием не 117 К у д и н о в И. В. только временного, но и пространственного микромасштаба релаксации. Для его оценки в некоторых теориях вводится «поток теплового потока», характеризующий влияние пространственного микромасштаба системы путем учета длины свободного пробега микрочастиц. Несмотря на наличие большого числа версий локально-неравновесных процессов переноса тепла, массы, импульса, в настоящее время отсутствует их единая теоретическая база, на основе которой можно было бы построить непротиворечивые теории по любому из направлений переноса. В последнее время в связи с необходимостью проведения различных физических процессов в экстремальных условиях возник интерес к изучению локально-неравновесных процессов. Интерес к таким процессам обусловлен возможностями познания новых закономерностей и свойств материалов, которые могут быть широко использованы для создания новых технологий получения наноматериалов и покрытий с уникальными физико-химическими свойствами. Однако для описания процессов, протекающих в локально-неравновесных условиях, классическая термодинамика неравновесных процессов не может быть использована, так как она основана на принципе локального равновесия, справедливого для относительно медленных процессов. Возникает проблема создания адекватных физических и соответствующих математических моделей описания локально-неравновесных процессов на основе современной версии термодинамики третьего типа - расширенной необратимой термодинамики, выходящей за рамки принципа локального равновесия. 1. Основные положения предлагаемого метода исследования. Известные математические постановки краевых задач математической физики основаны на бесконечной скорости распространения потенциалов исследуемых полей (задачи теплопроводности, теплообмена в жидкостях, динамической термоупругости, колебаний упругих тел и жидкостей и проч.), которая заложена в соответствующих диффузионных уравнениях: закон Фурье (для теплового потока), Фика (для задач диффузии), Гука (для нормального напряжения), Ньютона (для касательного напряжения в жидкости), Ома (для распространения тока в проводниках) и др. Полученные на основе этих законов дифференциальные уравнения являются локальными, то есть в них не учитывается пространственно-временная нелокальность. Например, математические постановки классических задач нестационарного тепломассопереноса основаны на параболических уравнениях, вывод которых базируется на уравнении теплового баланса (законе сохранения энергии) и феноменологическом законе Фурье для теплового потока (гипотезе Фурье): q = -λ ∂T . ∂x (1) Здесь q - тепловой поток; λ - коэффициент теплопроводности; T - температура; x - координата. Исследования аналитических решений полученных таким путем уравнений приводят к заключению о том, что скорость распространения тепла в изучаемых средах является бесконечной. Для большинства практических задач нестационарной теплопроводности классическая феноменология Фурье приводит к удовлетворительным результатам, подтверждаемым опытными данными. Однако существует большой класс задач, для которых классическая 118 Разработка математических моделей и исследование неравновесных явлений. . . модель теплопроводности оказывается несостоятельной. Это, например, краевые задачи при малых и сверхмалых значениях временной и пространственной переменных, задачи для быстропротекающих процессов, время протекания которых сопоставимо со временем релаксации, и др. Основная причина неприменимости классической модели к указанным процессам в том, что закон Фурье является лишь приближенным описанием процесса теплопроводности, не учитывающим инерционность теплового потока и скалярной величины градиента температуры. Первыми, кто ввел инерционность в уравнения переноса, были Максвелл, Вернотт, Каттанео. В частности, ими была предложена формула для теплового потока, известная как формула Максвелла- Каттанео, путем использования в законе Фурье релаксационного слагаемого, учитывающего конечную скорость распространения теплоты [6-8, 22-24]: q = -λ ∂T ∂q - τ1 . ∂x ∂t (2) Здесь t - время; τ1 - коэффициент релаксации теплового потока q. Подставлением (2) в уравнение теплового баланса cρ ∂q ∂T = , ∂t ∂x (3) где c - теплоемкость; ρ - плотность, выводится классическое гиперболическое уравнение теплопроводности в безразмерном виде: ∂Θ(ξ, Fo) ∂ 2 Θ(ξ, Fo) ∂ 2 Θ(ξ, Fo) + Fo1 = . ∂ Fo ∂ξ 2 ∂ Fo2 (4) Здесь at x T - Tw , Fo = 2 , ξ = ; T0 - Tw δ δ T - температура, Tw - температура стенки, T0 - начальная температура, t - время, a - коэффициент температуропроводности, δ - определяющий размер, например, толщина пластины, если уравнение (4) рассматривать как уравнение теплопроводности для пластины; Θ - безразмерная температура, Fo - безразмерное время, ξ - безразмерная координата. Точное аналитическое решение уравнения (4) при однородных граничных условиях первого рода Θ= Θ(0, Fo) = 0, Θ(1, Fo) = 0 и начальном условии Θ(ξ, 0) = 1 имеет вид ∞ Θ(ξ, Fo) = exp(s Fo) B1k cos(βk Fo)- k=1 - B2k sin(βk Fo) cos 2k - 1 (1 - 2ξ)π , (5) 2 119 К у д и н о в И. В. где βk B2k , s B2k = 4s , (2k - 1)βk π 4 Fo1 νk - 1 , 2 Fo1 s=- 1 , 4 Fo1 B1k = √ βk = Fo1 = aτ1 ; δ2 τ1 - коэффициент релаксации теплового потока; Fo1 - безразмерный коэффициент релаксации теплового потока. Многочисленные исследования решения (5) позволили заключить, что уравнение (4), так же как и параболическое уравнение (то есть при Fo1 = 0), не обеспечивает достаточно адекватного описания исследуемых процессов. И, в частности, решения гиперболического уравнения (4) приводят к возникновению скачков температуры, что равнозначно возникновению изотерм внутри тела и бесконечных тепловых потоков на фронте тепловой волны, а также к появлению отрицательных температур в обратной волне. Эти результаты свидетельствуют о нарушении принципа максимального значения в краевых задачах и закона сохранении энергии (рис. 1). Анализ причин получаемых из решения (5) результатов позволил сделать заключение о необходимости учета инерционности не только теплового потока (временная нелокальность), но и скалярной величины градиента температуры (пространственная нелокальность). Полученная таким путем формула (1) закона Фурье для теплового потока принимает вид q + τ1 ∂q ∂2q ∂T ∂2T ∂3T + τ12 2 + · · · = -λ + τ2 + τ22 + ··· , ∂t ∂t ∂x ∂x∂t ∂x∂t2 (6) где τ1 - коэффициент релаксации теплового потока; τ2 - коэффициент релаксации скалярной величины градиента температуры. Аналогичные урав- Рис. 1. Распределение температуры Θ для различных чисел Fo (значения чисел - метки у кривых) в пластине при n = 10 000, где n - число членов ряда точного аналитического решения, для Fo1 = 0.0015 (слева) и Fo1 = 0.1 (справа) [Figure 1. The temperature distribution in the plate for different Fourier numbers (Fo) (the values of the Fourier numbers are the labels of the curves) when n = 10 000 (n is the number of terms in the series of the exact analytical solution) for Fo1 = 0.0015 (left) and Fo1 = 0.1 (right)] 120 Разработка математических моделей и исследование неравновесных явлений. . . нения для законов Ньютона τ = µ∂ϑ/∂y и Гука σ = E∂U/∂x имеют вид 2 ∂τ ∂ϑ ∂2ϑ ∂3ϑ 2 ∂ τ 2 + τ1N + · · · = µ + τ + τ + · · · ; (7) 2N 2N ∂t ∂t2 ∂y ∂y∂t ∂y∂t2 2 3 ∂U ∂2U ∂σ 2 ∂ σ 2 ∂ U + · · · = E + · · · , (8) + τ1G + τ + τ σ + τ1G 2G 2G ∂t ∂t2 ∂x ∂x∂t ∂x∂t2 где τ , σ - касательное и нормальное напряжения; ϑ - скорость; U - перемещение; µ - динамическая вязкость; E - модуль упругости; τ1N , τ1G - коэффициенты релаксации напряжений τ и σ; τ2N , τ2G - коэффициенты релакса∂U ции градиентов скорости и деформации ∂ϑ ∂y и ∂x . Выражая q, τ , σ, ограничиваясь двумя членами правой и левой части соотношений (6)-(8), находим τ + τ1N ∂T ∂2T ∂q - λτ2 - τ1 ; (9) ∂x ∂x∂t ∂t ∂ϑ ∂2ϑ ∂τ τ =µ + µτ2N - τ1N ; (10) ∂y ∂y∂t ∂t ∂U ∂2U ∂σ σ=E + Eτ2G - τ1G . (11) ∂x ∂x∂t ∂t Отметим, что соотношения (9)-(11) совпадают с формулами, полученными из сформулированной А. В. Лыковым системы дифференциальных уравнений Онзагера (Lars Onsager) (исходя из гипотезы о конечной скорости диффузии массы и теплоты) [24]: q = -λ Ji = (r) ∂Ji Li ∂t N + Lik Xk + Lik k=1 ∂Xk , ∂t (12) где Ji - поток субстанции (теплоты, массы и т.д.); Xk - термодинамические (r) движущие силы (градиенты соответствующих величин); Li , Lik , Lik - постоянные феноменологические коэффициенты переноса (Lik = Lki ). Если положить L = -λ, L(r) = -τ1 , L = -λτ2 , Xk = ∂T ∂x , то из (12) получаем формулу (9). Если положить L = µ, L = µτ2N , L(r) = -τ1N , Xk = ∂ϑ ∂x , то уравнение (12) оказывается тождественным уравнению (10) [24]. При L = E, L = Eτ2G , L(r) = -τ1G из (12) получаем формулу (11). Соотношения (10), (11) полностью совпадают с усложненными моделями вязкоупругого тела, известными как модели Максвелла и Кельвина-Фойхта [25]. Эти модели отличаются лишь различными формулами для коэффициентов релаксации. Их физический смысл в том, что в них учитывается временная зависимость напряжений и деформаций и их взаимное влияние друг на друга. Совпадение моделей (10), (11) с моделями Максвелла, Кельвина-Фойхта (с точностью до постоянных) свидетельствует об использовании одинаковых исходных принципов. 2. Математические модели теплопроводности с учетом пространственно-временной нелокальности. Рассмотрим последовательность построения локально-неравновесных моделей теплопроводности на основе формулы (6) для теплового потока в случае бесконечной пластины с толщиной δ. 121 К у д и н о в И. В. Подставляя (6) в уравнение теплового баланса (3), при учете трех членов правой и левой части соотношения (6) получим уравнение теплопроводности для бесконечной пластины в безразмерном виде [26-29]: ∂Θ ∂2Θ ∂3Θ ∂2Θ ∂3Θ ∂4Θ + Fo1 + Fo2 = + Fo1 2 + Fo2 2 , 2 3 2 ∂ Fo ∂ξ ∂ξ ∂ Fo ∂ Fo ∂ Fo ∂ξ ∂ Fo2 (13) где T - Tw x aτ1 a2 τ 2 at , ξ = , Fo = 2 , Fo1 = 2 , Fo2 = 42 ; T0 - Tw δ δ δ δ T - температура, Tw - температура стенки пластины, T0 - начальная температура, t - время, a - коэффициент температуропроводности, δ - толщина пластины, τ1 - коэффициент релаксации теплового потока, τ2 - коэффициент релаксации скалярной величины градиента температуры; Θ - безразмерная температура, Fo - безразмерное время, ξ - безразмерная координата, Fo1 - безразмерный коэффициент релаксации скалярной величины градиента температуры, Fo2 - безразмерный коэффициент релаксации скалярной величины градиента температуры. Если положить Fo1 = Fo2 = 0 (при τ1 = τ2 = 0), то уравнение (13) приводится к классическому параболическому уравнению теплопроводности. Точное аналитическое решение уравнения (13) для бесконечной пластины при граничных условиях первого рода Θ= Θ(0, Fo) = 0, Θ(1, Fo) = 0 и начальном условии Θ(ξ, 0) = 1 имеет вид ∞ Θ(ξ, Fo) = D1k exp(z1k Fo) + D2k exp(z2k Fo)+ k=1 + D3k exp(z3k Fo) cos 2k - 1 (1 - 2ξ)π , (14) 2 где z1k , z2k , z3k - корни характеристического уравнения; D1k , D2k , D3k - константы интегрирования, определяемые из начальных условий краевой задачи. Результаты расчетов по формуле (14) даны на рис. 2. Их анализ позволяет заключить, что при малых значениях чисел Fo1 и Fo2 получаемые решения в диапазоне чисел 10-6 Fo < ∞ полностью совпадают с решением аналогичной задачи для параболического уравнения нестационарной теплопроводности. При Fo < 10-6 распределение температуры во времени имеет следующие особенности. При значениях величин Fo1 = 10-6 , Fo2 = 10-12 , например, при Fo = 10-4 температура тела в бесконечно малых окрестностях точек ξ = 0 и ξ = 1 скачкообразно изменяется от величины Θ(0, 10-4 ) = Θ(1, 10-4 ) = 0 непосредственно на границах (выполнение граничных условий первого рода) до Θ(0.0001, 10-4 ) = Θ(0.9999, 10-4 ) = 0.91 в точках ξ = 0.0001 и ξ = 0.9999, т. е. на некотором малом расстоянии от точек ξ = 0 и ξ = 1 (см. рис. 2 сверху). 122 Разработка математических моделей и исследование неравновесных явлений. . . При увеличении чисел Fo1 , Fo2 температурные кривые становятся более пологими и при Fo1 > 5 - практически параллельными оси ξ. Скачок температуры на стенках (в окрестностях точек ξ = 0 и ξ = 1) в этом случае имеет место практически для всех чисел Fo вплоть до наступления стационарного состояния (см. рис. 2 снизу слева). Изменение безразмерного теплового потока q¯ = qδ/(λ(T0 -Tw )), найденное по формулам (1), (2), (9), приведено на рис. 3. Анализ результатов расчетов позволяет сделать следующие выводы. Тепловой поток, определяемый по формуле (1), непосредственно на стороне пластины ξ = 1 при Fo = 0 равен бесконечно большой величине (кривая 1 на рис. 3). С увеличением времени величина q¯(Fo) уменьшается по экспоненциальной зависимости и при Fo ≈ 3 становится практически равной нулю. Тепловой поток, определяемый по формуле, представляющей решение уравнения (9), при ξ = 1 и Fo = 0 равен нулю. С увеличением времени он возрастает и при Fo = Fo∗ ≈ 0.072 достигает максимального значения q¯ ≈ 0.37 (кривая 2 на рис. 3). Изменение теплового потока, определяемого формулой (2), иллюстрируется кривой 3 на рис. 3. В данном случае тепловой поток на стороне ξ = 1, имея нулевое значение при Fo = 0, с увеличением числа Fo возрастает до некоторого максимального значения, равного q¯ ≈ 1.83 (при Fo = Fo ≈ 0.08). При дальнейшем увеличении числа Fo тепловой поток изменяется скачкообразно, принимая в некоторых диапазонах числа Fo отрицательные значения. Такое изменение теплового потока приводит к появлению скачков температуры и к их отрицательным значениям. Рис. 2. Распределение температуры Θ для различных чисел Fo (значения чисел - метки у кривых) в пластине (см. формулу (14)) для Fo1 = 10-6 , Fo2 = 10-12 , n = 100 000 (сверху) и Fo1 = 5, Fo2 = 25, n = 1 000 (снизу слева); n - число членов ряда (14) [Figure 2. The temperature distribution in the plate (see Eq. (14)) for different Fourier numbers(Fo) (the values of the Fourier numbers are the labels of the curves) for Fo1 = 10-6 , Fo2 = 10-12 , n = 100 000 (top) and Fo1 = 5, Fo2 = 25, n = 1 000 (bottom left); n is the number of terms in the series (14)] 123 К у д и н о в И. В. Рис. 3. Изменение теплового потока q¯ на стороне пластины ξ = 1 при Fo1 = 0.3 по формуле (1) (кривая 1), по формуле (9) (кривая 2) и по формуле (2) (кривая 3) [Figure 3. The change in the heat flux (¯ q ) on the side ξ = 1 of the plate when Fo1 = 0.3 by Eq. (1) (curve 1), by Eq. (9) (curve 2), and by Eq. (2) (curve 3)] Обнаруженный из анализа точного аналитического решения (14) факт задержки во времени принятия граничного условия первого рода свидетельствует о том, что ввиду сопротивления тела, оказываемого процессу проникновения теплоты, его мгновенный прогрев на границе невозможен ни при каких условиях теплообмена с окружающей средой. Следовательно, коэффициент теплоотдачи на стенке зависит не только от условий теплообмена (скорость среды, вязкость и проч.), но и от физических свойств тела, и он, во-первых, является переменной во времени величиной и, во-вторых, он не может превысить некоторой предельной для каждого конкретного случая величины. Конкретный пример определения коэффициента теплоотдачи для стальной пластины (λ = 49 Вт/(м · К)) толщиной δ = 0.01 м при Fo1 = 0.3, Fo2 = 0, который оказался равным α = 4005 Вт/(м2 · К), приведен в [29]. Таким образом, механизм теплопроводности с учетом релаксационных свойств материалов ограничивает количество тепла, отводимого от системы. Этот факт необходимо учитывать при рассмотрении всех быстропротекающих процессов - при решении задач теплового воспламенения (теплового взрыва), задач динамической термоупругости, где в условиях теплового удара классическая теория приводит к бесконечным динамическим напряжениям, при лазерной обработке и упрочнении материалов и во многих других процессах. 3. Математическая модель колебаний давления в несжимаемой жидкости с учетом пространственно-временной нелокальности. С учетом соотношения (9) авторами работы [27] было получено дифференциальное уравнение и сформулирована краевая задача, описывающая колебание давления упругой жидкости в трубопроводе в условиях гидравлического удара. Математическая постановка задачи в данном случае имеет вид For ∂ 3 Θ(ξ, Fo) ∂Θ(ξ, Fo) ∂ 2 Θ(ξ, Fo) + + Fo1 = 2 ∂ Fo ∂ Fo ∂ Fo3 ∂ 2 Θ(ξ, Fo) ∂ 3 Θ(ξ, Fo) = + Fo , 2 ∂ξ 2 ∂ξ 2 ∂ Fo Θ(ξ, Fo) 124 Fo 0, ∂Θ(ξ, Fo) ∂ 2 Θ(ξ, Fo) = 0, Fo=0 ∂ Fo Fo=0 ∂ Fo2 ∂Θ(ξ, Fo) = 0, Θ(ξ, Fo) ξ=1 = 0. ∂ξ ξ=0 = 1, 0 Fo=0 ξ = 0; 1; (15) (16) (17) Разработка математических моделей и исследование неравновесных явлений. . . В данной модельной задаче используются следующие обозначения: Θ= Fo1 = cs τ1 , l p - p1 , p0 - p1 Fo2 = Fo = cs τ2 , l cs t , l For = ξ =1- 2bl , cs x , l 2b = 32ν ; d2 p - давление, p0 - начальное давление, p1 - давление на входе в трубопровод, x - координата, l - длина трубы, cs - скорость звука в жидкости; ν - вязкость жидкости, d - диаметр трубы; τ1 - коэффициент релаксации касательного напряжения, τ2 - коэффициент релаксации градиента скорости; Θ - безразмерное давление, Fo - безразмерное время, ξ - безразмерная координата, Fo1 - безразмерный коэффициент релаксации касательного напряжения, Fo2 - безразмерный коэффициент релаксации градиента скорости, For - безразмерный коэффициент сопротивления среды. Краевая задача (15)-(17) описывает распространение скачка давления в трубопроводе с неподвижной в исходном состоянии жидкостью для случая, когда в сечении ξ = 1 произошло скачкообразное изменение давления, а сечение ξ = 0 перекрыто (скорость течения жидкости равна нулю). Точное аналитическое решение задачи (15)-(17) находится в виде ∞ Θ(ξ, Fo) = C1k exp (z1k Fo) + C2k exp (z2k Fo) + k=1 + C3k exp (z3k Fo) cos 2k - 1 πξ , (18) 2 где C1k , C2k , C3k - константы интегрирования, определяемые из начальных условий (16); z1k , z2k , z3k - корни характеристического уравнения. Результаты расчетов по формуле (18) приведены на рис. 4 (для ξ = 0). Здесь даны результаты решения задачи (15)-(17) в квазистационарной постановке (при Fo1 = Fo2 = 0), а также для случая, когда учитывается конечная скорость распространения потенциалов исследуемой среды (Fo1 = 0, Fo2 = 0), то есть с учетом влияния нестационарности процесса. На этих рисунках представлены также результаты расчетов, выполненных автором работы [31], и результаты экспериментальных исследований гидравлического удара для сильновязкой жидкости (масло) при ламинарном режиме течения [32]. Теоретические исследования в [31] выполнялись по двум методикам - с учетом влияния нестационарности процесса на силу трения и без ее учета (в квазистатической постановке). Анализ результатов расчетов по формуле (18) позволяет сделать вывод о том, что решение задачи (15)-(17) в квазистатической постановке (при Fo1 = Fo2 = 0) незначительно отличается от решения аналогичной задачи, полученного в [31] с использованием метода Бернулли-Фурье. Однако результаты двух теоретических методов значительно отличаются от экспериментальных данных. Учет нестационарности градиента скорости и касательного напряжения приводит к существенному приближению расчетных данных к результатам эксперимента. На основе экспериментальных данных, 125 К у д и н о в И. В. Рис. 4. Распределение давления в трубопроводе при гидравлическом ударе (ξ = 0), когда For = 0.0538 и Fo1 = Fo2 = 0 (слева) и For = 0.0538, Fo1 = 0.4, Fo2 = 0.043 (справа): 1 - расчет по формуле (18) при n = 100; 2 - по методу [31]; 3 - эксперимент [32] [Figure 4. The pressure distribution in a pipeline under the conditions of a hydraulic shock (ξ = 0) when For = 0.0538, Fo1 = Fo2 = 0 (left), and For = 0.0538, Fo1 = 0.4, Fo2 = 0.043 (right): curve 1 - calculation by Eq. (18) when n = 100; curve 2 - by method [31]; curve 3 - experimental data [32]] используя соотношение (18), путем решения обратной задачи были найдены коэффициенты релаксации: τ1 = 0.0108 с, τ2 = 0.0116 с. Отметим, что в классическом гиперболическом уравнении, описывающем распространение скачка давления при гидравлическом ударе, когда Fo1 = = Fo2 = 0, зависимости между значениями скорости, плотности и давления потока жидкости принимаются осредненными по его сечению. Кроме этого, они содержат среднюю во времени величину касательного напряжения на стенке трубопровода, то есть используется гипотеза квазистационарности, в соответствии с которой считается, что закон Ньютона для касательного напряжения выполняется и при нестационарном процессе изменения скорости. Предполагается также, что соответствующие соотношения сохраняются и для коэффициента гидравлического сопротивления. Строгое обоснование этого допущения отсутствует, поэтому его справедливость может быть подтверждена лишь путем сравнения результатов расчета с экспериментальными исследованиями. Полученные двумя различными методами результаты позволяют заключить, что гипотеза квазистационарности, вообще говоря, неверна, что подтверждается приведенными выше сравнениями решений квазистационарной и нестационарной задач с результатами натурного эксперимента. Из двух рассмотренных выше направлений учета нестационарности касательных напряжений первое связано с рассмотрением дифференциальных уравнений Навье-Стокса, описывающих изменение местных (локальных - неосредненных) величин [31]. Второе направление, рассмотренное в данной работе, связано с введением временн´oй зависимости между градиентом скорости и касательным напряжением в формуле (10) закона Ньютона. В данном случае устраняется проблема бесконечной скорости распространения потенциалов исследуемого поля, а также проблема бесконечных величин касательных напряжений на стенке трубопровода при Fo → 0. Анализ безразмерного касательного напряжения τ¯, определяемого на основе формулы (10), которая при τ1 = τ2 = 0 приводится к формуле закона Ньютона, позволяет заключить, что τ¯ → ∞ при Fo → 0 (см. рис. 5). При τ1 126 Разработка математических моделей и исследование неравновесных явлений. . . Рис. 5. Изменение касательного напряжения на стенке трубопровода при For = 2al/c = 0.3; решение уравнения (10) при τ1 = τ2 = 0 (линия 1) и при τ1 = 0.0108, τ2 = 0.0116 (линия 2) [Figure 5. The change of the tangential stress (¯ τ ) distribution on the pipeline wall, which was found as solution of Eq. (10) when τ1 = τ2 = 0 (curve 1), and τ1 = 0.0108, τ2 = 0.0116 (curve 2); here For = 2al/c = 0.3] и τ2 , отличных от нуля, касательное напряжение τ¯, имея в начальный момент времени нулевое значение, с увеличением времени возрастает, принимая при некотором Fo = Fo∗ максимальное значение τ¯ = τ¯∗ . При дальнейшем возрастании времени τ¯ уменьшается и при Fo → ∞ τ¯ → 0. Полученные результаты применительно к решению уравнения (10) позволяют сделать следующие выводы: на стенке трубопровода касательное напряжение τ¯ → 0 при Fo → 0 и не может превысить некоторой максимальной для данных конкретных условий (определяемых физическими свойствами среды) величины; применение формулы (10) позволяет избежать бесконечных значений касательных напряжений. Изменение давления в случае Fo1 = Fo2 = 0 (квазистационарный режим) приведено на левом рис. 4 (линия 1). Анализ результатов позволяет сделать заключение о скачкообразном изменении давления в любой точке трубы во времени, свидетельствующем о бесконечной скорости передачи импульса. Отметим, что скорость перемещения волны при этом вполне конечна и равна скорости звука в жидкости. Бесконечной является скорость передачи импульса, заложенная в классической (нерелаксированной) формуле закона Ньютона для касательного напряжения. Полученные результаты позволяют сделать вывод об ограниченности математических моделей колебаний упругих жидкостей, в которых не учитываются их релаксационные свойства, ввиду существенного отклонения получаемых данных от результатов экспериментальных исследований (см. рис. 4 слева). Учет релаксационных свойств позволяет через модификацию закона Ньютона для касательного напряжения учесть нестационарность изменения скорости и касательного напряжения. 4. Математическая модель теплообмена в движущейся жидкости с учетом ее релаксационных свойств. В предложенной выше математической модели (15)-(17) исследуется изменение давления в жидкости в условиях гидравлического удара. Рассмотрим вывод уравнения применительно к нестационарному теплообмену при ламинарном течении вязкой несжимаемой жидкости в плоскопараллельном канале толщиной 2δ (см. рис. 6). Формулы для теплового потока в направлении координатных осей x и y по аналогии с (9) будут иметь вид [33] ∂T ∂qx ∂2T - τ1 - λτ2 + ρωx h; ∂x ∂t ∂x∂t (19) ∂qy ∂T ∂2T qy = - λ - τ1 - λτ2 + ρωy h, ∂y ∂t ∂y∂t где qx , qy - тепловые потоки по осям x и y; ωx , ωy - скорости течения жидqx = - λ 127 К у д и н о в И. В. Рис. 6. Схема теплообмена при ламинарном течении жидкости в плоскопараллельном канале [Figure 6. A heat exchange scheme for laminar flow of liquid in the plane-parallel channel; Ti is the initial liquid temperature; T0 is the liquid temperature at the channel inlet; Tw is the liquid temperature on the channel wall; 2δ is the channel width; ωx is the fluid flow velocity in the direction of the x-axis] кости в направлении координатных осей x и y; h - теплосодержание (энтальпия жидкости); λ - коэффициент теплопроводности; x, y - продольная и поперечная координаты плоскопараллельного канала; τ1 , τ2 - коэффициенты релаксации соответственно теплового потока и градиента температуры. Отличие формул (19) от соотношения (9) состоит в том, что при движении жидкости теплота переносится не только теплопроводностью, но и конвекцией, учитываемой последними слагаемыми этих формул. Уравнение теплового баланса в данном случае имеет вид ρ ∂T ∂qx ∂qy ∂h = cρ =- + . ∂t ∂t ∂x ∂y (20) Подставляя (19) в (20), пренебрегая конвективным переносом теплоты в поперечном направлении (ρωy ∂h ∂y = 0) и теплопроводностью - в продольном 2 (λ ∂∂xT2 = 0), получаем ∂T ∂T ∂2T ∂2T ∂ + ωx = a 2 - τ1 2 + aτ2 ∂t ∂x ∂y ∂t ∂t ∂2T ∂2T + ∂x2 ∂y 2 , (21) где λ 3 y2 , ωx = ωav 1 - 2 ; cρ 2 δ ωav - средняя скорость течения жидкости. В данной модели будем использовать следующие обозначения: Ti - начальная температура жидкости; T0 - температура жидкости на входе в канал; Tw - температура жидкости на стенке; 2δ - ширина плоского канала, а также следующие безразмерные переменные и параметры: a= T - Tw ; T0 - Tw aτ2 Fo∗2 = 2 ; δ Θ= y at aτ1 ; Fo = 2 ; Fo∗1 = 2 ; δ δ δ 2 ax 2 a 2 ξ= ; Pe = . 3 δ 2 ωav 3 δωav η= (22) С учетом обозначений (22) уравнение (21) приводится к виду ∂2Θ ∂Θ ∂2Θ ∂2Θ ∂Θ ∂2Θ ∗ ∂ + Fo∗1 + (1 - η) = + Fo Pe + , 2 ∂ Fo ∂ξ ∂η 2 ∂ Fo ∂ξ 2 ∂η 2 ∂ Fo2 Fo > 0, ξ > 0, 0 < η < 1. 128 (23) Разработка математических моделей и исследование неравновесных явлений. . . Начальные и граничные условия имеют вид Θ(ξ, η, 0) = Θi , Θ(0, η, Fo) = 1, ∂Θ(ξ, η, Fo) ∂η η=0 ∂Θ(ξ, η, Fo) ∂ Fo ∂Θ(ξ, η, Fo) ∂ξ = 0, Fo=0 ξ=∞ = 0; = 0; (24) Θ(ξ, 1, Fo) = 0, где Θi = (Ti - Tw )/(T0 - Tw ) - безразмерная начальная температура. Для решения задачи (23), (24) использовался метод конечных разностей. Результаты расчетов (для Θi = 0.5) представлены на рис. 7, 8. На рис. 7 даны результаты, полученные при равенстве нулю последних слагаемых правой и левой частей уравнения (23), которое в этом случае приводится к уравнению (4). Анализ результатов расчетов позволяет сделать заключение о практическом совпадении расчетов численным методом с точным аналитическим решением уравнения (4). На рис. 8 дано распределение температуры жидкости при малых значениях величин поперечной координаты 0.99 η 1.0 и чисел Фурье из диапазона 2 · 10-8 Fo 104 . Анализ полученных результатов позволяет отметить интересную особенность теплообмена при задании граничных условий первого рода (тепловой удар) на стенке в случае, когда учитываются релакРис. 7. Распределение температуры в плоском канале с неподвижной жидкостью (Fo∗1 = Fo∗2 = = 10-7 ) для различных чисел Fo (значения чисел - метки у кривых): сплошные линии - точное решение [26-30]; ◦ - метод конечных разностей [Figure 7. The temperature distribution in a flat channel channel with a fixed liquid volume (when Fo∗1 = Fo∗2 = 10-7 ) for different Fourier numbers (Fo) (the values of the Fourier numbers are the labels of the curves); the solid lines are the exact solution [26-30] and lines of points (◦) are the simulation results by the finite difference method] Рис. 8. Распределение температуры по ширине канала при ξ = 5 · 10-5 , Fo∗1 = Fo∗2 = 10-7 , Pe = 10-7 в диапазоне безразмерного времени 2 · 10-8 Fo 10-4 : сплошные линии - точное решение [26-30]; ◦ - метод конечных разностей [Figure 8. The temperature distribution along the channel width (when ξ = 5 · 10-5 , Fo∗1 = Fo∗2 = 10-7 , Pe = 10-7 ) for 2 · 10-8 Fo 10-4 (the values of the Fourier numbers are the labels of the curves); the solid lines are the exact solution [26-30] and lines of points (◦) are the simulation results by the finite difference method] 129 К у д и н о в И. В. сационные свойства среды. Так, несмотря на выполнение граничного условия первого рода Θ(ξ, 1, Fo) = 0 в точке ξ = 1, в бесконечно малой окрестности этой точки, например, ξ = 0.999, температура стенки не принимает значение Θ(ξ, 1, Fo) = 0 мгновенно, то есть непосредственно при Fo = 0. И, в частности, для того чтобы температура стенки приняла значение Θ(ξ, 1, Fo) = 0, требуется некоторый промежуток времени, равный Fo = 10-4 (при Fo∗1 = Fo∗2 = 10-7 ). Полученные результаты полностью согласуются с результатами точного аналитического решения вида (14) [26-30]. Отметим, что при некоторых больших значениях числа Фурье теплообмен стабилизируется и наступает квазистационарный режим, когда температура жидкости в любой точке канала принимает неизменное во времени значение. 5. Локально-неравновесные модели процессов колебаний упругих твердых тел. Используя модифицированную формулу закона Гука вида (11), рассмотрим последовательность построения локально-неравновесной модели колебаний упругого стержня длины δ с учетом сопротивления, оказываемого материалом процессу колебаний. Уравнение равновесия (движения) в данном случае имеет вид ρ ∂2U ∂σ ∂U = - γρ , ∂t2 ∂x ∂t (25) где γ - коэффициент сопротивления; ρ - плотность. Отметим, что второе слагаемое правой части уравнения (25) представляет силу сопротивления среды, приходящуюся на единицу объема. Подставляя (11) в (25), с учетом того, что ∂σ ∂2U =ρ 2 , ∂x ∂t находим τ1c 2 3 ∂3U ∂2U ∂U 2∂ U 2 c ∂ U + = e + e τ -γ , s s 2 ∂t3 ∂t2 ∂x2 ∂x2 ∂t ∂t (26) где es = E/ρ - скорость распространения волны (скорость звука в материале стержня). Очевидно, что при τ1c = τ2c = γ = 0 уравнение (26) приводится к волновому уравнению незатухающих колебаний: ∂2U ∂2U = e . s ∂t2 ∂x2 Найдем точное аналитическое решение краевой задачи о колебаниях стержня, один торец которого жестко закреплен, а второй находится в свободном состоянии. В начальный момент времени стержень деформирован по линейному закону, согласно которому максимальное перемещение имеет свободный торец стержня. Воспользовавшись безразмерными переменными Θ= 130 U , U0 ξ= x , δ Fo = es t , δ Foc1 = es τ1c , δ Foc2 = es τ2c , δ Foc3 = δγ , es Разработка математических моделей и исследование неравновесных явлений. . . рассматриваемую задачу запишем в виде Foc3 3 ∂ 2 Θ(ξ, Fo) ∂Θ(ξ, Fo) ∂ 3 Θ(ξ, Fo) ∂ 2 Θ(ξ, Fo) c ∂ Θ(ξ, Fo) + = + Fo + Foc1 , 2 3 2 ∂ Fo ∂ξ 2 ∂ξ 2 ∂ Fo ∂ Fo ∂ Fo Fo > 0, 0 < ξ < 1; (27) Θ(ξ, 0) = 1, ∂ 2 Θ(ξ, Fo) ∂Θ(ξ, Fo) = 0, ∂ Fo Fo=0 ∂ Fo2 ∂Θ(ξ, Fo) = 0, Θ(1, Fo) = 0. ∂ξ ξ=0 Fo=0 = 0; (28) (29) Здесь Θ - безразмерное перемещение; ξ - безразмерная координата; Fo - число Фурье (безразмерное время); U0 - начальное перемещение стержня; Foc1 , Foc2 - безразмерные коэффициенты релаксации; Foc3 - безразмерный коэффициент сопротивления среды; δ - длина стержня. Точное аналитическое решение задачи (27)-(29), найденное методом разделения переменных, записывается в виде ∞ exp(γk Fo) B1k cos(βk Fo) - B2k sin(βk Fo) + Θ(ξ, Fo) = k=1 + B3k exp(z3k Fo) cos 2k - 1 πξ , (30) 2 где γk , βk , z3k - корни характеристического уравнения; B1k , B2k , B3k - константы интегрирования, определяемые из начальных условий (28). На рис. 9, 10 приведены результаты расчетов перемещений по формуле (30) при Foc1 = Foc2 = 0.1, Foc3 = 0.5 (рис. 9) и Foc1 = Foc2 = 0.1, Foc3 = 30 (рис. 10). Из их анализа можно заключить, что с увеличением числа Foc3 (при неизменных Foc1 и Foc2 ) колебания смещаются в область положительных значений перемещений и при некоторых больших значениях величины Foc3 (Foc3 > 5) возврат стержня в исходное состояние происходит практически без колебаний (критическое затухание (рис. 10)). На рис. 11 приведены результаты расчетов перемещений по формуле (30) при Foc1 = Foc2 = 10, Foc3 = 0. Их анализ позволяет заключить, что с увеличением Foc1 и Foc2 (при Foc3 = 0) амплитуда колебаний уменьшается. При этом в каждой точке стержня до момента его возвращения в исходное (недеформированное) состояние совершаются многократные колебательные процессы в области лишь положительных значений перемещений. С увеличением времени (Fo 35) процесс колебаний стабилизируется вблизи нулевого значения перемещения. Найдем распределение напряжений в стержне по формуле закона Гука без учета (τ1c = τ2c = 0) и с учетом (τ1c = 0, τ2c = 0) релаксационных слагаемых. Закон Гука (соотношение (11) при τ1c = τ2c = 0) с учетом обозначений Θ = = U/U0 и ξ = x/δ будет иметь следующий вид: σ ¯ (ξ, Fo) = ∂Θ(ξ, Fo) , ∂ξ (31) 131 К у д и н о в И. В. Рис. 9. Изменение перемещений стержня во времени для различных ξ при Foc1 = Foc2 = 0.1, Foc3 = 0.5, n = 100; n - число членов ряда (30) [Figure 9. The change of the rod displacement over time for various ξ when Foc1 = Foc2 = 0.1, Foc3 = 0.5, n = 100; n is the number of terms in the series (30)] Рис. 10. Изменение перемещений стержня во времени для различных ξ при Foc1 = Foc2 = 0.1, Foc3 = 30, n = 100; n - число членов ряда (30) [Figure 10. The change of the rod displacement over time for various ξ when Foc1 = Foc2 = 0.1, Foc3 = 30, n = 100; n is the number of terms in the series (30)] Рис. 11. Изменение перемещений стержня во времени для различных ξ при Foc1 = = Foc2 = 0.1, Foc3 = 0.5, n = 100; n - число членов ряда (30) [Figure 11. The change of the rod displacement over time for various ξ when Foc1 = Foc2 = 0.1, Foc3 = 0, n = 100; n is the number of terms in the series (30)] где σ ¯ (ξ, Fo) = σ(x, t) δ , E U0 Θ(ξ, Fo) - точное аналитическое решение уравнения незатухающих колебаний: ∞ 8 1 2k - 1 2k - 1 Θ(ξ, Fo) = 2 cos π Fo cos πξ . (32) π (2k - 1)2 2 2 k=1 Результаты расчетов безразмерного напряжения по формуле (31) с учетом 132 Разработка математических моделей и исследование неравновесных явлений. . . формулы (32) приведены на рис. 12. Их анализ позволяет заключить, что колебания σ ¯ происходят скачкообразно во времени от σ ¯ = -0.1 до σ ¯ = 0.1. Выпишем соотношение (11) в случае, когда τ1c = 0 и τ2c = 0, в безразмерном виде: ∂σ ¯ (ξ, Fo) + A1 σ ¯ (ξ, Fo) = A2 (ξ, Fo), (33) ∂ Fo где 1 ∂Θ(ξ, Fo) Foc2 ∂ 2 Θ(ξ, Fo) 1 + c . A1 = c , A2 (ξ, Fo) = c Fo1 Fo1 ∂ξ Fo1 ∂ξ∂ Fo Интегрируя уравнение (33), находим σ ¯ (ξ, Fo) = exp(-A1 Fo) A2 (ξ, Fo) exp(A1 Fo)d Fo +C , (34) где C - постоянная интегрирования, определяемая из начального условия σ ¯ (ξ, 0) = 0. В качестве функции Θ(ξ, Fo) в формуле (34) использовалось решение (30). Анализ результатов расчетов напряжения по формуле (34) позволяет заключить, что напряжения различных точек стержня изменяются во времени по гармонической зависимости с затухающей амплитудой (рис. 13). Таким образом, учет релаксационных слагаемых позволяет устранить скачкообразное изменение напряжений во времени. Рис. 12. (онлайн в цвете) Изменение напряжений в стержне во времени при различных значениях ξ, когда Foc1 = Foc2 = 0 (расчет по формуле (31)) [Figure 12. (Color online) The stress change in the rod depending on dimensionless time (Fo) for various ξ when Foc1 = Foc2 = 0 (calculation by Eq. (31))] Рис. 13. (онлайн в цвете) Изменение напряжений в стержне во времени при различных значениях ξ, когда Foc1 = Foc2 = 0.1 (расчет по формуле (34)) [Figure 13. (Color online) The stress change in the rod depending on dimensionless time (Fo) for various ξ when Foc1 = Foc2 = 0.1 (calculation by Eq. (34))] 133 К у д и н о в И. В. 6. Резонансные колебания стержня с учетом релаксационных свойств материалов. Если в задаче (27)-(29) к свободному торцу стержня приложить нагрузку, изменяющуюся по гармоническому закону, то первое граничное условие из (29) примет вид ∂Θ(ξ, Fo) = Fo4 cos(Fo5 Fo), (35) ∂ξ ξ=0 где Fo4 = δ/U0 ; Fo5 = ωδ/es ; ω = 2πνl - круговая частота; νl - частота колебаний нагрузки. Далее будет показано, что при совпадении частоты собственных колебаний стержня с частотой Fo5 внешней нагрузки наблюдаются резонансные колебания стержня, в том числе и бифуркационные (флаттерные). Получение точного аналитического решения краевой задачи (27)-(29) с учетом граничного условия (35) существенно усложняется, и для ее решения был использован конечно-разностный метод. Результаты выполненных исследований позволяют сделать вывод, что при Foc1 = Foc2 = Foc3 = 0 колебания незатухающие. При Foc3 = 0.1 и Foc1 = = Foc2 = Foc4 = 0 колебания становятся затухающими с экспоненциально уменьшающейся амплитудой. При Foc1 = Foc2 = Foc3 = 0 и Fo4 = 0.5, Fo5 = = 1.575 безразмерная частота собственных колебаний стержня νl = 1.575 совпадает с частотой вынужденных колебаний, возникающих под действием внешней нагрузки вида (35). При этом наблюдается неограниченное возрастание амплитуды колебаний (рис. 14). С увеличением коэффициента сопротивления среды Foc3 амплитуда колебаний внутренних точек стержня уменьшается (рис. 15) и при некотором большом его значении (Foc3 100) возврат стержня в исходное состояние происходит практически при отсутствии колебательного процесса внутренних точек стержня (критическое затухание) при неизменной амплитуде колебаний внешней нагрузки. Результаты расчетов резонансных колебаний для Foc1 = Foc2 = 2, Foc3 = 1, Fo4 = 0.1, Fo5 = 1.575 приведены на рис. 16. При дальнейшем увеличении коэффициентов релаксации (Foc1 = Foc2 = 10, Foc3 = 0.3, Fo4 = 0.1, Fo5 = 1.575) в процессе колебаний наблюдается бифуркационный резонанс, при котором амплитуда колебаний возрастает периодически. При этом в зависимости от соотношения величин Foc1 , Foc2 и Foc3 процесс колебаний может быть как затухающим, так и незатухающим. Отметим, что бифуркационное изменение амплитуды колебаний наблюдается также на частотах колебаний внешней нагрузки, близких к частоте собственных колебаний стержня, но не совпадающих с ней полностью. Результаты расчетов для такого типа колебаний приведены на рис. 17, 18. Обнаруженные в результате исследований бифуркационные (флаттерные) колебания как в резонансных, так и (что особенно важно) в нерезонансных частотах представляют серьезную проблему при работе различных технических устройств и приводят к значительному сокращению их ресурса. Результаты исследований позволяют найти такое соотношение частот собственных и вынужденных колебаний, при которых бифуркационные колебания могут быть исключены, то есть эти результаты позволят разработать рекомендации по установке демпфирующих устройств, способствующих затуханию колебаний, минуя бифуркационный режим. 134 Разработка математических моделей и исследование неравновесных явлений. . . Рис. 14. (онлайн в цвете) Изменение перемещений стержня в резонансных частотах при различных значениях ξ, когда Foc1 = Foc2 = Foc3 = 0, Fo4 = 0.5, Fo5 = 1.575 [Figure 14. (Color online) The change of the rod displacement at the resonant frequencies for various ξ when Foc1 = Foc2 = Foc3 = 0, Fo4 = 0.5, Fo5 = 1.575] Рис. 15. (онлайн в цвете) Изменение перемещений стержня в резонансных частотах при различных значениях ξ, когда Foc1 = Fo2 = 0, Foc3 = 100, Fo4 = 1, Fo5 = 1.575 [Figure 15. (Color online) The change of the rod displacement at the resonant frequencies for various ξ when Foc1 = Fo2 = 0, Foc3 = 100, Fo4 = 1, Fo5 = 1.575] Рис. 16. (онлайн в цвете) Изменение перемещений стержня в резонансных частотах при различных значениях ξ, когда Foc1 = Foc2 = 2, Foc3 = 1, Fo4 = 0.1, Fo5 = 1.575 [Figure 16. (Color online) The change of the rod displacement at the resonant frequencies for various ξ when Foc1 = Foc2 = 2, Foc3 = 1, Fo4 = 0.1, Fo5 = 1.575] 135 К у д и н о в И. В. Рис. 17. (онлайн в цвете) Изменение перемещений стержня в нерезонансных частотах при различных значениях ξ, когда Foc1 = Foc2 = 10, Foc3 = 0.3, Fo4 = 0.1, Fo5 = 0.1 [Figure 17. (Color online) The change of the rod displacement at the nonresonant frequencies for various ξ when Foc1 = Foc2 = 10, Foc3 = 0.3, Fo4 = 0.1, Fo5 = 0.1] Рис. 18. (онлайн в цвете) Изменение перемещений стержня в нерезонансных частотах при различных значениях ξ, когда Foc1 = Fo2 = 10, Foc3 = 0.3, Fo4 = 0.1, Fo5 = 1.5 [Figure 18. (Color online) The change of the rod displacement at the nonresonant frequencies for various ξ when Foc1 = Fo2 = 10, Foc3 = 0.3, Fo4 = 0.1, Fo5 = 1.5] 7. Математическая модель динамической термоупругости с учетом пространственно-временной нелокальности. При больших скоростях изменения температуры температурные напряжения оказываются зависимыми от распределения напряжений в предшествующие моменты времени. Время в данном случае становится независимой переменной и поэтому в уравнениях термоупругости появляются инерционные члены, учитывающие изменение напряжений во времени [34-37]. Классическое одномерное уравнение динамической термоупругости для незакрепленной пластины имеет вид [34-36] ∂ 2 σx (x, t) 1 ∂ 2 σx (x, t) 1 + ν ∂ 2 T (x, t) - = ρα , ∂x2 c2w ∂t2 1-ν ∂t2 t > 0, 0 < x < δ, (36) где c2w = E(1 - ν)/[(1 + ν)(1 - 2ν)ρ], cw - скорость распространения упругой волны, δ - толщина пластины, σx - нормальное напряжение в направлении координаты x, T - температура, T0 - температура в ненапряженном состоянии, t - время, E - модуль упругости, ν - коэффициент Пуассона, α - коэффициент линейного расширения, ρ - плотность. Распределение температуры как функции пространственной переменной и времени предполагается известным из решения соответствующей задачи теплопроводности. 136 Разработка математических моделей и исследование неравновесных явлений. . . Уравнение (35) является волновым, описывающим незатухающие колебания температурных напряжений. Отсутствие затухания связано с тем, что данное уравнение не содержит членов, учитывающих силы внутреннего сопротивления среды процессу изменения ее формы в результате температурных деформаций. Для учета сопротивления среды будем использовать уравнение равновесия (движения) вида ∂σx ∂U ∂2U - γρ =ρ 2 , ∂x ∂t ∂t (37) где второе слагаемое в левой части представляет силу сопротивления, приходящуюся на единицу объема среды, γ - коэффициент сопротивления. Уравнение (35) выведено без учета пространственно-временной локальной неравновесности. Для ее учета при выводе уравнения динамической термоупругости будем использовать формулу закона Гука вида (11). Подставляя (11) в (36) с учетом соотношения между напряжением и деформацией, находим следующее уравнение динамической термоупругости (в безразмерном виде): 3 ¯ (ξ, Fo) ∂2σ ¯ (ξ, Fo) ∂σ ¯ (ξ, Fo) ∂2σ ¯ (ξ, Fo) p ∂ σ + Fo - - Fop3 = 1 2 2 2 ∂ξ ∂ξ ∂ Fo ∂ Fo ∂ Fo ∂ 3 Θ(ξ, Fo) ∂ 2 Θ(ξ, Fo) ∂Θ(ξ, Fo) + , = Fop2 + Fop3 2 2 ∂ξ ∂ Fo ∂ Fo ∂ Fo Fo > 0, 0 < ξ < 1, 0 Θ 1. (38) Здесь Fop1 = c τ1 E τ2 - 2 , 2δ c ρ Fop2 = 1 τ1 E , 2 δcρ Fop3 = 2cργδ (1 + ν)(1 - 2ν) , E 1-ν T - Tw x ct (1 - 2ν)σx , ξ = , Fo = , σ ¯= , ∆T = T0 - Tw ; ∆T 2δ 2δ Eα∆T σ ¯ , Θ, ξ - соответственно безразмерные напряжение, температура, координата; Fo - число Фурье (безразмерное время); T0 - начальная температура; Tw - температура стенки пластины; E - модуль упругости; δ - толщина пластины. Начальные и граничные условия к уравнению (37) для бесконечной пластины при отсутствии внешней нагрузки имеют вид Θ= ∂σ(ξ, Fo) = 0; ∂ Fo Fo=0 σ(0, Fo) = 0, σ(1, Fo) = 0. σ(ξ, 0) = 0, (39) (40) Очевидно, что при Fop1 = Fop2 = Fop3 = 0 задача (37)-(40) приводится к задаче динамической термоупругости без учета пространственно-временной нелокальности и сопротивления среды. В качестве решения температурной задачи использовались приведенные выше точные аналитические решения задачи теплопроводности с учетом пространственно-временной нелокальности. 137 К у д и н о в И. В. Точное аналитическое решение задачи (37)-(40) имеет вид ∞ C1k e-ν1k Fo + C2k e-ν2k Fo - σ ¯ (ξ, Fo) = k=1 - 1 A2 H3 ez1k m Fo + B2 H4 ez2k m Fo H3 H4 cos (2k - 1) (1 - 2ξ)π , (41) 2 где ν1,2 k = 1 Fo4k ± Fo24k -4µk , 2 z1k = -µk = (2k - 1)2 π 2 , z2k = - 1 ; Fop2 C1k , C2k - константы интегрирования, определяемые из начальных условий (39); m, A2 , B2 , H3 , H4 - заданные постоянные; Fo4k = µk Fop1 Fop2 . Результаты расчетов температурных напряжений по формуле (41) приведены на рис. 19-22. В частности, на рис. 19 даны результаты расчетов при Fop1 = Fop2 = Fop3 = 0 для значений ξ = 0.1, 0.3, 0.5. Их анализ позволяет заключить, что напряжения в этих точках с течением времени изменяются скачкообразно, периодически изменяя знак. Волновой процесс изменения напряжений при Fop1 = Fop2 = Fop3 = 0 является незатухающим во времени, что объясняется отсутствием в уравнении (37) членов, отвечающих затуханию колебаний. Отметим, что при Fop3 = 0 в уравнении (37) появляется слагаемое, которое учитывает силы внутреннего трения. На рис. 20 приведены результаты расчетов распределения напряжений при Fop1 = Fop2 = Fop3 = 0 по координате ξ для различных моментов времени. Их анализ позволяет заключить, что изменение напряжений по форме эквивалентно случаю колебаний закрепленной с двух сторон струны, имеющей изломы (в виде прямоугольных скачков напряжений), перемещающиеся по координате ξ во времени. Отметим, что результаты расчетов, приведенные на рис. 19, 20, совпадают с результатами, полученными в работе [37]. Результаты расчетов напряжений в точке ξ = 0.5 для различных значений Fop1 = Fop2 и при Fop3 = 0.1 приведены на рис. 21. Их анализ позволяет заключить, что процесс колебаний является затухающим во времени. При увеличении значений Fop1 = Fop2 (при неизменном Fop3 = 0.1) скачки напряжений сглаживаются и при некоторых больших значениях Fop1 = Fop2 возникающие от теплового удара напряжения монотонно уменьшаются до нулевого значения при практическом отсутствии колебательного процесса. На рис. 22 приведены результаты расчетов напряжений по координате ξ для 0.1 Fo 1.8. Их анализ позволяет заключить, что колебания напряжений по форме аналогичны колебаниям закрепленной на концах струны, не имеющей изломов. Если в качестве решения температурной задачи использовать точное аналитическое решение классического параболического уравнения теплопроводности, то решение динамической задачи термоупругости приводит к бесконечным температурным напряжениям в точках задания граничного условия первого рода (тепловой удар), что связано со следующими обстоятельствами. При выполнении условий теплового удара предполагается мгновенное во времени принятие поверхностью пластины температуры, задаваемой граничным условием первого рода. Так как в динамической задаче термоупругости 138 Разработка математических моделей и исследование неравновесных явлений. . . Рис. 19. (онлайн в цвете) Изменение напряжений в пластине для различных ξ (расчет по формуле (41)), когда Fop1 = = Fop2 = Fop3 = 0, m = 0.001, For = 0.001 - коэффициент релаксации из решения температурной задачи; n = 1000 - число членов ряда [Figure 19. (Color online) The change of stresses in the plate for various ξ (calculation by Eq. (41)) when Fop1 = Fop2 = Fop3 = 0, m = 0.001, For = 0.001, n = 1000; For is the relaxation coefficient from the solution of the temperature problem; n is the number of terms in the series (41)] Рис. 20. (онлайн в цвете) Распределение напряжений в пластине для различных чисел Fo (расчет по формуле (41)), когда Fop1 = Fop2 = Fop3 = 0, m = 0.001, For = 0.001 - коэффициент релаксации из решения температурной задачи; n = 1000 - число членов ряда [Figure 20. (Color online) The stress distribution in the plate for different Fourier numbers (Fo) (calculation by Eq. (41)) when Fop1 = Fop2 = Fop3 = 0, m = 0.001, For = 0.001, n = 1000; For is the relaxation coefficient from the solution of the temperature problem; n is the number of terms in the series (41)] Рис. 21. (онлайн в цвете) Распределение напряжений в точке ξ = 0.5 для различных Fop1 = Fop2 (расчет по формуле (41)), когда Fop3 = 0.1, m = 0.001, For = = 0.001 - коэффициент релаксации из решения температурной задачи; n = 1000 - число членов ряда [Figure 21. (Color online) The stress distribution at the point ξ = 0.5 for different Fop1 = Fop2 (calculation by Eq. (41)) when Fop3 = 0.1, m = 0.001, For = 0.001, n = 1000; For is the relaxation coefficient from the solution of the temperature problem; n is the number of terms in the series (41)] 139 К у д и н о в И. В. Рис. 22. (онлайн в цвете) Распределение напряжений в пластине для различных чисел Fo (расчет по формуле (41)), когда Fop1 = Fop2 = Fop3 = 0.01, m = 0.001, For = 0.001 - коэффициент релаксации из решения температурной задачи; n = 1000 - число членов ряда [Figure 22. (Color online) The stress distribution in the plate for different Fourier numbers (Fo) (calculation by Eq. (41)) when Fop1 = Fop2 = Fop3 = 0.01, m = 0.001, For = 0.001, n = 1000; For is the relaxation coefficient from the solution of the temperature problem; n is the number of terms in the series (41)] учитывается изменение термических напряжений в зависимости от скорости изменения температуры во времени, то мгновенное ее задание означает бесконечную скорость изменения температуры, что и приводит к бесконечным значениям температурных напряжений. Если использовать в качестве решения температурной задачи точное аналитическое решение краевой задачи теплопроводности с учетом пространственно-временномй нелокальности, то напряжения в точке задания граничных условий первого рода имеют конечную величину. Этот факт можно объяснить тем, что температура стенки в данном случае не может быть задана мгновенно, а только лишь в течение некоторого конечного промежутка времени. 8. Исследование локально-неравновесных процессов теплового воспламенения. В теории теплового воспламенения изучается протекание химической реакции в условиях прогрессивного самоускорения, приводящего к неограниченному возрастанию температуры (тепловому воспламенению) [38, 39]. В уравнениях известных математических моделей теплового воспламенения не учитываются релаксационные свойства материалов (локальная неравновесность процесса), что, как показали приведенные ниже исследования, приводит к существенному отклонению критических условий и времени задержки теплового воспламенения по сравнению со случаем их учета. Данный факт можно объяснить тем, что в основе известных моделей лежат параболические уравнения, выведенные из условия бесконечной скорости распространения теплоты. В работах [26-29] приводятся результаты построения моделей, в которых учитывается ускорение во времени как теплового потока, так и градиента температуры. В настоящей работе теоретические положения построения таких моделей были распространены применительно к решению краевой задачи теплового воспламенения для бесконечной пластины с экспоненциально зависящим от температуры внутренним источником теплоты ω(T ) = Qρk0 exp - 140 E RT Разработка математических моделей и исследование неравновесных явлений. . . в следующей математической постановке: ∂Θ ∂2Θ ∂3Θ ∂2Θ γ = + Fo × + Fo1 - 2 2 2 2 ∂ Fo ∂ξ ∂ξ ∂ Fo (Θ∆T + Tw )2 ∂ Fo E ∂Θ E × exp - + τ p exp - , (42) R(Θ∆T + Tw ) ∂ Fo R(Θ∆T + Tw ) Fo > 0, 0 < ξ < 1; ∂Θ(ξ, Fo) Θ(ξ, Fo) = 0; = 0, Θ(1, Fo) = 0, (43) Θ(ξ, 0) = 1, ∂ ∂ Fo Fo=0 ∂ξ ξ=0 где T - Tw , ∆T aτ1 aτ2 , Fo2 = 2 , 2 δ δ 2 δ Qk0 τ1 k0 E τp = , ∆T = T0 - Tw , γ = ; ca∆T cR Θ, Fo, ξ - соответственно безразмерные температура, время, координата; τ p - безразмерный параметр; Fo1 , Fo2 - безразмерные коэффициенты релаксации; T - температура; x - координата; a - коэффициент температуропроводности; t - время; ω(T ) - мощность внутреннего источника теплоты; ρ - плотность; c - теплоемкость; Q - тепловой эффект реакции; k0 - предэкспоненциальный множитель; T0 - начальная температура; Tw - температура стенки; τ1 , τ2 - коэффициенты релаксации теплового потока и градиента температуры в формуле закона Фурье; E - энергия активации; R - универсальная газовая постоянная. Вывод уравнения (42) базируется на формуле закона Фурье (9), в которой учитывается ускорение во времени теплового потока и градиента температуры, и уравнении теплового баланса Θ= Fo = cρ at , δ2 ξ= x , δ Fo1 = ∂q ∂T + + ω(T ) = 0, ∂t ∂x (44) где λ - коэффициент теплопроводности. Для решения задачи (42), (43) применялся метод конечных разностей [40] со следующими исходными данными: τ p = 4.06 · 108 ; E = 127 кДж/моль; R = 8.31 Дж/(моль · К); ∆T = 10 К; Tw = 850 К; k0 = 1.925 · 106 Гц. Результаты расчетов приведены на рис. 23-25. Их анализ позволяет заключить, что величина безразмерных коэффициентов релаксации существенно влияет на время задержки теплового воспламенения. Так, при Fo1 = Fo2 = = 0.15 оно возрастает примерно на 1/3 по сравнению со случаем, когда релаксационные свойства материала не учитываются (Fo1 = Fo2 = 0) (см. рис. 23). Факт задержки теплового воспламенения объясняется релаксационными свойствами материалов. Их учет приводит к невозможности мгновенного принятия граничного условия первого рода на стенке (при ξ = 1) - процесс их установления занимает некоторый начальный диапазон времени 0 Fo Fo∗ . И, в частности, при Fo1 = Fo2 = 0.15 время принятия граничного условия 1-го рода Fo∗ = 0.5 (см. рис. 24). Отсюда следует, что граничные условия первого рода (тепловой удар) не могут быть установлены мгновенно ни при каких условиях теплообмена с окружающей средой. Следовательно, 141 К у д и н о в И. В. Рис. 23. (онлайн в цвете) Изменение температуры Θ, когда ξ = 0.15, γ = 0 [Figure 23. (Color online) The change of the temperature (Θ) when ξ = 0.15, γ = 0] Рис. 24. Распределение температуры Θ для различных чисел Fo, когда Fo1 = Fo2 = 0.15, γ = 0 [Figure 24. The temperature distribution for different Fourier numbers (Fo) when Fo1 = Fo2 = 0.15, γ = 0] Рис. 25. (онлайн в цвете) Изменение температуры Θ, когда ξ = 0.9 [Figure 25. (Color online) The change of the temperature (Θ) when ξ = 0.9] 142 Разработка математических моделей и исследование неравновесных явлений. . . в реальных процессах коэффициент теплоотдачи ограничен некоторой постоянной величиной, определяемой релаксационными свойствами материала. На рис. 25 приведены результаты расчетов в случае когда γ = 0. Их анализ позволяет заключить, что при малых значениях Fo1 = Fo2 величина значений γ оказывает незначительное влияние на время задержки воспламенения. С увеличением значений коэффициентов релаксации Fo1 и Fo2 его влияние значительно возрастает. 9. Расчет критических условий теплового взрыва с учетом релаксационных свойств материалов. Известные математические модели с нелинейными источниками теплоты, допускающие возникновение условий неограниченного возрастания температуры (теплового взрыва), основаны на параболических уравнениях, при выводе которых не учитывается пространственно-временная нелокальность (релаксационные свойства материалов). Такие модели в области малых значений временной и пространственной переменных неадекватно описывают реальные физические процессы, что связано с локальностью описывающих эти процессы дифференциальных уравнений, выведенных исходя из принципа локального термодинамического равновесия и гипотезы сплошной среды [1-6]. Подобные модели допускают бесконечные скорости распространения теплоты и бесконечные значения тепловых потоков (например, при тепловом ударе на поверхности). В работах [7-6, 26-30] приводятся принципы построения моделей, свободных от указанных недостатков, что связано с учетом локальной неравновесности процессов. В настоящей работе указанные принципы распространены на модели с нелинейным внутренним источником теплоты. В качестве конкретного примера рассмотрим задачу теплопроводности для бесконечной пластины при неоднородных граничных условиях третьего рода с линейно зависящим от температуры внутренним источником теплоты ω(T ) = ω0 (1 + βT ) в следующей математической постановке: ∂2Θ ∂2Θ ∂3Θ ∂Θ + Fo1 + Fo + Po1 Θ + Po2 ; = 2 ∂ Fo ∂ξ 2 ∂ξ 2 ∂ Fo ∂ Fo2 ∂Θ(ξ, Fo) Θ(ξ, 0) = 0, = 0; ∂ Fo Fo=0 ∂ 2 Θ(ξ, Fo) ∂Θ(ξ, Fo) + Fo2 - ∂ξ ∂ξ∂ Fo ∂Θ(ξ, Fo) -Bi1 ∆T1 + Θ(ξ, Fo) + Fo1 = 0; ∂ Fo ξ=0 ∂Θ(ξ, Fo) ∂ 2 Θ(ξ, Fo) + Fo2 - ∂ξ ∂ξ∂ Fo ∂Θ(ξ, Fo) -Bi2 ∆T2 - Θ(ξ, Fo) - Fo1 = 0, ∂ Fo ξ=1 (1 + γ) (45) (46) (47) (48) где Θ= T - T0 , T1 - T0 Fo = at , δ2 ξ= x , δ Po1 = ω0 δ 2 β , acρ Po2 = Po1 1 + βT0 , β∆T 143 К у д и н о в И. В. ∆T = T1 - T0 , γ= Fo1 = τ1 ω0 β , cρ aτ1 , δ2 aτ2 α1 δ α2 δ , Bi1 = , Bi2 = , 2 δ λ λ T0 - T1 T2 - T0 ∆T1 = , ∆T2 = ; ∆T ∆T Fo2 = T - температура; T0 - начальная температура; T1 , T2 - температуры сред; α1 , α2 - коэффициенты теплоотдачи; δ - толщина пластины; ρ - плотность; c - теплоемкость; a - коэффициент температуропроводности; ω0 - мощность источника теплоты при начальной температуре; β - коэффициент; τ1 , τ2 - коэффициенты релаксации соответственно теплового потока и градиента температуры в формуле закона Фурье; t - время; x - координата. Вывод уравнения (45) базируется на формуле закона Фурье (9), в которой учитывается ускорение во времени теплового потока и градиента температуры, и уравнении теплового баланса (44). Задача (45)-(48) решалась методом конечных разностей [40]. Анализ полученных результатов (при Bi1 = 5; Bi2 = 10; γ = 0) позволяет заключить, что при β = 0.001 для любых Fo1 и Fo2 неограниченного возрастания температуры (ситуация теплового взрыва) не происходит и при Fo 0.42 наблюдается установившееся (стационарное) состояние, при котором теплоприток от источника теплоты компенсируется теплоотводом в окружающую среду (см. рис. 26 и кривую 3 на рис. 27). С увеличением β (например при β = 0.1) тепловой взрыв возникает при любых Fo1 и Fo2 , причем время его задержки существенно зависит от их величины. Так, при увеличении Fo1 = Fo2 от их значения, равного 10-3 , до величины 0.05 время задержки теплового взрыва в точке ξ = 0.9 возрастает почти в два раза (см. кривые 1, 2 на рис. 27). Причиной его задержки является тот факт, что при учете релаксационных свойств материала тепловые потоки на границах в некотором диапазоне времени возрастают от нуля до значения, определяемого граничными условиями третьего рода. В силу учета релаксационных свойств материала граничные условия (независимо от их вида - первого, второго, третьего) могут быть выполнены лишь по истечении некоторого начального диапазона времени 0 Fo Fo∗ , где Fo∗ - время, при котором начинает выполняться заданное граничное условие краевой задачи. Следовательно, в определенном диапазоне начального временного участка прогрев пластины происходит в условиях ограниченного теплообмена с окружающей средой, участвующей в прогреве конструкции (при T1 > T0 , T2 > T0 ), что и является сдерживающим фактором процесса теплового взрыва. Было выполнено также исследование влияния релаксационных параметров на время задержки теплового взрыва в случае полного отсутствия теплообмена с окружающей средой, то есть, когда температура тела определяется действием лишь внутреннего источника теплоты. Полученные результаты позволяют заключить, что время задержки теплового взрыва в данном случае существенно увеличивается, что объясняется учетом ускорения во времени теплового потока и градиента температуры. На рис. 28 приведены результаты расчетов изменения температуры во времени в точке ξ = 0.9 при γ = 0 и γ = 2.5 (β = 0.1, Fo1 = Fo2 = 0.1, Bi1 = 5, Bi2 = 10). Их анализ позволяет заключить, что величина γ оказывает определенное влияние на время задержки теплового взрыва, которое возрастает с увеличением γ. 144 Разработка математических моделей и исследование неравновесных явлений. . . Рис. 26. Распределение температуры Θ для различных чисел Fo, когда β = 0.001, Fo1 = Fo2 = 0.1 [Figure 26. The temperature distribution for different Fourier numbers (Fo) when β = 0.001, Fo1 = Fo2 = 0.1] Рис. 27. Изменение температуры Θ, когда ξ = 0.9 и β = 0.1, Fo1 = Fo2 = 10-3 (линия 1), β = = 0.1, Fo1 = Fo2 = 0.05 (линия 2), β = 0.001, Fo1 = Fo2 = 0.1 (линия 3) [Figure 27. The change temperature (Θ) when ξ and β = 0.1, Fo1 = 10-3 (curve 1), and β Fo1 = Fo2 = 0.05 (curve β = 0.001, Fo1 = Fo2 (curve 3)] of the = 0.9, Fo2 = = 0.1, 2), and = 0.1 Рис. 28. Изменение температуры Θ, когда ξ = 0.9, Fo1 = Fo2 = = 0.1, β = 0.1 и γ = 0 (линия 1), γ = 2.5 (линия 2) [Figure 28. The change of the temperature (Θ) when ξ = 0.9, Fo1 = Fo2 = 0.1, and γ = 0 (curve 1), and γ = 2.5 (curve 2)] 10. Экспериментальные исследования продольных колебаний стержня. Экспериментальные исследования продольных колебаний стержня были выполнены на специальном стенде. В процессе испытаний закрепленный на одном из торцов стержень длиной 30 см растягивался на 0.1 мм. После сброса нагрузки с помощью специальной измерительной системы фиксировались колебания свободного торца стержня во времени. Результаты исследований приведены на рис. 29. Их анализ позволяет заключить, что концевое сечение стержня участвует в двух колебательных процессах - с малой амплитудой и большой частотой (временной участок a на рис. 29) и одновременно - с большой амплитудой и малой частотой (временной участок b на рис. 29). 145 К у д и н о в И. В. Рис. 29. Перемещения свободного торца стержня [Figure 29. The displacements of free end of the rod] Полученные экспериментальные данные качественно полностью согласуются с результатами теоретических исследований по разработанным моделям колебаний упругих тел с учетом пространственно-временной нелокальности. Количественное сравнение экспериментальных и теоретических исследований показывает, что при Fo1 = Fo2 = 10 расчетная частота колебаний составляет 2 380 Гц, а экспериментальная - 2 400 Гц. Расчетная амплитуда колебаний в диапазоне от 0 до 0.01 с равна 0.0033 мм, средняя экспериментальная амплитуда в этом же диапазоне времени составляет 0.0028 мм. Выводы. Сформулируем основные выводы по работе. 1. С учетом пространственно-временн´oй нелокальности разработана математическая модель теплопереноса, основанная на представлении классической формулы закона Фурье для теплового потока в ряд по степеням малых параметров коэффициентов релаксации теплового потока и скалярной величины градиента температуры. На основе этой формулы и уравнения теплового баланса получено дифференциальное уравнение локально-неравновесного теплообмена, включающее производные высокого порядка, в том числе и смешанные производные. 2. Разработаны методы получения точных аналитических решений краевых задач локально-неравновесного теплообмена, основанные на использовании метода разделения переменных и операционных методов. Исследование полученных решений позволило сделать вывод об отсутствии скачков температуры внутри рассматриваемых областей, а также отрицательных значений температур в обратной тепловой волне. Скачки температур, наблюдающиеся лишь в бесконечно малой окрестности граничных точек на начальном участке временн´oй переменной, свидетельствуют о физической невозможности мгновенного задания граничных условий первого рода. 3. С учетом пространственно-временн´oй нелокальности разработана математическая модель нестационарного неравновесного теплообмена при ламинарном течении жидкости в цилиндрической трубе. Анализ численного решения задачи (метод прогонки) позволил сделать заключение об отсутствии скачков температуры внутри движущейся жидкости, а также отрицательных ее значений во всем диапазоне времени нестационарного процесса. Выполненные исследования показали невозможность мгновенного охлаждения жидкости на стенке - температура жидкости принимает температуру стенки за некоторый конечный интервал 146 Разработка математических моделей и исследование неравновесных явлений. . . времени. 4. С учетом пространственно-временн´oй нелокальности разработана математическая модель затухающих колебаний упругих тел (стержней, пружин, упругих жидкостей и проч.), основанная на представлении классической формулы закона Гука для нормального напряжения по степеням малых параметров коэффициентов релаксации напряжения и деформации (скалярной величины градиента скорости). На основе этой формулы и уравнения равновесия (уравнение второго закона Ньютона) получено дифференциальное уравнение локально-неравновесных колебаний упругих тел. Исследование его точных аналитических решений позволило заключить, что использование релаксационных составляющих напряжений и деформаций в формуле Гука приводит к устранению скачкообразного изменения напряжений в процессе колебаний среды, выведенной из состояния равновесия. 5. Результаты теоретических исследований колебаний закрепленного на одном из торцов стержня сравнивались с данными натурных экспериментов. Стержень подвергался растяжению на 0.1 мм и после сброса нагрузки с помощью лазерного интерферометра определялись его продольные колебания во времени. Результаты экспериментальных исследований подтвердили результаты точного аналитического решения, согласно которому каждая точка стержня участвует в двух колебательных процессах. В одном из них совершаются высокочастотные и низкоамплитудные колебания, а в другом - низкочастотные и высокоамплитудные. На основе полученного точного аналитического решения и экспериментальных данных путем решения обратной задачи были найдены коэффициенты релаксации напряжения и градиента перемещения. 6. Сравнение результатов теоретических исследований с экспериментальными данными было выполнено также и для жидкости (масло) в условиях гидравлического удара. С использованием точного аналитического решения и результатов экспериментов, выполненных Е. Л. Холмбоу и В. Т. Руло [32], путем решения обратной задачи были найдены коэффициенты релаксации касательного напряжения и градиента перемещения, которые оказались равными τ1 = 0.0108 c; τ2 = 0.0116 c. Удовлетворительное совпадение результатов теоретических и экспериментальных исследований при найденных таким путем коэффициентах релаксации как для стержня, так и для несжимаемой жидкости свидетельствует о физической правомерности представления формул законов Гука и Ньютона с учетом ускорения во времени напряжений (нормального и касательного) и градиентов соответствующих величин (перемещений и скорости). 7. С учетом пространственно-временн´oй нелокальности разработана математическая модель динамической термоупругости, основанная на разложении классической формулы закона Гука в ряд по степеням малых параметров коэффициентов релаксации напряжений и деформаций при использовании локально-неравновесной модели теплопереноса. На основе метода разделения переменных получено точное аналитическое решение краевой задачи локально-неравновесной динамической термоупругости для свободной бесконечной пластины при симметрич147 К у д и н о в И. В. ных граничных условиях первого рода температурной задачи. Конкурирующие интересы. У меня нет конкурирующих интересов. Авторская ответственность. Я несу полную ответственность за предоставление окончательной версии рукописи в печать. Окончательная версия рукописи мною одобрена.

About the authors

Igor V Kudinov

Samara State Technical University

Email: igor-kudinov@bk.ru
Cand. Techn. Sci.; Associate Professor; Dept. of Theoretical Fundamentals of Heat-Engineering and Hydromechanics 244, Molodogvardeyskaya st., Samara, 443100, Russian Federation

References

Supplementary files