FPPowerSystemComputeDiffNumber.c File Reference

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

Go to the source code of this file.

Functions
double	FPPowerSystemDissipationFunction (int, void *, double)
double	FPPowerSystemComputeDiffNumber (void *s, void *m, double dt, double t)

Function Documentation

FPPowerSystemDissipationFunction()

double FPPowerSystemDissipationFunction	(	int	,
		void *	,
		double
	)

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);
}

FPPowerSystemComputeDiffNumber()

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

Definition at line 11 of file FPPowerSystemComputeDiffNumber.c.

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

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

  double  max_diff = 0;
  int     index[ndims];
  int done = 0; _ArraySetValue_(index,ndims,0);
  while (!done) {
    double dxinv; _GetCoordinate_(0,index[0],dim,ghosts,solver->dxinv,dxinv);
    double dyinv; _GetCoordinate_(1,index[1],dim,ghosts,solver->dxinv,dyinv);
    double dissp_x= FPPowerSystemDissipationFunction(0,params,t);
    double dissp_y= FPPowerSystemDissipationFunction(1,params,t);

    double local_diff_x = absolute(dissp_x) * dt * dxinv * dxinv;
    double local_diff_y = absolute(dissp_y) * dt * dyinv * dyinv;

    if (local_diff_x > max_diff) max_diff = local_diff_x;
    if (local_diff_y > max_diff) max_diff = local_diff_y;

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

  return(max_diff);
}