[hjasak]

Dear All,

I have recently implemented a cyclic jump boundary condition during work with a client. This allows you to have a cyclic geometry and specify a jump between the patches across it from the "left" to the "right" side.

In order to parallelise this properly, I need to go through some pain and re-organisation (the cyclic condition itself was trivial). So, before I do that, I'd like to hear a word from the Forum - sounds very Roman

Do you think this is useful/do you need it or intend to use it?

In case of decent response, I'll spend 2-3 days to wrap this up properly and sort out the parallelisation, domain decomposition etc.

Looking forward to hearing from you all,

Hrv

[stephan]

hi,

great! i am certainly interested and would appreciate it a lot.
would it be posibble to define the "jump" as a function of calculated fields?
thanks in advance for your work!
best regards
stephan

[alberto]

Add a vote from me :-)

Thanks!
Alberto

[olesen]

This could be useful for simulating porous baffles. It'd be interesting to see how much the solution oscillates with higher resistances.

[gschaider]

Hi Hrv!

Yes. That would be very interesting. In fact I'm working on something similar (porous jump) for a customer of ours. Would be interested how you did it (==how it should be done)

Bernhard

[nishant_hull]

I need it for some repeating structure for LES simulation.
take my vote for it as well.

Nishant

[suredross]

i will need it for my thesis work as well.i'm all up for it!

davey

[suredross]

Hi Jasak,
are you done with the cyclic jump boundary condition?if yes,where can i access it?i am in dire need of some help!

thanx
davey

[mou_mi]

Hi

I am intersted to test the jump cyclic boundary condition in OF1.5 and I did a test. here is what I did;
I made a block by blockMesh:
patches
(
wall maxY
(
(3 7 6 2)
)
wall minY
(
(1 5 4 0)
)
cyclic X
(
(2 6 5 1)
(0 4 7 3)
)
wall Z
(
(0 3 2 1)
(4 5 6 7)
)
);

Then I changed the boundary file and adjust the initial files like:
--------U
maxY
{
type fixedValue;
value uniform (0 0 0);
}
minY
{
type fixedValue;
value uniform (0 0 0);
}
X
{
type jumpCyclic;
value uniform (0 0 0);
}
Z
{
type fixedValue;
value uniform (0 0 0);
}

------------p
maxY
{
type zeroGradient;
}
minY
{
type zeroGradient;
}
X
{
type jumpCyclic;
value uniform 0;
}
Z
{
type zeroGradient;
}

then here if the run with "icoFoam" solver:

Create time

Create mesh for time = 0

Reading transportProperties

Reading field p

Reading field U

Reading/calculating face flux field phi


Starting time loop

Time = 0.02

Courant Number mean: 0.009308026 max: 0.02843602



gradientInternalCoeffs cannot be called for a genericFvPatchField (actual type jumpCyclic)
on patch X of field U in file "OpenFOAM/mou-1.5/run/tutorials/test1/0/U"
You are probably trying to solve for a field with a generic boundary condition.

From function genericFvPatchField<type>::gradientInternalCoeffs( ) const
in file fields/fvPatchFields/basic/generic/genericFvPatchField.C at line 692.

FOAM exiting



I did not have any problem for "cyclic" BC, but with the "jumpCyclic" it give me this error, Would you help if you have any suggestion about that?

Thank you
mou

[jens_klostermann]

Hi Hrv,

I just wonder if the jumpCyclic got checked in into svn OF-1.4.1-dev, because I didn't find it yet?! Or is the one you were talking at the beginning of the thread the fanBC in OF-1.5?

Best regards Jens

[leejc]

Hi Dr. Jasak,

I was able to run dieselFoam with jumpcyclic but I could not figure out how to specify the value of the jump. Say, I like to have delta(U) = 1 along the normal of the surface, how do I specify this in the boundary condition file?

Thank you.

Jerry