[shash]

Hi openfoamers ,

i was trying to create a parabolic velocity profile in oscillationFixedValue , as i require to generate v(y,t) where its parabolic in space and oscillating in time.i am new to openfoam ,just tried modifying the oscillatingfixedvalue code but its giving errors.can somebody suggest if my modifications r correct i hav attached my header and source.
[ATTACH][ATTACH]13684[/ATTACH][/ATTACH]

[gschaider]

Originally Posted by shash

Hi openfoamers ,

i was trying to create a parabolic velocity profile in oscillationFixedValue , as i require to generate v(y,t) where its parabolic in space and oscillating in time.i am new to openfoam ,just tried modifying the oscillatingfixedvalue code but its giving errors.can somebody suggest if my modifications r correct i hav attached my header and source.
[ATTACH][ATTACH]13684[/ATTACH][/ATTACH]

if your main interest is getting the boundary condition and not learning C++ then you might be interested in groovyBC.

Or, if you don't want to install 3rd-party software, have a look at the coded-boundary condition. There you only have to program a small C++-code-snipplet that sets the boundary condition but not the whole bunch of constructors and other methods.

The advantage of both methods is that the definition is where it belongs: in the case and not some extra library.

[shash]

hi gschaider,
i implemented the codedfixedvalue bc with:
inlet
{
type codedFixedValue;
value uniform (0 0 0); //$internalField;
redirectType awsum_bc;
code
#{
fixedValueFvPatchVectorField myPatch (*this);
forAll(myPatch, celli)
{
myPatch[celli] =(.1*(1-(this->patch().Cf()[celli](2)*this->patch().Cf()[celli](2)/.01))*(1 + sin(this->db().time().value())));
}
operator==(myPatch);
#};
}

but its giving errors, can you suggest what is wrong ,my velocity profile is v=v0(1+ sin(t)) ;v0 being parabolic f(y).

[gschaider]

Originally Posted by shash

hi gschaider,
i implemented the codedfixedvalue bc with:
inlet
{
type codedFixedValue;
value uniform (0 0 0); //$internalField;
redirectType awsum_bc;
code
#{
fixedValueFvPatchVectorField myPatch (*this);
forAll(myPatch, celli)
{
myPatch[celli] =(.1*(1-(this->patch().Cf()[celli](2)*this->patch().Cf()[celli](2)/.01))*(1 + sin(this->db().time().value())));
}
operator==(myPatch);
#};
}

but its giving errors, can you suggest what is wrong ,my velocity profile is v=v0(1+ sin(t)) ;v0 being parabolic f(y).

Sorry. My time is limited. I support "my" stuff (groovyBC) but I don't do C++-lectures

[shash]

i already solved my problem , thank you for all ur help.

[vitors]

How shash? I'm interested in your solution.

[shash]

hi vitor,

read this thread http://www.cfd-online.com/Forums/ope...condition.html , sorry for such late reply.