Невязкий аналог задачи Пуазейля

Полный текст

1. Введение Задача о движении жидкости в канале является одной из важных задач теоретической гидромеханики. Эта задача имеет не только теоретическое, но и прикладное значение, так как в основе лежит важный практический вопрос - каков закон движения жидкости в канале при заданных условиях на границах [1-3]. В прикладной гидромеханике движения такого типа получили названия напорных течений. Необходимость рассмотрения такого типа течений возникает при расчётах агрегатов и механизмов на трубопроводном транспорте, в машиностроении, энергетике, гидротехнике. На сегодняшний день многие вопросы, связанные с задачами такого типа, прояснены не до конца и требуют дополнительного изучения. Например, не ясно, по какому закону происходит падение давления вдоль оси канала, каков характер рассеяния энергии потока, при каких условиях одномерное движение разрушается и заменяется на более сложное двумерное с наличием выраженных вихревых зон, каковы закономерности перехода от ламинарного режима движения к турбулентному [4, 5]. Простейшая постановка предполагает рассмотрение плоской задачи и основывается на модели идеальной несжимаемой жидкости, движение которой описывается уравнениями Эйлера. Для случая 2D установившегося движения и в предположении об отсутствии внешних массовых сил, уравнения Эйлера имеют вид [1-3] u¶u ¶x + v¶u ¶y = -1 r ¶p ¶x, (1) u¶v ¶x + v¶v ¶y = -1 r ¶p ¶y, (2) ¶u ¶x + ¶v ¶y = 0, (3) где p, u, v обозначают основные неизвестные - давление и компоненты вектора скорости, соответственно продольной и поперечной; r есть плотность жидкости, которая в рамках модели несжимаемой среды неизменна и представляет заданный положительный параметр. Каждая из величин p, u, v, является некоторой неизвестной функцией координат x, y. Основная задача состоит в определении этих неизвестных. Важнейшей особенностью уравнений (??)-(??) является наличие нелинейных членов, которые присутствуют в левых частях первых двух уравнений. Эти члены уравнений создают основные сложности при исследовании и решении. Но именно эти члены являются неотъемлемой частью уравнений Эйлера и должны быть учтены при решении по возможности более строго. Ранее было обнаружено, что уравнения (??)-(??) допускают простое одномерное решение. Действительно, если предположить, что течение имеет лишь одну ненулевую компоненту скорости вдоль оси OX, тогда v = 0,¶u ¶x = 0,¶p ¶x = 0,¶p ¶y = 0. (4) В результате получаем p = const и u = u0(y). То есть давление не изменяется и профиль продольной скорости, который был во входном сечении канала, будет сохраняться без каких-либо изменений и при движении вдоль всего канала. Решение (??), которое соответствует такому простейшему движению, можно назвать тривиальным. Однако экспериментальные результаты и простые визуальные наблюдения говорят о том, что, как правило, существуют зоны течения, где движение заведомо не одномерно и v№0. Кроме того, в реальных течениях всегда существует некоторый ненулевой перепад давления, который является основной причиной движения. Все это заставляет вновь обращаться к этой известной задаче теоретической гидромеханики и сосредоточиться на поиске нетривиальных решений, которые соответствуют наблюдаемым на практике заведомо не одномерным движениям. 2. Постановка задачи Будем считать, что движение жидкости происходит в канале между двумя параллельными плоскостями, расположенными на расстоянии 2H друг от друга. Выберем систему координат так, что ось OY перпендикулярна ограничивающим поверхностям, а ось OX им параллельна и отстоит от каждой из поверхностей на расстояние H. Ограничимся рассмотрением плоской задачи, считая, что движение одинаково во всех плоскостях, перпендикулярных ограничивающим поверхностям и параллельных оси OX. Считаем, что канал имеет конечную протяжённость вдоль оси OX. Пусть для простоты - L < x < L, где величину L полагаем положительной и изначально заданной. Будем считать, что на концах канала при x = ±L создаётся перепад давления Dp, который и является основной причиной движения. Таким образом, постановка задачи аналогична известной задаче Пуазейля для движения вязкой жидкости в канале [1-4] с тем лишь различием, что вместо вязкой жидкости в нашем случае рассматривается идеальная и канал представляет не цилиндрическую трубу, как обычно предполагается в задаче Пуазейля, а имеет более простую конфигурацию. Основной вопрос состоит в том, чтобы определить гидродинамические характеристики потока при заданном перепаде давления на концах канала. В данной работе предлагается подход, основная идея которого состоит в том, чтобы строить решение задачи основываясь не на уравнениях Эйлера (??)-(??) непосредственно, а на основе первого интеграла этих уравнений. Первый интеграл уравнений (??)-(??) есть следствие интеграла уравнений Навье - Стокса, полученного автором для общего случая 3D неустановившегося движения вязкой несжимаемой жидкости [6-8] при сохранении нелинейных членов. Для рассматриваемого в данной работе более простого случая 2D движения при 1 Re = 0 и ¶ ¶t = 0 вывод соотношений, представляющих первый интеграл, сводится к следующему. Преобразуем нелинейные члены уравнений (??)-(??) с учётом известных формул дифференцирования произведения и уравнения неразрывности (??) как v¶u ¶y = ¶uv ¶y - u¶v ¶y = ¶uv ¶y + ¶ ¶x u2 2 , u¶v ¶x = ¶uv ¶x - v¶u ¶x = ¶uv ¶x + ¶ ¶y v2 2 . (5) С учётом выражений для нелинейных членов в виде (??) каждое из уравнений системы (??)-(??) представляется в дивергентном виде ¶Pi ¶x + ¶Qi ¶y = 0, (6) где Pi и Qi есть некоторые комбинации неизвестных. Уравнение неразрывности (??) изначально имеет вид (??), а для уравнений (??)-(??) нужные представления следующие: ¶ ¶x u2 + p r + ¶uv ¶y = 0, ¶uv ¶x + ¶ ¶y v2 + p r = 0. (7) Известно, что решение уравнения неразрывности (??) задаётся выражениями [1-3] u = ¶Y1 ¶y ,v = -¶Y1 ¶x , (8) где Y1(x,y) представляет функцию тока, хорошо известную в гидромеханике величину. Аналогично можно построить решение всякого уравнения вида (??). В частности, решения уравнений (??) можно представить в виде u2 + p r = ¶Y1,2 ¶y ,uv = -¶Y1,2 ¶x , uv = ¶Y1,3 ¶y ,v2 + p r = -¶Y1,3 ¶x , (9) где Y1,2 и Y1,3 представляют некоторые функции переменных x и y. Обозначения этих функций соответствуют названиям, которые были предложены в [6-8]. В указанных работах эти функции были названы псевдофункциями тока первого порядка с номерами 2 и 3 соответственно. Из уравнений (??) следуют соотношения p r = -u2 + v2 2 + 1 2 ¶Y1,2 ¶y -¶Y1,3 ¶x ,u2 - v2 = ¶Y1,2 ¶y + ¶Y1,3 ¶x , ¶Y1,3 ¶y + ¶Y1,2 ¶x = 0. (10) Последнее из уравнений (??) представляет уравнение вида (??), и его решение можно представить как Y1,2 = ¶Y2,1 ¶y ,Y1,3 = -¶Y2,1 ¶x , (11) где Y2,1 есть псевдофункция тока второго порядка с номером 1. Она представляет некоторую функцию переменных x и y. Вводя эту функцию в качестве нового ассоциированного неизвестного и используя для неё более простое обозначение Y2, из (??) получаем определяющие соотношения в виде p - p0 r = -U2 2 - d, (12) u2 - v2 = ¶2Y 2 ¶y2 -¶2Y2 ¶x2 , (13) uv = -¶2Y2 ¶x¶y. (14) В уравнениях (??)-(??) p0 есть аддитивная постоянная давления; U = u2 + v2 - модуль вектора скорости; Y2 - новое ассоциированное неизвестное, возникающее при интегрировании; величина d в правой части (??) есть приведённая плотность вихревой энергии, вычисляемая по формуле d = -1 2 ¶2Y2 ¶x2 + ¶2Y2 ¶y2 . (15) Если выполнены (??)-(??), то исходные уравнения Эйлера (??)-(??) также будут удовлетворены. Соотношения (??)-(??) представляют первый интеграл уравнений Эйлера для случая 2D установившегося движения идеальной несжимаемой среды. Уравнения (??)-(??) имеют нулевой порядок производных относительно основных неизвестных u, v, p и значит, решение этих уравнений сводится к более простой математической задаче. Именно уравнения (??)-(??) вместе с уравнением неразрывности (??) возьмём в качестве основных для дальнейшего решения задачи. Задачу будем рассматривать в безразмерных переменных при следующем выборе масштабов. В нашем случае специфика задачи подсказывает выбрать для координат x и y разные масштабы. В качестве масштаба координаты y возьмём H, в качестве масштаба координаты x - величину L, масштабом продольной скорости выберем величину U0, а масштабом поперечной скорости - величину kU0, где k = H L . В качестве масштаба давления полагаем rU02. Для всех безразмерных величин сохраним те же обозначения, что и для соответствующих им размерным величинам. С учётом указанного выбора масштабов область течения в нашем случае будет задаваться неравенствами - 1 < x < 1 и - 1 < y < 1. Кроме того, некоторые из определяющих уравнений несколько изменят свой вид. Так с учётом выбора масштабов уравнения (??) и (??) будут иметь вид p - p0 = -U2 2 - d, (16) u2 - k2v2 = ¶2Y 2 ¶y2 - k2¶2Y 2 ¶x2 . (17) По той же причине выражение для модуля скорости примет вид U = u2 + k2 v2, а уравнение (??) должно быть записано как d = -1 2 k2¶2Y 2 ¶x2 + ¶2Y2 ¶y2 . (18) Граничные условия зададим следующим образом. На ограничивающих поверхностях считаем выполненными условия непротекания [2-4]. Для нашего случая они сводятся к равенству нулю поперечной составляющей скорости вдоль границ канала при y = ±1. Таким образом, потребуем выполнения равенств v(x,±1) = 0. (19) Будем считать, что на входе в канал в сечении x = -1 скорость имеет только продольную составляющую, тогда как поперечная составляющая равна нулю: v(-1,y) = 0. (20) Считаем также, что на концах канала при x = ±1 перепад давления задан. Пусть в безразмерных переменных его значение равно величине θ, где θ представляет заданный положительный параметр, θ = Dp∕(rU02). Тогда данное данное граничное условие задаётся равенством p(-1,y) - p(1,y) = θ. (21) Таким образом, в рамках рассматриваемой постановки задача сводится к решению уравнений (??), (??)-(??), (??) при трёх граничных условиях (??)-(??). Величины k и θ представляют изначально заданные положительные параметры. 3. Решение определяющих уравнений В уравнениях (??), (??)-(??) отсутствуют производные основных неизвестных u, v, p. То есть по отношению к этим неизвестных имеем чисто алгебраические уравнения. Предлагается подход к решению, который был использован автором в работах [9, 10]. На первом этапе зададим общий вид выражений для u(x,y) и v(x,y) так, чтобы заведомо удовлетворить и уравнению неразрывности (??), и граничным условиям (??)-(??). На втором этапе уточним полученные выражения, чтобы удовлетворить уравнениям (??), (??). На последних этапах составим выражение для p - p0 в соответствии с равенством (??) и удовлетворим граничному условию (??). Таким образом, всё последующее решение можно представить, как решение отдельных более простых задач. 3.1. Первый этап Чтобы изначально удовлетворить уравнению неразрывности (??), введём функцию тока Y1. Тогда для неизвестных u и v справедливы соотношения (??). Функцию тока зададим в виде разложения по степеням xnym, где показатели степеней ограничены неравенством 0 Ј n + m Ј N. Таким образом, положим Y1(x,y) = е n=0N е m=0N-na nmxnym, (22) где anm представляют некоторые, пока не определённые коэффициенты. Предварительный анализ показывает, что нетривиальные решения получаются уже, начиная с пятого приближения N і 5. В дальнейшем полагаем N = 5 и будем искать неизвестное Y1(x,y) в виде многочлена пятой степени относительно xnym. Без учёта аддитивной постоянной a00, значение которой произвольно, на данном этапе имеем 20 неопределённых коэффициентов. При этом u(x,y) и v(x,y),, вычисленные согласно (??), представляют многочлены четвёртой степени относительно xnym, и они заведомо удовлетворяют уравнению неразрывности (??). 3.2. Второй этап Подберём коэффициенты anm так, чтобы удовлетворить граничному условию (??). Рассмотрим выражение для поперечной скорости, согласно второму из (??) и с учётом (??). Получаем выражение v = -(a10 + 2a20x + a11y + 3a30x2 + 2a 21xy + a12y2 + 4a 40x3 + 3a 31x2y+ 2a22xy2 + a 13y3 + 5a 50x4 + 4a 41x3y + 3a 32x2y2 + 2a 23xy3 + a 14y4).(23) Условие (??) требует, чтобы при y = ±1 поперечная скорость v обратилась бы в нуль. А значит, в рамках рассматриваемого приближения должно быть выполнено v = (1 - y2)(h 00 + h10x + h01y + h20x2 + h 11xy + h02y2), (24) где hij некоторые вспомогательные коэффициенты. Составляя равенства правых частей (??) и (??), и приравнивая коэффициенты при одинаковых степенях xnym, получаем выражения для вспомогательных коэффициентов hij через основные коэффициенты anm h00 = -a10,h10 = -2a20,h01 = -a11,h20 = -3a30, h11 = -2a21,h02 = -a12 - a10. Кроме того, в результате указанной процедуры получаются и дополнительные ограничения на коэффициенты anm a40 = 0,a31 = 0,a22 = -a20,a13 = -a11,a23 = -a21, a50 = 0,a41 = 0,a32 = -a30,a14 = -a12 - a10. (25) Последние равенства гарантируют выполнимость граничного условия (??). 3.3. Третий этап Наложим новые ограничения на anm, чтобы удовлетворить граничному условию (??). Это условие требует, чтобы при x = -1 поперечная скорость v была нулевая. Обратимся к выражению (??) с учётом (??), и положим x = -1. В результате получаем v = -((a10 - 2a20 + 3a30) + y(a11 - 2a21) + y2(a 12 + 2a20 - 3a30)+ y3(a 13 + 2a21) - y4(a 12 + a10)). Правая часть обращается в нуль при условии, что все коэффициенты при степенях ym будут нулевыми. Рассматривая полученные равенства как систему уравнений, получаем три новых ограничения на коэффициенты a30 = 1 3(2a20 - a10),a21 = 1 2a11,a12 = -a10. (26) Таким образом, из двадцати коэффициентов, которые присутствовали в выражениях изначально, уже определено двенадцать. На данном этапе имеем выражения для скоростей, которые удовлетворяют и уравнению неразрывности, и граничным условиям (??)-(??). Эти выражения следующие u = a01 + a11x + 2a02y + 1 2a11x2 - 2a 10xy + 3a03y2 - 2a 20x2y - 3a 11xy2+ + 4a04y3 + 2 3(a10 - 2a20)x3y -3 2a11x2y2 + 5a 05y4, (27) v = -a10 - 2a20x - a11y + (a10 - 2a20)x2 - a 11xy + a10y2 + 2a 20xy2+ + a11y3 + (2a 20 - a10)x2y2 + a 11xy3.(28) Введём ещё одно предположение. Будем считать, что в правой части выражения (??) для u первое слагаемое равно единице, то есть a01 = 1. Это предположение соответствует тому, что в размерных единицах первое слагаемое совпадает с величиной U0, выбранной в качестве масштаба продольной скорости. 3.4. Четвёртый этап На этом этапе нужно обеспечить выполнимость уравнений (??), (??), которые связывают основные неизвестные u, v и производные ассоциированного неизвестного Y2. Ассоциированное неизвестное Y2(x,y), так же как и функцию тока Y1(x,y), представляем в виде разложения по степеням Y2(x,y) = е n=0N е m=0N-nb nmxnym, (29) где bnm есть некоторые, пока не определённые коэффициенты. Считаем, что N = 5 и ограничиваемся рассмотрением в правой части (??) степеней не выше пятой 0 Ј n + m Ј 5. Неизвестное Y2(x,y) полностью определено, если все коэффициенты bnm будут найдены. Обратимся к уравнениям (??), (??) и найдём разложения всех членов уравнений по степеням xnym. Сравнивая коэффициенты при одинаковых степенях xnym в левых и правых частях, получаем неоднородные системы уравнений относительно неизвестных bnm. Для n + m = 2 и n + m = 3 имеем решения только при условии выполнимости условий совместности. Эти условия задаются следующими тремя уравнениями 6a03a10 - 4k2a 20 + 2(1 + k2)a 10 + k2a 11a10 = 0, 6a03a20 + 2a20 - (2 - k2)a 11a10 - k2a 20a11 = 0, (30) 12a04a10 + 2(1 + k2)a 02a10 - 4k2a 20a02 + 3a03a11 + 3a11 - 2k2a 20a10 + k2 2 a112 = 0. Если выполнены соотношения (??), то коэффициенты bnm находятся однозначно за исключением b20, который может быть выбран произвольно. Вычисления приводят к следующим выражениям b02 = k2b 20 + 1 2 1 - k2a 102 ,b 11 = a10, b30 = 1 3 1 k2a02a10 - 1 2k2a11 + 2a20a10 ,b21 = a20 + 1 2a11a10, b12 = a02a10 + 1 2a11,b03 = 1 3 2a02 -k2 2 a11a10 + k2a 20 , b40 = 1 3a202 + 1 3a20a10 -1 + 2k2 12k2 a102 + 1 6k2a20a02 - 1 24k2 a112 + a 11 , b31 = 1 3 2a20 + 2a20a11 - a10 + 1 2a11a10 , b22 = 1 4 4a20a02 - 2a102 + a 11 + a112 ,b 13 = 1 3 2a11a02 - a10 + 3a03a10 , b04 = 1 2 a03 + 2 3a022 + k2 3 a20a02 + k2 6 a102 -k2 12a112 + k2 12a11 , b50 = 2 5a202 - 1 + 2k2 10k2 a20a10 + 1 15k2a20a02 - 1 30k2a02a10 - 1 40k2a112, b41 = 3 4a20a11 - 1 4a11a10,b32 = 2 3a20a02 - 1 3a02a10 - a20a10 + 1 4a112, b23 = a03a20 - a11a10 + 1 3 a02a11 - a20 , b14 = a04a10 + 3 4a03a11 - 1 2a02a10 - 1 4a11, b05 = k2 10 a03a20 + a11a10 + 1 5 3a03a02 + 2a04 + k2 30 a02a11 - a20 . (31) Выражения (??) полностью определяют неизвестное Y2(x,y) в рассматриваемом приближении. Уравнения (??), (??) также удовлетворены. 3.5. Пятый этап Осталось удовлетворить граничному условию (??). Рассмотрим правую часть уравнения (??). В правой части присутствуют вторые производные неизвестного Y2, поскольку через них определяется величина d согласно (??). В рассматриваемом приближении неизвестное Y2 определялось с точностью до пятых степеней, а значит, для неизвестного p должны получить разложение с точностью до третьих степеней. То есть неизвестное p будет определено равенством: p = p00 + p10x + p01y + p20x2 + p 11xy + p02y2 + p 30x3 + p 21x2y+ + p30x3 + p 12xy2 + p 03y3, (32) где p00 является произвольно выбираемой аддитивной постоянной давления, а остальные коэффициенты pnm должны быть определены. Рассматривая разложения по степеням для всех членов правой части (??), приходим к следующим выражениям p10 = 2a02a10 - a11,p01 = 2k2a 20 - k2a 11a10, p20 = 2a20a02 - a102 -1 2 a11 + a112 , p11 = 2k2a 20 + a10 1 - k2 + 3a 03 -k2 2 a11 , p02 = k2 a 102 + 2a 20a02 + 1 2 a11 - a112 , p30 = -1 2a112 - 2a 20a10 + 2 3a02 2a20 - a10 , (33) p21 = a20 1 + 3a03 + k2 2 a11 -1 + k2 2 a11a10, p12 = 3a11 2 1 + a03 + a02a10 1 - k2 + 6a 04a10 + k2 2a 20a02 + a20a10 -1 4a112 , p03 = 2k2 a 11a10 + a03a20 + 1 3a02a11 -1 3a20 . Условие (??) в нашем случае сводится к равенству p10 + p30 + p11y + p12y2 = -θ 2. (34) Это равенство должно выполняться тождественно при любых - 1 < y < 1. Для этого необходимо и достаточно p10 + p30 = -θ 2,p11 = 0,p12 = 0. (35) Выражения (??) позволяют уточнить левые части (??). В результате несложных вычислений приходим к трём уравнениям 4 3 a02a10 + a20a02 - a11 - 2a20a10 - 1 2a112 = -θ 2, 3a03a10 + 2k2a 20 - k2 2 a11a10 + 1 - k2 a 10 = 0, 6a04a10 + a02a10 + 3 2a11 a03 + 1 + k2 2a 20a02 + a20a10 - a02a10 -1 4a112 = 0. (36) Уравнения (??) вместе с уравнениями (??), определяющими условия совместности, представляют дополнительные ограничения на оставшиеся коэффициенты aij. Если они выполнены, то все уравнения и граничные условия удовлетворены и задача полностью решена. Уравнения (??), (??) образуют систему шести нелинейных уравнений с шестью неизвестными a10, a20, a11, a02, a03, a04. Для завершения решения задачи осталось разрешить эту систему и найти указанные определяющие коэффициенты. Все остальные интересующие нас величины, включая и основные гидродинамические характеристики u(x,y), v(x,y), p(x,y), могут быть найдены через определяющие коэффициенты по приведённым выше формулам. Для u(x,y), v(x,y) такими формулами являются (??), (??), а для p(x,y) это формулы (??) с учётом (??). Аналитическое решение системы (??), (??) затруднительно, поэтому решение было реализовано численно с помощью программ стандартного пакета Maple. 4. Результаты Решение было реализовано для случая θ = 0,1 и k = 0,5. Выбранные значения параметров соответствуют случаю сравнительно небольшого перепада давления и относительно малой протяжённости канала, для которого L = 2H. Численное решение системы приводит к следующим результатам. Значения шести определяющих коэффициентов получились как a10 = ±4,113;a20 = ∓28,789;a11 = -30,00; a02 = ∓5,567;a03 = -0,333;a04 = ±2,144. Имеем два набора коэффициентов, а значит, имеем два различных решения. Для обоих решений коэффициенты a11 и a03 совпадают, а значения остальных отличаются знаком. Ниже представлены таблицы значений основных гидродинамических характеристик, как функции координат x и y, которые соответствуют каждому из этих решений. Структура таблиц следующая. Каждая вертикальная колонка таблицы соответствует определённому значению координаты x, начиная с x = -0,8. Горизонтальные строки таблиц соответствуют значениям y, начиная с y = -0,9. Первыми представлены характеристики для решения 1, для которого выбирались верхние знаки коэффициентов aij. В обозначении характеристик для решения 1 и решения 2 добавляется соответствующий верхний индекс. Таблица 1 Решение 1: Профиль продольной скорости u(1)(x,y) y \ x - 0,8 - 0,4 0,0 0,4 0,8 - 0,9 - 36,77 - 18,65 3,96 16,85 5,81 - 0,6 - 9,11 - 1,22 5,47 1,49 - 22,64 - 0,3 7,82 8,07 4,02 - 9,06 - 35,91 0,0 15,40 10,60 1,00 - 13,40 - 32,60 0,3 15,03 7,77 - 2,20 - 10,14 - 11,33 0,6 8,08 0,97 - 4,19 2,10 29,30 0,9 - 4,04 - 8,42 - 3,58 24,71 90,66 Таблица 2 Решение 1: Профиль поперечной скорости v(1)(x,y) y \ x - 0,8 - 0,4 0,0 0,4 0,8 - 0,9 - 3,06 - 6,36 - 5,91 - 1,71 6,24 - 0,6 - 9,15 - 17,97 - 14,15 2,30 31,38 - 0,3 - 11,37 - 20,63 - 11,93 14,73 59,36 0,0 - 10,69 - 17,28 - 4,11 28,79 81,43 0,3 - 8,09 - 10,81 4,45 37,66 88,85 0,6 - 4,54 - 4,14 8,88 34,55 72,85 0,9 - 1,01 - 0,20 4,35 12,65 24,71 Таблица 3 Решение 1: Давление p(1)(x,y) -p 0, как функция координат y \ x - 0,8 - 0,4 0,0 0,4 0,8 - 0,9 - 244,63 - 74,21 - 22,96 - 84,82 - 253,76 - 0,6 - 190,41 - 55,52 - 16,11 - 66,13 - 199,54 - 0,3 - 134,06 - 34,71 - 7,15 - 45,33 - 143,20 0,0 - 79,53 - 15,72 0,00 - 26,33 - 88,67 0,3 - 30,74 - 2,46 1,41 - 13,08 - 39,88 0,6 8,37 1,12 - 6,85 - 9,49 - 0,76 0,9 33,89 - 8,90 - 28,71 - 19,51 24,76 Таблица 4 Решение 2: Профиль продольной скорости u(2)(x,y) y \ x - 0,8 - 0,4 0,0 0,4 0,8 - 0,9 - 4,04 - 8,42 - 3,58 24,71 90,66 - 0,6 8,08 0,97 - 4,19 2,10 29,30 - 0,3 15,03 7,77 - 2,20 - 10,14 - 11,33 0,0 15,40 10,60 1,00 - 13,40 - 32,60 0,3 7,82 8,07 4,02 - 9,06 - 35,91 0,6 - 9,11 - 1,22 5,47 1,49 - 22,64 0,9 - 36,77 - 18,65 3,96 16,85 5,81 Таблица 5 Решение 2: Профиль поперечной скорости v(2)(x,y) y \ x - 0,8 - 0,4 0,0 0,4 0,8 - 0,9 1,01 0,20 - 4,35 - 12,65 - 24,71 - 0,6 4,54 4,14 - 8,88 - 34,55 - 72,85 - 0,3 8,09 10,81 - 4,45 - 37,66 - 88,85 0,0 10,69 17,28 4,11 - 28,79 - 81,43 0,3 11,37 20,63 11,93 - 14,73 - 59,36 0,6 9,15 17,97 14,15 - 2,30 - 31,38 0,9 3,06 6,36 5,91 1,71 - 6,24 Таблица 6 Решение 2: Давление p(2)(x,y) -p 0, как функция координат y \ x - 0,8 - 0,4 0,0 0,4 0,8 - 0,9 33,89 - 8,90 - 28,71 - 19,51 24,76 - 0,6 8,37 1,12 - 6,85 - 9,49 - 0,76 - 0,3 - 30,74 - 2,46 1,41 - 13,08 - 39,88 0,0 - 79,53 - 15,72 0,00 - 26,33 - 88,67 0,3 - 134,06 - 34,71 - 7,15 - 45,33 - 143,20 0,6 - 190,41 - 55,52 - 16,11 - 66,13 - 199,54 0,9 - 244,63 - 74,21 - 22,96 - 84,82 - 253,76 5. Анализ результатов Представленные решения, как функции переменных x и y, обладают симметрией. При переходе от одного решения к другому и при одновременной замене y на - y, продольная скорость и давление не изменятся, тогда как поперечная скорость изменяет знак на противоположный u(1)(x,-y) = u(2)(x,y),p(1)(x,-y) = p(2)(x,y),v(1)(x,-y) = -v(2)(x,y). Например, для x = -0,4 и y = ±0,6 по табл. ?? находим значение продольной скорости для первого решения u(1)(-0,4;0,6) = 0,97. По табл. ?? для второго решения с противоположным знаком координаты y определяем u(2)(-0,4;-0,6) = 0,97. Имеем совпадающие значения u(1)(-0,4;0,6) = u(2)(-0,4;-0,6). Для поперечной скорости для первого решения по табл. ?? находим v(1)(-0,4;0,6) = -4,14, тогда как для второго решения по табл. ?? получаем v(2)(-0,4;-0,6) = 4,14. Указанные значения отличаются знаком. Для давления в случае первого решения из табл. ?? находим p(1)(-0,4;0,6) - p0 = 1,12, тогда как для второго решения из табл. ?? получаем такое же значение p(2)(-0,4;-0,6) - p0 = 1,12. Таким образом, оба решения в некотором смысле симметричны по отношению друг к другу и для выявления качественных особенностей достаточно проанализировать одно из них, например первое, которое представлено табл. ??, ??, ??. Рассмотрим значения продольной скорости u(1)(x;y) по табл. ??. Основным направлением движения является движение слева направо. При движении в этом направлении продольная скорость положительна, а при движении в обратном - отрицательна. Из табл. ?? следует, что положительное направление движения имеет основная масса жидкости, сосредоточенной в области - 0,3 < y < 0,6. И лишь в области вблизи границ y ± 1 продольная скорость отрицательна и движение происходит в обратном направлении. Однако при x > 0 направление движения изменяется. В области, прилегающей к центральной линии y = 0, движение происходит в обратном направлении, тогда как в области, прилегающей к границам y ± 1, движение происходит в основном направлении при положительной продольной скорости. Получается, что при движении вдоль всего канала - 1 < x < 1 существуют значительные области движения в обоих направлениях. Обнаруживаются довольно интенсивные возвратные движения, и основное направление движения теряется. Рассмотрим значения поперечной скорости v(1)(x;y) по табл. ??. При этом не забудем, что масштабы поперечной и продольной скоростей отличаются в k раз. В нашем случае k = 0,5 и для того, чтобы сравнение поперечной и продольной скоростей было адекватным, нужно все значения табл. ?? умножать на 0,5. С учётом этого заключаем, что по модулю поперечная скорость меньше продольной, но она ненулевая. Движение во всей области значительно отличается от одномерного. Расчёты показывают, что одномерное движение имело бы место только в случае θ = 0. Для этого случая система (??), (??) имела бы решения a10 = a20 = a11 = 0, и основные гидродинамические характеристики совпадали бы с (??). Однако для нашего случая θ = 0,1, и даже это сравнительно небольшое отличие параметра θ от нуля исключает наличие одномерного движения. Так что в рамках рассматриваемой постановки движение двумерно. Кроме того, из табл. ?? следует, что существуют значительные области, как с положительными, так и с отрицательными значениями продольной скорости v(x;y). Значит, имеются зоны с интенсивными вихревыми движениями. Значения, приведённые в табл. ??, подтверждают отмеченные выше закономерности. Из табл. ?? следует, что если говорить об области течения в целом, то величина p(x;y) - p0 изменяется как по знаку, так и значительно изменяется по абсолютной величине. Этот факт соответствует наличию и интенсивных возвратных течений, и интенсивных вихревых движений.

Об авторах

Александр Владимирович Коптев

Государственный университет морского и речного флота имени адмирала С. О. Макарова

Автор, ответственный за переписку.
Email: Alex.Koptev@mail.ru

доцент, кандидат физико-математических наук, доцент кафедры математики Института водного транспорта Государственного университета морского и речного флота имени адмирала С. О. Макарова

ул. Двинская, д. 5/7, г. Санкт-Петербург, Россия, 198035

Список литературы