FPPowerSystem1BusDiffusion.c File Reference

#include <stdlib.h>
#include <basic.h>
#include <arrayfunctions.h>
#include <physicalmodels/fppowersystem1bus.h>
#include <hypar.h>

Go to the source code of this file.

Functions
double	FPPowerSystem1BusDissipationFunction (int, int, void *, double)
int	FPPowerSystem1BusDiffusionLaplacian (double *f, double *u, int dir, void *s, double t)
int	FPPowerSystem1BusDiffusionGeneral (double *f, double *u, int dir1, int dir2, void *s, double t)

Function Documentation

double FPPowerSystem1BusDissipationFunction	(	int	,
		int	,
		void *	,
		double
	)

Definition at line 19 of file FPPowerSystem1BusFunctions.c.

{
  FPPowerSystem1Bus *params = (FPPowerSystem1Bus*) p;

  double sigma  = params->sigma;
  double omegaS = params->omegaS;
  double omegaB = params->omegaB;
  double lambda = params->lambda;
  double H      = params->H;
  double D      = params->D;

  double dissp = 0;
  if (dir1 == _YDIR_) {

    double term = (sigma*sigma*omegaS*omegaS) / (4.0*H*H);
    double expterm = exp(-t/lambda);

    if (dir2 == _XDIR_) {
      dissp = term * (lambda*omegaB) * (lambda*(1-expterm) - t*expterm);
      /* dissp = term * lambda*omegaB*lambda; */
    } else if (dir2 == _YDIR_) {
      double gamma = D*omegaS / (2.0*H);
      dissp = term * (lambda*(1-expterm) + (gamma*lambda*(t*expterm-lambda*(1-expterm))) );
      /* dissp = term * params->lambda * (1 - params->lambda*gamma); */
    }

  }

  return(dissp);
}

int FPPowerSystem1BusDiffusionLaplacian	(	double *	f,
		double *	u,
		int	dir,
		void *	s,
		double	t
	)

Definition at line 9 of file FPPowerSystem1BusDiffusion.c.

{
  HyPar             *solver = (HyPar*)             s;
  FPPowerSystem1Bus *params = (FPPowerSystem1Bus*) solver->physics;

  int *dim    = solver->dim_local;
  int ghosts  = solver->ghosts;
  int ndims   = solver->ndims;
  int nvars   = solver->nvars;

  /* calculate total size of arrays */
  int bounds[ndims]; _ArrayAddCopy1D_(dim,(2*ghosts),bounds,ndims);
  int size;          _ArrayProduct1D_(bounds,ndims,size); size *= nvars;

  /* calculate dissipation coefficient  -- constant in x and y */
  double dissipation = FPPowerSystem1BusDissipationFunction(dir,dir,params,t);

  /* f = dissipation * u */
  _ArrayScaleCopy1D_(u,dissipation,f,size);

  return(0);
}

int FPPowerSystem1BusDiffusionGeneral	(	double *	f,
		double *	u,
		int	dir1,
		int	dir2,
		void *	s,
		double	t
	)

Definition at line 32 of file FPPowerSystem1BusDiffusion.c.

{
  HyPar             *solver = (HyPar*)             s;
  FPPowerSystem1Bus *params = (FPPowerSystem1Bus*) solver->physics;

  int *dim    = solver->dim_local;
  int ghosts  = solver->ghosts;
  int ndims   = solver->ndims;
  int nvars   = solver->nvars;

  /* calculate total size of arrays */
  int bounds[ndims]; _ArrayAddCopy1D_(dim,(2*ghosts),bounds,ndims);
  int size;          _ArrayProduct1D_(bounds,ndims,size); size *= nvars;

  /* calculate dissipation coefficient  -- constant in x and y */
  double dissipation = FPPowerSystem1BusDissipationFunction(dir1,dir2,params,t);

  /* f = dissipation * u */
  _ArrayScaleCopy1D_(u,dissipation,f,size);

  return(0);
}