On 09/12/2018 08:28 PM, Daniel Wheeler wrote:
So, I've implemented the elasticity part of the problem using Sfepy, see,
for the problem definition. My implementation is contained in
Basically, I'm using both the dw_lin_elastic and dw_lin_prestress to implement the two terms in the stress / strain equation. Does what I've done make sense to you with respect to how those terms are being used?
If I understand it correctly, dw_lin_elastic computes the elastic strain/stress parts for eta-dependent strain/stress in dw_lin_prestress? It looks OK to me (after a cursory glance).
I'm not currently solving the phase field equation for $\eta$. I'm simply imposing an $\eta$ field. Initially, I'm going to solve the phase field equation using FiPy and then use Sfepy for the elasticity. Eventually, I'd like to learn how to do it all in Sfepy.
The two-equations approach from that article seems to me to be a viable route.
PS: I needed to adjust myself a bit to the functional style of your code and I am starting liking it. What was your motivation for using toolz and what practical experience have you got with using it? Just curious :) (I am trying to use a more functional style in another Python-based project.)
Any advice would be greatly appreciated.
On Thu, Sep 6, 2018 at 7:58 AM Robert Cimrman firstname.lastname@example.org... wrote: >
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.