===affil2: Argonne National Laboratory ===firstname: Jungho ===firstname4: Lois ===firstname3: Mihai ===lastname2: Abhyankar ===lastname: Lee ===firstname5: Todd ===affil6: Argonne National Laboratory ===lastname3: Anitescu ===email: julee@mcs.anl.gov ===lastname6: Smith ===affil5: Argonne National Laboratory ===otherauths: ===lastname4: McInnes ===affil4: Argonne National Laboratory ===lastname7: Wang ===affil7: Argonne National Laboratory ===firstname7: Lei ===postal: Argonne National Laboratory 9700 South Cass Avenue Building 240 Argonne, IL 60439 ===firstname6: Barry ===ABSTRACT: Recent progress on the development of scalable differential variational inequality multigrid-based solvers for the phase-field approach to mesoscale materials modeling is described. We have developed an active-set method for variational inequalities in PETSc, leveraging experience by the optimization community in TAO. A geometric multigrid solver in PETSc is used to solve the resulting linear systems. We present strong and weak scaling results for 2D coupled Allen-Cahn/Cahn-Hilliard systems. ===affil3: Argonne National Laboratory ===lastname5: Munson ===affilother: ===title: Scalable DVI and Multigrid for phase-field equations ===firstname2: Shrirang