[TCS]

Hi,

I'm editing the InterFoam solver to add an additional transport equation (for a passive scalar). To construct the compressive term (equiv. in alpha) in the scalar transport equation, I create a new normal (equiv. src/transportProperties/transportProperties.H).

This is fine, except when I run the solver, I get the following error:

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 (core dumped)

I strongly believe that memory is not the root cause of the issue. My added code:

// Overall normal for c
surfaceVectorField nHatfv_forc(gradcf/( mag(gradcf) + ( 1e-8/pow(average(c.mesh().V()),1.0/3.0) ) ));

// Cell surface vector
const surfaceVectorField& Sf = mesh.Sf();

// nHat dot Sf 
const surfaceScalarField& nHat_forc = nHatfv_forc & Sf;

// Standard face-flux compression coefficient
surfaceScalarField Aphic(cCzbja.value()*mag(phi/mesh.magSf()));

// Phi r (THIS LINE)
surfaceScalarField Aphir(Aphic*nHat_forc);

// 1 - c
surfaceScalarField c2(1.0 - fvc::interpolate(c));

// Explicit flux of c
surfaceScalarField phicc(fvc::interpolate(c)*( phi + c2*Aphir ));

If I comment out the line of code labelled THIS LINE, and everything after, it's fine. If not, it's not fine...

Any help appreciated. Thanks in advance

[TCS]

Ok, I found the obvious mistake I made! In case anyone else is struggling with a similar issue (now or in the future):

Clearly nHat_forc and (in general) Sf need to change during the simulation, so using const to define them as not modifiable is a bad idea!

Replace

const surfaceVectorField& Sf = mesh.Sf();

with

surfaceVectorField Sf = mesh.Sf();

and similarly for nHat_forc.

My bad!

[alimea]

Hi foamers

I have created a new viscosity model and compiled it successfully. But when I want to test it in a case it gives me this error:

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 (core dumped)

And here is the code that is want to calculate "etaP_T" :

Foam::PTTlinearTempDep::PTTlinearTempDep
(
    const word& name,
    const volVectorField& U,
    const surfaceScalarField& phi,
    const dictionary& dict
)
:
    constitutiveEq(name, U, phi),
    tau_
    (
        IOobject
        (
            "tau" + name,
            U.time().timeName(),
            U.mesh(),
            IOobject::MUST_READ,
            IOobject::AUTO_WRITE
        ),
        U.mesh()
    ),
    lambda_(dict.lookup("lambda")),
    etaP_T
    (
        IOobject
        (
            name,
            U.time().timeName(),
            U.db(),
            IOobject::NO_READ,
            IOobject::AUTO_WRITE
        ),
        calcEtaP_T()
    )

{
 checkForStab(dict);
}

What's wrong with that?

Thanks