[DuMuX] Source Term

[Gerardo Zegers R]

Hi,

Iam doing a modeling in dumux. I try to include a extraction well but it
dosnt work.

I try two diferent metods ( in the file .hh)

First using source AtPos

void sourceAtPos(PrimaryVariables &values,
                const GlobalPosition &globalPos) const
{
 if (globalPos[1] < 20  && globalPos[1] > 10 && globalPos[0] > 5 &&
globalPos[0] < 10)
        values = -1.0;
 else
 values= 0.0;
    }

Second using source

 void source(PrimaryVariables &values,
                const Element &element,
                const FVElementGeometry &fvGeometry,
                int scvIdx) const
    {
GlobalPosition globalPos = element.geometry().corner(scvIdx);
 if (globalPos[1] < 20  && globalPos[1] > 10 && globalPos[0] > 5 &&
globalPos[0] < 10)
        values[Indices::contiWEqIdx] = -1.0;
values[Indices::contiNEqIdx]= 0.0;
else
values[Indices::contiWEqIdx] = 0.0;
values[Indices::contiNEqIdx]= 0.0;
    }


Could you help me with this?

Thanks

Regards

Gerardo Zegers
-- 
Gerardo Zegers R
89850305

_______________________________________________
Dumux mailing list
Dumux@listserv.uni-stuttgart.de
https://listserv.uni-stuttgart.de/mailman/listinfo/dumux