Related papers: An efficient numerical method for the anisotropic …
It is well known that the seven-step backward difference formula (BDF) is unstable for the parabolic equations, since it is not even zero-stable. However, a linear combination of two non zero-stable schemes, namely the seven-step BDF and…
We apply a simple dynamical density functional theory, the phase-field crystal (PFC) model of overdamped conservative dynamics, to address polymorphism, crystal nucleation, and crystal growth in the diffusion-controlled limit. We refine the…
We consider the numerical approximations for a phase field model consisting of incompressible Navier--Stokes equations with a generalized Navier boundary condition, and the Cahn-Hilliard equation with a dynamic moving contact line boundary…
In this paper, we consider the numerical solution of a binary fluid-surfactant phase field model, in which the free energy contains a nonlinear coupling entropy, a Ginzburg-Landau double well potential, and a logarithmic Flory-Huggins…
Interface energy and kinetic coefficient of crystal growth strongly depend on the face of the crystalline lattice. To investigate the kinetic anisotropy and velocity of different crystallographic faces we use the hyperbolic (modified) phase…
We propose an efficient and accurate parametric finite element method (PFEM) for solving sharp-interface continuum models for solid-state dewetting of thin films with anisotropic surface energies. The governing equations of the…
We present a numerical method for approximating the solutions of degenerate parabolic equations with a formal gradient flow structure. The numerical method we propose preserves at the discrete level the formal gradient flow structure,…
In this paper, we present a high order finite difference solver for anisotropic diffusion problems based on the first-order hyperbolic system method. In particular, we demonstrate that the construction of a uniformly accurate fifth-order…
We propose some finite element schemes to solve a class of fourth-order nonlinear PDEs, which include the vector-valued Landau--Lifshitz--Baryakhtar equation, the Swift--Hohenberg equation, and various Cahn--Hilliard-type equations with…
Facetted growth of snow crystals leads to a rich diversity of forms, and exhibits a remarkable sixfold symmetry. Snow crystal structures result from diffusion limited crystal growth in the presence of anisotropic surface energy and…
A nonlinear diffusion equation, interpreted as a Wasserstein gradient flow, is numerically solved in one space dimension using a higher-order minimizing movement scheme based on the BDF (backward differentiation formula) discretization. In…
We study dendritic microstructure evolution using an adaptive grid, finite element method applied to a phase-field model. The computational complexity of our algorithm, per unit time, scales linearly with system size, rather than the…
In this paper, we compute the stationary states of the multicomponent phase-field crystal model by formulating it as a block constrained minimization problem. The original infinite-dimensional non-convex minimization problem is approximated…
In this paper, we propose a numerical method to approximate the solution of partial differential equations in irregular domains with no-flux boundary conditions by means of spectral methods. The main features of this method are its…
In this paper, we present two multiple scalar auxiliary variable (MSAV)-based, finite element numerical schemes for the Abels-Garcke-Gr{\"u}n (AGG) model, which is a thermodynamically consistent phase field model of two-phase incompressible…
We consider an anisotropic version of Baxter's model of `sticky hard spheres', where a nonuniform adhesion is implemented by adding, to an isotropic surface attraction, an appropriate `dipolar sticky' correction (positive or negative,…
The numerical analysis of stochastic time fractional evolution equations presents considerable challenges due to the limited regularity of the model caused by the nonlocal operator and the presence of noise. The existing time-stepping…
We investigate how to simulate anisotropic heat conduction in a stable manner in Smoothed Particle Hydrodynamics. We show that the requirement for stability is that entropy must increase. From this, we deduce that methods involving direct…
High-order time-stepping schemes are crucial for simulating incompressible fluid flows due to their ability to capture complex turbulent behavior and unsteady motion. In this work, we propose a third-order accurate numerical scheme for the…
In this paper, we develop a first order (in time) numerical scheme for the binary fluid surfactant phase field model. The free energy contains a double-well potential, a nonlinear coupling entropy and a Flory-Huggins potential. The…