[arnolnic]

Hi everyone,

I am trying to implement a radial gravitation to my solver, which is based on the buoyantBussinesqPimpleFoam solver.

Therefore I added a volVectorField "Ge" to createField.H as follows

PHP Code:

    //Definition vom Ursprungsvektor
    vector origin (0,0,0);

   //Definition der Konstante gc
    scalar gc (10);  


    //Definition of radial vectorfield Ge
    volVectorField Ge
    (
        IOobject
        (
       "Ge",
        runTime.timeName(),
        mesh,
        IOobject::NO_READ,
        IOobject::AUTO_WRITE
        ),
        mesh,
        dimensionedVector("Ge", dimAcceleration, vector::zero)
    );

    //Definition für Feldschleife
    forAll(mesh.C(),i)
    {
    //Info<<mesh.V()[i]<<endl;
    vector dif (origin-mesh.C()[i]);                    //
    Ge[i]=(pow(mag(origin-mesh.C()[i]),-1))*dif*gc;     //unit vector*gc
    } 


This part works just fine, as the attached picture proves.
Afterwards I simply added a volScalarField "Geh" and a surfaceScalarField "Gehf" like it's done with the normal gravity g:

PHP Code:

    //Definition vom neuen Feld Geh & Gehf
    Info<< "Calculating field g.h\n" << endl;
    volScalarField gh("gh", g & mesh.C());
    surfaceScalarField ghf("ghf", g & mesh.Cf());
    volScalarField Geh("Geh", Ge & mesh.C());
    surfaceScalarField Gehf("Gehf", g & mesh.Cf());

    volScalarField p
    (
        IOobject
        (
            "p",
            runTime.timeName(),
            mesh,
            IOobject::NO_READ,
            IOobject::AUTO_WRITE
        ),
        p_rgh + rhok*(gh+Geh)
    ); 


Unfortunately it's not that simple... I can't substitute Ge for g in the definition of the surfaceScalarField Gehf or otherwise an error occurs while compiling.

PHP Code:

error: no matching constructor for initialization of
      'surfaceScalarField' (aka 'GeometricField<scalar, fvsPatchField,
      Foam::surfaceMesh>')
    surfaceScalarField Gehf("Gehf", Ge & mesh.Cf()); 


I guess that's the cracking point why the "new gravity Field" has no influence on the convective heat transfer inside my geometry!

I appreciate any kind of help.
Thanks in advance