Разработка алгоритмов случайного блуждания для решения нестационарных задач математической физики тема автореферата и диссертации по математике, 01.01.07 ВАК РФ

ВВЕДЕНИЕ.

ГЛАВА I. АЛГОРИТМ БЛУЖДАНИЯ ПО ГРАНИЦЕ ДЛЯ РЕШЕНИЯ НЕСТАЦИОНАРНЫХ КРАЕВЫХ ЗАДАЧ.

§1.1. Постановка задачи и основные определения

§1.2. Алгоритмы случайного блуждания по границе для решения интегральных уравнений теории тепловых потенциалов.

§1.3. Учет неоднородностей в правой части и в началь ном условии уравнения теплопроводности

§1.4. Доказательство конечности дисперсии и исслет-. . дование трудоемкости. Метод установления

§1.5. Вычисление производных на границе области от решения первой краевой задачи с однородным граничным условием.

ГЛАВА 2. МЕТОЛУ БЛУЖДАНИЯ ВНУТРИ ОБЛАСТИ И ИССЛЕДОВАНИЕ

ИХ ТРУДОЕМКОСТИ.

§2.1. Соотношение о. среднем, для. уравнения. теплопрот-, водности.*.

§2.2. Алгоритм блуждания по сферам для решения урав . нения теплопроводности - &и Ся^).

§2.3. Алгоритм блуждания по сферам для решения, урав , нения и^- Дисх^-ьесхуь)*^ £с*^)

§2.4. Оценка среднего числа шагов до попадания в е--окрестность границы.в алгоритмах,блуждания. . внутри области

ГЛАВА 3. АЛГОРИТМЫ СЛУЧАЙНОГО БЛУЖДАНИЯ ДЛЯ РЕШЕНИЯ ПЕРВОЙ КРАЕВОЙ ЗАДАЧИ ДЛЯ ПАРАБОЛИЧЕСКОГО УРАВНЕНИЯ С ПЕРЕМЕННЫМИ КОЭФВДИЕНТАМИ ПРИ МЛАДЦИХ ЧЛЕНАХ.

§3.1. Соотношение о среднем для параболического уравнения второго порядка.с переменными коэффициентами. при младших членах.

§3.2. Решение задачи Коши для уравнения (I.I) методом Монте-Карло.

§3.3. Решение первой краевой задачи.для,уравнения.(IЛ) . методом блуждания по сферам.

§3.4. Использование сопряженных задач для оценивания. многих функционалов методом Монте-Карло

ГЛАВА 4. РЕШЕНИЕ НЕКОТОРЫХ' МОДЕЛЬНЫХ И ПРИКЛАДНЫХ' ДЦШУШСШ

НЫХ ЗАДАЧ.

§4.1. Численные эксперименты по сравнению алгоритмов,. , блуждания по сферам и границе

§4.2. Численное решение задачи о диффузионном осаждении полидисперсных.аэрозолей.в трубах.методом блуждав ния по сферам.

В последнее время методы статистического моделирования, традиционно используемые в задачах теории переноса и статистической физики (см., например, [15, 16, 32] )получили интенсивное развитие при решении многомерных краевых задач математической физики. Большой интерес к данному подходу вызван рядом причин, связанных как с перспективой создания новых эффективных численных методов для решения практических задач, так и с наличием глубоких связей между дифференциальными уравнениями и случайными процессами, требующих всестороннего изучения. Такая связь известна уже давно: теория дифференциальных уравнений широко использовалась в теории вероятностей. Например, А.Н.Колмогоров [21] показал в 1931г., что переходная функция РС'Ц'эс^Г) - вероятность того, что броуновская частица, вышедшая из точки х , через время "Ь попадает в множество Р, удовлетворяет некоторому параболическому дифференциальному уравнению.

С другой стороны, аппарат теории вероятностей с успехом применялся при исследовании краевых задач; отметим, в частности, работы Р.Куранта, К.Фридрихса, Я.Леви [58], К.Дуба[59], Е.Б.Дын-кина [ю], А.Д.Вентцеля [з] , и др.

В методах Монте-Карло также использовались вероятностные представления решений краевых задач/в виде континуальных интегралов, при этом соответствующий стахостический процесс аппроксимировался дискретным случайным блужданием (см. например,[7,8]).

