Source author record

Ludmil T. Zikatanov

Ludmil T. Zikatanov 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)

preprint2020arXiv

Discrete Trace Theorems and Energy Minimizing Spring Embeddings of Planar Graphs

Tutte's spring embedding theorem states that, for a three-connected planar graph, if the outer face of the graph is fixed as the complement of some convex region in the plane, and all other vertices are placed at the mass center of their neighbors, then this results in a unique embedding, and this embedding is planar. It also follows fairly quickly that this embedding minimizes the sum of squared edge lengths, conditional on the embedding of the outer face. However, it is not at all clear how to embed this outer face. We consider the minimization problem of embedding this outer face, up to some normalization, so that the sum of squared edge lengths is minimized. In this work, we show the connection between this optimization problem and the Schur complement of the graph Laplacian with respect to the interior vertices. We prove a number of discrete trace theorems, and, using these new results, show the spectral equivalence of this Schur complement with the boundary Laplacian to the one-half power for a large class of graphs. Using this result, we give theoretical guarantees for this optimization problem, which motivates an algorithm to embed the outer face of a spring embedding.

preprint2020arXiv

On the unisolvence for the quasi-polynomial spaces of differential forms

We consider quasi-polynomial spaces of differential forms defined as weighted (with a positive weight) spaces of differential forms with polynomial coefficients. We show that the unisolvent set of functionals for such spaces on a simplex in any spatial dimension is the same as the set of such functionals used for the polynomial spaces. The analysis in the quasi-polynomial spaces, however, is not standard and requires a novel approach. We are able to prove our results without the use of Stokes' Theorem, which is the standard tool in showing the unisolvence of functionals in polynomial spaces of differential forms. These new results provide tools for studying exponentially-fitted discretizations stable for general convection-diffusion problems in Hilbert differential complexes.

preprint2020arXiv

Robust preconditioners for a new stabilized discretization of the poroelastic equations

In this paper, we present block preconditioners for a stabilized discretization of the poroelastic equations developed in [45]. The discretization is proved to be well-posed with respect to the physical and discretization parameters, and thus provides a framework to develop preconditioners that are robust with respect to such parameters as well. We construct both norm-equivalent (diagonal) and field-of-value-equivalent (triangular) preconditioners for both the stabilized discretization and a perturbation of the stabilized discretization that leads to a smaller overall problem after static condensation. Numerical tests for both two- and three-dimensional problems confirm the robustness of the block preconditioners with respect to the physical and discretization parameters.

preprint2016arXiv

A Nonconforming Finite Element Method for the Biot's Consolidation Model in Poroelasticity

A stable finite element scheme that avoids pressure oscillations for a three-field Biot's model in poroelasticity is considered. The involved variables are the displacements, fluid flux (Darcy velocity), and the pore pressure, and they are discretized by using the lowest possible approximation order: Crouzeix-Raviart finite elements for the displacements, lowest order Raviart-Thomas-Nedelec elements for the Darcy velocity, and piecewise constant approximation for the pressure. Mass lumping technique is introduced for the Raviart-Thomas-Nedelec elements in order to eliminate the Darcy velocity and, therefore, reduce the computational cost. We show convergence of the discrete scheme which is implicit in time and use these types of elements in space with and without mass lumping. Finally, numerical experiments illustrate the convergence of the method and show its effectiveness to avoid spurious pressure oscillations when mass lumping for the Raviart-Thomas-Nedelec elements is used.

preprint2016arXiv

Arbitrary Dimension Convection-Diffusion Schemes for Space-Time Discretizations

This note proposes embedding a time dependent PDE into a convection-diffusion type PDE (in one space dimension higher) with singularity, for which two discretization schemes, the classical streamline-diffusion and the EAFE (edge average finite element) one, are investigated in terms of stability and error analysis. The EAFE scheme, in particular, is extended to be arbitrary order which is of interest on its own. Numerical results, in combined space-time domain demonstrate the feasibility of the proposed approach.

preprint2016arXiv

On the Approximation of Laplacian Eigenvalues in Graph Disaggregation

Graph disaggregation is a technique used to address the high cost of computation for power law graphs on parallel processors. The few high-degree vertices are broken into multiple small-degree vertices, in order to allow for more efficient computation in parallel. In particular, we consider computations involving the graph Laplacian, which has significant applications, including diffusion mapping and graph partitioning, among others. We prove results regarding the spectral approximation of the Laplacian of the original graph by the Laplacian of the disaggregated graph. In addition, we construct an alternate disaggregation operator whose eigenvalues interlace those of the original Laplacian. Using this alternate operator, we construct a uniform preconditioner for the original graph Laplacian.

preprint2016arXiv

Robust Solvers for Maxwell's Equations with Dissipative Boundary Conditions

