ausmdv.c File Reference

Wada and Liou's flux calculator. More...

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

Defines
#define	K_SWITCH   10.0
#define	C_EFIX   0.125
Functions
int	ausmdv (struct flow_state Lft[], struct flow_state Rght[], struct flux_data_1D IF[], int first, int last)
	Compute the fluxes across an interface.

---

Detailed Description

Wada and Liou's flux calculator.

Implemented from details in their AIAA paper with hints from Ian Johnston.

 * Y. Wada and M. -S. Liou (1994)
 * A flux splitting scheme with high-resolution and 
 * robustness for discontinuities.
 * AIAA-94-0083.
 * 

Author:

P. A. Jacobs Department of Mechanical Engineering The University of Queensland

Version:

 * 26-Jun-97: Initial coding.
 * 08-Jul-97: Fix entropy fix (ironic, eh?)
 *            And finally debugged by Ian J.
 * 12-Oct-97: data structure for fluxes, array data
 * 15-Oct-97: vectorise for Cray
 * 

Todo:

Really should get rid of the vector loop some time.

---

Function Documentation

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

Compute the fluxes across an interface. Parameters: Lft[]  : IN : array Left flow states (with velocities in local frame of reference) Rght[]  : IN : array of Right flow states IF[]  : OUT : array for interface flux data first  : IN : index of first interface to compute last  : IN : index of last interface to compute * IF[i].mass : flux of mass across the interface (kg/s/m**2) * IF[i].n_mom : flux of normal momentum * IF[i].p_mom : flux of tangential momentum * IF[i].q_mom * IF[i].energy : flux of energy * IF[i].species[] : array of species fluxes * e_vib[] : vibrational energies * e_e : free-electron energy *

---