@NilsBerglund

@NilsBerglund

Nils Berglund | Exciting resonant modes in a circle @NilsBerglund | Uploaded September 2024 | Updated October 2024, 1 hour ago.
This video is inspired by some recent simulations on this channel, showing resonators shaped like magnetrons. These resonators feature several circular cavities, in which the waves displayed interesting standing wave patterns. In this simulation, the source of the waves is located in a smaller circular cavity, and the waves propagate through a channel to the larger cavity, where they produce interference patterns. Note that some of these patterns are reminiscent of heart-shaped reflections that you can sometimes see in a cup of coffee.
This video has two parts, showing the same evolution with two different color gradients:
Wave height: 0:00
Averaged wave energy: 1:39
In the first part, the color hue depends on the height of the wave. In the second part, it depends on the energy of the wave, slightly averaged over a sliding time window.

Render time: 21 minutes 47 seconds
Compression: crf 23
Color scheme: Part 1 - Twilight by Bastian Bechtoldgithub.com/bastibe/twilight
Part 2 - Magma by Nathaniel J. Smith and Stefan van der Walt
github.com/BIDS/colormap

Music: "Highway 5" by TrackTribe@TrackTribe

See also
https://images.math.cnrs.fr/des-ondes-dans-mon-billard-partie-i/ for more explanations (in French) on a few previous simulations of wave equations.

The simulation solves the wave equation by discretization. The algorithm is adapted from the paper hplgit.github.io/fdm-book/doc/pub/wave/pdf/wave-4print.pdf
C code: github.com/nilsberglund-orleans/YouTube-simulations
https://www.idpoisson.fr/berglund/software.html
Many thanks to Marco Mancini and Julian Kauth for helping me to accelerate my code!

#wave #resonator

Exciting resonant modes in a circle

