[ssa_cfd]

Hi all,



I created a PtrList to create multiple volScalarField in createFields.H


The code looks like this: (in createFields.H)

label N = 325;
    
    PtrList<volScalarField> kf(N);

for(int i = 0; i < N; i++)
{
word kfi ("kf" + name(i));
    kf.set
    (
        i,
        new volScalarField 
    (
            IOobject
            (
                kfi,
                runTime.timeName(),
                mesh,
                IOobject::NO_READ,
                IOobject::AUTO_WRITE
            ),
            mesh,
    dimensionedScalar(kfi, dimless, 0.0)
        )
    );
}

Now I want to access this volScalarField from chemistryModel.c.
currently my code looks like below in chemistryModel.c

word kfi ("kf" + name(i));
volScalarField tkf = mesh1.lookupObject<volScalarField>(kfi);

tkf[celli] = R.kf(pkf, Tkf, c2);

Now it stores the caluclated values only in "tkf[celli]" but doesn't change the value of "kfi" globally.



What should I do to modify the "kfi" values globally.

[ybapat]

Hello,


You have created local variable volScalarField. You need to use reference to original volScalarField then only it will be updated globally.



Regards,
-Yogesh

[ssa_cfd]

I tried like this

const volScalarField& tkf = mesh1.lookupObject<volScalarField>(kfi);
tkf[celli] = R.kf(pkf, Tkf, c2);

and I got the below error

lnInclude/chemistryModel.C: In instantiation of ‘Foam::scalar Foam::chemistryModel<CompType, ThermoType>::solve(const DeltaTType&) [with DeltaTType = Foam::Field<double>; CompType = Foam::psiChemistryModel; ThermoType = Foam::constTransport<Foam::species::thermo<Foam::hConstThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleEnthalpy> >; Foam::scalar = double]’:
lnInclude/chemistryModel.C:1079:43:   required from ‘Foam::scalar Foam::chemistryModel<CompType, ThermoType>::solve(const scalarField&) [with CompType = Foam::psiChemistryModel; ThermoType = Foam::constTransport<Foam::species::thermo<Foam::hConstThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleEnthalpy> >; Foam::scalar = double; Foam::scalarField = Foam::Field<double>]’
chemistrySolver/chemistrySolver/makeChemistrySolvers.C:91:1:   required from here
lnInclude/chemistryModel.C:929:13: error: assignment of read-only location ‘(& tkf)->Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>::<anonymous>.Foam::DimensionedField<double, Foam::volMesh>::<anonymous>.Foam::Field<double>::<anonymous>.Foam::List<double>::<anonymous>.Foam::UList<T>::operator[]<double>(celli)’

Then I tried like this

volScalarField& tkf = mesh1.lookupObject<volScalarField>(kfi);

and then I got the below error

lnInclude/chemistryModel.C: In instantiation of ‘Foam::scalar Foam::chemistryModel<CompType, ThermoType>::solve(const DeltaTType&) [with DeltaTType = Foam::Field<double>; CompType = Foam::psiChemistryModel; ThermoType = Foam::constTransport<Foam::species::thermo<Foam::hConstThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleEnthalpy> >; Foam::scalar = double]’:
lnInclude/chemistryModel.C:1079:43:   required from ‘Foam::scalar Foam::chemistryModel<CompType, ThermoType>::solve(const scalarField&) [with CompType = Foam::psiChemistryModel; ThermoType = Foam::constTransport<Foam::species::thermo<Foam::hConstThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleEnthalpy> >; Foam::scalar = double; Foam::scalarField = Foam::Field<double>]’
chemistrySolver/chemistrySolver/makeChemistrySolvers.C:91:1:   required from here
lnInclude/chemistryModel.C:917:62: error: invalid initialization of reference of type ‘Foam::volScalarField& {aka Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&}’ from expression of type ‘const Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>’

[einstein_zee]

lookupObjectRef

[ssa_cfd]

I think lookupObjectRef is available only in latest versions. I use version 2.4 because the solver was previously written in this version 2.4

[einstein_zee]

oh, that is too old. There might be a better solution but I think you can use const_cast. sth like this const_cast<volScalarField&>(tkf)[celli] ....

[ssa_cfd]

Thanks Hosein. const_cast solved the problem.!!