scholarly journals Minimization and Eulerian Formulation of Differential Geormetry Based Nonpolar Multiscale Solvation Models

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.

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.


2013 ◽  
Vol 13 (4) ◽  
pp. 1045-1065 ◽  
Author(s):  
Jie Shen ◽  
Xiaofeng Yang ◽  
Qi Wang

AbstractThe commonly used incompressible phase field models for non-reactive, binary fluids, in which the Cahn-Hilliard equation is used for the transport of phase variables (volume fractions), conserve the total volume of each phase as well as the material volume, but do not conserve the mass of the fluid mixture when densities of two components are different. In this paper, we formulate the phase field theory for mixtures of two incompressible fluids, consistent with the quasi-compressible theory [28], to ensure conservation of mass and momentum for the fluid mixture in addition to conservation of volume for each fluid phase. In this formulation, the mass-average velocity is no longer divergence-free (solenoidal) when densities of two components in the mixture are not equal, making it a compressible model subject to an internal con-straint. In one formulation of the compressible models with internal constraints (model 2), energy dissipation can be clearly established. An efficient numerical method is then devised to enforce this compressible internal constraint. Numerical simulations in confined geometries for both compressible and the incompressible models are carried out using spatially high order spectral methods to contrast the model predictions. Numerical comparisons show that (a) predictions by the two models agree qualitatively in the situation where the interfacial mixing layer is thin; and (b) predictions differ significantly in binary fluid mixtures undergoing mixing with a large mixing zone. The numerical study delineates the limitation of the commonly used incompressible phase field model using volume fractions and thereby cautions its predictive value in simulating well-mixed binary fluids.


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):  
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.


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.


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):  
yangxin Tang ◽  
Wei Gao

In this paper we present the phase-field models to describe nonisothermal solidification of ideal multicomponent and multiphase alloy systems. Governing equations are developed for the temporal and spatial variation of three phase-field functions, as well as the temperature field. The global existence of weak solutions to parabolic differential equations in three dimension was proved by the Galerkin method. The existence of a maximum theorem are also extensively studied.


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.


Sign in / Sign up

Export Citation Format

Share Document