[Nek5000-users] Conjugate Heat Transfer

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

Dear Dr. Fischer,

thank you for the reply.
We are currently testing the following procedure that would not make it
necessary to change prenek:
1) generate a mesh in prenek with all elements that will be solids
2) generate a mesh in prenek (matching the coordinate system of 1) with only the
fluid elements
3) merge all elements from 1) into 2), which should ensure that the solid
elements come behind the fluid ones in the order in the rea file. Then we would
only need to change the number of elements in the header lines of the element
section, correct?

Michael Meador, another student in Dr. Duggleby's group, is currently
investigating how to set up a basic conjugate heat transfer experiment (hot jet
impinging on conducting plate with finite thickness). We would like to compare
that data against a similar problem to be computed in nek5000. In case we'll
get something worth sharing, it'll be posted.

Thanks,
Markus



Quoting nek5000-users at lists.mcs.anl.gov:

>
> Hi Markus,
>
> We're currently not using the group distinction to segragate fluid/solids.
> (Though perhaps we should...)
>
> It is simply imposed by the fact that the first nelv elements must be
> fluid, and that the sum total must be thermal (fluid+solid).
>
> Thus, nelv \le nelt in the .rea file.
>
> A further consequence is that nek will look for only nelv*2*ndim bcs
> for the fluid.
>
> At present, prenek is not set up to supply the correct number of bcs -
> it requires intervention by another code or editing by hand.   This is
> all subject to change, but the truth is we just resurrected conjugate
> heat transfer for parallel processing in the past 12 months and, so far,
> have been the only customers.   I believe, however, that if a .rea file
> is set up this way, that genmap and nek5000 will work properly.  (I've
> used it both for 2D and 3D runs and all seems to be working ---- I never
> trust these things, however, until I've done it for at least 10
> different cases...).
>
> We can do some work on prenek to get it into shape for you, but I'm
> presently not very happy with the interface and need to kick around
> a few ideas with users (like you) and developers.
>
> I hope this helps.
>
> Paul
>
>
>
>
> On Mon, 1 Feb 2010, nek5000-users at lists.mcs.anl.gov wrote:
>
> > Hi,
> >
> > I have experienced some difficulties with setting up a 2D conjugate heat
> > transfer problem. When generating the .rea file with prenek, I create a
> mesh
> > of, say 60 elements and define, say, 20 of those as solids (material group
> 1
> > instead of 0). Trying to run this .rea file in genmap (even when it is
> > converted to re2) does not work since the fluid boundary condition section
> of
> > the .rea file contains entries for elements (the solids) without
> > connectivity/boundary information. Deleting those entries (4*20 lines for
> 2D)
> > lets genmap accept the rea file and a map is created with one section for
> the
> > fluids elements and one for the solids.
> > Running these files with nek5000, however, resulted in element mismatch
> > errors, no matter what combination of solid elements vs. number of
> processors
> > I tried. The newest nek5000 svn was used.
> >
> > The solution I found is the following: when one takes the .rea file for
> > conjugate heat transfer (as it came out of prenek) and changes all solid
> > elements back to fluids, a map can be generated from this only-fluids-.rea
> > file with genmap that is essentially left in ignorance of the conjugate
> heat
> > transfer problem/solid elements. Running this "fluids-only" map with the
> > original "conjugate heat"-.rea file (having deleted the solid elements
> > boundary/connectivity information in its fluids section) will work and
> > produces physically correct looking results.
> >
> > My questions now are the following:
> > -Is this an appropriate way to go or am I missing some parameters that need
> > to be set? The rea files come directly out of prenek, so I didn't manually
> > change any parameters.
> > -By generating a map file that does not resemble what's truly going on, do
> we
> > loose parallel scalability and if so, is that a considerable penalty?
> > -I had some trouble defining solid elements in 3D with prenek. Whenever I
> > click on an element to change its material group, it does not do it. Is
> there
> > a fix?
> > -Expanding a 2D conjugate heat transfer problem with n2to3 does not work.
> All
> > upper level elements become solids. Is there a fix for this?
> >
> > I apologize for my lengthy postings,
> > Markus
> > _______________________________________________
> > 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
>