On 09/05/2018 06:50 PM, Daniel Wheeler wrote:

Thanks for the references. I'll let you know how it goes. I might solve the elasticity part with Sfepy and the CH part with FiPy for now.

What is the best example in Sfepy to understand how to solve multiple PDEs?

Try looking at the examples in examples/multi_physics/ - it is usually elasticity + something.

r.

On Tue, Sep 4, 2018 at 5:33 PM Robert Cimrman <cimrman3@ntc.zcu.cz> wrote:

Hi Daniel,

On 09/04/2018 05:06 PM, Daniel Wheeler wrote:

Hi,

I'd like to try and solve the problem outlined here,

https://pages.nist.gov/pfhub/benchmarks/benchmark4.ipynb/

Is this a problem that Sfepy could handle? I'm assuming the elasticity part of the problem would be straightforward, but what about a fourth order phase field equation? Given that it's tractable are there examples that match the form of the elasticity equation and/or any Cahn-Hilliard type examples for Sfepy?

I have no experience with that kind of equation. However, I have found an article mentioning that the fourth order equation can be solved by splitting into two second order equations [1]. That should be doable with sfepy, talking about assembling the problem and getting the discretized system (if H1 elements are OK). Again, I have no experience with solving such systems.

Cheers, r.

[1] http://lib.dr.iastate.edu/cgi/viewcontent.cgi?article=1237&context=me_pubs or https://www.sciencedirect.com/science/article/pii/S0021999111002361

SfePy mailing list -- sfepy@python.org To unsubscribe send an email to sfepy-leave@python.org https://mail.python.org/mm3/mailman3/lists/sfepy.python.org/