Sx формула физика: Формулы кинематики пожалуйста – Школьные Знания.com

Содержание

Задачи на равноускоренное движение | Образовательный портал EduContest.Net — библиотека учебно-методических материалов


Чтобы посмотреть презентацию с картинками, оформлением и слайдами, скачайте ее файл и откройте в PowerPoint на своем компьютере.
Текстовое содержимое слайдов презентации:

Решение задач по теме «Прямолинейное равноускоренное движение» Урок физики в 9 – 10 классе Цели урока: Повторить основные формулы по теме «Прямолинейное равноускоренное движение».Сформировать навыки решения задач по данной теме. Основные формулы: Vx – V0x 1. аx = – ускорение t 2. Vx = Vox + aхt – скорость Vx + Vox 3. Sx = t 2 ax t2 4. Sx = Vox t + перемещение 2 Vx2 –Vox2 5. Sx = 2ax axt2 6. X = Xo + Vox t + – уравнение прямолинейного 2 равноускоренного движения Задача №1.С каким ускорением движется гоночный автомобиль, если его скорость за 6 с увеличивается со 144 до 216 км/ч? Задача №1.С каким ускорением движется гоночный автомобиль, если его скорость за 6 с увеличивается со 144 до 216 км/ч? Дано “СИ” Решение: Vo=144 км/ч 40 м/с V – Vo V = 216 км/ч 60 м/с а = t = 6 с t (60 – 40) м/с а – ? а = = 3,33 м/с2. 6 с Ответ: а = 3,33 м/с2.
км 144 · 1000 м м 144 = = 40 ч 3600 с с км 216 · 1000 м м 216 = = 60 ч 3600 с с Задача №2За какое время ракета приобретает первую космическую скорость 7,9 км/с, если она будет двигаться с ускорением 50 м/с2? Задача №2За какое время ракета приобретает первую космическую скорость 7,9 км/с, если она будет двигаться с ускорением 50 м/с2? Дано: “СИ” Решение. V = 7,9 км/с 7900 м/с V – Vo V Vo= 0 а = , т.к. Vo = 0, то а = а = 50 м/с2 t t V t = .
a 7900 м/с t = = 158 с. 50 м/с2 Ответ: t = 158 с. t – ? Задача №3Рассчитайте длину взлетной полосы, если скорость самолета 300 км/ч, а время разгона 40 с. Задача №3Рассчитайте длину взлетной полосы, если скорость самолета 300 км/ч, а время разгона 40 с. Дано: “СИ” Решение. V = 300 км/ч 83,3 м/с V + Vo Vo = 0 S = t t = 40с 2 (83,3 + 0) м/с S – ? S = · 40 с = 1666 м 2 Ответ: S = 1666 м ≈ 1,7 км.
Задача №4Скорость гоночного автомобиля в момент начала разгона 10 м/с, ускорение 5 м/с2. Определите путь, пройденный автомобилем за 10 с после начала движения. Какова скорость автомобиля в конце десятой секунды разгона? Задача №4Скорость гоночного автомобиля в момент начала разгона 10 м/с, ускорение 5 м/с2. Определите путь, пройденный автомобилем за 10 с после начала движения. Какова скорость автомобиля в конце десятой секунды разгона? Дано: Решение. Vo= 10 м/с a t2 5 м/с2 · (10 с)2 а = 5 м/с2 S = Vot + ; S = 10 м/с · 10 с + = 350 м. t = 10 с 2 2 S – ? V = Vo + a t ; V = 10 м/с + 5 м/с2 · 10 с = 60 м/с. V – ? Ответ: S = 350 м; V = 60 м/с. Задача №5Тормозной путь автомобиля, движущегося со скоростью 50 км/ч, равен 10 м.
Чему равен тормозной путь этого же автомобиля при скорости 100 км/ч? Задача №5Тормозной путь автомобиля, движущегося со скоростью 50 км/ч, равен 10 м. Чему равен тормозной путь этого же автомобиля при скорости 100 км/ч? Дано: “СИ” Решение. V = 0 Vo12 Vo12 Vo1 = 50 км/ч 13,9 м/с S1 = a = Vo2 = 100 км/ч 27,8 м/с 2a 2S1 S1 =10 м Vo22 Vo22 · 2S1 Vo22 S2 = = = S1 S2 – ? 2a 2 Vo12 Vo12 (27,8 м/с)2 772,84 S2 = 10 м —————– = 10 ————- = 40 м.
(13,9 м/с)2 193,21 Ответ: S2 = 40 м. Задача №6Какова длинна пробега самолета при посадке, если его посадочная скорость 140 км/ч, а ускорение при торможении 2 м/с2? Задача №6Какова длинна пробега самолета при посадке, если его посадочная скорость 140 км/ч, а ускорение при торможении 2 м/с2? Дано: “СИ” Решение. Vо = 140 км/ч 38,9 м/с V2 – Vo2 а = 2 м/с2 S = ; ax = – 2 м/с2. V = 0 2 ax S – ? Vo2 S = 2а км 140 · 1000 м м 140 = = 38,9 (38,9 м/с)2 ч 3600 с с S = ≈ 378 м.
2 · 2 м/с2 Ответ: S = 378 м. Задача №7Автомобиль, имея начальную скорость 54 км/ч, при торможении по сухой дороге проходит 30 м, а по мокрой – 90 м. Определите для каждого случая ускорение и время торможения. Задача №7Автомобиль, имея начальную скорость 54 км/ч, при торможении по сухой дороге проходит 30 м, а по мокрой – 90 м. Определите для каждого случая ускорение и время торможения. Дано: “СИ” Решение. V = 0 V + Vo Vot 2 SVo = 54 км/ч 15 м/с S = t S = t = . S1 = 30 м 2 2 Vo S2 = 90 м t1 = ? t2 = ? а – ? t – ? Vo2 Vo2 S = a = .
2 a 2 S a1 = ? a2 = ? Ответ: a1 = 3,75 м/с2; t1 = 4 с; a2 = 1,25 м/с2; t1 = 12 с. Задача №8При равноускоренном движении с начальной скоростью 5 м/с тело за 3 с прошло 20 м. С каким ускорением двигалось тело? Какова его скорость в конце третьей секунды? Задача №8При равноускоренном движении с начальной скоростью 5 м/с тело за 3 с прошло 20 м. С каким ускорением двигалось тело? Какова его скорость в конце третьей секунды? Дано: Решение. Vо = 5м/с a t2 a t2 2(S – Vot) t = 3 c S = Vot + = S – Vot a = . S = 20 м 2 2 t2 V = Vо + at .
a – ? V – ? 2 · (20 м – 5 м/с · 3с) a = ≈ 1,1 м/с2; 9 с2 V = 5 м/с + 1,1 м/с2 · 3 с = 8,3 м/с. Ответ: а = 1,1 м/с2; V = 8,3 м/с. Задача №9Два велосипедиста едут навстречу друг другу. Первый, имея начальную скорость 9 км/ч, спускается с горы с ускорением 0,4 м/с2. Второй поднимается в гору с начальной скоростью 18 км/ч и ускорением 0,2 м/с2. Через какое время встретятся велосипедисты, если начальное расстояние между ними 200 м? Задача №9Два велосипедиста едут навстречу друг другу.
Первый, имея начальную скорость 9 км/ч, спускается с горы с ускорением 0,4 м/с2. Второй поднимается в гору с начальной скоростью 18 км/ч и ускорением 0,2 м/с2. Через какое время встретятся велосипедисты, если начальное расстояние между ними 200 м? Дано: “СИ” Решение.Vo1x=9 км/ч 2,5 м/с а1х= 0,4 м/с2Vo2x= – 18 км/ч – 5 м/са2х= – 0,2 м/с2Хо2= 200 мХо1 = 0 м а1х t2t – ? Х1 = Хо1 + Vo1xt + 2 а2х t2 Х2 = Хо2 + Vo2xt + 2Место встречи Х1 = Х2 0,4 t2 0,2 t2 Х1 = 2,5 t + и X2 = 200 – 5 t – или 2 2 2,5 t + 0,2 t2 = 200 – 5 t – 0,1 t2 0,3 t2 + 7,5 t – 200 = 0 t ≈ 16,2 c. 200 м Х О Ответ: t = 16,2 c Задача №10Уравнение координаты имеет вид Х = 4 + 1,5t + t2. Какое это движение? Напишите формулу зависимости скорости тела от времени. Чему равны скорость и координата тела через 6 с? Задача №10Уравнение координаты имеет вид Х = 4 + 1,5t + t2. Какое это движение? Напишите формулу зависимости скорости тела от времени. Чему равны скорость и координата тела через 6 с? Дано: Решение. х = 4 + 1,5t + t2 Запишем уравнение равноускоренного движения в t = 6c общем виде: а t2V -? Х = Хо + Voxt + X -? 2 Сравним с данным уравнением: х = 4 + 1,5t +1t2 Х0 = 4 м а Vox = 1,5 м/с = 1 а = 2 м/с2 > 0 2 движение равноускоренное Запишем уравнение скорости: V = Vo + a t V = 1,5 + 2 t Вычисляем: V = 1,5 м/с + 2 м/с2 · 6с = 13,5 м/с. Х = 4м + 1,5 м/с · 6 с + 1м/с2 (6 с)2 = 49 мОтвет: V = 1,5 + 2 t ; V = 13,5 м/с; Х = 49 м. Желаем успеха в самостоятельном решении задач! Мякишев Г.Я., Буховцев Б.Б.,Сотский Н.Н. Физика 10 класс. – М.: Просвещение, 2007. – 365 с. Касьянов В.А. Физика 10 класс. – М.: Дрофа, 2006. – 410 с. Волков В.А. Поурочные разработки по физике. 10 класс. – М: Вако, 2006. – 400 с. Касаткина И.Л., Ларцева Н.А., Шкиль Т.В. Репетитор по физике. В 2-х томах. Том 1. – Ростов-на-Дону: Феникс, 1995. – 863 с.5. Рымкевич А.П. Задачник 10 – 11 классы. – М.: Дрофа, 2004. – 188 с.6. Степанова Г.Н. Сборник задач по физике 10 – 11 классы. – М: Просвещение, 2003. – 287 с. Литература

Приложенные файлы

Педагогическая мозаика. Персональный сайт преподавателя

Механика


1. Кинематика
2. Динамика
3. Законы сохранения

Кинематика

1.5 Прямолинейное равномерное  движение

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

Всплывают пузырьки в воде

 

Вращается планета 

Девочка едет на роликах

Скачет баскетбольный мяч


Посмотрите на рисунки и ответьте на вопрос: какие тела движутся прямолинейно и равномерно?

Для характеристики быстроты движения тела существует понятие скорости.

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




Отсюда уравнение прямолинейного движения выглядит так:Прибор для измерения скорости называется  

Определите  в каждом случае по спидометру скорость в км/ч.

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

Проекция скорости Vx может быть как положительной, так и отрицательной –  в зависимости от того, в каком направлении оси х  движется тело.


Перемещение при прямолинейном равномерном движении:  Sx = Vxt,

Координата при прямолинейном равномерном движении:    Х = Хо + Vxt.



Мотоциклист движется равномерно со скоростью 36 км/ч. Запишите уравнение координаты для мотоциклиста в СИ. И определите его координату через 5 с.

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

Так как таких тестов  в дальнейшем будет много, то вам нужно скачать и установить на своем компьютере программу MyTest X с сайта mytest.klyaksa.net.После установки программы скачайте  Тест №1 и пройдите тестирование.

                           

Как найти перемещение в физике формула?

Отсюда следует, что проекция перемещения выражается формулой: sx = (vx + v0x)t. При равноускоренном прямолинейном движении скорость тела в любой момент времени равна vx = v0x + axt, следовательно, sx = (2v0x + axt)t.

Как найти перемещение в физике формулы?

Перемещение s→- перемещение; вектор, соединяющий начальную и конечную точки движения. Проекция вектора перемещения на данную ось: sx = x — x0. Путь (l) — длина траектории.

Как найти перемещение по уравнению?

Уравнение перемещения — зависимость перемещения тела от времени s = s(t).

Как найти T по физике?

Время обозначается как t. Единица измерения времени – с (секунды). Самая простая формула при равномерном прямолинейном движении. Время, необходимое для прохождения пути равняется частному от деления пути на скорость равномерного прямолинейного движения: t = S / v.

Как найти скорость равномерного движения формула?

Формула выглядит следующим образом: s = v 0 t + a t 2 2 , где а — это ускорение.

Как найти путь формула?

1) Формула пути:

  1. Чтобы найти расстояние, нужно скорость умножить на время.
  2. Чтобы найти время, надо расстояние разделить на скорость.
  3. Чтобы найти скорость, надо расстояние разделить на время.

Как найти перемещение при равноускоренном движении?

равноускоренном движении: Sх = V0x t + ахt2/2.

Как найти уравнение зависимости координаты от времени?

зависимость координаты движущегося тела от времени имеет вид: x = x 0 + v 0 x t + a x t 2 2 . Последняя формула выражает кинематический закон равнопеременного движения.

Как найти перемещение тела при равномерном прямолинейном движении?

При прямолинейном равномерном движении модуль перемещения равен пройденному пути: |vec{s}|=s, если направление скорости совпадает с направлением вектора перемещения. Модуль вектора перемещения s, совершённого телом за промежуток времени t 1 , в данном случае определяется по формуле: s = v 1 ⋅ t 1 .

Как найти перемещение без начальной скорости?

Перемещение при прямолинейном равноускоренном движении без начальной скорости Проекция вектора перемещения при прямолинейном равноускоренном движении рассчитывается по следующей формуле: Sx=V0x*t+(ax*t^2)/2.

Как определить координаты тела в любой момент времени?

х=х0+vхt. Это уравнение есть уравнение равномерного прямолинейного движения точки, записанное в координатной форме. Оно позволяет найти координату х тела при этом движении в любой момент времени, если известны проекция его скорости на ось ОX и его начальная координата х0.

Как вычислить период ротации?

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

Что бы найти массу?

Масса тела равна произведению плотности вещества на объем этого тела: m = ρ · V. Чтобы вычислить объем тела, нужно массу тела разделить на его плотность: V = m : p.

Перемещение при прямолинейном равноускоренном движении

 

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

Для этого обратимся к графику зависимости проекции скорости прямолинейного равноускоренного движения от времени.

График зависимости проекции скорости прямолинейного равноускоренного движения от времени

Ниже на рисунке представлен график, для проекции скорости некоторого тела, которое движется с начальной скорость V0 и постоянным ускорением а.

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

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

Найдем площадь фигуры ограниченную осью оt, отрезками АО и ВС, а также отрезком АС.

 

Выделим на оси ot малый промежуток времени db. Проведем через эти точки перпендикуляры к оси времени, до их пересечения с графикос проекции скорости. Отметим точки пересечения a и c. За этот промежуток времени скорость тела поменяется от Vax до Vbx.

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

