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

Shiva Rudraraju

Shiva Rudraraju 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-scimath-phmath.MPphysics.comp-phcond-mat.softmath.NANumerical Analysisphysics.chem-phTissues and Organs

Catalog footprint

What is connected

10works
9topics
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

10 published item(s)

preprint2022arXiv

A numerical investigation of dimensionless numbers characterizing meltpool morphology of the laser powder bed fusion process

Microstructure evolution in metal additive manufacturing (AM) is a complex multi-physics and multi-scale problem. Understanding the impact of AM process conditions on the microstructure evolution and the resulting mechanical properties of the printed part is an active area of research. At the meltpool scale, the thermo-fluidic governing equations have been extensively modeled in the literature to understand the meltpool conditions and the thermal gradients in its vicinity. In many phenomena governed by partial differential equations, dimensional analysis and identification of important dimensionless numbers can provide significant insights into the process dynamics. In this context, a novel strategy using dimensional analysis and the method of linear least squares regression to numerically investigate the thermo-fluidic governing equations of the Laser Powder Bed Fusion AM process is presented in this work. First, the governing equations are solved using the Finite Element Method, and the model predictions are validated by comparing with experimentally estimated cooling rates, and with numerical results from the literature. Then, through dimensional analysis, an important dimensionless quantity - interpreted as a measure of heat absorbed by the powdered material and the meltpool, is identified. This dimensionless measure of heat absorbed, along with classical dimensionless quantities such as Peclet, Marangoni, and Stefan numbers, is used to investigate advective transport in the meltpool for different alloys. Further, the framework is used to study the variations of thermal gradients and the solidification cooling rate. Important correlations linking meltpool morphology and microstructure evolution related variables with classical dimensionless numbers are the key contribution of this work.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

DFT-FE -- A massively parallel adaptive finite-element code for large-scale density functional theory calculations

We present an accurate, efficient and massively parallel finite-element code, DFT-FE, for large-scale ab-initio calculations (reaching $\sim 100,000$ electrons) using Kohn-Sham density functional theory (DFT). DFT-FE is based on a local real-space variational formulation of the Kohn-Sham DFT energy functional that is discretized using a higher-order adaptive spectral finite-element (FE) basis, and treats pseudopotential and all-electron calculations in the same framework, while accommodating non-periodic, semi-periodic and periodic boundary conditions. We discuss the main aspects of the code, which include, the various strategies of adaptive FE basis generation, and the different approaches employed in the numerical implementation of the solution of the discrete Kohn-Sham problem that are focused on significantly reducing the floating point operations, communication costs and latency. We demonstrate the accuracy of DFT-FE by comparing the energies, ionic forces and periodic cell stresses on a wide range of problems with popularly used DFT codes. Further, we demonstrate that DFT-FE significantly outperforms widely used plane-wave codes---both in CPU-times and wall-times, and on both non-periodic and periodic systems---at systems sizes beyond a few thousand electrons, with over $5-10$ fold speedups in systems with more than 10,000 electrons. The benchmark studies also highlight the excellent parallel scalability of DFT-FE, with strong scaling demonstrated on up to 192,000 MPI tasks.

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

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

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