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

Krishna Garikipati

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

ResearcherUnclaimed source record
cond-mat.mtrl-sciphysics.comp-phBiological PhysicsMachine Learningmath.NANumerical Analysismath-phmath.MPSubcellular ProcessesTissues and OrgansCell BehaviorComputational Engineering, Finance, and Sciencecond-mat.softphysics.app-phphysics.chem-phPopulations and Evolution

Catalog footprint

What is connected

19works
16topics
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

19 published item(s)

preprint2022arXiv

A fourth-order phase-field fracture model: Formulation and numerical solution using a continuous/discontinuous Galerkin method

Modeling crack initiation and propagation in brittle materials is of great importance to be able to predict sudden loss of load-carrying capacity and prevent catastrophic failure under severe dynamic loading conditions. Second-order phase-field fracture models have gained wide adoption given their ability to capture the formation of complex fracture patterns, e.g. via crack merging and branching, and their suitability for implementation within the context of the conventional finite element method. Higher-order phase-field models have also been proposed to increase the regularity of the exact solution and thus increase the spatial convergence rate of its numerical approximation. However, they require special numerical techniques to enforce the necessary continuity of the phase field solution. In this paper, we derive a fourth-order phase-field model of fracture in two independent ways; namely, from Hamilton's principle and from a higher-order micromechanics-based approach. The latter approach is novel, and provides a physical interpretation of the higher-order terms in the model. In addition, we propose a continuous/discontinuous Galerkin (C/DG) method for use in computing the approximate phase-field solution. This method employs Lagrange polynomial shape functions to guarantee $C^0$-continuity of the solution at inter-element boundaries, and enforces the required $C^1$ regularity with the aid of additional variational and interior penalty terms in the weak form. The phase-field equation is coupled with the momentum balance equation to model dynamic fracture problems in hyper-elastic materials. Two benchmark problems are presented to compare the numerical behavior of the C/DG method with mixed finite element methods.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

A heteroencoder architecture for prediction of failure locations in porous metals using variational inference

In this work we employ an encoder-decoder convolutional neural network to predict the failure locations of porous metal tension specimens based only on their initial porosities. The process we model is complex, with a progression from initial void nucleation, to saturation, and ultimately failure. The objective of predicting failure locations presents an extreme case of class imbalance since most of the material in the specimens do not fail. In response to this challenge, we develop and demonstrate the effectiveness of data- and loss-based regularization methods. Since there is considerable sensitivity of the failure location to the particular configuration of voids, we also use variational inference to provide uncertainties for the neural network predictions. We connect the deterministic and Bayesian convolutional neural networks at a theoretical level to explain how variational inference regularizes the training and predictions. We demonstrate that the resulting predicted variances are effective in ranking the locations that are most likely to fail in any given specimen.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Machine Learning in Heterogeneous Porous Materials

The "Workshop on Machine learning in heterogeneous porous materials" brought together international scientific communities of applied mathematics, porous media, and material sciences with experts in the areas of heterogeneous materials, machine learning (ML) and applied mathematics to identify how ML can advance materials research. Within the scope of ML and materials research, the goal of the workshop was to discuss the state-of-the-art in each community, promote crosstalk and accelerate multi-disciplinary collaborative research, and identify challenges and opportunities. As the end result, four topic areas were identified: ML in predicting materials properties, and discovery and design of novel materials, ML in porous and fractured media and time-dependent phenomena, Multi-scale modeling in heterogeneous porous materials via ML, and Discovery of materials constitutive laws and new governing equations. This workshop was part of the AmeriMech Symposium series sponsored by the National Academies of Sciences, Engineering and Medicine and the U.S. National Committee on Theoretical and Applied Mechanics.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Bayesian neural networks for weak solution of PDEs with uncertainty quantification

