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_pu... or https://www.sciencedirect.com/science/article/pii/S0021999111002361