$Waves escaping a ring of obstacles: Polar grid This is a simulation of waves originating from a point source crossing a set of circular obstacles placed on a polar grid, meaning that the locations form a regular grid in polar coordinates. Here there are 60 radial lines with 25 discs each. As one would expect, waves pass this lattice easier than for other arrangements seen on this channel. Still, diffraction has a large impact on the outgoing waves. The fact that the 60-fold symmetry is not preserved is probably due to round-off errors coming from the discretization. This video has two parts, showing the same evolution with two different color gradients: Averaged wave energy: 0:00 Wave height: 1:14 In the first part, the color hue depends on the energy of the wave, averaged over a sliding time window. In the second part, it depends on the height of the wave. The contrast has been enhanced by a shading procedure, similar to the one I have used on videos of reaction-diffusion equations. The process is to compute the normal vector to a surface in 3D that would be obtained by using the third dimension to represent the field, and then to make the luminosity depend on the angle between the normal vector and a fixed direction. The color in the central region has been brightened to white, because the waves tend to make large-amplitude oscillations there, which would not be very pleasant to watch. There are absorbing boundary conditions on the borders of the simulated rectangle. The display at the right shows a time-averaged version of the signal near the right boundary of the simulated rectangular area. More precisely, it shows the square root of an average of squares of the respective field value (wave height or energy). Render time: 30 minutes 38 seconds Compression: crf 23 Color scheme: Part 1 - Magma by Nathaniel J. Smith and Stefan van der Walt https://github.com/BIDS/colormap Part 2 - Twilight by Bastian Bechtold https://github.com/bastibe/twilight Music: Desert Planet by Quincas Moreira@QuincasMoreira See also https://images.math.cnrs.fr/des-ondes-dans-mon-billard-partie-i/ for more explanations (in French) on a few previous simulations of wave equations. The simulation solves the wave equation by discretization. The algorithm is adapted from the paper https://hplgit.github.io/fdm-book/doc/pub/wave/pdf/wave-4print.pdf C code: https://github.com/nilsberglund-orleans/YouTube-simulations https://www.idpoisson.fr/berglund/software.html Many thanks to Marco Mancini and Julian Kauth for helping me to accelerate my code! #wave #diffraction$ $The double peak disappears when the source filmed with a gradient index lens is slower The recent video https://youtu.be/idB6qn802lM showed a source of light moving at about 33% of the speed of light being filmed by a gradient index lens. It showed a somewhat strange double peak, or echo being formed by the lens, making the source appear doubled. The question arose in the comments whether this was due to the very fast speed of the source. This simulation shows a similar situation, with a source moving at a slower speed, about 17% of the speed of light. It turns out that there is no double peak here, confirming that in the other simulation it must have been an effect of the higher speed. The refractive index n(r) of the gradient index lens used here depends on the distance r to the horizontal symmetry axis of the lens in a quadratic way, like n(r) = n0 - a*r², with n0 = 1.25 and a = 0.4375, making the refractive index smaller than 1 at the outer boundary of the lens, where r = 1. In order to visualize the refractive index as well, the luminosity of the background depends on the refractive index. The videos are inspired by Huygens Optics recent short https://youtube.com/shorts/VGd3Ajnp6e0 showing the principle of a gradient index lens. Lenses focus incoming rays of light by delaying them more near the center of the lens than at its periphery. This is often done with a material of constant index of refraction, by making the lens thicker near the center, as shown for instance in the simulation https://youtu.be/rrJJBh9ubUE . However, one can also build lenses of constant thickness, by making the index of refraction of their material depend on the location in the lens. In this simulation, the index decreases like the square of the distance to the center (that is, it is of the form n0 - a*r², where r is the distance to the axis of symmetry). This results in the incoming waves being focused at two points in the (estimated) focal plane, marked by a vertical line. The plot to the right shows a time-averaged value of the field along that plane. This video has two parts, showing the same evolution with two different color gradients: Wave height: 0:00 Averaged wave energy: 2:16 In the first part, the color hue depends on the height of the wave. In the second part, it depends on the energy of the wave, averaged over a sliding time window. There are absorbing boundary conditions on the borders of the simulated rectangle. The display at the right shows the signal along the focal plane, which is indicated by a vertical line. Render time: 1 hour 6 minutes Compression: crf 23 Color scheme: Part 1 - Viridis, by Nathaniel J. Smith, Stefan van der Walt and Eric Firing Part 2 - Plasma by Nathaniel J. Smith and Stefan van der Walt https://github.com/BIDS/colormap Music: Lazy Boy Blues by the Unicorn Heads@UnicornHeads See also https://images.math.cnrs.fr/Des-ondes-dans-mon-billard-partie-I.html for more explanations (in French) on a few previous simulations of wave equations. The simulation solves the wave equation by discretization. The algorithm is adapted from the paper https://hplgit.github.io/fdm-book/doc/pub/wave/pdf/wave-4print.pdf C code: https://github.com/nilsberglund-orleans/YouTube-simulations https://www.idpoisson.fr/berglund/software.html Many thanks to Marco Mancini and Julian Kauth for helping me to accelerate my code! #wave #lens #gradient_index$

Vortex dynamics on the sphere in a strong Coriolis force

Have you ever questioned the nature of your reality? First successful DNA replication in silico

Within cells interlinked: Replicated DNA curling up in a box

Density, speed and vorticity of unstable vortices on the sphere, featuring music by State Azure

