1 / 45

Introduction to computational plasma physics

Introduction to computational plasma physics. 雷奕安 62755208 , yalei@pku.edu.cn. 课程概况. http://www.phy.pku.edu.cn/~fusion/forum/viewtopic.php?t=77 上机 成绩评定为期末大作业. Related disciplines. Computation fluid dynamics (CFD) Applied mathematics, PDE, ODE Computational algorithms

Download Presentation

Introduction to computational plasma physics

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript


  1. Introduction to computational plasma physics 雷奕安 62755208,yalei@pku.edu.cn

  2. 课程概况 • http://www.phy.pku.edu.cn/~fusion/forum/viewtopic.php?t=77 • 上机 • 成绩评定为期末大作业

  3. Related disciplines • Computation fluid dynamics (CFD) • Applied mathematics, PDE, ODE • Computational algorithms • Programming language, C, Fortran • Parallel programming, OpenMP, MPI • Plasma physics, space, fusion, … • Unix, Linux, …

  4. 大规模数值模拟的特殊性

  5. Contents • What is plasma • Basic properties of plasma • Plasma simulation challenges • Simulation principles

  6. What is plasma • Partially ionized gas, quasi-neutral • Widely existed • Fire, lightning, ionosphere, polar aurora • Stars, solar wind, interplanetary (stellar, galactic) medium, accretion disc, nebula • Lamps, neon signs, ozone generator, fusion energy, electric arc, laser-material interaction • Basic properties • Density, degree of ionization, temperature, conductivity, quasi-neutrality • magnetization

  7. Plasma vs gas

  8. Basic properties • Temperature • Quasi-neutrality • Thermal speed • Plasma frequency • Plasma period

  9. U→0 λD Debye length • System size and time • Debye shielding

  10. Debye lengths

  11. Plasma parameter • Strong coupling • Weak coupling

  12. Weakly coupled plasmas

  13. Collision frequency • Mean-free-path • Collisional plasma • (Collisionless) • Collisioning frequency

  14. Magnetized plasma • Anisotropic • Gyroradius • Gyrofrequency • Magnetization parameter • Plasma beta

  15. Simulation challenges • Problem size: 1014 ~ 1024 particles • Debye sphere size: 102 ~ 106 particles • Time steps: 104 ~ 106 • Point particle, computational unstable, sigularities

  16. Solution • No details, essence of the plasma • One or two dimension to reduce the size • No high frequency phenomenon, increase time step length • Reduce ND, mi / me • Smoothing particle charge, clouds • Fluidal approaches, single or double • Kinetic approaches, df/f

  17. Simple Simulation • Electrostatic 1 dimensional simulation, ES1 • Self and applied electrostatic field • Applied magnetic field • Couple with both theory and experiment, and complementing them

  18. Basic model

  19. Basic model

  20. Basic model • Field -> force -> motion -> field -> … • Field: Maxwell's equations • Force: Newton-Lorentz equations • Discretized time and space • Finite size particle • Beware of nonphysical effects

  21. Computational cycle

  22. Equation of motion • vi, pi, trajectory • Integration method, fastest and least storage • Runge-Kutta • Leap-frog

  23. Planet Problem x0 = 1; vx0 = 0; y0 = 0; vy0 = 1 read (*,*) dt N = 30/dt do i = 0, N+3 x1 = x0 + vx0*dt y1 = y0 + vy0*dt r = sqrt(x0*x0 + y0*y0) fx = -x0/r**3 fy = -y0/r**3 vx1 = vx0 + fx*dt vy1 = vy0 + fy*dt ! if(mod(i,N/10).eq.2) write(*,*) x0, y0, -1/r+(vx0*vx0+vy0*vy0)/2 x0 = x1; y0 = y1; vx0 = vx1; vy0 = vy1 enddo end Forward differencing

  24. Planet Problem ./a.out > data 0.1 $ gnuplot Gnuplot> plot “data” u 1:2

  25. Planet Problem ./a.out > data 0.01 $ gnuplot Gnuplot> plot “data” u 1:2

  26. Planet Problem x0 = 1; vx0 = 0; y0 = 0; vy0 = 1 read (*,*) dt N = 30/dt x1 = x0 + vx0*dt y1 = y0 + vy0*dt xh0 = (x0+x1)/2; yh0 = (y0+y1)/2 do i = 0, N xh1 = xh0+vx0*dt; yh1 = yh0 + vy0*dt; r = sqrt(xh0*xh0 + yh0 *yh0 ) fx = -xh1/r**3 fy = -yh1/r**3 vx1 = vx0 + fx*dt vy1 = vy0 + fy*dt ! if(mod(i,N/100).eq.0) write(*,*) xh0, yh0, -1/r+(vx0*vx0+vy0*vy0)/2 xh0 = xh1; yh0 = yh1; vx0 = vx1; vy0 = vy1 enddo end Leap Frog

  27. Planet Problem ./a.out > data 0.1 $ gnuplot Gnuplot> plot “data” u 1:2

  28. Planet Problem ./a.out > data 0.01 $ gnuplot Gnuplot> plot “data” u 1:2

  29. Field equations • Poisson’s equation

  30. Field equations • Poisson’s equation is solvable • In periodic boundary conditions, fast Fourier transform (FFT) is used, filtering the high frequency part (smoothing), is easy to calculate

  31. Particle and force weighting • Particle positions are continuous, but fields and charge density are not, interpolating • Zero-order weighting • First-order weighting, cloud-in-cell

  32. Higher order weighting • Quadratic or cubic splines, rounds of roughness, reduces noise, more computation

  33. Initial values • Number of particles and cells • Weighting method • Initial distribution and perturbation • The simplest case: perturbed cold plasma, with fixed ions. • Warm plasma, set velocities

  34. Initial values

  35. Diagnostics • Graphical snapshots of the history • x, v, r, f, E, etc. • Not all ti • For particle quantities, phase space, velocity space, density in velocity • For grid quantities, charge density, potential, electrical field, electrostatic energy distribution in k space

  36. Tests • Compare with theory and experiment, with answer known • Change nonphysical initial values (NP, NG, Dt, Dx, …) • Simple test problems

  37. Server connection SshHost: 162.105.23.110, protocol: ssh2 Your username & password Vnc connectionIn ssh shell: “vncserver”, input vnc passwd, remember xwindow number Tightvnc: 162.105.23.110:xx (the xwindow number) Kill vncserver: “vncserver –kill :xx” (x-win no.)

  38. Xes1 Xes1 document Xgrafix already compiled in /usr/local Xes1 makefile make ./xes1 -i inp/ee.inp LIBDIRS = -L/usr/local/lib -L/usr/lib -L/usr/X11R6/lib64

  39. Clients Sshputty.exe Vncviewerhttp://www.phy.pku.edu.cn/~lei/vncviewer.exe Pscp: http://www.phy.pku.edu.cn/~lei/pscp.exe

More Related