Numerical scheme for solving the nonuniformly forced cubic and quintic Swift-Hohenberg equations strictly respecting the Lyapunov functional

DOI10.1016/J.CAM.2021.114005zbMATH Open1503.65166arXiv2007.16080OpenAlexW3046883219MaRDI QIDQ2075950FDOQ2075950

Authors: D. L. Coelho, E. Vitral, José Pedro Pontes, N. Mangiavacchi 

Publication date: 16 February 2022

Published in: Journal of Computational and Applied Mathematics (Search for Journal in Brave)

Abstract: Computational modeling of pattern formation in nonequilibrium systems is a fundamental tool for studying complex phenomena in biology, chemistry, materials science and engineering. The pursuit for theoretical descriptions of some among those physical problems led to the Swift-Hohenberg equation (SH3) which describes pattern selection in the vicinity of instabilities. A finite differences scheme, known as Stabilizing Correction (Christov & Pontes; 2001 DOI: 10.1016/S0895-7177(01)00151-0), developed to integrate the cubic Swift-Hohenberg equation in two dimensions, is reviewed and extended in the present paper. The original scheme features Generalized Dirichlet boundary conditions (GDBC), forcings with a spatial ramp of the control parameter, strict implementation of the associated Lyapunov functional, and second order representation of all derivatives. We now extend these results by including periodic boundary conditions (PBC), forcings with gaussian distributions of the control parameter and the quintic Swift-Hohenberg (SH35) model. The present scheme also features a strict implementation of the functional for all test cases. A code verification was accomplished, showing unconditional stability, along with second order accuracy in both time and space. Test cases confirmed the monotonic decay of the Lyapunov functional and all numerical experiments exhibit the main physical features: highly nonlinear behaviour, wavelength filter and competition between bulk and boundary effects.

Full work available at URL: https://arxiv.org/abs/2007.16080

Recommendations

• Numerical scheme for Swift-Hohenberg equation with strict implementation of Lyapunov functional.
• A new space-time discretization for the Swift-Hohenberg equation that strictly respects the Lyapunov functional
• Numerical simulation and analysis of the Swift-Hohenberg equation by the stabilized Lagrange multiplier approach
• Numerical simulation of Swift-Hohenberg equation by the fourth-order compact scheme
• A semi-analytical Fourier spectral method for the Swift-Hohenberg equation

zbMATH Keywords

nonlinear systemspattern formationSwift-Hohenberg equationfinite difference methods

Mathematics Subject Classification ID

Finite difference methods applied to problems in fluid mechanics (76M20) Homoclinic and heteroclinic orbits for dynamical systems (37C29) Finite difference methods for initial value and initial-boundary value problems involving PDEs (65M06) Stability and convergence of numerical methods for initial value and initial-boundary value problems involving PDEs (65M12) Free convection (76R10)

Cites Work

• On the Numerical Solution of Heat Conduction Problems in Two and Three Space Variables
• Instabilities and self-organisation in materials. Vol. 1: Fundamentals of nanoscience. Vol. 2: Applications in materials design and nanotechnology.
• Pattern formation outside of equilibrium
• Review of code and solution verification procedures for computational simulation
• To Snake or Not to Snake in the Planar Swift–Hohenberg Equation
• Title not available (Why is that?)
• Implicit time splitting for fourth-order parabolic equations
• Numerical scheme for Swift-Hohenberg equation with strict implementation of Lyapunov functional.
• Multiscale finite-difference-diffusion-Monte-Carlo method for simulating dendritic solidification
• An energy stable method for the Swift-Hohenberg equation with quadratic-cubic nonlinearity
• On dislocation patterning: multiple slip effects in the rate equation approach
• A simple and efficient finite difference method for the phase-field crystal equation on curved surfaces
• Boundary conditons on the envelope function of convective rolls close to onset
• On the Onset of Rayleigh-Bénard Convection in a Fluid Layer of Slowly Increasing Depth
• A simple and efficient scheme for phase field crystal simulation
• An efficient and stable compact fourth-order finite difference scheme for the phase field crystal equation
• Stripe patterns orientation resulting from nonuniform forcings and other competitive effects in the Swift-Hohenberg dynamics
• The onset of cellular convection in a shallow two-dimensional container of fluid heated non-uniformly from below

Cited In (5)

• A new space-time discretization for the Swift-Hohenberg equation that strictly respects the Lyapunov functional
• Rigorous numerics for localized patterns to the quintic Swift-Hohenberg equation
• Numerical scheme for Swift-Hohenberg equation with strict implementation of Lyapunov functional.
• A non-iterative and unconditionally energy stable method for the Swift-Hohenberg equation with quadratic-cubic nonlinearity
• Linear and energy stable schemes for the Swift-Hohenberg equation with quadratic-cubic nonlinearity based on a modified scalar auxiliary variable approach

This page was built for publication: Numerical scheme for solving the nonuniformly forced cubic and quintic Swift-Hohenberg equations strictly respecting the Lyapunov functional

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q2075950)