FPPowerSystemFunctions.c File Reference

#include <math.h>
#include <physicalmodels/fppowersystem.h>

Go to the source code of this file.

Functions
double	FPPowerSystemDriftFunction (int dir, void *p, double x, double y, double t)
double	FPPowerSystemDissipationFunction (int dir, void *p, double t)

Function Documentation

FPPowerSystemDriftFunction()

double FPPowerSystemDriftFunction	(	int	dir,
		void *	p,
		double	x,
		double	y,
		double	t
	)

Definition at line 4 of file FPPowerSystemFunctions.c.

{
  FPPowerSystem *params = (FPPowerSystem*) p;

  double drift = 0;
  if (dir == 0)   drift = params->O_s * (y - 1.0);
  else if (dir == 1) {
    if      (t < params->tf)    params->Pmax = params->E*params->V/params->g1;
    else if (t < params->tcl)   params->Pmax = 0.0;
    else                        params->Pmax = params->E*params->V/params->g2;
    drift = (1.0/(2.0*params->H)) * (params->Pm - params->Pmax*sin(x) - params->D*(y-1.0));
  }

  return drift;
}

FPPowerSystemDissipationFunction()

double FPPowerSystemDissipationFunction	(	int	dir,
		void *	p,
		double	t
	)

Definition at line 20 of file FPPowerSystemFunctions.c.

{
  FPPowerSystem *params = (FPPowerSystem*) p;

  double dissp = 0;
  if (dir == 1) {
    double gamma = params->D / (2.0*params->H);
    dissp =   (1.0/(2.0*params->H)) * (1.0/(2.0*params->H)) 
            * ((params->l*params->q*params->q)/(params->l*gamma+1.0))
            * (1.0 - exp(-(gamma+1.0/params->l)*t));
  }

  return(dissp);
}