Flux.cpp File Reference

Computes the numerical fluxes HLL and HLLD. More...

#include "Carmen.h"

Include dependency graph for Flux.cpp:

Functions
Vector	Flux (Cell &Cell1, Cell &Cell2, Cell &Cell3, Cell &Cell4, int AxisNo)
	Returns the flux at the interface between Cell2 and Cell3. Here a 4-point space scheme is used. Cell2 and Cell3 are the first neighbours on the left and right sides. Cell1 and Cell4 are the second neighbours on the left and right sides. More...

Detailed Description

Computes the numerical fluxes HLL and HLLD.

Function Documentation

Vector Flux	(	Cell &	Cell1,
		Cell &	Cell2,
		Cell &	Cell3,
		Cell &	Cell4,
		int	AxisNo
	)

Returns the flux at the interface between Cell2 and Cell3. Here a 4-point space scheme is used. Cell2 and Cell3 are the first neighbours on the left and right sides. Cell1 and Cell4 are the second neighbours on the left and right sides.

Parameters

Cell1	second neighbour on the left side
Cell2	first neighbour on the left side
Cell3	first neighbour on the right side
Cell4	second neighbour on the right side
AxisNo	Axis of interest

Returns

Vector

{
    // --- Local variables ---

    Vector Result(QuantityNb);

    Cell C1, C2, C3, C4;

    int BoundaryCell1 = BoundaryRegion(Cell1.center());
    int BoundaryCell2 = BoundaryRegion(Cell2.center());
    int BoundaryCell3 = BoundaryRegion(Cell3.center());
    int BoundaryCell4 = BoundaryRegion(Cell4.center());

    bool UseBoundaryCells = (UseBoundaryRegions && (BoundaryCell1!=0 || BoundaryCell2!=0 || BoundaryCell3!=0 || BoundaryCell4!=0));
    
    // --- Take into account boundary conditions ---

            if (UseBoundaryCells)
                GetBoundaryCells(Cell1, Cell2, Cell3, Cell4, C1, C2, C3, C4, AxisNo);

    switch(SchemeNb)
    {
        case 1:
        default:
            if (UseBoundaryCells)
                Result = SchemeHLL(C1, C2, C3, C4, AxisNo);
            else
                Result = SchemeHLL(Cell1, Cell2, Cell3, Cell4, AxisNo);
            break;

        case 2:
            if (UseBoundaryCells)
                Result = SchemeHLLD(C1, C2, C3, C4, AxisNo);
            else
                Result = SchemeHLLD(Cell1, Cell2, Cell3, Cell4, AxisNo);
            break;
        break;

    }


    return Result;
}

Here is the caller graph for this function: