New scalar in phaseModel.c

[bestucan]

@kimy yes, declaration in .H file and definition in .C file.

declaration of function myBeta() in phaseModel.H

//- Return the maximum phase-fraction (e.g. packing limit)
        scalar alphaMax() const;
        scalar myBeta() const;

and definition of function myBeta() in phaseModel.C

Foam::scalar Foam::phaseModel::alphaMax() const
{
    return alphaMax_;
}

Foam::scalar Foam::phaseModel::myBeta() const
{
    return myBeta_;
}

[kimy]

@bestucan But why alphaMax() was originally defined in the phaseModel.H without any error?

        scalar alphaMax() const
        {
            return alphaMax_;
        }

After I defined only scalar myBeta() in phaseModel.H, the following error occurs:

In file included from phaseModel/phaseModel.C:26:0:
phaseModel/phaseModel.H:158:29: error: invalid declarator before 'd'
         tmp<volScalarField> d() const;
                             ^
phaseModel/phaseModel.H:162:26: error: non-member function 'const Foam::PhaseCompressibleTurbulenceModel<Foam::phaseModel>& Foam::turbulence()' cannot have cv-qualifier
             turbulence() const;

[bestucan]

@kimy because the position of this part

 scalar alphaMax() const
        {
            return alphaMax_;
        }

in phaseModel.H belong to the declaration of class phaseModel.
so you do have to name alphaMax() as phaseModel::alphaMax()
as for "Foam::", just for rigorous.

"declaration in H and define in C" is not something like ban. just a rule free to follow. and it's a good habit.

---

you had edit your last post and delete the important part. I've seen it. you lost a ";". and error also said "invalid declarator before 'd'"

[kimy]

@bestucan Many thanks. The same error occurs, same as original when I run my case.

Could you please upload your phaseModel.C and H files here?

[bestucan]

pasheModel_and_systemsolve.zip

systemsolve is where I call myBeta()

[kimy]

@bestucan I am a little bit crazy because my files are totally same as yours.......Still error after run. Note that I did not change the default library into user's library (All of them I leave it as defualt "LIB = $(FOAM_LIBBIN)/liblibcompressibleTwoPhaseSystem"
).

[kimy]

@bestucan

Selecting default blending method: none
Selecting dragModel for (solids in water): SchillerNaumann
Selecting swarmCorrection for (solids in water): none
#0  Foam::error::printStack(Foam::Ostream&) at ??:?
#1  Foam::sigSegv::sigHandler(int) at ??:?
#2  ? in /lib64/libc.so.6
#3  ? at ??:?
#4  __libc_start_main in /lib64/libc.so.6
#5  ? at ??:?
Segmentation fault (core dumped)

[bestucan]

@kimy do it slowly.

1 find a normal case which works well. ensure basement is good
2 add myBeta() only, and not use it. ensure compile works
3 use cout << myBeta() << "somg flag". ensure function myBeta() works.
4 use myBeta_ as a constant value in calculation. ensure framework is well.
5 use myBeta_ as variable

divide int more step if needed

then use it as your wish. special case, special solver. If something wrong. there's nothing to do with myBeta(). maybe case? maybe solver? maybe variable overflow?

you combine all steps into one step. it's hardly to debug.
have you print myBeta_ value? just same as the value in case file?

the error message seems not a special error. something about cout (Ostream&)

[kimy]

@bestucan Thanks. Your suggestions are very helpful. However, the solver met the same error mentioned before. I uploaded my solver and case files. Hope you can see them if you wish (I cannot upload). I modified the solver based on twophaseeulerfoam.

[kimy]

2020twoPhaseEulerFoam.zip

[bestucan]

I've tried. seems like your OpenFOAM version is different with mine. my version is 8

[kimy]

@bestucan Thanks a lot. What you mean is that you did not meet any problem in your OpenFOAM version? Could you please upload the whole modified solver files here (twophaseeulerfoam)? The files I uploaed above was copied from blueCFD (maybe openfoam 5). Additionally, I am also using HPC where openfoam-v1906 is loaded.

[kimy]

This is the solver file of openfoam_v1906

slurrypipebendsolver.zip

[kimy]

I have already reached my goal in another way. Thanks for all helps.

[bestucan]

coagulation. complex for me to compile your solver. downloading and compiling openfoam 5 and openfoam v1906 cost lot of time. and several source file end with ".T.C" hadn't been upload.

emmm..... not a good way to debug. we have different OS and different environment.

anyway, you've solved the problem.

[李东岳]

@bestucan Thanks for tracking this issue

[bestucan]

@东岳

[kimy]

@bestucan Thanks for your attention.