DivCleaning.cpp File Reference

#include "Carmen.h"

Include dependency graph for DivCleaning.cpp:

Functions
Vector	DivCleaning (const Cell &Cell0, const Cell &Cell1, const Cell &Cell2, int AxisNo)

Function Documentation

Vector DivCleaning	(	const Cell &	Cell0,
		const Cell &	Cell1,
		const Cell &	Cell2,
		int	AxisNo
	)

                                                                                      {

    Vector  LeftAvg(QuantityNb);
    Vector  RightAvg(QuantityNb);
    Vector  Avg(QuantityNb);
    Vector  p8wave(QuantityNb);
    real    Bdiv=0;
    real    udotB=0;
    real    psiGrad=0;
    real    dx=0;

    LeftAvg  = Cell1.average();
    RightAvg = Cell2.average();
    Avg      = Cell0.average();

    dx = Cell0.size(AxisNo);
    dx = 2*dx;
    Bdiv = (RightAvg.value(AxisNo+6) - LeftAvg.value(AxisNo+6))/dx;
    p8wave.setZero();

    if(DivClean==1){
        udotB = Avg.value(2)*Avg.value(7) + Avg.value(3)*Avg.value(8) + Avg.value(4)*Avg.value(9);

        p8wave.setValue(2, -Bdiv*Avg.value(7));
        p8wave.setValue(3, -Bdiv*Avg.value(8));
        p8wave.setValue(4, -Bdiv*Avg.value(9));
        p8wave.setValue(7, -Bdiv*Avg.value(2)/Avg.value(1));
        p8wave.setValue(8, -Bdiv*Avg.value(3)/Avg.value(1));
        p8wave.setValue(9, -Bdiv*Avg.value(4)/Avg.value(1));
        p8wave.setValue(5, -Bdiv*udotB/Avg.value(1));

        return p8wave;

    }else if(DivClean==2){

        psiGrad = (RightAvg.value(6) - LeftAvg.value(6))/dx;

        p8wave.setValue(2, -Bdiv*Avg.value(7));
        p8wave.setValue(3, -Bdiv*Avg.value(8));
        p8wave.setValue(4, -Bdiv*Avg.value(9));
        p8wave.setValue(5, -Avg.value(AxisNo+6)*psiGrad);
        p8wave.setValue(6, Avg.value(6)*exp(-(cr*ch*TimeStep/SpaceStep)));

        return p8wave;
    }else if(DivClean==3){
        if(AxisNo==3)
            p8wave.setValue(6, Avg.value(6)*exp(-(cr*ch*TimeStep/SpaceStep)));
        return p8wave;
    }else
        return p8wave;


}