Researcher profile

Nilima Nigam

Nilima Nigam contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
7works
0followers
8topics
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

7 published item(s)

preprint2025arXiv

At the intersection of Numerical Analysis and Spectral Geometry

How do the geometric properties of a domain impact the spectrum of an operator defined on it? How do we compute accurate and reliable approximations of these spectra? The former question is studied in spectral geometry, and the latter is a central concern in numerical analysis. In this short expository survey we revisit the process of eigenvalue approximation, from the perspective of computational spectral geometry. Over the years a multitude of methods -- for discretizing the operator and for the resultant discrete system -- have been developed and analyzed in the field of numerical analysis. High-accuracy and provably convergent discretization approaches can be used to examine the interplay between the spectrum of an operator and the geometric properties of the spatial domain or manifold it is defined on. While computations have been used to guide conjectures in spectral geometry, in recent years approximation-theoretic tools and validated computations are also being used as part of proof strategies in spectral geometry. Given a particular spectral feature of interest, should we discretize the original problem, or seek a reformulation? Of the many possible approximation strategies, which should we choose? These choices are inextricably linked to the objective: on the one hand, rapid, specialized methods are often ideal for conjecture formulation (prioritizing efficiency and accuracy), whereas schemes with guaranteed, computable error bounds are needed when computation is incorporated into a proof strategy. We also review instances where the demanding requirements of spectral geometry -- the need for rigorous error control or the robust calculation of higher eigenvalues -- motivate new developments in numerical analysis.

preprint2025arXiv

On the dynamical stability of skeletal muscle

There has been debate for over 70-years about whether active skeletal muscle is dynamically stable at lengths greater than its optimal length. The stability of computational muscle models is a critical issue, as it directly affects our ability to simulate muscle deformation across different operating lengths, especially at lengths where muscles are known to remain functional despite model-predicted instabilities. In this study, we revisit the question of dynamical stability of ODE-based models of skeletal muscle. In particular, we investigate whether activation-independent tissue properties can provide stability to contractions along the dip region of the total force-length curve. First, using a combination of analytical tools (eigenvalue analysis and non-dimensionalization) and numerical simulations, we confirm that traditional Hill-type muscle models can display divergent dynamics in this region. Then, we propose a stabilized version of a 1D Hill-type muscle model that incorporates the 3D nature of skeletal muscle deformation. This results in a completely convex force-length relationship that can bring robustness to numerical simulations, while preserving the computational efficiency of 1D models. Our findings suggest that activation-independent intrinsic mechanical properties of muscle are sufficient to stabilize contractions even in the dip region, offering new insight into how muscles maintain functional integrity during active stretch.

preprint2022arXiv

A three-dimensional model of skeletal muscle tissues

Skeletal muscles are living tissues that can undergo large deformations in short periods of time and that can be activated to produce force. In this paper we use the principles of continuum mechanics to propose a dynamic, fully non-linear, and three-dimensional model to describe the deformation of these tissues. We model muscles as a fibre-reinforced composite and transversely isotropic material. We introduce a flexible computational framework to approximate the deformations of skeletal muscle to provide new insights into the underlying mechanics of these tissues. The model parameters and mechanical properties are obtained through experimental data and can be specified locally. A semi-implicit in time, conforming finite element in space scheme is used to approximate the solutions to the governing nonlinear dynamic model. We provide a series of numerical experiments demonstrating the application of this framework to relevant problems in biomechanics, and also discuss questions around model validation.

preprint2020arXiv

A simple extrapolation method for clustered eigenvalues

This paper introduces a simple variant of the power method. It is shown analytically and numerically to accelerate convergence to the dominant eigenvalue/eigenvector pair; and, it is particularly effective for problems featuring a small spectral gap. The introduced method is a one-step extrapolation technique that uses a linear combination of current and previous update steps to form a better approximation of the dominant eigenvector. The provided analysis shows the method converges exponentially with respect to the ratio between the two largest eigenvalues, which is also approximated during the process. An augmented technique is also introduced, and is shown to stabilize the early stages of the iteration. Numerical examples are provided to illustrate the theory and demonstrate the methods.

preprint2019arXiv

Optimization of Steklov-Neumann eigenvalues

This paper examines the Laplace equation with mixed boundary conditions, the Neumann and Steklov boundary conditions. This models a container with holes in it, like a pond filled with water but partly covered by immovable pieces on the surface. The main objective is to determine the right extent of the covering pieces, so that any shock inside the container yields a resonance. To this end, an algorithm is developed which uses asymptotic formulas concerning perturbations of the partitioning of the boundary pieces. Proofs for these formulas are established. Furthermore, this paper displays some results concerning bounds and examples with regards to the governing problem.

preprint2019arXiv

Revisiting the Jones eigenproblem in fluid-structure interaction

The Jones eigenvalue problem first described by D.S. Jones in 1983 concerns unusual modes in bounded elastic bodies: time-harmonic displacements whose tractions and normal components are both identically zero on the boundary. This problem is usually associated with a lack of unique solvability for certain models of fluid-structure interaction. The boundary conditions in this problem appear, at first glance, to rule out {\it any} non-trivial modes unless the domain possesses significant geometric symmetries. Indeed, Jones modes were shown to not be possible in most $C^\infty$ domains (see article by T. Hargé 1990). However, we should in this paper that while the existence of Jones modes sensitively depends on the domain geometry, such modes {\it do} exist in a broad class of domains. This paper presents the first detailed theoretical and computational investigation of this eigenvalue problem in Lipschitz domains. We also analytically demonstrate Jones modes on some simple geometries.

preprint2010arXiv

Numerical integration for high order pyramidal finite elements

We examine the effect of numerical integration on the convergence of high order pyramidal finite element methods. Rational functions are indispensable to the construction of pyramidal interpolants so the conventional treatment of numerical integration, which requires that the finite element approximation space is piecewise polynomial, cannot be applied. We develop an analysis that allows the finite element approximation space to include rational functions and show that despite this complication, conventional rules of thumb can still be used to select appropriate quadrature methods on pyramids. Along the way, we present a new family of high order pyramidal finite elements for each of the spaces of the de Rham complex.