Hello Kyle,<br><br>I do not know much about the set up for your case, however I see an error in the way the .box file is setup.<br><br>In a .box file, the points should be in an ascending order. Hence for the following lines in the your .box file the points in y-direction should be in the reverse order, i.e. start with -3.0 to -1.0.<br>
<br>Box lower fuel<br>
-8 8 nelx,nely for Box<br>
0 8 1. x0,x1,gain (rescaled in usrdat)<br>
-1.00 -1.04 -1.25 -1.55 -2.00 -2.45 -2.75 -2.96 -3.00 y0,y1,...,yn<br>
, , , bc's (3 chars each!)<br>
I ,I ,I ,E bc's (3 chars each!)<br>
<br>Hope this helps.<br><br>Harish.<br><br><div class="gmail_quote">On Tue, Jun 19, 2012 at 6:16 PM, <span dir="ltr"><<a href="mailto:nek5000-users@lists.mcs.anl.gov" target="_blank">nek5000-users@lists.mcs.anl.gov</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Thanks Justin. I did look at that example, and it is helpful (especially for setting up the boundary conditions in the .box file, which I have now adjusted to be blank for the fluid bc's in the solid domain).<br>
<br>
My issue is that I don't understand the general steps necessary to setup an arbitrary domain. Simply changing the dimensions in the .box file, and using the conj_ht.rea as the base .rea file like this:<br>
<br>
conj_ht.rea<br>
2 spatial dimension ( < 0 --> generate .rea/.re2 pair)<br>
2 number of fields<br>
<br>
Box fluid<br>
-8 8 nelx,nely for Box<br>
0 8 1. x0,x1,gain (rescaled in usrdat)<br>
-1.0 -0.96 -0.75 -0.45 0.00 0.45 0.75 0.96 1.00 y0,y1,...,yn<br>
P ,P ,W ,W bc's (3 chars each!)<br>
P ,P ,E ,E bc's (3 chars each!)<br>
Box upper fuel<br>
-8 8 nelx,nely for Box<br>
0 8 1. x0,x1,gain (rescaled in usrdat)<br>
1.00 1.04 1.25 1.55 2.00 2.45 2.75 2.96 3.00 y0,y1,...,yn<br>
, , , bc's (3 chars each!)<br>
I ,I ,E ,I bc's (3 chars each!)<br>
Box lower fuel<br>
-8 8 nelx,nely for Box<br>
0 8 1. x0,x1,gain (rescaled in usrdat)<br>
-1.00 -1.04 -1.25 -1.55 -2.00 -2.45 -2.75 -2.96 -3.00 y0,y1,...,yn<br>
, , , bc's (3 chars each!)<br>
I ,I ,I ,E bc's (3 chars each!)<br>
<br>
and then using genbox, changing nelv, deleting the fluid bc's for elements greater than 64, using genmap results in the peculiar error when running nek:<br>
<br>
setup mesh topology<br>
<br>
WARNINGa: Detected non-right-handed element.<br>
Number 129 C1-4: -0.4000E-01 -0.4000E-01 -0.4000E-01 -0.4000E-01<br>
<br>
This is with the latest repo (843). Maybe I am missing some necessary changes to the .rea file?<br>
<br>
Thanks,<br>
<br>
Kyle<br>
<br>
----- Original Message -----<br>
From: <a href="mailto:nek5000-users@lists.mcs.anl.gov">nek5000-users@lists.mcs.anl.gov</a><br>
To: <a href="mailto:nek5000-users@lists.mcs.anl.gov">nek5000-users@lists.mcs.anl.gov</a><br>
Sent: Tuesday, June 19, 2012 3:17:21 PM<br>
Subject: Re: [Nek5000-users] Setting up solid elements with genbox<br>
<br>
<br>
Hi Kyle,<br>
<br>
Did you take a look at the conj_ht case in the examples folder? It has something pretty close to what you want.<br>
<br>
Justin<br>
<br>
<br>
On Tue, Jun 19, 2012 at 1:18 PM, < <a href="mailto:nek5000-users@lists.mcs.anl.gov">nek5000-users@lists.mcs.anl.gov</a> > wrote:<br>
<br>
<br>
Hi all,<br>
<br>
I am attempting to setup a simple 2D conjugate heat transfer problem using genbox. I use the example "turbChannel" as my model. I want to setup a fluid box in between two solid boxes. Therefore, my 2DturbChannel.box file appears:<br>
<br>
base.rea<br>
2 spatial dimension ( < 0 --> generate .rea/.re2 pair)<br>
2 number of fields<br>
Box<br>
-8 8 nelx,nely for Box<br>
0 8 1. x0,x1,gain (rescaled in usrdat)<br>
-1.0 -0.96 -0.75 -0.45 0.00 0.45 0.75 0.96 1.00 y0,y1,...,yn<br>
P ,P ,W ,E bc's (3 chars each!)<br>
P ,P ,W ,E bc's (3 chars each!)<br>
Box<br>
-8 8 nelx,nely for Box<br>
0 8 1. x0,x1,gain (rescaled in usrdat)<br>
1.00 1.04 1.25 1.55 2.00 2.45 2.75 2.96 3.00 y0,y1,...,yn<br>
P ,P ,E ,W bc's (3 chars each!)<br>
P ,P ,E ,W bc's (3 chars each!)<br>
Box<br>
-8 8 nelx,nely for Box<br>
0 8 1. x0,x1,gain (rescaled in usrdat)<br>
-1.00 -1.04 -1.25 -1.55 -2.00 -2.45 -2.75 -2.96 -3.00 y0,y1,...,yn<br>
P ,P ,W ,E bc's (3 chars each!)<br>
P ,P ,W ,E bc's (3 chars each!)<br>
<br>
<br>
Running genbox on this file produces a .rea file with 192 fluid elements. I first try to change "nelv" in the .rea file to 64. Running genmap on this file produces a "PERIODIC MISMATCH" error, likely due to the 192 boundary conditions for the fluid field. So I then delete all of the boundary conditions for elements greater than 64 in the fluid (leaving thermal boundary conditions for all 192 elements). Genmap now works, but running nek with this .rea/.map combination gives the error:<br>
<br>
ERROR READING INITIAL CONDITION/DRIVE FORCE DATA<br>
ABORTING IN ROUTINE RDICDF.<br>
<br>
<br>
I thought that this procedure would work, given previous discussion on the problem of conjugate heat transfer in Nek5000 ( <a href="http://lists.mcs.anl.gov/pipermail/nek5000-users/2010-February/000121.html" target="_blank">http://lists.mcs.anl.gov/pipermail/nek5000-users/2010-February/000121.html</a> ). Perhaps I generated the .map file from the incorrect form of the .rea file? My question: is there a simple way to deal with solid elements in this configuration without using prenek?<br>
<br>
Thanks,<br>
<br>
Kyle<br>
_______________________________________________<br>
Nek5000-users mailing list<br>
<a href="mailto:Nek5000-users@lists.mcs.anl.gov">Nek5000-users@lists.mcs.anl.gov</a><br>
<a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users" target="_blank">https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users</a><br>
<br>
<br>
_______________________________________________<br>
Nek5000-users mailing list<br>
<a href="mailto:Nek5000-users@lists.mcs.anl.gov">Nek5000-users@lists.mcs.anl.gov</a><br>
<a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users" target="_blank">https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users</a><br>
_______________________________________________<br>
Nek5000-users mailing list<br>
<a href="mailto:Nek5000-users@lists.mcs.anl.gov">Nek5000-users@lists.mcs.anl.gov</a><br>
<a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users" target="_blank">https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users</a><br>
</blockquote></div><br>