[Nek5000-users] 答复: Nek5000-users Digest, Vol 93, Issue 6

Wed Nov 9 20:11:56 CST 2016

Hello, Nicolas
I read the documentation again carefully and very happy that my problem has
been solved. Thank you very much for advice. I would like to share the curve
definition parts in my mesh converting commands for a ring below:

Curves={4:{'type':'C','depth':4,'circle_center':(0,0),
'radius':1},3:{'type':'C','depth':4,'circle_center':(0,0),'radius':0.5}}

For 4 and 3 are outer and inner circle.

----------------------------------------------------------------
Bolun Xu
University of Science and Technology of China
Hefei, Anhui, China

-----邮件原件-----
发件人: nek5000-users-bounces at lists.mcs.anl.gov
[mailto:nek5000-users-bounces at lists.mcs.anl.gov] 代表
nek5000-users-request at lists.mcs.anl.gov
发送时间: 2016年11月9日 23:27
收件人: nek5000-users at lists.mcs.anl.gov
主题: Nek5000-users Digest, Vol 93, Issue 6

Send Nek5000-users mailing list submissions to
	nek5000-users at lists.mcs.anl.gov

To subscribe or unsubscribe via the World Wide Web, visit
	https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
or, via email, send a message with subject or body 'help' to
	nek5000-users-request at lists.mcs.anl.gov

You can reach the person managing the list at
	nek5000-users-owner at lists.mcs.anl.gov

When replying, please edit your Subject line so it is more specific than
"Re: Contents of Nek5000-users digest..."

Today's Topics:

   1. Re: A problem in Prenek (nek5000-users at lists.mcs.anl.gov)
   2. 2D circular ring mesh converting by mshconvert.py
      (nek5000-users at lists.mcs.anl.gov)
   3. Re: 2D circular ring mesh converting by	mshconvert.py
      (nek5000-users at lists.mcs.anl.gov)
   4. Re: Read Nek output in Matlab (nek5000-users at lists.mcs.anl.gov)
   5. Name of the variable in userf and userq
      (nek5000-users at lists.mcs.anl.gov)

----------------------------------------------------------------------

Message: 1
Date: Tue, 8 Nov 2016 21:32:25 +0000
From: nek5000-users at lists.mcs.anl.gov
To: "nek5000-users at lists.mcs.anl.gov"
	<nek5000-users at lists.mcs.anl.gov>
Subject: Re: [Nek5000-users] A problem in Prenek
Message-ID:
	<mailman.5887.1478640760.3602.nek5000-users at lists.mcs.anl.gov>
Content-Type: text/plain; charset="utf-8"

Hi Saikat,

I think this is what is expected. If we look at postpro.f in nek?s source
code (subroutine gen_rea_bc), there are 4 different formats in which the
boundary conditions will be written depending on the element count.
In your case, since the element count is greater than 1000 but less than
100,000, we expect the format to be: first 4 character spaces for boundary
condition, next 5 characters for element number, 1 character for edge/face
number depending on 2d/3d and so on. As a result, your boundary condition
line shows as:

 E   30724   2731.00       2.00000       0.00000       0.00000       0.00000

Thanks,
Ketan
-------------- next part --------------
An HTML attachment was scrubbed...
URL:
<http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20161108/6d5c2
cd9/attachment-0001.html>

------------------------------

Message: 2
Date: Wed, 9 Nov 2016 10:16:20 +0800
From: nek5000-users at lists.mcs.anl.gov
To: <nek5000-users at lists.mcs.anl.gov>
Subject: [Nek5000-users] 2D circular ring mesh converting by
	mshconvert.py
Message-ID:
	<mailman.5892.1478657813.3602.nek5000-users at lists.mcs.anl.gov>
Content-Type: text/plain; charset="us-ascii"

Hi, neks

I tried to convert a 2D circular ring mesh generated by gambit to .rea file
using the code mshconvert.py, but failed. I converted the example mesh file
sirkel4.msh following the step in README, then I have a question that what
command I should input when my mesh has more than one curved boundary just
like a ring mesh. For sirkel4.msh the command "7:{..}" means that the zone 7
is curved, if I have two curved zones, for example, zone 3 and 4, should I
input "3:{....},4:{....}" in one command line? 

Thank you in advance.

Regards,

Xu

----------------------------------------------------------------

Bolun Xu

University of Science and Technology of Chnia

Hefei, Anhui, China

-------------- next part --------------
An HTML attachment was scrubbed...
URL:
<http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20161109/51b79
3d8/attachment-0001.html>

------------------------------

Message: 3
Date: Wed, 9 Nov 2016 13:45:32 +0100
From: nek5000-users at lists.mcs.anl.gov
To: nek5000-users at lists.mcs.anl.gov
Subject: Re: [Nek5000-users] 2D circular ring mesh converting by
	mshconvert.py
Message-ID:
	<mailman.5908.1478698778.3602.nek5000-users at lists.mcs.anl.gov>
Content-Type: text/plain; charset="utf-8"

Dear Bolun,

Yes, you should specify the different curved zones separately. There is a
bit of documentation at the end of the file "mshconvert.py". Here is what is
said about the curves input:

