FPPowerSystem1BusComputeCFL.c File Reference

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

Go to the source code of this file.

Functions
double	FPPowerSystem1BusDriftFunction (int, void *, double, double, double)
double	FPPowerSystem1BusComputeCFL (void *s, void *m, double dt, double t)

Function Documentation

FPPowerSystem1BusDriftFunction()

double FPPowerSystem1BusDriftFunction	(	int	,
		void *	,
		double	,
		double	,
		double
	)

Definition at line 4 of file FPPowerSystem1BusFunctions.c.

{
  FPPowerSystem1Bus *params = (FPPowerSystem1Bus*) p;

  double drift = 0;
  if (dir == _XDIR_) {
    drift = params->omegaB * (y - params->omegaS);
  } else if (dir == _YDIR_) {
    drift = (params->omegaS/(2*params->H)) 
          * (params->Pm_avg - params->Pmax*sin(x) - params->D*(y-params->omegaS));
  }

  return drift;
}

FPPowerSystem1BusComputeCFL()

double FPPowerSystem1BusComputeCFL	(	void *	s,
		void *	m,
		double	dt,
		double	t
	)

Definition at line 11 of file FPPowerSystem1BusComputeCFL.c.

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

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

  double  max_cfl = 0;
  int     index[ndims];
  int done = 0; _ArraySetValue_(index,ndims,0);
  while (!done) {
    double x;     _GetCoordinate_(_XDIR_,index[_XDIR_],dim,ghosts,solver->x,x);
    double y;     _GetCoordinate_(_YDIR_,index[_YDIR_],dim,ghosts,solver->x,y);
    double dxinv; _GetCoordinate_(_XDIR_,index[_XDIR_],dim,ghosts,solver->dxinv,dxinv);
    double dyinv; _GetCoordinate_(_YDIR_,index[_YDIR_],dim,ghosts,solver->dxinv,dyinv);
    double drift_x= FPPowerSystem1BusDriftFunction(_XDIR_,params,x,y,t);
    double drift_y= FPPowerSystem1BusDriftFunction(_YDIR_,params,x,y,t);

    double local_cfl_x = absolute(drift_x) * dt * dxinv;
    double local_cfl_y = absolute(drift_y) * dt * dyinv;

    if (local_cfl_x > max_cfl) max_cfl = local_cfl_x;
    if (local_cfl_y > max_cfl) max_cfl = local_cfl_y;

    _ArrayIncrementIndex_(ndims,dim,index,done);
  }

  return(max_cfl);
}