[chandra shekhar pant]

Dear Foamers,
For a couple of weeks I am trying to write a time dependent boundary condition for pressure, after looking multiple post in the CFD Online, I found that codeStream boundary condition could do that (also swak4FOAM could do that, but unfortunately I am unable to compile that). This code runs and compiles well with OF v 1906/2112 but it's actually not been able to calculate the time variable, instead it is treating "t" as variable with value 0.

source
{
type fixedValue;
value #codeStream
{
codeInclude
#{
#include "fvCFD.H"
#};

codeOptions
#{
-I$(LIB_SRC)/finiteVolume/lnInclude \
-I$(LIB_SRC)/meshTools/lnInclude
#};

codeLibs
#{
-lmeshTools \
-lfiniteVolume
#};

code
#{
const IOdictionary& d = static_cast<const IOdictionary&>
(
dict.parent().parent()
);
const fvMesh& mesh = refCast<const fvMesh>(d.db());
const label id = mesh.boundary().findPatchID("source");
const fvPatch& patch = mesh.boundary()[id];
scalarField p(patch.size(), scalar(0));

const scalar t = d.db().time().value();
scalar t2 = t-tp*(floor(t/tp)); //
const scalar tp = 0.035;
const scalar tf = 0.0116;
{
forAll(p, i)
{
{
if (t2<tf)
{
p[i] = 1000000.0*t;
}
if (t2>=tf && t2<=tp)
{
p[i] = 1.0*(t-0.01);
}
}
}
}

p.writeEntry("", os);


#};
};

 }

Any comment/suggestion would be a great help, thanks in advance.

[chandra shekhar pant]

I am still struggling with this issue, could any one please help.

[chandra shekhar pant]

Hello again,

After struggling for a very quite long time, finally I managed to get the things working using the instructions given in the link http://www.wolfdynamics.com/wiki/pro...streamINIT.pdf



Previously maybe I was messing up with the initial condition and boundary condition. I have done the following to make it work (could be useful to someone):

 source
    {
      type             codedFixedValue;
      value             uniform 0;
      name            inletProfile2;
       
    code
    #{

      const fvPatch& boundaryPatch = patch();
        const vectorField& Cf = boundaryPatch.Cf();
        scalarField& p = *this;
        p = patchInternalField();
        scalar t = this->db().time().value();
      const scalar tf = 0.5;
    const scalar tp = 1.0;

  forAll(p, i)
        {
        
        {
            if (t<=tf)
            {
                    p[i] = 1000000.0*t;
                    
                    
            }
            if (t>=tf && t<=tp)
            {
                p[i] = 1.0*(t-0.01);    
                
            }
            
            }   
        
        }
    #}; 
   }