Тогда можно считать отрезок ac горизонтальным, а abcd – прямоугольником. Площадь abcd будет численно равна проекции вектора перемещения, за промежуток времени db. Мы можем разбить на такие малые промежутки времени всю площадь фигуры OACB. 

То есть мы получили, что проекция вектора перемещения Sx за промежуток времени, соответствующий отрезку ОВ, будет численно равна площади S трапеции ОACB, и будет определяться по той же формуле, что и эта площадь.

Следовательно,

Так как Vx=V0x+ax*t и S=Sx, полученная формула примет следующий вид:

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

В случае равнозамедленного движения формула примет следующий вид:

Нужна помощь в учебе?



Предыдущая тема: Скорость прямолинейного равноускоренного движения: график скорости
Следующая тема:&nbsp&nbsp&nbspПеремещение при прямолинейном равноускоренном движении без начальной скорости

Движение с постоянным ускорением

При равноускоренном движении справедливы следующие уравнения, которые мы приводим без вывода:

      =   o t  +  ½  t²      ⇒sx  =  υox t  +  ½ ax t²
sy  =  υoy t  +  ½ ay t²

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

Зная, что  sx = x – xo  и   sy = y – yo  (см. § 12-е), из двух скалярных формул из правой верхней колонки получим уравнения для координат:

x  =  xo + υox t  + ½ ax t²    и y  =  yo + υoy t  + ½ ay t²

Поскольку ускорение при равноускоренном движении тела постоянно, то координатные оси всегда можно расположить так, чтобы вектор ускорения был направлен параллельно одной оси, например оси Y. Следовательно, уравнение движения вдоль оси X заметно упростится:

x  =  xo + υox t  + ( 0 )       и       y  =  yo + υoy t  + ½ ay t²

Обратите внимание, что левое уравнение совпадает с уравнением равномерного прямолинейного движения (см. § 12-ж). Это означает, что равноускоренное движение может «складываться» из равномерного движения вдоль одной оси и равноускоренного движения вдоль другой. Подтверждением этому служит опыт с ядром на яхте (см. § 12-б).

Задача. Вытянув руки, девочка подбросила шар. Он поднялся на 80 cм и вскоре упал к ногам девочки, пролетев 180 cм. С какой скоростью шар был подброшен и какую скорость шар имел при ударе о землю?

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

Возведём в квадрат обе части уравнения для проекции на ось Y мгновенной скорости: υy  =  υoy + ay t  (см. § 12-и). Получим равенство:

υy²  =  ( υoy + ay t )²  =  υoy² + 2 υoy ay t + ay² t²

Вынесем за скобки множитель  2 ay  только для двух правых слагаемых:

υy²  =  υoy² + 2 ay ( υoy t + ½ ay t² )

Заметим, что в скобках получилась формула для вычисления проекции перемещения:  sy = υoy t + ½ ay t².  Заменяя её на sy , получим:

υy²  =  υoy²  + 2 ay sy

Решение. Сделаем чертёж: ось Y направим вверх, а начало координат поместим на земле у ног девочки. Применим выведенную нами формулу для квадрата проекции скорости сначала в верхней точке подъёма шара:

0 = υoy² + 2·(–g)·(+h)     ⇒     υoy = ±√¯2gh = +4 м/с

Затем при начале движения из верхней точки вниз:

υy² = 0 + 2·(–g)·(–H)     ⇒     υy = ±√¯2gh = –6 м/с

Ответ: шар был брошен вверх со скоростью 4 м/с, а в момент приземления имел скорость 6 м/с, направленную против оси Y.

Примечание. Надеемся, вы понимаете, что формула для квадрата проекции мгновенной скорости будет верна по аналогии и для оси X:

υx²  =  υox²  + 2 ax sx

Если движение одномерное, то есть происходит только вдоль одной оси, можно пользоваться любой из двух формул в рамках.

Графики равномерного, равноускоренного движения, сравнение. Линейная, квадратная зависимость. Правила определения параметров

Тестирование онлайн

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

Графики равномерного движения

Зависимость ускорения от времени. Так как при равномерном движении ускорение равно нулю, то зависимость a(t) – прямая линия, которая лежит на оси времени.

Зависимость скорости от времени. Скорость со временем не изменяется, график v(t) – прямая линия, параллельная оси времени.

Правило определения пути по графику v(t): Численное значение перемещения (пути) – это площадь прямоугольника под графиком скорости.

Зависимость пути от времени. График s(t) – наклонная линия.

Правило определения скорости по графику s(t): Тангенс угла наклона графика к оси времени равен скорости движения.

Графики равноускоренного движения

Зависимость ускорения от времени. Ускорение со временем не изменяется, имеет постоянное значение, график a(t) – прямая линия, параллельная оси времени.

Зависимость скорости от времени. При равномерном движении путь изменяется, согласно линейной зависимости . В координатах . Графиком является наклонная линия.

Правило определения пути по графику v(t): Путь тела – это площадь треугольника (или трапеции) под графиком скорости.

Правило определения ускорения по графику v(t): Ускорение тела – это тангенс угла наклона графика к оси времени. Если тело замедляет движение, ускорение отрицательное, угол графика тупой, поэтому находим тангенс смежного угла.

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

График движения при . График движения при

График движения при . График движения при

Сравнительная таблица графиков

Урок-соревнование “Восхождение на пик Механики”. 10-й класс

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

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

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

Игровое поле состоит из красочного планшета, на котором изображен горный пейзаж с нанесенным на него маршрутом восхождения и привалами. Привалы пронумерованы – их 10; старт обозначен флажком. Сбоку на планшете находятся карманы (тоже пронумерованы), в которых помещены карточки с заданиями для каждого привала. На некоторых привалах команду ожидает сюрприз – неудача. Так на карточке написано: “Туман, снегопад, команде вернуться на базу” или “Ожидается сход лавины, срочно спуститься на один переход”.

Класс делится на две команды: “Эрудит” и “Квант”, и капитаны команд, которые представляют свои команды и принимают важные решения (замену игроков или задания выполняют сами).

“Эрудит” – эти ребята умные, достойные и талантливые. “Квант” команда веселых, активных, находчивых, толковых.

Игра проводится под девизом: “Один за всех и все за одного”

Пусть кипит работа
Сильней соревнованье.
Успех решает не судьба,
А ваше знанье!
Мы, играя, проверяем,
Что умеем и что знаем.

Урок проводится с привалами:

1. Блиц-турнир

2. Найди ошибку(формулы)

3. Реши задачу

4. Поиск(экспериментальные задачи)

5. Графики

6. Задачи-шутки

7. Погружение

8. Эрудит

9. Средняя скорость

10. Тесты

11. Подведение итогов урока

На каждом привале выполняется задание участниками. Если участники выполняют задание правильно, то получают “5”, если с помощью команды, то получают жетон. Команда, выполнившая задание первой – тоже жетон. Побеждает та команда, которая первой доставит флажок на пик.

Учитель: “Механика – древнейшая наука, но не устаревшая”

Рекорды Гинесса на уроках физики “Основы кинематики”

Самые быстрые и самые медленные птицы

  1. Наиболее быстрыми из летающих птиц считаются утки и гуси, а такие виды, как красногрудый крохаль, гага, нырок североамериканский и гусь со шпорами, могут иногда развивать скорость до 108 км/ч, спасаясь от опасности.
  2. Белогрудый стриж, живущий в районах Азии, и альпийский стриж развивают чрезвычайно высокую скорость во время брачных танцев в воздухе и боев с соперниками – до 170 км/ч.
  3. Однако можно сказать, что в прямом полете скорость у половины летающих птиц не превышает 64 км/ч.
  4. Минимальную скорость в горизонтальном полете, не превышающую 8 км/ч, имеет американский лесной петух.

Самые быстрые из летающих насекомых.

В 1926 г. американским ученым было сделано сообщение, что он наблюдал полет оленьего овода со скоростью 1316 км/ч на высоте 3657 м. Однако это маловероятно.

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

Современные исследования показали, что средняя скорость полета любого насекомого (включая оленьего овода) составляет 39 км/ч, а максимальная – на коротких дистанциях – 58 км/ч.

1. Блиц-турнир (капитаны)

1. Сколько м в 1 км?

2. Сколько с в 1 часе?

З. Что больше: 10 м/с или 36 км/ч?

4. Ускорение: величина векторная или скалярная?

5. Чему равна длина траектории?

6. Чему равно ускорение свободного падения?

7. Как называется движение, если скорость не меняется?

8. Что показывает счетчик пути в автомобиле: путь или перемещение?

9. Чему равна скорость, если тело покоится?

10. Чему равно ускорение, если тело движется равномерно?

2. Формулы (исправить ошибку-1 человек из команды)

Не зная формул, пропадешь,
Когда в мир знаний ты войдешь.

3. Реши задачу. (1 человек из команды)

Энрико Ферми: “Знать физику – это значит уметь решать задачи”.

Автомобиль “Москвич”, двигаясь из состояния покоя, через 250 м развивает скорость 20 м/с. С каким ускорением он движется? Сколько времени длится разгон?

Дано: Решение:
V0x=0 Sx=(Vx2-V0x2)/2ax
Sx=250 м
Vx=20 м/с ax=(Vx2-V0x2)/2Sx
ах – ? ах=(400-0)/2*250=400/500=0,8
t – ? Sx=axt2/2
t=25 с

Ответ: 0,8м/с; 25 с.

4. Поиск (экспериментальные задачи – команда)

1. Мальчик, стоя на балконе пятиэтажного дома уронил орех и решил определить ускорение свободного падения. И очень скоро определил. Как?

(Ответ: измерив высоту и время падения по формуле S=gt2/2).

2.Скала наконец была покорена. Робинзон стоял на ее вершине, внизу шумело море. “Какова высота этой скалы? Надо ее определить”, – подумал он. Как?

(Ответ: измерив t по часам падения камня до уровня моря и по формуле S=gt2/2 вычислить высоту скалы).

3.С помощью установки (желоб, шарик, измерительная лента, часы) определить ускорение, с которым скатывается шарик?

a = 2S/t2

5. Графики. (команда)

Задание:

а) обозначить оси;

б) записать формулы

6. Задачи – шутки (капитаны)

1. Белку с орехами в лапках посадили на гладкий горизонтальный стол и толкнули вдоль него. Приблизившись к краю стола, белка почувствовала опасность. Она знала законы движения Ньютона и, пользуясь одним из них, предотвратила свое падение на пол. Каким образом?

(ответ: выбросила вперед орехи)

2. “Если он не погасил фары, то я не погашу свои”, – так рассуждают некоторые водители. Аналогию с каким законом механики здесь можно провести?

(ответ: с третьим законом Ньютона.Действие равно противодействию)

3. Часто на борту грузовика имеется надпись: “Не уверен – не обгоняй”. Имеет ли эта фраза физический смысл.

7. Комбинированные задачи “Погружение” (1 человек из команды)

1. Автомобиль, массой 6 т проезжает с выключенным мотором до полной остановки 15 м, под действием силы тормозящей 300Н. Сколько времени длилось торможение?

2. Велосипедист, прилагая силу в 5 Н, до полного разгона проезжает путь 50 метров. Какую скорость он имеет в конце пути, если его масса вместе с велосипедом 100 кг?

8. Эрудит (команда)

1. Выехав рано из города на ровное и пустынное шоссе, шофер решил устроить первую остановку ровно через час. Как ему выполнить свое намерение, не глядя на часы? Радиоприемник в машине отсутствует.

(ответ: по спидометру и счетчику пути, то есть шофер должен поддерживать скорость постоянной и дождаться увеличения километража на величину, численно равную этой скорости).

2. “Мировой рекорд” по прыжкам в высоту среди животных принадлежит маленькой южноафриканской антилопе. На какую высоту прыгнет антилопа, если она отталкивается от Земли вверх со скоростью 12 м/с? (ответ: 7,2 м)

3. Установлено, что человек может безопасно для жизни прыгнуть с высоты 2 м. С какой скоростью он достигнет при этом земли? (ответ: 6,3 м/с)

9. Средняя скорость

(задачи барона Мюнхаузена – капитаны)

Когда моя лошадь подворачивает ногу, я взваливаю лошадь на плечо, и мы продолжаем движение, но медленнее: когда я на лошади, мы движемся со скоростью Vi=120 км/ч, а когда лошадь на мне, со скоростью V2=30 км/ч. Чему равна наша средняя скорость, если:

1. Я еду полпути, а потом несу лошадь.

2. Я еду половину времени, а потом несу лошадь?

Тесты 10 класс Вариант I (1 чел.)

10. Основы кинематики

1. Какая единица длины принята основной в Международной системе?

а) 1с; б) 1 мин; в) 1ч; г) 1 сут; д) 1 год

2. Какие из перечисленных ниже величин векторные?

1) скорость 2) ускорение 3) путь

а) только 1; б) только 2; в) только 3; г) 1 и 2; д) 1 и 3; е) 1,2,3

3. В какой из двух задач, приведенных ниже, можно считать шар материальной точкой?

1) измерить время свободного падения шара радиусом I см с высоты 100м.

2) рассчитать архимедову силу, действующую на этот шар, погруженный в воду.

а) только в первой задаче;

б) только во второй задаче;

в) в обеих задачах;

г) ни в первой, ни во второй задачах.

4.По графику, представленному на рисунке, определите скорость движения велосипедиста через 3 с после начала движения.

а) 3 м/с; б) 9 м/с; в) 0; г) 27 м/с; д) среди ответов а - г нет правильного.

Тесты 10 класс Вариант II (1 чел).

Основы кинематики

1. Какая единица длины принята основной в Международной системе?

а) 1 мм; б) 1 см; в) 1м; г) 1км; д) 300000км.

2. Какие из перечисленных ниже величин скалярные:

1) скорость 2) путь 3) перемещение

а) только 1; б) только 2; в) только 3; г) 1 и 2; д) 2 и 3; е) 1,2,3.

3. В какой из двух задач, приведенных ниже, можно рассматривать Землю как материальную точку?

1) рассчитать период обращения Земли вокруг Солнца.

2) рассчитать линейную скорость движения точек поверхности Земли в результате ее суточного вращения.

а) только в первом случае; б) только во втором случае; в) в обоих случаях; г) ни в первом, ни во втором случаях.

4. По графику, представленному на рисунке, определите скорость движения велосипедиста через 2 с после начала движения.

а) 2 м/с; б) 6 м/с; в) 3 м/с; г) 12 м/с; д) среди ответов а – г нет правильного.

11. Подведение итогов урока.

Физика! Какая емкость слова!

Физика – для нас не просто звук.

Физика – опора и основа всех без исключения наук.

12. Домашнее задание.

В одной системе координат постройте по два графика скорости:

1) а) V0x=2 м/с, аx=1,5 м/с2

б) V0x=4 м/с, аx= – 2 м/с2

2) а) V0x=3 м/с, аx=1,5 м/с2

б) V0x=5 м/с, аx= 0 м/с2