"
curves = Dictionary of curve information. Keys are curve zones and value is
either
    {'type': 'C', 'radius': radius, 'circle_center': (x, y), 'depth': depth}
for a circle or
    {'type': 'm'}
for midpoint with nek5000 or
    {'type': 'spline'}
for a curved side with semtex. Here a .geom file containing the spline
information will be created.

The circle may provide the radius or the center of the circle through
'circle_center'. The curvature may also be used in the internal elements
inside the surface through specifying the depth. depth=4 means that the
curvature is used throughout the first four elements inside that surface.
This is necessary to get good quality meshes in, e.g., a cylinder. The
radius for an internal face is always computed as the distance to the
circle_center. Not for FEniCS.
"

I hope that helps.
Best,
Nicolas Offermans

On Wed, Nov 9, 2016 at 3:16 AM, <nek5000-users at lists.mcs.anl.gov> wrote:

> Hi, neks
>
> I tried to convert a 2D circular ring mesh generated by gambit to .rea 
> file using the code mshconvert.py, but failed. I converted the example 
> mesh file sirkel4.msh following the step in README, then I have a 
> question that what command I should input when my mesh has more than 
> one curved boundary just like a ring mesh. For sirkel4.msh the command 
> ?7:{??}? means that the zone 7 is curved, if I have two curved zones, 
> for example, zone 3 and 4, should I input ?3:{??..},4:{??..}? in one
command line?
>
> Thank you in advance.
>
> Regards,
>
> Xu
>
>
>
> ----------------------------------------------------------------
>
> Bolun Xu
>
> University of Science and Technology of Chnia
>
> Hefei, Anhui, China
>
>
>
>
>
> _______________________________________________
> Nek5000-users mailing list
> 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/20161109/63f8f
5a1/attachment-0001.html>

------------------------------

Message: 4
Date: Wed, 9 Nov 2016 13:50:59 +0100
From: nek5000-users at lists.mcs.anl.gov
To: nek5000-users at lists.mcs.anl.gov
Subject: Re: [Nek5000-users] Read Nek output in Matlab
Message-ID:
	<mailman.5909.1478698778.3602.nek5000-users at lists.mcs.anl.gov>
Content-Type: text/plain; charset=windows-1252; format=flowed

Hi Alessandro,

alternatively you can use the functions that I and a colleague of mine
wrote.
They can read and write nek's files.

You can find them here
https://github.com/nfabbiane/nekmatlab

Cheers,
Jacopo

On 2016-11-02 15:20, nek5000-users at lists.mcs.anl.gov wrote:
> Alessandro,
>
> Set, in usrdat2 in your .usr file the following lines:
>
>        param(66)=0
>        param(67)=0
>
> this will give you an ascii file output.  blah.fld01  blah.fld02 etc
>
> You can edit these files or write a short script to take off the 
> header or whatever and then load the data into matlab.
>
> Usual data lay out (for 2D) is:
>
> x y u v p T
>
> for fld01
>
> and
>
> u v p T
>
> for all other fields.
>
> Paul
>
> ________________________________________
> From: nek5000-users-bounces at lists.mcs.anl.gov 
> [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: Wednesday, November 02, 2016 7:44 AM
> To: nek5000-users at lists.mcs.anl.gov
> Subject: [Nek5000-users] Read Nek output in Matlab
>
> Hi Neks,
>
> I would like to make some postprocessing with the Nek output. Can I 
> read that files in Matlab (or C)?
>
> Just to make it clear: If my .rea file is called whatever.rea and the 
> output is whatever0.f00****
>
> I would like to read whatever0.f00**** and make some postprocessing. 
> It will be also very helpful to learn how to extrapolate the mesh so 
> that I can plot combinations of the solutions for instance.
>
> Thank you so much,
>
> Alessandro
>
>
> --
> Dr Alessandro Alla
>
> Florida State University
> Department of Scienfitic Computing
> 400 Dirac Science Library
> Tallahassee FL 32306-4120
>
> Room:   462
> Phone:  +1 8506441010
>
> E-Mail: aalla at fsu.edu
> Skype: alessandro.alla
> www.alessandroalla.com
>
> _______________________________________________
> 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

------------------------------

Message: 5
Date: Wed, 9 Nov 2016 16:27:12 +0100
From: nek5000-users at lists.mcs.anl.gov
To: nek5000-users at lists.mcs.anl.gov
Subject: [Nek5000-users] Name of the variable in userf and userq
Message-ID:
	<mailman.5924.1478705236.3602.nek5000-users at lists.mcs.anl.gov>
Content-Type: text/plain; charset=utf-8; format=flowed

Hi Neks,

I am trying to set volume forcing in my equations through userf and userq.
However it is not for the baseflow but for perturbations.

Does anyone know the name of the perturbations variables which are
equivalent to ffx, ffy, ffz, and qvol? So I can write something like :

ffxp = ...
ffyp = ...
ffzp = ...
qvolp = ...

Thank you for your time,

Best wishes,

Arnold

------------------------------

_______________________________________________
Nek5000-users mailing list
Nek5000-users at lists.mcs.anl.gov
https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users

End of Nek5000-users Digest, Vol 93, Issue 6
********************************************