ausm.c File Reference

One-dimensional flux splitting according to Liou and Steffen. More...

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

Defines
#define	FIRST_ORDER   1
#define	SECOND_ORDER   2
#define	P_SPLIT   FIRST_ORDER
Functions
int	ausm (struct flow_state QL[], struct flow_state QR[], struct flow_state QIF[], double WSL[], double WSR[], int first, int last)
	Compute interface flux of mass, momentum and energy.

---

Detailed Description

One-dimensional flux splitting according to Liou and Steffen.

Todo:

Really should get rid of the vector loop some time.

---

Function Documentation

int ausm (  struct flow_state  QL[], struct flow_state  QR[], struct flow_state  QIF[], double  WSL[], double  WSR[], int  first, int  last )

Compute interface flux of mass, momentum and energy. Given the initial LEFT and RIGHT states (QL and QR) either side of an interface, compute the interface state. This solution is returned as values of the flow quantities at the interface (QIF) and the velocities of the left and right waves. Parameters: *QL  : IN : pointer to the array of LEFT states see "flux_calc.h" for a definition of the components *QR  : IN : RIGHT states first  : IN : first element to be computed last  : IN : last element to be computed QIF  : OUT : Flow state at the interface WSL  : OUT : Left wave velocity WSR  : OUT : Right wave velocity Version: 1.0, 20 Aug 91 1.1, 18-Oct-96, second-order polynomial for pressure splitting as per Liou and Steffen's J. Comput. Phys paper Author: P.A. Jacobs ICASE Mail Stop 123C NASA Langley Rearch Centre Hampton VA 23665. * References ... * M. -S. Liou and C. J. Steffen Jr. * A new flux splitting scheme. * NASA Technical Memorandum 104404 *

---