[jp279]

Hello everyone,
I am trying to develop a new boundary condition which is a modified version of a fixed BC. Here, I also read a scalarField H as follows :

top
{
type newBC;
value uniform 1;
H uniform 0.0; // scalarField
}

After reading this scalarField H, I update it (by adding some constant to it) for each time step in the newBC.C file as follows:

H += someConstantScalarField; ... (1)

And then I use this H to compute the value at this boundary.

The problem is that, since I am using a PIMPLE solver, if I use (1) it updates it for all the iterations of the PIMPLE loop of that particular time step. In other words, I want to perform (1) just once for a time step and use that updated H for all the PIMPLE iterations within this time step. Something like -

// if(new time iteration?) OR if(last innermost iteration of pimple?)
{
H += someConstantScalarField; ... (1)
}

How can I access the "pimple" object (defined in the solver) from within the updateCoeffs() member function of a boundary condition?

Kindly, let me know if you have any suggestions.

[syavash]

Originally Posted by jp279

Hello everyone,
I am trying to develop a new boundary condition which is a modified version of a fixed BC. Here, I also read a scalarField H as follows :

top
{
type newBC;
value uniform 1;
H uniform 0.0; // scalarField
}

After reading this scalarField H, I update it (by adding some constant to it) for each time step in the newBC.C file as follows:

H += someConstantScalarField; ... (1)

And then I use this H to compute the value at this boundary.

The problem is that, since I am using a PIMPLE solver, if I use (1) it updates it for all the iterations of the PIMPLE loop of that particular time step. In other words, I want to perform (1) just once for a time step and use that updated H for all the PIMPLE iterations within this time step. Something like -

// if(new time iteration?) OR if(last innermost iteration of pimple?)
{
H += someConstantScalarField; ... (1)
}

How can I access the "pimple" object (defined in the solver) from within the updateCoeffs() member function of a boundary condition?

Kindly, let me know if you have any suggestions.

Hi,

It's a pretty old post, however I have faced the same issue. How can I access pimple control in my own boundary condition?

Kind regards,
syavash