ИСТИНА

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