Применительно к тем уравнениям, которые можно представить в форме (2.1), классификация ДУЧП может быть проведена путем проверки с использованием критериев (2.2). Когда же это невозможно, например в случае системы ДУЧП, для надлежащей классификации приходится обычно проверять поведение характеристик (см. п. 2.1.3).

Грубо говоря, различные категории ДУЧП могут ассоциироваться с различными типами гидроаэродинамических задач. В общем случае задачи, содержащие зависимость от времени сводятся к решению либо параболических, либо гиперболических ДУЧП. Параболические ДУЧП определяют собой течения, обладающие тем или иным механизмом диссипации, например такие, где существенны вязкие напряжения или теплопроводность. В этом случае решение будет гладким, а градиенты с течением времени будут уменьшаться, если только граничные условия не зависят от времени. Если же диссипативные механизмы отсутствуют, то решение будет сохранять постоянную амплитуду в случае линейности ДУЧП и будет даже возрастать в случае его нелинейности. Такое решение характерно для течений, определяемых гиперболическими ДУЧП. Эллиптические ДУЧП обычно определяют собой задачи, относящиеся к установившимся или к равновесным состояниям. Однако при исследовании некоторых установившихся течений приходится иметь дело с параболическими ДУЧП (установившееся течение в пограничном слое) и с гиперболическими ДУЧП (установившееся невязкое сверхзвуковое течение).

2.1.1. Природа корректно поставленной задачи

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

1) решение задачи существует,

2) решение является единственным,

3) решение непрерывно зависит от вспомогательных данных.

Вопрос о существовании решения обычно не создает каких-либо трудностей. Исключение составляет введение точных решений уравнения Лапласа (см. § 11.3), когда в некоторых изолированных точках решение может не существовать. Так,

например, оно не существует в центре источника, при г = и при использовании (11.53). Зачастую на практике удается избежать столкновения с этой проблемой путем расположения источника за пределами вычислительной области, например внутри тела, показанного на рис. 11.7.

Обычная причина неединственности связана с несоответствием вспомогательных условий типу определяющего ДУЧП. Для уравнения потенциала, определяющего невязкие безвихревые течения, а также для уравнений пограничного слоя надлежащая форма начальных и граничных условий хорошо установлена. Для уравнений Навье - Стокса хорошо известны надлежащие условия на твердой поверхности, однако существует некоторая неопределенность относительно надлежащего выбора граничных условий вдали от тела. В общем случае недоопределение граничных условий приводит к неединственности, тогда как их переопределение - к получению нефизических решений вблизи соответствующей границы.

Можно указать на некоторые проблемы обтекания, для которых множественность решений может ожидаться в силу физических причин. К таким проблемам уже нельзя применять вышеприведенные критерии корректности математической постановки. Именно такая ситуация нередко возникает при рассмотрении течений, испытывающих переход от ламинарного состояния движения к турбулентному. Однако при надлежащем представлении о гидроаэродинамике обычно удается идентифицировать такие виды течений, расчет которых может быть усложнен из-за необходимости заботы о корректной постановке математической задачи.

Третий из перечисленных выше критериев требует, чтобы малое изменение начальных или граничных условий влекло за собой лишь малое изменение решения. Довольно часто вспомогательные условия вводятся в типовой вычислительный алгоритм в приближенной форме. Следовательно, если третье условие не выполняется, то ошибки, введенные вместе со вспомогательными данными, будут распространяться во внутреннюю область, вызывая тем самым быстрый рост решения, особенно для гиперболических ДУЧП.

Вышеуказанные критерии обычно приписывают Адамару [Garabedian, 1964]. В дополнение к этому мы могли бы провести простую параллель и потребовать, чтобы для корректно поставленного расчета выполнялись следующие условия:

1) вычислительное решение существует,

2) вычислительное решение единственно,

3) вычислительное решение непрерывно зависит от приближенных вспомогательных данных.

Вводимые данные, а

Вычислительное решение, U

Рис. 2.1. Построение вычислительного решения

Таким образом, для проведения корректно поставленных вычислений необходимо, чтобы не только задача для исходного ДУЧП со вспомогательными условиями была корректно поставлена, но чтобы и алгоритм был также корректно поставлен (устойчив). При этом подразумевается, что приближенное решение, полученное в результате корректно поставленного расчета, будет в определенном смысле близко к точному решению корректно поставленной задачи. Этот вопрос будет рассматриваться в § 4.1.

2.7.2. Граничные и начальные условия

В результате проведенного в п. 2.1.1 обсуждения корректно поставленных задач и корректно поставленных вычислений становится ясно, что вспомогательные данные являются в определенном смысле отправной точкой для получения внутреннего решения, в особенности по отношению к задачам о распространении. Если мы не делаем различия между временем и пространством как независимыми переменными, то вспомогательные данные, вводимые на контуре dR (рис. 2.2), через посредство вычислительного алгоритма (основанного на ДУЧП) экстраполируются , обеспечивая тем самым построение решения во внутренней области R.

Вспомогательные условия определяются как относящиеся к одному из трех классов:

1) условие Дирихле, например u = f нг dR,

2) условие Неймана (для производной), например ди/дп = = f или du/ds = g на dR,

Процесс построения вычислительного решения может быть представлен схематически так, как это сделано на рис. 2.1. Здесь термин вводимые данные соответствует приближенному представлению начальных и граничных условий. Если граничные условия формулируются для производных функции и, то некоторая ошибка будет вноситься за счет приближенного лредставления этих условий. Вычислительный алгоритм, как правило, строится на основе определяющего ДУЧП (см. § 3.1) и должен быть устойчивым (см. § 4.3), чтобы выполнялись все три вышеприведенных условия.