Dynamics of the viscous Cahn-Hilliard equation

We study generalized viscous Cahn-Hilliard problems with nonlinearities satisfying critical growth conditions in W-0(1,p)(Omega), where Omega is a bounded smooth domain in R-n, n >= 3. In the critical growth case, we prove that the problems are locally well posed and obtain a bootstrapping procedure showing that the solutions are classical. For p = 2 and almost critical dissipative nonlinearities we prove global well posedness, existence of global attractors in H-0(1)(Omega) and, uniformly with respect to the viscosity parameter, L-infinity(Omega) bounds for the attractors. Finally, we obtain a result on continuity of regular attractors which shows that, if n = 3, 4, the attractor of the Cahn-Hilliard problem coincides (in a sense to be specified) with the attractor for the corresponding semilinear heat equation. (C) 2008 Elsevier Inc. All rights reserved.

The persistence of phase-separation in LiFePO4 with two-dimensional Li+ transport : the Cahn-Hilliard-reaction equation and the role of defects

We examine the solution of the two-dimensional Cahn-Hilliard-reaction (CHR) equation in the xy plane as a model of Li+ intercalation into LiFePO4 material. We validate our numerical solution against the solution of the depth-averaged equation, which has been used to model intercalation in the limit of highly orthotropic diffusivity and gradient penalty tensors. We then examine the phase-change behaviour in the full CHR system as these parameters become more isotropic, and find that as the Li+ diffusivity is increased in the x direction, phase separation persists at high currents, even in small crystals with averaged coherency strain included. The resulting voltage curves decrease monotonically, which has previously been considered a hallmark of crystals that fill homogeneously.

A least squares based finite volume method for the Cahn-Hilliard and Cahn-Hilliard-reaction equations

A vertex-centred finite volume method (FVM) for the Cahn-Hilliard (CH) and recently proposed Cahn-Hilliard-reaction (CHR) equations is presented. Information at control volume faces is computed using a high-order least-squares approach based on Taylor series approximations. This least-squares problem explicitly includes the variational boundary condition (VBC) that ensures that the discrete equations satisfy all of the boundary conditions. We use this approach to solve the CH and CHR equations in one and two dimensions and show that our scheme satisfies the VBC to at least second order. For the CH equation we show evidence of conservative, gradient stable solutions, however for the CHR equation, strict gradient-stability is more challenging to achieve.

An extended Cahn-Hilliard model for interfaces with cubic anisotropy

For studying systems with a cubic anisotropy in interfacial energy sigma, we extend the Cahn-Hilliard model by including in it a fourth-rank term, namely, gamma (ijlm) [partial derivative (2) c/(partial derivativex(i) partial derivativex(j))] [partial derivative (2) c/(partial derivativex(l) partial derivativex(m))]. This term leads to an additional linear term in the evolution equation for the composition parameter field. It also leads to an orientation-dependent effective fourth-rank coefficient gamma ([hkl]) in the governing equation for the one-dimensional composition profile across a planar interface. The main effect of a non-negative gamma ([hkl]) is to increase both sigma and interfacial width w, each of which, upon suitable scaling, is related to gamma ([hkl]) through a universal scaling function. In this model, sigma is a differentiable function of interface orientation (n) over cap, and does not exhibit cusps; therefore, the equilibrium particle shapes (Wulff shapes) do not contain planar facets. However, the anisotropy in the interfacial energy can be large enough to give rise to corners in the Wulff shapes in two dimensions. In particles of finite sizes, the corners become rounded, and their shapes tend towards the Wulff shape with increasing particle size.

Simulação tridimensional adaptativa da separação das fases de uma mistura bifásica usando a equação de Cahn-Hilliard

Simulamos a separação dos componentes de uma mistura bifásica com a equação de Cahn-Hilliard. Esta equação contém intrincados termos não lineares e derivadas de alta ordem. Além disso, a delgada região de transição entre os componentes da mistura requer muita resolução. Assim, determinar a solução numérica da equação de Cahn-Hilliard não é uma tarefa fácil, principalmente em três dimensões. Conseguimos a resolução exigida no tempo usando uma discretização semi-implícita de segunda ordem. No espaço, obtemos a precisão requerida utilizando malhas refinadas localmente com a estratégia AMR. Essas malhas se adaptam dinamicamente para recobrir a região de transição. O sistema linear proveniente da discretização é solucionado por intermédio de técnicas multinível-multigrid.

A QUADRATIC C degrees INTERIOR PENALTY METHOD FOR LINEAR FOURTH ORDER BOUNDARY VALUE PROBLEMS WITH BOUNDARY CONDITIONS OF THE CAHN-HILLIARD TYPE

We develop a quadratic C degrees interior penalty method for linear fourth order boundary value problems with essential and natural boundary conditions of the Cahn-Hilliard type. Both a priori and a posteriori error estimates are derived. The performance of the method is illustrated by numerical experiments.

Fast iterative solvers for Cahn-Hilliard problems

