Phase field methods have been widely used to study phase transitions and polarization switching in ferroelectric thin films. In this paper, we develop an efficient numerical scheme for the variational phase field model based on variational forms of the electrostatic energy and the relaxation dynamics of the polarization vector. The spatial discretization combines the Fourier spectral method with the finite difference method to handle three-dimensional mixed boundary conditions. It allows for an efficient semi-implicit discretization for the time integration of the relaxation dynamics. This method avoids explicitly solving the electrostatic equilibrium equation (a Poisson equation) and eliminates the use of associated Lagrange multipliers. We present several numerical examples including phase transitions and polarization switching processes to demonstrate the effectiveness of the proposed method.
The study of polycrystalline materials requires theoretical and computational techniques enabling multiscale investigations. The amplitude expansion of the phase field crystal model (APFC) allows for describing crystal lattice properties on diffusive timescales by focusing on continuous fields varying on length scales larger than the atomic spacing. Thus, it allows for the simulation of large systems still retaining details of the crystal lattice. Fostered by the applications of this approach, we present here an efficient numerical framework to solve its equations. In particular, we consider a real space approach exploiting the finite element method. An optimized preconditioner is developed in order to improve the convergence of the linear solver. Moreover, a mesh adaptivity criterion based on the local rotation of the polycrystal is used. This results in an unprecedented capability of simulating large, three-dimensional systems including the dynamical description of the microstructures in polycrystalline materials together with their dislocation networks.
After decades of searching for robust nanoscale ferroelectricity that could enable integration into the next generation memory and logic devices, hafnia-based thin films have appeared as the ultimate candidate because their ferroelectric (FE) polarization becomes more robust as the size is reduced. This exposes a new kind of ferroelectricity, whose mechanism still needs to be understood. Towards this end, thin films with increased crystal quality are needed. We report the epitaxial growth of Hf0.5Zr0.5O2 (HZO) thin films on (001)-oriented La0.7Sr0.3MnO3/SrTiO3 (STO) substrates. The films, which are under epitaxial compressive strain and are predominantly (111)-oriented, display large FE polarization values up to 34 {mu}C/cm2 and do not need wake-up cycling. Structural characterization reveals a rhombohedral phase, different from the commonly reported polar orthorhombic phase. This unexpected finding allows us to propose a compelling model for the formation of the FE phase. In addition, these results point towards nanoparticles of simple oxides as a vastly unexplored class of nanoscale ferroelectrics.
We develop a phenomenological thermodynamic theory of ferroelectric BaTiO3 (BT) thin films epitaxially grown on cubic substrates using the Landau-Devonshire eight-order potential. The constructed misfit-temperature phase diagram is asymmetrical. We found that, overall view of the phase diagram depends on the values of compliances used in calculations and provide two qualitatively different diagrams. A thermodynamic path for BT film grown onto particular substrate can be found using a plot of the room-temperature tetragonal distortion (c-a)/a as a function of misfit strain.
We present a finite element based variational interface-preserving and conservative phase-field formulation for the modeling of incompressible two-phase flows with surface tension dynamics. The preservation of the hyperbolic tangent interface profile of the convective Allen-Cahn phase-field formulation relies on a novel time-dependent mobility model. The mobility coefficient is adjusted adaptively as a function of gradients of the velocity and the order parameter in the diffuse interface region in such a way that the free energy minimization properly opposes the convective distortion. The ratio of the convective distortion to the free energy minimization is termed as the convective distortion parameter, which characterizes the deviation from the hyperbolic tangent shape due to the convection effect. The mass conservation is achieved by enforcing a Lagrange multiplier with both temporal and spatial dependence on the phase-field function. We integrate the interface-preserving and conservative phase-field formulation with the incompressible Navier-Stokes equations and the continuum surface tension force model for the simulation of incompressible two-phase flows. A positivity preserving scheme designed for the boundedness and stability of the solution is employed for the variational discretization using unstructured finite elements. We examine the convergence and accuracy of the Allen-Cahn phase-field solver through a generic one-dimensional bistable diffusion-reaction system in a stretching flow. We quantify and systematically assess the relative interface thickness error and the relative surface tension force error with respect to the convective distortion parameter. Two- and three-dimensional rising bubble cases are further simulated to examine the effectiveness of the proposed model on the volume-preserving mean curvature flow and the interface-preserving capability.
In this paper, we develop an efficient lattice Boltzmann (LB) model for simulating immiscible incompressible $N$-phase flows $(N geq 2)$ based on the Cahn-Hilliard phase field theory. In order to facilitate the design of LB model and reduce the calculation of the gradient term, the governing equations of the $N$-phase system are reformulated, and they satisfy the conservation of mass, momentum and the second law of thermodynamics. In the present model, $(N-1)$ LB equations are employed to capture the interface, and another LB equation is used to solve the Navier-Stokes (N-S) equations, where a new distribution function for the total force is delicately designed to reduce the calculation of the gradient term. The developed model is first validated by two classical benchmark problems, including the tests of static droplets and the spreading of a liquid lens, the simulation results show that the current LB model is accurate and efficient for simulating incompressible $N$-phase fluid flows. To further demonstrate the capability of the LB model, two numerical simulations, including dynamics of droplet collision for four fluid phases and dynamics of droplets and interfaces for five fluid phases, are performed to test the developed model. The results show that the present model can successfully handle complex interactions among $N$ ($N geq 2$) immiscible incompressible flows.
Ruotai Li
,Qiang Du
,Lei Zhang
.
(2020)
.
"Numerical Discretization of Variational Phase Field Model for Phase Transitions in Ferroelectric Thin Films"
.
Ruotai Li
هل ترغب بارسال اشعارات عن اخر التحديثات في شمرا-اكاديميا