[Nek5000-users] Mesh with two blocks

nek5000-users at lists.mcs.anl.gov nek5000-users at lists.mcs.anl.gov
Fri Jun 1 14:34:14 CDT 2018


Hi Marco,


The way I normally do this is as follows:


  1.  build a valid fluid mesh ( "f.rea" )
  2.  build a valid solid mesh ("s.rea")
  3.  run pretex and pick option "3"
  4.  Enter session name (say, "c" for "combined")
  5.  Enter fluid session name, f
  6.  Enter solid session name, s


Hopefully, it produces a valid "c.rea" file ready for conjugate heat transfer.


Let me know if that doesn't work for you.


Best,


Paul


________________________________
From: Nek5000-users <nek5000-users-bounces at lists.mcs.anl.gov> on behalf of nek5000-users at lists.mcs.anl.gov <nek5000-users at lists.mcs.anl.gov>
Sent: Friday, June 1, 2018 2:10:04 PM
To: nek5000-users at lists.mcs.anl.gov
Subject: Re: [Nek5000-users] Mesh with two blocks

Paul,

thank you again for your reply.

When I generated a mesh in Cubit that has two blocks and then run it with Nek5000 I get nelv = nelt. How does Nek5000 know that block 1 is the fluid mesh and block 2 is the solid mesh?

I think I missing a step in the setup of my problem. I wonder how the mesh used in conj_ht was generated?

Thanks,

Marco

On Jun 1, 2018, at 3:04 PM, nek5000-users at lists.mcs.anl.gov<mailto:nek5000-users at lists.mcs.anl.gov> wrote:


Depends on which routine.  Suppose userbc:

      if (ifield.eq.1) then ! fluid field

      else   ! thermal field

      endif

Or... if in userchk:

      integer e

      do e=1,nelt

           if (e.le.nelv) then ! fluid + temperature

           else    !  solid==temperature only

           endif

     enddo

Local elements are sorted:  1...nelv is fluid

1...nelt is thermal

nelv is <= nelt

Paul

________________________________
From: Nek5000-users <nek5000-users-bounces at lists.mcs.anl.gov<mailto:nek5000-users-bounces at lists.mcs.anl.gov>> on behalf of nek5000-users at lists.mcs.anl.gov<mailto:nek5000-users at lists.mcs.anl.gov> <nek5000-users at lists.mcs.anl.gov<mailto:nek5000-users at lists.mcs.anl.gov>>
Sent: Friday, June 1, 2018 1:53:03 PM
To: nek5000-users at lists.mcs.anl.gov<mailto:nek5000-users at lists.mcs.anl.gov>
Subject: Re: [Nek5000-users] Mesh with two blocks

Thank you for your replies.

I am actually interested in generating a mesh for a conjugate gradient solve. In the conj_ht example, the solid block is differentiated from the liquid block by testing the element index, i.e., ieg > nelgt. Assuming that I generate a mesh with two blocks that represent the fluid and solid regions, respectively, how do I access these blocks in the usr file?

Marco

> On Jun 1, 2018, at 2:35 PM, nek5000-users at lists.mcs.anl.gov<mailto:nek5000-users at lists.mcs.anl.gov> wrote:
>
> Why don't you use the sidesets to set your BC?
>
> -----Original message-----
>> From:nek5000-users at lists.mcs.anl.gov<mailto:nek5000-users at lists.mcs.anl.gov> <nek5000-users at lists.mcs.anl.gov<mailto:nek5000-users at lists.mcs.anl.gov>>
>> Sent: Friday 1st June 2018 20:06
>> To: nek5000-users at lists.mcs.anl.gov<mailto:nek5000-users at lists.mcs.anl.gov>
>> Subject: [Nek5000-users] Mesh with two blocks
>>
>> Hi,
>>
>> I would like to get some insights on how Nek5000 deals with meshes that contains more than one block. I generated a 2-D mesh with Cubit and was able to convert it to a re2 format with exo2nek but I do not know how to differentiate each block in the usr file to set boundary conditions for instance.
>>
>> I went through the examples and did not find anything relevant to my question. Does anybody have experience with multi-block mesh? What is the proper approach when dealing with multi-block meshes in Nek5000?
>>
>> Thanks,
>>
>> Marco
>> _______________________________________________
>> Nek5000-users mailing list
>> Nek5000-users at lists.mcs.anl.gov<mailto: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<mailto: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<mailto: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<mailto: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/20180601/6fb8d56a/attachment.html>


More information about the Nek5000-users mailing list