Continuum thermodynamics of ferroelectric domain evolution: Theory, finite element implementation, and application to domain wall pinning

A continuum thermodynamics framework is devised to model the evolution of ferroelectric domain structures. The theory falls into the class of phase-field or diffuse-interface modeling approaches. Here a set of micro-forces and governing balance laws are postulated and applied within the second law of thermodynamics to identify the appropriate material constitutive relationships. The approach is shown to yield the commonly accepted Ginzburg-Landau equation for the evolution of the polarization order parameter. Within the theory a form for the free energy is postulated that can be applied to fit the general elastic, piezoelectric and dielectric properties of a ferroelectric material near its spontaneously polarized state. Thereafter, a principle of virtual work is specified for the theory and is implemented to devise a finite element formulation. The theory and numerical methods are used to investigate the fields near straight 180° and 90° domain walls and to determine the electromechanical pinning strength of an array of line charges on 180° and 90° domain walls.