Journal of Applied Mathematics and Mechanics (ZAMM) / The polarization process of ferroelectric materials in the framework of variational inequalities

We are concerned with the mathematical modeling of the polarization process in ferroelectric media. We assume that this dissipative process is governed by two constitutive functions, which are the free energy function and the dissipation function. The dissipation function, which is closely connected to the dissipated energy, is usually nondifferentiable. Thus, a minimization condition for the overall energy includes the subdifferential of the dissipation function. This condition can also be formulated by way of a variational inequality in the unknown fields strain, dielectric displacement, remanent polarization and remanent strain. We analyze the mathematical wellposedness of this problem. We provide an existence and uniqueness result for the timediscrete update equation. Under stronger assumptions, we can prove existence of a solution to the timedependent variational inequality. To solve the discretized variational inequality, we use mixed finite elements, where mechanical displacement and dielectric displacement are unknowns, as well as polarization (and, if included in the model, remanent strain). It is then possible to satisfy Gauss' law of zero free charges exactly. We propose to regularize the dissipation function and solve for all unknowns at once in a single Newton iteration. We present numerical examples gained in the open source software package Netgen/NGSolve. (VLID)5630319 Version of record