Этот подход прост в реализации, не требует большой памяти ЭВМ, но является сравнительно трудоемким из-за необходимости моделирования "длинных" случайных траекторий. Существуют и другие, отличные от алгоритма»основанного на моделировании дискретных случайных блужданий, приближенные методы вычисления континуальных интегралов - как детерминированные (см., например,[б, 7,57] ), так и методы Монте-Карло, использующие, например, другие аппроксимации винеровского процесса (см. например,[44, 54] ).

В случае уравнения Лапласа специфика винеровских интегралов позволяет для их оценки строить траектории менее подробно, например, достаточно моделировать выход винеровского процесса на границу области. Обоснование такого метода, известного под названием алгоритма блуждания по сферам, впервые описанном, по-видимому, Брауном в [2] , было дано в [61] . После появления этих работ начаты интенсивные исследования по построению и изучению алгоритмов метода Монте-Карло, основанных на использовании соотношения о среднем и формул Грина (см., например, [12,13,14,17,22,24,28, 35,37,42,43,48,50,51,52.] ).

Г. А.Михайловым и Б.С.Елеповым [12] впервые предложен подход, основанный на сведении исходной дифференциальной задачи к специальному интегральному уравнению с обобщенным ядром, что дало возможность использовать развитый аппарат методов Монте-Карло, для решения линейных и нелинейных интегральных уравнений 2-го рода [15,16] . Именно этот подход позволил разработать методы статистического моделирования для решения методом Монте-Карло, более общих уравнений и систем (см., например,[43,47,48,52] ). Первая монография [14] на эту тему опубликована в 1980 г.

Важный вопрос исследования трудоемкости методов статистического моделирования рассматривался в ряде работ. Как известно см., например, [16] ), трудоемкость оценок метода Монте-Карло определяется дисперсией этой оценки и количеством операций, необходимом для получения одного выборочного значения. Это количество операций пропорционально среднему числу блужданий до попадания марковской цепи в поглощающее, состояние.

В [бо! получена логарифмическая оценка среднего числа блуждания по сферам до выхода в ¿-окрестность границы. Отметим, что в работе [61] подробно исследовался вопрос трудоемкое-ти в зависимости от размерности пространства, однако логарифмическая оценка в ней не приводится.

Независимо от [60] в [13 ] Б.С.Елеповым и Г.А.Михайловым логарифмическая оценка и конечность дисперсии были получены на основе полуэвристической оценки дяя плотности среднего числа, сфер вблизи границы. Обобщение логарифмической.оценки на случай невыпуклых областей осуществлено.в работах Г.А.Михайлова [16], А.А.Кронберга [23] и автора [26].

Следует отметить, что во всех вышеперечисленных работах алгоритмы метода Монте-Карло строились в рамках предположения о сходимости ряда Неймана, в связи с чем не было попыток использовать представление решений стационарных задач в виде потенциалов простого и двойного слоя, плотность которых удовлетворяет интегральным уравнениям со спектральным радиусом, равным единице. На преодоление этих трудностей направлены две следующие группы работ, использующие два различных подхода. В заметке К.К.Сабельфельда [42] для решения граничных интегральных уравнений теории потенциала предложен алгоритм "блуждания по границеV основанный на аналитическом продолжении резольвенты этих уравнений. В работах К.К. Сабельфельда и Н.А.Симонова [46,47] дана оценка трудоемкости, проведено обобщение на случай системы Ламе. Второй подход связан с переходом к нестационарным задачам и применением метода установления (см., например,[зо]): дело в том, что при использовании тепловых потенциалов интегральные уравнения душ плотности имеют'нулевой спектральный радиус (уравнения Вольтерра по времени), и здесь можно применять развитую технику решения интегральных уравнений методом Монте-Карло. Данный подход был использован в работах [27,30,45] , где ососбое внимание было уделено разработке алгоритмов статистического моделирования для решения нестационарных задач.

Остановимся вкратце на особенностях алгоритма блуждания по границе. Во - первых, важным свойством этого метода является его "глобальность", т.е. возможность на одной траектории случайного блуждания получать оценки решения одновременно в произвольном множестве точек. Следует отметить, что существует и "глобальный" метод блуждания по сферам, предложенный в [48) , однако он применим лишь для простейших уравнений с однородными граничными условиями и симметричной функцией Грина.

Вторая особенность алгоритма "блуждания по границе" заключается в том, что он позволяет решать как внутренние, так и внешние задачи Дирихле и Неймана, а также смешанные краевые задачи.

В третьих, в алгоритме "блуждания по границе" существенно проще строить оценки производных, поскольку прямо из соотношения о среднем получить выражение для производных удается лишь в простейших случаях.

Наконец, в отличие от метода блуждания по сферам, здесь не приходиться прибегать к "сносу" граничных условий, что может оказаться существенным в задачах с большими градиентами вблизи границы.

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

Следует отметить, что методы статистического моделирования применялись в основном для решения стационарных задач математической физики; имеется небольшое число работ, посвященных развитию метода Монте-Карло для решения нестационарных задач, отметим здесь работы А.Хаджи-Шейха, Спэрроу, А.З.Веселовской, А.С.Сипина, Е.О.Яжетываева и К.К.Сабелъфельда [4, 9, 17, 55] . В [55] автор исходил из вероятностного представления решения и обобщил метод блуждания по сферам для уравнения теплопроводности. В работах |~4, 9] применяются интегральные преобразования типа Лапласа для перехода к стационарному случаю, а затем строятся специальные коррелированные оценки при обращении образов.

