Рекурсия при задании граничного условия

[dolivanov 7]

Всем привет.

При моделировании течения крови по артериям нужно задать давление на выходе из сосуда, которое пересчитывается через ряд параметров на каждом временном шаге. Давление зависит, в первую очередь, от расхода на этом выходе.

Модель, которая позволяет это реализовать, называется Windkessel.

Для внедрения этой модели надо на каждом временном шаге рассчитывать значение рекурсивной функции по формуле. Формулу прилагаю. В ней слева функция на текущем шаге, справа есть параметры на текущем шаге, параметры на предыдущем шаге и сама функция на текущем шаге.

 

Как такую рекурсивную функцию реализовать? Если в виде expression, то тогда непонятно как инициализировать эту переменную.

Спасибо.

 

[sportsfun 1]

В cfx это в явном виде не реализовано, насколько я знаю. Можно попробовать через subroutine, см. файл

Time operations in CFX.doc