BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Debasish Roy

Debasish Roy contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
Methodologymath-phmath.MPApplicationsComputationComputational Engineering, Finance, and Sciencecond-mat.softcond-mat.stat-mechmath.NANumerical Analysis

Trust snapshot

Quick read

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

11 published item(s)

preprint2022arXiv

A novel four-field mixed variational approach to Kirchhoff rods implemented with finite element exterior calculus

A four-field mixed variational principle is proposed for large deformation analysis of Kirchhoff rods with a $C^0$ mixed FE approximations. The core idea behind the approach is to introduce a one-parameter family of points (the centerline) and a separate one-parameter family of orthonormal frames (the Cartan moving frame) that are specified independently. The curvature and torsion of the curve are related to the relative rotation of neighboring frames. The relationship between the frame and the centerline is then enforced at the solution step using a Lagrange multiplier (which plays the role of section force). Well known frames like the Frenet-Serret are defined only using the centerline, which demands higher-order smoothness for the centerline approximation. Decoupling the frame from the position vector of the base curve leads to a description of torsion and curvature that is independent of the position information, thus allowing for simpler interpolations. This approach is then cast in the language of exterior calculus. In this framework, the strain energy may be interpreted as a differential form leading to the natural force-displacement (velocity) pairing. The four-field mixed variational principle we propose has frame, differential forms, and position vector as input arguments. While the position vector is interpolated linearly, the frames are interpolated as piecewise geodesics on the rotation group. Similarly, consistent interpolation schemes are also adopted to obtain finite dimensional approximations for other differential forms aswell. Using these discrete approximations, a discrete mixed variational principle is laid out which is then numerically extremized. The discrete approximation is then applied to benchmark problems, our numerical studies reveal an impressive performance of the proposed method without numerical instabilities or locking.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Geometric Thermodynamics of Strain-Induced Crystallization in Polymers

Going beyond the classical Gaussian approximation of Einstein's fluctuation theory, Ruppeiner gave it a Riemannian geometric structure with an entropic metric. This yielded a fundamental quantity - the Riemannian curvature, which was used to extract information on the nature of interactions between molecules in fluids, ideal gases and other open systems. In this article, we examine the implications of this curvature in a non-equilibrium thermodynamic system where relaxation is sufficiently slow so as not to invalidate the local equilibrium hypothesis. The non-equilibrium system comprises of a rubbery polymer undergoing strain induced crystallization. The curvature is found to impart information on a spurious isochoric energy arising from the conformational stretching of already crystallized segments. This unphysical component perhaps arises as the crystallized manifold is considered Euclidean with the stretch measures defined via the Euclidean metric. The thermodynamic state associated with curvature is the key to determine the isochoric stretch and hence the spurious energy. We determine this stretch and propose a form for the spurious free energy that must be removed from the total energy in order that the correct stresses are recovered.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Geometrically adapted Langevin dynamics for Markov chain Monte Carlo simulations

Markov Chain Monte Carlo (MCMC) is one of the most powerful methods to sample from a given probability distribution, of which the Metropolis Adjusted Langevin Algorithm (MALA) is a variant wherein the gradient of the distribution is used towards faster convergence. However, being set up in the Euclidean framework, MALA might perform poorly in higher dimensional problems or in those involving anisotropic densities as the underlying non-Euclidean aspects of the geometry of the sample space remain unaccounted for. We make use of concepts from differential geometry and stochastic calculus on Riemannian manifolds to geometrically adapt a stochastic differential equation with a non-trivial drift term. This adaptation is also referred to as a stochastic development. We apply this method specifically to the Langevin diffusion equation and arrive at a geometrically adapted Langevin dynamics. This new approach far outperforms MALA, certain manifold variants of MALA, and other approaches such as Hamiltonian Monte Carlo (HMC), its adaptive variant the no-U-turn sampler (NUTS) implemented in Stan, especially as the dimension of the problem increases where often GALA is actually the only successful method. This is evidenced through several numerical examples that include parameter estimation of a broad class of probability distributions and a logistic regression problem.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

A mixed variational principle in nonlinear elasticity using Cartan's moving frame and implementation with finite element exterior calculus

