ramrunner / mhd-lisp Goto Github PK

MHD simulation by DsP v.03-alpha

spherical front simulation:
  Currently implemented algorithm :
  a) at first we create a 2d vector, we can scale it by parameter a
  b) we use it to create a set of parallel vectors with different magnitudes
  c) we rotate the whole set in order to create a disk (we also can do it by rotation matrix)
  d) now we project stereographicly to create a spherical front

On the Sagdeev 1d plasma simulation technique:
  We implement a safe Newton Raphson method to calculate the zeros of the S(phi, M) potential
  then using a discrete scheme based on the second derivative of phi with respect to X we
  calculate Phi as function of X

  Todo:
    Use the Phi(X) functions as initials conditions on a multigrid solver to test the stability
    of those pulses.

Install: 
get a scheme distro for your system (pref chicken scheme).
have gnuplot installed to view data output.
will soon provide a makefile ;)