[chusma]

Hello everyone,

I am trying to include a source term in the k-epsilon model. The source term is related with the porosity applying the Darcy-Forchheimer Law.

I have created a Field called Darcy which my program will read. With that field I calculate other two terms.

There are not compilation errors, however when I launch porousSimpleFoam the following error occurs:

new cannot satisfy memory request.
This does not necessarily mean you have run out of virtual memory.
It could be due to a stack violation caused by e.g. bad use of pointers or an out of date shared library
Aborted

Despite I have not created any new pointer to the functions the error still appearing.

Somebody could help me with this? I will appreciate any help.

Here is how I define the functions in my code:

tmp<volScalarField> kEpsilonFurow :: permeability()//New
{
volScalarField permeability_("permeability", porosity()/Darcy_);

return permeability_;
}

tmp<volScalarField> kEpsilonFurow :: porosity()//New
{




volScalarField porosity_("porosity", sqrt(1/(0.00450799*Darcy_ +1 )));

return porosity_;

}
tmp <volScalarField> kEpsilonFurow ::Gk()//New
{

volScalarField Gk("Gk", Ck_*k_*porosity()/sqrt(permeability())* sqrt(U_ & U_) );
volScalarField Zero("Zero", 0.0*k_);
if (Darcy_ > Zero)
return Gk;
else
return Zero;


}

And in the k equation:

tmp<fvScalarMatrix> kEqn
(
fvm::ddt(k_)
+ fvm::div(phi_, k_)
- fvm::laplacian(DkEff(), k_)
==
G + Gk()

- fvm::Sp(epsilon_/k_, k_)
);

[agustinvo]

Hi,


I know it's been a while, but did you find an answer to this issue? I'm facing a similar problem. I have a function which returns a tmp field, but if I create the field inside the function I get also this error message when running.