-
Linearly implicit exponential integrators for damped Hamiltonian PDEs
Authors:
Murat Uzunca,
Bülent Karasözen
Abstract:
Structure-preserving linearly implicit exponential integrators are constructed for Hamiltonian partial differential equations with linear constant damping. Linearly implicit integrators are derived by polarizing the polynomial terms of the Hamiltonian function and portioning out the nonlinearly of consecutive time steps. They require only a solution of one linear system at each time step. Therefor…
▽ More
Structure-preserving linearly implicit exponential integrators are constructed for Hamiltonian partial differential equations with linear constant damping. Linearly implicit integrators are derived by polarizing the polynomial terms of the Hamiltonian function and portioning out the nonlinearly of consecutive time steps. They require only a solution of one linear system at each time step. Therefore they are computationally more advantageous than implicit integrators. We also construct an exponential version of the well-known one-step Kahan's method by polarizing the quadratic vector field. These integrators are applied to one-dimensional damped Burger's, Korteweg-de-Vries, and nonlinear Schr{ö}dinger equations. Preservation of the dissipation rate of linear and quadratic conformal invariants and the Hamiltonian is illustrated by numerical experiments.
△ Less
Submitted 18 March, 2024; v1 submitted 25 September, 2023;
originally announced September 2023.
-
Reduced-order modeling for Ablowitz-Ladik equation
Authors:
Murat Uzunca,
Bülent Karasözen
Abstract:
In this paper, reduced-order models (ROMs) are constructed for the Ablowitz-Ladik equation (ALE), an integrable semi-discretization of the nonlinear Schrödinger equation (NLSE) with and without damping. Both ALEs are non-canonical conservative and dissipative Hamiltonian systems with the Poisson matrix, depending quadratically on the state variables and with quadratic Hamiltonian. The full-order s…
▽ More
In this paper, reduced-order models (ROMs) are constructed for the Ablowitz-Ladik equation (ALE), an integrable semi-discretization of the nonlinear Schrödinger equation (NLSE) with and without damping. Both ALEs are non-canonical conservative and dissipative Hamiltonian systems with the Poisson matrix, depending quadratically on the state variables and with quadratic Hamiltonian. The full-order solutions are obtained with the energy preserving midpoint rule for the conservative ALE and exponential midpoint rule for the dissipative ALE. The reduced-order solutions are constructed intrusively by preserving the skew-symmetric structure of the reduced non-canonical Hamiltonian system by applying proper orthogonal decomposition with the Galerkin projection. For an efficient offline-online decomposition of the ROMs, the quadratic nonlinear terms of the Poisson matrix are approximated by the discrete empirical interpolation method. The computation of the reduced-order solutions is further accelerated by the use of tensor techniques. Preservation of the Hamiltonian and momentum for the conservative ALE, and preservation of dissipation properties of the dissipative ALE, guarantee the long-term stability of soliton solutions.
△ Less
Submitted 23 June, 2022;
originally announced July 2022.
-
Global energy preserving model reduction for multi-symplectic PDEs
Authors:
Murat Uzunca,
Bülent Karasözen,
Ayhan Aydın
Abstract:
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 a…
▽ More
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.
△ Less
Submitted 29 August, 2022; v1 submitted 13 February, 2022;
originally announced March 2022.
-
Nonintrusive model order reduction for cross-diffusion systems
Authors:
Bulent Karasozen,
Murat Uzunca,
Gulden Mulayim
Abstract:
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…
▽ More
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.
△ Less
Submitted 5 September, 2021;
originally announced September 2021.
-
Intrusive and non-intrusive reduced order modeling of the rotating thermal shallow water equation
Authors:
Süleyman Yıldız,
Murat Uzunca,
Bülent Karasözen
Abstract:
In this paper, we investigate projection-based intrusive and data-driven non-intrusive model order reduction methods in numerical simulation of rotating thermal shallow water equation (RTSWE) in parametric and non-parametric form. Discretization of the RTSWE in space with centered finite differences leads to Hamiltonian system of ordinary differential equations with linear and quadratic terms. The…
▽ More
In this paper, we investigate projection-based intrusive and data-driven non-intrusive model order reduction methods in numerical simulation of rotating thermal shallow water equation (RTSWE) in parametric and non-parametric form. Discretization of the RTSWE in space with centered finite differences leads to Hamiltonian system of ordinary differential equations with linear and quadratic terms. The full-order model (FOM) is obtained by applying linearly implicit Kahan's method in time. Applying proper orthogonal decomposition with Galerkin projection (POD-G), we construct the intrusive reduced-order model (ROM). We apply operator inference (OpInf) with re-projection for non-intrusive reduced-order modeling. In the parametric case, we make use of the parameter dependency at the level of the PDE without interpolating between the reduced operators. The least-squares problem of the OpInf is regularized with the minimum norm solution. Both ROMs behave similar and are able to accurately predict the test and training data and capture system behavior in the prediction phase with several orders of computational speedup over the FOM. The preservation of system physics such as the conserved quantities of the RTSWE by both ROMs enables that the models fit better to data and stable solutions are obtained in long-term predictions, which are robust to parameter changes.
△ Less
Submitted 14 April, 2021; v1 submitted 31 March, 2021;
originally announced April 2021.
-
Energy preserving reduced-order modelling of thermal shallow water equation
Authors:
Suleyman Yildiz,
Murat Uzunca,
Bulent Karasozen
Abstract:
In this paper, Hamiltonian and energy preserving reduced-order models are developed for the rotating thermal shallow water equation (RTSWE) in the non-canonical Hamiltonian form with the state-dependent Poisson matrix. The high fidelity full solutions are obtained by discretizing the RTSWE in space with skew-symmetric finite-differences, that preserve the Hamiltonian structure. The resulting skew-…
▽ More
In this paper, Hamiltonian and energy preserving reduced-order models are developed for the rotating thermal shallow water equation (RTSWE) in the non-canonical Hamiltonian form with the state-dependent Poisson matrix. The high fidelity full solutions are obtained by discretizing the RTSWE in space with skew-symmetric finite-differences, that preserve the Hamiltonian structure. The resulting skew-gradient system is integrated in time with the energy preserving average vector field (AVF) method. The reduced-order model (ROM) is constructed in the same way as the full order model (FOM), preserving the reduced skew-symmetric structure and integrating in time with the AVF method. Relying on structure-preserving discretizations in space and time and applying proper orthogonal decomposition (POD) with the Galerkin projection, an energy preserving reduced order model (ROM) is constructed. The nonlinearities in the ROM are computed by applying the discrete empirical interpolation (DEIM) method to reduce the computational cost. The computation of the reduced-order solutions is accelerated further by the use of tensor techniques. The overall procedure yields a clear separation of the offline and online computational cost of the reduced solutions. The accuracy and computational efficiency of the ROMs are demonstrated for a numerical test problem. Preservation of the energy (Hamiltonian), and other conserved quantities, i.e. mass, buoyancy, and total vorticity show that the reduced-order solutions ensure the long-term stability of the solutions while exhibiting several orders of magnitude computational speedup over the FOM.
△ Less
Submitted 16 February, 2021; v1 submitted 1 November, 2020;
originally announced November 2020.
-
Data-Driven Learning of Reduced-order Dynamics for a Parametrized Shallow Water Equation
Authors:
Süleyman Yıldız,
Pawan Goyal,
Peter Benner,
Bülent Karasözen
Abstract:
This paper discusses a non-intrusive data-driven model order reduction method that learns low-dimensional dynamical models for a parametrized shallow water equation. We consider the shallow water equation in non-traditional form (NTSWE). We focus on learning low-dimensional models in a non-intrusive way. That means, we assume not to have access to a discretized form of the NTSWE in any form. Inste…
▽ More
This paper discusses a non-intrusive data-driven model order reduction method that learns low-dimensional dynamical models for a parametrized shallow water equation. We consider the shallow water equation in non-traditional form (NTSWE). We focus on learning low-dimensional models in a non-intrusive way. That means, we assume not to have access to a discretized form of the NTSWE in any form. Instead, we have snapshots that are obtained using a black-box solver. Consequently, we aim at learning reduced-order models only from the snapshots. Precisely, a reduced-order model is learnt by solving an appropriate least-squares optimization problem in a low-dimensional subspace. Furthermore, we discuss computational challenges that particularly arise from the optimization problem being ill-conditioned. Moreover, we extend the non-intrusive model order reduction framework to a parametric case where we make use of the parameter dependency at the level of the partial differential equation. We illustrate the efficiency of the proposed non-intrusive method to construct reduced-order models for NTSWE and compare it with an intrusive method (proper orthogonal decomposition). We furthermore discuss the predictive capabilities of both models outside the range of the training data.
△ Less
Submitted 4 August, 2020; v1 submitted 28 July, 2020;
originally announced July 2020.
-
Reduced order modelling of nonlinear cross-diffusion systems
Authors:
Bülent Karasözen,
Gülden Mülayim,
Murat Uzunca,
Süleyman Yıldız
Abstract:
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…
▽ More
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.
△ Less
Submitted 8 February, 2021; v1 submitted 9 June, 2020;
originally announced June 2020.
-
Structure-preserving reduced-order modelling of Korteweg de Vries equation
Authors:
Bulent Karasozen,
Murat Uzunca,
Suleyman Yildiz
Abstract:
Computationally efficient, structure-preserving reduced-order methods are developed for the Korteweg de Vries (KdV) equations in Hamiltonian form. The KdV equation is discretized in space by finite differences. 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…
▽ More
Computationally efficient, structure-preserving reduced-order methods are developed for the Korteweg de Vries (KdV) equations in Hamiltonian form. The KdV equation is discretized in space by finite differences. 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 the full-order model (FOM) is preserved by the reduced-order model (ROM). The quadratic nonlinear terms of the KdV equation are evaluated efficiently by the use of tensorial methods, clearly separating the offline-online cost of the FOMs and ROMs. The accuracy of the reduced solutions, preservation of the Hamiltonian, momentum and mass, and computational speed-up gained by ROMs are demonstrated for the one-dimensional KdV equation, coupled KdV equations and two-dimensional Zakharov-Kuznetzov equation with soliton solutions
△ Less
Submitted 7 February, 2021; v1 submitted 17 April, 2020;
originally announced April 2020.
-
Model Order Reduction in Neuroscience
Authors:
Bülent Karasözen
Abstract:
The human brain contains approximately $10^9$ neurons, each with approximately $10^3$ connections, synapses, with other neurons. Most sensory, cognitive and motor functions of our brains depend on the interaction of a large population of neurons. In recent years, many technologies are developed for recording large numbers of neurons either sequentially or simultaneously. An increase in computation…
▽ More
The human brain contains approximately $10^9$ neurons, each with approximately $10^3$ connections, synapses, with other neurons. Most sensory, cognitive and motor functions of our brains depend on the interaction of a large population of neurons. In recent years, many technologies are developed for recording large numbers of neurons either sequentially or simultaneously. An increase in computational power and algorithmic developments have enabled advanced analyses of neuronal population parallel to the rapid growth of quantity and complexity of the recorded neuronal activity. Recent studies made use of dimensionality and model order reduction techniques to extract coherent features which are not apparent at the level of individual neurons. It has been observed that the neuronal activity evolves on low-dimensional subspaces. The aim of model reduction of large-scale neuronal networks is an accurate and fast prediction of patterns and their propagation in different areas of the brain. Spatiotemporal features of the brain activity are identified on low dimensional subspaces with methods such as dynamic mode decomposition (DMD), proper orthogonal decomposition (POD), discrete empirical interpolation (DEIM) and combined parameter and state reduction. In this paper, we give an overview of the currently used dimensionality reduction and model order reduction techniques in neuroscience.
This work will be featured as a chapter in the upcoming Handbook on Model Order Reduction,(P. Benner, S. Grivet-Talocia, A. Quarteroni, G. Rozza, W. H. A. Schilders, L. M. Silveira, eds, to appear on DE GRUYTER)
△ Less
Submitted 17 April, 2020; v1 submitted 11 March, 2020;
originally announced March 2020.
-
Structure-preserving Reduced Order Modeling of non-traditional Shallow Water Equation
Authors:
Süleyman Yıldız,
Murat Uzunca,
Bülent Karasözen
Abstract:
An energy preserving reduced order model is developed for the nontraditional shallow water equation (NTSWE) with full Coriolis force. The NTSWE in the noncanonical Hamiltonian/Poisson form is discretized in space by finite differences. The resulting system of ordinary differential equations is integrated in time by the energy preserving average vector field (AVF) method. The Poisson structure of t…
▽ More
An energy preserving reduced order model is developed for the nontraditional shallow water equation (NTSWE) with full Coriolis force. The NTSWE in the noncanonical Hamiltonian/Poisson form is discretized in space by finite differences. The resulting system of ordinary differential equations is integrated in time by the energy preserving average vector field (AVF) method. The Poisson structure of the NTSWE in discretized exhibits a skew-symmetric matrix depending on the state variables. An energy preserving, computationally efficient reduced-order model (ROM) is constructed by proper orthogonal decomposition with Galerkin projection. The nonlinearities are computed for the ROM efficiently by discrete empirical interpolation method. Preservation of the semi-discrete energy and the enstrophy are shown for the full order model, and for the ROM which ensures the long term stability of the solutions. The accuracy and computational efficiency of the ROMs are shown by two numerical test problems
△ Less
Submitted 26 February, 2020;
originally announced February 2020.
-
Structure Preserving Model Order Reduction of Shallow Water Equations
Authors:
Bülent Karasözen,
Süleyman Yıldız,
Murat Uzunca
Abstract:
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 me…
▽ More
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.
△ Less
Submitted 22 July, 2020; v1 submitted 10 July, 2019;
originally announced July 2019.
-
Structure preserving reduced order modeling for gradient systems
Authors:
Tuğba Akman Yıldız,
Murat Uzunca,
Bülent Karasözen
Abstract:
Minimization of energy in gradient systems leads to formation of oscillatory and Turing patterns in reaction-diffusion systems. These patterns should be accurately computed using fine space and time meshes over long time horizons to reach the spatially inhomogeneous steady state. In this paper, a reduced order model (ROM) is developed which preserves the gradient dissipative structure. The coupled…
▽ More
Minimization of energy in gradient systems leads to formation of oscillatory and Turing patterns in reaction-diffusion systems. These patterns should be accurately computed using fine space and time meshes over long time horizons to reach the spatially inhomogeneous steady state. In this paper, a reduced order model (ROM) is developed which preserves the gradient dissipative structure. The coupled system of reaction-diffusion equations are discretized in space by the symmetric interior penalty discontinuous Galerkin (SIPG) method. The resulting system of ordinary differential equations (ODEs) are integrated in time by the average vector field (AVF) method, which preserves the energy dissipation of the gradient systems. The ROMs are constructed by the proper orthogonal decomposition (POD) with Galerkin projection. The nonlinear reaction terms are computed efficiently by discrete empirical interpolation method (DEIM). Preservation of the discrete energy of the FOMs and ROMs with POD-DEIM ensures the long term stability of the steady state solutions. Numerical simulations are performed for the gradient dissipative systems with two specific equations; real Ginzburg-Landau equation and Swift-Hohenberg equation. Numerical results demonstrate that the POD-DEIM reduced order solutions preserve well the energy dissipation over time and at the steady state.
△ Less
Submitted 27 November, 2018; v1 submitted 13 June, 2018;
originally announced June 2018.
-
Energy preserving model order reduction of the nonlinear Schrödinger equation
Authors:
Bülent Karasözen,
Murat Uzunca
Abstract:
An energy preserving reduced order model is developed for two dimensional nonlinear Schrödinger equation (NLSE) with plane wave solutions and with an external potential. The NLSE is discretized in space by the symmetric interior penalty discontinuous Galerkin (SIPG) method. The resulting system of Hamiltonian ordinary differential equations are integrated in time by the energy preserving average v…
▽ More
An energy preserving reduced order model is developed for two dimensional nonlinear Schrödinger equation (NLSE) with plane wave solutions and with an external potential. The NLSE is discretized in space by the symmetric interior penalty discontinuous Galerkin (SIPG) method. The resulting system of Hamiltonian ordinary differential equations are integrated in time by the energy preserving average vector field (AVF) method. The mass and energy preserving reduced order model (ROM) is constructed by proper orthogonal decomposition (POD) Galerkin projection. The nonlinearities are computed for the ROM efficiently by discrete empirical interpolation method (DEIM) and dynamic mode decomposition (DMD). Preservation of the semi-discrete energy and mass are shown for the full order model (FOM) and for the ROM which ensures the long term stability of the solutions. Numerical simulations illustrate the preservation of the energy and mass in the reduced order model for the two dimensional NLSE with and without the external potential. The POD-DMD makes a remarkable improvement in computational speed-up over the POD-DEIM. Both methods approximate accurately the FOM, whereas POD-DEIM is more accurate than the POD-DMD.
△ Less
Submitted 14 February, 2018; v1 submitted 30 May, 2017;
originally announced June 2017.
-
Reduced Order Optimal Control of the Convective FitzHugh-Nagumo Equation
Authors:
Bülent Karasözen,
Murat Uzunca,
Tuğba Küçükseyhan
Abstract:
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 mo…
▽ More
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.
△ Less
Submitted 20 August, 2019; v1 submitted 28 February, 2017;
originally announced March 2017.
-
Reduced-Order Modeling for Heston Stochastic Volatility Model
Authors:
Sinem Kozpınar,
Murat Uzunca,
Bülent Karasözen
Abstract:
In this paper, we compare the intrusive proper orthogonal decomposition (POD) with Galerkin projection and the data-driven dynamic mode decomposition (DMD), for Heston's option pricing model. The full order model is obtained by discontinuous Galerkin discretization in space and backward Euler in time. Numerical results for butterfly spread, European and digital call options reveal that in general…
▽ More
In this paper, we compare the intrusive proper orthogonal decomposition (POD) with Galerkin projection and the data-driven dynamic mode decomposition (DMD), for Heston's option pricing model. The full order model is obtained by discontinuous Galerkin discretization in space and backward Euler in time. Numerical results for butterfly spread, European and digital call options reveal that in general DMD requires more modes than the POD modes for the same level of accuracy. However, the speed-up factors are much higher for DMD than POD due to the non-intrusive nature of the DMD.
△ Less
Submitted 17 March, 2024; v1 submitted 29 October, 2016;
originally announced November 2016.
-
Energy stable model order reduction for the Allen-Cahn equation
Authors:
Murat Uzunca,
Bülent Karasözen
Abstract:
The Allen-Cahn equation is a gradient system, where the free-energy functional decreases monotonically in time. We develop an energy stable reduced order model (ROM) for a gradient system, which inherits the energy decreasing property of the full order model (FOM). For the space discretization we apply a discontinuous Galerkin (dG) method and for time discretization the energy stable average vecto…
▽ More
The Allen-Cahn equation is a gradient system, where the free-energy functional decreases monotonically in time. We develop an energy stable reduced order model (ROM) for a gradient system, which inherits the energy decreasing property of the full order model (FOM). For the space discretization we apply a discontinuous Galerkin (dG) method and for time discretization the energy stable average vector field (AVF) method. We construct ROMs with proper orthogonal decomposition (POD)-greedy adaptive sampling of the snapshots in time and evaluating the nonlinear function with greedy discrete empirical interpolation method (DEIM). The computational efficiency and accuracy of the reduced solutions are demonstrated numerically for the parametrized Allen-Cahn equation with Neumann and periodic boundary conditions.
△ Less
Submitted 5 February, 2017; v1 submitted 25 June, 2016;
originally announced June 2016.
-
Pricing European and American Options under Heston Model using Discontinuous Galerkin Finite Elements
Authors:
Sinem Kozpınar,
Murat Uzunca,
Bülent Karasözen
Abstract:
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 Amer…
▽ More
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.
△ Less
Submitted 15 March, 2020; v1 submitted 30 May, 2016;
originally announced June 2016.
-
Moving Mesh Discontinuous Galerkin Methods for PDEs with Traveling Waves
Authors:
Murat Uzunca,
Bülent Karasözen,
Tuğba Küçükseyhan
Abstract:
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…
▽ More
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.
△ Less
Submitted 31 January, 2016;
originally announced February 2016.
-
Model Order Reduction for Pattern Formation in FitzHugh-Nagumo Equation
Authors:
Bülent Karasözen,
Murat Uzunca,
Tuğba Küçükseyhan
Abstract:
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 (FO…
▽ More
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.
△ Less
Submitted 15 December, 2015;
originally announced December 2015.
-
Time-Space Adaptive Method of Time Layers for the Advective Allen-Cahn Equation
Authors:
Murat Uzunca,
Bülent Karasözen,
Ayşe Sarıaydın Filibelioğlu
Abstract:
We develop an adaptive method of time layers with a linearly implicit Rosenbrock method as time integrator and symmetric interior penalty Galerkin method for space discretization for the advective Allen-Cahn equation with non-divergence-free velocity fields. Numerical simulations for convection dominated problems demonstrate the accuracy and efficiency of the adaptive algorithm for resolving the s…
▽ More
We develop an adaptive method of time layers with a linearly implicit Rosenbrock method as time integrator and symmetric interior penalty Galerkin method for space discretization for the advective Allen-Cahn equation with non-divergence-free velocity fields. Numerical simulations for convection dominated problems demonstrate the accuracy and efficiency of the adaptive algorithm for resolving the sharp layers occurring in interface problems with small surface tension.
△ Less
Submitted 15 December, 2015;
originally announced December 2015.
-
Optimal Control of Convective FitzHugh-Nagumo Equation
Authors:
Murat Uzunca,
Tuğba Küçükseyhan,
Hamdullah Yücel,
Bülent Karasözen
Abstract:
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 nu…
▽ More
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.
△ Less
Submitted 28 February, 2017; v1 submitted 9 September, 2015;
originally announced September 2015.
-
Structure preserving integration and model order reduction of skew-gradient reaction-diffusion systems
Authors:
Bülent Karasözen,
Tuğba Küçükseyhan,
Murat Uzunca
Abstract:
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 dissip…
▽ More
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.
△ Less
Submitted 22 October, 2015; v1 submitted 31 July, 2015;
originally announced August 2015.
-
Energy stable interior penalty discontinuous Galerkin finite element method for Cahn-Hilliard equation
Authors:
Bülent Karasözen,
Ayşe Sarıaydın Filibelioğlu,
Murat Uzunca
Abstract:
An energy stable conservative method is developed for the Cahn--Hilliard (CH) equation with the degenerate mobility. The CH equation is discretized in space with the mass conserving symmetric interior penalty discontinuous Galerkin (SIPG) method. The resulting semi-discrete nonlinear system of ordinary differential equations are solved in time by the unconditionally energy stable average vector fi…
▽ More
An energy stable conservative method is developed for the Cahn--Hilliard (CH) equation with the degenerate mobility. The CH equation is discretized in space with the mass conserving symmetric interior penalty discontinuous Galerkin (SIPG) method. The resulting semi-discrete nonlinear system of ordinary differential equations are solved in time by the unconditionally energy stable average vector field (AVF) method. We prove that the AVF method preserves the energy decreasing property of the CH equation. Numerical results confirm the theoretical convergence rates and the performance of the proposed approach.
△ Less
Submitted 4 June, 2015;
originally announced June 2015.
-
Time-space adaptive discontinuous Galerkin method for advection-diffusion equations with non-linear reaction mechanism
Authors:
Bülent Karasözen,
Murat Uzunca
Abstract:
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 discontinuo…
▽ More
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.
△ Less
Submitted 17 May, 2015;
originally announced May 2015.
-
A Matlab Tutorial for Diffusion-Convection-Reaction Equations using DGFEM
Authors:
Murat Uzunca,
Bülent Karasözen
Abstract:
We present a collection of MATLAB routines using discontinuous Galerkin finite elements method (DGFEM) for solving steady-state diffusion-convection-reaction equations. The code employs the sparse matrix facilities of MATLAB with "vectorization" and uses multiple matrix multiplications {\it "MULTIPROD"} to increase the efficiency of the program.
We present a collection of MATLAB routines using discontinuous Galerkin finite elements method (DGFEM) for solving steady-state diffusion-convection-reaction equations. The code employs the sparse matrix facilities of MATLAB with "vectorization" and uses multiple matrix multiplications {\it "MULTIPROD"} to increase the efficiency of the program.
△ Less
Submitted 10 February, 2015;
originally announced February 2015.
-
Adaptive discontinuous Galerkin methods for non-linear diffusion-convection-reaction equations
Authors:
Murat Uzunca,
Bülent Karasözen,
Murat Manguoğlu
Abstract:
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 DGAFE…
▽ More
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.
△ Less
Submitted 15 September, 2014;
originally announced September 2014.
-
Energy Stable Discontinuous Galerkin Finite Element Method for the Allen-Cahn Equation
Authors:
Bülent Karasözen,
Ayşe Sarıaydın Filibelioğlu,
Murat Uzunca
Abstract:
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 fo…
▽ More
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.
△ Less
Submitted 16 May, 2015; v1 submitted 13 September, 2014;
originally announced September 2014.
-
Model Order Reduction for Nonlinear Schrödinger Equation
Authors:
Bülent Karasözen,
Canan Akkoyunlu,
Murat Uzunca
Abstract:
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,…
▽ More
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.
△ Less
Submitted 13 September, 2014;
originally announced September 2014.
-
Space-Time Discontinuous Galerkin Solution of Convection Dominated Optimal Control Problems
Authors:
Tuğba Akman,
Bülent Karasözen
Abstract:
In this paper, a space-time discontinuous Galerkin finite element method for distributed optimal control problems governed by unsteady diffusion-convection-reaction equations with control constraints is studied. Time discretization is performed by discontinuous Galerkin method with piecewise constant and linear polynomials, while symmetric interior penalty Galerkin with upwinding is used for space…
▽ More
In this paper, a space-time discontinuous Galerkin finite element method for distributed optimal control problems governed by unsteady diffusion-convection-reaction equations with control constraints is studied. Time discretization is performed by discontinuous Galerkin method with piecewise constant and linear polynomials, while symmetric interior penalty Galerkin with upwinding is used for space discretization. The numerical results presented confirm the theoretically observed convergence rates.
△ Less
Submitted 8 August, 2013;
originally announced August 2013.
-
Destruction of the family of steady states in the planar problem of Darcy convection
Authors:
V. G. Tsybulin,
B. Karasözen
Abstract:
The natural convection of incompressible fluid in a porous medium causes for some boundary conditions a strong non-uniqueness in the form of a continuous family of steady states. We are interested in the situation when these boundary conditions are violated. The resulting destruction of the family of steady states is studied via computer experiments based on a mimetic finite-difference approach.…
▽ More
The natural convection of incompressible fluid in a porous medium causes for some boundary conditions a strong non-uniqueness in the form of a continuous family of steady states. We are interested in the situation when these boundary conditions are violated. The resulting destruction of the family of steady states is studied via computer experiments based on a mimetic finite-difference approach. Convection in a rectangular enclosure is considered under different perturbations of boundary conditions (heat sources, infiltration). Two scenario of the family of equilibria are found: the transformation to a limit cycle and the formation of isolated convective patterns.
△ Less
Submitted 1 March, 2008;
originally announced March 2008.
-
Staggered grids discretization in three-dimensional Darcy convection
Authors:
B. Karasozen,
A. D. Nemtsev,
V. G. Tsybulin
Abstract:
We consider three-dimensional convection of an incompressible fluid saturated in a parallelepiped with a porous medium. A mimetic finite-difference scheme for the Darcy convection problem in the primitive variables is developed. It consists of staggered nonuniform grids with five types of nodes, differencing and averaging operators on a two-nodes stencil. The nonlinear terms are approximated usi…
▽ More
We consider three-dimensional convection of an incompressible fluid saturated in a parallelepiped with a porous medium. A mimetic finite-difference scheme for the Darcy convection problem in the primitive variables is developed. It consists of staggered nonuniform grids with five types of nodes, differencing and averaging operators on a two-nodes stencil. The nonlinear terms are approximated using special schemes. Two problems with different boundary conditions are considered to study scenarios of instability of the state of rest. Branching off of a continuous family of steady states was detected for the problem with zero heat fluxes on two opposite lateral planes.
△ Less
Submitted 17 February, 2008; v1 submitted 28 October, 2007;
originally announced October 2007.
-
Poisson integrators
Authors:
B. Karasözen
Abstract:
An overview of Hamiltonian systems with noncanonical Poisson structures is given. Examples of bi-Hamiltonian ode's, pde's and lattice equations are presented. Numerical integrators using generating functions, Hamiltonian splitting, symplectic Runge-Kutta methods are discussed for Lie-Poisson systems and Hamiltonian systems with a general Poisson structure. Nambu-Poisson systems and the discrete…
▽ More
An overview of Hamiltonian systems with noncanonical Poisson structures is given. Examples of bi-Hamiltonian ode's, pde's and lattice equations are presented. Numerical integrators using generating functions, Hamiltonian splitting, symplectic Runge-Kutta methods are discussed for Lie-Poisson systems and Hamiltonian systems with a general Poisson structure. Nambu-Poisson systems and the discrete gradient methods are also presented.
△ Less
Submitted 13 May, 2005;
originally announced May 2005.
-
Poisson integrators for Volterra lattice equations
Authors:
T. Ergenç,
B. Karasözen
Abstract:
The Volterra lattice equations are completely integrable and possess bi-Hamiltonian structure. They are integrated using partitioned Lobatto IIIA-B methods which preserve the Poisson structure. Modified equations are derived for the symplectic Euler and second order Lobatto IIIA-B method. Numerical results confirm preservation of the corresponding Hamiltonians, Casimirs, quadratic and cubic inte…
▽ More
The Volterra lattice equations are completely integrable and possess bi-Hamiltonian structure. They are integrated using partitioned Lobatto IIIA-B methods which preserve the Poisson structure. Modified equations are derived for the symplectic Euler and second order Lobatto IIIA-B method. Numerical results confirm preservation of the corresponding Hamiltonians, Casimirs, quadratic and cubic integrals in the long-term with different orders of accuracy.
△ Less
Submitted 13 May, 2005;
originally announced May 2005.