We model grain level phenomena linked with contact driven mechanical deformation. The three key phenomena being modeled are microstructure evolution, elasto-plastic deformation of the multi-grain aggregates, and contact modeling of the interface friction.
Graduate student: Prakarsh Pandey and Rahul Gulati. In collaboration with Prof. Izabela Szlufarska (MSE). Funded by ARO.
Friction Stir Welding is a relatively new solid-state welding process that has significant advantages over conventional fusion processes. We model plastic flow induced by the underlying thermo-mechanical processes and seek to gain a fundamental understanding of how and when defects (voids) are formed, and identify process parameters that allow for optimal control over the properties of the resulting weld.
Graduate student: Mohammad Ali Ansari. In collaboration with Prof. Frank Pfefferkorn and Prof. Mike Zinn. Funded by NSF.
The formation and constriction of membrane necks as a precursor to membrane fission is a fundamental step for a variety of biological processes including endocytosis, cytokinesis, etc. Using a curvilinear coordinates based Kirchhoff-Love thin shell kinematics framework, we model the deformational response and instabilities in various membrane geometries, study the suitability of Helfrich like material models, and determine the boundary conditions needed to induce preferred deformation modes.
Graduate student: Debabrata Auddya and Prakarsh Pandey. In collaboration with Prof. Padmini Rangamani (UCSD) and Prof. Krishna Garikipati (UM).
R. Vasan, S. Rudraraju, M. Akamatsu, K. Garikipati, P. Rangamani, ``A mechanical model reveals that non-axisymmetric buckling lowers the energy barrier associated with membrane neck constriction'', [arXiv](under review)
We model dendritic and grain level microstructure evolution, driven by alloy specific solidification models and thermo-mechanics. The numerical framework couples phase field based solidification models to mesoscale crystal plasticity and continuum plasticity models, and is geared towards handling realistic length scales and microstructures with multi-million degrees of freedom.
Graduate student: Kunal Bhagat. In collaboration with Prof. Dan Thoma (GIE). Funded by WARF.
Formation of multi-cell aggregates is a foundational process in the evolution of multicellular organisms. Beginning with a single cell or a small cluster, the growth of aggregates is driven by cell division, differentiation, migration and cell-cell interactions. Understanding the processes underlying the formation of these aggregates is central to many phenomena in cellular biology and physiology, including embryogenesis, wound healing and tissue engineering. In this work, we develop a model for cell growth, division and packing under soft constraints that arise from the deformability of the cells as well as of a membrane that encloses them. Our treatment falls within the framework of diffuse interface methods, under which each cell is represented by a scalar phase field and the zero level set of the phase field represents the cell membrane.
Graduate student: Debabrata Auddya. In collaboration with Prof. Krishna Garikipati (UM).
J. Jiang, K. Garikipati, S. Rudraraju``A diffuse interface framework for modelling the evolution of multi-cell aggregates as a soft packing problem due to growth and division of cells'', Bulletin of Mathematical Biology , 2019. [journal][arXiv]
Molluskan shells are an important model system for understanding mechano-morphological basis of the evolution of exoskeleton in invertebrates. Here we propose a three-dimensional computational framework coupling morphology, incremental surface growth by accretion, and morpho-elastic volume growth, to enable an improved representation of the growth and structural parameters controlling the evolution of shell structural ornamentations. We study the effects of growth rates, and of previously deposited shell geometries on the resulting modes of mantle deformation, and determine a ``phase diagram'' of morphogenesis in molluskan shells. Our main motivation for focusing on generic physical processes involved in development is that they may shape living beings in a predictive way and partly determine the spectrum of forms that have been and could have been generated during evolution.
In collaboration with Derek Moulton (Oxford), Alain Goriely (Oxford) and Krishna Garikipati (UM).
S. Rudraraju, D.E. Moulton, R. Chirat, A. Goriely, K. Garikipati, ''A computational framework for the morpho-elastic development of molluskan shells by surface and volume growth'', PLOS Computational Biology, 15(7): e1007213, 2019. [journal][arXiv]
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 mechanochemical 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. Specifically, phase transformations in shape memory alloys, electrode materials for Li-ion batteries and certain high-temperature ceramics are of interest.
In collaboration with Krishna Garikipati (UM) and Anton Van der Ven (UCSB).
S. Rudraraju, A van der Ven, K. Garikipati, ``Mechano-chemical spinodal decomposition: A phenomenological theory of phase transformations in multi-component, crystalline solids'', Nature npj Computational Materials, 2016, doi:10.1038/npjcompumats.2016.12. [journal] [arXiv]
Cubic spinel LiTiO is a promising electrode material because it exhibits a high lithium diffusivity and undergoes minimal changes in lattice parameters during lithiation and delithiation, thereby ensuring favorable cyclability. This work is a multiphysics and multiscale study of LiTiO that combines first-principles computations of thermodynamic and kinetic properties with continuum scale modeling of lithiation−delithiation kinetics. Several case studies explore the temporal evolution of peak stresses and the potential for crack initiation, during lithiation and delithiation.
In collaboration with Krishna Garikipati (UM) and Michael Falk (JHU).
T. Jiang, S. Rudraraju, A. Roy, A. Van der Ven, K. Garikipati, M. L. Falk, ``Multi-physics simulations of lithiation-induced stress in LiTiO electrode particles'', Journal of Physical Chemistry C, 2016, [journal] [arXiv]
Classical elasticity is scale invariant: it does not admit length scales intrinsic to the material. The only length scale which may manifest itself arises from the problem geometry. However, real materials do posses intrinsic length scales (i.e., interatomic distance, dislocation cells, grain sizes, etc). If the deformation varies on these scales, size effects may arise. One framework that is appropriate for treating such problems is the strain gradient formulation of elasticity. We have developed a comprehensive numerical framework for Toupin’s theory of nonlinear gradient elasticity and obtained the first complete three-dimensional numerical solutions to a broad range of boundary value problems in gradient elasticity at finite strain. The governing equation is a fourth order partial differential equation with complex boundary conditions. Studies conducted using the framework include gradient length scale effects on microscale/nanoscale structural problems, regularization of crack tip solutions and modeling of mechanically induced microstructure.
In collaboration with Krishna Garikipati (UM) and Anton Van der Ven (UCSB).
S. Rudraraju, A van der Ven, K. Garikipati, ``Three-dimensional iso-geometric solutions to general boundary value problems of Toupin's gradient elasticity theory at finite strains'', Computer Methods in Applied Mechanics and Engineering (CMAME), Vol 278, Pages 705-728, 2014. [journal][arXiv]
Z. Wang, S. Rudraraju, K. Garikipati,``A three dimensional field formulation, and isogeometric solutions to point and line defects using Toupin's theory of gradient elasticity at finite strains'', Journal of the Mechanics and Physics of Solids (JMPS), Vol. 94: 336-361, 2016, doi:10.1016/j.jmps.2016.03.028.[journal][arXiv]
The dynamical processes of chemistry, transport and mechanics that govern tumor growth can be broadly classified into three distinct spatial scales: the tumor scale, the cell-extracellular matrix (ECM) interactions and the sub-cellular processes. In this work, we focus on tumor scale investigations, and model the emergent morphology of tumors to gain system-level insights into their progression. The evolution of various constituents (cancer cells, stem cells, necrotic cells, nutrients, byproducts, etc) is modelled with reaction-advection-diffusion equations coupled with finite strain mechanics.
In collaboration with Kristen Mills (Max Planck Institute for Intelligent Systems, now at RPI), Ralf Kemkemer (Max Planck Institute for Intelligent Systems) and Krishna Garikipati (UM)
S. Rudraraju, K.L.Mills, R Kemkemer, K. Garikipati, ``Multiphysics Modeling of Reactions, Mass Transport and Mechanics of Tumor Growth'' Computer Models in Biomechanics, Pages 293-303, 2013. [journal]
Modeling crack propagation in fiber reinforced composites is more challenging than conventional brittle/ductile materials due to the presence of large process zones behind the crack tip. These process zones, also known as bridging zones or cohesive zones, are regions of complex fiber-matrix micromechanics which significantly enhance the fracture toughness. To address this modeling problem, we developed a numerical framework for cohesive crack propagation based on the variational multiscale method. It treats cracks as displacement discontinuities and represents them as fine scale variations with compact support on a coarse scale displacement field. This treatment leads to a multiscale model with a coarse scale problem enforcing the conservation of linear momentum and a fine scale problem representing the micromechanics (cohesive surface laws).
S. Rudraraju, A. Salvi, K. Garikipati, A. M. Waas, ``In-Plane fracture of laminated fiber reinforced composites with varying fracture resistance: experimental observations and numerical crack propagation simulations'', International Journal of Solids and Structures (IJSS), Vol 47, Issues 7-8, Pages 901-911, 2010. [journal] [arXiv]
S. Rudraraju, A. Salvi, K. Garikipati, A. M. Waas, ``Predictions of crack propagation using a variational multiscale approach and its application to fracture in laminated fiber reinforced composites'', Composite Structures, Vol 94, Issue 11, Pages 3336-3346, 2012. [journal]
S. Rudraraju, A. Salvi, K. Garikipati, A. M. Waas, ``Experimental observations and numerical simulations of curved crack propagation in laminated fiber composites'' Composites Science and Technology (CST), Vol 72, Issue 10, Pages 1064-1074, 2012. [journal]
S. Rudraraju, K. Garikipati, A.M. Waas, B.A. Bednarcyk, ``On the Theory and Numerical Simulation of Cohesive Crack Propagation with Application to Fiber-Reinforced Composites'', NASA Technical Publication (peer reviewed) (NASA/TP-2013-217431, E-18123), 2013. [NASA Archive]