Utilizing Block Restriction To Generate Two GlobalStrain Objects #27885
Replies: 2 comments 1 reply
-
Block restriction seems to be implemented so this looks like the expected syntax |
Beta Was this translation helpful? Give feedback.
-
@lkuna24 your zz epitaxy should be minus :) I just want to (re)iterate that this is part of a long standing issue within Ferret of computing ferroelectric thin film phase diagrams: See We get many requests (approx 1-2/month) to be able to solve this problem and I think its time we finally pin this down. We are slowly refactoring things to make them use more MOOSE systems but we still can't find what is going on. As soon as we throw in a second eigenstrain (proportional to the square of the phase field variable), we find that it dominates the calculation whereas instead we expect that at the end of the solve, the global strain of the film block should be very much equal to the eigenstrain "epitaxy" that we apply above. As such, the epitaxial condition (eigenstrain) coupling to the phase field problem is only a weak influence. |
Beta Was this translation helpful? Give feedback.
-
Hi All,
I am trying to carry out a simulation of a thin film with an applied epitaxy in the form of an eigenstrain. In the simulation the thin film sits on top of a substrate which imposes a uniform mismatch strain throughout the film (so long as the film is sufficiently thin). The domain (both film and substrate) is periodic in the x and y direction. Applying GlobalStrain to both the film and substrate with just the epitaxy (in the form of an eigenstrain) acting on the film yields some unexpected results, but block restricting and applying GlobalStrain to only the film with the epitaxy applied gives the expected result. I wanted to attempt to separate out GlobalStrain utilizing block restriction as shown in the input file below, but I get the error:
As I need the whole domain to be periodic, I believe I need a GlobalStrain on both blocks. According to the documentation here GlobalStrain, it seems as though I should be able to separate GlobalStrain Objects utilizing block restriction, or is this actually not possible?
Ultimately there will be another eigenstrain contributing to the thin film which will be coupled to the evolution of polarization throughout the volume of the film. But for now, I am just trying to see if I can establish the mechanical properties of the film and substrate.
Thanks for you time.
Beta Was this translation helpful? Give feedback.
All reactions