Обобщение алгоритма блуждания по сферам осуществлено А.С.Сипиным [17] и основано на соотношениях о среднем, полученных в [24] . В работе [28] проведено дальнейшее обобщение этого метода и его усовершенствование.

Итак,можно выделить три подхода, используемые в методах Монте-Карло для решения краевых задач:

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

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

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

Следует дать некоторое поясление к данной классификации. Во втором подходе интегральное уравнение также обычно записывается во всей области, однако эта запись лишь формальна, и ядро этого уравнения сосредоточено на указанных локальных областях. Как следствие, в методах П, в отличие от Ш, как правило , . приходится вводить £- границу и сносить граничные условия, т.е. прибегать к аппроксимации граничных условий.

Далее для полноты следует упомянуть еще об одном подходе, широко применяемом в практических расчетах в силу своей простоты и сравнительной универсальности. Речь идет о методе блуждания по решетке. Однако его можно отнести к подходу I, поскольку этот метод представляет собой частный случай дискретного приближения континуальных интегралов.

Безусловно, предлагаемая классификация носит условный характер, и приводится для удобства дальнейшего изложения результатов диссертации.

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

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

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

Перейдем к изложению содержания диссертации по главам.

В диссертации формулы имеют замкнутую нумерацию внутри каждой главы. В тех случаях, когда делается ссылка на формулы из другой главы, указываются номер формулы и номер главы.

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

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

Следующий, третий параграф, посвящен обобщению алгоритма из §1.2 на неоднородные уравнения с нулевыми начальными данными на основе решения задачи Коши или функции Грина с последующим использованием двойной рандомизации.

В четвертом параграфе содержится доказательство конечности дисперсии оценок, предложенных в §1.2, и рассматривается метод установления душ решения стационарных задач, оценивается трудоемкость алгоритма блуждания по границе.

В пятом параграфе первой главы на основе одного метода представления решения однородной задачи Дирихле в виде потенциала простого слоя (см., [31 ] стр. 464), строится оценка как самого решения, так и ее производных. Причем дисперсия этих оценок конечна и в граничных точках, что не справедливо для обычных оценок, рассмотренных в предыдущих параграфах.

Вторая глава посвящена разработке и обоснованию алгоритмов блуждания внутри области, для.решения нестационарных и некоторых стационарных краевых задач.

В первом параграфе этой главы выводится соотношение о среднем для уравнения теплопроводности, которое является обобщением соответствующей теоремы из работы [24] .

Во втором параграфе предлагается модификация алгоритма блуждания по шароидам А.С.Еипина, основанного на соотношении о среднем, полученном в §2Л данной главы. В этой модификации случайное блуждание в ряде случаев заканчивается существенно быстрее.

Третий параграф содержит обобщение вышеупомянутого алгоритма на случай уравнения гц —Ди-*сс*с ограниченным коэффициентом

В четвертом параграфе даются достаточные условия,при которых справедлива логарифмическая оценка для среднего числа блужданий N £ до попадания в £ - окрестность границы для процессов блуждания внутри области.

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

ГГ - М.дис*»-Ь) + ссх,-ии = (I)

В первом параграфе выводится соотношение о среднем для уравнений типа (I).

Во втором параграфе излагаются два алгоритма для решения задачи Коши для уравнения (I).

В третьем параграфе данной главы строится алгоритм блуждания, по сферам для решения первой краевой задачи для уравнения (I).

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

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

В первом параграфе данной главы численно решается одна задача с известным точным решением, методами блуждания по границе и сферам.