Приложение.

Как рассчитать линейность | Sciencing

Способность вычислять линейность (или корреляцию, как ее часто называют) – очень ценный навык. Линейность – это количественная оценка того, насколько сильно связан набор данных. Линейность варьируется от 0 (совсем не связана) до 1 (полностью связана) и дает полезный числовой индикатор, который можно использовать вместе с числовым графиком. Для наших расчетов будет использоваться следующая выборка пар (x, y): x: 2.4, 3.4, 4.6, 3.7, 2.2, 3.3, 4.0, 2.2). Вы получите ответ 0,5124, который мы называем Sy.

Вычисление Sxy

    Вычислите x_y, умножив каждое значение x на соответствующее значение y. Ваши значения x_y будут 3,192, 7,208, 8,280, 6,105, 4,400, 5,808, 8,440, 3,423.

    Сложите вместе все ваши значения x_y, и вы получите сумму (x_y) = 46,856.

    Умножьте сумму (x) на сумму (y), и вы получите ответ 370,08.

    Разделите 370,08 на 8 (общее количество пар данных в нашем примере данных).Вы получите ответ 46,26.

    Вычтите 46,26 из суммы (x * y) (из шага 2), и вы получите ответ 0,5960, который мы называем Sxy.

Объединяем

    Извлеките квадратный корень из Sx, и ответ будет 2,398.

    Извлеките квадратный корень из Sy, и ответ будет 0,716.

    Умножьте свои ответы из шагов 1 и 2, и вы получите ответ 1,717.

    Разделите Sxy на 1,717 (из шага 3), чтобы рассчитать окончательную линейность 0.347. Такой низкий уровень линейности предполагает, что данные слабо связаны и лишь немного линейны.

Формула кинематических уравнений

Кинематика – это исследование движущихся объектов и их взаимосвязей. Есть четыре (4) кинематических уравнения, которые относятся к смещению D, скорости v, времени t и ускорению a.

a) D = v i t + 1/2 при 2 b) (v i + v f ) / 2 = D / t

c) a = (v f – v i ) / t d) v f 2 = v i 2 + 2aD

D = смещение

a = ускорение

t = время

v f = конечная скорость

v i = начальная скорость

Формула кинематических уравнений.

1) Боб едет на велосипеде в магазин со скоростью 4 м / с, когда перед ним выбегает кошка. Он быстро тормозит до полной остановки, с ускорением – 2м / с 2 . Какое у него перемещение?

Ответ: Поскольку Боб остановлен, конечная скорость v f = 0. Его начальная скорость v i = 4 м / с. Ускорение, a = -2 м / с 2 . Время не указано, поэтому используйте уравнение (d) для смещения, D, потому что оно не зависит от времени.

v f 2 = v i 2 + 2aD

(0) 2 = (4 м / с) 2 +2 (- 2 м / с 2 ) D

0 = 16 м 2 / с 2 + (- 4 м / с 2 ) D

-16 м 2 / с 2 = (- 4 м / с 2 ) D

16 м 2 / с 2 = 4 м / с 2 ) D

(16 м 2 / с 2 ) / (4 м / с 2 ) = D

Водоизмещение полное 4 м.

2) Вы путешествуете с постоянной скоростью 11 м / с в течение 5 минут. Как далеко вы уехали?

Ответ: При постоянной скорости v i = v f = 11 м / с. Время t = 5 мин или t = (60 сек / мин x 5 мин) = 300 сек. Теперь используйте уравнение (b), чтобы найти смещение D.

(v i + v f ) / 2 = D / t

D = [(v i + v f ) / 2] t

D = [(11 м / с + 11 м / с) / 2] x 300 с

D = (22 м / с) / 2 x 300 с

D = 11 м / с x 300 с

D = 3300 м. Водоизмещение полное 3,300 м.

3) Каково ускорение автомобиля, который разгоняется с 11 до 40 м / с за 10 секунд?

Ответ: V i = 11 м / с. V f = 40 м / с. Время, t = 10 с. Используйте кинематическое уравнение c), чтобы найти ускорение.

a = (v f – v i ) / t

a = (40 м / с – 11 м / с) / 10 с

a = (29 м / с) / 10 с = 2,9 м / с 2

4) Если автомобиль разгоняется на 3.0 м / с 2 от полной остановки, сколько времени потребуется, чтобы проехать 3000 м?

Ответ: Ускорение a = 2,9 м / с 2 и перемещение D = 3000 м. Автомобиль был неподвижен, поэтому v i = 0. Используйте уравнение a), чтобы найти время.

D = v i t + 1/2 при 2

3000 м = 0т + 1/2 (3,0 м / с 2 ) т 2

3000 м = 1/2 (3,0 м / с 2 ) / т 2

3000 м / 1.5 м / с 2 = t 2

2000 с 2 = t 2

t = 44,72 с

Учебное пособие по линейной регрессии в Excel

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


Пример данных.

Допустим, у нас есть набор данных, , показано слева. Если у нас есть основания полагать, что существует линейная зависимость между переменными x и y , мы можем построить данные и провести по ним «наиболее подходящую» прямую . Конечно, это соотношение регулируется знакомым уравнением . Затем мы можем найти угол наклона , м, , и , пересечение по оси Y, b , для данных, которые показаны на рисунке ниже.

Давайте введем вышеуказанные данные в электронную таблицу Excel, Постройте данные, создайте линия тренда и отображение ее наклона, пересечение по оси Y и значение R-квадрат. Напомним, что значение R-квадрат – это квадрат коэффициента корреляции. (Наиболее статистические тексты показывают коэффициент корреляции как « r », но Excel показывает коэффициент как « R ». Пишете ли вы как r или R , коэффициент корреляции дает нам меру надежности линейная зависимость между значениями x и y .(Значения близкое к 1 означает отличную линейную надежность.))

Введите свои данные, как мы это делали в столбцах B и C. Причина этого строго косметический, как вы скоро увидите.


Уравнения линейной регрессии.

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

Учитывая набор данных с n точками данных, наклоном, пересечением по оси y и коэффициентом корреляции, r , можно определить с помощью следующего:


(Обратите внимание, что пределы суммирования: от i до n , а индексы суммирования для x и y опущены.)

Неявное применение регрессии к выборке данных.

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

  1. Сначала добавьте три столбца, которые будут использоваться для определения количества xy , x 2 и y 2 , для каждой точки данных.
  2. Далее используйте Excel, чтобы оценить следующее: Sx , Sy , S (ху) , S (х 2 ) , S (y 2 ) , (Sx) 2 , (Sy) 2 . Напомним, что символ S означает «суммирование». Кроме того, термин xy является произведением x и y , то есть: x * y .Также термин S (x 2 ) сильно отличается от термина (Sx) 2 . Будьте осторожны с вашим заказом операций!
  3. Теперь используйте Excel для подсчета количества точек данных, n . (Сделать это, используйте функцию Excel COUNT (). Синтаксис COUNT () в этом примере: = СЧЁТ (B3: B8) и отображается в строке формул на снимке экрана ниже.
  4. Наконец, используйте указанные выше компоненты и уравнения линейной регрессии. данные в предыдущем разделе для расчета наклон (м) , пересечение оси y (b) и коэффициент корреляции (r) данных.Если вы будете осторожны, ваш таблица должна выглядеть как наша. Обратите внимание, что наши уравнения для наклона, Y-точка пересечения и коэффициент корреляции выделены желтым цветом.

Линейная регрессия со встроенными функциями.

Очевидно, что вычисленные значения наклона и точки пересечения по оси Y с использованием методов линейной регрессии идентичны значениям более знакомая линия тренда из графика в первом разделе; а именно т = 0.5842 и b = 1,6842. Кроме того, Excel можно использовать для отображения значения R-квадрат. Опять таки, 2 = 2 . Из графика видим, что R 2 = 0,9488. Из нашего анализа линейной регрессии мы находим, что r = 0,9741, поэтому r 2 = 0,9488, что согласуется с графиком.

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

В Excel есть три встроенные функции, позволяющие использовать третий метод определения значения наклона, точки пересечения по оси Y, коэффициента корреляции и R-квадрата набора данных. Функции являются SLOPE (), INTERCEPT (), CORREL () и RSQ (), а также рассматриваются в разделе статистики этого руководства.

Синтаксис для каждого следующий:

  • Уклон, м: = НАКЛОН (известные_y, известные_x)
  • Y-пересечение, b: = ПЕРЕСЕЧЕНИЕ (известные_y, известные_x)
  • Коэффициент корреляции, r: = КОРРЕЛЬ (известные_y, известные_x)
  • R-квадрат, r 2 : = RSQ (известные_y, известные_x)

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

Итак, повторюсь, мы можем определить наклон, точку пересечения по оси y и коэффициент корреляции любого набора данных с помощью трех методов Excel:

  1. Постройте данные и добавьте линию тренда
  2. Неявно использовать методы линейной регрессии
  3. Используйте встроенные функции Excel


Если у вас есть вопрос или комментарий, отправьте электронное письмо по адресу .


E = mc2 Calculator

Вы можете поспорить, что любой, кого вы встретите, хотя бы слышал об этом знаменитом уравнении. Вы, наверное, уже сто раз сталкивались с этим. Но что на самом деле означает E равняется mc в квадрате ? Что такое загадочный принцип эквивалентности массы и энергии? При чем здесь Эйнштейн? Продолжайте читать, чтобы узнать!

Теория относительности Эйнштейна

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

Формула эквивалентности энергии массы

Хорошо, вы, наверное, знаете, что:

E = mc²

Но что именно обозначают буквы в этом знаменитом уравнении?

  • м – это просто масса объекта.

  • c – это скорость света – постоянное значение 299 792 458 м / с.

  • E – энергия покоя объекта.

Последствия E = mc2

Если задуматься, последствия формулировки уравнения E = mc² были поразительными. Эйнштейн предложил мир, в котором масса – это просто энергия, ожидающая «освобождения». И не просто энергия – ее огромное количество.

Средний вес взрослого человека 62 кг. Такой человек, согласно Эйнштейну, имеет энергию покоя 5,6 × 10 12 мегаджоулей (МДж). Вы можете проверить это с помощью нашего калькулятора E mc2. Для сравнения: бомба, сброшенная на Нагасаки, имела энергию 8,4 × 10 7 МДж. По сути, если вам удастся взорвать и освободить всю свою энергию покоя (что еще не достижимо … пока), вы вызовете те же разрушения, что и более 66 000 ядерных бомб.Сказать, что это много, – определенно ничего не сказать.

Эта формула была спасением для всех, кто пытался понять, как устроена Вселенная. Это объяснило, почему радиоактивные материалы не «плавятся», испуская потоки частиц (это из-за чрезвычайной эффективности преобразования массы в энергию). Это также объяснило, почему у звезд не заканчивается водород, их основное топливо, в течение миллиардов лет. Внезапно возраст Вселенной оказался более достоверным.

Энергия из длины волны Пример задачи

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

Ключевые выводы: найти энергию фотона по длине волны

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

Энергия из проблемы длины волны – Энергия лазерного луча

Красный свет гелий-неонового лазера имеет длину волны 633 нм. Что такое энергия одного фотона?

Для решения этой проблемы вам нужно использовать два уравнения:

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

E = hν

где
E = энергия
h = постоянная Планка = 6,626 x 10 -34 Дж · с
ν = частота

Второе уравнение – это волновое уравнение, которое описывает скорость света через длину волны и частоту. Вы используете это уравнение, чтобы найти частоту для включения в первое уравнение. Волновое уравнение:
c = λν

где
c = скорость света = 3 x 10 8 м / сек
λ = длина волны
ν = частота

Перепишите уравнение для определения частоты:
ν = c / λ

Затем замените частоту в первом уравнении на c / λ, чтобы получить формулу, которую вы можете использовать:
E = hν
E = hc / λ

Другими словами, энергия фотографии прямо пропорциональна его частоте и обратно пропорциональна длине волны.

Осталось только подставить значения и получить ответ:
E = 6,626 x 10 -34 Дж · сек 3 x 10 8 м / сек / (633 нм x 10 -9 м / 1 нм )
E = 1,988 x 10 -25 Дж · м / 6,33 x 10 -7 м E = 3,14 x -19 Дж
Ответ:
Энергия одиночного фотона красного света от гелия. -неоновый лазер 3,14 x -19 Дж.

Энергия одного моля фотонов

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

Источник света излучает излучение с длиной волны 500,0 нм. Найдите энергию одного моля фотонов этого излучения. Выразите ответ в кДж.

Обычно требуется выполнить преобразование единиц измерения длины волны, чтобы заставить его работать в уравнении. Сначала преобразуйте нм в м. Nano- это 10 -9 , поэтому все, что вам нужно сделать, это переместить десятичный разряд на 9 точек или разделить на 10 9 .

500.0 нм = 500.0 x 10 -9 м = 5.000 x 10 -7 м

Последнее значение – это длина волны, выраженная в экспоненциальном представлении и правильном количестве значащих цифр.

Вспомните, как уравнение Планка и волновое уравнение были объединены, чтобы дать:

E = hc / λ

E = (6,626 x 10 -34 Дж · с) (3.000 x 10 8 м / с) / (5.000 x 10 -17 м)
E = 3,9756 x 10 -19 Дж

Однако это энергия одиночного фотона.Умножьте значение на число Авогадро для энергии моля фотонов:

энергия моля фотонов = (энергия одиночного фотона) x (число Авогадро)

энергия моля фотонов = (3,9756 x 10 -19 Дж) (6,022 x 10 23 моль -1 ) [подсказка: умножьте десятичные числа, а затем вычтите показатель знаменателя из показателя числителя, чтобы получить степень 10)

энергия = 2,394 x 10 5 Дж / моль

на один моль энергия равна 2.394 х 10 5 Дж

Обратите внимание, как значение сохраняет правильное количество значащих цифр. Для окончательного ответа его еще нужно преобразовать из Дж в кДж:

энергия = (2,394 x 10 5 Дж) (1 кДж / 1000 Дж)
энергия = 2,394 x 10 2 кДж или 239,4 кДж

Помните, что если вам нужно выполнить дополнительные преобразования единиц, следите за своими значащими цифрами.

Источники

  • Френч, А.П., Тейлор, Э.Ф. (1978). Введение в квантовую физику . Ван Ностранд Рейнхольд. Лондон. ISBN 0-442-30770-5.
  • Гриффитс, Д.Дж. (1995). Введение в квантовую механику . Прентис Холл. Верхняя река Сэдл, штат Нью-Джерси. ISBN 0-13-124405-1.
  • Landsberg, P.T. (1978). Термодинамика и статистическая механика . Издательство Оксфордского университета. Оксфорд, Великобритания. ISBN 0-19-851142-6.

Калькулятор потенциального магнитного поля для приложений физики Солнца с использованием шахматной сетки

A&A 625, A47 (2019)

Калькулятор потенциального магнитного поля для приложений физики Солнца с использованием шахматной сетки

Каллум М.Букок и Давид Циклаури

Школа физики и астрономии Лондонского университета королевы Марии, здание G.O. Jones, 327 Mile End Road, London E1 4NS, UK
электронная почта: [email protected]

Поступило: 19 Ноябрь 2018 г.
Принято: 18 маршировать 2019 г.

Аннотация

Программа была разработана для точной генерации потенциального магнитного поля на шахматной сетке путем экстраполяции магнитного поля, нормального к фотосферной поверхности.Программа сначала вычисляет магнитный потенциал, используя метод функции Грина, а затем использует схему конечных разностей для вычисления магнитного поля из потенциала. Была получена новая формула конечных разностей, учитывающая смещение сетки; показано, что эта формула дает численное приближение, наиболее близкое к реальному потенциальному полю. Также показано, что расширение области, в которой задается нормальное фотосферное поле, может повысить точность создаваемого потенциального поля.Программа представляет собой код FORTRAN 90, который можно использовать для генерации входных сигналов потенциального магнитного поля для Lare3d и других МГД-решателей, которые используют шахматную сетку для компонентов магнитного поля. Программа может быть распараллелена для быстрой работы на нескольких вычислительных ядрах. Код и вспомогательное описание приведены в приложениях.

Ключевые слова: Солнце: магнитные поля / методы: численная / магнитогидродинамика (МГД)

© ESO 2019

1. Введение

Проблема нагрева короны, вопрос о том, почему корона Солнца намного горячее (~ 1 МК), чем фотосфера (~ 6000 К), является постоянной проблемой в физике Солнца.Исследователи согласны с тем, что энергия переносится в корону за счет нетеплового переноса энергии через магнитное поле Солнца (Arregui 2015), хотя доминирующий механизм нагрева короны обсуждается (Parnell & De Moortel 2012). Вычислительные трехмерные магнитогидродинамические (МГД) модели часто используются для решения проблемы нагрева короны (Климчук, 2015).

Вычислительные модели МГД требуют задания начального магнитного поля. Для многих проблем, в первую очередь тех, которые касаются МГД-волн, полезно иметь статическое равновесие, на которое можно воздействовать возмущениями скорости и магнитных полей (Goossens 2003).В случае статического равновесия вектор скорости v и его производная по времени ∂ v / ∂ t должны равняться нулю.

Аналитические выражения существуют для ряда структур, которые описывают магнитостатические равновесия как в 2D, так и в 3D, например (Смит и др., 2007; Петрухин и др., 2018; Оливер и др., 1998; Куперман и др., 1989; Гент и др., 2013). ). Для более детальных магнитных структур нам нужна информация из наблюдений за Солнцем. Карты прямой видимости и векторные магнитограммы от фотосферы могут быть измерены с помощью спектрополяриметрических методов, таких как эффект Зеемана, т.е.е. расщепление спектральных линий в присутствии магнитного поля (Beckers 1971). Однако гораздо труднее напрямую измерить магнитное поле в солнечной короне (Руан и др., 2008). К счастью, фотосферное магнитное поле можно использовать для восстановления коронального поля путем экстраполяции. Экстраполяция фотосферного магнитного поля в настоящее время является основным инструментом моделирования корональных магнитных полей (Tadesse et al. 2014).