Solving partial differential equations (PDEs) is the canonical approach for understanding the behavior of physical systems. However, large scale solutions of PDEs using state of the art discretization techniques remains an expensive proposition. In this work, a new physics-constrained neural network (NN) approach is proposed to solve PDEs without labels, with a view to enabling high-throughput solutions in support of design and decision-making. Distinct from existing physics-informed NN approaches, where the strong form or weak form of PDEs are used to construct the loss function, we write the loss function of NNs based on the discretized residual of PDEs through an efficient, convolutional operator-based, and vectorized implementation. We explore an encoder-decoder NN structure for both deterministic and probabilistic models, with Bayesian NNs (BNNs) for the latter, which allow us to quantify both epistemic uncertainty from model parameters and aleatoric uncertainty from noise in the data. For BNNs, the discretized residual is used to construct the likelihood function. In our approach, both deterministic and probabilistic convolutional layers are used to learn the applied boundary conditions (BCs) and to detect the problem domain. As both Dirichlet and Neumann BCs are specified as inputs to NNs, a single NN can solve for similar physics, but with different BCs and on a number of problem domains. The trained surrogate PDE solvers can also make interpolating and extrapolating (to a certain extent) predictions for BCs that they were not exposed to during training. Such surrogate models are of particular importance for problems, where similar types of PDEs need to be repeatedly solved for many times with slight variations. We demonstrate the capability and performance of the proposed framework by applying it to steady-state diffusion, linear elasticity, and nonlinear elasticity.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Multiscale modeling meets machine learning: What can we learn?

Machine learning is increasingly recognized as a promising technology in the biological, biomedical, and behavioral sciences. There can be no argument that this technique is incredibly successful in image recognition with immediate applications in diagnostics including electrophysiology, radiology, or pathology, where we have access to massive amounts of annotated data. However, machine learning often performs poorly in prognosis, especially when dealing with sparse data. This is a field where classical physics-based simulation seems to remain irreplaceable. In this review, we identify areas in the biomedical sciences where machine learning and multiscale modeling can mutually benefit from one another: Machine learning can integrate physics-based knowledge in the form of governing equations, boundary conditions, or constraints to manage ill-posted problems and robustly handle sparse and noisy data; multiscale modeling can integrate machine learning to create surrogate models, identify system dynamics and parameters, analyze sensitivities, and quantify uncertainty to bridge the scales and understand the emergence of function. With a view towards applications in the life sciences, we discuss the state of the art of combining machine learning and multiscale modeling, identify applications and opportunities, raise open questions, and address potential challenges and limitations. We anticipate that it will stimulate discussion within the community of computational mechanics and reach out to other disciplines including mathematics, statistics, computer science, artificial intelligence, biomedicine, systems biology, and precision medicine to join forces towards creating robust and efficient models for biological systems.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

System inference for the spatio-temporal evolution of infectious diseases: Michigan in the time of COVID-19

We extend the classical SIR model of infectious disease spread to account for time dependence in the parameters, which also include diffusivities. The temporal dependence accounts for the changing characteristics of testing, quarantine and treatment protocols, while diffusivity incorporates a mobile population. This model has been applied to data on the evolution of the COVID-19 pandemic in the US state of Michigan. For system inference, we use recent advances; specifically our framework for Variational System Identification (Wang et al., Comp. Meth. App. Mech. Eng., 356, 44-74, 2019; arXiv:2001.04816 [cs.CE]) as well as Bayesian machine learning methods.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Variational Extrapolation of Implicit Schemes for General Gradient Flows

We introduce a class of unconditionally energy stable, high order accurate schemes for gradient flows in a very general setting. The new schemes are a high order analogue of the minimizing movements approach for generating a time discrete approximation to a gradient flow by solving a sequence of optimization problems. In particular, each step entails minimizing the associated energy of the gradient flow plus a movement limiter term that is, in the classical context of steepest descent with respect to an inner product, simply quadratic. A variety of existing unconditionally stable numerical methods can be recognized as (typically just first order accurate in time) minimizing movement schemes for their associated evolution equations, already requiring the optimization of the energy plus a quadratic term at every time step. Therefore, our approach gives a painless way to extend these to high order accurate in time schemes while maintaining their unconditional stability. In this sense, it can be viewed as a variational analogue of Richardson extrapolation.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Second Order Threshold Dynamics Schemes for Two Phase Motion by Mean Curvature

