SHARE
TWEET

Untitled

a guest Jun 13th, 2019 47 Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
  1. f = 1.3 10^9; (*Hz*)
  2. [Omega] = 2. Pi f;
  3. Trf = 1./f;
  4.  
  5. Tc = 9.25; (*K*)
  6. T = 1.5; (*K*)
  7.  
  8. [Mu]0 = 4. Pi 10^-7; (*H/m*)
  9. [Phi]0 = 2.07 10^-15; (*Wb*)
  10. n = 5.56 10^28; (*m^-3*)
  11. e = 1.6 10^-19; (*C*)
  12. m = 9.1 10^-31; (*kg*)
  13.  
  14. [Tau] = 0.5 10^-9; (*s*)
  15. vf = 1.37 10^6; (*m/s*)
  16. Bc2 = 410. 10^-3; (*T*)
  17. l = 100 10^-9; (*m*)
  18. [Lambda] = 49.5 10^-9; (*m*)
  19. [Xi] = 28 10^-9; (*m*)
  20. [Kappa] = [Lambda]/[Xi];
  21.  
  22. [Gamma] = ([Phi]0)/([Mu]0 [Lambda] );
  23.  
  24. v0 = Sqrt[(l vf Sqrt[14 Zeta[3] (1 - T/Tc)])/(3 Pi [Tau])];
  25.  
  26. [Rho] = (m vf)/(n e^2 l);
  27. [Eta] = ([Phi]0 Bc2)/[Rho];
  28.  
  29. g = Log[[Kappa]]+0.5+Exp[-0.4-0.8 Log[[Kappa]]-0.1(Log[[Kappa]])^2];
  30. [Epsilon] = ([Phi]0^2 g)/(4. Pi [Mu]0 [Lambda]^2);
  31.  
  32. Tmax = 3. Trf;
  33. Zmax = 1. 10^-6; (*m*)
  34.  
  35. B = 100. 10^-3; (*T*)
  36.  
  37. sol = NDSolve[{
  38.      (v0^2 [Eta])/(v0^2 + D[x[t, z], t]^2)D[x[t, z],t] == [Epsilon] D[x[t, z], z, z] + [Gamma] B Cos[[Omega] t] Exp[-(z/[Lambda])],
  39.      x[0., z] == 0.,
  40.      x[t, Zmax] == 0.,
  41.      (D[x[t, z], z] /. z -> 0) == 0.
  42.      }, x, {t, 0., Tmax}, {z, 0., Zmax},
  43.     AccuracyGoal -> 13,
  44.     PrecisionGoal -> 2,
  45.     MaxSteps -> Infinity];
  46. u[t_, z_] = Evaluate[x[t, z] /. sol][[1]];
  47. du[t_, z_] = D[u[t, z],t];
  48.  
  49. phaseSpace =
  50.    ParallelTable[{Re[u[t, 0.]] 10^6, Re[du[t, 0.]] 10^-3}, {t, Trf,
  51.      Tmax, Tmax/500.}];
  52. ListPlot[phaseSpace, Joined -> True, PlotRange -> All]
RAW Paste Data
We use cookies for various purposes including analytics. By continuing to use Pastebin, you agree to our use of cookies as described in the Cookies Policy. OK, I Understand
 
Top