Otto-von-Guericke-Universität Magdeburg, Fakultät für Mathematik, Dissertation, 2016

A model for diffusive and displacive phase transitions: Thermo-chemo-mechanical coupling effects

A diffuse interface phase field model is proposed for the unified analysis of diffusive and displacive phase transitions under nonisothermal conditions. Two order parameters are used for the description of the phenomena: one is related to the solute mass fraction and the other to the strain. The model governing equations come from the balance of linear momentum, the solute mass balance (which will lead to the Cahn-Hilliard equation) and the balance of internal energy. Thermodynamic restrictions allow to define constitutive relations for the thermodynamic forces and for the mechanical and chemical dissipations. Numerical tests carried out at different values of the initial temperature show that the model is able to describe the main features of both the displacive and the diffusive phase transitions, as well as their effect on the temperature. © 2010, Advanced Engineering Solutions.

Three-dimensional, fully adaptive simulations of phase-field fluid models

We present an efficient numerical methodology for the 31) computation of incompressible multi-phase flows described by conservative phase-field models We focus here on the case of density matched fluids with different viscosity (Model H) The numerical method employs adaptive mesh refinements (AMR) in concert with an efficient semi-implicit time discretization strategy and a linear, multi-level multigrid to relax high order stability constraints and to capture the flow`s disparate scales at optimal cost. Only five linear solvers are needed per time-step. Moreover, all the adaptive methodology is constructed from scratch to allow a systematic investigation of the key aspects of AMR in a conservative, phase-field setting. We validate the method and demonstrate its capabilities and efficacy with important examples of drop deformation, Kelvin-Helmholtz instability, and flow-induced drop coalescence (C) 2010 Elsevier Inc. All rights reserved

Incorporating memory effects in phase separation processes

We consider the modification of the Cahn-Hilliard equation when a time delay process through a memory function is taken into account. We then study the process of spinodal decomposition in fast phase transitions associated with a conserved order parameter. Finite-time memory effects are seen to affect the dynamics of phase transition at short times and have the effect of delaying, in a significant way, the process of rapid growth of the order parameter that follows a quench into the spinodal region. These effects are important in several systems characterized by fast processes, like non-equilibrium dynamics in the early universe and in relativistic heavy-ion collisions. (C) 2006 Elsevier B.V. All rights reserved.

Memory effect and fast spinodal decomposition

We consider the modification of the Cahn-Hilliard equation when a time delay process through a memory function is taken into account. We then study the process of spinodal decomposition in fast phase transitions associated with a conserved order parameter. The introduced memory effect plays an important role to obtain a finite group velocity. Then, we discuss the constraint for the parameters to satisfy causality. The memory effect is seen to affect the dynamics of phase transition at short times and have the effect of delaying, in a significant way, the process of rapid growth of the order parameter that follows a quench into the spinodal region.

A model for diffusive and displacive phase transitions in steels

Heat treatment of steels is a process of fundamental importance in tailoring the properties of a material to the desired application; developing a model able to describe such process would allow to predict the microstructure obtained from the treatment and the consequent mechanical properties of the material. A steel, during a heat treatment, can undergo two different kinds of phase transitions [p.t.]: diffusive (second order p.t.) and displacive (first order p.t.); in this thesis, an attempt to describe both in a thermodynamically consistent framework is made; a phase field, diffuse interface model accounting for the coupling between thermal, chemical and mechanical effects is developed, and a way to overcome the difficulties arising from the treatment of the non-local effects (gradient terms) is proposed. The governing equations are the balance of linear momentum equation, the Cahn-Hilliard equation and the balance of internal energy equation. The model is completed with a suitable description of the free energy, from which constitutive relations are drawn. The equations are then cast in a variational form and different numerical techniques are used to deal with the principal features of the model: time-dependency, non-linearity and presence of high order spatial derivatives. Simulations are performed using DOLFIN, a C++ library for the automated solution of partial differential equations by means of the finite element method; results are shown for different test-cases. The analysis is reduced to a two dimensional setting, which is simpler than a three dimensional one, but still meaningful.

Decomposition driven interface evolution for layers of binary mixtures: I. Model dirivation and base states

A dynamical model is proposed to describe the coupled decomposition and profile evolution of a free surfacefilm of a binary mixture. An example is a thin film of a polymer blend on a solid substrate undergoing simultaneous phase separation and dewetting. The model is based on model-H describing the coupled transport of the mass of one component (convective Cahn-Hilliard equation) and momentum (Navier-Stokes-Korteweg equations) supplemented by appropriate boundary conditions at the solid substrate and the free surface. General transport equations are derived using phenomenological nonequilibrium thermodynamics for a general nonisothermal setting taking into account Soret and Dufour effects and interfacial viscosity for the internal diffuse interface between the two components. Focusing on an isothermal setting the resulting model is compared to literature results and its base states corresponding to homogeneous or vertically stratified flat layers are analyzed.