Researcher profile

Alexander Zlotnik

Alexander Zlotnik contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
12works
0followers
5topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

12 published item(s)

preprint2024arXiv

On properties of a semi-explicit in time fourth-order vector compact scheme for the multidimensional acoustic wave equation

We deal with an initial-boundary value problem for the multidimensional acoustic wave equation, with the variable speed of sound. For a three-level semi-explicit in time higher-order vector compact scheme, we prove stability and derive 4th order error bound in the enlarged energy norm. This scheme is three-point in each spatial direction, and it exploits additional sought functions which approximate 2nd order non-mixed spatial derivatives of the solution to the equation. At the first time level, a similar two-level in time scheme is applied, with no derivatives of the data. No iterations are required to implement the scheme. We also present results of various 3D numerical experiments that demonstrate a very high accuracy of the scheme for smooth data, its advantages in the error behavior over the classical explicit 2nd order scheme for nonsmooth data as well and an example of the wave in a layered medium initiated by the Ricker-type wavelet source function.

preprint2021arXiv

On Properties of Compact 4th order Finite-Difference Schemes for the Variable Coefficient Wave Equation

We consider an initial-boundary value problem for the $n$-dimensional wave equation with the variable sound speed, $n\geq 1$. We construct three-level implicit in time and compact in space (three-point in each space direction) 4th order finite-difference schemes on the uniform rectangular meshes including their one-parameter (for $n=2$) and three-parameter (for $n=3$) families. We also show that some already known methods can be converted into such schemes. In a unified manner, we prove the conditional stability of schemes in the strong and weak energy norms together with the 4th order error estimate under natural conditions on the time step. We also transform an unconditionally stable 4th order two-level scheme suggested for $n=2$ to the three-level form, extend it for any $n\geq 1$ and prove its stability. We also give an example of a compact scheme for non-uniform in space and time rectangular meshes. We suggest simple fast iterative methods based on FFT to implement the schemes. A new effective initial guess to start iterations is given too. We also present promising results of numerical experiments.

preprint2018arXiv

Finite element error analysis for measure-valued optimal control problems governed by a 1D wave equation with variable coefficients

This work is concerned with the optimal control problems governed by a 1D wave equation with variable coefficients and the control spaces $\mathcal M_T$ of either measure-valued functions $L_{w^*}^2(I,\mathcal M(Ω))$ or vector measures $\mathcal M(Ω,L^2(I))$. The cost functional involves the standard quadratic tracking terms and the regularization term $α\|u\|_{\mathcal M_T}$ with $α>0$. We construct and study three-level in time bilinear finite element discretizations for this class of problems. The main focus lies on the derivation of error estimates for the optimal state variable and the error measured in the cost functional. The analysis is mainly based on some previous results of the authors. The numerical results are included.

preprint2017arXiv

A "converse" stability condition is necessary for a compact higher order scheme on non-uniform meshes

The stability bounds and error estimates for a compact higher order Numerov-Crank-Nicolson scheme on non-uniform space meshes for the 1D time-dependent Schrödinger equation have been recently derived. This analysis has been done in $L^2$ and $H^1$ mesh norms and used the non-standard "converse" condition $h_ω\leq c_0τ$, where $h_ω$ is the mean space step, $τ$ is the time step and $c_0>0$. Now we prove that such condition is necessary for some families of non-uniform meshes and any space norm. Also numerical results show unacceptably wrong behavior of numerical solutions (their dramatic mass non-conservation) when this condition is violated.

preprint2017arXiv

On fast Fourier solvers for the tensor product high-order FEM for a generalized Poisson equation

We present direct logarithmically optimal in theory and fast in practice algorithms to implement the tensor product high order finite element method on multi-dimensional rectangular parallelepipeds for solving PDEs of the Poisson kind. They are based on the well-known Fourier approaches. The key new points are the fast direct and inverse FFT-based algorithms for expansion in eigenvectors of the 1D eigenvalue problems for the high order FEM. The algorithms can further be used for numerous applications, in particular, to implement the tensor product high order finite element methods for various time-dependent PDEs. Results of numerical experiments in 2D and 3D cases are presented.

preprint2016arXiv

A direct fast FFT-based implementation for high order finite element method on rectangular parallelepipeds for PDE

We present a new direct logarithmically optimal in theory and fast in practice algorithm to implement the high order finite element method on multi-dimensional rectangular parallelepipeds for solving PDEs of the Poisson kind. The key points are the fast direct and inverse FFT-based algorithms for decomposition in eigenvectors of the 1D eigenvalue problems for the high order FEM. The algorithm can further be used for numerous applications, in particular, to implement the high order finite element methods for various time-dependent PDEs.

preprint2015arXiv