$Waves branching through Poisson disc obstacles This is a simulation of waves originating from a point source crossing a set of obstacles forming a Poisson disc process, also known as blue noise. This is obtained by placing the obstacles randomly, but with an imposed minimal distance. Note the nice standing wave patterns in the wave height part, and the filamentation in the energy representation, similar to a branched flow, illustrated for instance in the video https://youtu.be/kOR6OiRlL8Y This video has two parts, showing the same evolution with two different color gradients: Wave height: 0:00 Averaged wave energy: 1:27 In the first part, the color hue depends on the height of the wave. In the second part, it depends on the energy of the wave, averaged over a sliding time window. The contrast has been enhanced by a shading procedure, similar to the one I have used on videos of reaction-diffusion equations. The process is to compute the normal vector to a surface in 3D that would be obtained by using the third dimension to represent the field, and then to make the luminosity depend on the angle between the normal vector and a fixed direction. There are absorbing boundary conditions on the borders of the simulated rectangle. The display at the right shows a time-averaged version of the signal near the right boundary of the simulated rectangular area. More precisely, it shows the square root of an average of squares of the respective field value (wave height or energy). Render time: 37 minutes 6 seconds Compression: crf 23 Color scheme: Part 1 - Parula, originally from Matlab https://www.mathworks.com/help/matlab/ref/colormap.html Part 2 - Turbo, by Anton Mikhailov https://gist.github.com/mikhailov-work/6a308c20e494d9e0ccc29036b28faa7a Music: Desert Catharsis by Asher Fulero@AsherFulero See also https://images.math.cnrs.fr/des-ondes-dans-mon-billard-partie-i/ for more explanations (in French) on a few previous simulations of wave equations. The simulation solves the wave equation by discretization. The algorithm is adapted from the paper https://hplgit.github.io/fdm-book/doc/pub/wave/pdf/wave-4print.pdf C code: https://github.com/nilsberglund-orleans/YouTube-simulations https://www.idpoisson.fr/berglund/software.html Many thanks to Marco Mancini and Julian Kauth for helping me to accelerate my code! #wave #diffraction$

Mega-tetris: Coagulating falling squares with alignment

Youngs double-slit experiment for a quantum particle

$Viewing two light sources through a gradient index lens Several recent videos on this channel, such as https://youtu.be/CdpTq3dSrXA and https://youtu.be/aYvsBgyJgIE , showed planar waves crossing a gradient index lens. In this video, we look instead at two light sources that are located quite close to the lens. The lens is indeed able to focus the light on two spots located in the focal plane. In practice, one seldom has light sources located so close to the lens (except perhaps in macro photography), but this set-up was better suited to the simulation method used here. To be able to get such a small focal length, the refractive index has been chosen rather large in the center of the lens, with a value of 1.82. The videos are inspired by Huygens Optics recent short https://youtube.com/shorts/VGd3Ajnp6e0 showing the principle of a gradient index lens. Lenses focus incoming rays of light by delaying them more near the center of the lens than at its periphery. This is often done with a material of constant index of refraction, by making the lens thicker near the center, as shown for instance in the simulation https://youtu.be/rrJJBh9ubUE . However, one can also build lenses of constant thickness, by making the index of refraction of their material depend on the location in the lens. In this simulation, the index decreases like sqrt(n0² - a*r²), where r is the distance to the axis of symmetry. This results in the incoming waves being focused at two points in the (estimated) focal plane, marked by a vertical line. The plot to the right shows a time-averaged value of the field along that plane. This video has two parts, showing the same evolution with two different color gradients: Wave height: 0:00 Averaged wave energy: 1:12 In the first part, the color hue depends on the height of the wave. In the second part, it depends on the energy of the wave, averaged over a sliding time window. There are absorbing boundary conditions on the borders of the simulated rectangle. The display at the right shows the signal along the focal plane, which is indicated by a vertical line. Render time: 38 minutes 12 seconds Compression: crf 23 Color scheme: Part 1 - Twilight by Bastian Bechtold https://github.com/bastibe/twilight Part 2 - Inferno by Nathaniel J. Smith and Stefan van der Walt https://github.com/BIDS/colormap Music: Jupiter One by Riot See also https://images.math.cnrs.fr/Des-ondes-dans-mon-billard-partie-I.html for more explanations (in French) on a few previous simulations of wave equations. The simulation solves the wave equation by discretization. The algorithm is adapted from the paper https://hplgit.github.io/fdm-book/doc/pub/wave/pdf/wave-4print.pdf C code: https://github.com/nilsberglund-orleans/YouTube-simulations https://www.idpoisson.fr/berglund/software.html Many thanks to Marco Mancini and Julian Kauth for helping me to accelerate my code! #wave #lens #gradient_index$

Growing quasicrystals from pentagons: Orientation and number of neighbors

Exciting resonant modes in a circle @NilsBerglund