В конечных разностях схемную вязкость оценивают следующим образом: берут дифференциальный оператор, представляют его разностным оператором. Затем от разностного оператора отнимают слагаемые из дифференциального уравнения, которые разложены в ряд Тейлора. Полученный остаток и называют схемной, а кто-то и аппроксимационной вязкостью (тут тоже вопрос в терминологии). Как подступиться к подобным оценкам в МКЭ пока не представляю. Может кто-то делал подобные вещи?
|