Существует множество методов экстраполяции структуры магнитного поля из измерений поверхности.Обычно используемые методы основываются на различных предположениях (Neukirch 2005). Чаще всего пренебрегают немагнитными силами, такими как градиенты давления и гравитация; это хорошо оправдано в солнечной короне из-за низкого бета-излучения плазмы (Wiegelmann, 2008) и с учетом масштабов, меньших, чем гидростатический масштаб высоты (Peter et al. 2015). Для магнитостатического равновесия магнитная сила Лоренца должна быть равна нулю Дж × B = 0. Это определяет бессиловое поле (Wiegelmann & Sakurai 2012).Другой способ выразить это условие дается ниже в формуле. (1) где α может быть либо нулем, либо константой, либо переменной, постоянной вдоль силовых линий. Эти случаи соответствуют потенциальному, линейному и нелинейному бессиловым полям соответственно (Aschwanden 2004). Это выражение записывается как

(1)

Различные методы экстраполяции используются для восстановления потенциальных полей, линейных бессиловых полей (Gary 1989) и нелинейных бессиловых полей (Wiegelmann 2008).Существуют также методы экстраполяции для восстановления несиловых полей, хотя они могут потребовать дополнительных данных (Wiegelmann 2004). Более того, некоторые из этих методов могут быть расширены от декартовых до сферических координат, так что солнечное корональное поле можно рассчитать глобально (Wiegelmann 2007; Wiegelmann et al. 2007).

В частном случае потенциального бессилового поля существует нулевой ток Дж = × B = 0. Поле называется потенциальным магнитным полем, потому что магнитное поле может быть выражено через скалярный потенциал B = ∇ ϕ .Потенциальные магнитные поля полезны, потому что их можно использовать, например, в качестве начальной топологии для МГД-моделирования (Gudiksen & Nordlund 2005; Masson et al. 2009; Bingert & Peter 2011, 2013; Bourdin et al. 2013) для изучения МГД-волн. явлений (Thackray & Jain 2017; Ofman & Selwa 2008; Smith et al. 2007; Petrukhin et al. 2018), чтобы внедрить непотенциальные поля (Browning et al. 2008) или создать несиловые МГД-равновесия ( Гордовский и др. 2014; Соланки и Штайнер 1990; Хоменко и Колладос 2006; Пиццо 1986; Вигельманн и Нойкирх 2006; Иноуэ 2016; Иноуэ и др.2013).

Чтобы вычислить потенциальное магнитное поле из нормального магнитного поля B n в фотосфере, можно использовать множество методов, таких как разложение Фурье, разложения по сферическим гармоникам и функция Грина. Сравнение метода Фурье и метода функций Грина проведено Сакураи (1982). Экстраполированные поля требуют, чтобы граничные условия для экстраполяции соответствовали условиям МГД-моделирования (Отто и др., 2007).

В этой работе мы рассматриваем метод функции Грина, который использует дискретную аппроксимацию для B n в пределах каждой сетки.Этот метод лучше подходит для рассмотрения непериодических изолированных магнитных областей, но имеет недостаток, заключающийся в предположении, что B n = 0 за пределами интересующей области. Этот метод был впервые использован Шмидтом (1964), позже был разработан Семелом (1967) для наклонных фотосферных данных, а затем адаптирован Сакураи, который скорректировал метод для практических приложений на конечной числовой сетке (Сакураи, 1982).

В этой статье описан калькулятор потенциального поля, основанный на модифицированном методе функции Грина.Целью данной статьи является демонстрация улучшенного метода конечных разностей, подходящего для решателя МГД, использующего шахматное магнитное поле. Это также демонстрирует, что за счет увеличения объема входных данных относительно вычислительной области можно уменьшить ошибки, связанные с предположением B n = 0 за пределами интересующей области.

Объясняются различные методы численного дифференцирования, которые рассматривались для вычислителя потенциального поля, и сравниваются с наиболее точными методами, выбранными для окончательной реализации.Калькулятор был разработан для создания начальных входных значений магнитного поля для Lare3d (Arber et al. 2001), но работает с любым решателем MHD, который использует шахматные сетки для компонентов магнитного поля.

2. Метод функции Грина

Принимая Дж = × B = 0 в качестве начального условия, закон магнетизма Гаусса может быть выражен как уравнение Лапласа для магнитного потенциала (Ашванден, 2004).

(2)

Решая уравнение Лапласа для скалярного потенциала ϕ , можно вычислить потенциальное магнитное поле B , которое приводит к статическому равновесию.Рассмотрим трехмерную область с вертикальной координатой z и границей фотосфер на z = 0. Обозначим магнитное поле, нормальное к границе фотосферы, как B z ( x , y , z ), а нормальное поле на границе фотосферы тогда будет B z 0 ( x , y ) = B z ( x , y , 0 ).

Уравнение Лапласа с фотосферным граничным условием B z 0 = ∂ ϕ / ∂ z может быть решено численно на дискретной сетке точек ( i , j , k ) с разделением сетки d , используя модифицированный метод функции Грина, приведенный ниже в формуле.(3) (Sakurai 1982), т.е.

(3)

, где i 1 и j 1 – фиктивные переменные, используемые для суммирования вкладов от каждой точки фотосферы. Нормальное магнитное поле вносит вклад в потенциал посредством дискретной аппроксимации функции Грина. В силу выбранной функции Грина потенциал удовлетворяет уравнению Лапласа и граничным условиям при z = 0 и r → ∞ следующим образом:

(4)

Этот конкретный выбор функции Грина представляет «подводный магнитный источник» на глубине под каждой точкой сетки на z = 0.Этот выбор функции Грина вызывает небольшое неудобство, поскольку вычисленное нормальное магнитное поле не совсем соответствует исходному нормальному полю при z = 0. Хотя есть и другие возможные варианты, которые дают более точное совпадение, один из которых описан. в Sakurai (1982) эта конкретная функция Грина была выбрана для удобства вычислений, потому что в этой статье основное внимание уделяется не выбору функции Грина, а, скорее, методу конечных разностей, используемому для дифференцирования магнитного потенциала.

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

(5)

Если для дифференцирования магнитного потенциала используется схема конечных разностей, то потенциал необходимо рассчитывать в дополнительных «фантомных» ячейках на границах домена. Количество дополнительных ячеек зависит от порядка разностной схемы. Призрачные ячейки также требуются во входных данных для Lare3d, как описано в Разд.3.

Поскольку потенциал должен быть рассчитан в ячейках-призраках, потенциал магнитного источника, представленный функцией Грина, также должен быть расположен под этими ячейками и вне области. Это может быть сделано либо путем увеличения индексов в направлении z так, чтобы фантомные ячейки поднимались над источником магнитного поля, либо путем изменения уравнения. 3) эффективно опустить магнитный источник на необходимое количество ячеек.

В любом случае общий эффект состоит в том, чтобы рассматривать самую низкую фантомную ячейку как новую позицию исходного нормального поля.Магнитное поле, рассматриваемое в МГД-решателе, затем смещается вверх от фотосферы на длину дополнительных фантомных ячеек. Хотя эта неточность в положении неудобна, и есть способы избежать этого, например, вручную исправить граничные условия в Lare3d или не дифференцировать для определения магнитного поля в двух нижних слоях, мы решили внести эту модификацию на основании того, что любая ручная фиксация граничных условий влияет на равновесие потенциального поля способами, которые не зависят от исследуемых методов конечной разности.Таким образом, модифицированный метод функции Грина, используемый в этой статье, представляет собой небольшую модификацию уравнения (3), который погружает магнитный источник, представленный функцией Грина, на дополнительные n ячеек, где n – максимальное необходимое количество фантомных ячеек, т. Е.

(6)

3. Сетка в шахматном порядке

Расчетная область в Lare3d имеет размер ( nx , ny , nz ), однако магнитное поле определяется на шахматной сетке, так что компоненты магнитного поля расположены в центрах граней каждой ячейки.Например, значения B x смещены от центра ячейки на половину ячейки в направлении x , как показано на рисунке 1, и аналогично для B y и B z в направлениях y и z соответственно. Таким образом, компоненты магнитного поля располагаются в шахматном порядке, чтобы сохранить соленоидальное свойство магнитного поля ∇ ⋅ B = 0, поскольку оно определяется уравнением индукции (Balsara & Kim 2004).Смещение сетки изменяет положение, в котором компоненты магнитного поля должны быть указаны по отношению к магнитному потенциалу. Из-за шахматной сетки размеры B x , B y и B z в вычислительной области записываются как

(7)

Рисунок 1.

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

Открыть с помощью DEXTER

Числовая сетка, на которой должно быть задано начальное магнитное поле в Lare3d, однако, больше. Это связано с тем, что Lare3d имеет второй порядок точности в пространстве и, следовательно, требует двух фантомных ячеек на каждой границе. Значения B x , B y и B z , которые должны быть указаны в начальном магнитном поле, поэтому записываются как

(8)

4.Методика сравнения

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

4.1. Тестовые примеры

Для этого исследования мы использовали два аналитически определенных тестовых случая.Первый тестовый случай представляет собой область, содержащую бипольную петлю, а второй – область, содержащую униполярное магнитное поле. Для обоих тестовых случаев магнитные поля рассчитываются для области размером nx = ny = nz = 200 с шагом сетки d = 0,01. Определяя расширение нашей области как [−1: 1, −1: 1,0: 2] с его началом в (0,0,0), тестовые случаи могут быть определены в терминах входного магнитного поля B z 0 следующим образом:

Контрольный пример 1 .Формула для бипольной петли определяется путем взятия аналитического выражения для бессилового бипольного поля, приведенного в Cuperman et al. (1989). Мы берем текущую бесплатную версию этих уравнений и устанавливаем координаты магнитных полюсов равными ( x 01 , y 01 , z 01 ) = (−0,5, 0, −1) и ( x 02 , y 02 , z 02 ) = (0,5, 0, −1). Наконец, мы берем компонент B z и устанавливаем z = 0, чтобы получить

(9)

Контрольный пример 2 .Формула для униполярного магнитного поля определяется исходя из предположения, что сила нормального магнитного поля B z спадает как гауссово по мере удаления от источника магнитного поля. Мы принимаем положение источника магнитного источника в начале координат на единицу расстояния ниже фотосферы, чтобы получить

.

(10)

Потенциальные поля, создаваемые этими тестовыми примерами, представлены диаграммами силовых линий на рис. 2 и 3; затенение у основания доменов на этих изображениях представляет силу и направление магнитного поля в фотосфере.

Рис. 2.

Линейная диаграмма магнитного поля, используемого в тестовом примере 1. Фотосферная поверхность в основании домена представляет собой цветной контур B z 0 .

Открыть с помощью DEXTER
Рис 3.

Линейная диаграмма магнитного поля, используемого в тестовом примере 2. Фотосферная поверхность в основании домена представляет собой цветной контур B z 0 .

Открыть с помощью DEXTER

4.2. Метрики

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

Затем мы установили начальную скорость на ноль, начальное давление на постоянное значение P = 5 × 10 −9 и удельное сопротивление на η = 5 × 10 −5 . Мы установили γ = 5/3, чтобы жидкость была адиабатической, и установили ударную вязкость ν 1 = 0.1 и ν 2 = 0,5. Мы устанавливаем граничные условия так, чтобы скорость и магнитное поле имели фиксированные значения на границах.

Затем мы позволяем моделированию развиваться для 100 τ A , где τ A – время Альвена. Количество требуемых временных шагов варьируется, поскольку Lare3d использует ограниченный временной интервал CFL на каждой итерации. Это не является фазой релаксации, а скорее проверкой того, насколько стабильна конфигурация потенциального поля с течением времени.Наконец, мы измеряем ток, скорость и разницу между конечным и начальным магнитными полями в домене.

