scholarly journals Improved phase-field models of melting and dissolution in multi-component flows

Author(s):  
Eric W. Hester ◽  
Louis-Alexandre Couston ◽  
Benjamin Favier ◽  
Keaton J. Burns ◽  
Geoffrey M. Vasil

We develop and analyse the first second-order phase-field model to combine melting and dissolution in multi-component flows. This provides a simple and accurate way to simulate challenging phase-change problems in existing codes. Phase-field models simplify computation by describing separate regions using a smoothed phase field. The phase field eliminates the need for complicated discretizations that track the moving phase boundary. However, standard phase-field models are only first-order accurate. They often incur an error proportional to the thickness of the diffuse interface. We eliminate this dominant error by developing a general framework for asymptotic analysis of diffuse-interface methods in arbitrary geometries. With this framework, we can consistently unify previous second-order phase-field models of melting and dissolution and the volume-penalty method for fluid–solid interaction. We finally validate second-order convergence of our model in two comprehensive benchmark problems using the open-source spectral code Dedalus.

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.


2019 ◽  
Vol 2019 ◽  
pp. 1-10
Author(s):  
Darae Jeong ◽  
Yibao Li ◽  
Chaeyoung Lee ◽  
Junxiang Yang ◽  
Yongho Choi ◽  
...  

In this paper, we propose a verification method for the convergence rates of the numerical solutions for parabolic equations. Specifically, we consider the numerical convergence rates of the heat equation, the Allen–Cahn equation, and the Cahn–Hilliard equation. Convergence test results show that if we refine the spatial and temporal steps at the same time, then we have the second-order convergence rate for the second-order scheme. However, in the case of the first-order in time and the second-order in space scheme, we may have the first-order or the second-order convergence rates depending on starting spatial and temporal step sizes. Therefore, for a rigorous numerical convergence test, we need to perform the spatial and the temporal convergence tests separately.


2016 ◽  
Vol 4 (1) ◽  
Author(s):  
Zhan Chen

AbstractIn this work, the existence of a global minimizer for the previous Lagrangian formulation of nonpolar solvation model proposed in [1] has been proved. One of the proofs involves a construction of a phase field model that converges to the Lagrangian formulation. Moreover, an Eulerian formulation of nonpolar solvation model is proposed and implemented under a similar parameterization scheme to that in [1]. By doing so, the connection, similarity and difference between the Eulerian formulation and its Lagrangian counterpart can be analyzed. It turns out that both of them have a great potential in solvation prediction for nonpolar molecules, while their decompositions of attractive and repulsive parts are different. That indicates a distinction between phase field models of solvation and our Eulerian formulation.


2017 ◽  
Vol 44 (2) ◽  
pp. 189-214 ◽  
Author(s):  
M. Dehsara ◽  
H. Fu ◽  
S.Dj Mesarovic ◽  
D.P. Sekulic ◽  
M. Krivilyov

Phase field (diffuse interface) models accommodate diffusive triple line motion with variable contact angle, thus allowing for the no-slip boundary condition without the stress singularities. We consider two commonly used classes of phase field models: the compositionally compressible (CC) model with compressibility limited to the fluid mix within the diffuse interface, and the incompressible (IC) model. First, we show that the CC model applied to fluids with dissimilar mass densities exhibits the computational instability leading to the breakup of the triple line. We provide a qualitative physical explanation of this instability and argue that the compositional compressibility within the diffuse interface is inconsistent with the global incompressible flow. Second, we derive the IC model as a systematic approximation to the CC model, based on a suitable choice of continuum velocity field. Third, we benchmark the IC model against sharp interface theory and experimental kinetics. The triple line kinetics is well represented by the triple line mobility parameter. Finally, we investigate the effects of the bulk phase field diffusional mobility parameter on the kinetics of the wetting process and find that within a wide range of magnitudes the bulk mobility does not affect the flow.


