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

Nilima Nigam

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

ResearcherUnclaimed source record
math.NANumerical AnalysisTissues and Organsmath.DSQuantitative Methodsmath-phmath.APmath.FAmath.MPmath.PRmath.SPphysics.med-ph

Catalog footprint

What is connected

20works
12topics
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

20 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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Nearly radial Neumann eigenfunctions on symmetric domains

We study the existence of Neumann eigenfunctions which do not change sign on the boundary of some special domains. We show that eigenfunctions which are strictly positive on the boundary exist on regular polygons with at least 5 sides, while on equilateral triangles and cubes it is not even possible to find an eigenfunction which is nonnegative on the boundary. We use analytic methods combined with symmetry arguments to prove the result for polygons with six or more sides. The case for the regular pentagon is harder. We develop a validated numerical method to prove this case, which involves iteratively bounding eigenvalues for a sequence of subdomains of the triangle. We use a learning algorithm to find and optimize this sequence of subdomains, making it straightforward to check our computations with standard software.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

On regularizations of the delta distribution

In this article we consider regularizations of the Dirac delta distribution with applications to prototypical elliptic and hyperbolic partial differential equations (PDEs). We study the convergence of a sequence of distributions $\mathcal{S}_H$ to a singular term $\mathcal{S}$ as a parameter $H$ (associated with the {support size} of $\mathcal{S}_H$) shrinks to zero. We characterize this convergence in both the weak-$\ast$ topology of distributions, as well as in a weighted Sobolev norm. These notions motivate a framework for constructing regularizations of the delta distribution that includes a large class of existing methods in the literature. This framework allows different regularizations to be compared. The convergence of solutions of PDEs with these regularized source terms is then studied in various topologies such as pointwise convergence on a deleted neighborhood and weighted Sobolev norms. We also examine the lack of symmetry in tensor product regularizations and effects of dissipative error in hyperbolic problems.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Structural Changes of Active Skeletal Muscles: Modelling, Validation and Numerical Experiments

The purpose of this study was to report numerical validation of a 3D finite element model of contracting muscle. The model was based on continuum theory for fibre-reinforced composite materials. Here we simulated contractions for an idealized medial gastrocnemius muscle in man, using the model. Simulations were performed to test the force-length relation of the whole muscle, to evaluate the changes in internal fascicle geometry during contractions, and to assess the importance of material formulations for the aponeurosis and tendon. The simulation results were compared to previously published experimental values. The force-length profile for the whole muscle showed a realistic profile. As the muscle contracted the fascicles curved into S-shaped trajectories and curled around 3D paths, both of which matched previous experimental findings. As the fascicles shortened they increased in their cross-sectional area, but this increase was asymmetric with the smaller increase occurring within the fascicle-plane: the Poisson's ratio in this plane matched that previously shown from ultrasound imaging. The distribution of strains in the aponeurosis and tendon was shown to be a function of their material properties. This study demonstrated that the model could replicate realistic patterns of whole muscle-force, and changes to the internal muscle geometry, and so will be useful for testing mechanisms that affect the structural changes within contracting muscle.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A boundary integral algorithm for the Laplace Dirichlet-Neumann mixed eigenvalue problem

We present a novel integral-equation algorithm for evaluation of Zaremba eigenvalues and eigenfunctions}, that is, eigenvalues and eigenfunctions of the Laplace operator with mixed Dirichlet-Neumann boundary conditions; of course, (slight modifications of) our algorithms are also applicable to the pure Dirichlet and Neumann eigenproblems. Expressing the eigenfunctions by means of an ansatz based on the single layer boundary operator, the Zaremba eigenproblem is transformed into a nonlinear equation for the eigenvalue $μ$. For smooth domains the singular structure at Dirichlet-Neumann junctions is incorporated as part of our corresponding numerical algorithm---which otherwise relies on use of the cosine change of variables, trigonometric polynomials and, to avoid the Gibbs phenomenon that would arise from the solution singularities, the Fourier Continuation method (FC). The resulting numerical algorithm converges with high order accuracy without recourse to use of meshes finer than those resulting from the cosine transformation. For non-smooth (Lipschitz) domains, in turn, an alternative algorithm is presented which achieves high-order accuracy on the basis of graded meshes. In either case, smooth or Lipschitz boundary, eigenvalues are evaluated by searching for zero minimal singular values of a suitably stabilized discrete version of the single layer operator mentioned above. (The stabilization technique is used to enable robust non-local zero searches.) The resulting methods, which are fast and highly accurate for high- and low-frequencies alike, can solve extremely challenging two-dimensional Dirichlet, Neumann and Zaremba eigenproblems with high accuracies in short computing times---enabling, in particular, evaluation of thousands of eigenvalues and corresponding eigenfunctions for a given smooth or non-smooth geometry with nearly full double-precision accuracy

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Integral equation methods for the Yukawa-Beltrami equation on the sphere

