Related papers: Method of lines transpose: Energy gradient flows u…
The $L^2$ gradient flow of the Ginzburg-Landau free energy functional leads to the Allen Cahn equation that is widely used for modeling phase separation. Machine learning methods for solving the Allen-Cahn equation in its strong form suffer…
We present Cahn-Hilliard and Allen-Cahn numerical integration algorithms that are unconditionally stable and so provide significantly faster accuracy-controlled simulation. Our stability analysis is based on Eyre's theorem and unconditional…
An efficient direct solver for volume integral equations with O(N) complexity for a broad range of problems is presented. The solver relies on hierarchical compression of the discretized integral operator, and exploits that off-diagonal…
We introduce and analyze various Regularized Combined Field Integral Equations (CFIER) formulations of time-harmonic Navier equations in media with piece-wise constant material properties. These formulations can be derived systematically…
This work aims to construct an efficient and highly accurate numerical method to address the time singularity at $t=0$ involved in a class of time-fractional parabolic integro-partial differential equations in one and two dimensions. The…
In this paper, we propose a variable time-step linear relaxation scheme for time-fractional phase-field equations with a free energy density in general polynomial form. The $L1^{+}$-CN formula is used to discretize the fractional…
We present a novel numerical method and algorithm for the solution of the 3D axially symmetric time-dependent Schr\"odinger equation in cylindrical coordinates, involving singular Coulomb potential terms besides a smooth time-dependent…
This study proposes an algorithm for modeling compressible flows in spherical shells in nearly incompressible and weakly compressible regimes based on an implicit direction splitting approach. The method retains theoretically expected…
We present a numerical scheme for solving a sixth-order Cahn-Hilliard type equation that captures the dynamics of phase transitions in a ternary mixture consisting of two immiscible fluids and a surface active molecule that is amphiphilic.…
We present a second-order-in-time finite difference scheme for the Cahn-Hilliard-Hele-Shaw equations. This numerical method is uniquely solvable and unconditionally energy stable. At each time step, this scheme leads to a system of…
Time-implicit schemes are attractive since they allow numerical time steps that are much larger than those permitted by the Courant-Friedrich-Lewy criterion characterizing time-explicit methods. This advantage comes, however, with a cost:…
In this paper we analyze a fully discrete scheme for a general Cahn-Hilliard equation coupled with a nonsteady Magneto-hydrodynamics flow, which describes two immiscible, incompressible and electrically conducting fluids with different…
In the present work, a high order finite element type residual distribution scheme is designed in the framework of multidimensional compressible Euler equations of gas dynamics. The strengths of the proposed approximation rely on the…
A novel numerical strategy is introduced for computing approximations of solutions to a Cahn-Hilliard model with degenerate mobilities. This model has recently been introduced as a second-order phase-field approximation for surface…
In this paper, we propose a quadratic reformulation theory for rational-like functions. Based on this theory, we develop the Quadratic Conservation Elevation (QCE) method, which combines the Scalar Auxiliary Variable (SAV) method with the…
This paper describes the first steps of development of a new multidimensional time implicit code devoted to the study of hydrodynamical processes in stellar interiors. The code solves the hydrodynamical equations in spherical geometry and…
In this article we present a refined convergence analysis for a second order accurate in time, fourth order finite difference numerical scheme for the 3-D Cahn-Hilliard equation, with an improved convergence constant. A modified backward…
A numerical method for solving elliptic PDEs with variable coefficients on two-dimensional domains is presented. The method is based on high-order composite spectral approximations and is designed for problems with smooth solutions. The…
The discrete gradient structure and the positive definiteness of discrete fractional integrals or derivatives are fundamental to the numerical stability in long-time simulation of nonlinear integro-differential models. We build up a…
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…