[Lewis Liang]

Currently, I am trying to mesh a vertical cylinder by using the m4 macro provides by Ehsan Madadi. The script is shown below:

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     |                                                 |
|   \\  /    A nd           | Copyright (C) 2016 Ehsan Madadi-Kandjani        |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    `format'      ascii;
    class       dictionary;
    object      blockMeshDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
// General macros to create cylinder mesh

changecom(//)changequote([,])
define(calc, [esyscmd(perl -e 'use Math::Trig; print ($1)')])
define(VCOUNT, 0)
define(vlabel, [[// ]Vertex $1 = VCOUNT define($1, VCOUNT)define([VCOUNT], incr(VCOUNT))])

define(hex2D, hex ($1b $2b $3b $4b $1t $2t $3t $4t))
define(btQuad, ($1b $2b $2t $1t))
define(topQuad, ($1t $4t $3t $2t))
define(bottomQuad, ($1b $2b $3b $4b))

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

convertToMeters 0.001;

// Inner square side half
define(s, 0.35)

// Inner square side curvature
define(sc, 0.4)

// cylinder radius
define(r, 0.7)

// Height of cylinder
define(z, 20.0)

// Base z
define(Zb, 0)

// Outlet z
define(Zt, calc(Zb + z))

// Number of cells at inner square
define(Ns, 20)

// Number of cells between inner square and circle
define(Ni, 10)

// Number of cells in the cylinder height
define(Nz, 100)

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

As introduced from the http://www.ehsanmadadi.com/cylinder-mesh/ and the OpenFOAMWiki https://openfoamwiki.net/index.php/H...ckMesh_with_m4. I just can mesh a cylinder with a centre of bottom plane located in (0 0 0). It is straightforward to set up the length and original point in the Z direction. However, I am confusing to modify original point in XY plane, for example, to build a cylinder whose central point of bottom plane located in (10 10 0). Could anyone face this same problem?

Thanks in advance!

Zuodong

[Antimony]

Hi,

One possible solution would be to use transformPoints and translate by the vector (10 10 0) after you have created the mesh using blockMesh

Cheers,
Antimony

[Lewis Liang]

Hi Antimony! It works now after typing this command:

transformPoints -translate '(10 10 0)'

Many thanks!

Zuodong