The threshold dynamics algorithm of Merriman, Bence, and Osher is only first order accurate in the two-phase setting. Its accuracy degrades further to half order in the multi-phase setting, a shortcoming it has in common with other related, more recent algorithms such as the equal surface tension version of the Voronoi implicit interface method. As a first, rigorous step in addressing this shortcoming, we present two different second order accurate versions of two-phase threshold dynamics. Unlike in previous efforts in this direction, we present careful consistency calculations for both of our algorithms. The first algorithm is consistent with its limit (motion by mean curvature) up to second order in any space dimension. The second achieves second order accuracy only in dimension two, but comes with a rigorous stability guarantee (unconditional energy stability) in any dimension -- a first for high order schemes of its type.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A comparison of Redlich-Kister polynomial and cubic spline representations of the chemical potential in phase field computations

Free energies play a central role in many descriptions of equilibrium and non-equilibrium properties of solids. Continuum partial differential equations (PDEs) of atomic transport, phase transformations and mechanics often rely on first and second derivatives of a free energy function. The stability, accuracy and robustness of numerical methods to solve these PDEs are sensitive to the particular functional representations of the free energy. In this communication we investigate the influence of different representations of thermodynamic data on phase field computations of diffusion and two-phase reactions in the solid state. First-principles statistical mechanics methods were used to generate realistic free energy data for HCP titanium with interstitially dissolved oxygen. While Redlich-Kister polynomials have formed the mainstay of thermodynamic descriptions of multi-component solids, they require high order terms to fit oscillations in chemical potentials around phase transitions. Here we demonstrate that high fidelity fits to rapidly fluctuating free energy functions are obtained with spline functions. Spline functions that are many degrees lower than Redlich-Kister polynomials provide equal or superior fits to chemical potential data and, when used in phase field computations, result in solution times approaching an order of magnitude speed up relative to the use of Redlich-Kister polynomials.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A three dimensional field formulation, and isogeometric solutions to point and line defects using Toupin's theory of gradient elasticity at finite strains

We present a field formulation for defects that draws from the classical representation of the cores as force dipoles. We write these dipoles as singular distributions. Exploiting the key insight that the variational setting is the only appropriate one for the theory of distributions, we arrive at universally applicable weak forms for defects in nonlinear elasticity. Remarkably, the standard, Galerkin finite element method yields numerical solutions for the elastic fields of defects that, when parameterized suitably, match very well with classical, linearized elasticity solutions. The true potential of our approach, however, lies in its easy extension to generate solutions to elastic fields of defects in the regime of nonlinear elasticity, and even more notably for Toupin's theory of gradient elasticity at finite strains(Arch. Rat. Mech. Anal., 11, 385, 1962). In computing these solutions we adopt recent numerical work on an isogeometric analytic framework that enabled the first three-dimensional solutions to general boundary value problems of Toupin's theory (Rudraraju et al. Comp. Meth. App. Mech. Engr., 278, 705, 2014). We first present exhaustive solutions to point defects, edge and screw dislocations, and a study on the energetics of interacting dislocations. Then, to demonstrate the generality and potential of our treatment, we apply it to other complex dislocation configurations, including loops and low-angle grain boundaries.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A variational treatment of material configurations with application to interface motion and microstructural evolution

We present a unified variational treatment of evolving configurations in crystalline solids with microstructure. The crux of our treatment lies in the introduction of a vector configurational field. This field lies in the material, or configurational, manifold, in contrast with the traditional displacement field, which we regard as lying in the spatial manifold. We identify two distinct cases which describe (a) problems in which the configurational field's evolution is localized to a mathematically sharp interface, and (b) those in which the configurational field's evolution can extend throughout the volume. The first case is suitable for describing incoherent phase interfaces in polycrystalline solids, and the latter is useful for describing smooth changes in crystal structure and naturally incorporates coherent (diffuse) phase interfaces. For sharp interfaces that are out-of-equilibrium, the second law of thermodynamics furnishes restrictions on the kinetic law for the interface velocity. The class of problems in which the material undergoes configurational changes between distinct, stable crystal structures are characterized by free energy density functions that are non-convex with respect to configurational strain. For physically meaningful solutions and mathematical well-posedness, it becomes necessary to incorporate interfacial energy. This we have done by introducing a configurational strain gradient dependence in the free energy density function following ideas laid out by Toupin (Arch. Rat. Mech. Anal. 11, 1962, 385-414). The variational treatment leads to a system of partial differential equations governing the configuration that is coupled with the traditional equations of nonlinear elasticity. Numerical examples are presented to demonstrate interface motion as well as evolving microstructures of crystal structures.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Mechano-chemical spinodal decomposition: A phenomenological theory of phase transformations in multi-component, crystalline solids

