Файл:Mie resonances vs Radius.gif
Материал из Википедии — свободной энциклопедии
Перейти к навигации
Перейти к поиску
Mie_resonances_vs_Radius.gif (243 × 243 пкс, размер файла: 5,67 МБ, MIME-тип: image/gif, закольцованный, 202 фрейма, 30 с)
Этот файл находится на Викискладе. Сведения о нём показаны ниже.
Викисклад — централизованное хранилище для свободных файлов, используемых в проектах Викимедиа.
Сообщить об ошибке с файлом |
Краткое описание
ОписаниеMie resonances vs Radius.gif |
English: Effect of a dielectric sphere (technically a disk, as the simulation is in 2D) on an incident plane wave as a function of the radius.
The patterns you see flashing in are the Mie resonances. The incident plane wave is coming from the bottom. |
Дата | |
Источник | https://twitter.com/j_bertolotti/status/1388821332989775875 |
Автор | Jacopo Bertolotti |
Права (Повторное использование этого файла) |
https://twitter.com/j_bertolotti/status/1030470604418428929 |
Mathematica 12.0 code
\[Lambda]0 = 1.; k0 =
N[(2 \[Pi])/\[Lambda]0]; (*The wavelength in vacuum is set to 1, so all lengths are now in units of wavelengths*)
\[Delta] = \[Lambda]0/20; \[CapitalDelta] = 20*\[Lambda]0; (*Parameters for the grid*)
sourcef[x_, y_] := E^(I k0 y); (*Important! The source MUST be a solution of the Helmholtz equation in vacuum*)
\[Phi]in = Table[sourcef[x, y], {x, -\[CapitalDelta]/2, \[CapitalDelta]/ 2, \[Delta]}, {y, -\[CapitalDelta]/2, \[CapitalDelta]/ 2, \[Delta]}]; (*Discretized source*)
d = \[Lambda]0/1; (*typical scale of the absorbing layer*)
imn = Table[5 I (E^-((x + \[CapitalDelta]/2)/d) + E^((x - \[CapitalDelta]/2)/d) + E^-((y + \[CapitalDelta]/2)/d) + E^((y - \[CapitalDelta]/2)/d)), {x, -\[CapitalDelta]/2, \[CapitalDelta]/2, \[Delta]}, {y, -\[CapitalDelta]/ 2, \[CapitalDelta]/ 2, \[Delta]}]; (*Imaginary part of the refractive index (used to emulate absorbing boundaries)*)
dim = Dimensions[\[Phi]in][[1]];
L = -1/\[Delta]^2*KirchhoffMatrix[GridGraph[{dim, dim}]]; (*Discretized Laplacian*)
r[t_] := (\[Lambda]0*3 - \[Delta]*5) Sin[\[Pi]/ 2 t]^2 + \[Delta]*5; (*The radius changes with the parameter t*)
frames = Table[
n = Table[
If[y^2 + x^2 <= r[t]^2, 2, 1], {x, -\[CapitalDelta]/ 2, \[CapitalDelta]/2, \[Delta]}, {y, -\[CapitalDelta]/ 2, \[CapitalDelta]/2, \[Delta]}] + imn; (*Matrix with the refractive index*)
b = -(Flatten[n]^2 - 1) k0^2 Flatten[\[Phi]in]; (*Right-hand side of the equation we want to solve*)
M = L + DiagonalMatrix[SparseArray[Flatten[n]^2 k0^2]]; (*Operator on the left-hand side of the equation we want to solve*)
\[Phi]s = Partition[LinearSolve[M, b], dim]; (*Solve the linear system*)
ImageAdd[
ArrayPlot[
Transpose[Abs[\[Phi]in + \[Phi]s]^2/ Max[Abs[\[Phi]in + \[Phi]s]^2]][[(4 d)/\[Delta] ;; (-4 d)/\[Delta], (4 d)/\[Delta] ;; (-4 d)/\[Delta]]], ColorFunction -> "AvocadoColors" , DataReversed -> True,
Frame -> False, PlotRange -> {0, 0.8}], ArrayPlot[Transpose@Re[(n - 1)/10] , DataReversed -> True , ColorFunctionScaling -> False, ColorFunction -> GrayLevel, Frame -> False]
](*Plot everything*)
, {t, 0, 1, 1./100}];
ListAnimate[frames]
Лицензирование
Я, владелец авторских прав на это произведение, добровольно публикую его на условиях следующей лицензии:
Этот файл доступен на условиях Creative Commons CC0 1.0 Универсальной передачи в общественное достояние (Universal Public Domain Dedication). | |
Лица, связанные с работой над этим произведением, решили передать данное произведение в общественное достояние, отказавшись от всех прав на произведение по всему миру в рамках закона об авторских правах (а также связанных и смежных прав), в той степени, которую допускает закон. Вы можете копировать, изменять, распространять, исполнять данное произведение в любых целях, в том числе в коммерческих, без получения на это разрешения автора.
http://creativecommons.org/publicdomain/zero/1.0/deed.enCC0Creative Commons Zero, Public Domain Dedicationfalsefalse |
Элементы, изображённые на этом файле
изображённый объект
У этого свойства есть некоторое значение без элемента в
image/gif
История файла
Нажмите на дату/время, чтобы посмотреть файл, который был загружен в тот момент.
Дата/время | Миниатюра | Размеры | Участник | Примечание | |
---|---|---|---|---|---|
текущий | 08:51, 4 мая 2021 | 243 × 243 (5,67 МБ) | Berto | Uploaded own work with UploadWizard |
Использование файла
Следующая страница использует этот файл:
Глобальное использование файла
Данный файл используется в следующих вики:
- Использование в en.wikipedia.org
Метаданные
Файл содержит дополнительные данные, обычно добавляемые цифровыми камерами или сканерами. Если файл после создания редактировался, то некоторые параметры могут не соответствовать текущему изображению.
Примечание GIF-файла | Created with the Wolfram Language : www.wolfram.com |
---|