examples.meshing.inputGrid2D

To run this example from the base FiPy directory, type:

$ python examples/meshing/inputGrid2D.py --numberOfElements=X

This example demonstrates how to build a 1D mesh and obtain basic mesh information. The command line argument, X, controls the number of elements on the mesh. Firstly parse the command line argument for numberOfElements, with the default set at 100.

>>> from fipy.tools.parser import parse
>>> numberOfElements = parse('--numberOfElements', action = 'store', type = 'int', default = 100)

A Grid2D object is invoked in the following way,

>>> from fipy import Grid2D, CellVariable, Viewer
>>> from fipy.tools import numerix

>>> nx = int(numerix.sqrt(numberOfElements))
>>> ny = nx
>>> dx = 1.
>>> dy = 1.
>>> mesh = Grid2D(nx = nx, ny = nx, dx = dx, dy = dy)

Once the mesh has been built information about the mesh can be obtained. For example the mesh volumes can be obtained with the getCellVolumes() method.

>>> vols = mesh.cellVolumes
>>> numerix.allclose(dx * dy * numerix.ones(nx * ny), vols)
1

Obtain the number of cells in the mesh

>>> N = mesh.numberOfCells
>>> numerix.allclose(N, numberOfElements)
1

Obtain all the left exterior faces, this is equal to ny.

>>> faces = mesh.facesLeft
>>> len(faces) == ny
1

One can view the mesh with the following code,

>>> if __name__ == '__main__':
...     viewer = Viewer(CellVariable(value = 0, mesh = mesh))
...     viewer.plot()

Last updated on Jun 26, 2024. Created using Sphinx 7.1.2.