BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

N. Sukumar

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

ResearcherUnclaimed source record
math.NANumerical Analysiscond-mat.mtrl-sciphysics.comp-phComputational Engineering, Finance, and ScienceNeural and Evolutionary Computing

Catalog footprint

What is connected

9works
6topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

9 published item(s)

preprint2022arXiv

Virtual element method for modeling the deformation of multiphase composites

In this paper, we study applications of the virtual element method (VEM) for simulating the deformation of multiphase composites. The VEM is a Galerkin approach that is applicable to meshes that consist of arbitrarily-shaped polygonal and polyhedral (simple and nonsimple) elements. In the VEM, the basis functions are defined as the solution of a local elliptic partial differential equation, and are never explicitly computed in the implementation of the method. The stifness matrix of each element is built by using the elliptic projection operator of the internal virtual work (bilinear form) and it consists of two terms: a consistency term that is exactly computed (linear patch test is satisfied) and a correction term (ensures stability) that is orthogonal to affine displacement fields and has the right scaling. The VEM simplifies mesh generation for a multiphase composite: a stiff inclusion can be modeled using a single polygonal or polyhedral element. Attributes of the virtual element approach are highlighted through comparisons with Voronoi-cell lattice models, which provide discrete representations of material structure. The comparisons involve a suite of two-dimensional linear elastic problems: patch test, axisymmetric circular inclusion problem, and the deformation of a three-phase composite. The simulations demonstrate the accuracy and flexibility of the virtual element method.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Virtual elements on agglomerated finite elements to increase the critical time step in elastodynamic simulations

