With the amazing functionality offered by the LEoPart library I am able to track “strain-history” of particles inside a complex domain.
In the picture above, one can distinguish the particles inside the red box that will leave through Outlet1 and the particles inside the blue box leaving through Outlet2. These particles will carry information of their own “strain-history”.
For my application I am interested in the cumulative strain-history of particles leaving a specific boundary. For example, I would like to know the average of strain-histories of particles leaving through outlet1 and compare the average strain-histories of particles leaving through outlet2. With the effect of the stenosis, I expect higher strain-histories carried by particles inside the red box.
LEoPart functionality enables specification of “open” boundaries, which means particles will be deleted if they “escaped” through a facet on an open boundary. To this end, one should specify an “open boundary” by marking it with integer “2” in a MeshFunction object.
Currently, I have a model with 1 inlet, 2 outlets and a wall. From a .h5-file I obtain a MeshFunction object carrying the markers for these boundaries: the inlet is marked as “1”, outlet 1 as “2”, outlet 3 as “3” and wall is marked as “4”. I use these markers to assign boundary conditions and calculate flow at the specific outlets, for example.
To enable LEoPart “open” boundary functionality, I would require a second MeshFunction object with different markers.
Is it possible to create a second meshfunction object and to reassign markers conveniently?
For my application I would need to reassign facets carrying markers “2” and “3”, so that they both carry marker “2” to become “open” in LEoPart.
With kind regards,