In this paper, we design robust and efficient linear solvers for the numerical approximation of solutions to Maxwell's equations with dissipative boundary conditions. We consider a structure-preserving finite-element approximation with standard Nedelec--Raviart--Thomas elements in space and a Crank--Nicolson scheme in time to approximate the electric and magnetic fields. We focus on two types of block preconditioners. The first type is based on the well-posedness results of the discrete problem. The second uses an exact block factorization of the linear system, for which the structure-preserving discretization yields sparse Schur complements. We prove robustness and optimality of these block preconditioners, and provide supporting numerical tests.

preprint2015arXiv

On the Maximal Error of Spectral Approximation of Graph Bisection

Spectral graph bisections are a popular heuristic aimed at approximating the solution of the NP-complete graph bisection problem. This technique, however, does not always provide a robust tool for graph partitioning. Using a special class of graphs, we prove that the standard spectral graph bisection can produce bisections that are far from optimal. In particular, we show that the maximum error in the spectral approximation of the optimal bisection (partition sizes exactly equal) cut for such graphs is bounded below by a constant multiple of the order of the graph squared.

preprint2014arXiv

A Cascadic Multigrid Algorithm for Computing the Fiedler Vector of Graph Laplacians

In this paper, we develop a cascadic multigrid algorithm for fast computation of the Fiedler vector of a graph Laplacian, namely, the eigenvector corresponding to the second smallest eigenvalue. This vector has been found to have applications in fields such as graph partitioning and graph drawing. The algorithm is a purely algebraic approach based on a heavy edge coarsening scheme and pointwise smoothing for refinement. To gain theoretical insight, we also consider the related cascadic multigrid method in the geometric setting for elliptic eigenvalue problems and show its uniform convergence under certain assumptions. Numerical tests are presented for computing the Fiedler vector of several practical graphs, and numerical results show the efficiency and optimality of our proposed cascadic multigrid algorithm.

preprint2014arXiv

A uniform additive Schwarz preconditioner for the $hp$-version of Discontinuous Galerkin approximations of elliptic problems

In this paper we design and analyze a uniform preconditioner for a class of high order Discontinuous Galerkin schemes. The preconditioner is based on a space splitting involving the high order conforming subspace and results from the interpretation of the problem as a nearly-singular problem. We show that the proposed preconditioner exhibits spectral bounds that are uniform with respect to the discretization parameters, i.e., the mesh size, the polynomial degree and the penalization coefficient. The theoretical estimates obtained are supported by several numerical simulations.

preprint2013arXiv

Numerical Approximation of Asymptotically Disappearing Solutions of Maxwell's Equations

This work is on the numerical approximation of incoming solutions to Maxwell's equations with dissipative boundary conditions whose energy decays exponentially with time. Such solutions are called asymptotically disappearing (ADS) and they play an importarnt role in inverse back-scatering problems. The existence of ADS is a difficult mathematical problem. For the exterior of a sphere, such solutions have been constructed analytically by Colombini, Petkov and Rauch [7] by specifying appropriate initial conditions. However, for general domains of practical interest (such as Lipschitz polyhedra), the existence of such solutions is not evident. This paper considers a finite-element approximation of Maxwell's equations in the exterior of a polyhedron, whose boundary approximates the sphere. Standard Nedelec-Raviart-Thomas elements are used with a Crank-Nicholson scheme to approximate the electric and magnetic fields. Discrete initial conditions interpolating the ones chosen in [7] are modified so that they are (weakly) divergence-free. We prove that with such initial conditions, the approximation to the electric field is weakly divergence-free for all time. Finally, we show numerically that the finite-element approximations of the ADS also decay exponentially with time when the mesh size and the time step become small.

preprint2012arXiv

An exponential fitting scheme for general convection-diffusion equations on tetrahedral meshes

This paper contains construction and analysis a finite element approximation for convection dominated diffusion problems with full coefficient matrix on general simplicial partitions in $R^d$, $d=2,3$. This construction is quite close to the scheme of Xu and Zikatanov (Math. Comp. 1999) where a diagonal coefficient matrix has been considered. The scheme is of the class of exponentially fitted methods that does not use upwind or checking the flow direction. It is stable for sufficiently small discretization step-size assuming that the boundary value problem for the convection-diffusion equation is uniquely solvable. Further, it is shown that, under certain conditions on the mesh the scheme is monotone. Convergence of first order is derived under minimal smoothness of the solution.

preprint2012arXiv

Polynomial of best uniform approximation to $x^{-1}$ and smoothing in two-level methods

We derive a three-term recurrence relation for computing the polynomial of best approximation in the uniform norm to $x^{-1}$ on a finite interval with positive endpoints. As application, we consider two-level methods for scalar elliptic partial differential equation (PDE), where the relaxation on the fine grid uses the aforementioned polynomial of best approximation. Based on a new smoothing property of this polynomial smoother that we prove, combined with a proper choice of the coarse space, we obtain as a corollary, that the convergence rate of the resulting two-level method is uniform with respect to the mesh parameters, coarsening ratio and PDE coefficient variation.