Source author record

Murat Uzunca

Murat Uzunca appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

13works
3topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

13 published item(s)

preprint2022arXiv

Global energy preserving model reduction for multi-symplectic PDEs

Many Hamiltonian systems can be recast in multi-symplectic form. We develop a reduced-order model (ROM) for multi-symplectic Hamiltonian partial differential equations (PDEs) that preserves the global energy. The full-order solutions are obtained by finite difference discretization in space and the global energy preserving average vector field (AVF) method. The ROM is constructed in the same way as the full-order model (FOM) applying proper orthogonal decomposition (POD) with the Galerkin projection. The reduced-order system has the same structure as the FOM, and preserves the discrete reduced global energy. Applying the discrete empirical interpolation method (DEIM), the reduced-order solutions are computed efficiently in the online stage. A priori error bound is derived for the DEIM approximation to the nonlinear Hamiltonian. The accuracy and computational efficiency of the ROMs are demonstrated for the Korteweg de Vries (KdV) equation, Zakharov-Kuznetzov (ZK) equation, and nonlinear Schr{ö}dinger (NLS) equation in multi-symplectic form. Preservation of the reduced energies shows that the reduced-order solutions ensure the long-term stability of the solutions.

preprint2021arXiv

Nonintrusive model order reduction for cross-diffusion systems

In this paper, we investigate tensor based nonintrusive reduced-order models (ROMs) for parametric cross-diffusion equations. The full-order model (FOM) consists of ordinary differential equations (ODEs) in matrix or tensor form resulting from finite-difference discretization of the differential operators by taking the advantage of Kronecker structure. The matrix/tensor differential equations are integrated in time with the implicit-explicit (IMEX) Euler method. The reduced bases, relying on a finite sample set of parameter values, are constructed in form of a two-level approach by applying higher-order singular value decomposition (HOSVD) to the space-time snapshots in tensor form, which leads to a large amount of computational and memory savings. The nonintrusive reduced approximations for an arbitrary parameter value are obtained through tensor product of the reduced basis by the parameter dependent core tensor that contains the reduced coefficients. The reduced coefficients for new parameter values are computed using radial basis function (RBF) interpolation. The efficiency of the proposed method is illustrated through numerical experiments for two-dimensional Schnakenberg and three-dimensional Brusselator cross-diffusion equations. The spatiotemporal patterns are accurately predicted by the reduced-order models with speed-up factors of orders two and three over the full-order models.

preprint2021arXiv

Reduced order modelling of nonlinear cross-diffusion systems

In this work, we present a reduced-order model for a nonlinear cross-diffusion problem from population dynamics, for the Shigesada-Kawasaki-Teramoto (SKT) equation with Lotka-Volterra kinetics. The finite-difference discretization of the SKT equation in space results in a system of linear--quadratic ordinary differential equations (ODEs). The reduced order model (ROM) has the same linear-quadratic structure as the full order model (FOM). Using the linear-quadratic structure of the ROM, the reduced-order solutions are computed independent of the full solutions with the proper orthogonal decomposition (POD). The computation of the reduced solutions is further accelerated by applying tensorial POD. The formation of the patterns of the SKT equation consists of a fast transient phase and a long steady-state phase. Reduced order solutions are computed by separating the time, into two-time intervals. In numerical experiments, we show for one-and two-dimensional SKT equations with pattern formation, the reduced-order solutions obtained in the time-windowed form, i.e., principal decomposition framework (P-POD), are more accurate than the global POD solutions (G-POD) obtained in the whole time interval. Furthermore, we show the decrease of the entropy numerically by the reduced solutions, which is important for the global existence of nonlinear cross-diffusion equations such as the SKT equation.

preprint2020arXiv

Adaptive Discontinuous Galerkin Finite Elements for Advective Allen-Cahn Equation

We apply a space adaptive interior penalty discontinuous Galerkin method for solving advective Allen-Cahn equation with expanding and contracting velocity fields. The advective Allen-Cahn equation is first discretized in time and the resulting semi-linear elliptic PDE is solved by an adaptive algorithm using a residual-based a posteriori error estimator. The a posteriori error estimator contains additional terms due to the non-divergence-free velocity field. Numerical examples demonstrate the effectiveness and accuracy of the adaptive approach by resolving the sharp layers accurately.

preprint2020arXiv

Pricing European and American Options under Heston Model using Discontinuous Galerkin Finite Elements

