[vivek05]

Hi Everyone,
I am trying to impose a particular condition for source term in fvScalarMatrix Eqn using if - else condition with the help of forAll loop to initialize in all computational cells. Below I mentioned the condition. Now the problem is after 18th time steps the code got crashed down. Log file also has been attached with this thread.

PHP Code:

const volScalarField turb_nut = turbulence->nut(); 
const volScalarField turb_k = turbulence->k(); 

    Info << "Source term calculation starting" <<endl;
        forAll(mesh.cells(),celli)
        {
        if(((ALPHA[counter][celli])*(1.0-(ALPHA[counter][celli]))) > 0.001)
        {
    SourceElsa[celli] = ((CTurbElsaSigmaPrime*(turb_k[celli]/(ElsaTinyJap+turb_nut[celli]))*tmpElsaSigma[celli]*(1-(tmpElsaSigma[celli]/(tmpElsaSigmaEquil[celli]+sigmaTiny)))/(1 - runTime.deltaT().value()*CTurbElsaSigmaPrime*(turb_k[celli])/(ElsaTinyJap+turb_nut[celli]))*(1-2*((tmpElsaSigma[celli]/(tmpElsaSigmaEquil[celli]+sigmaTiny))))));
        }
        else
        {
        SourceElsa[celli] = SourceZero[celli];    
        }
        } 


Could someone figure out where did I make a mistake here??

[agustinvo]

Hi


coud you post the equations you want to implement? Did you check in the fvOptions if there is something you could reuse, or in other solvers?

[vivek05]

Originally Posted by agustinvo

Hi


coud you post the equations you want to implement? Did you check in the fvOptions if there is something you could reuse, or in other solvers?

Hi Augustin Villa,
I am very much obliged to hear from you!!

PHP Code:

 fvScalarMatrix ElsaSigmaPrimeEqn
        (
            fvm::ddt(SIGMA[counter])
            + fvm::div(phi,SIGMA[counter],sigmaScheme)
            - fvm::laplacian
            (
            Dab + alphatab*turbulence->nut(), SIGMA[counter],
            "laplacian(Dab,ElsaSigmaPrime)"
            )
            -Foam::sqrt((2/qWiener)*alphatab*turbulence->nut())*(fvc::grad(SIGMA[counter]) & dWdt[counter])
            ==  SourceElsa
        ); 


I want to calculate the source term in RHS side mentioned in the above equation. Now I am calculating source term in each cell using forAll loop if X(1-X)>0.001 and zero value for source term in the remaining portion of the domain. Whether is this correct way of assigning a source to all cells??

Best Wishes,
Vivek