Closed RadostW closed 2 weeks ago
This looks great. You can add it as ex50.py, implement a test for it under tests/test_examples.py, and use the following as a template for listofexamples.rst (under fluid mechanics category, perhaps?):
Example 36: Nearly incompressible hyperelasticity
-------------------------------------------------
This example demonstrates the implementation of a two field mixed formulation
for nearly incompressible Neo-Hookean solids.
.. figure:: https://user-images.githubusercontent.com/22624037/91212007-4055aa80-e6d5-11ea-8572-f27986887331.png
The displacement contour of Example 36.
The figure was created using `ParaView <https://www.paraview.org/>`__.
See the source code of :exlink:`36` for more information.
Simply copy the image url from this issue and replace it in the template.
I suggest that you save the mesh to file through mesh.save(”docs/examples/meshes/cylinder.msh”), include the mesh in the commit and avoid using pygmsh by loading the mesh in the beginning of the example, via MeshTri.load(…). I think pygmsh is not included in our CI.
Finally, we accept nowadays only examples that have a matching license (BSD-3-clause) so hopefully you are able to use such a license instead of GPL.
We've managed to solve an advection-diffusion problem
$$0 = \Delta \phi - \mathrm{Pe} (u \cdot \nabla \phi) $$
in cylindrical coordinates with $U$ coming from Stokes flow around a sphere using
scikit-fem
package -- code hereI'd like to contribute it as an use example -- is there some clever way this code should be structured to work with
listofexamples.rst
?