In this paper, we use the first-order virtual element method (VEM) to investigate the effect of shape quality of polyhedra in the estimation of the critical time step for explicit three-dimensional elastodynamic finite element (FE) simulations. Low-quality finite elements are common when meshing realistic complex components, and while tetrahedral meshing technology is generally robust, meshing algorithms cannot guarantee high-quality meshes for arbitrary geometries or for non-water-tight computer-aided design models. For reliable simulations on such meshes, we consider FE meshes with tetrahedral and prismatic elements that have badly-shaped elements$-$tetrahedra with dihedral angles close to $0^\circ$ and $180^\circ$, and slender prisms with triangular faces that have short edges$-$and agglomerate such `bad' elements with neighboring elements to form a larger polyhedral virtual element. On each element, the element-eigenvalue inequality is used to estimate the critical time step. For a suite of illustrative finite element meshes with $ε$ being a mesh-coordinate parameter that leads to poor mesh quality, we show that adopting VEM on the agglomerated polyhedra yield critical time steps that are insensitive as $ε\rightarrow 0$. The significant reduction in solution time on meshes with agglomerated virtual elements vis-$à$-vis tetrahedral meshes is demonstrated through explicit dynamics simulations on a tapered beam.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Exact imposition of boundary conditions with distance functions in physics-informed deep neural networks

In this paper, we introduce a new approach based on distance fields to exactly impose boundary conditions in physics-informed deep neural networks. The challenges in satisfying Dirichlet boundary conditions in meshfree and particle methods are well-known. This issue is also pertinent in the development of physics informed neural networks (PINN) for the solution of partial differential equations. We introduce geometry-aware trial functions in artifical neural networks to improve the training in deep learning for partial differential equations. To this end, we use concepts from constructive solid geometry (R-functions) and generalized barycentric coordinates (mean value potential fields) to construct $ϕ$, an approximate distance function to the boundary of a domain. To exactly impose homogeneous Dirichlet boundary conditions, the trial function is taken as $ϕ$ multiplied by the PINN approximation, and its generalization via transfinite interpolation is used to a priori satisfy inhomogeneous Dirichlet (essential), Neumann (natural), and Robin boundary conditions on complex geometries. In doing so, we eliminate modeling error associated with the satisfaction of boundary conditions in a collocation method and ensure that kinematic admissibility is met pointwise in a Ritz method. We present numerical solutions for linear and nonlinear boundary-value problems over domains with affine and curved boundaries. Benchmark problems in 1D for linear elasticity, advection-diffusion, and beam bending; and in 2D for the Poisson equation, biharmonic equation, and the nonlinear Eikonal equation are considered. The approach extends to higher dimensions, and we showcase its use by solving a Poisson problem with homogeneous Dirichlet boundary conditions over the 4D hypercube. This study provides a pathway for meshfree analysis to be conducted on the exact geometry without domain discretization.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Extended virtual element method for two-dimensional linear elastic fracture

In this paper, we propose an eXtended Virtual Element Method (X-VEM) for two-dimensional linear elastic fracture. This approach, which is an extension of the standard Virtual Element Method (VEM), facilitates mesh-independent modeling of crack discontinuities and elastic crack-tip singularities on general polygonal meshes. For elastic fracture in the X-VEM, the standard virtual element space is augmented by additional basis functions that are constructed by multiplying standard virtual basis functions by suitable enrichment fields, such as asymptotic mixed-mode crack-tip solutions. The design of the X-VEM requires an extended projector that maps functions lying in the extended virtual element space onto a set spanned by linear polynomials and the enrichment fields. An efficient scheme to compute the mixed-mode stress intensity factors using the domain form of the interaction integral is described. The formulation permits integration of weakly singular functions to be performed over the boundary edges of the element. Numerical experiments are conducted on benchmark mixed-mode linear elastic fracture problems that demonstrate the sound accuracy and optimal convergence in energy of the proposed formulation.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

A nodal integration scheme for meshfree Galerkin methods using the virtual element decomposition

In this paper, we present a novel nodal integration scheme for meshfree Galerkin methods that draws on the mathematical framework of the virtual element method. We adopt linear maximum-entropy basis functions for the discretization of field variables, although the proposed scheme is applicable to any linear meshfree approximant. In our approach, the weak form integrals are nodally integrated using nodal representative cells that carry the nodal displacements and state variables such as strains and stresses. The nodal integration is performed using the virtual element decomposition, wherein the bilinear form is decomposed into a consistency part and a stability part that ensure consistency and stability of the method. The performance of the proposed nodal integration scheme is assessed through benchmark problems in linear and nonlinear analyses of solids for small displacements and small-strain kinematics. Numerical results are presented for linear elastostatics and linear elastodynamics, and viscoelasticity. We demonstrate that the proposed nodally integrated meshfree method is accurate, converges optimally, and is more reliable and robust than a standard cell-based Gauss integrated meshfree method.

0 citations0 reviewsNo code linkedOpen access
preprint2018arXiv

Orbital-enriched Flat-top Partition of Unity Method for the Schrödinger Eigenproblem

Quantum mechanical calculations require the repeated solution of a Schrödinger equation for the wavefunctions of the system. Recent work has shown that enriched finite element methods significantly reduce the degrees of freedom required to obtain accurate solutions. However, time to solution has been adversely affected by the need to solve a generalized eigenvalue problem and the ill-conditioning of associated systems matrices. In this work, we address both issues by proposing a stable and efficient orbital-enriched partition-of-unity method to solve the Schrödinger boundary-value problem in a parallelepiped unit cell subject to Bloch-periodic boundary conditions. In our proposed PUM, the three-dimensional domain is covered by overlapping patches, with a compactly-supported, non-negative weight function, that is identically equal to unity over some finite subset of its support associated with each patch. This so-called flat-top property provides a pathway to devise a stable approximation over the whole domain. On each patch, we use $p$-th degree orthogonal polynomials that ensure $p$-th order completeness, and in addition include eigenfunctions of the radial solution of the Schrödinger equation. Furthermore, we adopt a variational lumping approach to construct a block-diagonal overlap matrix that yields a standard eigenvalue problem and demonstrate accuracy, stability and efficiency of the method.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Convergence analysis of a locally accelerated preconditioned steepest descent method for Hermitian-definite generalized eigenvalue problems

By extending the classical analysis techniques due to Samokish, Faddeev and Faddeeva, and Longsine and McCormick among others, we prove the convergence of preconditioned steepest descent with implicit deflation (PSD-id) method for solving Hermitian-definite generalized eigenvalue problems. Furthermore, we derive a nonasymptotic estimate of the rate of convergence of the \psdid method. We show that with the proper choice of the shift, the indefinite shift-and-invert preconditioner is a locally accelerated preconditioner, and is asymptotically optimal that leads to superlinear convergence. Numerical examples are presented to verify the theoretical results on the convergence behavior of the \psdid method for solving ill-conditioned Hermitian-definite generalized eigenvalue problems arising from electronic structure calculations. While rigorous and full-scale convergence proofs of preconditioned block steepest descent methods in practical use still largely eludes us, we believe the theoretical results presented in this paper sheds light on an improved understanding of the convergence behavior of these block methods.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Gradient bounds for Wachspress coordinates on polytopes

We derive upper and lower bounds on the gradients of Wachspress coordinates defined over any simple convex d-dimensional polytope P. The bounds are in terms of a single geometric quantity h_*, which denotes the minimum distance between a vertex of P and any hyperplane containing a non-incident face. We prove that the upper bound is sharp for d=2 and analyze the bounds in the special cases of hypercubes and simplices. Additionally, we provide an implementation of the Wachspress coordinates on convex polyhedra using Matlab and employ them in a 3D finite element solution of the Poisson equation on a non-trivial polyhedral mesh. As expected from the upper bound derivation, the H^1-norm of the error in the method converges at a linear rate with respect to the size of the mesh elements.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Hybrid preconditioning for iterative diagonalization of ill-conditioned generalized eigenvalue problems in electronic structure calculations

The iterative diagonalization of a sequence of large ill-conditioned generalized eigenvalue problems is a computational bottleneck in quantum mechanical methods employing a nonorthogonal basis for {\em ab initio} electronic structure calculations. We propose a hybrid preconditioning scheme to effectively combine global and locally accelerated preconditioners for rapid iterative diagonalization of such eigenvalue problems. In partition-of-unity finite-element (PUFE) pseudopotential density-functional calculations, employing a nonorthogonal basis, we show that the hybrid preconditioned block steepest descent method is a cost-effective eigensolver, outperforming current state-of-the-art global preconditioning schemes, and comparably efficient for the ill-conditioned generalized eigenvalue problems produced by PUFE as the locally optimal block preconditioned conjugate-gradient method for the well-conditioned standard eigenvalue problems produced by planewave methods.

0 citations0 reviewsNo code linkedOpen access