Латиноамериканская модель глобального развития - Математическое и компьютерное моделирование в естествознании

В 1974г. группа аргентинских ученых во главе с профессором А. Эррерой получила предварительные результаты работы над латиноамериканской моделью глобального развития. Предпосылки для выполнения работы при обсуждении модели "Мир-3" послужил тезис о том, что основные преграды на пути гармонического развития человечества заключаещися главным образом в неравномерном распределении богатства между различными странами.

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

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

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

Компьютерное моделирование в естествознании: возможности, достижения, перспективы

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

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

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

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

Результативность компьютерной модели в значительной степени определяется качеством используемого программного обеспечения. Основные требования, предъявляемые к программам - это, конечно, простота ввода и корректировки исходных данных, а также визуализация (наглядность) результатов счета. Сегодня имеются и мощные специализированные системы программирования (MAPLE, SolidWorks, AutoCAD и др.) и специальные программы, в которых реализуется удобные графические пользовательские возможности.

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

Приведем примеры естественнонаучных задач, смоделированных при помощи компьютера:

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

В начале работы необходимо уточнить формулировку и условия задачи:

По заданным условиям (начальные координаты малого тела, начальная скорость малого тела, масса большого тела) построить траекторию движения малого тела, проверить является ли она эллипсом, если да, то проверить выполнение второго закона Кеплера для данных тел. Считать, что система, в которой производится моделирование, состоит из 2х тел. Считать гравитационную постоянную. Считать центральное тело неподвижным. Производить расчеты для тел, удаленных от центра системы не более чем на 100 а. е. Считать а. е. = 149 597 890 000м, а центром системы - центр большого тела.

Далее опишем методы исследования процесса:

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

(1)

Здесь G = 6,67-10-11 м3/кг-с2) - гравитационная постоянная.

выбор системы координат при решении задачи двух тел

Рис. 1 Выбор системы координат при решении задачи двух тел

Знак "минус" в формуле (1) связан с тем, что гравитационная сила является силой притяжения, т. е. стремится уменьшить расстояние г между телами.

Рассмотрим случай, когда центральным телом является Солнце, масса которого

М=1,9891Ч1030 кг, а обращающимся Земля, т. е. за расстояние примем 1 астрономическую единицу, равную 149 597 890 000м, пренебрегая при этом относительно небольшими силами притяжения от всех прочих небесных тел. Разумеется, мы тем самым произвели ранжирование факторов, и наши последующие действия имеют отношение к реальности лишь в меру соблюдения определенных условий.

Уравнение, описывающее движение тела m в указанной системе координат, имеет вид

Или в проекциях на оси х, у

(2)

Интересующая нас орбита сильно зависит от "начальной скорости" тела т и "начального расстояния. При моделировании нам придется принять некоторое положение условно за начало, а затем изучать движение дальше. Очень часто космические тела движутся практически с постоянной скоростью по орбитам, близким к круговым. Для таких орбит легко найти элементарное соотношение между скоростью и радиусом. В этом случае сила тяготения выступает в роли центростремительной, а центростремительная сила при постоянной скорости выражается известной из начального курса физики формулой mv2/r. Таким образом, имеем

Или

    (3) - искомое соотношение.

Период движения по такой орбите

Если соотношение (3) нарушено, то орбита не будет круговой. Выяснить, какой она будет, можно в ходе численного моделирования. Сведем (2) к системе четырех дифференциальных уравнений первого порядка:

(4).

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

Метод Эйлера дает возможность построить таблично заданную функцию на некотором интервале ее значений. По методу Эйлера напишем систему соотношений (6) для нахождения таблицы значений функций описанных уравнением (3) и удовлетворяющих условию (5). Подробнее процесс построения таблицы значений функций по методу Эйлера описан в [6].

(5)

(6)

Где шаг сетки значений h.

Для проверки закона Кеплера нам необходимо сравнить площади секторов, заметаемых радиус-вектором за одинаковые промежутки времени. Если площади будут равные в достаточной степени точности, то мы сможем утверждать, что второй закон Кеплера выполняется. Для этого возьмем достаточно малое приращение по времени, для того, чтобы площадь данного сектора стремилась к площади треугольника, со сторонами а, b, c (рис.2).

Рис.2

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

(7)

Где p - полупериметр треугольника, изображенного на рис. 2, a, b, c - длины его сторон.

Из системы уравнений (7) мы можем выразить координаты точек эллипса, удовлетворяющих заданной степени точности. Координаты начала векторов a, b будут равны (0,0), т. к. за начало координат принято Солнце. Длины векторов вычисляются по известной из алгебры формуле:

