Structure preserving integration and model order reduction of skew-gradient reaction-diffusion systems

Karasözen, Bülent
Uzunca, Murat
Activator-inhibitor FitzHugh-Nagumo (FHN) equation is an example for reaction-diffusion equations with skew-gradient structure. We discretize the FHN equation using symmetric interior penalty discontinuous Galerkin (SIPG) method in space and average vector field (AVF) method in time. The AVF method is a geometric integrator, i.e. it preserves the energy of the Hamiltonian systems and energy dissipation of the gradient systems. In this work, we show that the fully discrete energy of the FHN equation satisfies the mini-maximizer property of the continuous energy for the skew-gradient systems. We present numerical results with traveling fronts and pulses for one dimensional, two coupled FHN equations and three coupled FHN equations with one activator and two inhibitors in skew-gradient form. Turing patterns are computed for fully discretized two dimensional FHN equation in the form of spots and labyrinths. Because the computation of the Turing patterns is time consuming for different parameters, we applied model order reduction with the proper orthogonal decomposition (POD). The nonlinear term in the reduced equations is computed using the discrete empirical interpolation (DEIM) with SIPG discretization. Due to the local nature of the discontinuous Galerkin method, the nonlinear terms can be computed more efficiently than for the continuous finite elements. The reduced solutions are very close to the fully discretized ones. The efficiency and accuracy of the POD and POD-DEIM reduced solutions are shown for the labyrinth-like patterns.


Energy Stable Discontinuous Galerkin Finite Element Method for the Allen-Cahn Equation
Karasözen, Bülent; Sariaydin-Filibelioglu, Ayse; Yücel, Hamdullah (2018-05-01)
In this paper, we investigate numerical solution of Allen-Cahn equation with constant and degenerate mobility, and with polynomial and logarithmic energy functionals. We discretize the model equation by symmetric interior penalty Galerkin (SIPG) method in space, and by average vector field (AVF) method in time. We show that the energy stable AVF method as the time integrator for gradient systems like the Allen-Cahn equation satisfies the energy decreasing property for fully discrete scheme. Numerical result...
Structure-preserving reduced-order modeling of Korteweg–de Vries equation
Uzunca, Murat; Karasözen, Bülent; Yıldız, Süleyman (2021-10-01)
Computationally efficient, structure-preserving reduced-order methods are developed for the Korteweg–de Vries (KdV) equations in Hamiltonian form. The semi-discretization in space by finite differences is based on the Hamiltonian structure. The resulting skew-gradient system of ordinary differential equations (ODEs) is integrated with the linearly implicit Kahan's method, which preserves the Hamiltonian approximately. We have shown, using proper orthogonal decomposition (POD), the Hamiltonian structure of t...
Exact Solutions of Effective-Mass Dirac-Pauli Equation with an Electromagnetic Field
Arda, Altug; Sever, Ramazan (Springer Science and Business Media LLC, 2017-01-01)
The exact bound state solutions of the Dirac-Pauli equation are studied for an appropriate position-dependent mass function by using the Nikiforov-Uvarov method. For a central electric field having a shifted inverse linear term, all two kinds of solutions for bound states are obtained in closed forms.
Structure preserving model order reduction of shallow water equations
Karasözen, Bülent; UZUNCA, MURAT (2020-07-01)
In this paper, we present two different approaches for constructing reduced-order models (ROMs) for the two-dimensional shallow water equation (SWE). The first one is based on the noncanonical Hamiltonian/Poisson form of the SWE. After integration in time by the fully implicit average vector field method, ROMs are constructed with proper orthogonal decomposition(POD)/discrete empirical interpolation method that preserves the Hamiltonian structure. In the second approach, the SWE as a partial differential eq...
Numerical solution of nonlinear reaction-diffusion and wave equations
Meral, Gülnihal; Tezer, Münevver; Department of Mathematics (2009)
In this thesis, the two-dimensional initial and boundary value problems (IBVPs) and the one-dimensional Cauchy problems defined by the nonlinear reaction- diffusion and wave equations are numerically solved. The dual reciprocity boundary element method (DRBEM) is used to discretize the IBVPs defined by single and system of nonlinear reaction-diffusion equations and nonlinear wave equation, spatially. The advantage of DRBEM for the exterior regions is made use of for the latter problem. The differential quad...
Citation Formats
B. Karasözen and M. Uzunca, “Structure preserving integration and model order reduction of skew-gradient reaction-diffusion systems,” ANNALS OF OPERATIONS RESEARCH, pp. 79–106, 2017, Accessed: 00, 2020. [Online]. Available: