Solving Allen-Cahn equations with periodic and nonperiodic boundary conditions using mimetic finite-difference
Saulo Orizaga1, Gilberto González-Parra1,2, Logan Forman1
1Department of Mathematics, New Mexico Tech, Socorro, NM, USA.
None:
In this paper, we investigate and implement a numerical method that is based on the mimetic finite difference operator in order to solve the nonlinear Allen-Cahn equation with periodic and non-periodic boundary conditions. In addition, we also analyze the performance of this mimetic-based method by using the classical heat equation with a variety of boundary conditions. We assess the performance of the mimetic-based numerical method by comparing the errors of its solutions with those obtained by a classical finite difference method and the pdepde built-in Matlab function. We compute the errors by using the exact solutions when they are available or with reference solutions. We adapt and implement the mimetic-based numerical method by using the MOLE (Mimetic Operators Library Enhanced) library that includes some built-in functions that return representations of the curl, divergence and gradient operators, in order to deal with the Allen-Cahn and heat equations. We present several results with regard to errors and numerical convergence tests in order to provide insight into the accuracy of the mimetic-based numerical method. The results show that the numerical method based on the mimetic difference operator is a reliable method for solving the Allen-Cahn and heat equations with periodic and non-periodic boundary conditions. The numerical solutions generated by the mimetic-based method are relatively accurate. We also proposed a new method based on the mimetic finite difference operator and the convexity splitting approach to solve Allen-Cahn equation in 2D. We found that, for small time step sizes the solutions generated by the mimetic-based method are more accurate than the ones generated by the pdepe Matlab function and similar to the solutions given by a finite difference method.
Related Concept Videos
Transmission-Line Differential Equations
Line Section Model
A circuit representing a line section of length Δx helps in understanding the transmission line parameters. The voltage V(x) and current i(x) are measured...
Linear Approximation in Time Domain
For a simple pendulum with a mass evenly distributed along its length and the center of mass located at half the pendulum's length,...
RLC Circuit as a Damped Oscillator
Consider a series RLC circuit. Here, the presence of resistance in the circuit leads to energy loss due to joule heating in the resistance. Therefore, the total electromagnetic energy in the circuit is no longer constant and decreases with time. Since the magnitude of charge, current, and potential difference continuously decreases, their oscillations are said to be damped. This is...
Linear Approximation in Frequency Domain
In contrast, nonlinear systems do not inherently possess these properties. However, for small deviations around an operating point, a nonlinear system can often be approximated as linear....
Difference Equation Solution using z-Transform
The z-transform facilitates handling delayed signals by shifting the signal in the z-domain, which corresponds to delaying the signal in the time domain, and advancing signals by similarly shifting in the...
Magnetostatic Boundary Conditions


