[Tobi]

Dear all,

I want to change the boundary types of a special pointScalarField.
As I can figure out with doxygen, I need the GeometricBoundaryField for my further manipulations. Sofar so good, done:

pointScalarField::GeometricBoundaryField& gBF = pointMotionU_.boundaryField();

I also can get the types of the BC via .types(). But now the problem is that I don't be successful to assign a new type to the corresponding patch. With doxygen I can figure out that there is the assignment operator=() but till now I failed.

Any suggestions are appreciated.
Thanks in advance.


PS: To explain the problem more precisely. I want to use ACMI with moving meshes (not solidBodyMotion). Hence we create a pointScalarField called pointMotionU_. The problem now is that the patch types in that field for ACMI boundaries are of type ACMI. The problem is that these patches should become zeroGradient.

[epi_c]

Hi Tobias,

When you say operator=() did you mean:

//- Assignment to FieldField<PatchField, Type>
void operator=(const FieldField<PatchField, Type>&);

To change an individual patch's type, I think this is the best choice.

B.R.

[Tobi]

Hey,

thanks for the replay. I made a beautiful mistake. I always checked out Doxygen for 2.3.1 but I use 3.0+. I checked the online version and could not get the same classes that I have in 2.3.1. Tomorrow I will compile Doxygen for 3.0+ and recheck the stuff.

I will check your funtion and give a replay. Thanks in advance.
Tobi

[Tobi]

Dear Bryan (?), B.R.

The problem at the moment (may be I am too focused on the patch type change) is how to do exactly what you have. Here the code that I already have (nothing special):

    Info<< "Manipulating patch fields for ACMI\n" << endl;

    pointScalarField::GeometricBoundaryField& gBF =
        pointMotionU_.boundaryField();

    const wordList& typeNames = gBF.types();
    Info<< gBF << endl;
    
    forAll(typeNames, n)
    {
        if (typeNames[n] == "cyclicACMI")
        {
            //- Get patchField
            //patchField& pF = gBF[](n);

            //- Set new patchField with new patch type
            //gBF[](n)=
        }
    }

The problem that I have up to know. To get the patch field is not working because patchField is not defined. I already included the fvPatch.H but still it is not defined. Furthermore I can not figure out how to change the patch type of a given patch field.

Any suggestions are appreciated.

Thanks in advance,
Tobi

[epi_c]

Hi Tobias,

As far as I know:

typedef GeometricField<scalar, pointPatchField, pointMesh> pointScalarField;

So you need to include pointPatchFields.H (pointPatchField.H) rather than fvPatch.H.

Can you try this:

Info<< "Manipulating patch fields for ACMI\n" << endl;

pointScalarField::GeometricBoundaryField& gBF = pointMotionU_.boundaryField();
Info<< gBF << endl;

const wordList& typeNames = gBF.types();
forAll(typeNames, i)
{
	if (typeNames[i] == "cyclicACMI")
	{
		FieldField<pointPatchField, scalar>& xxx = gBF;
		PtrList<pointPatchField<scalar> >& yyy = xxx;
		const pointPatch& pp = yyy[i].patch();
		const DimensionedField<scalar, pointMesh>& pif = yyy[i].dimensionedInternalField();
		yyy.set(i, new zeroGradientPointPatchField<scalar>(pp, pif);
	}
}

I'm not familiar with pointMotion, so It's hard for me to test it.

Best Regard

[Tobi]

Hi,

thanks for the replay. I also tested out something different:

pointScalarField tmp
(
    IOobject
    (
        "tmp",
        runTime.timeName(),
        mesh,
        IOobject::NO_READ,
        IOobject::NO_WRITE
    ),
    point::New(mesh),
    dimensionedScalar("tmp", dimensionSet(0,1,-1,0,0,0,0), scalar(0)),
    newTypeNames,
    oldTypeNames
);

Info<< tmp.boundaryField() << endl;

pointMotionU_ = tmp;

Info<< pointMotionU_.boundaryField() << endl;

This I based on the constructor:

GeometricField (const IOobject &, const Mesh &, const dimensioned< Type > &, const wordList &wantedPatchTypes, const wordList &actualPatchTypes=wordList())  
Constructor given IOobject, mesh, dimensioned<Type> and patch types.

For the tmp field it was working. But the copy was not successful. I will check out your solution tomorrow. Thanks in advance. Tobi

[Tobi]

Dear all,

first, thank you for the replay. The code you mentioned is working and I get the modified patch types for the ACMI. The solution looks better but still I get interpolation somehow. Hence I did the same for the volMotionU_ field. After that the stuff looked even better but I still get interpolations in between the point field of the ACMI patches (if they get connected). I do not know why but it is based on the ACMI type. See the pictures. Therefore, I need more investigations how to manipulate the ACMI patch values. If this is done, I think the ACMI will also work fine with dynamic meshes motion (not zone motion). Maybe I will focus on the manipulation on the pointScalarFields but there is only less time doing this. I thought after changing the ACMI patch types it is done, but not as I saw, not really (I was too enthusiastic) . Maybe it would be even better to manipulate the ACMI patch class somehow to forbid interpolation for the motion fields. Thanks epi_c for the kind help. Again I learned more about programming in FOAM.

Up to now, I think I got a good insight into what we have to manipulated but it is not finished . if I get any progress I will let you know.

[epi_c]

Hi Tobias,

I'm glad the code is working and the result improved.

Best Regards