A molecular sieve: sorting particles by size Nils Berglund 2024-09-01 | In this simulation, I tried to come up with a system that is able to sort particles according to their size. It consists in three grids of obstacles with decreasing space between them, acting as particle sieves. The obstacles rotate and exert a tangential force on the particles, in order to decrease clogging of the sieves.The conveyor belt effect results from the segments forming the belt exerting a tangential force on the polygons, in addition to the normal force. The tangential force is proportional to the difference between the tangential speed of the polygon and the speed of the belt. To compute the force and torque of polygon j on polygon i, the code computes the distance of each vertex of polygon j to the faces of polygon i. If this distance is smaller than a threshold, the force increases linearly with a large spring constant. In addition, radial forces between the vertices of the polygons have been added, whenever a vertex of polygon j is not on a perpendicular to a face of polygon i. This is important, because otherwise triangles can approach each other from the vertices, and when one vertex moves sideways, it is suddenly strongly accelerated, causing numerical instability. A weak Lennard-Jones interaction between polygons has been added, as it seems to increase numerical stability. Unlike in some previous videos involving interacting polygons, there is no thermostat in this simulation. Instead, friction forces (both linear and angular) have been added for numerical stability. In addition, the particles are subject to a gravitational force directed downwards. The color of the polygons depends on their size.To save on computation time, particles are placed into a "hash grid", each cell of which contains between 3 and 10 particles. Then only the influence of other particles in the same or neighboring cells is taken into account for each particle. The Lennard-Jones potential is strongly repulsive at short distance, and mildly attracting at long distance. It is widely used as a simple yet realistic model for the motion of electrically neutral molecules. The force results from the repulsion between electrons due to Pauli's exclusion principle, while the attractive part is a more subtle effect appearing in a multipole expansion. For more details, see en.wikipedia.org/wiki/Lennard-Jones_potentialRender time: 50 minutes 6 secondsCompression: crf 23Color scheme: Turbo, by Anton Mikhailovgist.github.com/mikhailov-work/6a308c20e494d9e0ccc29036b28faa7a Music: "Crystal" by Vibe TracksCurrent version of the C code used to make these animations: github.com/nilsberglund-orleans/YouTube-simulationshttps://www.idpoisson.fr/berglund/software.htmlSome outreach articles on mathematics:https://images.math.cnrs.fr/auteurs/nils-berglund/(in French, some with a Spanish translation)#molecular_dynamics #polygons #conveyor
Trajectories in a mass spectrometer Nils Berglund 2024-10-18 | This remake of the video youtu.be/O837373mq-Y with particle trajectories was requested by viewers. The simulation illustrates the principle of a mass spectrometer, which is a device allowing to sort charged particles according the their mass/charge ratio.The particles injected on the left side have an random radius, a mass proportional to their radius squared, and a constant positive charge. There is a constant electric field, directed from left to right, that accelerates the particles. In the middle of the simulation region, there is a constant magnetic field, perpendicular to the simulation plane. The particles interact via a Lennard-Jones potential, and are subject to a viscous drag, in order to avoid numerical instability. The Coulomb interaction between particles has been turned off, however, since otherwise they would not stay in the bins at the right (in practice, one uses more elaborate ways to collect the separated particles). The particles' acceleration, resulting from the Lorentz force, is proportional to the ratio between their charge and their mass. Therefore, lighter particles are deflected more, and tend to land in lower bins than heavier particles. The video has two parts, showing the same simulation with two different representations.With trajectories: 0:00Without trajectories: 2:44The particles' color depends on their size. To save on computation time, particles are placed into a "hash grid", each cell of which contains between 3 and 10 particles. Then only the influence of other particles in the same or neighboring cells is taken into account for each particle. The Lennard-Jones potential is strongly repulsive at short distance, and mildly attracting at long distance. It is widely used as a simple yet realistic model for the motion of electrically neutral molecules. The force results from the repulsion between electrons due to Pauli's exclusion principle, while the attractive part is a more subtle effect appearing in a multipole expansion. For more details, see en.wikipedia.org/wiki/Lennard-Jones_potentialRender time: 3 minutes 7 seconds with tracers2 minutes 15 seconds without tracersCompression: crf 23Color scheme: Turbo, by Anton Mikhailovgist.github.com/mikhailov-work/6a308c20e494d9e0ccc29036b28faa7a Music: "Dusty Finger" by DJ Williams@djwmusicCurrent version of the C code used to make these animations: github.com/nilsberglund-orleans/YouTube-simulationshttps://www.idpoisson.fr/berglund/software.htmlSome outreach articles on mathematics:https://images.math.cnrs.fr/auteurs/nils-berglund/(in French, some with a Spanish translation)#molecular_dynamics #mass_spectrometer
Creating rotating waves in a circle with fifteen angled secondary cavities Nils Berglund 2024-10-17 | This simulation of a resonating cavity implements a suggestion made by viewers of a previous simulation. It is a variant of the simulation youtu.be/cimKdtw04hg using 15 instead of 10 secondary cavities. By placing the channels connecting the main cavity to the secondary ones at an angle, instead of radially, the waves excited in the main cavity should have a tendency to rotate. It turns out that this works indeed quite well. Another difference with the previous video is the order in which the sources fire: instead of oscillating in the same order they are placed around the central cavity, the sources at odd positions fire before those at even positions. This different phase shift creates a larger angular velocity.This video has two parts, showing the same evolution with two different color gradients:Wave height: 0:00Averaged wave energy: 1:43In 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: 20 minutes 26 secondsCompression: crf 23Color scheme: Part 1 - Viridis by Nathaniel J. Smith, Stefan van der Walt and Eric Firinggithub.com/BIDS/colormapPart 2 - Turbo, by Anton Mikhailovgist.github.com/mikhailov-work/6a308c20e494d9e0ccc29036b28faa7a Music: "Electrosamba" by Quincas Moreira@QuincasMoreiraSee 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-simulationshttps://www.idpoisson.fr/berglund/software.htmlMany thanks to Marco Mancini and Julian Kauth for helping me to accelerate my code!#wave #resonator
Simulation of an asteroid impact between Ireland and Newfoundland Nils Berglund 2024-10-16 | This new simulation of an asteroid impact assumes the asteroid falls into the North Atlantic ocean about halfway between Ireland and Newfoundland, more to the north than in the previous simulation youtu.be/cr3rCxjWL4k of a North Atlantic impact. The simulation includes a respawning of tracer particles, which gives a better idea of the velocity field. The simulation is based on a shallow water equation. The model includes lunar forcing accounting for tides, and allows for some flooding of land areas. The vertical scale of both mountains and waves has been exaggerated. The main difference between the nonlinear equation and the linear wave equation is that the wave speed becomes larger when the water is shallower, which can lead to build-up of the waves in coastal and other shallow regions. This causes problems in the simulation, because land masses can lead to blow-up of the solution. This problem is circumvented here by replacing the continents by a repelling force field, plus a dissipative term. The initial state features velocities radiating outward from the impact point of the asteroid. The shallow water equations are nonlinear equations, which give a better description of the motion of water than the linear wave equation. In particular, unlike the linear wave equation, they conserve the total volume of water. The linear equation gives an approximation of the solutions, when the wave height remains close to its average over space. The equations used here include viscosity and dissipation, as described for instance in en.wikipedia.org/wiki/Shallow_water_equations#Non-conservative_form , including the Coriolis force. One difficulty is to model the wetting boundary, which separates regions that are under water and those which are not. This difficulty has been circumvented here by replacing the continents by a repulsive force field, directed downslope, instead of a sharp boundary. The video has four parts, showing simulations at two different speeds and with two different visualizations:Time lapse, 3D: 0:00Time lapse, 2D: 0:16Original speed, 3D: 0:32Original speed, 2D: 1:38The color hue and radial coordinate show the height of the water, on an exaggerated radial scale. The 2D parts use a projection in equirectangular coordinates. In the 3D parts, the point of view is slowly rotating around the Earth in a plane containing its center. In parts 1 and 2, the animation has been speeded up by a factor 4.The velocity field is materialized by 2000 tracer particles that are advected by the flow. The tracers are randomly "respawned", by moving them to a random new location at random times. This is to avoid that particles concentrate in some areas, as has happened on some previous simulations. Render time: 3D parts - 2 hours 26 minutes2D parts - 2 hours 18 minutesColor scheme: Viridis by Nathaniel J. Smith, Stefan van der Walt and Eric Firinggithub.com/BIDS/colormapMusic: "The Arid Land" by Sir Cubworth@SirCubworthSee 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 2D shallow water equation by discretization (finite differences).C code: github.com/nilsberglund-orleans/YouTube-simulationshttps://www.idpoisson.fr/berglund/software.htmlMany thanks to Marco Mancini and Julian Kauth for helping me to accelerate my code!#shallowwater #waves #Earth
Winnowing 3: Throwing the particles from the ground Nils Berglund 2024-10-15 | This variant of the video youtu.be/jd5dmbsgqic shows again the process of sorting particles by size using wind, the difference being that here the particles are thrown from near the ground instead of being dropped from above. This is closer to the real process used for winnowing, or separating grain from chaff by throwing it into the air, after it has been threshed. There are a few more sorting errors than in the previous simulation, which seem to be due to collisions between particles. In this simulation, the particles have different sizes, but the same mass. They are subjected to gravity, and to a horizontal force proportional to their radius, modeling the effect of wind. The wind stops a short distance above the bins at the bottom, to avoid that the particles bounce on the bins' walls.The particles' color depends on their size. To save on computation time, particles are placed into a "hash grid", each cell of which contains between 3 and 10 particles. Then only the influence of other particles in the same or neighboring cells is taken into account for each particle. The Lennard-Jones potential is strongly repulsive at short distance, and mildly attracting at long distance. It is widely used as a simple yet realistic model for the motion of electrically neutral molecules. The force results from the repulsion between electrons due to Pauli's exclusion principle, while the attractive part is a more subtle effect appearing in a multipole expansion. For more details, see en.wikipedia.org/wiki/Lennard-Jones_potentialRender time: 4 minutes 12 secondsCompression: crf 23Color scheme: Turbo, by Anton Mikhailovgist.github.com/mikhailov-work/6a308c20e494d9e0ccc29036b28faa7a Music: "Two Sides" by SefChol@SefCholCurrent version of the C code used to make these animations: github.com/nilsberglund-orleans/YouTube-simulationshttps://www.idpoisson.fr/berglund/software.htmlSome outreach articles on mathematics:https://images.math.cnrs.fr/auteurs/nils-berglund/(in French, some with a Spanish translation)#molecular_dynamics #winnowing
Creating rotating waves in a circle with ten angled secondary cavities Nils Berglund 2024-10-14 | This simulation of a resonating cavity implements a suggestion made by viewers of a previous simulation: By placing the channels connecting the main cavity to the secondary ones at an angle, instead of radially, the waves excited in the main cavity should have a tendency to rotate. It turns out that this works indeed quite well.This video has two parts, showing the same evolution with two different color gradients:Wave height: 0:00Averaged wave energy: 2:04In 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: 22 minutes 52 secondsCompression: crf 23Color scheme: Part 1 - Cividis by Jamie R. Nuñez, Christopher R. Anderton, Ryan S. Renslowjournals.plos.org/plosone/article?id=10.1371/journal.pone.0199239Part 2 - Inferno by Nathaniel J. Smith and Stefan van der Waltgithub.com/BIDS/colormapMusic: "September Pass" by Asher Fulero@AsherFuleroSee 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-simulationshttps://www.idpoisson.fr/berglund/software.htmlMany thanks to Marco Mancini and Julian Kauth for helping me to accelerate my code!#wave #resonator
Viewer request: Simulation of an asteroid impact off the coast of California Nils Berglund 2024-10-13 | This simulation has been requested by a viewer, who wanted to see the effect of an asteroid impact "500 km off the coast of Santa Monica". They didn't say why, maybe it's because their mother in law lives there, but since I have the code ready, I am happy to oblige. The simulation is based on a shallow water equation. The model includes lunar forcing accounting for tides, and allows for some flooding of land areas. The vertical scale of both mountains and waves has been exaggerated. The main difference between the nonlinear equation and the linear wave equation is that the wave speed becomes larger when the water is shallower, which can lead to build-up of the waves in coastal and other shallow regions. This causes problems in the simulation, because land masses can lead to blow-up of the solution. This problem is circumvented here by replacing the continents by a repelling force field, plus a dissipative term. The initial state features velocities radiating outward from the impact point of the asteroid. The shallow water equations are nonlinear equations, which give a better description of the motion of water than the linear wave equation. In particular, unlike the linear wave equation, they conserve the total volume of water. The linear equation gives an approximation of the solutions, when the wave height remains close to its average over space. The equations used here include viscosity and dissipation, as described for instance in en.wikipedia.org/wiki/Shallow_water_equations#Non-conservative_form , including the Coriolis force. One difficulty is to model the wetting boundary, which separates regions that are under water and those which are not. This difficulty has been circumvented here by replacing the continents by a repulsive force field, directed downslope, instead of a sharp boundary. The video has four parts, showing simulations at two different speeds and with two different visualizations:Time lapse, 3D: 0:00Time lapse, 2D: 0:14Original speed, 3D: 0:29Original speed, 2D: 1:26The color hue and radial coordinate show the height of the water, on an exaggerated radial scale. The 2D parts use a projection in equirectangular coordinates. In the 3D parts, the point of view is slowly rotating around the Earth in a plane containing its center. In parts 1 and 2, the animation has been speeded up by a factor 4.The velocity field is materialized by 2000 tracer particles that are advected by the flow. The tracers are randomly "respawned", by moving them to a random new location at random times. This is to avoid that particles concentrate in some areas, as has happened on some previous simulations. Render time: 3D parts - 2 hours 2 minutes2D parts - 1 hour 50 minutesColor scheme: Viridis by Nathaniel J. Smith, Stefan van der Walt and Eric Firinggithub.com/BIDS/colormapMusic: "Frequency" by Silent PartnerSee 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 2D shallow water equation by discretization (finite differences).C code: github.com/nilsberglund-orleans/YouTube-simulationshttps://www.idpoisson.fr/berglund/software.htmlMany thanks to Marco Mancini and Julian Kauth for helping me to accelerate my code!#shallowwater #waves #Earth
Winnowing 2: Separating triangles by size using wind Nils Berglund 2024-10-12 | This is a variant of the video youtu.be/jd5dmbsgqic , illustrating the process of winnowing, but with triangular instead of circular particles. Using wind to separate particles according to their size is an old method, going back to winnowing, or separating grain from chaff by throwing it into the air after it has been threshed. In this simulation, the particles have different sizes, but the same mass. They are subjected to gravity, and to a horizontal force proportional to their radius, modeling the effect of wind. The wind stops a short distance above the bins at the bottom, to avoid that the particles bounce on the bins' walls.The video has two parts, showing the same simulation with two different color gradients.Particle size: 0:00Kinetic energy: 2:36In part 1, the particles' color depends on their size. In part 2, in depends on the direction of their velocity. To save on computation time, particles are placed into a "hash grid", each cell of which contains between 3 and 10 particles. Then only the influence of other particles in the same or neighboring cells is taken into account for each particle. Render time: 4 minutes 12 secondsCompression: crf 23Color scheme: Turbo, by Anton Mikhailovgist.github.com/mikhailov-work/6a308c20e494d9e0ccc29036b28faa7a Music: "Sunshine on Sand" by the Unicorn Heads@UnicornHeadsCurrent version of the C code used to make these animations: github.com/nilsberglund-orleans/YouTube-simulationshttps://www.idpoisson.fr/berglund/software.htmlSome outreach articles on mathematics:https://images.math.cnrs.fr/auteurs/nils-berglund/(in French, some with a Spanish translation)#molecular_dynamics #winnowing
Exciting resonant modes in a circle, with a single source and sixteen secondary cavities Nils Berglund 2024-10-11 | This is a variant of the video youtu.be/XmJ_Ydldqvk , showing waves in a circular cavity excited by pulsing sources in other, connected cavities. In this case, there are sixteen secondary cavities, but only one source. One can notice that after a while, this source excites standing waves in all sixteen cavities. This video has two parts, showing the same evolution with two different color gradients:Wave height: 0:00Averaged wave energy: 2:12In 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: 25 minutes 54 secondsCompression: crf 23Color scheme: Part 1 - Cividis by Jamie R. Nuñez, Christopher R. Anderton, Ryan S. Renslowjournals.plos.org/plosone/article?id=10.1371/journal.pone.0199239Part 2 - Inferno by Nathaniel J. Smith and Stefan van der Waltgithub.com/BIDS/colormapMusic: "Lazy Laura" by Quincas Moreira@QuincasMoreiraSee 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-simulationshttps://www.idpoisson.fr/berglund/software.htmlMany thanks to Marco Mancini and Julian Kauth for helping me to accelerate my code!#wave #resonator
A tide simulation with a more realistic lunar forcing Nils Berglund 2024-10-10 | The simulation youtu.be/QYqjmapPKsc tried to reproduce tides on the Earth, but the lunar forcing was way too strong. In this variant, the lunar forcing has been reduced to about 40% of its value in the previous video, with a more realistic result.The effect of the Moon is modeled here by a force acting on the water height, which is maximal at the part of Earth closest to the Moon and at its antipodal point. The motion of the water is modeled by a shallow water equation, taking the varying ocean depth into account.The main difference between the nonlinear equation and the linear wave equation is that the wave speed becomes larger when the water is shallower, which can lead to build-up of the waves in coastal and other shallow regions. This causes problems in the simulation, because land masses can lead to blow-up of the solution. This problem is circumvented here by replacing the continents by a repelling force field, plus a dissipative term. The shallow water equations are nonlinear equations, which give a better description of the motion of water than the linear wave equation. In particular, unlike the linear wave equation, they conserve the total volume of water. The linear equation gives an approximation of the solutions, when the wave height remains close to its average over space. The equations used here include viscosity and dissipation, as described for instance in en.wikipedia.org/wiki/Shallow_water_equations#Non-conservative_form , including the Coriolis force. One difficulty is to model the wetting boundary, which separates regions that are under water and those which are not. This difficulty has been circumvented here by replacing the continents by a repulsive force field, directed downslope, instead of a sharp boundary. The video has four parts, showing simulations at two different speeds and with two different visualizations:Original speed, 3D: 0:00Original speed, 2D: 1:13Time lapse, 3D: 2:27Time lapse, 2D: 2:46The color hue and radial coordinate show the height of the water, on an exaggerated radial scale. The 2D parts use a projection in equirectangular coordinates. The longitude closest to the Moon is indicated by a vertical line in the 2D parts. In the 3D parts, the point of view is slowly rotating around the Earth in a circular orbit. In parts 3 and 4, the animation has been speeded up by a factor 4.The velocity field is materialized by 2000 tracer particles that are advected by the flow. A modification compared to previous simulations is that the particles are "respawned" at random times, meaning that their location is changed randomly. This avoids the concentration of tracers seen in some previous simulations.Render time: 3D parts - 2 hours 12 minutes2D parts - 2 hour 46 minutesColor scheme: Viridis by Nathaniel J. Smith, Stefan van der Walt and Eric Firinggithub.com/BIDS/colormapMusic: "Boom Bap Flick" by Quincas Moreira@QuincasMoreiraSee 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 2D shallow water equation by discretization (finite differences).C code: github.com/nilsberglund-orleans/YouTube-simulationshttps://www.idpoisson.fr/berglund/software.htmlMany thanks to Marco Mancini and Julian Kauth for helping me to accelerate my code!#shallowwater #waves #Earth #tsunami
Winnowing: Separating particles by size using wind Nils Berglund 2024-10-09 | Using wind to separate particles according to their size is an old method, going back to winnowing, or separating grain from chaff by throwing it into the air after it has been threshed. In this simulation, the particles have different sizes, but the same mass. They are subjected to gravity, and to a horizontal force proportional to their radius, modeling the effect of wind. The wind stops a short distance above the bins at the bottom, to avoid that the particles bounce on the bins' walls.The video has two parts, showing the same simulation with two different color gradients.Particle size: 0:00Kinetic energy: 2:24In part 1, the particles' color depends on their size. In part 2, in depends on the direction of their velocity. To save on computation time, particles are placed into a "hash grid", each cell of which contains between 3 and 10 particles. Then only the influence of other particles in the same or neighboring cells is taken into account for each particle. The Lennard-Jones potential is strongly repulsive at short distance, and mildly attracting at long distance. It is widely used as a simple yet realistic model for the motion of electrically neutral molecules. The force results from the repulsion between electrons due to Pauli's exclusion principle, while the attractive part is a more subtle effect appearing in a multipole expansion. For more details, see en.wikipedia.org/wiki/Lennard-Jones_potentialRender time: 4 minutes 12 secondsCompression: crf 23Color scheme: Turbo, by Anton Mikhailovgist.github.com/mikhailov-work/6a308c20e494d9e0ccc29036b28faa7a Music: John Stockton Slow Drag by Chris Zabriskie is licensed under a Creative Commons Attribution 4.0 licence. creativecommons.org/licenses/by/4.0/Source: http://chriszabriskie.com/uvp/Artist: http://chriszabriskie.comCurrent version of the C code used to make these animations: github.com/nilsberglund-orleans/YouTube-simulationshttps://www.idpoisson.fr/berglund/software.htmlSome outreach articles on mathematics:https://images.math.cnrs.fr/auteurs/nils-berglund/(in French, some with a Spanish translation)#molecular_dynamics #winnowing
A more symmetric version of resonances in a circle excited by five out of phase sources Nils Berglund 2024-10-08 | This is a variant of the video youtu.be/4Qjj9gTLm4k showing waves excited by five out of phase sources in a circular cavity. As was pointed out in the comments, the result lacked five-fold symmetry, which seemed to be due to a phase shift of a tenth, rather than a fifth full turn. Indeed the additional factor two stemmed from the fact that the excitations changed sign after each period. This has been corrected here, resulting in somewhat more symmetric patterns. One can still notice some imperfections, though, that probably stem from the simulation lattice.This video has two parts, showing the same evolution with two different color gradients:Wave height: 0:00Averaged wave energy: 1:32In 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: 15 minutes 13 secondsCompression: crf 23Color scheme: Part 1 - Cividis by Jamie R. Nuñez, Christopher R. Anderton, Ryan S. Renslowjournals.plos.org/plosone/article?id=10.1371/journal.pone.0199239Part 2 - Twilight by Bastian Bechtoldgithub.com/bastibe/twilightMusic: Spooky Ride by Twin Musicom is licensed under a Creative Commons Attribution 4.0 licence. creativecommons.org/licenses/by/4.0/Source: http://www.twinmusicom.org/song/250/spooky-rideArtist: http://www.twinmusicom.orgSee 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-simulationshttps://www.idpoisson.fr/berglund/software.htmlMany thanks to Marco Mancini and Julian Kauth for helping me to accelerate my code!#wave #resonator
How tides could look if the Moon were much closer to Earth Nils Berglund 2024-10-07 | Like in the video youtu.be/DNPFUomUFyw , this simulation uses a shallow water equation as a simplified model for the tides. The main difference with the previous simulation is the initial state, which already takes the influence of the Moon into account.The effect of the Moon is modeled by a force acting on the water height, which is maximal at the position of the Moon and at its antipode. The parameter tuning this effect has been chosen too large, however, so that in effect we get an impression of what tides would look like if the Moon were much closer to the Earth than it actually is. The main difference between the nonlinear equation and the linear wave equation is that the wave speed becomes larger when the water is shallower, which can lead to build-up of the waves in coastal and other shallow regions. This causes problems in the simulation, because land masses can lead to blow-up of the solution. This problem is circumvented here by replacing the continents by a repelling force field, plus a dissipative term. The shallow water equations are nonlinear equations, which give a better description of the motion of water than the linear wave equation. In particular, unlike the linear wave equation, they conserve the total volume of water. The linear equation gives an approximation of the solutions, when the wave height remains close to its average over space. The equations used here include viscosity and dissipation, as described for instance in en.wikipedia.org/wiki/Shallow_water_equations#Non-conservative_form , including the Coriolis force. One difficulty is to model the wetting boundary, which separates regions that are under water and those which are not. This difficulty has been circumvented here by replacing the continents by a repulsive force field, directed downslope, instead of a sharp boundary. The video has four parts, showing simulations at two different speeds and with two different visualizations:Original speed, 3D: 0:00Original speed, 2D: 1:05Time lapse, 3D: 2:12Time lapse, 2D: 2:28The color hue and radial coordinate show the height of the water, on an exaggerated radial scale. The 2D parts use a projection in equirectangular coordinates. The longitude closest to the Moon is indicated by a vertical line in the 2D parts. In the 3D parts, the point of view is slowly rotating around the Earth in a circular orbit. In parts 3 and 4, the animation has been speeded up by a factor 4.The velocity field is materialized by 2000 tracer particles that are advected by the flow. Render time: 3D parts - 2 hours 29 minutes2D parts - 1 hour 51 minutesColor scheme: Viridis by Nathaniel J. Smith, Stefan van der Walt and Eric Firinggithub.com/BIDS/colormapMusic: "You're Not Wrong" by roljui@roljui1445See 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 2D shallow water equation by discretization (finite differences).C code: github.com/nilsberglund-orleans/YouTube-simulationshttps://www.idpoisson.fr/berglund/software.htmlMany thanks to Marco Mancini and Julian Kauth for helping me to accelerate my code!#shallowwater #waves #Earth #tsunami
A particle sorter with rollers of constant radius and increasing gaps Nils Berglund 2024-10-06 | This simulation uses a similar device sorting particles by size as the videos youtu.be/HNBTewymZa4 and youtu.be/U13YpO-nepw . The main difference is that instead of using a sieve with rollers having equally spaced centers, and decreasing radius, here the rollers have constant radius, while the distance between their axes increases. This modification was suggested in a comment to a previous video.The obstacles rotate and exert a tangential force on the particles, in order to decrease clogging of the sieves. The angular speed of the obstacles varies periodically in time, in order to reduce the chance of particles getting stuck. In addition, the obstacles have a circular motion of small amplitude, which reduces the chances of particles getting stuck even more. The conveyor belt effect results from the segments forming the belt exerting a tangential force on the polygons, in addition to the normal force. The tangential force is proportional to the difference between the tangential speed of the polygon and the speed of the belt. To compute the force and torque of polygon j on polygon i, the code computes the distance of each vertex of polygon j to the faces of polygon i. If this distance is smaller than a threshold, the force increases linearly with a large spring constant. In addition, radial forces between the vertices of the polygons have been added, whenever a vertex of polygon j is not on a perpendicular to a face of polygon i. This is important, because otherwise triangles can approach each other from the vertices, and when one vertex moves sideways, it is suddenly strongly accelerated, causing numerical instability. A weak Lennard-Jones interaction between polygons has been added, as it seems to increase numerical stability. Unlike in some previous videos involving interacting polygons, there is no thermostat in this simulation. Instead, friction forces (both linear and angular) have been added for numerical stability. In addition, the particles are subject to a gravitational force directed downwards. The color of the polygons depends on their size.To save on computation time, particles are placed into a "hash grid", each cell of which contains between 3 and 10 particles. Then only the influence of other particles in the same or neighboring cells is taken into account for each particle. The Lennard-Jones potential is strongly repulsive at short distance, and mildly attracting at long distance. It is widely used as a simple yet realistic model for the motion of electrically neutral molecules. The force results from the repulsion between electrons due to Pauli's exclusion principle, while the attractive part is a more subtle effect appearing in a multipole expansion. For more details, see en.wikipedia.org/wiki/Lennard-Jones_potentialRender time: 1 hour 15 minutesCompression: crf 23Color scheme: Turbo, by Anton Mikhailovgist.github.com/mikhailov-work/6a308c20e494d9e0ccc29036b28faa7a Music: Chee Zee Beach - Latinesque by Kevin MacLeod is licensed under a Creative Commons Attribution 4.0 licence. creativecommons.org/licenses/by/4.0Source: http://incompetech.com/music/royalty-free/index.html?isrc=USUAN1100686Artist: http://incompetech.comCurrent version of the C code used to make these animations: github.com/nilsberglund-orleans/YouTube-simulationshttps://www.idpoisson.fr/berglund/software.htmlSome outreach articles on mathematics:https://images.math.cnrs.fr/auteurs/nils-berglund/(in French, some with a Spanish translation)#molecular_dynamics #polygons #conveyor
Exciting resonant modes in a circle, with fewer sources than cavities Nils Berglund 2024-10-05 | This video is based on a suggestion in a comment. Like several previous videos on this channel , this one shows waves in a circular cavity excited by pulsing sources in other, connected cavities. In this case, there are sixteen cavities, but only six sources, pulsing with a phase shift of one twelfth of the period. Source are located in every third cavity, where the spacing of 3 has been chosen because it is coprime with 16. One can notice that after a while, these sources excite standing waves in all sixteen cavities. This video has two parts, showing the same evolution with two different color gradients:Wave height: 0:00Averaged wave energy: 1:40In 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: 18 minutes 53 secondsCompression: crf 23Color scheme: Part 1 - Viridis by Nathaniel J. Smith, Stefan van der Walt and Eric Firinggithub.com/BIDS/colormapPart 2 - Cividis by Jamie R. Nuñez, Christopher R. Anderton, Ryan S. Renslowjournals.plos.org/plosone/article?id=10.1371/journal.pone.0199239Music: Living Voyage by Kevin MacLeod is licensed under a Creative Commons Attribution 4.0 licence. creativecommons.org/licenses/by/4.0/Source: http://incompetech.com/music/royalty-free/index.html?isrc=USUAN1100594Artist: http://incompetech.comSee 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-simulationshttps://www.idpoisson.fr/berglund/software.htmlMany thanks to Marco Mancini and Julian Kauth for helping me to accelerate my code!#wave #resonator
Trying to model tides with a shallow water equation Nils Berglund 2024-10-04 | This is an attempt at modeling tides with a shallow water equation. A version using the linear wave equation previously appeared in the video youtu.be/SqwvZX20Q2Y . The initial state has been chosen somewhat arbitrarily, just to have some currents in the oceans. Then the effect of the Moon is modeled by a force acting on the water height, which is maximal at the position of the Moon and at its antipode. The longitude closest to the Moon is indicated by a vertical line in the 2D parts.The main difference between the nonlinear equation and the linear wave equation is that the wave speed becomes larger when the water is shallower, which can lead to build-up of the waves in coastal and other shallow regions. This causes problems in the simulation, because land masses can lead to blow-up of the solution. This problem is circumvented here by replacing the continents by a repelling force field, plus a dissipative term. The shallow water equations are nonlinear equations, which give a better description of the motion of water than the linear wave equation. In particular, unlike the linear wave equation, they conserve the total volume of water. The linear equation gives an approximation of the solutions, when the wave height remains close to its average over space. The equations used here include viscosity and dissipation, as described for instance in en.wikipedia.org/wiki/Shallow_water_equations#Non-conservative_form , including the Coriolis force. One difficulty is to model the wetting boundary, which separates regions that are under water and those which are not. This difficulty has been circumvented here by replacing the continents by a repulsive force field, directed downslope, instead of a sharp boundary. The video has four parts, showing simulations at two different speeds and with two different visualizations:Time lapse, 3D: 0:00Time lapse, 2D: 0:17Original speed, 3D: 0:35Original speed, 2D: 1:45The color hue and radial coordinate show the height of the water, on an exaggerated radial scale. The 2D parts use a projection in equirectangular coordinates. In the 3D parts, the point of view is slowly rotating around the Earth in a circular orbit. In parts 1 and 2, the animation has been speeded up by a factor 4.The velocity field is materialized by 2000 tracer particles that are advected by the flow. Render time: 3D parts - 1 hour 58 minutes2D parts - 2 hours 42 minutesColor scheme: Viridis by Nathaniel J. Smith, Stefan van der Walt and Eric Firinggithub.com/BIDS/colormapMusic: "5 Star" by Causmic@CausmicSee 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 2D shallow water equation by discretization (finite differences).C code: github.com/nilsberglund-orleans/YouTube-simulationshttps://www.idpoisson.fr/berglund/software.htmlMany thanks to Marco Mancini and Julian Kauth for helping me to accelerate my code!#shallowwater #waves #Earth #tsunami
A mass spectrometer Nils Berglund 2024-10-03 | This simulation illustrates the principle of a mass spectrometer, which is a device allowing to sort charged particles according the their mass/charge ratio.The particles injected on the left side have an random radius, a mass proportional to their radius squared, and a constant positive charge. There is a constant electric field, directed from left to right, that accelerates the particles. In the middle of the simulation region, there is a constant magnetic field, perpendicular to the simulation plane. The particles interact via a Lennard-Jones potential, and are subject to a viscous drag, in order to avoid numerical instability. The Coulomb interaction between particles has been turned off, however, since otherwise they would not stay in the bins at the right (in practice, one uses more elaborate ways to collect the separated particles). The particles' acceleration, resulting from the Lorentz force, is proportional to the ratio between their charge and their mass. Therefore, lighter particles are deflected more, and tend to land in lower bins than heavier particles. The video has two parts, showing the same simulation with two different color gradients.Particle size: 0:00Kinetic energy: 2:44In part 1, the particles' color depends on their size. In part 2, in depends on their kinetic energy, which is proportional to their mass times their speed squared. To save on computation time, particles are placed into a "hash grid", each cell of which contains between 3 and 10 particles. Then only the influence of other particles in the same or neighboring cells is taken into account for each particle. The Lennard-Jones potential is strongly repulsive at short distance, and mildly attracting at long distance. It is widely used as a simple yet realistic model for the motion of electrically neutral molecules. The force results from the repulsion between electrons due to Pauli's exclusion principle, while the attractive part is a more subtle effect appearing in a multipole expansion. For more details, see en.wikipedia.org/wiki/Lennard-Jones_potentialRender time: 4 minutes 15 secondsCompression: crf 23Color scheme: Turbo, by Anton Mikhailovgist.github.com/mikhailov-work/6a308c20e494d9e0ccc29036b28faa7a Music: "Trips" by Topher Mohr and Alex ElenaCurrent version of the C code used to make these animations: github.com/nilsberglund-orleans/YouTube-simulationshttps://www.idpoisson.fr/berglund/software.htmlSome outreach articles on mathematics:https://images.math.cnrs.fr/auteurs/nils-berglund/(in French, some with a Spanish translation)#molecular_dynamics #mass_spectrometer
Refraction and reflection of a shock wave - Source in medium with refractive index 3 Nils Berglund 2024-10-02 | Like the video youtu.be/196um7yJzf8 , this one shows a simulation of shock waves created by a source traveling faster than the wave speed in a medium, near an interface to a medium with different wave speed. The difference is that here, the lower medium has relative refractive index 0.333 (the wave speed is 3 times faster in the lower part). As a result, the waves in the lower medium sometimes propagate faster than the source. As before, part of the shock wave is refracted, changing its angle in the second medium, while part of it is reflected. The relative intensity of the reflected and refracted parts depends on the angle between shock wave and interface.This video has six parts, showing the evolution with three different speeds of the source, and two different color gradients for each speed:Speed 1, wave height: 0:00Speed 1, averaged wave energy: 0:52Speed 2, wave height: 1:52Speed 2, averaged wave energy: 2:44Speed 3, wave height: 3:44Speed 3, averaged wave energy: 4:36In parts 1, 3 and 5, the color hue depends on the height of the wave. In the other parts, 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.Render time: Part 1 - 21 minutes 26 secondsPart 2 - 20 minutes 26 secondsPart 3 - 21 minutes 10 secondsCompression: crf 23Color scheme: Part 1 - Viridis by Nathaniel J. Smith, Stefan van der Walt and Eric FiringPart 2 - Inferno by Nathaniel J. Smith and Stefan van der Waltgithub.com/BIDS/colormapMusic: Smooth Sailing (with Guitar) by Audionautix is licensed under a Creative Commons Attribution 4.0 licence. creativecommons.org/licenses/by/4.0/Artist: http://audionautix.comSee 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-simulationshttps://www.idpoisson.fr/berglund/software.htmlMany thanks to Marco Mancini and Julian Kauth for helping me to accelerate my code!#wave #shock_wave #sonic_boom
Modeling the 2023 tsunami due to a landslide and glacier collapse in eastern Greenland Nils Berglund 2024-10-01 | On 16th of September, 2023, seismologists around the world recorded an unusual signal, with an oscillation period of 90 seconds, and that lasted for 9 days. As explained in the video youtu.be/60T9TKuuujs , scientists found out that the origin of the signal was a landslide in uninhabited Dickson fjord in eastern Greenland, probably due to the collapse of a melting glacier. This landslide caused a tsunami with 200 meter high waves in the fjord. Fortunately, the wave height quickly decreased with distance from the landslide. The seismic signal was shown to have been caused by a standing wave in Dickson fjord, called a seiche, that took 9 days to dissipate.This very simplified simulation uses a shallow water equation to model the propagation of the tsunami through the Atlantic Ocean. Height of land masses and waves have been greatly exaggerated. In reality, unlike what this simulation might suggest, the tsunami had dissipated too much when reaching the coasts of Norway and Iceland to affect them. However, the seismic signal did travel much further. The main difference between the nonlinear equation and the linear wave equation is that the wave speed becomes larger when the water is shallower, which can lead to build-up of the waves in coastal and other shallow regions. This causes problems in the simulation, because land masses can lead to blow-up of the solution. This problem is circumvented here by replacing the continents by a repelling force field, plus a dissipative term. The initial state features velocities radiating outward from the impact point of the asteroid. The shallow water equations are nonlinear equations, which give a better description of the motion of water than the linear wave equation. In particular, unlike the linear wave equation, they conserve the total volume of water. The linear equation gives an approximation of the solutions, when the wave height remains close to its average over space. The equations used here include viscosity and dissipation, as described for instance in en.wikipedia.org/wiki/Shallow_water_equations#Non-conservative_form , including the Coriolis force. One difficulty is to model the wetting boundary, which separates regions that are under water and those which are not. This difficulty has been circumvented here by replacing the continents by a repulsive force field, directed downslope, instead of a sharp boundary. The video has four parts, showing simulations at two different speeds and with two different visualizations:Time lapse, 3D: 0:00Time lapse, 2D: 0:15Original speed, 3D: 0:29Original speed, 2D: 1:30The color hue and radial coordinate show the height of the water, on an exaggerated radial scale. The 2D parts use a projection in equirectangular coordinates. In the 3D parts, the point of view is slowly rotating around the Earth in a polar orbit. In parts 1 and 2, the animation has been speeded up by a factor 4.The velocity field is materialized by 2000 tracer particles that are advected by the flow. Thanks to Florentin Bieder from the University of Basel for letting me know about these findings. Render time: 3D parts - 2 hours 4 minutes2D parts - 2 hours 8 minutesColor scheme: Viridis by Nathaniel J. Smith, Stefan van der Walt and Eric Firinggithub.com/BIDS/colormapMusic: "Spring" by HOVATOFF@hovatoff7997See 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 2D shallow water equation by discretization (finite differences).C code: github.com/nilsberglund-orleans/YouTube-simulationshttps://www.idpoisson.fr/berglund/software.htmlMany thanks to Marco Mancini and Julian Kauth for helping me to accelerate my code!#shallowwater #waves #Earth #tsunami
Sorting heptagons with a linear sieve Nils Berglund 2024-09-30 | This is a variant of the video youtu.be/HNBTewymZa4 of a device sorting particles by size. The particles are heptagons instead of pentagons, and the friction is 20% lower as in the previous video. The particle sorter in this simulation is inspired by a comment to a previous video. Instead of using several sieves one above the other, it uses a single sieve, with rotating obstacles of decreasing radius and increasing gaps between them.The obstacles rotate and exert a tangential force on the particles, in order to decrease clogging of the sieves. The angular speed of the obstacles varies periodically in time, in order to reduce the chance of particles getting stuck. In addition, the obstacles have a circular motion of small amplitude, which reduces the chances of particles getting stuck even more. The conveyor belt effect results from the segments forming the belt exerting a tangential force on the polygons, in addition to the normal force. The tangential force is proportional to the difference between the tangential speed of the polygon and the speed of the belt. To compute the force and torque of polygon j on polygon i, the code computes the distance of each vertex of polygon j to the faces of polygon i. If this distance is smaller than a threshold, the force increases linearly with a large spring constant. In addition, radial forces between the vertices of the polygons have been added, whenever a vertex of polygon j is not on a perpendicular to a face of polygon i. This is important, because otherwise triangles can approach each other from the vertices, and when one vertex moves sideways, it is suddenly strongly accelerated, causing numerical instability. A weak Lennard-Jones interaction between polygons has been added, as it seems to increase numerical stability. Unlike in some previous videos involving interacting polygons, there is no thermostat in this simulation. Instead, friction forces (both linear and angular) have been added for numerical stability. In addition, the particles are subject to a gravitational force directed downwards. The color of the polygons depends on their size.To save on computation time, particles are placed into a "hash grid", each cell of which contains between 3 and 10 particles. Then only the influence of other particles in the same or neighboring cells is taken into account for each particle. The Lennard-Jones potential is strongly repulsive at short distance, and mildly attracting at long distance. It is widely used as a simple yet realistic model for the motion of electrically neutral molecules. The force results from the repulsion between electrons due to Pauli's exclusion principle, while the attractive part is a more subtle effect appearing in a multipole expansion. For more details, see en.wikipedia.org/wiki/Lennard-Jones_potentialRender time: 1 hour 30 minutesCompression: crf 23Color scheme: Turbo, by Anton Mikhailovgist.github.com/mikhailov-work/6a308c20e494d9e0ccc29036b28faa7a Music: "Digifunk" by DivKid@DivKidCurrent version of the C code used to make these animations: github.com/nilsberglund-orleans/YouTube-simulationshttps://www.idpoisson.fr/berglund/software.htmlSome outreach articles on mathematics:https://images.math.cnrs.fr/auteurs/nils-berglund/(in French, some with a Spanish translation)#molecular_dynamics #polygons #conveyor
Exciting spirals in a circle with sixteen out of phase sources Nils Berglund 2024-09-29 | Like several previous videos on this channel , this one shows waves in a circular cavity excited by pulsing sources in other, connected cavities. In this case, there are sixteen sources, pulsing a fraction 1/32 of a period out of phase. The resulting wave patterns start to resemble spirals, as suggested in comments to a previous video. This video has two parts, showing the same evolution with two different color gradients:Wave height: 0:00Averaged wave energy: 1:27In 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. The color dependence on wave height has been chosen quite important, so that the darker regions concentrate near the nodal lines of the waves, that is, the curves where the wave height is zero. Render time: 17 minutes 51 secondsCompression: crf 23Color scheme: Part 1 - Twilight by Bastian Bechtoldgithub.com/bastibe/twilightPart 2 - Inferno by Nathaniel J. Smith and Stefan van der Waltgithub.com/BIDS/colormapMusic: "Asleep" by Verified Picasso@verifiedpicasso7260See 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-simulationshttps://www.idpoisson.fr/berglund/software.htmlMany thanks to Marco Mancini and Julian Kauth for helping me to accelerate my code!#wave #resonator
An asteroid impact in the South Pacific Ocean Nils Berglund 2024-09-28 | This is an attempt at modeling an asteroid impact in the South Pacific Ocean using a nonlinear wave equation. This model allows for some flooding of land areas. The vertical scale of both mountains and waves has been exaggerated. The main difference between the nonlinear equation and the linear wave equation is that the wave speed becomes larger when the water is shallower, which can lead to build-up of the waves in coastal and other shallow regions. This causes problems in the simulation, because land masses can lead to blow-up of the solution. This problem is circumvented here by replacing the continents by a repelling force field, plus a dissipative term. The initial state features velocities radiating outward from the impact point of the asteroid. The shallow water equations are nonlinear equations, which give a better description of the motion of water than the linear wave equation. In particular, unlike the linear wave equation, they conserve the total volume of water. The linear equation gives an approximation of the solutions, when the wave height remains close to its average over space. The equations used here include viscosity and dissipation, as described for instance in en.wikipedia.org/wiki/Shallow_water_equations#Non-conservative_form , including the Coriolis force. One difficulty is to model the wetting boundary, which separates regions that are under water and those which are not. This difficulty has been circumvented here by replacing the continents by a repulsive force field, directed downslope, instead of a sharp boundary. The video has four parts, showing simulations at two different speeds and with two different visualizations:Time lapse, 3D: 0:00Time lapse, 2D: 0:17Original speed, 3D: 0:34Original speed, 2D: 1:45The color hue and radial coordinate show the height of the water, on an exaggerated radial scale. The 2D parts use a projection in equirectangular coordinates. In the 3D parts, the point of view is slowly rotating around the Earth in a plane containing its center. In parts 1 and 2, the animation has been speeded up by a factor 4.The velocity field is materialized by 2000 tracer particles that are advected by the flow. Render time: 3D parts - 2 hours 19 minutes2D parts - 2 hours 25 minutesColor scheme: Viridis by Nathaniel J. Smith, Stefan van der Walt and Eric Firinggithub.com/BIDS/colormapMusic: It Is Lost by Kevin MacLeod is licensed under a Creative Commons Attribution 4.0 licence. creativecommons.org/licenses/by/4.0Source: http://incompetech.com/music/royalty-free/index.html?isrc=USUAN1100190Artist: http://incompetech.comSee 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 2D shallow water equation by discretization (finite differences).C code: github.com/nilsberglund-orleans/YouTube-simulationshttps://www.idpoisson.fr/berglund/software.htmlMany thanks to Marco Mancini and Julian Kauth for helping me to accelerate my code!#shallowwater #waves #Earth
A particle sorter with linearly increasing gaps Nils Berglund 2024-09-27 | The particle sorter in this simulation is inspired by a comment to a previous video. Instead of using several sieves one above the other, it uses a single sieve, with rotating obstacles of decreasing radius and increasing gaps between them.The obstacles rotate and exert a tangential force on the particles, in order to decrease clogging of the sieves. The angular speed of the obstacles varies periodically in time, in order to reduce the chance of particles getting stuck. In addition, the obstacles have a circular motion of small amplitude, which reduces the chances of particles getting stuck even more. The conveyor belt effect results from the segments forming the belt exerting a tangential force on the polygons, in addition to the normal force. The tangential force is proportional to the difference between the tangential speed of the polygon and the speed of the belt. To compute the force and torque of polygon j on polygon i, the code computes the distance of each vertex of polygon j to the faces of polygon i. If this distance is smaller than a threshold, the force increases linearly with a large spring constant. In addition, radial forces between the vertices of the polygons have been added, whenever a vertex of polygon j is not on a perpendicular to a face of polygon i. This is important, because otherwise triangles can approach each other from the vertices, and when one vertex moves sideways, it is suddenly strongly accelerated, causing numerical instability. A weak Lennard-Jones interaction between polygons has been added, as it seems to increase numerical stability. Unlike in some previous videos involving interacting polygons, there is no thermostat in this simulation. Instead, friction forces (both linear and angular) have been added for numerical stability. In addition, the particles are subject to a gravitational force directed downwards. The color of the polygons depends on their size.To save on computation time, particles are placed into a "hash grid", each cell of which contains between 3 and 10 particles. Then only the influence of other particles in the same or neighboring cells is taken into account for each particle. The Lennard-Jones potential is strongly repulsive at short distance, and mildly attracting at long distance. It is widely used as a simple yet realistic model for the motion of electrically neutral molecules. The force results from the repulsion between electrons due to Pauli's exclusion principle, while the attractive part is a more subtle effect appearing in a multipole expansion. For more details, see en.wikipedia.org/wiki/Lennard-Jones_potentialRender time: 1 hour 9 minutesCompression: crf 23Color scheme: Turbo, by Anton Mikhailovgist.github.com/mikhailov-work/6a308c20e494d9e0ccc29036b28faa7a Music: "Fairy Meeting" by Emily A. Sprague@emilysprague6983Current version of the C code used to make these animations: github.com/nilsberglund-orleans/YouTube-simulationshttps://www.idpoisson.fr/berglund/software.htmlSome outreach articles on mathematics:https://images.math.cnrs.fr/auteurs/nils-berglund/(in French, some with a Spanish translation)#molecular_dynamics #polygons #conveyor
Exciting resonant modes in a circle with ten out of phase sources Nils Berglund 2024-09-26 | Like several previous videos on this channel, this one shows waves in a circular cavity excited by pulsing sources in other, connected cavities. In this case, there are ten sources, and unlike in previous simulations, the ten sources are pulsing a twentieth of a period out of phase. Some commenters suggested trying this, in the hope that it might lead to spiraling waves. This is not quite the case, but the resulting interference patterns are interesting nonetheless. This video has two parts, showing the same evolution with two different color gradients:Wave height: 0:00Averaged wave energy: 1:44In 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: 15 minutes 30 secondsCompression: crf 23Color scheme: Part 1 - Viridis by Nathaniel J. Smith, Stefan van der Walt and Eric FiringPart 2 - Magma by Nathaniel J. Smith and Stefan van der Waltgithub.com/BIDS/colormapMusic: "Absolutely Nothing" by Jeremy Blake@RedMeansRecordingSee 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-simulationshttps://www.idpoisson.fr/berglund/software.htmlMany thanks to Marco Mancini and Julian Kauth for helping me to accelerate my code!#wave #resonator
An equatorial asteroid impact in the Pacific Ocean Nils Berglund 2024-09-25 | This is an attempt at modeling an equatorial asteroid impact in the Pacific Ocean using a nonlinear wave equation. This model allows for some flooding of land areas. The vertical scale of both mountains and waves has been exaggerated. The main difference between the nonlinear equation and the linear wave equation is that the wave speed becomes larger when the water is shallower, which can lead to build-up of the waves in coastal and other shallow regions. This causes problems in the simulation, because land masses can lead to blow-up of the solution. This problem is circumvented here by replacing the continents by a repelling force field, plus a dissipative term. The initial state features velocities radiating outward from the impact point of the asteroid. The shallow water equations are nonlinear equations, which give a better description of the motion of water than the linear wave equation. In particular, unlike the linear wave equation, they conserve the total volume of water. The linear equation gives an approximation of the solutions, when the wave height remains close to its average over space. The equations used here include viscosity and dissipation, as described for instance in en.wikipedia.org/wiki/Shallow_water_equations#Non-conservative_form , including the Coriolis force. One difficulty is to model the wetting boundary, which separates regions that are under water and those which are not. This difficulty has been circumvented here by replacing the continents by a repulsive force field, directed downslope, instead of a sharp boundary. The video has four parts, showing simulations at two different speeds and with two different visualizations:Time lapse, 3D: 0:00Time lapse, 2D: 0:16Original speed, 3D: 0:33Original speed, 2D: 1:39The color hue and radial coordinate show the height of the water, on an exaggerated radial scale. The 2D parts use a projection in equirectangular coordinates. In the 3D parts, the point of view is slowly rotating around the Earth in a plane containing its center. In parts 1 and 2, the animation has been speeded up by a factor 4.The velocity field is materialized by 2000 tracer particles that are advected by the flow. Render time: 3D parts - 1 hour 49 minutes2D parts - 2 hours 17 minutesColor scheme: Viridis by Nathaniel J. Smith, Stefan van der Walt and Eric Firinggithub.com/BIDS/colormapMusic: "Red Sea" by RiotSee 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 2D shallow water equation by discretization (finite differences).C code: github.com/nilsberglund-orleans/YouTube-simulationshttps://www.idpoisson.fr/berglund/software.htmlMany thanks to Marco Mancini and Julian Kauth for helping me to accelerate my code!#shallowwater #waves #Earth
Sorting magnetized bestagons by size Nils Berglund 2024-09-24 | This simulation uses the same sorting device as in the video youtu.be/EcCHSkKQ8ic , including sieves that "rattle" to avoid clogging, with hexagonal instead of square particles. The friction force is proportional to the radius of the hexagons, which helps stabillizing the larger hexagons in the rightmost bin. In addition, a torque between hexagons has been added, which tends to align their orientation.The sorting set-up consists in three grids of obstacles with decreasing space between them, acting as particle sieves. The obstacles rotate and exert a tangential force on the particles, in order to decrease clogging of the sieves. The angular speed of the obstacles varies periodically in time, in order to reduce the chance of particles getting stuck. In addition, the obstacles have a circular motion of small amplitude, which reduces the chances of particles getting stuck even more. The conveyor belt effect results from the segments forming the belt exerting a tangential force on the polygons, in addition to the normal force. The tangential force is proportional to the difference between the tangential speed of the polygon and the speed of the belt. To compute the force and torque of polygon j on polygon i, the code computes the distance of each vertex of polygon j to the faces of polygon i. If this distance is smaller than a threshold, the force increases linearly with a large spring constant. In addition, radial forces between the vertices of the polygons have been added, whenever a vertex of polygon j is not on a perpendicular to a face of polygon i. This is important, because otherwise triangles can approach each other from the vertices, and when one vertex moves sideways, it is suddenly strongly accelerated, causing numerical instability. A weak Lennard-Jones interaction between polygons has been added, as it seems to increase numerical stability. Unlike in some previous videos involving interacting polygons, there is no thermostat in this simulation. Instead, friction forces (both linear and angular) have been added for numerical stability. In addition, the particles are subject to a gravitational force directed downwards. The color of the polygons depends on their size.To save on computation time, particles are placed into a "hash grid", each cell of which contains between 3 and 10 particles. Then only the influence of other particles in the same or neighboring cells is taken into account for each particle. The Lennard-Jones potential is strongly repulsive at short distance, and mildly attracting at long distance. It is widely used as a simple yet realistic model for the motion of electrically neutral molecules. The force results from the repulsion between electrons due to Pauli's exclusion principle, while the attractive part is a more subtle effect appearing in a multipole expansion. For more details, see en.wikipedia.org/wiki/Lennard-Jones_potentialRender time: 25 minutes 45 secondsCompression: crf 23Color scheme: Turbo, by Anton Mikhailovgist.github.com/mikhailov-work/6a308c20e494d9e0ccc29036b28faa7a Music: "Islabonita" by An Jone@SOYEMILIACurrent version of the C code used to make these animations: github.com/nilsberglund-orleans/YouTube-simulationshttps://www.idpoisson.fr/berglund/software.htmlSome outreach articles on mathematics:https://images.math.cnrs.fr/auteurs/nils-berglund/(in French, some with a Spanish translation)#molecular_dynamics #polygons #conveyor
Exciting resonant modes in a circle with five out of phase sources Nils Berglund 2024-09-23 | Like several previous videos on this channel , this one shows waves in a circular cavity excited by pulsing sources in other, connected cavities. In this case, there are five sources, and unlike in previous simulations, the five sources are pulsing a fifth of a period out of phase. Some commenters suggested trying this, in the hope that it might lead to spiraling waves. This is not quite the case, but the resulting interference patterns are interesting nonetheless. This video has two parts, showing the same evolution with two different color gradients:Wave height: 0:00Averaged wave energy: 1:32In 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: 15 minutes 30 secondsCompression: crf 23Color scheme: Part 1 - Cividis by Jamie R. Nuñez, Christopher R. Anderton, Ryan S. Renslowjournals.plos.org/plosone/article?id=10.1371/journal.pone.0199239Part 2 - Twilight by Bastian Bechtoldgithub.com/bastibe/twilightMusic: "Doorway" by Slenderbeats@slenderbodiesSee 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-simulationshttps://www.idpoisson.fr/berglund/software.htmlMany thanks to Marco Mancini and Julian Kauth for helping me to accelerate my code!#wave #resonator
An asteroid impact in the North Pacific Ocean, with tsunami-induced flooding Nils Berglund 2024-09-22 | This is an attempt at modeling an asteroid impact in the North Pacific Ocean using a nonlinear wave equation. This model allows for some flooding of land areas. The vertical scale of both mountains and waves has been exaggerated. A mistake in the code has been corrected, which incorrectly computed the influence of the non-constant ocean depth - the east-west component of that force was using the north-south gradient of the depth, instead of the east-west gradient. The display of land masses has also been improved, by using a different criterion than before to distinguish land from sea. The main difference between the nonlinear equation and the linear wave equation is that the wave speed becomes larger when the water is shallower, which can lead to build-up of the waves in coastal and other shallow regions. This causes problems in the simulation, because land masses can lead to blow-up of the solution. This problem is circumvented here by replacing the continents by a repelling force field, plus a dissipative term. The initial state features velocities radiating outward from the impact point of the asteroid. The shallow water equations are nonlinear equations, which give a better description of the motion of water than the linear wave equation. In particular, unlike the linear wave equation, they conserve the total volume of water. The linear equation gives an approximation of the solutions, when the wave height remains close to its average over space. The equations used here include viscosity and dissipation, as described for instance in en.wikipedia.org/wiki/Shallow_water_equations#Non-conservative_form , including the Coriolis force. One difficulty is to model the wetting boundary, which separates regions that are under water and those which are not. This difficulty has been circumvented here by replacing the continents by a repulsive force field, directed downslope, instead of a sharp boundary. The video has four parts, showing simulations at two different speeds and with two different visualizations:Time lapse, 3D: 0:00Time lapse, 2D: 0:15Original speed, 3D: 0:30Original speed, 2D: 1:33The color hue and radial coordinate show the height of the water, on an exaggerated radial scale. The 2D parts use a projection in equirectangular coordinates. In the 3D parts, the point of view is slowly rotating around the Earth in a plane containing its center. In parts 1 and 2, the animation has been speeded up by a factor 4.The velocity field is materialized by 2000 tracer particles that are advected by the flow. Render time: 3D parts - 1 hour 44 minutes2D parts - 1 hours 38 minutesColor scheme: Viridis by Nathaniel J. Smith, Stefan van der Walt and Eric Firinggithub.com/BIDS/colormapMusic: "Wandering Gaze" by Chasms@felteSee 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 2D shallow water equation by discretization (finite differences).C code: github.com/nilsberglund-orleans/YouTube-simulationshttps://www.idpoisson.fr/berglund/software.htmlMany thanks to Marco Mancini and Julian Kauth for helping me to accelerate my code!#shallowwater #waves #Earth
Sorting triangles with rattling sieves Nils Berglund 2024-09-21 | This simulation uses the same sorting device as in the video youtu.be/EcCHSkKQ8ic , including sieves that "rattle" to avoid clogging, with triangular instead of square particles. Another difference is that larger particles experience a larger friction force. This reduces the instability of the largest particles, in the rightmost bin, to some degree, although they still have trouble settling down. The sorting set-up consists in three grids of obstacles with decreasing space between them, acting as particle sieves. The obstacles rotate and exert a tangential force on the particles, in order to decrease clogging of the sieves. The angular speed of the obstacles varies periodically in time, in order to reduce the chance of particles getting stuck. In addition, the obstacles have a circular motion of small amplitude, which reduces the chances of particles getting stuck even more. The conveyor belt effect results from the segments forming the belt exerting a tangential force on the polygons, in addition to the normal force. The tangential force is proportional to the difference between the tangential speed of the polygon and the speed of the belt. To compute the force and torque of polygon j on polygon i, the code computes the distance of each vertex of polygon j to the faces of polygon i. If this distance is smaller than a threshold, the force increases linearly with a large spring constant. In addition, radial forces between the vertices of the polygons have been added, whenever a vertex of polygon j is not on a perpendicular to a face of polygon i. This is important, because otherwise triangles can approach each other from the vertices, and when one vertex moves sideways, it is suddenly strongly accelerated, causing numerical instability. A weak Lennard-Jones interaction between polygons has been added, as it seems to increase numerical stability. Unlike in some previous videos involving interacting polygons, there is no thermostat in this simulation. Instead, friction forces (both linear and angular) have been added for numerical stability. In addition, the particles are subject to a gravitational force directed downwards. The color of the polygons depends on their size.To save on computation time, particles are placed into a "hash grid", each cell of which contains between 3 and 10 particles. Then only the influence of other particles in the same or neighboring cells is taken into account for each particle. The Lennard-Jones potential is strongly repulsive at short distance, and mildly attracting at long distance. It is widely used as a simple yet realistic model for the motion of electrically neutral molecules. The force results from the repulsion between electrons due to Pauli's exclusion principle, while the attractive part is a more subtle effect appearing in a multipole expansion. For more details, see en.wikipedia.org/wiki/Lennard-Jones_potentialRender time: 56 minutes 56 secondsCompression: crf 23Color scheme: Turbo, by Anton Mikhailovgist.github.com/mikhailov-work/6a308c20e494d9e0ccc29036b28faa7a Music: Chipper Doodle - Electronic Light by Kevin MacLeod is licensed under a Creative Commons Attribution 4.0 licence. creativecommons.org/licenses/by/4.0Source: http://incompetech.com/music/royalty-free/index.html?isrc=USUAN1100854Artist: http://incompetech.comCurrent version of the C code used to make these animations: github.com/nilsberglund-orleans/YouTube-simulationshttps://www.idpoisson.fr/berglund/software.htmlSome outreach articles on mathematics:https://images.math.cnrs.fr/auteurs/nils-berglund/(in French, some with a Spanish translation)#molecular_dynamics #polygons #conveyor
Using seven sources to excite resonant modes in a circle Nils Berglund 2024-09-20 | Like the videos youtu.be/HawUSO5QYl4 and youtu.be/3fmdJAltpKA , this one shows waves in a circular cavity excited by pulsing sources in other, connected cavities. The difference is that here, there are seven smaller cavities instead of six or three.This video has two parts, showing the same evolution with two different color gradients:Wave height: 0:00Averaged wave energy: 1:24In 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: 13 minutes 59 secondsCompression: crf 23Color scheme: Part 1 - Magma by Nathaniel J. Smith and Stefan van der WaltPart 2 - Plasma by Nathaniel J. Smith and Stefan van der Waltgithub.com/BIDS/colormapMusic: "Aurora Borealis Expedition" by Asher Fulero@AsherFuleroSee 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-simulationshttps://www.idpoisson.fr/berglund/software.htmlMany thanks to Marco Mancini and Julian Kauth for helping me to accelerate my code!#wave #resonator
An asteroid impact in the North Atlantic Ocean, with tsunami-induced flooding Nils Berglund 2024-09-19 | This is an attempt at modeling an asteroid impact in the North Atlantic Ocean using a nonlinear wave equation. This model allows for some flooding of land areas. The vertical scale of both mountains and waves has been exaggerated. Compared to the earlier simulation youtu.be/th_TzO6RsvI , a mistake in the code has been corrected, which incorrectly computed the influence of the non-constant ocean depth - the east-west component of that force was using the north-south gradient of the depth, instead of the east-west gradient. The display of land masses has also been improved, by using a different criterion than before to distinguish land from sea. The main difference between the nonlinear equation and the linear wave equation is that the wave speed becomes larger when the water is shallower, which can lead to build-up of the waves in coastal and other shallow regions. This causes problems in the simulation, because land masses can lead to blow-up of the solution. This problem is circumvented here by replacing the continents by a repelling force field, plus a dissipative term. The initial state features velocities radiating outward from the impact point of the asteroid. The shallow water equations are nonlinear equations, which give a better description of the motion of water than the linear wave equation. In particular, unlike the linear wave equation, they conserve the total volume of water. The linear equation gives an approximation of the solutions, when the wave height remains close to its average over space. The equations used here include viscosity and dissipation, as described for instance in en.wikipedia.org/wiki/Shallow_water_equations#Non-conservative_form , including the Coriolis force. One difficulty is to model the wetting boundary, which separates regions that are under water and those which are not. This difficulty has been circumvented here by replacing the continents by a repulsive force field, directed downslope, instead of a sharp boundary. The video has four parts, showing simulations at two different speeds and with two different visualizations:Time lapse, 3D: 0:00Time lapse, 2D: 0:14Original speed, 3D: 0:28Original speed, 2D: 1:27The color hue and radial coordinate show the height of the water, on an exaggerated radial scale. The 2D parts use a projection in equirectangular coordinates. In the 3D parts, the point of view is slowly rotating around the Earth in a plane containing its center. In parts 1 and 2, the animation has been speeded up by a factor 4.The velocity field is materialized by 2000 tracer particles that are advected by the flow. Render time: 3D parts - 1 hour 55 minutes2D parts - 1 hours 29 minutesColor scheme: Viridis by Nathaniel J. Smith, Stefan van der Walt and Eric Firinggithub.com/BIDS/colormapMusic: "Virtual Light" by Houses of Heaven@felteSee 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 2D shallow water equation by discretization (finite differences).C code: github.com/nilsberglund-orleans/YouTube-simulationshttps://www.idpoisson.fr/berglund/software.htmlMany thanks to Marco Mancini and Julian Kauth for helping me to accelerate my code!#shallowwater #waves #Earth
Rattling the sieves in a particle sorter avoids clogging Nils Berglund 2024-09-18 | In the previous simulation youtu.be/ZA3gcbNm1l8 of a device sorting particles according to their size, particles sometimes got stuck in the sieves, resulting in other particles ending up in the wrong bin. This simulation implements a suggestion made by several viewers, which is to "rattle" the sieves in order to dislodge stuck particles.The sorting set-up consists in three grids of obstacles with decreasing space between them, acting as particle sieves. The obstacles rotate and exert a tangential force on the particles, in order to decrease clogging of the sieves. The angular speed of the obstacles varies periodically in time, in order to reduce the chance of particles getting stuck. In addition, the obstacles have a circular motion of small amplitude, which reduces the chances of particles getting stuck even more. The conveyor belt effect results from the segments forming the belt exerting a tangential force on the polygons, in addition to the normal force. The tangential force is proportional to the difference between the tangential speed of the polygon and the speed of the belt. To compute the force and torque of polygon j on polygon i, the code computes the distance of each vertex of polygon j to the faces of polygon i. If this distance is smaller than a threshold, the force increases linearly with a large spring constant. In addition, radial forces between the vertices of the polygons have been added, whenever a vertex of polygon j is not on a perpendicular to a face of polygon i. This is important, because otherwise triangles can approach each other from the vertices, and when one vertex moves sideways, it is suddenly strongly accelerated, causing numerical instability. A weak Lennard-Jones interaction between polygons has been added, as it seems to increase numerical stability. Unlike in some previous videos involving interacting polygons, there is no thermostat in this simulation. Instead, friction forces (both linear and angular) have been added for numerical stability. In addition, the particles are subject to a gravitational force directed downwards. The color of the polygons depends on their size.To save on computation time, particles are placed into a "hash grid", each cell of which contains between 3 and 10 particles. Then only the influence of other particles in the same or neighboring cells is taken into account for each particle. The Lennard-Jones potential is strongly repulsive at short distance, and mildly attracting at long distance. It is widely used as a simple yet realistic model for the motion of electrically neutral molecules. The force results from the repulsion between electrons due to Pauli's exclusion principle, while the attractive part is a more subtle effect appearing in a multipole expansion. For more details, see en.wikipedia.org/wiki/Lennard-Jones_potentialRender time: 50 minutes 22 secondsCompression: crf 23Color scheme: Turbo, by Anton Mikhailovgist.github.com/mikhailov-work/6a308c20e494d9e0ccc29036b28faa7a Music: "Parakeet" by Quincas Moreira@UCL1zFMJb0sthwdAlGjGbdyCurrent version of the C code used to make these animations: github.com/nilsberglund-orleans/YouTube-simulationshttps://www.idpoisson.fr/berglund/software.htmlSome outreach articles on mathematics:https://images.math.cnrs.fr/auteurs/nils-berglund/(in French, some with a Spanish translation)#molecular_dynamics #polygons #conveyor
Resonances in a circle, excited by four sources Nils Berglund 2024-09-17 | In this simulation of the wave equation, resonances are excited in a circular cavity, by sources placed in four secondary circular cavities. The design is inspired by cavity magnetrons, which are used to generate microwaves in microwave ovens. The color hue depends on the height of the wave. Render time: 6 minutes 52 secondsCompression: crf 23Color scheme: Inferno by Nathaniel J. Smith and Stefan van der Waltgithub.com/BIDS/colormapMusic: "U Make Me Feel" by MK2See 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-simulationshttps://www.idpoisson.fr/berglund/software.htmlMany thanks to Marco Mancini and Julian Kauth for helping me to accelerate my code!#wave #resonator
An asteroid impact in the Indian Ocean, with tsunami-induced flooding Nils Berglund 2024-09-16 | This is an attempt at modeling an asteroid impact in the Indian ocean using a nonlinear wave equation. This model allows for some flooding of land areas. The vertical scale of both mountains and waves has been exaggerated. A mistake in the code has been corrected, which incorrectly computed the influence of the non-constant ocean depth - the east-west component of that force was using the north-south gradient of the depth, instead of the east-west gradient. A version using the linear wave equation appeared previously on this channel, see youtu.be/DRsKobzGCVk The main difference in the nonlinear equation is that the wave speed becomes larger when the water is shallower, which can lead to build-up of the waves in coastal and other shallow regions. This causes problems in the simulation, because land masses can lead to blow-up of the solution. This problem is circumvented here by replacing the continents by a repelling force field, plus a dissipative term. The initial state features velocities radiating outward from the impact point of the asteroid. The shallow water equations are nonlinear equations, which give a better description of the motion of water than the linear wave equation. In particular, unlike the linear wave equation, they conserve the total volume of water. The linear equation gives an approximation of the solutions, when the wave height remains close to its average over space. The equations used here include viscosity and dissipation, as described for instance in en.wikipedia.org/wiki/Shallow_water_equations#Non-conservative_form , including the Coriolis force. One difficulty is to model the wetting boundary, which separates regions that are under water and those which are not. This difficulty has been circumvented here by replacing the continents by a repulsive force field, directed downslope, instead of a sharp boundary. The video has four parts, showing simulations at two different speeds and with two different visualizations:Time lapse, 3D: 0:00Time lapse, 2D: 0:13Original speed, 3D: 0:27Original speed, 2D: 1:21The color hue and radial coordinate show the height of the water, on an exaggerated radial scale. The 2D parts use a projection in equirectangular coordinates. In the 3D parts, the point of view is slowly rotating around the Earth in a plane containing its center. In parts 1 and 2, the animation has been speeded up by a factor 4.The velocity field is materialized by 2000 tracer particles that are advected by the flow. Render time: 3D parts - 1 hour 58 minutes2D parts - 1 hours 50 minutesColor scheme: Viridis by Nathaniel J. Smith, Stefan van der Walt and Eric Firinggithub.com/BIDS/colormapMusic: "Mast" by Silent PartnerSee 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 2D shallow water equation by discretization (finite differences).C code: github.com/nilsberglund-orleans/YouTube-simulationshttps://www.idpoisson.fr/berglund/software.htmlMany thanks to Marco Mancini and Julian Kauth for helping me to accelerate my code!#shallowwater #waves #Earth
A conveyor belt with shovels at lower friction Nils Berglund 2024-09-15 | This simulation uses the same conveyor belt with shovels as in the video youtu.be/beBpMuZ6cn0 , but the friction is four times lower. As a consequence, particles accelerate more when they fall, and also rotate more. Other differences are that the particles are pentagons instead of hexagons, have different radii, and there is no torque aligning neighboring polygons.The conveyor belt effect results from a combination of two factors: the shovels, and the fact that the segments forming the belt exert a tangential force on the polygons, in addition to the normal force. The tangential force is proportional to the difference between the tangential speed of the polygon and the speed of the belt. To compute the force and torque of polygon j on polygon i, the code computes the distance of each vertex of polygon j to the faces of polygon i. If this distance is smaller than a threshold, the force increases linearly with a large spring constant. In addition, radial forces between the vertices of the polygons have been added, whenever a vertex of polygon j is not on a perpendicular to a face of polygon i. This is important, because otherwise triangles can approach each other from the vertices, and when one vertex moves sideways, it is suddenly strongly accelerated, causing numerical instability. A weak Lennard-Jones interaction between polygons has been added, as it seems to increase numerical stability. Unlike in some previous videos involving interacting polygons, there is no thermostat in this simulation. Instead, friction forces (both linear and angular) have been added for numerical stability. In addition, the particles are subject to a gravitational force directed downwards. This simulation has two parts, showing the evolution with two different color gradients:Initial x position: 0:00Velocity: 2:24In the first part, the particles' color depends on their initial x position. This allows to get a sense of how particles are mixed while falling through the funnel. In the second part, the due depends on the direction in which the particles are moving, while the luminosity depends on their speed. Both quantities are averaged over a time interval. To save on computation time, particles are placed into a "hash grid", each cell of which contains between 3 and 10 particles. Then only the influence of other particles in the same or neighboring cells is taken into account for each particle. The Lennard-Jones potential is strongly repulsive at short distance, and mildly attracting at long distance. It is widely used as a simple yet realistic model for the motion of electrically neutral molecules. The force results from the repulsion between electrons due to Pauli's exclusion principle, while the attractive part is a more subtle effect appearing in a multipole expansion. For more details, see en.wikipedia.org/wiki/Lennard-Jones_potentialRender time: 1 hour 32 minutesCompression: crf 23Color scheme: Part 1 - Turbo, by Anton Mikhailovgist.github.com/mikhailov-work/6a308c20e494d9e0ccc29036b28faa7a Part 2 - Twilight by Bastian Bechtoldgithub.com/bastibe/twilightMusic: Cold Funk - Funkorama by Kevin MacLeod is licensed under a Creative Commons Attribution 4.0 licence. creativecommons.org/licenses/by/4.0/Source: http://incompetech.com/music/royalty-free/index.html?isrc=USUAN1100499Artist: http://incompetech.comCurrent version of the C code used to make these animations: github.com/nilsberglund-orleans/YouTube-simulationshttps://www.idpoisson.fr/berglund/software.htmlSome outreach articles on mathematics:https://images.math.cnrs.fr/auteurs/nils-berglund/(in French, some with a Spanish translation)#molecular_dynamics #polygons #conveyor
Refraction and reflection of a shock wave - Source in the lower speed medium Nils Berglund 2024-09-14 | This variant of the simulation youtu.be/fQV6vgF9vSI shows again shock waves created by a source traveling faster than the wave speed in a medium, near an interface to a medium with different wave speed. However, here the wave speeds have been reversed. The lower medium has relative refractive index 0.667 (the wave speed is 1.5 times faster in the lower part). As a result, part of the shock wave is refracted, changing its angle in the second medium, while part of it is reflected. The relative intensity of the reflected and refracted parts depends on the angle between shock wave and interface.This video has six parts, showing the evolution with three different speeds of the source, and two different color gradients for each speed:Speed 1, wave height: 0:00Speed 1, averaged wave energy: 0:44Speed 2, wave height: 1:36Speed 2, averaged wave energy: 2:19Speed 3, wave height: 3:12Speed 3, averaged wave energy: 3:56In parts 1, 3 and 5, the color hue depends on the height of the wave. In the other parts, 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.Render time: Part 1 - 17 minutes 5 secondsPart 2 - 15 minutes 58 secondsPart 3 - 16 minutes 2 secondsCompression: crf 23Color scheme: Part 1 - Viridis by Nathaniel J. Smith, Stefan van der Walt and Eric FiringPart 2 - Inferno by Nathaniel J. Smith and Stefan van der Waltgithub.com/BIDS/colormapMusic: "Stardrive" by Jeremy Blake@RedMeansRecordingSee 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-simulationshttps://www.idpoisson.fr/berglund/software.htmlMany thanks to Marco Mancini and Julian Kauth for helping me to accelerate my code!#wave #shock_wave #sonic_boom
An asteroid impact in the South Atlantic, modeled with a shallow water equation Nils Berglund 2024-09-13 | This is an attempt at modeling an asteroid impact in the South Atlantic ocean using a nonlinear wave equation. The main difference with the linear wave equation is that the wave speed becomes larger when the water is shallower, which can lead to build-up of the waves in coastal and other shallow regions. This causes problems in the simulation, because land masses can lead to blow-up of the solution. This problem is circumvented here by replacing the continents by a repelling force field, plus a dissipative term. The initial state features velocities radiating outward from the impact point of the asteroid. The shallow water equations are nonlinear equations, which give a better description of the motion of water than the linear wave equation. In particular, unlike the linear wave equation, they conserve the total volume of water. The linear equation gives an approximation of the solutions, when the wave height remains close to its average over space. The equations used here include viscosity and dissipation, as described for instance in en.wikipedia.org/wiki/Shallow_water_equations#Non-conservative_form , including the Coriolis force. One difficulty is to model the wetting boundary, which separates regions that are under water and those which are not. This difficulty has been circumvented here by replacing the continents by a repulsive force field, directed downslope, instead of a sharp boundary. The color hue and radial coordinate show the height of the water, on an exaggerated radial scale. The velocity field is materialized by 2000 tracer particles that are advected by the flow. Render time: 1 hour 51 minutesColor scheme: Viridis by Nathaniel J. Smith, Stefan van der Walt and Eric Firinggithub.com/BIDS/colormapMusic: "Traversing" by Godmode@GODMODEMUSICSee 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 2D shallow water equation by discretization (finite differences).C code: github.com/nilsberglund-orleans/YouTube-simulationshttps://www.idpoisson.fr/berglund/software.htmlMany thanks to Marco Mancini and Julian Kauth for helping me to accelerate my code!#shallowwater #waves #Earth
How does the particle size sorter work at very low friction? Nils Berglund 2024-09-12 | Not too well. This simulation uses the same sorting set-up as in the video youtu.be/ZA3gcbNm1l8 , but the friction here is four times weaker. As a result, the particles bounce more, and many of them end up in the wrong bin. The sorting set-up consists in three grids of obstacles with decreasing space between them, acting as particle sieves. The obstacles rotate and exert a tangential force on the particles, in order to decrease clogging of the sieves. The angular speed of the obstacles varies periodically in time, in order to reduce the chance of particles getting stuck. The conveyor belt effect results from the segments forming the belt exerting a tangential force on the polygons, in addition to the normal force. The tangential force is proportional to the difference between the tangential speed of the polygon and the speed of the belt. To compute the force and torque of polygon j on polygon i, the code computes the distance of each vertex of polygon j to the faces of polygon i. If this distance is smaller than a threshold, the force increases linearly with a large spring constant. In addition, radial forces between the vertices of the polygons have been added, whenever a vertex of polygon j is not on a perpendicular to a face of polygon i. This is important, because otherwise triangles can approach each other from the vertices, and when one vertex moves sideways, it is suddenly strongly accelerated, causing numerical instability. A weak Lennard-Jones interaction between polygons has been added, as it seems to increase numerical stability. Unlike in some previous videos involving interacting polygons, there is no thermostat in this simulation. Instead, friction forces (both linear and angular) have been added for numerical stability. In addition, the particles are subject to a gravitational force directed downwards. The color of the polygons depends on their size.To save on computation time, particles are placed into a "hash grid", each cell of which contains between 3 and 10 particles. Then only the influence of other particles in the same or neighboring cells is taken into account for each particle. The Lennard-Jones potential is strongly repulsive at short distance, and mildly attracting at long distance. It is widely used as a simple yet realistic model for the motion of electrically neutral molecules. The force results from the repulsion between electrons due to Pauli's exclusion principle, while the attractive part is a more subtle effect appearing in a multipole expansion. For more details, see en.wikipedia.org/wiki/Lennard-Jones_potentialRender time: 37 minutes 43 secondsCompression: crf 23Color scheme: Turbo, by Anton Mikhailovgist.github.com/mikhailov-work/6a308c20e494d9e0ccc29036b28faa7a Music: "Fanta Mankane" by the Mini Vandals@theminivandals1840featuring Mamadou Kolta and LassoCurrent version of the C code used to make these animations: github.com/nilsberglund-orleans/YouTube-simulationshttps://www.idpoisson.fr/berglund/software.htmlSome outreach articles on mathematics:https://images.math.cnrs.fr/auteurs/nils-berglund/(in French, some with a Spanish translation)#molecular_dynamics #polygons #conveyor
Using three sources to excite resonant modes in a circle Nils Berglund 2024-09-11 | Like the video youtu.be/HawUSO5QYl4 , this one shows waves in a circular cavity excited by pulsing sources in other, connected cavities. The difference is that here, there are three smaller cavities instead of six.This video has two parts, showing the same evolution with two different color gradients:Wave height: 0:00Averaged wave energy: 1:30In 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: 16 minutes 44 secondsCompression: crf 23Color scheme: Part 1 - Viridis by Nathaniel J. Smith, Stefan van der Walt and Eric Firinggithub.com/BIDS/colormapPart 2 - Turbo, by Anton Mikhailovgist.github.com/mikhailov-work/6a308c20e494d9e0ccc29036b28faa7a Music: "Mantis" by Quincas Moreira@QuincasMoreiraSee 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-simulationshttps://www.idpoisson.fr/berglund/software.htmlMany thanks to Marco Mancini and Julian Kauth for helping me to accelerate my code!#wave #resonator
Bloopers 19: Asteroid impact in the Indian Ocean, with a too strong Coriolis force Nils Berglund 2024-09-10 | This attempt at simulating an asteroid impact in the Indean Ocean uses a strong Coriolis force. As it turns out, the Coriolis force is so strong that the waves created by the impact tend to move inward after initially moving outward, earning this simulation a spot in my list of bloopers youtube.com/playlist?list=PLAZp3rbgWLo0CUHGuFfPUUrjdRWNCelDrThe shallow water equations are nonlinear equations, which give a better description of the motion of water than the linear wave equation. In particular, unlike the linear wave equation, they conserve the total volume of water. The linear equation gives an approximation of the solutions, when the wave height remains close to its average over space. The initial state features velocities radiating outward from the impact point of the asteroid. The equations used here include viscosity and dissipation, as described for instance in en.wikipedia.org/wiki/Shallow_water_equations#Non-conservative_form , including the Coriolis force. One difficulty is to model the wetting boundary, which separates regions that are under water and those which are not. This difficulty has been circumvented here by replacing the continents by a repulsive force field, directed downslope, instead of a sharp boundary. The video has four parts, showing simulations at two different speeds and with two different visualizations:Time lapse, 3D: 0:00Time lapse, 2D: 0:12Original speed, 3D: 0:25Original speed, 2D: 1:15The color hue and radial coordinate show the height of the water, on an exaggerated radial scale. The 2D parts use a projection in equirectangular coordinates. In the 3D parts, the point of view is slowly rotating around the Earth in a plane containing its center. In parts 1 and 2, the animation has been speeded up by a factor 4.The velocity field is materialized by 2000 tracer particles that are advected by the flow. Render time: 3D parts - 2 hour 38 minutes2D parts - 3 hours 20 minutesColor scheme: Viridis by Nathaniel J. Smith, Stefan van der Walt and Eric Firinggithub.com/BIDS/colormapMusic: "Grand Avenue" by Text Me Records - Bobby Renz@socialxworkSee 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 2D shallow water equation by discretization (finite differences).C code: github.com/nilsberglund-orleans/YouTube-simulationshttps://www.idpoisson.fr/berglund/software.htmlMany thanks to Marco Mancini and Julian Kauth for helping me to accelerate my code!#shallowwater #waves #Earth
Sorting particles at lower friction Nils Berglund 2024-09-09 | This simulation uses a similar set-up for sorting particles as in the video youtu.be/-k5cxxNL2po , but the friction (both linear and angular) is about 4 times smaller. This allows the particles to accelerate more when falling, to bounce when hitting a conveyor or other particles, and to rotate more. The drawback is that the largest particles, in the rightmost bin, are more "jumpy", which is probably due to a bit of numerical instability. The sorting set-up consists in three grids of obstacles with decreasing space between them, acting as particle sieves. The obstacles rotate and exert a tangential force on the particles, in order to decrease clogging of the sieves. The angular speed of the obstacles varies periodically in time, in order to reduce the chance of particles getting stuck. The conveyor belt effect results from the segments forming the belt exerting a tangential force on the polygons, in addition to the normal force. The tangential force is proportional to the difference between the tangential speed of the polygon and the speed of the belt. To compute the force and torque of polygon j on polygon i, the code computes the distance of each vertex of polygon j to the faces of polygon i. If this distance is smaller than a threshold, the force increases linearly with a large spring constant. In addition, radial forces between the vertices of the polygons have been added, whenever a vertex of polygon j is not on a perpendicular to a face of polygon i. This is important, because otherwise triangles can approach each other from the vertices, and when one vertex moves sideways, it is suddenly strongly accelerated, causing numerical instability. A weak Lennard-Jones interaction between polygons has been added, as it seems to increase numerical stability. Unlike in some previous videos involving interacting polygons, there is no thermostat in this simulation. Instead, friction forces (both linear and angular) have been added for numerical stability. In addition, the particles are subject to a gravitational force directed downwards. The color of the polygons depends on their size.To save on computation time, particles are placed into a "hash grid", each cell of which contains between 3 and 10 particles. Then only the influence of other particles in the same or neighboring cells is taken into account for each particle. The Lennard-Jones potential is strongly repulsive at short distance, and mildly attracting at long distance. It is widely used as a simple yet realistic model for the motion of electrically neutral molecules. The force results from the repulsion between electrons due to Pauli's exclusion principle, while the attractive part is a more subtle effect appearing in a multipole expansion. For more details, see en.wikipedia.org/wiki/Lennard-Jones_potentialRender time: 1 hour 28 minutesCompression: crf 23Color scheme: Turbo, by Anton Mikhailovgist.github.com/mikhailov-work/6a308c20e494d9e0ccc29036b28faa7a Music: "Parasail" by Silent PartnerCurrent version of the C code used to make these animations: github.com/nilsberglund-orleans/YouTube-simulationshttps://www.idpoisson.fr/berglund/software.htmlSome outreach articles on mathematics:https://images.math.cnrs.fr/auteurs/nils-berglund/(in French, some with a Spanish translation)#molecular_dynamics #polygons #conveyor
Using six sources to excite resonant modes in a circle Nils Berglund 2024-09-08 | Like the video youtu.be/HawUSO5QYl4 , this one shows waves in a circular cavity excited by pulsing sources in other, connected cavities. The difference is that here, there are six smaller cavities instead of only one. In this way, the waves keep a six-fold symmetry (or rather, for the first part, a three-fold symmetry with an inversion). This video has two parts, showing the same evolution with two different color gradients:Wave height: 0:00Averaged wave energy: 1:20In 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: 15 minutes 24 secondsCompression: crf 23Color scheme: Part 1 - Twilight by Bastian Bechtoldgithub.com/bastibe/twilightPart 2 - Viridis by Nathaniel J. Smith, Stefan van der Walt and Eric Firinggithub.com/BIDS/colormapMusic: "Alone" by Emmit Fenn@emmitfennSee 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-simulationshttps://www.idpoisson.fr/berglund/software.htmlMany thanks to Marco Mancini and Julian Kauth for helping me to accelerate my code!#wave #resonator
An improved device sorting particles by size Nils Berglund 2024-09-07 | This variant of the videos youtu.be/qozGp_KR_DA and youtu.be/feEDk9e4HPM of a device sorting particles by size contains two improvements. The first one, suggested in a comment, is to have a conveyor belt bringing the particles from the opposite side, and dropping them more to the left, in order to use a longer stretch of the sieves. The second improvement is that the circular obstacles forming the sieves do not rotate at constant speed, but rotate sometimes faster and sometimes slower. This reduces the number of particles getting stuck in the sieves, in a similar way as "rattling" the sieves now and then would. The sorting set-up consists in three grids of obstacles with decreasing space between them, acting as particle sieves. The conveyor belt effect results from the segments forming the belt exerting a tangential force on the polygons, in addition to the normal force. The tangential force is proportional to the difference between the tangential speed of the polygon and the speed of the belt. To compute the force and torque of polygon j on polygon i, the code computes the distance of each vertex of polygon j to the faces of polygon i. If this distance is smaller than a threshold, the force increases linearly with a large spring constant. In addition, radial forces between the vertices of the polygons have been added, whenever a vertex of polygon j is not on a perpendicular to a face of polygon i. This is important, because otherwise triangles can approach each other from the vertices, and when one vertex moves sideways, it is suddenly strongly accelerated, causing numerical instability. A weak Lennard-Jones interaction between polygons has been added, as it seems to increase numerical stability. Unlike in some previous videos involving interacting polygons, there is no thermostat in this simulation. Instead, friction forces (both linear and angular) have been added for numerical stability. In addition, the particles are subject to a gravitational force directed downwards. The color of the polygons depends on their size.To save on computation time, particles are placed into a "hash grid", each cell of which contains between 3 and 10 particles. Then only the influence of other particles in the same or neighboring cells is taken into account for each particle. The Lennard-Jones potential is strongly repulsive at short distance, and mildly attracting at long distance. It is widely used as a simple yet realistic model for the motion of electrically neutral molecules. The force results from the repulsion between electrons due to Pauli's exclusion principle, while the attractive part is a more subtle effect appearing in a multipole expansion. For more details, see en.wikipedia.org/wiki/Lennard-Jones_potentialRender time: 1 hour 17 minutesCompression: crf 23Color scheme: Turbo, by Anton Mikhailovgist.github.com/mikhailov-work/6a308c20e494d9e0ccc29036b28faa7a Music: "Running Errands" by TrackTribe@TrackTribeCurrent version of the C code used to make these animations: github.com/nilsberglund-orleans/YouTube-simulationshttps://www.idpoisson.fr/berglund/software.htmlSome outreach articles on mathematics:https://images.math.cnrs.fr/auteurs/nils-berglund/(in French, some with a Spanish translation)#molecular_dynamics #polygons #conveyor
An asteroid impact in the North Atlantic, modeled with a shallow water equation Nils Berglund 2024-09-06 | This is an attempt at modeling an asteroid impact in the North Atlantic ocean using a nonlinear wave equation. A version using the linear wave equation appeared previously on this channel, see youtu.be/mjSFAVcpglA The main difference in the nonlinear equation is that the wave speed becomes larger when the water is shallower, which can lead to build-up of the waves in coastal and other shallow regions. This causes problems in the simulation, because land masses can lead to blow-up of the solution. This problem is circumvented here by replacing the continents by a repelling force field, plus a dissipative term. The initial state features velocities radiating outward from the impact point of the asteroid. One issue of the simulation is that waves tend to slow down as time goes on. This may be related to parameter values used in the equation. The shallow water equations are nonlinear equations, which give a better description of the motion of water than the linear wave equation. In particular, unlike the linear wave equation, they conserve the total volume of water. The linear equation gives an approximation of the solutions, when the wave height remains close to its average over space. The equations used here include viscosity and dissipation, as described for instance in en.wikipedia.org/wiki/Shallow_water_equations#Non-conservative_form , including the Coriolis force. One difficulty is to model the wetting boundary, which separates regions that are under water and those which are not. This difficulty has been circumvented here by replacing the continents by a repulsive force field, directed downslope, instead of a sharp boundary. The video has four parts, showing simulations at two different speeds and with two different visualizations:Time lapse, 3D: 0:00Time lapse, 2D: 0:16Original speed, 3D: 0:33Original speed, 2D: 1:39The color hue and radial coordinate show the height of the water, on an exaggerated radial scale. The 2D parts use a projection in equirectangular coordinates. In the 3D parts, the point of view is slowly rotating around the Earth in a plane containing its center. In parts 1 and 2, the animation has been speeded up by a factor 4.The velocity field is materialized by 2000 tracer particles that are advected by the flow. Render time: 3D parts - 2 hour 38 minutes2D parts - 3 hours 20 minutesColor scheme: Viridis by Nathaniel J. Smith, Stefan van der Walt and Eric Firinggithub.com/BIDS/colormapMusic: "Sink Whole Dream" by the 129ersSee 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 2D shallow water equation by discretization (finite differences).C code: github.com/nilsberglund-orleans/YouTube-simulationshttps://www.idpoisson.fr/berglund/software.htmlMany thanks to Marco Mancini and Julian Kauth for helping me to accelerate my code!#shallowwater #waves #Earth
Exciting resonant modes in a circle Nils Berglund 2024-09-05 | 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:00Averaged wave energy: 1:39In 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 secondsCompression: crf 23Color scheme: Part 1 - Twilight by Bastian Bechtoldgithub.com/bastibe/twilightPart 2 - Magma by Nathaniel J. Smith and Stefan van der Waltgithub.com/BIDS/colormapMusic: "Highway 5" by TrackTribe@TrackTribeSee 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-simulationshttps://www.idpoisson.fr/berglund/software.htmlMany thanks to Marco Mancini and Julian Kauth for helping me to accelerate my code!#wave #resonator
Sorting hexagons by size Nils Berglund 2024-09-04 | This simulation uses the same set-up as seen in the video youtu.be/qozGp_KR_DA to sort particles according to their size, but with hexagons instead of squares. The set-up consists in three grids of obstacles with decreasing space between them, acting as particle sieves. The obstacles rotate and exert a tangential force on the particles, in order to decrease clogging of the sieves.The conveyor belt effect results from the segments forming the belt exerting a tangential force on the polygons, in addition to the normal force. The tangential force is proportional to the difference between the tangential speed of the polygon and the speed of the belt. To compute the force and torque of polygon j on polygon i, the code computes the distance of each vertex of polygon j to the faces of polygon i. If this distance is smaller than a threshold, the force increases linearly with a large spring constant. In addition, radial forces between the vertices of the polygons have been added, whenever a vertex of polygon j is not on a perpendicular to a face of polygon i. This is important, because otherwise triangles can approach each other from the vertices, and when one vertex moves sideways, it is suddenly strongly accelerated, causing numerical instability. A weak Lennard-Jones interaction between polygons has been added, as it seems to increase numerical stability. Unlike in some previous videos involving interacting polygons, there is no thermostat in this simulation. Instead, friction forces (both linear and angular) have been added for numerical stability. In addition, the particles are subject to a gravitational force directed downwards. The color of the polygons depends on their size.To save on computation time, particles are placed into a "hash grid", each cell of which contains between 3 and 10 particles. Then only the influence of other particles in the same or neighboring cells is taken into account for each particle. The Lennard-Jones potential is strongly repulsive at short distance, and mildly attracting at long distance. It is widely used as a simple yet realistic model for the motion of electrically neutral molecules. The force results from the repulsion between electrons due to Pauli's exclusion principle, while the attractive part is a more subtle effect appearing in a multipole expansion. For more details, see en.wikipedia.org/wiki/Lennard-Jones_potentialRender time: 50 minutes 6 secondsCompression: crf 23Color scheme: Turbo, by Anton Mikhailovgist.github.com/mikhailov-work/6a308c20e494d9e0ccc29036b28faa7a Music: "Crystal" by Vibe TracksCurrent version of the C code used to make these animations: github.com/nilsberglund-orleans/YouTube-simulationshttps://www.idpoisson.fr/berglund/software.htmlSome outreach articles on mathematics:https://images.math.cnrs.fr/auteurs/nils-berglund/(in French, some with a Spanish translation)#molecular_dynamics #polygons #conveyor
Video #1300: What if the Earths oceans were much shallower? Nils Berglund 2024-09-03 | This is the 1300th video published on this channel, not counting some videos published multiple times because of corruption issues. Thanks again to all you viewers for your fidelity, comments and suggestions.This is a first simulation of the shallow water equations on a model for the Earth. The wave speed in these equations becomes larger when the water is shallower, which can lead to build-up of the waves in coastal and other shallow regions. The most important parameter in shallow water equations is the ratio between the size of depth variation and the average depth. This videos uses two different values of this parameter, the last two parts having a larger ratio. In both cases, this ratio is quite large compared to what applies to the Earth, making the effect of shallower parts more important.The shallow water equations are nonlinear equations, which give a better description of the motion of water than the linear wave equation. In particular, unlike the linear wave equation, they conserve the total volume of water. The linear equation gives an approximation of the solutions, when the wave height remains close to its average over space. The equations used here include viscosity and dissipation, as described for instance in en.wikipedia.org/wiki/Shallow_water_equations#Non-conservative_form , including the Coriolis force. One difficulty is to model the wetting boundary, which separates regions that are under water and those which are not. This difficulty has been circumvented here by replacing the continents by a repulsive force field, directed downslope, instead of a sharp boundary. The video has four parts, showing simulations with two different parameter values and two different visualizations:Smaller ratio, 3D: 0:00Smaller ratio, 2D: 0:48Larger ratio, 3D: 1:36Larger ratio, 2D: 2:24The color hue and radial coordinate show the height of the water, on an exaggerated radial scale. The 2D parts use a projection in equirectangular coordinates. In the 3D parts, the point of view is slowly rotating around the Earth in a plane containing its center. The position of the sphere is represented by a line starting above the north pole, and pointing in a fixed direction. The velocity field is materialized by 2000 tracer particles that are advected by the flow. Render time: Part 1 - 1 hour 30 minutesPart 2 - 1 hour 14 minutesPart 3 - 1 hour 8 minutesPart 4 - 1 hour 36 minutesColor scheme: Viridis by Nathaniel J. Smith, Stefan van der Walt and Eric Firinggithub.com/BIDS/colormapMusic: "Gymnopédie No. 1" by Erik SatieSee 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 2D shallow water equation by discretization (finite differences).C code: github.com/nilsberglund-orleans/YouTube-simulationshttps://www.idpoisson.fr/berglund/software.htmlMany thanks to Marco Mancini and Julian Kauth for helping me to accelerate my code!#shallowwater #waves #Earth
Classics revisited: A hexagonal parabolic resonator Nils Berglund 2024-09-02 | This is a remake of the video youtu.be/lVDwH5PfEGY of a parabolic resonator with six sides, in much higher resolution, and with a pulsing source. The resonator is made of six parts of parabolas, which share the central point as a common focus. Whenever a circular wave centered at the focus hits a parabola, it is transformed into a linear wave. When this linear wave hits another parabola, it is transformed back into a circular wave, converging at the focal point. This pattern repeats almost periodically, except that diffraction effects (and probably also some numerical dispersion) slowly degrade the pattern. This video has two parts, showing the same evolution with two different color gradients:Wave height: 0:00Averaged wave energy: 1:52In 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: 22 minutes 14 secondsCompression: crf 23Color scheme: Part 1 - Twilight by Bastian Bechtoldgithub.com/bastibe/twilightPart 2 - Turbo, by Anton Mikhailovgist.github.com/mikhailov-work/6a308c20e494d9e0ccc29036b28faa7a Music: "Into the Void" by TrackTribe@TrackTribeSee 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-simulationshttps://www.idpoisson.fr/berglund/software.htmlMany thanks to Marco Mancini and Julian Kauth for helping me to accelerate my code!#wave #resonator
Refraction and reflection of a shock wave Nils Berglund 2024-08-31 | This simulation has been suggested in a comment to a previous video. It shows the shock wave created by a source traveling faster than the wave speed in a medium. The wave encounters an interface to a different medium of relative refractive index 1.5 (the wave speed is 1.5 times slower in the lower part). As a result, part of the shock wave is refracted, changing its angle in the second medium, while part of it is reflected.This video has six parts, showing the evolution with three different speeds of the source, and two different color gradients for each speed:Speed 1, wave height: 0:00Speed 1, averaged wave energy: 0:44Speed 2, wave height: 1:36Speed 2, averaged wave energy: 2:19Speed 3, wave height: 3:12Speed 3, averaged wave energy: 3:56In parts 1, 3 and 5, the color hue depends on the height of the wave. In the other parts, 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.Render time: Part 1 - 16 minutes 29 secondsPart 2 - 16 minutes 13 secondsPart 3 - 16 minutes 15 secondsCompression: crf 23Color scheme: Part 1 - Viridis by Nathaniel J. Smith, Stefan van der Walt and Eric FiringPart 2 - Inferno by Nathaniel J. Smith and Stefan van der Waltgithub.com/BIDS/colormapMusic: "Eye Do" by Jeremy Blake@RedMeansRecordingSee 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-simulationshttps://www.idpoisson.fr/berglund/software.htmlMany thanks to Marco Mancini and Julian Kauth for helping me to accelerate my code!#wave #shock_wave #sonic_boom
What happened to your luggage at the airport? Nils Berglund 2024-08-30 | Last time I took a plane, I had a short connection, and guess what? My luggage was delayed. I ended up recovering my luggage all right, but since I have been making videos on conveyor belts recently, I realized that airport luggage sorting systems are a good example of the use of conveyor belts (as we all know thanks to John McClane). The example shown here is of course not to be taken too seriously.The conveyor belt effect results from a combination of two factors: the shovels, and the fact that the segments forming the belt exert a tangential force on the polygons, in addition to the normal force. The tangential force is proportional to the difference between the tangential speed of the polygon and the speed of the belt. To compute the force and torque of polygon j on polygon i, the code computes the distance of each vertex of polygon j to the faces of polygon i. If this distance is smaller than a threshold, the force increases linearly with a large spring constant. In addition, radial forces between the vertices of the polygons have been added, whenever a vertex of polygon j is not on a perpendicular to a face of polygon i. This is important, because otherwise triangles can approach each other from the vertices, and when one vertex moves sideways, it is suddenly strongly accelerated, causing numerical instability. A weak Lennard-Jones interaction between polygons has been added, as it seems to increase numerical stability. Unlike in some previous videos involving interacting polygons, there is no thermostat in this simulation. Instead, friction forces (both linear and angular) have been added for numerical stability. In addition, the particles are subject to a gravitational force directed downwards. The color of the polygons depends on the time they are created, and changes slowly over time.To save on computation time, particles are placed into a "hash grid", each cell of which contains between 3 and 10 particles. Then only the influence of other particles in the same or neighboring cells is taken into account for each particle. The Lennard-Jones potential is strongly repulsive at short distance, and mildly attracting at long distance. It is widely used as a simple yet realistic model for the motion of electrically neutral molecules. The force results from the repulsion between electrons due to Pauli's exclusion principle, while the attractive part is a more subtle effect appearing in a multipole expansion. For more details, see en.wikipedia.org/wiki/Lennard-Jones_potentialRender time: 40 minutes 47 secondsCompression: crf 23Color scheme: Turbo, by Anton Mikhailovgist.github.com/mikhailov-work/6a308c20e494d9e0ccc29036b28faa7a Music: Dispersion Relation by Kevin MacLeod is licensed under a Creative Commons Attribution 4.0 licence. creativecommons.org/licenses/by/4.0/Source: http://incompetech.com/music/royalty-free/index.html?isrc=USUAN1100258Artist: http://incompetech.comCurrent version of the C code used to make these animations: github.com/nilsberglund-orleans/YouTube-simulationshttps://www.idpoisson.fr/berglund/software.htmlSome outreach articles on mathematics:https://images.math.cnrs.fr/auteurs/nils-berglund/(in French, some with a Spanish translation)#molecular_dynamics #polygons #conveyor