An integral equation method for solving the Yukawa-Beltrami equation on a multiply-connected sub-manifold of the unit sphere is presented. A fundamental solution for the Yukawa-Beltrami operator is constructed. This fundamental solution can be represented by conical functions. Using a suitable representation formula, a Fredholm equation of the second kind with a compact integral operator needs to be solved. The discretization of this integral equation leads to a linear system whose condition number is bounded independent of the size of the system. Several numerical examples exploring the properties of this integral equation are presented.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Fast integral equation methods for the Laplace-Beltrami equation on the sphere

Integral equation methods for solving the Laplace-Beltrami equation on the unit sphere in the presence of multiple "islands" are presented. The surface of the sphere is first mapped to a multiply-connected region in the complex plane via a stereographic projection. After discretizing the integral equation, the resulting dense linear system is solved iteratively using the fast multipole method for the 2D Coulomb potential in order to calculate the matrix-vector products. This numerical scheme requires only O(N) operations, where $N$ is the number of nodes in the discretization of the boundary. The performance of the method is demonstrated on several examples.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Convergence analysis of a multigrid algorithm for the acoustic single layer equation

We present and analyze a multigrid algorithm for the acoustic single layer equation in two dimensions. The boundary element formulation of the equation is based on piecewise constant test functions and we make use of a weak inner product in the multigrid scheme as proposed in \cite{BLP94}. A full error analysis of the algorithm is presented. We also conduct a numerical study of the effect of the weak inner product on the oscillatory behavior of the eigenfunctions for the Laplace single layer operator.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Boundary Integral Equations for the Laplace-Beltrami Operator

We present a boundary integral method, and an accompanying boundary element discretization, for solving boundary-value problems for the Laplace-Beltrami operator on the surface of the unit sphere $§$ in $\mathbb{R}^3$. We consider a closed curve ${\cal C}$ on ${\cal S}$ which divides ${\cal S}$ into two parts ${\cal S}_1$ and ${\cal S}_2$. In particular, ${\cal C} = \partial {\cal S}_1$ is the boundary curve of ${\cal S}_1$. We are interested in solving a boundary value problem for the Laplace-Beltrami operator in $§_2$, with boundary data prescribed on $\C$.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Geometric integration on spheres and some interesting applications

Geometric integration theory can be employed when numerically solving ODEs or PDEs with constraints. In this paper, we present several one-step algorithms of various orders for ODEs on a collection of spheres. To demonstrate the versatility of these algorithms, we present representative calculations for reduced free rigid body motion (a conservative ODE) and a discretization of micromagnetics (a dissipative PDE). We emphasize the role of isotropy in geometric integration and link numerical integration schemes to modern differential geometry through the use of partial connection forms; this theoretical framework generalizes moving frames and connections on principal bundles to manifolds with nonfree actions.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

On the well-posedness of the stochastic Allen-Cahn equation in two dimensions

White noise-driven nonlinear stochastic partial differential equations (SPDEs) of parabolic type are frequently used to model physical and biological systems in space dimensions d = 1,2,3. Whereas existence and uniqueness of weak solutions to these equations are well established in one dimension, the situation is different for d \geq 2. Despite their popularity in the applied sciences, higher dimensional versions of these SPDE models are generally assumed to be ill-posed by the mathematics community. We study this discrepancy on the specific example of the two dimensional Allen-Cahn equation driven by additive white noise. Since it is unclear how to define the notion of a weak solution to this equation, we regularize the noise and introduce a family of approximations. Based on heuristic arguments and numerical experiments, we conjecture that these approximations exhibit divergent behavior in the continuum limit. The results strongly suggest that a series of published numerical studies are problematic: shrinking the mesh size in these simulations does not lead to the recovery of a physically meaningful limit.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The cellular dynamics of bone remodeling: a mathematical model

The mechanical properties of vertebrate bone are largely determined by a process which involves the complex interplay of three different cell types. This process is called {\it bone remodeling}, and occurs asynchronously at multiple sites in the mature skeleton. The cells involved are bone resorbing osteoclasts, bone matrix producing osteoblasts and mechanosensing osteocytes. These cells communicate with each other by means of autocrine and paracrine signaling factors and operate in complex entities, the so-called bone multicellular units (BMU). To investigate the BMU dynamics in silico, we develop a novel mathematical model resulting in a system of nonlinear partial differential equations with time delays. The model describes the osteoblast and osteoclast populations together with the dynamics of the key messenger molecule RANKL and its decoy receptor OPG. Scaling theory is used to address parameter sensitivity and predict the emergence of pathological remodeling regimes. The model is studied numerically in one and two space dimensions using finite difference schemes in space and explicit delay equation solvers in time. The computational results are in agreement with in vivo observations and provide new insights into the role of the RANKL/OPG pathway in the spatial regulation of bone remodeling.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

High-order finite elements on pyramids. II: unisolvency and exactness

We present degrees of freedom to accompany the approximation spaces already presented in a companion paper and thus complete the definition of families of high-order conforming finite elements on pyramids for the spaces of the de Rham complex. We prove that the elements are unisolvent; are compatible with conventional tetrahedral and hexahedral elements; satisfy a commuting diagram property and contain high-degree polynomials. We also tabulate shape functions for each element.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access