Francesco Caimmi <francesco.caimmi(a)polimi.it> writes:
first of all my thanks to Andreas for making Hedge
I am looking for a solution to do some DG analysis in the field of solid
mechanics, and Hedge looked promising, as I am pretty accustomed to Python.
However by looking at the documentation I am not really sure if it is really
the tool for me, nor if the effort I would need to make it work as I wish
would be sensible or not.
I am interested in solid mechanics problems which have solutions with evolving
discontinuities across some subdomain, usually a surface or a line ( I am
thinking of crack propagation inside bodies).
I know there many techniques to address allow for such a feature in the
solution, so I am wondering if Hedge have some built-in support for this kind
Thanks for your message, and sorry about the delay in responding. Hedge
has some amount of built-in support for artificial viscosity, which lets
it--to some extent--deal with discontinuities (shocks, contact d.)
arising in high-speed flow problems.
The precise procedure is described here:
I don't believe that this approach is very suited to crack propagation
however--other things like XFEM may ultimately serve you better.
My medium-term plan for hedge is to continue to support its operator
language, but to replace its computational guts (JIT compilation and
PyCUDA) with a more modern approach based on Loopy  and meshmode
. I'm always looking for allies in making this happen, and I think
what you describe could fit well with what I'm imagining.
Let me know if you're interested, and I'd be happy to talk more.