[jlorenzo]

Dear all,

I'm trying to code a solver where some boundaries are moved according to the solution, and I need the user to specify in a file NumericalParameters the names of the patches whose nodes will be displaced.

I do this with the following piece of code in createFields.H:

IOdictionary numericalParameters
(
    IOobject
    (
	"numericalParameters",
	runTime.constant(),
	mesh,
	IOobject::MUST_READ,
	IOobject::NO_WRITE
    )
);

scalar gamma(readScalar(numericalParameters.lookup("gamma")));
label incount(readLabel(numericalParameters.lookup("incount")));
wordList movingPatches(numericalParameters.lookup("movingPatches"));

Then at some point in the main .C file I have something like this

forAll(movingPatches,i)
{
    const word& patchName = movingPatches[i];
    label patchID = mesh.boundaryMesh().findPatchID(patchName);

    // Extra code
}

It seems that the code fails in

label patchID = mesh.boundaryMesh().findPatchID(patchName);

As if the wordList movingPatches didn't store the actual list that is provided by the user. In fact, I have observed that when using the command Info to show the list elements everything seems OK, but when using the command cout I get the list elements followed by \n1 (new line and 1). I have also tried movingPatches[i].value but this doesn't work.

Thank you in advance.

[Jibran]

​ Hi,

It would be useful if you could post the error as well. My guess is that the word list is not stored properly. Can you access all of the patches?

forAll(movingPatches,i)
{
    const word& patchName = movingPatches[i];
    Info << "patch name = " << patchName << endl;
}

Moreover, I would recommend implementing boundary conditions outside of your main solver just like standard OpenFOAM BCs. You can easily read fields from your solver within the new boundary condition through 'lookupPatchField' or 'lookupObject' member functions (search the forum for its usage). ​