efm.c File Reference

Mike Macrossan's Equilibrium Flux Method via Paul Petrie. More...

#include <stdio.h>
#include <math.h>
#include "../../util/source/compiler.h"
#include "../../util/source/useful.h"
#include "flux_calc.h"

Defines
#define	PHI   1.0
#define	dtwspi   0.282094792
#define	P   0.327591100
#define	A1   0.254829592
#define	A2   -0.284496736
#define	A3   1.421413741
#define	A4   -1.453152027
#define	A5   1.061405429
#define	LIMIT   5.0
#define	EXLIM   0.138879e-10
#define	EFLIM   1.0
#define	DSIGN(val, sgn)   ( (sgn >= 0.0)? fabs(val): -fabs(val) )
Functions
int	efmflx (struct flow_state Lft[], struct flow_state Rght[], struct flux_data_1D IF[], int first, int last)
	Compute the fluxes across an interface using the Equilibrium Flux Method of Macrossan & Pullin.
int	exxef (double sn, double *ex, double *ef)
	Compute exp(-x**2) and erf(x) with a polynomial approximation.

---

Detailed Description

Mike Macrossan's Equilibrium Flux Method via Paul Petrie.

Author:

Paul Petrie

PJ

Version:

02-May-96: Code obtained from Paul Petrie and adapted for mb_cns.

03-Nov-96: Update species mass flux code. Have to fudge the gas constants.

26-Jun-97: density weighting of gas constants

12-Oct-97: flux_data_1D used to return the flux data array data

---

Function Documentation

int efmflx (  struct flow_state  Lft[], struct flow_state  Rght[], struct flux_data_1D  IF[], int  first, int  last )

Compute the fluxes across an interface using the Equilibrium Flux Method of Macrossan & Pullin. Parameters: Lft[]  : IN : array of Left flow states (with velocities in local frame of reference) Rght[]  : IN : array of Right flow state IF[]  : IN/OUT : array of interface flux data structures first  : IN : index of first interface to compute last  : IN : index of last interface to compute * interface data -- contains... * flux of mass across the interface (kg/s/m**2) * flux of normal momentum * flux of tangential momentum * flux of energy * array of species fluxes * vibrational energies * free-electron energy *

int exxef (  double  sn, double *  ex, double *  ef )

Compute exp(-x**2) and erf(x) with a polynomial approximation. Parameters: sn  : IN : x *ex  : OUT : exp(x**2) *ef  : OUT : erf(x) error function

---