Abstract:Phase interfaces in melting and solidification processes are strongly affected by the presence of convection in the liquid. One way of modeling their transient evolution is to couple an incompressible flow model to an energy balance in enthalpy formulation. Two strong nonlinearities arise, which account for the viscosity variation between phases and the latent heat of fusion at the phase interface. The resulting coupled system of PDE's can be solved by a single-domain semi-phase-field, variable viscosity, finite element method with monolithic system coupling and global Newton linearization. A robust computational model for realistic phase-change regimes furthermore requires a flexible implementation based on sophisticated mesh adaptivity. In this article, we present first steps towards implementing such a computational model into a simulation tool which we call Phaseflow. Phaseflow utilizes the finite element software FEniCS, which includes a dual-weighted residual method for goal-oriented adaptive mesh refinement. Phaseflow is an open-source, dimension-independent implementation that, upon an appropriate parameter choice, reduces to classical benchmark situations including the lid-driven cavity and the Stefan problem. We present and discuss numerical results for these, an octadecane PCM convection-coupled melting benchmark, and a preliminary 3D convection-coupled melting example, demonstrating the flexible implementation. Though being preliminary, the latter is, to our knowledge, the first published 3D result for this method. In our work, we especially emphasize reproducibility and provide an easy-to-use portable software container using Docker.

Computer simulations of a heterogeneous membrane with enhanced sampling techniques

Dynamic simulations of multicomponent lipid membranes over long length and time scales

TriMem: A parallelized hybrid Monte Carlo software for efficient simulations of lipid membranes

Identifying the Onset of Phase Separation in Quaternary Lipid Bilayer Systems from Coarse-Grained Simulations

Non-ideal mixing of lipids: a molecular dynamics perspective

Experimental validation of a phase-field model to predict coarsening dynamics of lipid domains in multicomponent membranes

Efficient Simulation of Tunable Lipid Assemblies Across Scales and Resolutions

Configurational Sampling of All-Atom Solvated Membranes Using Hybrid Nonequilibrium Molecular Dynamics Monte Carlo Simulations

Amphiphiles at Interfaces: Simulation of Structure and Phase Behavior

Phase separation and near-critical fluctuations in two-component lipid membranes: Monte Carlo simulations on experimentally relevant scales

Large-scale simulations of fluctuating biological membranes

Characterizing Microheterogeneity in Liquid Mixtures via Local Density Fluctuations

Lipid domain coarsening and fluidity in multicomponent lipid vesicles: A continuum based model and its experimental validation

Density Fluctuations, Solvation Thermodynamics and Coexistence Curves in Grand Canonical Molecular Dynamics Simulations

Binary bilayer simulations for partitioning within membranes

Efficient and realistic simulation of phase coexistence

Monolithic simulation of convection-coupled phase-change - verification and reproducibility

Conformational landscapes of membrane proteins delineated by enhanced sampling molecular dynamics simulations

Hybrid simulations of lateral diffusion in fluctuating membranes

Dynamic framework for large-scale modeling of membranes and peripheral proteins

A computational study of lateral phase separation in biological membranes