This article offers a new perspective for the mechanics of solids using moving Cartan's frame, specifically discussing a mixed variational principle in non-linear elasticity. We treat quantities defined on the co-tangent bundles of reference and deformed configurations as additional unknowns. Such a treatment invites compatibility of the fields with base-space (configurations of the body), so that the configuration can be realised as a subset of the Euclidean space. We first rewrite the metric and connection using differential forms, which are further utilised to write the deformation gradient and Cauchy-Green deformation tensor in terms of frame and co-frame fields. The geometric understanding of stress as a co-vector valued 2-form fits squarely within our overall program. We show that, for a hyperelastic solid, an equation similar to the Doyle-Erciksen formula may be written for the co-vector part of stress. Using these, we write a mixed energy functional in terms of differential forms, whose extremum leads to the compatibility of deformation, constitutive rules and equations of equilibrium. Finite element exterior calculus is then utilised to construct a finite dimensional approximation for the differential forms appearing in the variational principle. These approximations are then used to construct a discrete functional which is then numerically extremised. This discertization leads to a mixed method as it uses independent approximations for differential forms related to stress and deformation gradient. The mixed variational principle is then specialized for 2D case, whose discrete approximation is applied to problems in nonlinear elasticity. An important feature of our FE technique is the lack of additional stabilization. From the numerical study, it is found that the present discretization also does not suffer form locking and related convergence issues.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A geometric approach to Hu-Washizu variational principle in nonlinear elasticity

We discuss the Hu-Washizu (HW) variational principle from a geometric standpoint. The mainstay of the present approach is to treat quantities defined on the co-tangent bundles of reference and deformed configurations as primal. Such a treatment invites compatibility equations so that the base space (configurations of the solid body) could be realised as a subset of an Euclidean space. Cartan's method of moving frames and the associated structure equations establish this compatibility. Moreover, they permit us to write the metric and connection using 1-forms. With the mathematical machinery provided by differentiable manifolds, we rewrite the deformation gradient and Cauchy-Green deformation tensor in terms of frame and co-frame fields. The geometric understanding of stress as a co-vector valued 2-form fits squarely within our overall program. We also show that for a hyperelastic solid, an equation similar to the Doyle-Erciksen formula may be written for the co-vector part of the stress 2-form. Using this kinetic and kinematic understanding, we rewrite the HW functional in terms of frames and differential forms. Finally, we show that the compatibility of deformation, constitutive rules and equations of equilibrium are obtainable as Euler-Lagrange equations of the HW functional when varied with respect to traction 1-forms, deformation 1-forms and the deformation. This new perspective that involves the notion of kinematic closure precisely explicates the necessary geometrical restrictions on the variational principle, without which the deformed body may not be realized as a subset of the Euclidean space. It also provides a pointer to how these restrictions could be adjusted within a non-Euclidean setting.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Stochastic dynamical systems developed on Riemannian manifolds

We propose a method for developing the flows of stochastic dynamical systems, posed as Ito's stochastic differential equations, on a Riemannian manifold identified through a suitably constructed metric. The framework used for the stochastic development, viz. an orthonormal frame bundle that relates a vector on the tangent space of the manifold to its counterpart in the Euclidean space of the same dimension, is the same as that used for developing a standard Brownian motion on the manifold. Mainly drawing upon some aspects of the energetics so as to constrain the flow according to any known or prescribed conditions, we show how to expediently arrive at a suitable metric, thus briefly demonstrating the application of the method to a broad range of problems of general scientific interest. These include simulations of Brownian dynamics trapped in a potential well, a numerical integration scheme that reproduces the linear increase in the mean energy of conservative dynamical systems under additive noise and non-convex optimization. The simplicity of the method and the sharp contrast in its performance vis-a-vis the correspondent Euclidean schemes in our numerical work provide a compelling evidence to its potential.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A global optimization paradigm based on change of measures

A global optimization framework, acronymed COMBEO (Change OfMeasure Based Evolutionary Optimization), is proposed. An important aspect in the development is a set of derivative-free additive directional terms obtainable through a change of measures en route to the imposition of any stipulated conditions aimed at driving the realized design variables (particles) to the global optimum. The generalized setting offered by the new approach also enables several basic ideas, used with other global search methods such as the particle swarm or the differential evolution, to be rationally incorporated in the proposed setup via a change of measures. The global search may be further aided by imparting to the directional update terms additional layers of random perturbations such as scrambling and selection. Depending on the precise choice of the optimality conditions and the extent of random perturbation, the search can be readily rendered either greedy or more exploratory. As numerically demonstrated, the new proposal appears to provide for a more rational, more accurate and faster alternative to most available evolutionary optimization schemes, prominent amongst which are the differential evolution and the particle swarm methods.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

