spin:esc202_fs2022
Differences
This shows you the differences between two versions of the page.
| Both sides previous revision Previous revision Next revision | Previous revision | ||
|
spin:esc202_fs2022 [2022/03/14 14:52] stadel [Lectures] |
spin:esc202_fs2022 [2022/04/11 14:25] (current) stadel [Lectures] |
||
|---|---|---|---|
| Line 15: | Line 15: | ||
| 14. Mar. 2022: {{ :spin:sphkernels.pdf | SPH - Equations and Kernels}} | 14. Mar. 2022: {{ :spin:sphkernels.pdf | SPH - Equations and Kernels}} | ||
| + | |||
| + | 21. Mar. 2022: {{ :spin:sphmainequations.pdf | SPH - Main numerical equations}} | ||
| + | |||
| + | 28. Mar. 2022: {{ :spin:sphwindtunnel.pdf |SPH - Wind Tunnel}} | ||
| + | |||
| + | 4. Apr. 2022: {{ :spin:2-disingmodel.pdf |2-D Ising Model}} | ||
| + | |||
| + | 11. Apr. 2022: {{ :spin:travelling_salesperson.pdf |Traveling Salesperson}} | ||
| ====== Assignments ====== | ====== Assignments ====== | ||
| Line 69: | Line 77: | ||
| 3. For each particle calculate the “top-hat” density from the 32 nearest neighbors and plot it using a colormap. Use a priority queue with a heap algorithm to implement “replace” and “max” functions (you can make use of the heapq library for python). We will need the density to implement SPH so it needs to be well tested. Also, make sure it can work with periodic boundary conditions! **(Hand in by March 13 2022)** | 3. For each particle calculate the “top-hat” density from the 32 nearest neighbors and plot it using a colormap. Use a priority queue with a heap algorithm to implement “replace” and “max” functions (you can make use of the heapq library for python). We will need the density to implement SPH so it needs to be well tested. Also, make sure it can work with periodic boundary conditions! **(Hand in by March 13 2022)** | ||
| + | |||
| + | 4. Now, calculate the density using the Monaghan kernel defined in the lecture. Plot and compare to the density you get from the “top-hat” kernel (the Monaghan result should be a little smoother). **(Hand in by March 20 2022)** | ||
| + | |||
| + | 5. Implement the SPH code introduced in the lecture to simulate the movement of the particles in the box, following the laws of hydrodynamics. You can start with a random particle distribution, or even simulate a "blast" with initial conditions where a particle somewhere has a high energy. **(Hand in by March 27 2022)** | ||
| + | |||
| + | 6. 2D Ising Model: Using the Metropolis algorithm, plot the mean magnetization of a N by N grid of spins (+1 and -1) depending on the temperature. Visualise the spin state of the grid at different temperatures **(due April 11th)**! | ||
| + | |||
| + | 7. Metropolis algorithm II: Implement the traveling merchant problem **(due April 17th)**! Have a look at some of the data found at http://comopt.ifi.uni-heidelberg.de/software/TSPLIB95/tsp/ to test your algorithm! | ||
spin/esc202_fs2022.1647265971.txt.gz · Last modified: 2022/03/14 14:52 by stadel
Page Tools
Except where otherwise noted, content on this wiki is licensed under the following license: CC Attribution-Share Alike 4.0 International
