Revision of the inextensibility constraint implementation without a penalty parameter

[blackwer]

The implementation of the inextensibility constraint is through a penalty method which has a tuning parameter that pretty much depends on the flow parameters. Simulation failures are largely due to this, so fiber solutions should be re-done to be more robust.