2020 ◽  
pp. 2000162
Author(s):  
Chao Yang ◽  
Houbing Huang ◽  
Wenbo Liu ◽  
Junsheng Wang ◽  
Jing Wang ◽  
...  

Author(s):  
Naoki Takada

For interface-tracking simulation of two-phase flows in various micro-fluidics devices, the applicability of two versions of Navier-Stokes phase-field method (NS-PFM) was examined, combining NS equations for a continuous fluid with a diffuse-interface model based on the van der Waals-Cahn-Hilliard free-energy theory. Through the numerical simulations, the following major findings were obtained: (1) The first version of NS-PFM gives good predictions of interfacial shapes and motions in an incompressible, isothermal two-phase fluid with high density ratio on solid surface with heterogeneous wettability. (2) The second version successfully captures liquid-vapor motions with heat and mass transfer across interfaces in phase change of a non-ideal fluid around the critical point.


Author(s):  
Arne Claus Hansen-Dörr ◽  
Jörg Brummund ◽  
Markus Kästner

Abstract In this contribution, a variational diffuse modeling framework for cracks in heterogeneous media is presented. A static order parameter smoothly bridges the discontinuity at material interfaces, while an evolving phase-field captures the regularized crack. The key novelty is the combination of a strain energy split with a partial rank-I relaxation in the vicinity of the diffuse interface. The former is necessary to account for physically meaningful crack kinematics like crack closure, the latter ensures the mechanical jump conditions throughout the diffuse region. The model is verified by a convergence study, where a circular bi-material disc with and without a crack is subjected to radial loads. For the uncracked case, analytical solutions are taken as reference. In a second step, the model is applied to crack propagation, where a meaningful influence on crack branching is observed, that underlines the necessity of a reasonable homogenization scheme. The presented model is particularly relevant for the combination of any variational strain energy split in the fracture phase-field model with a diffuse modeling approach for material heterogeneities.


Author(s):  
T. Philippe ◽  
H. Henry ◽  
M. Plapp

At equilibrium, the shape of a strongly anisotropic crystal exhibits corners when for some orientations the surface stiffness is negative. In the sharp-interface problem, the surface free energy is traditionally augmented with a curvature-dependent term in order to round the corners and regularize the dynamic equations that describe the motion of such interfaces. In this paper, we adopt a diffuse interface description and present a phase-field model for strongly anisotropic crystals that is regularized using an approximation of the Willmore energy. The Allen–Cahn equation is employed to model kinetically controlled crystal growth. Using the method of matched asymptotic expansions, it is shown that the model converges to the sharp-interface theory proposed by Herring. Then, the stress tensor is used to derive the force acting on the diffuse interface and to examine the properties of a corner at equilibrium. Finally, the coarsening dynamics of the faceting instability during growth is investigated. Phase-field simulations reveal the existence of a parabolic regime, with the mean facet length evolving in t , with t the time, as predicted by the sharp-interface theory. A specific coarsening mechanism is observed: a hill disappears as the two neighbouring valleys merge.


Author(s):  
Rakesh Dhote ◽  
Kamran Behdinan

In this paper, we study the dynamic thermo-mechanical behaviors of 3D shape memory alloy (SMA) nanostructures using the phase-field (PF) model. The PF model is based on the Ginzburg-Landau theory and requires a non-convex free energy function for an adequate description of the cubic-to-tetragonal martensitic phase transformations. We have developed a model that includes domain walls, treated as a diffuse interface, which leads to a fourth-order differential equation in a strain-based order parameter PF model. Arising numerical challenges have been overcome based on an isogeometric analysis (IGA) framework. Microstructure morphology evolution and consequent thermo-mechanical properties have been studied on SMA nanostructures of different geometries. The numerical results are in agreement with experimental observations. The developed coupled dynamic model has provided a better understanding of underlying microstructures and behaviors, which can be used for development of better SMA-based devices.


Sign in / Sign up

Export Citation Format

Share Document