spin:esc202_fs2020

This shows you the differences between two versions of the page.

Both sides previous revision Previous revision | |||

spin:esc202_fs2020 [2020/05/04 08:11] stadel [Lectures] |
spin:esc202_fs2020 [2020/05/04 08:12] (current) stadel [List of assignments] |
||
---|---|---|---|

Line 97: | Line 97: | ||

4. **Sedov Taylor Explosion in SPH**: Set up a uniform grid of SPH particles in a 2-D unit cell. Make it so that there are an odd number of particle on a side, such that you have a particle centered exactly in the center of the unit cell. All particle should have v = 0 and mass = 1/N, where N is the total number of particles such that the mass of the fluid in the unit cell is 1. If we use gamma=2, then for e=1 we have c=sqrt(2) for all particles and the SPH code should maintain the uniform density as everything should be in pressure equilibrium. Test this. Now set the specific internal energy of the central particle to e = 100! Now the gas should react by producing the Sedov-Taylor blast wave. You will have to think about the appropriate timestep for the simulation. Think back to the Courant Condition from the first semester, or simply experiment with different values (hint: you know the sound speed, but what is the grid speed? What defines a resolution length scale here? h? So could we regard h/delta_t as a "grid speed"? | 4. **Sedov Taylor Explosion in SPH**: Set up a uniform grid of SPH particles in a 2-D unit cell. Make it so that there are an odd number of particle on a side, such that you have a particle centered exactly in the center of the unit cell. All particle should have v = 0 and mass = 1/N, where N is the total number of particles such that the mass of the fluid in the unit cell is 1. If we use gamma=2, then for e=1 we have c=sqrt(2) for all particles and the SPH code should maintain the uniform density as everything should be in pressure equilibrium. Test this. Now set the specific internal energy of the central particle to e = 100! Now the gas should react by producing the Sedov-Taylor blast wave. You will have to think about the appropriate timestep for the simulation. Think back to the Courant Condition from the first semester, or simply experiment with different values (hint: you know the sound speed, but what is the grid speed? What defines a resolution length scale here? h? So could we regard h/delta_t as a "grid speed"? | ||

+ | |||

+ | 5. **2-d Ising Model** | ||

+ | |||

+ | 6. **2-d Travelling Salesman Problem** | ||

+ |

spin/esc202_fs2020.txt ยท Last modified: 2020/05/04 08:12 by stadel

Except where otherwise noted, content on this wiki is licensed under the following license: CC Attribution-Share Alike 4.0 International