stadelweb

This page is read only. You can view the source, but not change it. Ask your administrator if you think this is wrong.
====== Simulations in the Natural Sciences I ======                                                                              
{{:spin:combined.jpg?800x100|}}                                                                                                                          
                                                                                                                                                           
**ESC201:**Fall 2021: Monday Lecture: 13:00-14:00 Exercises: 14:00-17:00 in **online only**

**TAs**: Stefan Schafroth, Sebastian Schulz (**sebastian.schulz@uzh.ch**), Noah Kubli (**noah.kubli@uzh.ch**)

===== Corona Virus Measures =====
**
Dear Students,
**

Video lectures will be available (see below).


Exercises will take place via Microsoft Teams :[[https://teams.microsoft.com/l/team/19%3aMRQsJW_dmDXZEhh1zumrUqh36j4D64ZgArPrkNMQxFM1%40thread.tacv2/conversations?groupId=127eb34a-b6e8-44da-bc4c-4ad8df3a1924&tenantId=c7e438db-e462-4c22-a90a-c358b16980b3|ESC 201]]



Joachim Stadel
----
===== Video Lectures =====

Videos for the lectures can be found here (they are too large for my Wiki):

[[http://www.ics.uzh.ch/~stadel/krone/public_downloads/esc201|ESC201 Video downloads]]


====== Lectures ======

20. Sep. 2021: {{ :spin:hs2021_week1.pdf |First Lecture}}

27. Sep. 2021: {{ :spin:hs2021_week2.pdf |Newton's Method and Kepler's Equation}}

4. Oct. 2021: {{ :spin:hs2021_week3.pdf |Population Growth, Chaos and Fractals}}

11. Oct. 2021: {{ :spin:hs2021_week4.pdf |Complex Fractals, Start of ODEs}}

18. Oct. 2021: {{ :spin:hs2021_week5.pdf |Ordinary Differential Equations: Predator-Prey}}

25. Oct. 2021: {{ :spin:hs2021_week6.pdf |Symplectic Integrators: Leap-frog and the Simple Pendulum}}

1. Nov. 2021: {{ :spin:hs2021_week7.pdf |Many Bodies: a digital Orrery}}

8. Nov. 2021: {{ :spin:hs2021_week8.pdf |Partial Differential Equations: Jacobi & SOR}}

15. Nov. 2021: {{ :spin:hs2021_week9.pdf |Bi-linear and Bi-cubic interpolation and Electron Beams}}

22. Nov. 2021: {{ :spin:hs2021_week10.pdf |Diffusion Equation and Numerical Stability}}

29. Nov. 2021: {{ :spin:hs2021_week11.pdf |Hyperbolic PDEs}}

6. Dec. 2021: {{ :spin:hs2021_week12.pdf |Finite Volume Method, Modified Equation, 2-D Advection}}

13. Dec. 2021: {{ :spin:hs2021_week13.pdf |1-D Hydrodynamics and Riemann Solvers}}




====== Assignments ======

Your solutions should be handed in 13 days after each assignment has been given, i.e. **Sunday night by 21:00** two weeks after the Monday lecture.
Assignments should be **individual** and should be in python and **provide a 
correct virtual environment!** (if you stick to standard libraries like numpy, matplotlib, scipy... **you can also just submit your python source code together with the output** of your program)

For help getting started with virtual environments, please read carefully [[https://towardsdatascience.com/all-you-need-to-know-about-python-virtual-environments-9b4aae690f97|Python Virtual Environments for Pip]] and [[https://uoa-eresearch.github.io/eresearch-cookbook/recipe/2014/11/20/conda/|Python Virtual Environments for Conda]].

You should email 3 things to Sebastian (**sebastian.schulz@uzh.ch**):

  - The working **python source code**
  - The **requirements.txt** file for your virtual environment (if you used a virtual environment)
  - A **.pdf** or **.png** image or animation of the output of your program

Template: {{template.zip}}

** Instructions: **

Please add the names of the people you work together (if you do) to the comment section of your python scripts.

Create a virtual environment using
    
Pip
    
- run **virtualenv yourenv_name** to create a virtual environment

- run **source yourenv_name/bin/activate** to activate yourenv_name

- install necessary libraries that you want using **pip install package_name**

- work in that directory, get your outputs (*.pdf, *.png, *jpeg, *.mp4, etc...)

- run **pip freeze > requirements.txt** to get your list of libraries
      
Conda

- run **conda create -n yourenvname python=x.x anaconda** to create a virtual environment

- run **source activate yourenvname** to activate yourenv_name

- install necessary libraries that you want using **conda install -n yourenv_name package_name**

- work in that directory, get your outputs (*.pdf, *.png, *jpeg, *.mp4, etc...)

- run **conda list --export > requirements.txt** to get your list of libraries

====== List of assignments ======

  - Implement the bisection root finding method presented in the first lecture for a function f of your choice (e.g. f(x) = ax³ + bx² + cx + d, or x<sup>x</sup> - 100 = 0)! (**to submit by 3 October 20201, 9pm**)
  - Plot (and/or animate) the elliptical orbit of a planet around the sun by repeatedly solving Kepler's equation with Newton's method (or the bisection method), as explained in the lecture! (**to submit by 10 October 2021, 9pm**) 
  - Draw a Feigenbaum diagram that results from solving the logistic equation (**to submit by 17 October, 2021, 9pm**). (Optional: Implement a function that allows you to zoom into the Feigenbaum diagram)
  - Fractals: Draw some Julia sets with various constants c (you can start with the Mandelbrot set as it was explained in the lecture and the exercise class)! (**to submit by 24 October, 2021, 9pm**)
  - Ordinary Differential Equations: Solve the Lotka-Volterra equation using the Euler method and the midpoint Runge-Kutta method (optional: 4th order Runge Kutta method) and compare the results. Make two plots: the time dependence of both populations (mice and foxes), and the phase diagram using different initial conditions (**to submit by 31 October, 2021, 9pm**).
  - Symplectic Integrators: Use the Leap-Frog method to make a phase plot (p vs q) of the harmonic oscillator for different total energies. Compare the results with what you get using the Forward Euler method and the midpoint Runge-Kutta method. Make the same plot for a simple pendulum (**to submit by 7 November, 2021, 9pm**).
  - Make a solar system orrery following the steps outlined in the lecture (**to submit by 14 November, 2021, 9pm**)!
  - Elliptical partial differential equations: Solve the Poisson equation for the electromagnetic potential using the SOR method described in the lecture, with boundary conditions given by a 1000 Volt stick in the center of a 0 Volt box (as depicted in the lecture notes). Plot the contours of the resulting potential (**to submit by 21 November, 2021, 9pm**).
  - Interpolation, Part 1: Trace the movement of electrons in an electromagnetic potential (e.g. the one from the last exercise) with Leapfrog or Runge-Kutta. Use bilinear or bicubic interpolation for the potential. (**not graded, to submit by 28 November, 2021, 9pm**).
  - Interpolation, Part 2 (**WIN A PRIZE**): Design an optimal electron detector (specifics in lecture materials) (**to submit by 5 December, 2021, 9pm**).
  - Hyperbolic PDEs: Solve the linear advection equation by evolving an initial waveform in a periodic grid. See how the waveform behaves after passing through the grid multiple times and compare the results you get when using various methods introduced in the lecture (e.g. the LAX method, upwind scheme, LAX-Wendroff method...) (you can get 0.5 bonus points if you implement all three variants provided in Stefan's Hyperbolic Hints.txt) (**to submit by 12 December, 2021, 9pm**).
  - 2D advection: Solve the 2D advection problem using two methods introduced in the lecture (CIR and CTU) and compare if and how your solution diffuses numerically (**to submit by 19 December, 2021, 9pm**).
  - **Optional (grade for this exercise can replace your worst grade so far):** 1D Hydrodynamics: Solve the “shock tube” problem and the Sedov-Taylor blast wave using the three methods provided in the lecture (detailed assignment can be found in the lecture notes) (**to submit by 19 December, 2021, 9pm **).
stadelweb

User Tools

Site Tools

Page Tools
