[Nek5000-users] Periodic BC on non-parallel boundaries
nek5000-users at lists.mcs.anl.gov
nek5000-users at lists.mcs.anl.gov
Tue Oct 11 12:51:31 CDT 2011
Hi Stefan,
I`ve got another question related to the periodic/cyclic BC`s. In addition
to the partial 2D annulus, I`m also trying to simulate a 3D annular cylinder
(a full 2π radians for now). I made the domain as I made the full 2D
annulus, by remapping a periodic-side domain in userdat2. This worked in
2D, but in 3D the periodic BC are not working properly. Is this a case
where cyclic BC would be needed, rather than periodic? I`m hoping not since
I have stress-free walls, so I`m using Pn/Pn-2 again. (I tried the same
setup with no-slip BC in the Pn/Pn formulation, and that did work properly.)
Thanks,
David
P.S. When I tried the 2D partial-annulus case with Pn/Pn, I accidentally had
IFCYCLIC=F one time, and it still worked, so I can`t tell whether the cyclic
feature is having an effect at all.
On Mon, Oct 10, 2011 at 2:20 PM, <nek5000-users at lists.mcs.anl.gov> wrote:
> Hmm, just had a look at the code. It looks to me that it might work
> out of the box. Why don't you give it a shot.
>
> Cheers
> Sefan
>
> On 10/10/11, S K <stgeke at googlemail.com> wrote:
> > Let me check what we need to change to make it happen.
> > -Stefan
> >
> > On 10/10/11, nek5000-users at lists.mcs.anl.gov
> > <nek5000-users at lists.mcs.anl.gov> wrote:
> >> Well then it seems I am stuck. Might a work-around be possible for a
> 1/4
> >> or
> >> 1/2 annulus? For instance with 1/2, I would just need to match (u,v) to
> >> (-u,-v).
> >>
> >> David
> >>
> >> On Mon, Oct 10, 2011 at 1:41 PM, <nek5000-users at lists.mcs.anl.gov>
> wrote:
> >>
> >>> btw: I am afraid, the stress formulation doesn't support cyclic BCs at
> >>> the moment!
> >>> -Stefan
> >>>
> >>> On 10/10/11, S K <stgeke at googlemail.com> wrote:
> >>> > You are right, it's an undocumented AND experimental feature. Just
> >>> > flag the boundaries as periodic. Yes, just add another IFCYCLIC line.
> >>> > Currently, there is no stress formulation support for Pn/Pn.
> >>> >
> >>> > -Stefan
> >>> >
> >>> > On 10/10/11, nek5000-users at lists.mcs.anl.gov
> >>> > <nek5000-users at lists.mcs.anl.gov> wrote:
> >>> >> Hi Stephan,
> >>> >>
> >>> >> I see, thanks. I can`t find any documentation about cyclic BC`s.
> >>> >> What
> >>> >> is
> >>> >> their symbol for the box file? Also, I don`t see a logical switch
> >>> >> for
> >>> >> IFCYCLIC. Do I just add another line?
> >>> >>
> >>> >> Also, I have stress-free BC on curved boundaries. A series of error
> >>> >> messages led me to believe that I had to use the stress formulation,
> >>> >> and
> >>> >> that this required me to use Pn/Pn-2. Is this not the case?
> >>> >>
> >>> >> Best,
> >>> >>
> >>> >> David
> >>> >>
> >>> >
> >>> _______________________________________________
> >>> Nek5000-users mailing list
> >>> Nek5000-users at lists.mcs.anl.gov
> >>> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
> >>>
> >>
> >
> _______________________________________________
> Nek5000-users mailing list
> Nek5000-users at lists.mcs.anl.gov
> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20111011/49abdcec/attachment.html>
More information about the Nek5000-users
mailing list