Второй параграф четвертой главы посвящен задаче о трансформации спектра полидисперсных аэрозолей при распространении и осаждении в трубе. Эта задача решается методом блуждания по сферам.

Укажем наконец,.что в приложение вынесены некоторые технические подробности, связанные.с оценкой трудоемкости алгоритмов блуждания внутри области.

Настоящая работа выполнена в Новосибирском Государственном Университете им.Ленинского комсомола. Изложенные в ней результаты были представлены и докладывались на У и У1 Всесоюзных совещаниях по методам Монте-Карло в вычислительной математике и математической физике (Новосибирск, 1976, 1979), на Всесоюзном семинаре по диффузии (Новосибирск, 1983) и на семинарах по методам Монте-Карло Вычислительного центра СО АН СССР. Основные результаты диссертации опубликованы в работах [25 - 30, 45].

Пользуясь случаем, автор.выражает искреннюю благодарность своим научным, руководителям д,ф.-^м.н., профессору Г.А.Михайлову и к.ф.-м.н. К.К.Сабельфельду за постоянное внимание и руководство работой.

ЗАКЛЮЧЕНИЕ

Сформулируем основные результаты диссертационной работы:

1. Разработаны и обоснованы алгоритмы случайного блуждания по границе области для решения внутренних и внешних первой, второй и третьей краевых задач для уравнения теплопроводности.

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

3. Получена точная асимптотическая оценка для трудоемкости алгоритмов, основанных на процессах блуждания внутри области для широкого класса областей и процессов.

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

1. Амелин А. Г. Теоретические основы образования тумана при конденсации пара. - М.: Химия, 1972, 304 с.

2. Браун Дж. Методы Монте-Карло В кн.: Современная математика для инженеров. - М.: И.Л., 1959, с. 275-301.

3. Вентцель А.Д. Общие граничные задачи, связанные с диффузионными процессами. Успехи мат. наук, I960, т.15, вкп.2 (92), с. 202-204.

4. Веселовская А.З. Несмещенные оценки.функционалов от реиений некоторых гиперболических уравнений. В сб.: Методы Монте-Карло в вычислительной математике и математической физике. Труды ВЦ СОМ,СССР, Новосибирск, 1976, с. I2I-I28. .

5. Владимиров B.C. Уравнения математической физики. М.: Наука,1981. 512 с.

6. Владимиров B.C. О приближенной вычислении винеровских интегралов. Успехи мат. наук, I960, т.15, вып. 4, с. 129-135.

7. Гелъфанд Н.М., Яглом А.М. Интегрирование в функциональных пространствах и его применение в квантовой физике. Успехи мат. наук,.1956, т.II, вып. I (67), е.77-114.

8. Гельфанд Н.М., Фролов A.C., Ченцов H.H. Вычисление континуальных интегралов методом Монте-Карло. Изв. высш. учебн. заведений. Математика, 1958, № 5(6), с. 32-45.

9. Джетыбаев Е.О. , Сабельфельд К.К. Решение смешанной задачи для уравнений параболического и гиперболического типа методом Монте-Карло. Ж. вычислит.матем. и матем.физ., 1984,т. 24, 15, с. 677-685. . ,

10. Дынкин Е.Б. Марковские процессы. М.: Физматгиз, 1963. 860 с.

11. Дынкин Е.Б,, Шкевич A.A. Теоремы и задачи о процессах Маркова. М.: Наука, 1967. 231 с. V

12. Елепов B.C., Михайлов Г.А. О решении задачи Дирихле для уравнения = моделированием "блужданий по сфере". Ж. вычислит, матем. и матем. физ., 1969, т.9, № 3, с. 647-654.

13. Елепов Б.С., Михайлов Г.А. Алгоритм "блуждания по сферам" для уравнения = Докл. АН СССР, 1973, т.212, М, с. 15-19.

14. Решение краевых задач методом Монте-Карло./ Елепов B.C., Кронберг A.A., Михайлов Г.А., Сабельфельд К.К. Новосибирск: Наука, 1980. - 173 с. .

15. Ерманов С.М. Метод Монте-Карло и смежные вопросы. М.: Наука, 1975. - 471 с.

16. Ермаков С.М., Михайлов Г.А. Статистическое моделирование. М.: Наука, 1982. 295 е. .

17. Ермаков С.М., Некруткин В.В., Сипин.A.C. Случайные процессы . решения, задач математической физики. М.: Наука, 1984. с.

18. Ильин A.M., Калашников A.C., Ошейник O.A. Линейные уравнения второго порядка параболического типа. Успехи мат. наук,1962, т.17, вып. 3, с. 3-147.

