Tuğba Küçükseyhan

Murat Uzunca, Tuğba Küçükseyhan, Hamdullah Yücel et al.

We investigate smooth and sparse optimal control problems for convective FitzHugh-Nagumo equation with travelling wave solutions in moving excitable media. The cost function includes distributed space-time and terminal observations or targets. The state and adjoint equations are discretized in space by symmetric interior point Galerkin (SIPG) method and by backward Euler method in time. Several numerical results are presented for the control of the travelling waves. We also show numerically the validity of the second order optimality conditions for the local solutions of the sparse optimal control problem for vanishing Tikhonov regularization parameter. Further, we estimate the distance between the discrete control and associated local optima numerically by the help of the perturbation method and the smallest eigenvalue of the reduced Hessian.

Murat Uzunca, Bülent Karasözen, Tuğba Küçükseyhan

In this paper, a moving mesh discontinuous Galerkin (dG) method is developed for nonlinear partial differential equations (PDEs) with traveling wave solutions. The moving mesh strategy for one dimensional PDEs is based on the rezoning approach which decouples the solution of the PDE from the moving mesh equation. We show that the dG moving mesh method is able to resolve sharp wave fronts and wave speeds accurately for the optimal, arc-length and curvature monitor functions. Numerical results reveal the efficiency of the proposed moving mesh dG method for solving Burgers', Burgers'-Fisher and Schlögl(Nagumo) equations.

Bülent Karasözen, Murat Uzunca, Tuğba Küçükseyhan

We developed a reduced order model (ROM) using the proper orthogonal decomposition (POD) to compute efficiently the labyrinth and spot like patterns of the FitzHugh-Nagumo (FNH) equation. The FHN equation is discretized in space by the discontinuous Galerkin (dG) method and in time by the backward Euler method. Applying POD-DEIM (discrete empirical interpolation method) to the full order model (FOM) for different values of the parameter in the bistable nonlinearity, we show that using few POD and DEIM modes, the patterns can be computed accurately. Due to the local nature of the dG discretization, the POD-DEIM requires less number of connected nodes than continuous finite element for the nonlinear terms, which leads to a significant reduction of the computational cost for dG POD-DEIM.

Bülent Karasözen, Tuğba Küçükseyhan, Murat Uzunca

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 (DG) 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.