This paper deals with pricing of European and American options, when the underlying asset price follows Heston model, via the interior penalty discontinuous Galerkin finite element method (dGFEM). The advantages of dGFEM space discretization with Rannacher smoothing as time integrator with nonsmooth initial and boundary conditions are illustrated for European vanilla options, digital call and American put options. The convection dominated Heston model for vanishing volatility is efficiently solved utilizing the adaptive dGFEM. For fast solution of the linear complementary problem of the American options, a projected successive over relaxation (PSOR) method is developed with the norm preconditioned dGFEM. We show the efficiency and accuracy of dGFEM for option pricing by conducting comparison analysis with other methods and numerical experiments.

preprint2020arXiv

Structure Preserving Model Order Reduction of Shallow Water Equations

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/discrete empirical interpolation method (POD/DEIM) that preserves the Hamiltonian structure. In the second approach, the SWE as a partial differential equation with quadratic nonlinearity is integrated in time by the linearly implicit Kahan's method and ROMs are constructed with the tensorial POD that preserves the linear-quadratic structure of the SWE. We show that in both approaches, the invariants of the SWE such as the energy, enstrophy, mass, and circulation are preserved over a long period of time, leading to stable solutions. We conclude by demonstrating the accuracy and the computational efficiency of the reduced solutions by a numerical test problem.

preprint2019arXiv

Reduced Order Optimal Control of the Convective FitzHugh-Nagumo Equation

In this paper, we compare three model order reduction methods: the proper orthogonal decomposition (POD), discrete empirical interpolation method (DEIM) and dynamic mode decomposition (DMD) for the optimal control of the convective FitzHugh-Nagumo (FHN) equations. The convective FHN equations consists of the semi-linear activator and the linear inhibitor equations, modeling blood coagulation in moving excitable media. The semilinear activator equation leads to a non-convex optimal control problem (OCP). The most commonly used method in reduced optimal control is POD. We use DEIM and DMD to approximate efficiently the nonlinear terms in reduced order models. We compare the accuracy and computational times of three reduced-order optimal control solutions with the full order discontinuous Galerkin finite element solution of the convection dominated FHN equations with terminal controls. Numerical results show that POD is the most accurate whereas POD-DMD is the fastest.

preprint2016arXiv

Moving Mesh Discontinuous Galerkin Methods for PDEs with Traveling Waves

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.

preprint2015arXiv

Energy Stable Discontinuous Galerkin Finite Element Method for the Allen-Cahn Equation

Allen--Cahn equation with constant and degenerate mobility, and with polynomial and logarithmic energy functionals is discretized using symmetric interior penalty discontinuous Galerkin (SIPG) finite elements in space. We show that the energy stable average vector field (AVF) method as the time integrator for gradient systems like the Allen-Cahn equation satisfies the energy decreasing property for the fully discrete scheme. The numerical results for one and two dimensional Allen-Cahn equation with periodic boundary condition, using adaptive time stepping, reveal that the discrete energy decreases monotonically, the phase separation and metastability phenomena can be observed and the ripening time is detected correctly.

preprint2015arXiv

Time-space adaptive discontinuous Galerkin method for advection-diffusion equations with non-linear reaction mechanism

In this work, we apply a time-space adaptive discontinuous Galerkin method using the elliptic reconstruction technique with a robust (in Péclet number) elliptic error estimator in space, for the convection dominated parabolic problems with non-linear reaction mechanisms. We derive a posteriori error estimators in the $L^{\infty}(L^2)+L^2(H^1)$-type norm using backward Euler in time and discontinuous Galerkin (symmetric interior penalty Galerkin (SIPG)) in space. Numerical results for advection dominated reactive transport problems in homogeneous and heterogeneous media demonstrate the performance of the time-space adaptive algorithm.

preprint2014arXiv

Adaptive discontinuous Galerkin methods for non-linear diffusion-convection-reaction equations

In this work, we apply the adaptive discontinuous Galerkin (DGAFEM) method to the convection dominated non-linear, quasi-stationary diffusion-convection-reaction equations. We propose an efficient preconditioner using a matrix reordering scheme to solve the sparse linear systems iteratively arising from the discretized non-linear equations. Numerical examples demonstrate effectiveness of the DGAFEM to damp the spurious oscillations and resolve well the sharp layers occurring in convection dominated non-linear equations.

preprint2014arXiv

Model Order Reduction for Nonlinear Schrödinger Equation

We apply the proper orthogonal decomposition (POD) to the nonlinear Schrödinger (NLS) equation to derive a reduced order model. The NLS equation is discretized in space by finite differences and is solved in time by structure preserving symplectic mid-point rule. A priori error estimates are derived for the POD reduced dynamical system. Numerical results for one and two dimensional NLS equations, coupled NLS equation with soliton solutions show that the low-dimensional approximations obtained by POD reproduce very well the characteristic dynamics of the system, such as preservation of energy and the solutions.