Метрики, используемые для анализа эффективности каждого метода, приведены ниже. Чем ближе эти значения к нулю, тем ближе поле к потенциальному полю и тем ближе наше равновесие к статическому равновесию. Средние плотности тока – это среднее значение плотности тока Дж, по всем точкам сетки. Эти показатели включают

MAJ 0 : максимальная абсолютная плотность тока в нулевой момент времени.

J avg0 : Средняя абсолютная плотность тока в нулевой момент времени.

MAJ f : максимальная абсолютная плотность тока после 100 τ A .

J avg f : Средняя абсолютная плотность тока после 100 τ A .

MAV: максимальная абсолютная скорость после 100 τ A .

Δ B max : максимальная разница между начальным и конечным магнитными полями.

5. Дифференциация магнитного потенциала

ϕ

5.1. Метод

Чтобы вычислить компоненты магнитного поля в нашей области, нам нужно дифференцировать потенциал магнитного поля в соответствии с формулой. (5). Эффект смещения сетки заключается в том, что составляющая магнитного поля в определенном направлении должна вычисляться на краях ячейки в этом направлении, а не в центре ячейки. Однако магнитный потенциал ϕ рассчитывается в центрах ячеек.Это показано на рис. 1.

В этом разделе представлены три метода численного дифференцирования, каждый из которых по-своему учитывает это расхождение. Расчет B x ( i , j , k ) взят в качестве примера. Индексы j и k считаются фиксированными, а B x ( i , j , k ) обозначается просто как.

Метод А . В этом методе используется формула центральной разности порядка 𝒪 ( d 4 ) без учета смещения сетки. Магнитное поле, рассчитанное по методу А, обозначено как. Формула, используемая для расчета, приведена в формуле. (11) ниже. Этот метод требует дополнительных четырех фантомных ячеек в каждом направлении, т. Е.

(11)

Метод B . В этом методе смещение сетки учитывается путем вычисления компонентов магнитного поля в центрах ячеек с использованием метода A и последующего усреднения этих значений для получения компонентов магнитного поля на краях ячейки.Магнитное поле, рассчитанное по методу B , обозначено как. Формула, используемая для расчета, приведена в формуле. (12) ниже. Этот метод требует дополнительных пяти ячеек-призраков в каждом направлении, т. Е.

(12)

Метод C . В этом методе смещение сетки учитывается путем расчета составляющих магнитного поля с использованием другой формулы центральной разности. Полный вывод этой формулы приведен в Приложении A. Магнитное поле, рассчитанное с использованием метода C, обозначено как.Формула, используемая для расчета, приведена в формуле. (13) ниже. Этот метод требует дополнительных трех призрачных ячеек в каждом направлении, то есть

(13)

5.2. Результатов

Для каждого метода дифференцирования мы сгенерировали потенциальное поле для каждого тестового примера. Потенциальные поля были сгенерированы из входных полей B z 0 , определенных на основе нашей расчетной области. Мы проанализировали и сравнили потенциальные генерируемые поля, и результаты нашего анализа представлены в таблицах 1 и 2.

Таблица 1.

Значение каждой метрики для метода A, метода B и метода C численного дифференцирования для контрольного примера 1.

Таблица 2.

Значение каждой метрики для метода A, метода B и метода C численного дифференцирования для контрольного примера 2.

Для Контрольного примера 1 потенциальное поле, созданное с использованием метода 1, не было доведено до 100 τ A . Токи, присутствующие в начальном магнитном поле, вызвали увеличение v , что значительно уменьшило ρ в некоторых частях домена.Это привело к значительному сокращению ограниченного временного шага CFL по мере продвижения моделирования, что привело к невероятно длительному времени выполнения моделирования, равному 100 τ A ; вместо этого приведены значения для 50 τ A .

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

6. Размер области ввода для

B z0

6.1. Метод

Теперь рассмотрим размер области ввода, то есть области на границе фотосфер, над которой задано поле ввода B z 0 .Размер вычислительной области для ϕ обозначен как ( Nx , Ny , Nz ), а размер входной области для B z 0 обозначен как ( Nx 1 , Нью-Йорк 2 ). До сих пор поле ввода B z 0 было определено только на основе нашей вычислительной области. Другими словами ( Nx 1 , Ny 2 ) было установлено точно равным ( Nx , Ny ).Однако поле ввода B z 0 может быть определено для большей области.

Мотивом для расширения нашей области ввода является то, что текущее значение J в наших предыдущих потенциальных полях было сконцентрировано на границах x / y . Причину этого можно увидеть на рис. 6; на этих границах резко меняется градиент напряженности магнитного поля. Расширяя область ввода, мы пытаемся уменьшить влияние этого изменения градиента и уменьшить ток на границах.Рассматриваем два варианта размера входной области:

Стандартная область ввода – ( Nx 1 , Ny 2 ) установлена ​​равной точно ( Nx , Ny ), а центры области ввода и вычислительной области выровнены. Затем входные данные покрывают точно такую ​​же область фотосферы, что и расчетная область для ϕ . Это показано на рис. 4.

Инжир.4.

Иллюстрация расчетной области для ϕ , показанной серым прямоугольником, и области ввода для B z 0 , показанной синим квадратом, когда ( Nx 1 , Ny 2 ) равно ( Nx , Ny ).

Открыть с помощью DEXTER

Расширенная область ввода – ( Nx 1 , Ny 2 ) установлена ​​равной (3 Nx , 3 Ny ), а центры области ввода и вычислительной области выровнены.Тогда входные данные покрывают гораздо большую область фотосферы, чем расчетная область для ϕ . Это показано на рис. 5.

Рис. 5.

Иллюстрация расчетной области для ϕ , показанной серым прямоугольником, и области ввода для B z 0 , показанной синим квадратом, когда ( Nx 1 , Ny 2 ) равно (3 Nx , 3 Ny ).

Открыть с помощью DEXTER

6.2. Результатов

Для каждого тестового примера мы создаем потенциальные поля, используя как стандартные, так и расширенные области ввода. Мы использовали формулу. (6) для расчета нашего магнитного потенциала ϕ и Метод C для выполнения численного дифференцирования. Каждое сгенерированное потенциальное поле прогонялось в Lare3d для 100 τ A перед анализом и сравнением. Результаты нашего анализа представлены в таблицах 3 и 4.Мы можем видеть, что магнитное поле, создаваемое с использованием большей области ввода, ближе к потенциальному полю. Кроме того, как видно из рис. 6 и 7, что использование большей области ввода создает более гладкие магнитные поля на границах.

Таблица 3.

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

Таблица 4.

Значение каждой метрики с использованием стандартных и расширенных областей ввода для тестового примера 2.

Инжир.6.

Закрашенная поверхностная диаграмма напряженности магнитного поля при z = 0 для потенциального поля, созданного с использованием стандартной области ввода. Расширяющиеся края возникают из-за неоднородности на краю входной области.

Открыть с помощью DEXTER
Рис. 7.

Закрашенная поверхностная диаграмма напряженности магнитного поля при z = 0 для потенциального поля, созданного с использованием расширенной области ввода.Края более гладкие из-за расширенной области ввода.

Открыть с помощью DEXTER

Стоит отметить, что из-за характера используемого метода функции Грина вычислитель потенциала ожидает, что значение B z 0 будет равно нулю за пределами области ввода. Если за пределами области ввода ожидаются значительные магнитные поля, то точность вычислителя можно повысить, увеличив размер области ввода; однако существует компромисс между точностью и временем вычисления, помимо которого могут быть доступны ограниченные фотосферные данные.

7. Сравнение с аналитической бипольной петлей

7.1. Метод

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

Используя выбранный нами метод дифференцирования и расширенную область ввода, мы применяем наш метод Грина для экстраполяции, как подробно описано в Разд.3, чтобы рассчитать потенциальное поле. Это поле создается с использованием нормального фотосферного поля для бипольной петли, указанного как Контрольный пример 1 в Разд. 4. Мы сравниваем это полевое решение с аналитическим решением для потенциального бипольного поля, данным Куперманом и др. (1989), который имеет такое же нормальное магнитное поле в фотосфере. Это дает

(14)

, где R j – это расстояния от каждого магнитного полюса для j = 1, 2, расположенных в ( x 01 , y 01 , z 01 ) = (-0.5, 0, −1) и ( x 02 , y 02 , z 02 ) = (0,5, 0, −1).

7.2. Результатов

Магнитное поле, создаваемое нашим вычислителем потенциала, и потенциальное поле, определяемое уравнением. (14) были запущены в Lare3d и проанализированы так же, как и для предыдущих магнитных полей. Это было сделано для сравнения как начальных магнитных полей, так и эволюции этих полей после 100 τ A времени моделирования в Lare3d.

Результаты этого анализа приведены в таблице 5. Токи, присутствующие в аналитическом поле, можно объяснить численной диффузией. Числовая диффузия вызвана центральной разностной схемой, используемой для расчета токов в Lare3d. Ошибка усечения или центральная разностная схема, используемая для расчета токов, имеет порядок 𝒪 ( d 2 ), где h – шаг сетки. Действительно, мы можем видеть, что максимальный ток в этом решении также имеет порядок 𝒪 ( d 2 ) с шагом сетки d = 0.01.

Таблица 5.

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

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

8.Данные солнечной магнитограммы

8.1. Метод

Определив предпочтительную методологию для нашего потенциального калькулятора и сравнив ее с аналитическим решением, мы завершаем наш анализ рассмотрением магнитного поля, созданного на основе данных солнечной магнитограммы. Используемые данные были взяты из данных магнитограмм прямой видимости, собранных прибором гелиосейсмической и магнитной визуализации (HMI) на обсерватории солнечной динамики (SDO), и доступ к ним был получен через онлайн-базу данных Объединенного научного операционного центра (JSOC).Данные были получены из активной области вблизи центра солнечного диска между 14:02 и 14:56 31 августа 2003 года. Поскольку данные были взяты из небольшой области в центральном диске, они были приняты для приблизительного представления данных нормального магнитного поля в этой области. Данные были нормализованы так, чтобы максимальная напряженность поля была равна единице, и поэтому результаты сопоставимы с результатами из предыдущих разделов.

Контурный график данных магнитограммы приведен на рис. 8, а потенциальное поле, созданное с использованием этих данных магнитограммы, представлено диаграммой силовых линий, показанной на рис.9. Используя эти данные, потенциальные поля были созданы с использованием каждого из трех методов численного дифференцирования, описанных в разд. 5.1 для стандартной области ввода, как описано в Разд. 6.1, и используя наш предпочтительный метод C для численного дифференцирования для расширенной области ввода.

Рис. 8.

Контурная диаграмма исходных данных магнитограммы нормального магнитного поля на фотосфере. Цвет указывает направление и силу магнитного поля.

Открыть с помощью DEXTER
Рис.9.

Схема силовых линий потенциального магнитного поля, созданного с использованием данных солнечной магнитограммы. Фотосферная поверхность в основании домена представляет собой цветной контур B z 0 .

Открыть с помощью DEXTER

8.2. Результатов

Создаваемые потенциальные поля обрабатывались в Lare3d и анализировались так же, как и для предыдущих магнитных полей.Потенциальное поле, созданное с использованием метода А для численного дифференцирования, который игнорирует влияние шахматной сетки, не будет работать даже для 5 τ A и будет иметь начальные максимальные токи на порядок больше, чем любой другой потенциал. поля. Результаты анализа остальных потенциальных полей приведены в таблице 6.

Таблица 6.

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

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

Сравнивая теперь результаты для использования стандартной или расширенной области ввода (оба используют метод C для дифференциации), мы видим удивительный результат.Максимальный ток одинаков в обоих случаях, а средний ток немного ниже при использовании стандартной области ввода. Хотя максимальная скорость при 100 τ A ниже, когда используется расширенная область ввода, эта разница меньше, чем на порядок величины. Причина в том, что использование расширенной области ввода имеет небольшой эффект, что можно ясно увидеть, сравнив напряженность магнитного поля для потенциального поля, созданного с помощью стандартной области ввода и расширенной области ввода.Напряженность поля в фотосфере показана для каждого из этих полей на рис. 10 и 11. Хорошо видно, что в отличие от рис. 7 в разд. 6.1 напряженность поля на рис. 11 намного меньше и более плоская на границах, поэтому предположение, что B z 0 = 0 вне области имеет гораздо меньший эффект.

Рис 10.

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

Открыть с помощью DEXTER
Рис 11.

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

Открыть с помощью DEXTER

Хотя расширение области ввода в целом повышает точность решения потенциального поля, эти результаты показывают, что преимущества увеличения размера области ввода иногда являются лишь предельными, в зависимости от поля ввода B z 0 .Следует помнить, что для достижения небольшого снижения максимальной скорости при 100 τ A , область ввода пришлось расширить примерно в девять раз. Это означает, что требуется в девять раз больше исходных данных, а также означает, что время выполнения увеличится в девять раз.

9. Заключение

Вычислитель потенциального поля был разработан для точного расчета потенциального магнитного поля в области с нормальным полем B z 0 на границе фотосферы.Дизайн кода был направлен на повышение точности за счет учета: расчета магнитного потенциала ϕ , дифференцирования магнитного потенциала ϕ и размера входной области для B z 0 as описано в разд. 2, 5.1 и 6.1 этой статьи соответственно.

Калькулятор потенциального поля использует формулу. (6) для расчета магнитного потенциала. Затем потенциал дифференцируется с использованием новой формулы центральной разности Eq.(13), полученное в Приложении А, для расчета компонент магнитного поля. Было показано, что игнорирование разнесенных магнитных полей во время дифференцирования магнитного потенциала приводит к нестабильному потенциальному полю и что, хотя простая интерполяция может быть эффективной, выбранный метод численного дифференцирования в конечном итоге более точен, быстрее вычисляется и требует меньшего количества паразитов. клетки для выполнения. Далее это было показано в Разд. 7.1 видно, что для аналитически заданных входных данных этот метод может создавать потенциальные поля, которые являются правильными с точностью до числового распространения сетки моделирования.

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

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

Хотя методы функций Грина, используемые в этой работе, были доступны в течение значительного периода времени (Schmidt 1964; Semel 1967; Sakurai 1982), есть надежда, что эта точная и простая реализация будет весьма полезна для тех, кто работает с потенциальными полями. на шахматной сетке. Код потенциального калькулятора написан на FORTRAN 90 и был разработан для создания начальных вводов магнитного поля для Lare3d, но будет работать с любым решателем MHD, который использует шахматную сетку для компонентов магнитного поля.Программа может быть распараллелена для быстрой работы на нескольких вычислительных ядрах. Этот код и сопроводительная документация предоставляются вместе с этим документом.

Благодарности

Каллум Букок благодарит UK STFC DISCnet за финансовую поддержку его докторантуры. Исследование Давида Циклаури при поддержке Национального научного фонда Грузии Шота Руставели в рамках конкурса совместных исследований (звонок соотечественников) 2016-2019. В этом исследовании использовалась установка Queen Mary’s Apocrita HPC при поддержке QMUL Research-IT (King et al.2017).