(8)

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

После того как модель готова остается только написать программу, которая будет реализовывать данный алгоритм. Наша программа написана на языке С++.

Входные данные программы:

Mg= 132672970000000000000.0 метров - произведение массы Солнца на константу g.

X = 149600000000.0 - начальная координата малого тела по оси х

Y = 0 - начальная координата малого тела по оси у

Vx = 0 м/с - начальная скорость по оси х

Vy = 29785 м/с - начальная скорость по оси у

Выходные данные программы:

SMax - площадь максимального из заметаемых радиус-вектором секторов

SMin - площадь минимального из заметаемых радиус-вектором секторов

Тестирование:

При заданных входных параметрах, площади треугольников, получаемых с помощью описанного выше алгоритма приближенно равны (не равны точно, т. к. в языке С++ переменные типа float при вычислениях округляются). sMax = 8.02207e+018, sMin = 8.01897e+018. Это позволяет говорить о том, что второй закон Кеплера справедлив.

2. Постановка задачи. Построить изолинии поля, созданного четырьмя одноименными зарядами, находящимися в вершинах прямоугольника. Значения зарядов (при последовательном обходе вершин) есть q, 2q, 3q, 4q.

Для того чтобы понять, какой способ условного изображения поля наиболее удобен и понятен было нам было необходимо понять суть самого предмета электродинамики, и на основе полученных знаний, предположить, какими могут быть альтернативные способы. Для этого нами были прочитаны: главы 36-40 [5], глава 3.8 [3], и раздел "Описание физических процессов в приближении сплошной среды" [4], в последствии чего было выявлено, что наиболее эффективным и наглядным способом является изображение поля линиями равного потенциала.

Описание метода исследования процесса. Для построения изолиний поля, созданного несколькими зарядами, можно воспользоваться принципом суперпозиции, то есть потенциал в каждой точке поля равен, где цI Создаются в этой точке i-м зарядом. Потенциал поля, созданного зарядом Q на расстоянии r от него, равен, где Е0 - электрическая постоянная.

Для решения воспользуемся методом сеток. Пусть поле создается системой зарядов Q1 ... QP с координатами соответственно (x1, y1),...,( xP, yP). Выберем по осям х и у и некоторые шаги hX и hY и покроем плоскость сеткой, образованной прямыми, параллельными осям х и у и отстоящими друг от друга на расстояниях hX и hY. Точки пересечения этих прямых - узлы сетки. Пронумеруем их так: начало координат (0, 0), следующий по оси х вправо - (0, 1), влево - (0, -1); по оси у вверх - (1, 0), вниз - (-1, 0) и т. д. Значения потенциала, создаваемого системой зарядов Q1 ... QР в узле (i, k), согласно принципу суперпозиции, таково (i - номер строки, k - столбца сетки): .

Ограничимся прямоугольной областью в плоскости ху: [-mhX, mhX] по оси х и [-nhY, nhY] по оси у. В этой области (2m + l) - (2n + l) узлов. Вычислим значения потенциала в каждом из них по указанным формулам.

Фиксируем некоторое значение потенциала Ф и построим изолинию, соответствующую этому значению. Для этого проходим, к примеру, по i-ой горизонтальной линии сетки и ищем среди ее узлов такие соседние значения потенциала, причем выполняется. Если такая пара узлов найдена, то координату точки, значение потенциала в которой равна Ф, найдем приближенно с помощью линейной интерполяции:

Найдя в данной горизонтали все такие точки, переходим к следующей горизонтали, пока не исчерпаем их все. Для этого надо совершить двойной циклический проход: во внешнем цикле перебирать i от - п до +п, во внутреннем перебирать k от - т до +т.

После этого надо найти нужные точки на вертикальных линиях сетки. Формулы аналогичны и имеют вид: .

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

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

Реализация программы: входными параметрами нашей программы являются координаты точек, в которых расположены заряды; величины зарядов; значения потенциалов.

Следовательно выходные и выходные данные будут представлены в следующем виде:

Входные данные:

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

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

Массив, который содержит значения зарядов. Обозначим q

Массив, который содержит значения потенциалов. Обозначим g

Выходные данные:

Семейство изолиний.

Далее нам остается лишь корректно реализовать предложенный нами в главе I алгоритм.

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

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

Похожие статьи




Латиноамериканская модель глобального развития - Математическое и компьютерное моделирование в естествознании

Предыдущая | Следующая