[Nek5000-users] Setting up solid elements with genbox

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

Hi all,

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:

base.rea
2                      spatial dimension  ( < 0 --> generate .rea/.re2 pair)
2                      number of fields
Box
-8   8                                           nelx,nely for Box
0 8 1.                                              x0,x1,gain  (rescaled in usrdat)
-1.0 -0.96 -0.75 -0.45 0.00 0.45 0.75 0.96 1.00     y0,y1,...,yn
P  ,P  ,W  ,E                                 bc's  (3 chars each!)
P  ,P  ,W  ,E                                 bc's  (3 chars each!)
Box
-8   8                                           nelx,nely for Box
0 8 1.                                              x0,x1,gain  (rescaled in usrdat)
1.00 1.04 1.25 1.55 2.00 2.45 2.75 2.96 3.00     y0,y1,...,yn
P  ,P  ,E  ,W                                 bc's  (3 chars each!)
P  ,P  ,E  ,W                                 bc's  (3 chars each!)
Box
-8   8                                           nelx,nely for Box
0 8 1.                                              x0,x1,gain  (rescaled in usrdat)
-1.00 -1.04 -1.25 -1.55 -2.00 -2.45 -2.75 -2.96 -3.00     y0,y1,...,yn
P  ,P  ,W  ,E                                 bc's  (3 chars each!)
P  ,P  ,W  ,E                                 bc's  (3 chars each!)


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:

  ERROR READING INITIAL CONDITION/DRIVE FORCE DATA
  ABORTING IN ROUTINE RDICDF.


I thought that this procedure would work, given previous discussion on the problem of conjugate heat transfer in Nek5000 (http://lists.mcs.anl.gov/pipermail/nek5000-users/2010-February/000121.html). 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?

Thanks,

Kyle