Список литературы

  1. Арбер Т., Лонгботтом А., Джеррард К. и Милн А. 2001, J. Comput. Физ., 171, 151 [НАСА ОБЪЯВЛЕНИЕ] [CrossRef] [Google Scholar]
  2. Арреги, И.2015, Филос. Пер. R. Soc. Лондон сер. А, 373, 20140261 [НАСА ОБЪЯВЛЕНИЕ] [CrossRef] [Google Scholar]
  3. Ашванден, М.J. 2004, Физика солнечной короны. Введение (Praxis Publishing Ltd) [Google Scholar]
  4. Балсара, Д. С., и Ким, Дж. 2004, ApJ, 602, 1079 [НАСА ОБЪЯВЛЕНИЕ] [CrossRef] [Google Scholar]
  5. Бекерс, Дж.М. 1971 г., в сб .: Solar Magnetic Fields, ed. Р. Ховард, IAU Symp., 43, 3 [НАСА ОБЪЯВЛЕНИЕ] [CrossRef] [Google Scholar]
  6. Бингер, С., И Питер, Г. 2011, A&A, 530, A112 [НАСА ОБЪЯВЛЕНИЕ] [CrossRef] [EDP Sciences] [Google Scholar]
  7. Бингер, С., И Питер, H. 2013, A&A, 550, A30 [НАСА ОБЪЯВЛЕНИЕ] [CrossRef] [EDP Sciences] [Google Scholar]
  8. Бурден, П.А., Бингерт, С., и Питер, Х., 2013 г., A&A, 555, A123 [НАСА ОБЪЯВЛЕНИЕ] [CrossRef] [EDP Sciences] [Google Scholar]
  9. Браунинг, П.К., Джеррард К., Худ А. В., Кевис Р. и ван дер Линден Р. А. М. 2008, A&A, 485, 837 [НАСА ОБЪЯВЛЕНИЕ] [CrossRef] [EDP Sciences] [Google Scholar]
  10. Куперман, С., Офман Л. и Семел М. 1989, A&A, 216, 265. [НАСА ОБЪЯВЛЕНИЕ] [Google Scholar]
  11. Гэри, Г.А. 1989, ApJS, 69, 323 [НАСА ОБЪЯВЛЕНИЕ] [CrossRef] [EDP Sciences] [Google Scholar]
  12. Гент, Ф.А., Федун, В., Мамфорд, С. Дж., И Эрдейи, Р. 2013, MNRAS, 435, 689 [НАСА ОБЪЯВЛЕНИЕ] [CrossRef] [Google Scholar]
  13. Гуссенс, М.2003, Введение в астрофизику плазмы и магнитогидродинамику (Springer), 294 [CrossRef] [Google Scholar]
  14. Гордовский, М., Браунинг, П. К., Контар, Э. П., и Биан, Н. Х. 2014, A&A, 561, A72 [НАСА ОБЪЯВЛЕНИЕ] [CrossRef] [EDP Sciences] [Google Scholar]
  15. Гудиксен, Б.В., & Нордлунд, А. 2005, ApJ, 618, 1020 [НАСА ОБЪЯВЛЕНИЕ] [CrossRef] [Google Scholar]
  16. Иноуэ, С.2016, Прогресс в науке о Земле и планетах, 3, 19 [НАСА ОБЪЯВЛЕНИЕ] [CrossRef] [Google Scholar]
  17. Иноуэ, С., Магара, Т., и Пандей, В. 2013, ApJ, 780, 101 [НАСА ОБЪЯВЛЕНИЕ] [CrossRef] [Google Scholar]
  18. Хоменко, Э., & Колладос, М. 2006, ApJ, 653, 739 [НАСА ОБЪЯВЛЕНИЕ] [CrossRef] [Google Scholar]
  19. Кинг, Т., Бутчер, С., и Залевски, Л. 2017, Apocrita– кластер высокопроизводительных вычислений для Лондонского университета королевы Марии, DOI 10.5281 / zenodo.438045 [Google Scholar]
  20. Климчук Ю.А. 2015, Философ. Пер. R. Soc. Лондон сер. А, 373, 20140256 [НАСА ОБЪЯВЛЕНИЕ] [CrossRef] [Google Scholar]
  21. Массон, С., Pariat, E., Aulanier, G., & Schrijver, C.J. 2009, ApJ, 700, 559 [НАСА ОБЪЯВЛЕНИЕ] [CrossRef] [Google Scholar]
  22. Нойкирх, Т.2005, в Хромосферных и корональных магнитных полях, ред. Д. Э. Иннес, А. Лагг, С. А. Соланки, ESA Spec. Опубл., 596, 12.1 [Google Scholar]
  23. Офман Л. и Селва М. 2008, Proc. Int. Astron. Союзная, 4, 151 [CrossRef] [Google Scholar]
  24. Оливер, Р., Муравски К. и Баллестер Дж. Л. 1998, A&A, 330, 726 [НАСА ОБЪЯВЛЕНИЕ] [Google Scholar]
  25. Отто, А., Бюхнер, Дж. И Никутовски, Б. 2007, A&A, 468, 313 [НАСА ОБЪЯВЛЕНИЕ] [CrossRef] [EDP Sciences] [Google Scholar]
  26. Парнелл, К.E., & De Moortel, I. 2012, Philos. Пер. R. Soc. Лондон сер. А, 370, 3217 [НАСА ОБЪЯВЛЕНИЕ] [CrossRef] [Google Scholar]
  27. Питер, Х., Варнеке, Дж., Читта, Л. П., и Кэмерон, Р. Х. 2015, A&A, 584, A68 [НАСА ОБЪЯВЛЕНИЕ] [CrossRef] [EDP Sciences] [Google Scholar]
  28. Петрухин, Н., Рудерман, М., & Диденкулова, Е. 2018, МНРАС, 474, 2289 [НАСА ОБЪЯВЛЕНИЕ] [CrossRef] [Google Scholar]
  29. Пиццо, В.J. 1986, ApJ, 302, 785 [НАСА ОБЪЯВЛЕНИЕ] [CrossRef] [Google Scholar]
  30. Руан, П., Wiegelmann, T., Inhester, B., et al. 2008, A&A, 481, 827 [НАСА ОБЪЯВЛЕНИЕ] [CrossRef] [EDP Sciences] [Google Scholar]
  31. Сакураи, Т.1982, Sol. Физ., 76, 301 [НАСА ОБЪЯВЛЕНИЕ] [CrossRef] [Google Scholar]
  32. Шмидт, Х.U.1964, НАСА Spec. Опубл., 50, 107 [НАСА ОБЪЯВЛЕНИЕ] [Google Scholar]
  33. Семель, М.1967, Ann. Астрофиз., 30, 513 [НАСА ОБЪЯВЛЕНИЕ] [Google Scholar]
  34. Смит, П.Д., Циклаури, Д., и Рудерман, М. С. 2007, A&A, 475, 1111 [НАСА ОБЪЯВЛЕНИЕ] [CrossRef] [EDP Sciences] [Google Scholar]
  35. Соланки, С.К. и Штайнер, О. 1990, A&A, 234, 519 [НАСА ОБЪЯВЛЕНИЕ] [Google Scholar]
  36. Тадессе, Т., Вигельманн, Т., Госаин, С., Макнейс, П., и Певцов, А.А. 2014, A&A, 562, A105 [НАСА ОБЪЯВЛЕНИЕ] [CrossRef] [EDP Sciences] [Google Scholar]
  37. Текрей, Х., & Джайн, р. 2017, A&A, 608, A108 [НАСА ОБЪЯВЛЕНИЕ] [CrossRef] [EDP Sciences] [Google Scholar]
  38. Вигельманн, Т.2004, Sol. Физ., 219, 87 [НАСА ОБЪЯВЛЕНИЕ] [CrossRef] [Google Scholar]
  39. Вигельманн, Т.2007, Sol. Физ., 240, 227 [НАСА ОБЪЯВЛЕНИЕ] [CrossRef] [Google Scholar]
  40. Вигельманн, Т.2008, J. Geophys. Res .: Space Phys., 113, A3. [CrossRef] [Google Scholar]
  41. Вигельманн, Т., & Neukirch, T. 2006, A&A, 457, 1053 [НАСА ОБЪЯВЛЕНИЕ] [CrossRef] [EDP Sciences] [Google Scholar]
  42. Вигельманн, Т., & Сакураи, Т. 2012, Лив. Преподобный Sol. Физ., 9, 5 [Google Scholar]
  43. Wiegelmann, T., Neukirch, T., Ruan, P., & Inhester, B. 2007, A&A, 475, 701 [НАСА ОБЪЯВЛЕНИЕ] [CrossRef] [EDP Sciences] [Google Scholar]

Приложение A: Вывод формулы численного дифференцирования

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

В качестве примера возьмем расчет B x ( i , j , k ), индексы j и k приняты фиксированными, а B x ( i , j , k ) обозначается просто как.Теперь рассмотрим положение относительно точек, где ϕ определено вдоль оси x . Если мы определим d как шаг сетки и h как половину шага сетки, то для определенного в точке x 0 , ϕ определяется на расстояниях h и 3 h в любую сторону из x 0 . Это показано на рис. A.1. Теперь мы вычисляем производную ϕ относительно x в точке x 0 , то есть ϕ ′ ( x 0 ).Начнем с разложения Тейлора ϕ i −1 , ϕ i , ϕ i + 1 и ϕ i + 2 около x они приведены ниже в уравнениях. (A.1) – (A.4) следующим образом:

(A.1)

(A.2)

(A.3)

(A.4)

Рис. А.1.

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

Открыть с помощью DEXTER

Сначала устраним четные степени h , взяв ϕ i + 2 ϕ i −1 и ϕ i + 1 ϕ i , как показано ниже в уравнениях. (A.5) и (A.6):

(A.5)

(A.6)

Затем нам нужно удалить член третьего порядка, взяв 27 ( ϕ i + 1 ϕ i ) – ( ϕ i + 2 ϕ i −1 ), как показано ниже в уравнении.(A.7):

(A.7)

Преобразуя это уравнение, мы получаем формулу центральной разности,

(A.8)

Приложение Б: Калькулятор потенциала – BooTsik.f90

В этом приложении содержится код FORTRAN 90 для потенциального калькулятора BooTsik.f90, который также можно найти по адресу https://github.com/calboo/Potential_Field_Calculator. Пользователь должен установить размер домена: nx , ny , nz ; шаг сетки d ; размер области ввода nx 0, nx 1, ny 0, ny 1, и должен указывать нормальное поле B z 0 в фотосфере.Существуют предустановленные параметры для полей ввода B z 0 , используемых в этой статье, и для импорта поля ввода из неформатированного файла данных, которые можно использовать, просто раскомментировав их. Код приведен ниже:

! BooTsik.f90
! Код FORTRAN 90, предназначенный для расчета потенциального магнитного поля.
! путем экстраполяции нормального поля, заданного в основании области.
!
! Авторы: Каллум Букок и Дэвид Циклаури
! Организация: QMUL
! Электронная почта: [email protected]
! Дата: 3 ноября 2018 г.
!
! Авторские права (C) 2018
!
! Это бесплатное программное обеспечение: вы можете распространять и / или изменять
! это в соответствии с условиями Стандартной общественной лицензии GNU, опубликованной
! Free Software Foundation, либо версия 3 Лицензии, либо
! любая более поздняя версия.
!
! Эта программа распространяется в надежде, что она будет полезна,
! но БЕЗ КАКИХ-ЛИБО ГАРАНТИЙ; без даже подразумеваемой гарантии
! КОММЕРЧЕСКАЯ ЦЕННОСТЬ или ПРИГОДНОСТЬ ДЛЯ КОНКРЕТНОЙ ЦЕЛИ. Увидеть
! Стандартная общественная лицензия GNU для более подробной информации:
! .
!

Константы MODULE
ПОДРАЗУМЕВАЕМЫЙ НЕТ

! Размер области, в которой необходимо рассчитать трехмерное потенциальное поле.
ЦЕЛОЕ, ПАРАМЕТР :: nx = 100, ny = 100, nz = 100
!
! Установите параметр для шага сетки.
ДВОЙНАЯ ТОЧНОСТЬ, ПАРАМЕТР :: d = 0,01
!
! Установить размер области ввода
! Раскомментируйте 1-ю строку для стандарта, 2-ю строку для расширенного
! ДВОЙНАЯ ТОЧНОСТЬ, ПАРАМЕТР :: nx0 = -1, nx1 = nx + 2, ny0 = -1, ny1 = ny + 2
! ДВОЙНАЯ ТОЧНОСТЬ, ПАРАМЕТР :: nx0 = -nx-3, nx1 = 2 * nx + 4, ny0 = -ny-3, ny1 = 2 * ny + 4
!
! Параметры для определения нормального поля Bz0 в основании области z = 0.ДВОЙНАЯ ТОЧНОСТЬ, ПАРАМЕТР :: x01 = 0,0, y01 = 0,0, x02 = 0,0, y02 = 0,5, сигма = 1,0, amp = 1,0
!
! Значение числа пи
ДВОЙНАЯ ТОЧНОСТЬ, ПАРАМЕТР :: pi = 3,141558979323
!
! Массивы и индексы
ДВОЙНАЯ ТОЧНОСТЬ, РАЗМЕР (-3: nx + 4, -3: ny + 4, -3: nz + 4) :: phi
ДВОЙНАЯ ТОЧНОСТЬ, РАЗМЕР (-2: nx + 2, -1: ny + 2, -1: nz + 2) :: bx
ДВОЙНАЯ ТОЧНОСТЬ, РАЗМЕР (-1: nx + 2, -2: ny + 2, -1: nz + 2) :: by
ДВОЙНАЯ ТОЧНОСТЬ, РАЗМЕР (-1: nx + 2, -1: ny + 2, -2: nz + 2) :: bz
РЕАЛЬНОЕ, РАЗМЕР (nx0: nx1, ny0: ny1) :: bz0
ЦЕЛОЕ :: i, j, k, i1, j1

Константы END MODULE
 
ПРОГРАММНЫЙ потенциал
Константы USE
ПОДРАЗУМЕВАЕМЫЙ НЕТ
фи = 0.0

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
! Ниже приведены три варианта поля Bz0, которое необходимо указать на нижней границе.
! Возможные варианты: а) бипольная петля, б) гауссов полюс, в) ввод из файла данных. Для вариантов (а) и
! (б) параметры могут быть установлены в верхней части этого кода. Раскомментируйте любой профиль Bz0
! хотите использовать или, в качестве альтернативы, указать собственное поле Bz0.

