OPTIMAL TEMPERATURE DISTRIBUTION FOR A NONISOTHERMAL CAHN-HILLIARD SYSTEM IN TWO DIMENSIONS WITH SOURCE TERM AND DOUBLE OBSTACLE POTENTIAL

Abstract

In this note, we study the optimal control of a nonisothermal phase field system of Cahn-Hilliard type that constitutes an extension of the classical Caginalp model for nonisothermal phase transitions with a conserved order parameter. It couples a Cahn-Hilliard type equation with source term for the order parameter with the universal balance law of internal energy. In place of the standard Fourier form, the con­stitutive law of the heat flux is assumed in the form given by the theory developed by Green and Naghdi, which accounts for a possible thermal memory of the evolution. This has the consequence that the balance law of internal energy becomes a second-order in time equation for the thermal displacement or freezing index, that is, a primitive with respect to time of the temperature. Another particular feature of our system is the presence of the source term in the equation for the order parameter, which entails further mathematical difficulties because the mass conservation of the order parameter is no longer satisfied. In this paper, we study the case that the double-well potential driving the evolution of the phase transition is given by the nondifferentiable dou­ble obstacle potential, thereby complementing recent results obtained for the differentiable cases of regular and logarithmic potentials. Be­sides existence results, we derive first-order necessary optimality condi­tions for the control problem. The analysis is carried out by employing the so-called deep quench approximation in which the nondifferentiable double obstacle potential is approximated by a family of potentials of logarithmic structure for which meaningful first-order necessary opti­mality conditions in terms of suitable adjoint systems and variational inequalities are available. Since the results for the logarithmic poten­tials crucially depend on the validity of the so-called strict separation property which is only available in the spatially two-dimensional situ­ation, our whole analysis is restricted to the two-dimensional case.