Статьи журнала - Вычислительная механика сплошных сред
Все статьи: 692
Нелинейные режимы конвекции бинарной смеси в двухслойной пористой среде различной конфигурации
Статья научная
Представлены результаты численного исследования возникновения конвекции и ее нелинейных режимов в бинарной смеси жидких углеводородов, заключенной в двухслойной пористой среде. Состав смеси, тепловые условия и геометрия позволяют интерпретировать задачу как простую модель литологически ограниченной залежи месторождения углеводородов. Расчетная область представляет собой вытянутый в горизонтальном направлении прямоугольник, разделенный на два слоя. В одном из рассмотренных случаев эти слои равной высоты, а в другом - разной, так как граница раздела слоев имеет форму дуги окружности, обращенной выпуклостью вниз, что имитирует синклинальную геологическую складку. Слои обладают равной пористостью, но разной проницаемостью. В области существует геотермальный градиент со средним значением температуры, характерным для глубины 2000 м, что соответствует средней глубине залегания нефти. Задача решается в рамках модели Дарси-Буссинеска с учетом эффекта термодиффузии. Прослежены процессы образования и установления нелинейных режимов конвекции, распределение концентрации компонентов смеси для различных значений проницаемости слоев и их зависимость от того, какой из этих слоев лучше проницаем. Найдено, что в слоях с близкими значениями устанавливается стационарный режим течения при любой из конфигураций границы раздела, а также формируется «крупномасштабная» конвекция. При значительной разности проницаемостей возникающая конвекция «локальна», могут наблюдаться либо квазипериодические колебания сложной формы, либо нерегулярные колебания. При этом проникновение течения в слой с меньшей проницаемостью незначительно.
Бесплатно
Нелинейные режимы конвекции трехкомпонентной смеси в двухслойной пористой среде
Статья научная
Численно исследованы возникновение и нелинейные режимы конвекции смеси, состоящей из метана (35%), этана (35%) и бутана (30%), в вытянутой в горизонтальном направлении прямоугольной области пористой среды под действием геотермального градиента. Область имеет твердые, непроницаемые для вещества границы и образована двумя горизонтальными слоями, высоты которых соотносятся как 1:3. Слои характеризуются равной пористостью, но отличаются значениями проницаемости. Величины пористости и проницаемости выбирались близкими к значениям реальных сред, таких как пески, песчаники или известняки. Компоненты анализируемой смеси принадлежат к основным группам химических соединений, присутствующих в почвах нефтегазовых месторождений. Таким образом, описанная конфигурация представляет собой модель углеводородной залежи. Рассмотрены случаи, когда значение проницаемости верхнего слоя выше, чем нижнего, и, наоборот, нижний слой более проницаем, чем верхний. Остальные параметры пористой среды считаются одинаковыми во всей расчетной области. Задача решается в рамках модели Дарси-Буссинеска с учетом эффекта термодиффузии. Прослежена временная эволюция локальных характеристик и структуры формирующегося течения и распределения компонентов смеси. В более проницаемом слое меньшей высоты показан «локальный» характер возникновения конвекции. Именно при таком сочетании высоты и проницаемости в этом слое зарождается тчение и, в процессе развития конвекции, оно начинает проникать в менее проницаемый слой, однако центры формирующихся вихрей заметно сдвинуты в сторону более проницаемого слоя. Аналогичное смещение вихрей наблюдается и в толстом слое, если он обладает большей проницаемостью. Однако в этом случае конвекция носит «крупномасштабный» характер.
Бесплатно
Нелинейный изгиб и устойчивость эллипсоидальной оболочки, взаимодействующей с жидкостью
Статья научная
Рассматриваются осесимметричный изгиб и устойчивость сплюснутой эллипсоидальной оболочки под действием давления жидкости с выпуклой стороны оболочки, а также при контакте оболочки с несжимаемой жидкостью. Используются геометрически и физически нелинейные соотношения теории оболочек. Давление жидкости изменяется по адиабатическому закону. Полученные результаты свидетельствуют о зависимости нелинейного деформирования от характеристик нагружающей среды. Также показано, что особенности напряженно-деформированного состояния оболочки определяются переменностью кривизны срединной поверхности.
Бесплатно
Нестационарная слоистая тепловая и концентрационная конвекция Марангони вязкой несжимаемой жидкости
Статья научная
Построены и проанализированы точные стационарные и нестационарные решения задачи слоистой конвекции Марангони, которая является переопределенной краевой задачей, а также ее численное решение, принадлежащее к классу решений Бириха. Переопределенность разрешающей системы уравнений возникает вследствие равенства нулю скорости, параллельной оси аппликат. Рассмотрены случаи тепловой и концентрационной конвекции вязкой несжимаемой жидкости. Для разрешимости краевой задачи предложено использовать класс точных решений, в котором скорости одномерны по координатам, а поля давления и температуры являются трехмерными. Характерная особенность этого класса - тождественное обращение в нуль конвективной производной в уравнении сохранения импульса. При этом конвективная производная присутствует в калорическом уравнении состояния. Показано, что рассматриваемая краевая задача, в отличие от классического решения Бириха и его многочисленных обобщений, не может быть сведена к одномерной при задании градиента температуры на обеих границах слоя жидкости. Найденные в данной работе стационарные и нестационарные решения имеют в профиле скоростей застойную точку, что говорит о наличии противотечений при движении жидкости. Методами локализации корней полиномов стационарных решений продемонстрировано, что существует такое значение толщины слоя, при котором касательное напряжение может обратиться в нуль на нижней границе слоя жидкости только при тепловой конвекции Марангони. Полученные методом граничных элементов нестационарные решения, которые можно трактовать как точные, с течением времени выходят на стационарные решения. Применение метода граничных элементов существенно расширяет класс точных нестационарных решений, поскольку позволяет изучать и те из них, которые не обладают свойством инвариантности.
Бесплатно
Нестационарный изгиб консольно-закрепленной балки Бернулли-Эйлера с учетом диффузии
Статья научная
Рассматривается нестационарная задача изгиба консольно-закрепленной упруго-диффузионной однородной изотропной балки Бернулли-Эйлера. Математическая постановка представляет собой замкнутую систему уравнений поперечных нестационарных колебаний балки с учетом диффузии. Разрешающая система уравнений изгиба балки получена из общей модели упругой диффузии в сплошной среде с помощью вариационного принципа Даламбера. При этом предполагается, что прогибы балки являются малыми и выполняется гипотеза плоских сечений. Повороты сечений удовлетворяют гипотезе Бернулли-Эйлера. Решение ищется методом эквивалентных граничных условий, который позволяет перейти от исходной формулировки с произвольными граничными условиями к задаче того же вида, с той же геометрией области, но с заведомо достижимым решением. В данной работе в качестве вспомогательной выступает такая задача, решение которой находится путем интегрального преобразования Лапласа по времени и разложением в тригонометрические ряды Фурье. Далее строятся соотношения, связывающие правые части граничных условий исходной и вспомогательной задач. Они являются интегральными уравнениями Вольтерра 1-го рода и образуют систему, которая разрешается численно с использованием квадратурных формул средних прямоугольников. В итоге решение исходной задачи представляется в виде сверток функций Грина вспомогательной задачи с функциями, определяемыми из решения системы интегральных уравнений Вольтерра. На примере двухкомпонентного материала выполнено численное исследование взаимодействия нестационарных полей - механического и диффузионного, в изотропной балке. Результаты вычислений представлены в виде графиков зависимости искомых полей перемещений балки и приращений концентраций компонентов среды от времени и координат. Исходя из их анализа сделан вывод о влиянии связанности механодиффузионных полей на напряженно-деформированное состояние и массоперенос в балке.
Бесплатно
Нестационарный износ двухслойного покрытия с учётом тепловыделения от трения
Статья научная
Применение двухслойных покрытий является одним из распространённых путей повышения работоспособности высоконагруженных узлов трения. Проектирование двухслойных покрытий и прогнозирование их эффективности в условиях сухого трения скольжения требуют разработки математических моделей скользящего контакта, учитывающих износ и разогрев вследствие трения. Для этого в статье рассматривается нестационарная квазистатическая контактная задача несвязанной термоупругости в случае скольжения с постоянной скоростью жёсткого штампа, имеющего вид полуплоскости, по поверхности двухслойного упругого покрытия, сцепленного по нижней грани с плоской недеформируемой подложкой. Во время процесса жёсткий штамп не только скользит по поверхности термоупругого покрытия, но и внедряется в него по нормали. Трение штампа о покрытие описывается законом Амонтона-Кулона. Возникающий тепловой поток от трения направлен вглубь покрытия. Фрикционный скользящий контакт приводит к изнашиванию двухслойного покрытия...
Бесплатно
Неупругое взаимодействие и расщепление солитонов деформации, распространяющихся в стержне
Статья научная
Дается краткий аналитический обзор публикаций, посвященных теоретическому изучению возможности формирования и особенностей распространения солитонов продольной деформации в нелинейно-упругих стержнях, а также экспериментальному обнаружению таких волн. Основное внимание уделяется изложению результатов оригинальных исследований взаимодействия солитонов. Путем численного моделирования показано, что качественно различные сценарии такого взаимодействия зависят от относительной скорости столкновения солитонов. Так, при малой скорости столкновение происходит по сценарию, свойственному классическим солитонам, описываемым уравнением Кортевега-де Вриза, то есть вторичные солитоны сохраняют те же скорость, амплитуду и ширину, что и первичные солитоны. При большей относительной скорости столкновение имеет неупругий характер: часть энергии солитонов теряется при воздействии их друг на друга и реализуется в пакеты квазигармонических волн, движущихся со скоростью линейной волны. Дальнейшее увеличение скорости столкновения приводит к эффекту расщепления солитонов, под которым авторы данной работы понимают образование большего, чем вступавшее в связь, количества вторичных солитонов. Кроме согласованного поведения солитонов между собой, исследовано их взаимодействие с границей стержня.
Бесплатно
Низкодиссипативная фильтрация в разностных методах вычислительной газовой динамики
Статья научная
Сопоставлены различные методики фильтрации решения, такие как: явная низкодиcсипативная фильтрация (DRP), фильтрация с дополнительным выборочным фильтром низкого порядка для сглаживания скачков (SC), вейвлет-фильтрация, компактная фильтрация с регулируемой пропускной способностью. Показано, что хотя 11-точечная DRP-фильтрация и является одним из самых простых и быстрых методов, тем не менее, имеет более выраженные диссипативные свойства по сравнению с вейвлет-фильтрацией на основе 6-коэффициентной функции Добешии и с 5-точечной трехдиагональной компактной схемой фильтрации. Для представленных задач наименее диссипативными получились вейвлет-фильтры. А среди фильтров, дающих наименьший уровень паразитных пульсаций решения в области скачка, помимо фильтров низкого порядка и схем высокого порядка с использованием SC-фильтров, выделен компактный фильтр с регулируемой пропускной способностью.
Бесплатно
Новый итерационный метод расчета многолистовой рессоры
Статья научная
Предложен новый итерационный метод решения контактной задачи, возникающей в классической модели многолистовой рессоры. Рессора представлена в виде пачки консольных балок (листов) постоянного поперечного сечения. Длины и толщины листов различны, ширина одинакова. В отсутствие нагрузки листы являются плоскими и плотно прилегают друг к другу; трение не учитывается. К свободному краю самого длинного листа приложена сосредоточенная сила, под действием которой листы испытывают слабый статический изгиб; контакт листов происходит с возможным отставанием. Задача состоит в отыскании сил взаимодействия листов. Новый метод основан на полученном авторами ранее аналитическом решении обобщения поставленной выше задачи (на листы действуют силы с произвольными заданными плотностями), но для двухлистовой рессоры. Для многолистовой рессоры это решение используется для каждой пары соседних листов при итерационном уточнении сил взаимодействия. Получен ряд ранее не известных картин контакта, которые не могут быть найдены традиционными методами.
Бесплатно
О вихревых течениях двухфазной жидкости в пористой среде
Статья научная
Представлены результаты численного моделирования двухфазного изотермического течения несмешивающихся несжимаемых жидкостей в пористой среде с учетом капиллярных и гравитационных сил. Задача решена методом контрольных объемов в переменных «скорость-насыщенность». Для решения уравнения гиперболического типа использована схема WENO третьего порядка точности в комбинации с методом Рунге-Кутты. Показано, что под действием капиллярных и гравитационных сил при определенных условиях могут возникать вихревые потоки в каждой из фаз. В частности, вихревые течения наблюдаются как при внедрении тяжелой жидкости в пористое тело, насыщенное легкой жидкостью, так и при сегрегации жидкостей в пласте, если в нем имеются включения другой проницаемости. Также показано, что под влиянием вихревых потоков жидкости меняется конфигурация фронта вытеснения. Ранее явление образования вихревых течений в пористой среде изучалось посредством введения в рассмотрение новой искомой переменной - функции завихренности. В такой постановке задачи изучено появление вихревых течений в пористых средах на границах разделов двух жидкостей разной плотности, на границах скачкообразного изменения проницаемости и в некоторых других случаях. В настоящей работе показано, что наличие завихренности при двухфазном течении в пористой среде может быть выявлено путем непосредственного вычисления поля скоростей каждой из фаз.
Бесплатно
Статья научная
Узел опорной части мостового сооружения состоит из верхней стальной плиты со сферическим сегментом, обеспечивающим подвижность конструкции, нижней стальной плиты, отвечающей за передачу вертикальной опорной реакции и равномерное распределение давления от сферического сегмента, и антифрикционных полимерных слоев скольжения (сферического и плоского). При изготовлении узла предусматривается разное местоположение слоев скольжения относительно стальных плит. В работе анализируется влияние места нанесения сферического слоя скольжения в опорном узле на фрикционное контактное взаимодействие элементов опорной части и на их деформационное поведение. Плоский слой скольжения в модели не учитывается. Авторами не найдены в литературе сведения о подобной оценке связи геометрии и конфигурации антифрикционных прослоек с поведением конструкции при деформировании. Рассмотрено два варианта размещения слоя скольжения: он нанесен или на сферический стальной сегмент, или заполняет сферический вырез нижней стальной плиты (классический вариант). В качестве материала слоя скольжения используется радиационно-модифицированный фторопласт Ф-4 без наполнения. Опорный узел с прослойкой, расположенной в вырезе нижней стальной плиты, соответствует модели опорной части Л-100 производства ООО «АльфаТех» (г. Пермь). Модель Л-100 рассчитана на нормативную вертикальную нагрузку 1000 кН и обладает конструктивными размерами: максимальной длиной 155 мм, максимальной высотой 54 мм, толщиной 4 мм. Опорный узел с прослойкой на сферическом сегменте имеет аналогичные модели Л-100 геометрические показатели. Стандартный угол наклона торца антифрикционной прослойки составляет 30°. Установлено, что при стандартном угле в опорной части с прослойкой на сферическом сегменте наблюдается расхождение (отлипание) поверхностей сопряжения более чем на 2% ее общей площади. Поэтому в рамках данного исследования оценивается роль наклона торца антифрикционной прослойки на деформирование опорных частей. В численных экспериментах установлен ряд преимуществ классической конструкции сферической опорной части, показавшей при моделировании большую площадь полного прилипания поверхностей сопряжения элементов, равномерное распределение параметров контакта по поверхностям, меньшее деформирование профиля торца слоя скольжения. Также установлены значения угла наклона торца слоя скольжения, при которых получается благоприятное для функционирования мостовой конструкции в зоне опорного узла распределение свойств зон контакта и характеристик деформированного состояния опорных частей при двух вариантах местоположения антифрикционной прослойки.
Бесплатно
О влиянии геотермальной системы на деформации земной поверхности во время вулканического извержения
Статья научная
Измерение деформаций земной поверхности, происходящие в процессе вулканического извержения, являются одним из основных методов мониторинга активных вулканов. Источником деформаций служат процессы в системе очаг-канал, а также в геотермальных системах, прогреваемых поднимающейся магмой. Роль очага и в меньшей степени канала при деформировании окружающих пород достаточно хорошо изучена теоретически, однако влияние геотермальной системы на измеряемые деформации в процессе вулканического извержения ранее не рассматривалось. В статье сравниваются деформации от двух инициирующих деформации источников: щелевого канала с заданным избыточным давлением и геотермальной системы, прогрев которой осуществлялся за счет течения магмы. Показано, что вертикальные деформации вследствие активности геотермальной системы могут в разы превосходить деформации, связанные с течением магмы. Пространственное распределение деформаций также существенно различается. При деформациях, вызванных геотермальной системой, максимум вертикального перемещения находится над щелевым каналом, в случае же изменения давления в канале над ним наблюдается локальное опускание поверхности, а ее максимальное поднятие располагается на расстоянии, примерно вдвое превышающем глубину залегания верхней части щелевого канала. Влияние геотермальной системы необходимо учитывать при интерпретации данных мониторинга активных вулканов.
Бесплатно
О влиянии силы тяжести на упругие волновые процессы в идеальном газе
Статья научная
Рассматривается влияние силы тяжести на упругие волновые процессы в идеальном газе. Для описания таких процессов используется простейшая линейная баротропная энергетически замкнутая модель. При этом в волновом уравнении присутствует линейный поправочный член, учитывающий как силу тяжести, так и обусловленное ею изменение плотности газа по высоте. На примере типовых волновых процессов (плоской гармонической, сферической гармонической и нестационарной сферической волн) показывается, какие изменения вносит эта поправка в аналитические представления таких процессов и как отражается на их характере.
Бесплатно
Статья научная
Интерес к конвекции в условиях фазовых переходов, свойственных испаряющимся жидкостям и парогазовым смесям, вызван их широким применением в различных технологиях и промышленных процессах. Альтернативой натурной отработке технологических методик и экспериментальным исследованиям конвективного тепломассообмена являются методы математического моделирования. В настоящей работе в рамках приближения Обербека-Буссинеска рассматривается задача испарительной конвекции в двухфазной системе жидкость-газ. Для описания установившихся адвективных течений при наличии диффузионного испарения используется частично-инвариантное точное решение определяющих уравнений, корректно учитывающее влияние внешней тепловой нагрузки и термодиффузионных эффектов в газопаровом слое. С помощью точного решения оценивается влияние толщины жидкого слоя на кинематические и температурные характеристики возникающих режимов, параметры фазового перехода и паросодержания в несущем газе. Установлено, что с ростом толщины слоя происходит последовательная смена течения чисто термокапиллярного типа течениями смешанного и пуазейлевского типов и имеют место качественные изменения в процессах массопереноса через межфазную границу. Методом нормальных мод анализируется линейная устойчивость точного решения относительно плоских и пространственных возмущений. Найдены пороговые характеристики устойчивости, демонстрируется эволюция топологии нейтральных кривых и форм неустойчивости при изменении геометрии системы. Показано, что увеличение толщины жидкого слоя дестабилизирует течение, при этом в системе всегда реализуется колебательная неустойчивость. Представлены зависимости фазовых скоростей возмущений для систем разной геометрии. Выявлено, что формы неустойчивости в слое испаряющейся жидкости, увлекаемой спутным газовым потоком, предсказываемые на основе точного решения, совпадают с наблюдаемыми в теплофизических экспериментах.
Бесплатно
О выделении разрывов в расчетах динамики несжимаемой упругой среды
Статья научная
Задача алгоритмического выделения близко отстоящих друг от друга разрывов деформаций решается на примере одномерных цилиндрических поверхностей разрывов (ударных волн). Ударные волны создаются в цилиндрическом слое несжимаемой нелинейной упругой среды посредством скручивающего ударного воздействия в присутствии предварительно полученных антиплоских деформаций. Деформации среды (как предварительные, так и приобретенные в результате внешнего нагружения) полагаются конечными, в качестве их меры используется тензор Альманси. Показано, что граничное ударное возмущение вызывает в среде две фронтальные поверхности разрыва деформаций: плоскополяризованную ударную волну нагрузки, увеличивающую предварительный антиплоский сдвиг, и нейтральную волну круговой поляризации, меняющую направление сдвига в соответствии с производимым воздействием. Найдены скорости распространения образующихся поверхностей разрыва в среде. Выявлено, что скорость распространения плоскополяризованной волны нагрузки зависит от предварительных деформаций среды и интенсивности граничного возмущения. Скорость распространения ударной волны круговой поляризации (нейтральной ударной волны) полностью определяется предварительными деформациями среды. Для целей алгоритмического выделения при вычислении положения поверхностей разрывов и интенсивностей разрывов на каждом временном шаге строятся специальные прифронтовые лучевые разложения решений, которые встраиваются в расчетные конечно-разностные схемы, записанные для узлов сетки, не входящих в прифронтовую область. Указан способ построения прифронтовых лучевых разложений за поверхностями разрывов деформаций, опирающийся на рекуррентность свойств геометрических и кинематических условий совместности разрывов. Построен численный алгоритм реализации поставленной задачи и создана программа для расчета полей перемещений и компонент тензора напряжений Коши. Осуществлен вычислительный эксперимент для резиноподобного материала с заданными свойствами.
Бесплатно
Статья научная
В рамках теории больших деформаций приведено решение одномерной краевой задачи определения НДС материала с упругими, вязкими и пластическими свойствами в окрестности одиночного сферического дефекта сплошности при нагрузке и последующей разгрузке. Установлена закономерность продвижения упругопластической границы, приведены расчетные уровни и распределение остаточных напряжений. Указаны условия возникновения повторного пластического течения в процессе разгрузки.
Бесплатно
О механизмах каскадного переноса энергии в конвективной турбулентности
Статья научная
В работе изучаются особенности каскадных процессов в развитой турбулентности, существующей на фоне градиента плотности (температуры), либо сонаправленного с вектором силы тяжести (турбулентность в устойчиво стратифицированной среде - далее УС), либо противонаправленного ему (конвективная турбулентность - КТ). Основное внимание уделено режиму Обухова-Болджиано (ОБ), подразумевающего баланс сил Архимеда и нелинейных сил в достаточно протяженной части инерционного интервала. Уверенного подтверждения того, что режим ОБ возможен, не получено до сих пор, хотя фрагменты спектров с наклонами, близкими к «-11/5» и «-7/5», были зарегистрированы в некоторых работах по численному моделированию конвективной турбулентности. Проводится критическое сравнение этих данных с результатами расчетов, выполненных авторами настоящей статьи с помощью каскадной модели конвективной турбулентности, позволившей рассмотреть широкий диапазон значений управляющих параметров. Каскадная модель является новой и получена путем обобщения класса спиральных каскадных моделей на случай турбулентной конвекции. Показано, что в режимах развитой турбулентности, признаком которой является интервал с постоянным спектральным потоком кинетической энергии, силы Архимеда не могут конкурировать с нелинейными взаимодействиями и не оказывают существенного влияния на динамику инерционного интервала. В случае КТ именно они обеспечивают каскадный процесс энергией, но только на максимальных масштабах турбулентности. При УС силы плавучести снижают энергию турбулентных пульсаций. Но ни в том, ни в другом случае режим ОБ не возникает, а на масштабах, попадающих в инерционный интервал, устанавливается колмогоровская турбулентность с законом «-5/3», в которой температура ведет себя как пассивная примесь. Наблюдаемые отклонения от спектра «-5/3», ошибочно интерпретируемые как режим ОБ, появляются при недостаточном разделении макромасштаба турбулентности и диссипативного масштаба.
Бесплатно
О моделировании динамических процессов в сферических и цилиндрических оболочках
Статья научная
Различия между результатами математического моделирования динамических процессов в сплошных средах и результатами физических экспериментов определяется тремя видами погрешности: погрешностью физической модели (∆Ф), погрешностью математической модели (∆М) и погрешностью методов измерений (∆И). В статье на примере cреды в сферически симметричной оболочке сравниваются ∆Φ и ∆М. Показано, что в некотором диапазоне скоростей диссипация энергии, определяемая погрешностью конечно-разностной аппроксимации, может превосходить диссипацию энергии, порождаемую моделью идеальной пластичности.
Бесплатно
Статья научная
Дан краткий обзор публикаций журнала «Вычислительная механика сплошных сред» за 2008-2011 годы с целью анализа направлений развития вычислительных технологий в механике сплошных сред. Выделены восемь направлений и приведена их краткая характеристика на основе опубликованных работ.
Бесплатно
О некоторых моделях нанотрибологии
Статья научная
Предложена новая модель трения на наноуровне, представляющая собой обобщение классических моделей Томлинсона и Френкеля-Конторовой на случай учета дальнодействия в рамках теории центральных парных взаимодействий. Построены численные решения задач об определении сил трения, распределений сил и перемещений в деформируемой цепочке для разных моделей трения: Томлинсона, Френкеля-Конторовой и предложенной в работе обобщенной одномерной модели. Проведено сравнение решений, выявлены некоторые важные механические эффекты и указаны области применимости классических моделей.
Бесплатно