! Бипольная петля
! делать j = ny0, ny1
! делать i = nx0, nx1
! bz0 (i, j) = amp * ((1.0 / ((((real (i) - (nx / 2.0)) / (nx / 2.0)) - x01) ** 2 + &
! (((real (j) - (ny / 2.0)) / (ny / 2.0)) - y01) ** 2 + 1.0 ** 2) ** (3.0 / 2.0)) &
! - (1.0 / ((((real (i) - (nx / 2.0)) / (nx / 2.0)) - x02) ** 2 + &
! (((real (j) - (ny / 2.0)) / (ny / 2.0)) - y02) ** 2 + 1.0 ** 2) ** (3.0 / 2.0)))
! конец делать
! конец делать
! Гауссов полюс
! делать j = ny0, ny1
! делать i = nx0, nx1
! bz0 (i, j) = amp * exp (- ((((real (i) - (nx / 2.0)) / (nx / 2.0)) - x01) ** 2.0+ &
! (((real (j) - (ny / 2.0)) / (ny / 2.0)) - y01) ** 2.0 + 1.0 ** 2.0) / (2.0 * sigma ** 2))
! конец делать
! конец делать

! Прочитать в Bz0 из неформатированного файла данных
! Убедитесь, что размер импортированного массива равен размеру bz0, как определено в константах.! OPEN (unit = 12, FORM = ’UNFORMATTED’, STATUS = ’OLD’, ACTION = ’READ’, FILE = ’bz0.dat’)
! ПРОЧИТАТЬ (12) bz0
! ЗАКРЫТЬ (12)

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
! Этот цикл вычисляет магнитный потенциал с использованием метода модифицированной функции Грина.
!
сделать k = -3, nz + 4
сделать j = -3, ny + 4
сделать я = -3, nx + 4
сделать i1 = nx0, nx1
сделать j1 = ny0, ny1
phi (i, j, k) = phi (i, j, k) + bz0 (i1, j1) * d / (2,0 * pi * sqrt ((i-i1) ** 2+ (j-j1) ** 2+ (k + 4 + 1 / sqrt (2.0 * pi)) ** 2))
конец делать
конец делать
конец делать
конец делать
print *, ’Расчет Phi при HIGHT =’, k + 4, ’out of’, nz + 8
конец делать

!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
! В этом разделе рассчитываются компоненты магнитного поля путем численного дифференцирования.
! магнитного потенциала.!
bx = 0,0
по = 0,0
bz = 0,0
сделать k = -1, nz + 2
сделать j = -1, ny + 2
сделать я = -2, nx + 2
bx (i, j, k) = - (phi (i-1, j, k) -27.0 * phi (i, j, k) + 27.0 * phi (i + 1, j, k) -phi (i + 2, к, к)) / (24,0 * г)
конец делать
конец делать
print *, ’Теперь вычисляем Bx на HEIGHT =’, k + 2, ’out of’, nz + 4
конец делать
сделать k = -1, nz + 2
сделать j = -2, ny + 2
сделать я = -1, nx + 2
по (i, j, k) = - (phi (i, j-1, k) -27.0 * phi (i, j, k) + 27.0 * phi (i, j + 1, k) -phi (i, j + 2, k)) / (24,0 * г)
конец делать
конец делать
print *, ’Теперь вычисляется By на HEIGHT =’, k + 2, ’out of’, nz + 4
конец делать
сделать k = -2, nz + 2
сделать j = -1, ny + 2
сделать я = -1, nx + 2
bz (i, j, k) = - (phi (i, j, k-1) -27.0 * phi (i, j, k) +27.0 * phi (i, j, k + 1) -phi (i, j, k + 2)) / (24.0 * d)
конец делать
конец делать
print *, ’Теперь вычисляем Bz на HEIGHT =’, k + 2, ’out of’, nz + 4
конец делать
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!


! Этот раздел записывает компоненты магнитного поля в три файла.
! bx.dat, by.dat и bz.dat
!
ОТКРЫТЬ (unit = 66, FORM = ’UNFORMATTED’, FILE = ’bx.dat’)
ЗАПИСАТЬ (66) bx
ЗАКРЫТЬ (66)
ОТКРЫТЬ (unit = 66, FORM = ’UNFORMATTED’, FILE = ’by.dat’)
ЗАПИСАТЬ (66) по
ЗАКРЫТЬ (66)
ОТКРЫТЬ (unit = 66, FORM = ’UNFORMATTED’, FILE = ’bz.dat ’)
ЗАПИСАТЬ (66) bz
ЗАКРЫТЬ (66)
КОНЕЦ ПРОГРАММЫ
 

Этот код можно запускать параллельно в кластере HPC. Для этого код можно автоматически распараллелить, используя, например, компилятор Intel Fortran. Используя компилятор Intel Fortran, следующая строка параллельно компилирует код, готовый к отправке в качестве задания на узел HPC. Это значительно улучшает время выполнения в зависимости от количества используемых ядер. Строка записывается как

ifort -o потенциал -O3 -параллельный потенциал. F90 -xHost.

Приложение C: Инструмент визуализации IDL – Boobox.pro

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

Для использования инструмента пользователь должен указать размеры домена: nx , ny , nz ; экстенты домена для меток галочки и формат метки галочки; углы, под которыми следует визуализировать область, ax , az , и, наконец, необходимо указать значение для переменной GRID.Если для параметра GRID задано целочисленное значение 1, то линии поля проводятся из начальных точек, равномерно распределенных по поверхности фотосферы в основании домена. Если для GRID установлено любое другое значение, то линии поля строятся из начальных точек, случайным образом распределенных по всей области, создавая график, на котором плотность линий поля примерно соответствует напряженности поля. Примеры графиков, построенных для биполярной петли с использованием Boobox.pro с GRID, установленным в 1 и 0, приведены на рис.С.1 и С.2.

Рис. С.1.

Линейная диаграмма потенциального бипольного поля, визуализированная с помощью Boobox.pro с GRID = 1. Мы можем ясно видеть силовые линии магнитного поля, соединенные с каждой точкой фотосферы.

Открыть с помощью DEXTER
Рис. С.2.

Линейная диаграмма потенциального бипольного поля, визуализированная с помощью Boobox.pro с GRID = 0. Мы можем ясно видеть, что напряженность поля больше ближе к полюсам на фотосфере.

Открыть с помощью DEXTER

Если инструмент используется с файлами неструктурированных данных компонентов магнитного поля bx.dat, by.dat и bz.dat, то соответствующий блок необходимо раскомментировать, и инструмент можно запустить в IDL с помощью .r boobox. Если инструмент используется для визуализации данных из моментального снимка данных Lare3d, ds, тогда первая строка pro boobox, ds и соответствующий блок должны быть раскомментированы, и инструмент может быть запущен в IDL с boobox, ds, где ds – это снимок данных.Скрипт IDL для Boobox.pro приведен ниже:

; Boobox.pro
;
; Скрипт IDL для визуализации магнитных полей либо из файлов с неструктурированными данными, либо из снимков данных из Lare3d.
Поля визуализируются в виде линий поля в трехмерном блоке с основанием блока, очерченным в соответствии с силой и направлением Bz на этой границе.
;
; Автор: Каллум Букок
; Организация: QMUL
; Почта: [email protected]
; Дата: 3 ноября 2018 г.
;
; Авторские права (C) 2018
;
; Это бесплатное программное обеспечение: вы можете распространять и / или изменять
; это в соответствии с условиями Стандартной общественной лицензии GNU, опубликованной
; Free Software Foundation, либо версия 3 Лицензии, либо
; любая более поздняя версия.;
; Эта программа распространяется в надежде, что она будет полезна,
; но БЕЗ КАКИХ-ЛИБО ГАРАНТИЙ; без даже подразумеваемой гарантии
; КОММЕРЧЕСКАЯ ЦЕННОСТЬ или ПРИГОДНОСТЬ ДЛЯ КОНКРЕТНОЙ ЦЕЛИ. Увидеть
; Стандартная общественная лицензия GNU для более подробной информации:
; .
;
; Раскомментируйте следующую строку только при использовании ds моментальных снимков данных из Lare3d
; pro boobox, ds

; Размеры домена
nx = 100,000
ny = 100,000
nz = 100,000
; Экстенты домена, для галочки
xmin = -1,0
xmax = 1,0
ymin = -1.0
ymax = 1.0
zmin = 0,0
zmax = 2,0
; Формат меток галочки
tikform = ’(f5.1)’
; Углы обзора для визуализации
топор = 30
аз = 30
; Расположение силовых линий
; 1 - начальные точки, равномерно распределенные на фотосфере
; Любое другое значение - начальные точки случайным образом распределяются по домену.
СЕТКА = 0

; Раскомментируйте этот блок для использования с файлами .dat
;
; bx0 = dblarr (nx + 5, ny + 4, nz + 4)
; by0 = dblarr (nx + 4, ny + 5, nz + 4)
; bz0 = dblarr (nx + 4, ny + 4, nz + 5)
; х = dblarr (nx + 4)
; y = dblarr (ny + 4)
; z = dblarr (nz + 4)
; OPENR, 1, ’bx.dat ’, / F77_UNFORMATTED
; READU, 1, bx0
; ЗАКРЫТЬ, 1
; OPENR, 1, ’by.dat’, / F77_UNFORMATTED
; READU, 1, by0
; ЗАКРЫТЬ, 1
; OPENR, 1, ’bz.dat’, / F77_UNFORMATTED
; READU, 1, bz0
; ЗАКРЫТЬ, 1
; для i = 0, nx + 3 do x (i) = i
; для i = 0, ny + 3 do y (i) = i
; для i = 0, nz + 3 do z (i) = i
; bx0 = bx0 (2: nx + 3,2: ny + 2,2: nz + 2)
; by0 = by0 (2: nx + 2,2: ny + 3,2: nz + 2)
; bz0 = bz0 (2: nx + 2,2: ny + 2,2: nz + 3)

; Раскомментируйте этот блок для использования с ds моментальных снимков данных из Lare3d
;
; bx0 = ds.bx
; by0 = ds.by
; bz0 = ds.bz
; x = ds.x
; y = ds.y
; z = ds.z

; Рассчитайте и отформатируйте метки галочки
x2 = (xmin + xmax) / 2.0
х1 = (хмин + х2) /2.0
x3 = (x2 + xmax) /2.0
y2 = (ymin + ymax) /2.0
y1 = (ymin + y2) /2.0
y3 = (y2 + ymax) /2.0
z2 = (zmin + zmax) /2.0
z1 = (zmin + z2) /2,0
z3 = (z2 + zmax) /2.0
xmin = строка (xmin, FORMAT = tikform)
x2 = строка (x2, FORMAT = tikform)
x1 = строка (x1, FORMAT = tikform)
x3 = строка (x3, FORMAT = tikform)
xmax = строка (xmax, FORMAT = tikform)
ymin = строка (ymin, FORMAT = tikform)
y2 = строка (y2, FORMAT = tikform)
y1 = строка (y1, FORMAT = tikform)
y3 = строка (y3, FORMAT = tikform)
ymax = строка (ymax, FORMAT = tikform)
zmin = строка (zmin, FORMAT = tikform)
z2 = строка (z2, FORMAT = tikform)
z1 = строка (z1, FORMAT = tikform)
z3 = строка (z3, FORMAT = tikform)
zmax = строка (zmax, FORMAT = tikform)
xt = [strtrim (xmin, 2), strtrim (x1,2), strtrim (x2,2), strtrim (x3,2), strtrim (xmax, 2)]
yt = [strtrim (ymin, 2), strtrim (y1,2), strtrim (y2,2), strtrim (y3,2), strtrim (ymax, 2)]
zt = [strtrim (zmin, 2), strtrim (z1,2), strtrim (z2,2), strtrim (z3,2), strtrim (zmax, 2)]

; Установить цвета
ТВЛКТ, 255, 255, 255, 254; белый цвет
TVLCT, 0, 0, 0, 253; Черный цвет
ТВЛКТ, 88, 88, 88, 252; Цвет линии
!П.Цвет = 253
! P.Background = 254

; Чистое окно
WINDOW, 1, XSIZE = 600, YSIZE = 500, TITLE = ’Boobox’
СТЕРЕТЬ

; Настройте систему 3D-масштабирования:
SCALE3, xr = [0, nx], yr = [0, ny], zr = [0, nz], ax = ax, az = az

; Нарисуйте контур основы:
КОНТУР, bz0 (*, *, 0), / fill, nlevels = 25, / t3d, / noerase, zvalue = 0.0, / noclip, $
XSTYLE = 1, YSTYLE = 1, XRANGE = [0, nx], YRANGE = [0, ny], CHARSIZE = 4, $
XTITLE = ’X’, YTITLE = ’Y’, POS = [0.1,0.1, nx, ny], $
XTICKS = 4, XTICKNAME = xt, $
YTICKS = 4, YTICKNAME = yt, $
ZTICKS = 4, ZTICKNAME = zt

; Задайте трехмерное координатное пространство с помощью осей.
ПОВЕРХНОСТЬ, [[0, nx], [0, ny]], / NODATA, / SAVE, / noerase, / t3d, $
XSTYLE = 1, YSTYLE = 1, / noclip, XRANGE = [0, nx], YRANGE = [0, ny], ZRANGE = [0, nz], $
CHARSIZE = 4, POS = [0.1,0.1, nx, ny], $
XTICKS = 4, XTICKNAME = xt, $
YTICKS = 4, YTICKNAME = yt, $
ZTICKS = 4, ZTICKNAME = zt

; Постройте векторное поле
! P.Color = 252
ЕСЛИ (СЕТКА уравн. 1) ТО НАЧАТЬ
к = 0
п = 0
sx = FLTARR (1000)
sy = FLTARR (1000)
sz = FLTARR (1000)
для j = 0, ny-1, (ny / 20.0) действительно начинаются
   для i = 0, nx-1, (nx / 20.0) действительно начинаются
      sx (n) = я
      sy (n) = j
      sz (n) = k
      п = п + 1
   конец
конец
ENDIF ELSE НАЧАТЬ
seed = 12345678
sx = FIX (nx * RANDOMU (начальное число, 250))
sy = FIX (ny * RANDOMU (семя; 250))
sz = FIX (nz * RANDOMU (начальное число, 250))
ENDELSE
SCALE3, zr = [0, nz], ax = ax, az = az
ПОТОК3, bx0, by0, bz0, ARROWSIZE = 0.01, LEN = 2.0, NSTEPS = 9999, sx = sx, sy = sy, sz = sz
FLOW3, -bx0, -by0, -bz0, ARROWSIZE = 0,01, LEN = 2,0, NSTEPS = 9999, sx = sx, sy = sy, sz = sz
! P.Color = 253

; Раскомментируйте для создания файла png
; write_png, ’Fields.png’, tvrd (/ true)

конец
 

Приложение D: Импорт данных потенциального поля в Lare3d

В этом приложении объясняется, как данные о потенциальном поле, созданные BooTsik.f90 и сохраненные в трех файлах, bx.dat, by.dat и bz.dat, были импортированы в Lare3d и использованы в качестве начальной конфигурации магнитного поля.Сначала файлы данных bx.dat, by.dat и bz.dat были скопированы в тот же каталог, из которого был запущен Lare3d. Затем были внесены изменения в четыре исходных файла Lare3d.