19. Карлин С., Стадцен В. Чебышевские системы и их применение в анализе и статистике. М.: Наука, 1976. - 568 с.

20. Кравченко И.И. ,Лехтмахер С.О., Рузер Л.С. Расчет диффузиш-ного осаждения частиц аэрозоля s логарифмически нормальным, распределением частиц по размерам в цилиндрических каналах. Коллоидный,журнал, 1971, т.33, № 6, с.

21. Колмогоров А.Н. üt&c dte алмДлНьеМлн. Wtbkocle+ь in du,haibbtoCit. Ah-hdie^j404, HiS- 4S~3.

22. Кронберг A.A. К решению двух краевых задач методом Монте-Карло. Ж. вычисл. матем. и матем. физ., 1976, т. 16,1. Л I, с. 162-169.

23. Кронберг A.A. Об асимптотике математического ожидания числа шагов ¿-сферического процесса. Ж. вычисл. матем. и матем. физ., 1980, т.20, £2, с. 528-531.

24. Купцов Л. П. Свойство среднего и принцип максимума для параболических уравнений второго порядка. Доклады АН СССР, 1978, т. 242, №3, с. 529-532.

25. Курбанмурадов 0. Эффективность алгоритма "блуждания по решетке". в сб.: Методы Монте-Карло в вычислительной математике .и математической физике. Тр. ВЦ СОАН СССР, Новосибирск, 1976. , с.163-170.

26. Курбанмурадов 0. Оценка математического ожидания числа шагов е- сферического процесса. В сб.: Методы Монте-Карло в вычислительной математике и математической физике. Тр.

27. ВЦ СОАН СССР, Новосибирск, 1979, с. 137-144.

28. Курбанмурадов 0. Решение второй краевой задачи для смешанного уравнения теплопроводности методом Монте-Карло. Вопросы вычислительной и прикладной математики. Тр. Уз ИПО "Кибернетика", Ташкент, 1980, вып. 59.-е. 24-31.

29. Курбанмурадов 0. Методы блуждания по шароидам для решения уравнения теплопроводности. Теория и алгоритмы статистического моделирования. Тр. ВЦ СОАН СССР, Новосибирск, 1984.

30. Курбанмурадов 0. Исследование дисперсии векторного алгоритма "блуждания по сферам" для полигармонического уравнения. Математические и имитационные модели систем. Тр.ВЦ СОАН СССР, 1983, с. 26-32.

31. Курбанмурадов 0., Сабельфельд K.K. Решения многомерных задач теории потенциала методом блуждания по границе. В сб.: Численные методы механики сплошной среды, Новосибирск, 1984, т.15, с. 77-102.

32. Ладыженская O.A., Солонников В.А., Уральцева H.H. Линейные квазилинейные уравнения параболического типа. М.: Наука, 1967. - 736 с.

33. Метод Монте-Карло в атмосферной оптике. / Марчук Г.И., Михайлов Г.А., Назаралиев М.А. и др. Новосибирск: Наука, 1976. - 283 с.

34. Марчук Г.И. Математическое моделирование в проблеме окружаю. щей среды.-.М.: Наука, 1982. 319 с.

35. Михайлов Г.А. Минимаксная теория весовых методов Монте-Карло. -.Ж. вычислит, матем. и матем. физ., 1984, т.24, №9,с. 1294-1302.

36. Михайлов Г.А., Сабельфельд К.К., Ченцов H.H. Векторнот-стохас-тические модели некоторых задач математической,физики. В кн.: Актуальные проблемы прикл. матем. и матем. моделирования. Новосибирск: Наука, 1982, с. 69-82.

37. Михайлов Г.А.Оптимизация векторных алгоритмов метода Монте-. Карло. Докл. АН ССОР, 1981, т. 260, №1, с.26-31.

38. Мишустин Б.А. О решении задачи Дирихле для уравнения Лапласа, методом статистических, испытаний. Ж. вычисл.матем. и матем. физ., 1967, т.7, №5, с. II79-II85. .

39. Мейер П.-А. Вероятность и потенциалы. М.: Мир, 1973, 325 с.

40. Прохоров Ю.В., Розанов Ю.А. Теория вероятностей. Основные понятия, предельные теоремы, случайные процессы. СМБ.-М.: Наука, 1967, 495 с.

41. Петров В.Б. Суммы независимых случайных величин. М.: Наука, 1972.

