upxo.meshing.conformal_mesher2d module

class upxo.meshing.conformal_mesher2d.confMesh2d(gtess=None)

Bases: object

An orchastrator class capable of multiple 2D conformal mehsing pipelines.

Usage

from upxo.meshing.conformal_mesher2d import confMesh2d as cm2d

ValidElTypesOptions = ('triangle', 'quad')

gtess

classmethod from_geometric_pxtal(gsgen_method='shapely_pxtal_load', pxtal=None, xbound=None, ybound=None)

Class method to mesh geometrified grain struture

femesh_pygmsh(elementShape='tri', elementOrder=1, meshingAlgorithmID=4, elsize_global=[1.0], intermediateFilename='femesh', intermediateFileformat='vtk')

Femesh pygmsh.

set_fids(geomObject)

Assign each cell feature, a unique ID (fid) starting from 1.

Parameters:

geomObject (list of shapely geometries Geometries that will be used to form feature IDs.)

Example

gsConfMesh.set_fids(gsConfMesh.gtess.L0.pxtal.geoms)

filter_mesh()

Filter mesh.

get_mesh_geometry()

Return the mesh geometry.

assess_quality(qualityMeasures=['aspect_ratio', 'skew', 'min_angle', 'area'])

Assess quality.

see_mesh_quality(mqm_data, mqm_dataframe, data_to_vis='mesh > quality > field', qualityMeasures=['aspect_ratio', 'skew', 'min_angle', 'area'])

See mesh quality.

see_femesh(*args, **kwargs)

See femesh.

find_meshCellFeatTypes()

Find meshCellFeatTypes.

find_lineFeatLocation()

Find lineFeatLocation.

find_GBlines()

Find GBlines.

find_availableFeatures()

Find availableFeatures.

find_availableElTypes()

Find availableElTypes.

find_availableElTypeID()

Find availableElTypeID.

build_nodes()

Build and return nodes.

rebuild_elConnectivity()

Rebuild elconnectivity.

get_elCentroids_singleElType(nodes, elConn)

Return the elCentroids singleElType.

form_elsets_elType()

Form elsets eltype.

get_elCentroids()

Return the elCentroids.

find_elID_ranges()

Find elID ranges.

form_elsets()

Consilidates elment sets of all element types into one single dictionary.

extract_gblines_for_grain(grain_name)

Extract gblines for grain.

extract_gblines_grains()

Extract gblines grains.

extract_gbnodes_grains()

Extract gbnodes grains.

extract_gbnodeCoords_grains()

Extract gbnodecoords grains.

extract_gbCoords_ckdTrees(preserve_ZDim=False)

Extract gbcoords ckdtrees.

calc_element_qualities(ar=True, saa=False, throw=False)

Return the element qualities.

find_elIDs_by_quality(**kwargs)

Find elIDs by quality.

plot_elements_geometric_grain(grain_name, **kwargs)

Visualise elements geometric grain using Matplotlib or PyVista.

plot_elements_by_elIDs(element_ids, **kwargs)

Visualise elements by elIDs using Matplotlib or PyVista.

plot_elements_geometric_grains(**kwargs)

Visualise elements geometric grains using Matplotlib or PyVista.

extract_gb_elements_for_grain(grain_name)

Extract gb elements for grain.

collect_gb_elements_for_grains(grain_ids=None, grain_names=None, prefix='grain')

Orchestrator to extract GB elements for multiple grains using extract_gb_elements_for_grain. Accepts either grain_ids or grain_names. If grain_ids are provided, uses prefix + ‘.’ + str(grainid) to form grain names.

see_gbElements_grains(grain_name, **kwargs)

See gbelements grains.

build_global_element_numbering()

Build and return global element numbering.

find_el_neigh(element_ids, n_order=1, eltype=None, include_self=False)

Find el neigh.

find_nthOrderNeigh(n_order, el_subset, include_self=False)

Find nthOrderNeigh.

resolve_eltypes(grain_name, grainElements, grainCoordinates, eltypes=None)

Resolve eltypes.

build_element_ids_by_band(grain_name, bands, eltypes, grainElements, nearest_dist_by_type)

Build and return element ids by band.

pick_contrasting_colours_from_cmap(n_colours, cmap_name='nipy_spectral')

Pick contrasting colours from cmap.

resolve_band_colours(bands, band_colours=None, auto_cmap='nipy_spectral')

Resolve band colours.

resolve_plot_eltype(eltypes, plot_eltype=None)

Resolve plot eltype.

plot_band_elements(element_ids_by_band, bands, plot_eltype, gbcoords, gblines_by_grain=None, colours_to_use=None, title='Selected elements by band', band_facecolors=False)

Visualise band elements using Matplotlib or PyVista.

select_elements_in_bands(grain_name, bands, gbnodeCoords, grainCoordinates, grainElements, eltypes=None, plot_eltype='quad', gblines_by_grain=None, plot=True, band_colours=None, auto_cmap='nipy_spectral', band_facecolors=False)

Select elements in bands.

fids

pxtal_mesh

grid

mesherTool

elementShape

elementOrder

meshingAlgorithmID

filtered_cells

filtered_mesh

nodes

elConn

meshCellFeatTypes

lineFeatLocation

GBlines

availableFeatures

availableElTypes

availableElTypeID

elsets_eltype

elsets

elID_ranges

class upxo.meshing.conformal_mesher2d.confMesh2dGMSH

Bases: object

Conformal 2D mesher using the raw gmsh API (no pygmsh dependency).

Designed to work directly with a dict of Shapely polygons (flat_cells) and an optional grain-id mapping (gid_map). This class duplicates analytical / post-processing methods from confMesh2d so that pygmsh can be phased out without breaking existing workflows.

Usage

from upxo.meshing.conformal_mesher2d import confMesh2dGMSH as cm2dg

Typical workflow

m = confMesh2dGMSH() m.femesh_gmsh(flat_cells, gid_map,

mesh_size_gb=0.75, mesh_size_bulk=4.5, mesh_algo=8, mesh_order=1, recombine_to_quads=True)

m.form_elsets_gmsh() m.build_boundary_nsets() m.build_grain_nsets() m.build_gb_nset()

ValidElTypesOptions = ('triangle', 'quad')

gtess

fids

gid_map

point_registry

surface_tags

physical_surface_tags

mesherTool

elementShape

elementOrder

meshingAlgorithmID

recombine_to_quads

nodes

elConn

GBlines

availableFeatures

availableElTypes

availableElTypeID

elsets_eltype

elsets

elID_ranges

nsets

grid

femesh_gmsh(flat_cells, gid_map=None, mesh_size_gb=1.0, mesh_size_bulk=4.0, mesh_algo=6, mesh_order=1, recombine_to_quads=False, out_dir=None, basename='gs_mesh', formats=None)

Full gmsh meshing pipeline.

Parameters:

• flat_cells (dict) – {flat_id: shapely.geometry.Polygon} — grain geometries.

• gid_map (dict, optional) – {flat_id: original_grain_id}. If None, flat_id is used as-is.

• mesh_size_gb (float) – Target element size on grain boundaries.

• mesh_size_bulk (float) – Target element size in grain interiors.

• mesh_algo (int) – gmsh mesh algorithm ID (e.g. 6=Frontal, 8=Frontal-Delaunay quads).

• mesh_order (int) – Element order (1=linear, 2=quadratic).

• recombine_to_quads (bool) – Whether to recombine triangles into quads after meshing.

• out_dir (str or None) – Directory to write exported mesh files. No export when None.

• basename (str) – Filename stem for exported files (extension appended per format).

• formats (list of str or None) – File format extensions to export, e.g. ['msh', 'inp', 'vtk'].

get_elCentroids_singleElType(nodes, elConn)

Return the elCentroids singleElType.

form_elsets_gmsh(flat_cells, prefix='grain.')

Build per-grain element sets from flat_cells polygons.

Parameters:

• flat_cells (dict) – {flat_id: shapely.geometry.Polygon} — same dict used in femesh_gmsh.

• prefix (str) – Prefix for element set names (default ‘grain.’).

find_elID_ranges()

Find elID ranges.

get_elCentroids()

Return the elCentroids.

build_boundary_nsets(tol=1e-06)

Build node sets for domain boundaries and corners.

Populates self.nsets with keys:

‘LEFT’, ‘RIGHT’, ‘BOTTOM’, ‘TOP’, ‘BOTTOM_LEFT’, ‘BOTTOM_RIGHT’, ‘TOP_LEFT’, ‘TOP_RIGHT’

using node coordinates from self.nodes.

build_grain_nsets()

Build per-grain node sets from self.elsets_eltype. Populates self.nsets with keys matching elset names (e.g. ‘grain.1’).

build_gb_nset()

Build a node set containing all grain-boundary nodes. Populates self.nsets[‘GB’].

extract_gblines_for_grain(grain_name)

Extract gblines for grain.

extract_gblines_grains()

Extract gblines grains.

extract_gbnodes_grains()

Extract gbnodes grains.

extract_gbnodeCoords_grains()

Extract gbnodecoords grains.

extract_gbCoords_ckdTrees(preserve_ZDim=False)

Extract gbcoords ckdtrees.

extract_gb_elements_for_grain(grain_name)

Extract gb elements for grain.

collect_gb_elements_for_grains(grain_ids=None, grain_names=None, prefix='grain')

Collect gb elements for grains.

get_mesh_geometry()

Return (points_2d, lines, triangles, quads) arrays for plotting.

see_femesh(*args, **kwargs)

See femesh.

see_gbElements_grains(grain_name, **kwargs)

See gbelements grains.

calc_element_qualities(ar=True, saa=False, throw=False)

Return the element qualities.

find_elIDs_by_quality(**kwargs)

Find elIDs by quality.

find_el_neigh(element_ids, n_order=1, eltype=None, include_self=False)

Find el neigh.

find_nthOrderNeigh(n_order, el_subset, include_self=False)

Find nthOrderNeigh.

build_global_element_numbering()

Build and return global element numbering.

resolve_eltypes(grain_name, grainElements, grainCoordinates, eltypes=None)

Resolve eltypes.

build_element_ids_by_band(grain_name, bands, eltypes, grainElements, nearest_dist_by_type)

Build and return element ids by band.

pick_contrasting_colours_from_cmap(n_colours, cmap_name='nipy_spectral')

Pick contrasting colours from cmap.

resolve_band_colours(bands, band_colours=None, auto_cmap='nipy_spectral')

Resolve band colours.

resolve_plot_eltype(eltypes, plot_eltype=None)

Resolve plot eltype.

plot_elements_geometric_grain(grain_name, **kwargs)

Visualise elements geometric grain using Matplotlib or PyVista.

plot_elements_by_elIDs(element_ids, **kwargs)

Visualise elements by elIDs using Matplotlib or PyVista.

plot_elements_geometric_grains(**kwargs)

Visualise elements geometric grains using Matplotlib or PyVista.

plot_band_elements(element_ids_by_band, bands, plot_eltype, gbcoords, gblines_by_grain=None, colours_to_use=None, title='Selected elements by band', band_facecolors=False)

Visualise band elements using Matplotlib or PyVista.

select_elements_in_bands(grain_name, bands, gbnodeCoords, grainCoordinates, grainElements, eltypes=None, plot_eltype='quad', gblines_by_grain=None, plot=True, band_colours=None, auto_cmap='nipy_spectral', band_facecolors=False)

Select elements in bands.