[Nek5000-users] Setting up a case from scratch
nek5000-users at lists.mcs.anl.gov
nek5000-users at lists.mcs.anl.gov
Wed Nov 28 08:47:17 CST 2012
Note that this should be relatively trivial to do with the MOAB-based version of Nek (MOAB can read gmsh meshes
already). I'm not sure whether you can specify higher-order hexes in gmsh, but if you can, start by generating those
(27-node hex). Building and running Nek on MOAB is described in the Nek user's guide and on the nek wiki. The
examples/moab problem is a MOAB-based one, showing how to do this for a fluid-only problem; moab_conjht demonstrates it
for a fluid/solid problem. Note that the user documentation on text input is a bit out of date, so refer to the moab
example rather than the user's guide for that (will take this as an indication I need to update the user's guide and wiki).
As for boundary conditions, I'm not sure how this is represented in the gmsh files (whether groups of nodes/faces are
indicated in the file, or whether a field on nodes/faces/sides is used to indicate this). If the former, these should
read into MOAB sets marked with NEUMANN_SET or DIRICHLET_SET tags; if the latter, you'll have to interpret this data in
the .usr file to set the boundary conditions.
MOAB's gmsh reader is not parallelized, so if you want to run this problem in parallel, you'll need to translate into
MOAB's .h5m file format, then partition, as described in the wiki instructions. Building MOAB builds also 'mbconvert',
which converts between any file formats read/written by MOAB (for that matter, partitioning with MOAB's tool would allow
you to read the gmsh file directly and output the partitioning info with the mesh to .h5m).
- tim
On 11/27/2012 03:04 AM, nek5000-users at lists.mcs.anl.gov wrote:
> Hello
> I want to setup a case of 2d cylinder in a channel. I plan to do this way
>
> 1. Make mesh in gmsh
> 2. Convert to vtx
> 3. Setup in prex and read vtx mesh.
> 3.1 make curved sides
> 3.2 set boundary conditions
>
> I have already tried all these steps partly. I want to know if there is any easier way to do it. Especially making
> curved sides and setting boundary conditions needs to be done per boundary edge. Is there a way to make it faster ?
>
> Thanks
> praveen
>
>
> _______________________________________________
> Nek5000-users mailing list
> Nek5000-users at lists.mcs.anl.gov
> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
>
--
================================================================
"You will keep in perfect peace him whose mind is
steadfast, because he trusts in you." Isaiah 26:3
Tim Tautges Argonne National Laboratory
(tautges at mcs.anl.gov) (telecommuting from UW-Madison)
phone (gvoice): (608) 354-1459 1500 Engineering Dr.
fax: (608) 263-4499 Madison, WI 53706
More information about the Nek5000-users
mailing list