42. Расулов М.Л. Применение метода контурного интеграла к решению задач для параболических систем второго порядка. М.: Наука, 1975. - 256 с.

43. Сабельфельд К.К. Векторные алгоритмы метода Монте-Карло для решения систем эллиптических уравнений 2-го порядка и уравнения Ламе. Докл. АН СССР, 1982, т.262, №5,с. 1076-1080.

44. Сабельфельд К.К. Решение одной краевой задачи для метагар-монического уравнения методом Монте-Карло. Ж. вычисл. ма- тем.* и матем.физ., 1979, т.19, №4, с. 961-969.

45. Сабельфельд К.К. О приближенном вычислении винеровских.континуальных интегралов методом.Монте-Карло. Ж. вычисл.матем. и матем. физ., 1979, т.19,.ЖЕ, с. 29-43.

46. Сабельфельд К.К., Курбанмурадов 0. Алгоритмы прямого и сопряженного блуждания для решения.задачи о распространении примеси в случайном поле скоростей: Препринт № 506. Новосибирск, 1984. -г 30 с. В надзаг.: ВЦ СОАН СССР.

47. Сабельфельд К.К., Симонов H.A. Алгоритмы, случайного блуждания по границе для решения краевых задач. В сб.: Численные методы,механики сплошной среды, Новосибирск, т.14, №1,с. I16-134. .

48. Сабельфельд К.К., Симонов H.A. Решение пространственных задач теории упругости в детерминированной и стохастической постановках методом Монте-Карло. Докл. АНСССР, 1984,т.275, № 4, с. 802-805.

49. Сабелъфельд K.K. Принцип выметания и усреднения при построении алгоритмов метода Монте-Карло для решения краевых задач. Ж. вычисл. матем. и матем. физ., 1983, т.23, $2,с. 366 379.

50. Сабелъфельд К.К., Цецохо C.B. Решение диффузионных задач методом статистического моделирования лагранжевых тракторий частиц. В сб.: Методы и алгоритмы статистического модели. рования, 1983, с. 140-148.

51. Симонов H.A. Алгоритмы случайного блуждания для решения краевых задач с разбиением на подобласти. В кн.: Методыи.алгоритмы статистического моделирования, Новосибирск, 1983, . с. 48-58.

52. Сипин A.C. О решении задачи Неймана методом Монте-Карло.- В сб.: Методы Монте-Карло в.вычислит, матем. и матем. физ., Новосибирск, 1976, с. 129-135.

53. Сипин A.C. Решение первой краевой задачи для уравнения эллиптического типа методом Монте-Карло. В сб. : Методы Монте-Карло в вычислит, матем. и матем. физике, Новосибирск, 1979, с. II3-II9.

54. Теория.тепломассообмена / Исаев С.И., Коршунов И.А., Ко. фанов В.й. и др. М. : Высшая школа, 1979, 495 с.

55. Татарский В.И. О приближенном вычислении винеровских континуальных интегралов. В сб. Методы Монте-Карло в вычисл.матем. и матем. физ., Новосибирск, 1976, с. 60-74.

56. Хаджи-Шейх А., Спэрроу. Решение задач теплопроводности вероятностными методами. Труды американского общества инженеров-механиков, русский перевод, 1967, т.2, с.1-18.

57. Исследование объемной конденсации при течении парагазовой смеси в трубе /Шутов С.А., Пащенко С.Э., Анкилов А.Н. и др. -Коллоидный журнал, 1983, 2, с. 387-390.

58. Янович Л.А. Приближенное вычисление континуальных интегралов по гауссовым мерам. Минск., Наука и техника, 1976, 383 с.

59. Соил-CuKi Л, ? TtLccLiU&ki N., Ltvy И. W&t cät

60. L^ол,€>п l ¿oh owcfMu c(e>i ynath-mjbl-icAi&K Physik. Haihzmi. A werfe*., 34-W.59. ßoo£ J.L, -ZctvUfnoJvtCHcjalt OmqL SuthoJL^nofUz fundte/t*. Ткомл . Amet. Hatk. Joe., J6-9f.

61. MlkcU M, BvaütatCow Coutitiiuotändkod ty Шсид рковщ. A-hnah7/ tU УнкШс HoJbtiiical Моикк&маЛСс*) ТоКдб^1. M. H} 43-54.61. ilu-f-foz, M.t. eoxtinuOul Ко»Ы'0#Ло tootrhooli {Ol tfuL ¡biiAt^id ptcd&M. Аня. MM. ШЬ('НссЯ)