В области численного моделирования на мультикомпьютерах стоит проблема разработки эффективных параллельных программ (ПП). ПП могут работать эффективно только на определенных классах вычислительных систем (ВС). Для использования других ВС, вообще говоря, потребуются разные ПП. Возникает вопрос о переносимости ПП с учетом различных особенностей ВС. Системы автоматизации конструирования ПП позволяют обеспечивать такую переносимость, в том числе система LuNA. Она реализует подход фрагментированного программирования и исполняет фрагментированные программы (ФП). Не любая ПП является ФП, поэтому требуется преобразовать ее в ФП. Развитие отрасли знаний, позволяющей выполнять такое преобразование, может происходить при помощи исследования ручной разработки ФП для частных случаев среди задач численного моделирования. Данная работа посвящена такому исследованию для задачи моделирования процесса фильтрации двухфазной жидкости на основе законов сохранения в интегральной форме.