We present a phenomenological treatment of diffusion-driven martensitic phase transformations in multi-component crystalline solids that arise from non-convex free energies in mechanical and chemical variables. The treatment describes diffusional phase transformations that are accompanied by symmetry breaking structural changes of the crystal unit cell and reveals the importance of a mechano-chemical spinodal, defined as the region in strain-composition space where the free energy density function is non-convex. The approach is relevant to phase transformations wherein the structural order parameters can be expressed as linear combinations of strains relative to a high-symmetry reference crystal. The governing equations describing mechano-chemical spinodal decomposition are variationally derived from a free energy density function that accounts for interfacial energy via gradients of the rapidly varying strain and composition fields. A robust computational framework for treating the coupled, higher order diffusion and nonlinear strain gradient elasticity problems is presented. Because the local strains in an inhomogeneous, transforming microstructure can be finite, the elasticity problem must account for geometric nonlinearity. An evaluation of available experimental phase diagrams and first-principles free energies suggests mechano-chemical spinodal decomposition should occur in metal hydrides such as ZrH$_{2-2c}$. The rich physics that ensues is explored in several numerical examples in two and three dimensions and the relevance of the mechanism is discussed in the context of important electrode materials for Li-ion batteries and high temperature ceramics.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Multi-physics simulations of lithiation-induced stress in Li$_{\rm 1+x}$Ti$_2$O$_4$ electrode particles

Cubic spinel Li$_{\rm 1+x}$Ti$_2$O$_4$ is a promising electrode material as it exhibits a high lithium diffusivity and undergoes minimal changes in lattice parameters during lithiation and delithiation, thereby ensuring favorable cycleability. The present work is a multi-physics and multi-scale study of Li$_{\rm 1+x}$Ti$_2$O$_4$ that combines first principles computations of thermodynamic and kinetic properties with continuum scale modeling of lithiation-delithiation kinetics. Density functional theory calculations and statistical mechanics methods are used to calculate lattice parameters, elastic coefficients, thermodynamic potentials, migration barriers and Li diffusion coefficients. These quantities then inform a phase field framework to model the coupled chemo-mechanical evolution of electrode particles. Several case studies accounting for either homogeneous or heterogeneous nucleation are considered to explore the temporal evolution of maximum principle stress values, which serve to indicate stress localization and the potential for crack initiation, during lithiation and delithiation.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Perspectives on the mathematics of biological patterning and morphogenesis

A central question in developmental biology is how size and position are determined. The genetic code carries instructions on how to control these properties in order to regulate the pattern and morphology of structures in the developing organism. Transcription and protein translation mechanisms implement these instructions. However, this cannot happen without some manner of sampling of epigenetic information on the current patterns and morphological forms of structures in the organism. Any rigorous description of space- and time-varying patterns and morphological forms reduces to one among various classes of spatio-temporal partial differential equations. Reaction-transport equations represent one such class. Starting from simple Fickian diffusion, the incorporation of reaction, phase segregation and advection terms can represent many of the patterns seen in the animal and plant kingdoms. Morphological form, requiring the development of three-dimensional structure, also can be represented by these equations of mass transport, albeit to a limited degree. The recognition that physical forces play controlling roles in shaping tissues leads to the conclusion that (nonlinear) elasticity governs the development of morphological form. In this setting, inhomogeneous growth drives the elasticity problem. The combination of reaction-transport equations with those of elasto-growth makes accessible a potentially unlimited spectrum of patterning and morphogenetic phenomena in developmental biology. This perspective communication is a survey of the partial differential equations of mathematical physics that have been proposed to govern patterning and morphogenesis in developmental biology. Several numerical examples are included to illustrate these equations and the corresponding physics, with the intention of providing physical insight wherever possible.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Unconditionally stable, second-order accurate schemes for solid state phase transformations driven by mechano-chemical spinodal decomposition

