[us7]

Hi foamers,

I am using Openfoam V 2.4 and trying to add a relative permeability into DarcyForchheimer Model (already available in OpenFoam src/finiteVolume/cfdTools/gernal/porosity model). I want to use this with twoPhaseEulerFoam by specifying an explicitPorositySource --> DarcyForchheimer in the fvOptions file.

-available source term for DarcyForchheimer is
Si = -(μ D + 1/2ρ|u|F)ui where D = 1/K , K= intrinsic permeability that is directly feed in /constant/fvOptions e.g. d d [2e106 -1000 -1000] m2

What i am trying to do:
-Trying to add is K--> KKri in case of two phase for relative permeability if Kri=alpha then source term will become
Si = -(μ D/alpha + 1/2ρ|u|F)ui
In short i just want to divide D with alpha and alpha value should be coming from the actual solver that is twoPhaseEulerfoam just like its getting value of U.

For this i added scalarField alpha by just follwing U in all related Files

>modified DarcyForchheimerTemplate.C by changing
part of "
(mu[cellI]*dZones[j]/alpha[cellI]) + (rho[cellI]*mag(U[cellI]))*fZones[j];

const scalar isoCd = tr(Cd);

Udiag[cellI] += V[cellI]*isoCd;
Usource[cellI] -= V[cellI]*((Cd - I*isoCd) & U[cellI]);"

Full files are attached.

When i compile this it gives me following error:


required from ‘Foam:orosityModel::addmeshConstructorToTable<po rosityModelType>::addmeshConstructorToTable(const Foam::word&) [with porosityModelType = Foam:orosityModels::fixedCoeff]’
cfdTools/general/porosityModel/fixedCoeff/fixedCoeff.C:37:9: required from here
/home/umer/OpenFOAM/OpenFOAM-2.4.0/src/OpenFOAM/lnInclude/runTimeSelectionTables.H:76:66: error: invalid new-expression of abstract class type ‘Foam:orosityModels::fixedCoeff’
return autoPtr< baseType >(new baseType##Type parList); \
^


I am only changing DarcyForchheimer but somehow it is also linked with fixedCoeff model. Anyone tell me what is this problem and how can i resolve this?

Regards,
Umer

[A CFD free user]

Hi Us7,
Have you managed to solve the problem by now?

Or in other way, is it possible to define relative permeability in interFoam with porous media?

Thanks,

[roshanak.rabiee878]

Hi
have you solved your problem?

[Major0412]

Originally Posted by roshanak.rabiee878

Hi
have you solved your problem?

Hi! roshanak
I meet the same question, have you solved your problem?
thanks! Best wishes!