Место издания:Издательский отдел физического факультета МГУ Отдел оперативной печати физического факультета МГУ Москва
Первая страница:78
Последняя страница:79
Аннотация:При изучении модели вязкости Педерсен, которая широко используется в современных гидродинамических симуляторах для моделирования вязкости многокомпонентных многофазных углеводородных растворов, были выявлены несколько сложностей, одна из которых это решение уравнения состояния BWR (Бенедикта — Вебба — Рубина), в виде предложенном McCarty. В этом уравнении содержится 32 члена, часть из которых экспоненциальные. Уравнение задает плотность метана с высокой точностью ~0.1% во всей области существования вещества. В некотором диапазоне температур само уравнение осциллирует, но по-прежнему существуют монотонные газовая и жидкая ветви, имеющие физический смысл. Необходимо было и в этом диапазоне температур найти правильный корень уравнения: при заданных давлении и температуре найти плотность нужной (жидкой или газовой) фазы. Задача решалась методом Ньютона, в качестве начальных приближений использовались более простые уравнения состояния, задающие плотность метана с меньшей точностью, а также уравнения, описывающие поведение метана на кривой насыщения (кривой кипения/конденсации). Показано, что более простые уравнения состояния не могут давать правильное начальное приближение во всем рассматриваемом диапазоне давлений и температур. Также показано, что уравнения, описывающие поведение метана на кипения/конденсации позволяют правильно решить задачу и сделать проверку полученного результата: принадлежит ли найденный корень нужной ветви уравнения.