On Error Estimates of the Crank-Nicolson-Polylinear Finite Element Method with the Discrete TBC for the Generalized Schrödinger Equation in an Unbounded Parallelepiped

We deal with an initial-boundary value problem for the generalized time-dependent Schrödinger equation with variable coefficients in an unbounded $n$--dimensional parallelepiped ($n\geq 1$). To solve it, the Crank-Nicolson in time and the polylinear finite element in space method with the discrete transpa\-rent boundary conditions is considered. We present its stability properties and derive new error estimates $O(τ^2+|h|^2)$ uniformly in time in $L^2$ space norm, for $n\geq 1$, and mesh $H^1$ space norm, for $1\leq n\leq 3$ (a superconvergence result), under the Sobolev-type assumptions on the initial function. Such estimates are proved for methods with the discrete TBCs for the first time.

preprint2014arXiv

Some remarks on discrete and semi-discrete transparent boundary conditions for solving the time-dependent Schrödinger equation on the half-axis

We consider the generalized time-dependent Schrödinger equation on the half-axis and a broad family of finite-difference schemes with the discrete transparent boundary conditions (TBCs) to solve it. We first rewrite the discrete TBCs in a simplified form explicit in space step $h$. Next, for a selected scheme of the family, we discover that the discrete convolution in time in the discrete TBC does not depend on $h$ and, moreover, it coincides with the corresponding convolution in the semi-discrete TBC rewritten similarly. This allows us to prove the bound for the difference between the kernels of the discrete convolutions in the discrete and semi-discrete TBCs (for the first time). Numerical experiments on replacing the discrete TBC convolutions by the semi-discrete one exhibit truly small absolute errors though not relative ones in general. The suitable discretization in space of the semi-discrete TBC for the higher-order Numerov scheme is also discussed.

preprint2013arXiv

A Numerov-Crank-Nicolson-Strang scheme with discrete transparent boundary conditions for the Schrödinger equation on a semi-infinite strip

We consider an initial-boundary value problem for a 2D time-dependent Schrödinger equation on a semi-infinite strip. For the Numerov-Crank-Nicolson finite-difference scheme with discrete transparent boundary conditions, the Strang-type splitting with respect to the potential is applied. For the resulting method, the uniqueness of a solution and the uniform in time $L^2$-stability (in particular, $L^2$-conservativeness) are proved. Due to the splitting, an effective direct algorithm using FFT in the direction perpendicular to the strip is developed to implement the splitting method for general potential. Numerical results on the tunnel effect for smooth and rectangular barriers together with the practical error analysis on refining meshes are included as well.

preprint2013arXiv

A splitting higher order scheme with discrete transparent boundary conditions for the Schrödinger equation in a semi-infinite parallelepiped

An initial-boundary value problem for the $n$-dimensional ($n\geq 2$) time-dependent Schrödinger equation in a semi-infinite (or infinite) parallelepiped is considered. Starting from the Numerov-Crank-Nicolson finite-difference scheme, we first construct higher order scheme with splitting space averages having much better spectral properties for $n\geq 3$. Next we apply the Strang-type splitting with respect to the potential and, third, construct discrete transparent boundary conditions (TBC). For the resulting method, the uniqueness of solution and the unconditional uniform in time $L^2$-stability (in particular, $L^2$-conservativeness) are proved. Owing to the splitting, an effective direct algorithm using FFT (in the coordinate directions perpendicular to the leading axis of the parallelepiped) is applicable for general potential. Numerical results on the 2D tunnel effect for a Pöschl-Teller-like potential-barrier and a rectangular potential-well are also included.

preprint2013arXiv

The splitting in potential Crank-Nicolson scheme with discrete transparent boundary conditions for the Schrödinger equation on a semi-infinite strip

We consider an initial-boundary value problem for a generalized 2D time-dependent Schrodinger equation (with variable coefficients) on a semi-infinite strip. For the Crank-Nicolson-type finite-difference scheme with approximate or discrete transparent boundary conditions (TBCs), the Strang-type splitting with respect to the potential is applied. For the resulting method, the unconditional uniform in time $L^2$-stability is proved. Due to the splitting, an effective direct algorithm using FFT is developed now to implement the method with the discrete TBC for general potential. Numerical results on the tunnel effect for rectangular barriers are included together with the detailed practical error analysis confirming nice properties of the method.

preprint2012arXiv

On a family of finite-difference schemes with discrete transparent boundary conditions for a parabolic equation on the half-axis

An initial-boundary value problem for the 1D self-adjoint parabolic equation on the half-axis is solved. We study a broad family of two-level finite-difference schemes with two parameters related to averagings both in time and space. Stability in two norms is proved by the energy method. Also discrete transparent boundary conditions are rigorously derived for schemes by applying the method of reproducing functions. Results of numerical experiments are included as well.