Abstract. Optimal control of magnetization in a ferromagnet is formulated as a mathematical program with evolutionary equilibrium constraints. To this purpose, we construct an evolutionary infinite-dimensional model which is discretized both in the space as well as in time variables.
The evolutionary nature of this equilibrium is due to the hysteresis behavior of the respective magnetization process. To solve the problem numerically, we adapted the implicit programming technique. The adjoint equations, needed to compute subgradients of the composite objective, are derived using the generalized differential calculus of B. Mordukhovich. We solve two test examples and discuss numerical results.