Mapping OpenFOAM results into dolfinx

[dokken]

If you can read it into a pyvista unstructured grid, you can extract the topology and geometry and create a dolfinx mesh. I have tried to explain this in the following tutorial:
http://jsdokken.com/FEniCS23-tutorial/src/mesh_generation.html

You should also then consider mesh.topology.original_cell_index and mesh.geometry.input_global_indices, as for instance explained in

Wrong order of dofs after creating mesh from lists of points and cells

You are making some major assumptions here that simply cannot be held by DOLFINx This is because: Only simple function spaces has the same number of dofs as the mesh geometry The input ordering does not make sense once one wants to run in parallel. More information about this has been discussed in: Consider the following MWE to map a vertex->mesh node->original_input_position import dolfinx import matplotlib.pyplot as plt import numpy as np import ufl from mpi4py import MPI pp = np.arra…

How to apply a boundary condition to a range of DOFs?

I don’t have time to read that user manual. Assuming that what is in NSET describes the original nodes in the mesh you would like to constrain, you need to: Create a map from the input node index to a vertex on your process, this is simply a matter of inverting mesh.geometry.input_global_indices. Once you have a map from the global indices (in your case a map from 97, 96, 95, ... to what is local on the process), you remove all the results that are -1. You then map from geometry to topology …