scholarly journals Phase-field modelling and computing for a large number of phases

2019 ◽  
Vol 53 (3) ◽  
pp. 805-832
Author(s):  
Élie Bretin ◽  
Roland Denis ◽  
Jacques-Olivier Lachaud ◽  
Édouard Oudet

We propose a framework to represent a partition that evolves under mean curvature flows and volume constraints. Its principle follows a phase-field representation for each region of the partition, as well as classical Allen–Cahn equations for its evolution. We focus on the evolution and on the optimization of problems involving high resolution data with many regions in the partition. In this context, standard phase-field approaches require a lot of memory (one image per region) and computation timings increase at least as fast as the number of regions. We propose a more efficient storage strategy with a dedicated multi-image representation that retains only significant phase field values at each discretization point. We show that this strategy alone is unfortunately inefficient with classical phase field models. This is due to non local terms and low convergence rate. We therefore introduce and analyze an improved phase field model that localizes each phase field around its associated region, and which fully benefits of our storage strategy. To demonstrate the efficiency of the new multiphase field framework, we apply it to the famous 3D honeycomb problem and the conjecture of Weaire–Phelan’s tiling.

2011 ◽  
Vol 421 ◽  
pp. 90-97 ◽  
Author(s):  
Yun Chen ◽  
Na Min Xiao ◽  
Xiu Hong Kang ◽  
Dian Zhong Li

To understand the dendrite formation during solidification phase-field model has become a powerful numerical method of simulating crystal growth in recent years. Two phase-field models due to Wheeler et al. and Karma et al., respectively, have been employed for modeling the dendrite growth worldwidely. The comparison of the two models was performed. Then using the adaptive finite element method, both models were solved to simulate a free dendrite growing from highly undercooled melts of nickel at various undercoolings. The simulated results showed that the discrepancy between the two phase-field models is negligible. Careful comparison of the phase-filed simulations with LKT(BCT) theory and experimental data were carried out, which demonstrated that the phase-field models are able to quantitatively simulate the dendrite growth of nickel at low undercoolings, however, at undercoolings above ten percent of the melting point (around 180K), the simulated velocities by Wheeler and Karma model as well as the analytical predictions overestimated the reported experiment results.


2020 ◽  
Author(s):  
Mengfei Li ◽  
Karel Matouš ◽  
Robert Nerenberg

<p>In many environmental systems, such as membrane filtration systems, biofilm control is essential, but costly and requiring harsh chemicals. More effective biofilm control may be obtained using a “materials science” approach.  Biofilms can be characterized as viscoelastic materials, and biofilm “disruptors” can be characterized for their weakening effect on biofilm mechanical strength. By using a novel mathematical model that incorporates biofilm mechanical properties, fluid flow, and diffusion and reaction of disruptors, better cleaning strategies can be devised.</p> <p> </p> <p>Phase-field models, where the biofilm is treated like a viscoelastic fluid, are one of the few types of models that can predicting deformation and detachment based on mechanical properties. While several related studies have proposed phase-field models for predicting biofilm deformation, there has not been any validation of these models with experimental data. As a first step towards developing a material science strategy for biofilm control, this study validated the ability of a phase-field model to capture biofilm viscoelastic behavior.</p> <p> </p> <p>In this study, a two-dimensional continuum biofilm model was implemented with finite element method (FEM) using COMSOL Multiphysics (COMSOL v5.4, Comsol Inc, Burlington, MA). We applied the phase-field model with the Cahn-Hilliard equation to simulate biofilm mechanical behavior under fluid flow. The Oldroyd-B model, the simplest viscoelastic constitutive model, was applied to capture biofilm viscoelasticity. The biofilm was modeled as an incompressible viscoelastic fluid, with EPS and a water solvent. The phase-field physics were adapted from previous studies and applied to biofilm-fluid interactions. Two types of incompressible, immiscible fluids (EPS and water solvent) were studied as two components of a single fluid, with a fluid-fluid interface between the two.</p> <p> </p> <p>Homogeneous alginate was used as a synthetic biofilm for the experimental validation. The viscoelastic parameters of alginate were obtained by shear rheometry using stress relaxation tests. In experimental tests, the deformation behavior was observed in real time using optical coherence tomography (OCT). By importing the 2-D geometry from OCT and viscoelastic parameters from rheometry, the model was simulated and compared with real deformation in the flow cell.</p> <p> </p> <p>With the applied constant flow (Re=6), biofilm demonstrated viscoelastic behavior.The same behavior was observed in modeling as well. By tracking the movements of several locations of the biofilm geometry, it was concluded that the deformation of alginate biofilms was consistent with the computational results of phase-field models. The relative error between experiment and model for this certain location were 12.8%. Heterotrophic counter-diffusional biofilms cultured in membrane-aerated biofilm reactors were also tested in this study, with a relative error of 22.2%.</p> <p> </p> <p>In conclusion, the phase-field model, coupled with Oldroyd-B equation, could properly capture biofilm viscoelastic behavior. In a complex system, the phase-field model could be used as a tool to characterize the viscoelastic parameters from the observed deformation. With this information, the model can be used to predict the required disruptor dose to achieve high amounts of biofilm removal with a minimal amount of chemical addition. This can reduce operating costs and minimize the use of harsh chemicals.</p>


