[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