[yummigummi]

Hello everyone,
I'm trying to modify the viscoelaticModel to read viscosities from a library. My constructor looks as follows:


viscoelasticModel::viscoelasticModel
(
const volVectorField& U,
const surfaceScalarField& phi,
const volScalarField& fluidIndicator
)
:
IOdictionary
(
IOobject
(
"viscoelasticProperties",
U.time().constant(),
U.db(),
IOobject::MUST_READ,
IOobject::NO_WRITE
)
),
lawPtr_(viscoelasticLaw::New(word::null, U, phi, subDict("rheology"))),
mu1_(this->lookup("mu1")),
mu2_(this->lookup("mu2"))
{}




When I compile the library I get the following error:


viscoelasticModel/viscoelasticModel.C: In constructor ‘Foam::viscoelasticModel::viscoelasticModel(c onst volVectorField&, const surfaceScalarField&, const volScalarField&)’:
viscoelasticModel/viscoelasticModel.C:60:29: error: no matching function for call to ‘Foam::GeometricField<Foam::SymmTensor<double >, Foam::fvPatchField, Foam::volMesh>::GeometricField()’
mu2_(this->lookup("mu2"))



edit: I'm using foam extend 4.1

Thanks a lot!


Paul

[Jo88]

Hi Paul,



did you manage to to compile your code? I would like to do something similar and still try to figure out what's the easiest way. I created a new library where I can calculate different temperature shift factors and use them in my viscosity library (viscosityModels) as an input variable to calculate nu.



By the way, I also use foam-extend-4.1


Thanks a lot!
Jochen

[yummigummi]

Hey Jochen,


I ended up with a totally different work around where I received my solvent related viscosities from a different class. But I think the problem was that the order of variable declaration in my header file was different, than it was in my .C file. That's why there was the "no matching function to call" error.


I hope it helps you.
Paul

[Jo88]

Hey Paul,

thanks for your reply. In the meantime I also switched to a different approach and created a new class (temperatureShiftModel). In this class I will calculate the temperature shift factor and would like to pass it to the viscosity class. Although the code compiles I always get a symbolLookup-Error.

I think I misused pointers or did something generally wrong in exchanging variables between classes. Could you please share a code-example of how managed to pass a variable to another class.

Thanks in advance.

Best
Jochen

[yummigummi]

Hey Jochen,


what I did is copying the methods tau(), divTau(volVectorField& U), correct() and read() into my new class. Make sure to copy lawPtr_(viscoelasticLaw::New(word::null, U, phi, subDict("rheology"))) in your .C file of the new class aswell. By including #include "viscoelasticLaw.H" and #include "autoPtr.H" and declaring autoPtr<viscoelasticLaw> lawPtr_; in your header file it should be possible to totally bypass viscoelasticModel.
I guess, you could also do it by inheritance and give a new sub class utilities of viscoelasticModel.


Hope this helps you!
Bye
Paul