[Nek5000-users] Periodic Mismatch 1

[nek5000-users at lists.mcs.anl.gov]

Hi, 

I was able to solve this by removing the element face count column. I now run into the following issue shown below. I was able to find the subroutine spec_bis_conn where this error occurs. I believe this has to do with the periodic faces I have set up. My question is whether the periodic faces need to be continuous: 

For example, the geometry is similar to fluid flow over a flat plate on both sides. The fluid is periodic, but separated by the solid thickness (this is a Conj HT problem). Essentially there is a "gap" in the periodic fluid portion if you can imagine. I saw a note in the routine that says something like "ensure that all graphs are connected", not sure if this is what it means. I tried making the fluid periodic on only one side of the plate and it worked, however when I make the other side (bottom of the plate) periodic with itself I get the following error. Any thoughts? Thanks - Michael M. 

start rec_bisect: 2006 
done: 0.0% 
done: 1.0% 
done: 2.0% 
done: 3.0% 
infinite loop 
not connected 1 63 1 0 
not connected 1 32 2 0 
not connected 1 16 3 0 
not connected 1 8 4 0 
not connected 1 4 5 0 
done: 4.0% 
not connected 13 31 6 0 
Segmentation fault (core dumped) 




----- Original Message ----- 
From: nek5000-users at lists.mcs.anl.gov 
To: "Nekton User List" <nek5000-users at lists.mcs.anl.gov> 
Sent: Tuesday, June 29, 2010 11:53:10 AM GMT -06:00 Guadalajara / Mexico City / Monterrey 
Subject: Re: [Nek5000-users] Periodic Mismatch 1 


Hi, 


I am getting a 'Periodic Mismatch 1' error: 

abort: PERIODIC MISMATCH 1: 

1 3 P ie 
1 6 E je 
6 1 E ke 

I grep'd to the error in genmap.f and from what I can tell this bit of code is checking to ensure that the elements and faces I have said are periodic, match up with the other corresponding elements and faces. From the error, it seems to think that face 3 of element 1 is Periodic, and then returns this error because of course it is not. 

I have attached the rea. I am not sure what the issue is here, and why it is not recognizing face 1 as periodic (as you will see in the rea, Elem 1 face 1 is periodic with elem 6 face 3). This error occurs at the first element. The values in the error message all seem to come from the periodic condition I have specified, but it is almost as if they are not in the correct order in the rea. 

>From my understanding: " P 1 1.00 6.00000 3.00000 0.00000 0.00000 0.00000 " states that element 1 face 1 is periodic with element 6 face 3. 

Thanks for any help with this, 
Michael M. 


_______________________________________________ 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/20100629/06d08311/attachment.html>