An Ensemble Kushner-Stratonovich (EnKS) Nonlinear Filter: Additive Particle Updates in Non-Iterative and Iterative Forms

Despite the cheap availability of computing resources enabling faster Monte Carlo simulations, the potential benefits of particle filtering in revealing accurate statistical information on the imprecisely known model parameters or modeling errors of dynamical systems, based on limited time series data, have not been quite realized. A major numerical bottleneck precipitating this under-performance, especially for higher dimensional systems, is the progressive particle impoverishment owing to weight collapse and the aim of the current work is to address this problem by replacing weight-based updates through additive ones. Thus, in the context of nonlinear filtering problems, a novel additive particle update scheme, in its non-iterative and iterative forms, is proposed based on manipulations of the innovation integral in the governing Kushner-Stratonovich equation. Numerical evidence for the identification of nonlinear and large dimensional dynamical systems indicates a substantively superior performance of the non- iterative version of the EnKS vis-à-vis most existing filters. The costlier iterative version, though conceptually elegant, mostly appears to effect a marginal improvement in the reconstruction accuracy over its non-iterative counterpart. Prominent in the reported numerical comparisons are variants of the Ensemble Kalman Filter (EnKF) that also use additive updates, albeit with many inherent limitations of a Kalman filter.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

An ensemble Kushner-Stratonovich-Poisson filter for recursive estimation in nonlinear dynamical systems

Despite the numerous applications that may be expeditiously modelled by counting processes, stochastic filtering strategies involving Poisson-type observations still remain somewhat poorly developed. In this work, we propose a Monte Carlo stochastic filter for recursive estimation in the context of linear/nonlinear dynamical systems with Poisson-type measurements. A key aspect of the present development is the filter-update scheme, derived from an ensemble approximation of the time-discretized nonlinear filtering equation, modified to account for Poisson-type measurements. Specifically, the additive update through a gain-like correction term, empirically approximated from the innovation integral in the filtering equation, eliminates the problem of particle collapse encountered in many conventional particle filters. Through a few numerical demonstrations, the versatility of the proposed filter is brought forth, first with application to filtering problems with diffusive or Poisson-type measurements and then to an automatic control problem wherein the extremization of the associated cost functional is achieved simply by an appropriate redefinition of the innovation process.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Evolutionary global optimization posed as a randomly perturbed martingale problem and applied to parameter recovery of chaotic oscillators

A new global stochastic search, guided mainly through derivative-free directional information computable from the sample statistical moments of the design variables within a Monte Carlo setup, is proposed. The search is aided by imparting to a directional update term, which parallels the conventional Gateaux derivative used in a local search for the extrema of smooth cost functionals, additional layers of random perturbations referred to as 'coalescence' and 'scrambling'. A selection scheme, constituting yet another avenue for random perturbation, completes the global search. The direction-driven nature of the search is manifest in the local extremization and coalescence components, which are posed as martingale problems that yield gain-like update terms upon discretization. As anticipated and numerically demonstrated, to a limited extent, against the problem of parameter recovery given the chaotic response histories of a couple of nonlinear oscillators, the proposed method apparently provides for a more rational, more accurate and faster alternative to most available evolutionary schemes, prominently the particle swarm optimization.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Iterated Gain-based Stochastic Filters for Dynamic System Identification: Annealing-type Iterations and the Filter Bank

A novel form of nonlinear stochastic filtering employing an annealing-type iterative update scheme, aided by the introduction of an artificial diffusion parameter and based on the Gaussian sum approximations of the prior and posterior densities, is presented. The proposed Monte Carlo filter bank conforms in structure to the parent nonlinear filtering (Kushner-Stratonovich) equation, as reflected in the additive gain-based updates, and possesses excellent mixing properties enabling better explorations of the phase space of the state vector. The performance of the filter bank, presently assessed against a few carefully chosen numerical examples, provide ample evidence of its substantively improved performance in terms of filter convergence and estimation accuracy vis-à-vis a few other competing filters especially in higher dimensional dynamic system identification problems including cases that may demand estimating relatively minor variations in the parameter values from their reference states.

0 citations0 reviewsNo code linkedOpen access