Researcher profile

N. Sukumar

N. Sukumar contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 17 - UnverifiedVerification L1Unclaimed author
4works
0followers
6topics
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

4 published item(s)

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.

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.

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.

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.