HyPar  1.0
Finite-Difference Hyperbolic-Parabolic PDE Solver on Cartesian Grids
fppowersystem1bus.h
Go to the documentation of this file.
1 
45 /*
46  dp d[mu(x,y)p] d[nu(x,y)p] d^2 p d^2 p
47  -- + ----------- + ----------- = D_yx ----- + D_yy -----
48  dt dx dy dy dx dy^2
49 
50  p(x,y): probability
51  x : angle between axis of generator and the magnetic field (theta in paper)
52  y : generator angular speed (omega in paper)
53  t : time
54 
55  mu(x,y) = omegaB * (y-omegaS)
56  nu(x,y) = (omegaS/2H)[Pm_avg - Pmax sin(x) - D(y-omegaS)]
57 
58  sigma^2 omegaS^2
59  D_yx = ---------------- lambda^2 omegaB
60  4 H^2
61 
62  sigma^2 omegaS^2 D omegaS
63  D_yy = ---------------- lambda ( 1 - lambda ---------- )
64  4 H^2 2 H
65 
66 */
67 
68 #define _FP_POWER_SYSTEM_1BUS_ "fp-power-system-1bus"
69 
70 /* define ndims and nvars for this model */
71 #undef _MODEL_NDIMS_
72 #undef _MODEL_NVARS_
73 #define _MODEL_NDIMS_ 2
74 #define _MODEL_NVARS_ 1
75 
76 /* define grid directions */
77 #define _XDIR_ 0
78 #define _YDIR_ 1
79 
80 typedef struct fp_power_system_1bus_parameters {
81 
82  /* input parameters */
83  double omegaB, omegaS, H, D, Pm_avg, Pmax, sigma, lambda;
84 
85  /* calculated, not an input */
86  double pdf_integral;
87 
89 
90 int FPPowerSystem1BusInitialize (void*,void*);
91 int FPPowerSystem1BusCleanup (void*);
int FPPowerSystem1BusCleanup(void *)
int FPPowerSystem1BusInitialize(void *, void *)