Первым редактируемым исходным файлом был shared_ data.F90, который можно найти по пути src / core / shared_ data.F90. Следующие строки были добавлены в модуль shared_data для создания размещаемых массивов:

! Дополнительные размещаемые массивы для:
! Сохранение входных данных магнитного поля
! Сохранение начальных значений для переменных моделирования
! Расчет магнитных полей в центрах ячеек

  REAL (число), DIMENSION (:,:, :), ALLOCATABLE :: bx_init, by_init, bz_init
  REAL (число), DIMENSION (:,:, :), ALLOCATABLE :: rho0, energy0, bx0, by0, bz0
  REAL (число), DIMENSION (:,:, :), ALLOCATABLE :: arrx, arry, arrz
 

Вторым редактируемым исходным файлом был control.f90, который можно найти по пути src / control.f90. В этот файл были внесены следующие изменения. В подпрограмме user_normalisation все нормализации были установлены на 1.0. В подпрограмме control_variables количество точек сетки и размер домена были изменены, чтобы соответствовать нашим входным данным, время выполнения было установлено на 100.0 (это измерено в τ A ), удельное сопротивление было установлено на 5 × 10 −5 , а граничные условия были заданы пользователем BC_USER.Затем в конец подпрограммы control_variables были добавлены следующие строки для считывания данных потенциального магнитного поля в выделяемые массивы bx_init, by_init и bz_init:

   ! Здесь мы размещаем массивы bx_init, by_init и bz_init
   ! Затем мы читаем файлы bx.dat, by.dat и bz.dat.
   ! и сохраните входные магнитные поля в этих массивах.

   ALLOCATE (bx_init (-2: nx_global + 2, -1: ny_global + 2, -1: nz_global + 2))
   ALLOCATE (by_init (-1: nx_global + 2, -2: ny_global + 2, -1: nz_global + 2))
   ALLOCATE (bz_init (-1: nx_global + 2, -1: ny_global + 2, -2: nz_global + 2))

   ОТКРЫТЬ (unit = 12, FORM = ’UNFORMATTED’, STATUS = ’OLD’, ACTION = ’READ’, FILE = ’bx.dat ’)
   ЧИТАТЬ (12) bx_init (-2: nx_global + 2, -1: ny_global + 2, -1: nz_global + 2)

   ОТКРЫТЬ (unit = 12, FORM = ’UNFORMATTED’, STATUS = ’OLD’, ACTION = ’READ’, FILE = ’by.dat’)
   ПРОЧИТАЙТЕ (12) by_init (-1: nx_global + 2, -2: ny_global + 2, -1: nz_global + 2)

   ОТКРЫТЬ (unit = 12, FORM = ’UNFORMATTED’, STATUS = ’OLD’, ACTION = ’READ’, FILE = ’bz.dat’)
   ЧИТАТЬ (12) bz_init (-1: nx_global + 2, -1: ny_global + 2, -2: nz_global + 2)
 

Наконец, в подпрограмме set_output_dumps время между моментальными снимками было установлено на 10.0 (это измерено в τ A ), а маски дампа 17-19 для токов были установлены на.ИСТИНА ..

Третьим редактируемым исходным файлом был initial_conditions.f90, который можно найти по пути src / initial_conditions.f90. Файл initial_conditions.f90 устанавливает начальные значения плотности ρ , энергии ε , поля скоростей v и магнитного поля B для каждой подобласти после декомпозиции области. В этом файле мы заменяем подпрограмму set_initial_conditions следующей:

  SUBROUTINE set_initial_conditions

  ЦЕЛОЕ :: ix, iy, iz
  REAL (число) :: бета
! Гравитация и бета

  грав = 0.0_num
  бета = 1.0e-8_num

! Скорости
! Статический домен.

  vx = 0.0_num
  vy = 0.0_num
  vz = 0.0_num

! Магнитное поле
! Компоненты магнитного поля считываются с массивов
! bx_init, by_init, bz_init.

  DO ix = -2, nx + 2
     DO iy = -1, ny + 2
        ДО iz = -1, nz + 2
           bx (ix, iy, iz) = bx_init (ix + n_global_min (1), iy + n_global_min (2), iz + n_global_min (3))
        КОНЕЦ ДЕЛАТЬ
     КОНЕЦ ДЕЛАТЬ
   КОНЕЦ ДЕЛАТЬ
   DO ix = -1, nx + 2
      DO iy = -2, ny + 2
         ДО iz = -1, nz + 2
            by (ix, iy, iz) = by_init (ix + n_global_min (1), iy + n_global_min (2), iz + n_global_min (3))
         КОНЕЦ ДЕЛАТЬ
      КОНЕЦ ДЕЛАТЬ
   КОНЕЦ ДЕЛАТЬ
   
   DO ix = -1, nx + 2
      DO iy = -1, ny + 2
         ДО iz = -2, nz + 2
            bz (ix, iy, iz) = bz_init (ix + n_global_min (1), iy + n_global_min (2), iz + n_global_min (3))
         КОНЕЦ ДЕЛАТЬ
      КОНЕЦ ДЕЛАТЬ
   КОНЕЦ ДЕЛАТЬ

! Плотность

! В этом случае плотность была установлена ​​равной
! квадрат напряженности магнитного поля.ALLOCATE (arrx (-1: nx + 2, -1: ny + 2, -1: nz + 2))
  ALLOCATE (arry (-1: nx + 2, -1: ny + 2, -1: nz + 2))
  ALLOCATE (arrz (-1: nx + 2, -1: ny + 2, -1: nz + 2))

  DO iy = -1, ny + 2
     ДО iz = -1, nz + 2
        DO ix = -1, nx + 2
            arrx (ix, iy, iz) = 0,5 * bx (ix-1, iy, iz) + bx (ix, iy, iz)
            arry (ix, iy, iz) = 0,5 * по (ix, iy-1, iz) + по (ix, iy, iz)
            arrz (ix, iy, iz) = 0,5 * bz (ix, iy, iz-1) + bz (ix, iy, iz)
            rho (ix, iy, iz) = (arrx (ix, iy, iz) ** 2.0) + &
                           (arry (ix, iy, iz) ** 2.0) + &
                           (arrz (ix, iy, iz) ** 2.0)
        КОНЕЦ ДЕЛАТЬ
     КОНЕЦ ДЕЛАТЬ
  КОНЕЦ ДЕЛАТЬ

! Энергия
! Энергия задана таким образом, чтобы давление везде было равным бета / 2.

  энергия = 0,5_число * (бета * 1.0_число) / ((rho) * (гамма-1.0_число))

! Сохранить начальные значения
! В этих массивах хранятся начальные значения.

  ALLOCATE (rho0 (-1: nx + 2, -1: ny + 2, -1: nz + 2))
  ALLOCATE (bx0 (-2: nx + 2, -1: ny + 2, -1: nz + 2))
  ALLOCATE (by0 (-1: nx + 2, -2: ny + 2, -1: nz + 2))
  ALLOCATE (bz0 (-1: nx + 2, -1: ny + 2, -2: nz + 2))
  ALLOCATE (энергия0 (-1: nx + 2, -1: ny + 2, -1: nz + 2))

  bx0 = bx
  по0 = по
  bz0 = bz
  rho0 = rho
  energy0 = энергия

  КОНЕЦ ПОДПРОГРАММЫ set_initial_conditions
 

Для каждой подобласти эта подпрограмма обнуляет гравитацию и начальное поле скорости.Подпрограмма устанавливает начальную внутреннюю энергию ε так, чтобы начальное давление P составляло 5 × 10 −9 поперек области, импортирует правильное подразделение потенциального магнитного поля и устанавливает плотность, равную квадрату величины магнитного поля. напряженность магнитного поля ρ = B 2 . Наконец, эти начальные значения сохраняются в массивах rho0, bx0, by0, bz0, energy0.

Четвертым и последним исходным файлом, подлежащим редактированию, был Border.f90, который можно найти по пути src / initial_boundary.f90. Эти граничные условия в этом файле были отредактированы так, чтобы на каждой границе плотность ρ , внутренняя энергия ε и компоненты магнитного поля bx, by, bz были установлены на свои начальные значения rho0, bx0, by0, bz0, energy0 . Пример этого приведен ниже, пример приведен для магнитных полей с минимальным значением x – граница:

  ЕСЛИ (proc_x_min == MPI_PROC_NULL .AND. Xbc_min == BC_USER) ТО
    bx (-1,:, :) = bx0 (-1,:, :)
    bx (-2,:, :) = bx0 (-2,:, :)
    by (0,:, :) = by0 (0,:, :)
    by (-1,:, :) = by0 (-1,:, :)
    bz (0,:, :) = bz0 (0,:, :)
    bz (-1,:, :) = bz0 (-1,:, :)
  КОНЕЦ ЕСЛИ
 

Наконец, скорости на всех границах как на полном, так и на половинном временном шаге были установлены на ноль.Пример этого приведен ниже, пример приведен для скоростей полного временного шага при минимальной границе x :

  ЕСЛИ (proc_x_min == MPI_PROC_NULL .AND. Xbc_min == BC_USER) ТО
    vx (-2: 0,:, :) = 0.0_num
    vy (-2: 0,:, :) = 0.0_num
    vz (-2: 0,:, :) = 0.0_num
  КОНЕЦ ЕСЛИ
 

Все таблицы

Таблица 1.

Значение каждой метрики для метода A, метода B и метода C численного дифференцирования для контрольного примера 1.

Таблица 2.

Значение каждой метрики для метода A, метода B и метода C численного дифференцирования для контрольного примера 2.

Таблица 3.

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

Таблица 4.

Значение каждой метрики с использованием стандартных и расширенных областей ввода для тестового примера 2.

Таблица 5.

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

Таблица 6.

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

Все фигуры

Инжир.1.

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

Открыть с помощью DEXTER
По тексту
Рис. 2.

Линейная диаграмма магнитного поля, используемого для Контрольного примера 1.Фотосферная поверхность в основании домена представляет собой цветной контур B z 0 .

Открыть с помощью DEXTER
По тексту
Рис 3.

Линейная диаграмма магнитного поля, используемого в тестовом примере 2. Фотосферная поверхность в основании домена представляет собой цветной контур B z 0 .

Открыть с помощью DEXTER
По тексту
Инжир.4.

Иллюстрация расчетной области для ϕ , показанной серым прямоугольником, и области ввода для B z 0 , показанной синим квадратом, когда ( Nx 1 , Ny 2 ) равно ( Nx , Ny ).

Открыть с помощью DEXTER
По тексту
Рис. 5.

Иллюстрация расчетной области для ϕ , показанной серым прямоугольником, и области ввода для B z 0 , показанной синим квадратом, когда ( Nx 1 , Ny 2 ) равно (3 Nx , 3 Ny ).

Открыть с помощью DEXTER
По тексту
Рис. 6.

Закрашенная поверхностная диаграмма напряженности магнитного поля при z = 0 для потенциального поля, созданного с использованием стандартной области ввода. Расширяющиеся края возникают из-за неоднородности на краю входной области.

Открыть с помощью DEXTER
По тексту
Инжир.7.

Закрашенная поверхностная диаграмма напряженности магнитного поля при z = 0 для потенциального поля, созданного с использованием расширенной области ввода. Края более гладкие из-за расширенной области ввода.

Открыть с помощью DEXTER
По тексту
Рис. 8.

Контурная диаграмма исходных данных магнитограммы нормального магнитного поля на фотосфере. Цвет указывает направление и силу магнитного поля.

Открыть с помощью DEXTER
По тексту
Рис.9.

Схема силовых линий потенциального магнитного поля, созданного с использованием данных солнечной магнитограммы. Фотосферная поверхность в основании домена представляет собой цветной контур B z 0 .

Открыть с помощью DEXTER
По тексту
Инжир.10.

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

Открыть с помощью DEXTER
По тексту
Рис 11.

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

Открыть с помощью DEXTER
По тексту
Рис. А.1.

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

Открыть с помощью DEXTER
По тексту
Инжир.С.1.

Линейная диаграмма потенциального бипольного поля, визуализированная с помощью Boobox.pro с GRID = 1. Мы можем ясно видеть силовые линии магнитного поля, соединенные с каждой точкой фотосферы.

Открыть с помощью DEXTER
По тексту
Рис. С.2.

Линейная диаграмма потенциального бипольного поля, визуализированная с помощью Boobox.pro с GRID = 0. Мы можем ясно видеть, что напряженность поля больше ближе к полюсам на фотосфере.

Открыть с помощью DEXTER
По тексту

Разгон

Ускорение – это скорость изменения скорости:

  • Ускорение
  • Замедление (также называется замедлением )
  • Изменение направления
  • и т. Д.

Обычно отображается как:

м / с 2
“метров на секунду в квадрате”

Что это за “ на секунду в квадрате “? Пример поможет:

Бегун ускоряется с

5 м / с (5 метров в секунду) до 6 м / с всего за одну секунду

Значит, они ускоряются на 1 метр в секунду в секунду

Видите, как «в секунду» используется дважды?

Его можно представить как (м / с) / с, но обычно пишут м / с 2

Значит, их ускорение 1 м / с 2

Формула:

Ускорение = Изменение скорости (м / с) Время (с)

Пример: велогонка!

Вы участвуете в велогонке со стабильной скоростью 10 метров в секунду (10 м / с).

Acceleration : Теперь вы начинаете ездить быстрее! Вы увеличиваете до 14 м / с в течение следующих 2 секунд (все еще движетесь в том же направлении):

Ваша скорость увеличивается на 4 м / с , за период времени 2 секунды , поэтому:

Ускорение = Изменение скорости (м / с) Время (с)

= 4 м / с 2 с = 2 м / с 2

Ваша скорость изменится на 2 метра в секунду в секунду .

Пример: вы бежите со скоростью 7 м / с, а занос остановился за 2 секунды.

Вы перешли с 7 м / с на 0, так что это уменьшение скорости:

Ускорение = Изменение скорости (м / с) Время (с)

= −7 м / с 2 с = −3,5 м / с 2

Мы не всегда это говорим, но ускорение имеет направление (что делает его вектором):

Автомобиль движется на

Запад со скоростью 16 м / с .

Водитель поворачивает руль, и в течение 4 секунд машина направляется на восток со скоростью 16 м / с . Какое ускорение?

Цифры такие же, но направление другое!

Ускорение = Изменение скорости (м / с) Время (с)

Ускорение = С 16 м / с на запад до 16 м / с на восток 4 с

С 16 м / с на запад до 16 м / с на восток общее изменение на 32 м / с к востоку.

Ускорение = 32 м / с Восток 4 с = 8 м / с 2 Восток

Для более сложных изменений направления читайте векторы.

.

Оставить комментарий

Ваш адрес email не будет опубликован. Обязательные поля помечены *