An operator splitting scheme for numerical simulation of spinodal decomposition and microstructure evolution of binary alloys
Abdullah Shah,
Sana Ayub,
Muhammad Sohaib,
Sadia Saeed,
Saher Akmal Khan,
Suhail Abbas,
Said Karim Shah
Affiliations
Abdullah Shah
Department of Mathematics, King Fahd University of Petroleum and Minerals, 31261, Dhahran, Saudi Arabia; Corresponding author.
Sana Ayub
Department of Mathematics, COMSATS University Islamabad, Park Road, Islamabad-45550, Pakistan
Muhammad Sohaib
Department of Mathematics, King Fahd University of Petroleum and Minerals, 31261, Dhahran, Saudi Arabia; Department of Mathematics, COMSATS University Islamabad, Park Road, Islamabad-45550, Pakistan; Department of Mathematics and Statistics, Bacha Khan University, Charsadda 24461, Pakistan
Sadia Saeed
Department of Applied Mathematics and Statistics, Institute of Space Technology, Islamabad, Pakistan
Saher Akmal Khan
Department of Mathematics, COMSATS University Islamabad, Park Road, Islamabad-45550, Pakistan
Suhail Abbas
Department of Mathematics, Karakoram International University, Ghizer Campus, Ghizer 15200, Pakistan
Said Karim Shah
Department of Physics, Abdul Wali Khan University Mardan, 23200, Khyber Pakhtunkhwa, Pakistan
This article compares the operator splitting scheme to linearly stabilized splitting and semi-implicit Euler's schemes for the numerical solution of the Cahn-Hilliard equation. For the purpose of validation, the spinodal decomposition phenomena have been simulated. The efficacy of the three schemes has been demonstrated through numerical experiments. The computed results show that the schemes are conditionally stable. It has been observed that the operator splitting scheme is computationally more efficient.