We consider solid state phase transformations that are caused by free energy densities with domains of non-convexity in strain-composition space; we refer to the non-convex domains as mechano-chemical spinodals. The non-convexity with respect to composition and strain causes segregation into phases with different crystal structures. We work on an existing model that couples the classical Cahn-Hilliard model with Toupin's theory of gradient elasticity at finite strains. Both systems are represented by fourth-order, nonlinear, partial differential equations. The goal of this work is to develop unconditionally stable, second-order accurate time-integration schemes, motivated by the need to carry out large scale computations of dynamically evolving microstructures in three dimensions. We also introduce reduced formulations naturally derived from these proposed schemes for faster computations that are still second-order accurate. Although our method is developed and analyzed here for a specific class of mechano-chemical problems, one can readily apply the same method to develop unconditionally stable, second-order accurate schemes for any problems for which free energy density functions are multivariate polynomials of solution components and component gradients. Apart from an analysis and construction of methods, we present a suite of numerical results that demonstrate the schemes in action.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Three-dimensional iso-geometric solutions to general boundary value problems of Toupin's gradient elasticity theory at finite strains

We present, to the best of our knowledge, the first complete three-dimensional numerical solutions to a broad range of boundary value problems for a general theory of finite strain gradient elasticity. We have chosen for our work, Toupin's theory [Arch. Rat. Mech. Anal., 11(1), 385-414, 1962]-one of the more general formulations of strain gradient elasticity. Our framework has three crucial ingredients: The first is iso-geometric analysis [Hughes et al., Comp. Meth. App. Mech. Engrg., 194(39-41), 4135-4195, 2005], which we have adopted for its straightforward and robust representation of C1-continuity. The second is a weak treatment of the higher-order Dirichlet boundary conditions in the formulation, which control the development of strain gradients in the solution. The third ingredient is algorithmic (automatic) differentiation, which eliminates the need for linearization "by hand" of the rather complicated geometric and material nonlinearities in gradient elasticity at finite strains. We present a number of numerical solutions to demonstrate that the framework is applicable to arbitrary boundary value problems in three dimensions. We discuss size effects, the role of higher-order boundary conditions, and perhaps most importantly, the relevance of the framework to problems with elastic free energy density functions that are non-convex in strain space.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A computational study of stress fiber-focal adhesion dynamics governing cell contractility

We apply a recently developed model of cytoskeletal force generation to study a cell intrinsic contractility, as well as its response to external loading. The model is based on a non-equilibrium thermodynamic treatment of the mechano-chemistry governing force in the stress fiber-focal adhesion system. Our computational study suggests that the mechanical coupling between the stress fibers and focal adhesions leads to a complex, dynamic, mechano-chemical response. We collect the results in response maps whose regimes are distinguished by the initial geometry of the stress fiber-focal adhesion system, and by the external load on the cell. The results from our model connect qualitatively with recent studies on the force response of smooth muscle cells on arrays of polymeric microposts (Mann et al., Lab. on a Chip, 12, 731-740, 2012).

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Elastic free energy drives the shape of prevascular solid tumors

It is well established that the mechanical environment influences cell functions in health and disease. Here, we address how the mechanical environment influences tumor growth, in particular, the shape of solid tumors. In an in vitro tumor model, which isolates mechanical interactions between tumor cells and a hydrogel, we find that tumors grow as ellipsoids, resembling the same, oft-reported observation of in vivo tumors. Specifically, an oblate ellipsoidal tumor shape robustly occurs when the tumors grow in hydrogels that are stiffer than the tumors, but when they grow in more compliant hydrogels they remain closer to spherical in shape. Using large scale, nonlinear elasticity computations we show that the oblate ellipsoidal shape minimizes the elastic free energy of the tumor-hydrogel system. Having eliminated a number of other candidate explanations, we hypothesize that minimization of the elastic free energy is the reason for predominance of the experimentally observed ellipsoidal shape. This result may hold significance for explaining the shape progression of early solid tumors in vivo and is an important step in understanding the processes underlying solid tumor growth.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The mechano-chemistry of cytoskeletal force generation

In this communication, we propose a model to study the non-equilibrium process by which actin stress fibers develop force in contractile cells. The emphasis here is on the non-equilibrium thermodynamics, which is necessary to address the mechanics as well as the chemistry of dynamic cell contractility. In this setting we are able to develop a framework that relates (a) the dynamics of force generation within the cell and (b) the cell response to external stimuli to the chemical processes occurring within the cell, as well as to the mechanics of linkage between the stress fibers, focal adhesions and extra-cellular matrix.

0 citations0 reviewsNo code linkedOpen access