Author(s):  
Xinfu Chen ◽  
Gunduz Caginalp ◽  
Emre Esenturk

2001 ◽  
Vol 12 (1) ◽  
pp. 39-56 ◽  
Author(s):  
KARL GLASNER

Rapid solidification fronts are studied using a phase field model. Unlike slow moving solutions which approximate the Mullins–Sekerka free boundary problem, different limiting behaviour is obtained for rapidly moving fronts. A time-dependent analysis is carried out for various cases and the leading order behaviour of solidification front solutions is derived to be one of several travelling wave problems. An analysis of these problems is conducted, leading to expressions for front speeds in certain limits. The dynamics leading to these travelling wave solutions is derived, and conclusions about stability are drawn. Finally, a discussion is made of the relationship to other solidification models.


Author(s):  
Kais Ammar ◽  
Benoît Appolaire ◽  
Georges Cailletaud ◽  
Samuel Forest

A general constitutive framework is proposed to incorporate linear and nonlinear mechanical behaviour laws into a standard phase field model. In the diffuse interface region where both phases coexist, two mixture rules for strain and stress are introduced, which are based on the Voigt/Taylor and Reuss/Sachs well-known homogenization schemes and compared to the commonly used mixture rules in phase field models. Finite element calculations have been performed considering an elastoplastic precipitate growing in an elastic matrix in order to investigate the plastic accommodation processes.


2012 ◽  
Vol 736 ◽  
pp. 1-12 ◽  
Author(s):  
Rajdip Mukherjee ◽  
T.A. Abinandanan ◽  
M.P. Gururajan

Phase field models are widely used for the study of microstructures and their evolution. They can also be used as computer experiments. As computer experiments, they serve two important roles: (a) theoretical results which are hard to verify/validate experimentally can be verified/validated on the computer using phase field models; and, (b) when severe assumptions are made in a theory, they can be relaxed in the phase field model, and hence, results with wider reach can be obtained. In this paper, we discuss some such computer experiments in general, and the growth kinetics of precipitates in systems with tetragonal and cubic interfacial anisotropies in particular.


2021 ◽  
Vol 3 (9) ◽  
Author(s):  
Masato Kimura ◽  
Takeshi Takaishi ◽  
Sayahdin Alfat ◽  
Takumi Nakano ◽  
Yoshimi Tanaka

AbstractThree new industrial applications of irreversible phase field models for crack growth are presented in this study. The phase field model for irreversible crack growth in an elastic material is derived as a unidirectional gradient flow of the Francfort–Marigo energy with the Ambrosio–Tortorelli regularization, which is known to be consistent with classic Griffith fracture theory. The obtained compact parabolic-elliptic system of PDEs including two regularization parameters for space and time enables us to simulate various kinds of crack behaviors with standard finite element software, without any geometric restriction on the crack path. We extend the irreversible phase field model to thermal cracking in solder and to cracking in a viscoelastic material, keeping the compact forms of the PDEs and the energy consistency. On the other hand, for hydrogen-assisted cracking in metal, we propose a compact phase field model focusing on a kinematic jamming effect of the hydrogen by a weak coupling approach. Several numerical experiments for these three models show not only their reasonableness and usefulness but also flexible extendability of the phase field approach in fracture mechanics.


Sign in / Sign up

Export Citation Format

Share Document