open TELEMAC-MASCARET The mathematically superior suite of solvers

Hello everyone,
is there a posibility to add a new porosity program in telemac2d sources. I know that telemac have already the corpor.f to model porosity, I'm just asking if I can use my own code.
Thank you in advance for any kind of help or advise.

Hello Mourad

The corpor.f subroutine is an empty one. It is intended to be programmed by the user. You can define the value of the porosity of each element.

There's an example that defines a quadrilateral and sets the porosity of the elements inside it as 0.5. You can use it as a guide to write your own code.

Best regards,

Phelype

Thanks Phelype for your quick reply,
Actually I know how to use corpor.f, I did it ùany times by the pasr, but I'm not realy stisfied by the result. The corpor code have a great advantage to other shemes I've found in the litterature cause the mass balance is fullfilled, however I can't understand why only a portion of the total volume is loaded in the porous medias, this volume corresponds to (n*V), n is the porosity.

My question is about wether is possible to couple a new code for Saint-Venant equations with porosity with the full finite element scheme of Telemeac2D.

Oh, I completely misunderstood your question.

In this case I really don't know how to answer your question. Sorry.

I'll leave it to the developers.

Best Regards,

Phelype

Hi,

It depends on how the other code is used.
If it is a library with a api you can add the call in telemac.
And link it with the library.
If it is just an executable it is more complicated.