Isosurface¶

Description¶

Perform an isosurface using VTK.

==>

On the left: initial unstructured multi-element mesh and \(\rho u_x\) field.

On the right: resulting unstructured mesh after Isosurface Treatment with value = 0.

Parameters¶

base: Base
The base on which the isosurfacing will be applied.

coordinates: list(str), default= <base>.coordinate_names
The variable names that define the set of coordinates. The default coordinates comes from the input base attribute coordinate_names.
variable: str,
Name of the variable to make the isosurface.
value: float or list(float),
Value of the isosurface or list of values for isosurfaces.
with_families: bool, default= False
If True, the output of the treatment contains rebuilt families based on the input base. All the families and all the levels of sub-families are rebuilt, but only attributes and Zone are transfered (not yet implemented for Boundary and Window).
output_element: str
Options dedicated to expert users. Values are [‘triangle’]. ‘triangle’: do triangulate the isosurface. All elements are triangulated. Only used when polyhedral input.

Preconditions¶

Coordinates and variables must be available at nodes.

Zones may be either structured or unstructured.

Zones may contain multiple instants.

Postconditions¶

The output base is always unstructured. If with_families is enabled, the output base contains the reconstructed families of base.

Example¶

This example shows how to extract an isosurface. To be able to use this functionnality you must have vtk installed.

import antares

myt = Treatment('isosurface')
myt['base'] = base
myt['variable'] = 'Ux'
myt['value'] = 0.
base_iso = myt.execute()

Main functions¶

class antares.treatment.TreatmentIsosurface.TreatmentIsosurface¶

Class for isosurfaces.

execute()¶

Create isosurfaces.

Returns:: an unstructured base or a list of unstructured bases sorted as the list of provided iso values.
Return type:: Base or list(Base)

Example¶

"""
This example shows how to extract an isosurface.
To be able to use this functionnality you must have vtk installed.
"""
import os
if not os.path.isdir('OUTPUT'):
    os.makedirs('OUTPUT')

from antares import Reader, Treatment, Writer

# ------------------
# Reading the files
# ------------------
reader = Reader('bin_tp')
reader['filename'] = os.path.join('..', 'data', 'ROTOR37', 'GENERIC', 'flow_<zone>_ite0.dat')
ini_base = reader.read()

# -----------
# Isosurface
# -----------
treatment = Treatment('isosurface')
treatment['base'] = ini_base
treatment['variable'] = 'rovx'
treatment['value'] = 0.
result = treatment.execute()

# -------------------
# Writing the result
# -------------------
writer = Writer('bin_tp')
writer['filename'] = os.path.join('OUTPUT', 'ex_isosurface.plt')
writer['base'] = result
writer.dump()

# -----------
# Isosurfaces
# -----------
treatment = Treatment('isosurface')
treatment['base'] = ini_base
treatment['variable'] = 'rovx'
treatment['value'] = [0., 0.01]
result = treatment.execute()

for num, base in enumerate(result):
    writer = Writer('bin_tp')
    writer['filename'] = os.path.join('OUTPUT', 'ex_isosurface{}.plt'.format(num))
    writer['base'] = base
    writer.dump()