Point vortex dynamics simulation

Versión 1.2.0.0 (8,92 KB) por Tom Ashbee
Integration of N 2D point vortices in a cylinder, using an adaptive 4th order Runge-Kutta scheme.
5.0
(2)
810 descargas
Actualizado 22 mar 2015

Ver licencia

The point vortex model was introduced by Helmholtz in 1867 and it was subsequently proved by Lin in 1942 that the motion of N vortices in a bounded domain is a Hamiltonian system (sometime called the `Kirchhoff-Routh path function'). The statistical mechanics of the N vortices is remarkable due to the fact that the system exhibits negative temperature.
The integrator used is a custom adaptive 4th order Runge-Kutta scheme which ensures the convergence of vortex positions to below a tolerance parameter before advancing to the next time step. In this way the energy and angular momentum of the system (the only two invariants) are conserved to high precision.

Citar como

Tom Ashbee (2024). Point vortex dynamics simulation (https://www.mathworks.com/matlabcentral/fileexchange/49103-point-vortex-dynamics-simulation), MATLAB Central File Exchange. Recuperado .

Compatibilidad con la versión de MATLAB
Se creó con R2013a
Compatible con cualquier versión
Compatibilidad con las plataformas
Windows macOS Linux
Categorías
Más información sobre Fluid Dynamics en Help Center y MATLAB Answers.

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

INT_VORTEX/

Versión Publicado Notas de la versión
1.2.0.0

Similar to the previous update, the double sum in the energy (energy.m) has been completely replaced by vector operations. Gives O(2) speed up in energy calc (if N>100) => only really useful for microcanonical statistical mechanics calculations.
1.1.0.0

The double sum in the equations of motion (eqns_of_motion.m) has been completely replaced by vector operations. This gives speed improvements of around O(6).
1.0.0.0