Source author record

Giacomo Po

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

14works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

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

14 published item(s)

preprint2026arXiv

Active Learning of A Crystal Plasticity Flow Rule From Discrete Dislocation Dynamics Simulations

Continuum-scale material deformation models, such as crystal plasticity, can significantly enhance their predictive accuracy by incorporating input from lower-scale (i.e., mesoscale) models. The procedure to generate and extract the relevant information is however typically complex and ad hoc, involving decision and intervention by domain experts, leading to long development times. In this study, we develop a principled approach for calibration of continuum-scale models using lower scale information by representing a crystal plasticity flow rule as a Gaussian process model. This representation allows for efficient parameter space exploration, guided by the uncertainty embedded in the model through a process known as Bayesian optimization. We demonstrate a semi-autonomous Bayesian optimization loop which instantiates discrete dislocation dynamics simulations whose initial conditions are automatically chosen to optimize the uncertainty of a model crystal plasticity flow rule. Our self-guided computational pipeline efficiently generated a dataset and corresponding model whose error, uncertainty, and physical feature sensitivities were validated with comparison to an independent dataset four times larger, demonstrating a valuable and efficient active learning implementation readily transferable to similar material systems.

preprint2023arXiv

Plasticity in irradiated FeCrAl nanopillars investigated using discrete dislocation dynamics

In this paper, we investigate plasticity in irradiated FeCrAl nanopillars using discrete dislocation dynamics simulations (DDD), with comparisons to transmission electron microscopic (TEM) in situ tensile tests of ion and neutron irradiated commercial FeCrAl alloy C35M. The effects of irradiation-induced defects, such as a/2 111 and a 100 type loops and composition fluctuations representative of phase separation in irradiated FeCrAl alloys, are investigated separately as well as superposed together in simulations. We explore the effects of defects on the stress-strain behavior, specifically yield strength and hardening response, of FeCrAl nanopillars. Our simulations confirm the widely accepted fact that irradiated alloys exhibit a stress-strain response with higher yield strength and hardening as compared to homogeneous alloys. However, our DDD calculations reveal an atypical superposition of the hardening contributions due to composition inhomogeneity and irradiation loops wherein hardening due composition inhomogeneity counteracts hardening due to irradiation loops at small scales. As a result, we observe that the yield strength in irradiated alloys, after taking into consideration the effects of both composition inhomogeneity and irradiation loops, is smaller than the yield strength of the alloys with only irradiation loops and is approximately same for the alloy with composition inhomogeneity alone. We identify this destructive interference in the superposition in our parallel TEM in situ tensile tests on unirradiated, ion irradiated, and neutron irradiated C35M FeCrAl alloy as well. This destructive interference in the hardening contributions contrasts with the widely utilized dispersed barrier hardening (DBH) models by the experimental community to model the hardening contributions due to different irradiation induced defects.

preprint2021arXiv

Discrete Dislocation Dynamics for Crystal RVEs. Part 1: Periodic Loop Kinematics

A novel implementation of the dislocation flux boundary condition in discrete dislocation dynamics is presented. The continuity of the individual dislocation loops in a periodic representative crystal volume (RVE) is enforced across the boundary of the RVE with the help of a dual topological description for representing dislocation line kinematics in two equivalent spaces representing the deforming crystal, the RVE and the unbounded crystal spaces. The former describes the motion of the dislocations in the simulated crystal RVE whereas the latter represents the motion of dislocations in an infinite space containing all replicas of the RVE. A mapping between the two spaces forms the basis of the implementation of flux boundary condition. The implementation details are discussed in the context of statistical homogeneity of bulk crystals undergoing macroscopically homogeneous plastic deformation. Some test cases are presented and discussed to establish the proposed approach. It has been observed that the boundary nodes associated with dislocation segments bear no relevance in the motion of the dislocation segments from the standpoint of bulk crystal plasticity modeling throughout the deformation history of the crystal.

preprint2020arXiv

Inference and Prediction of Nanoindentation Response in FCC Crystals: Methods and Discrete Dislocation Simulation Examples

Nanoindentation, a common technique for probing the mechanical properties of crystalline materials, exhibits both surface and bulk-dominated responses that are linked to the parent crystal's elasticity and plasticity. For FCC crystals, the nanoindentation response is primarily controlled by the indented-grain's crystalline orientation and dislocation density. However, the indentation properties, such as hardness, have ambiguous significance at the nanoscale due to tip-specific size effects. Through extensive discrete dislocation dynamics simulations of flat-punch nanoindentation on FCC single crystals, we show that a one-to-one correspondence between an indented location and the pre-existing dislocation density is possible with the help of statistical, machine-learning algorithms. We cluster and classify various experimentally plausible ensembles with varying dislocation density and/or crystalline orientation, using nanoindentation force-depth curves or post-indent surface displacement images. We use this classification to statistically predict the nanoindentation response of a given location, through the average response of classified load-displacement curves. In addition, we demonstrate that our simulation results qualitatively agree with experimental data on common FCC single crystals. Ultimately, we propose a novel pathway to uncover the depth dependence of local surface material properties for more complete, non-destructive, microstructural characterization, through nanoindentation.

preprint2020arXiv

Plasticity without phenomenology: a first step

A novel, concurrent multiscale approach to meso/macroscale plasticity is demonstrated. It utilizes a carefully designed coupling of a partial differential equation (pde) based theory of dislocation mediated crystal plasticity with time-averaged inputs from microscopic Dislocation Dynamics (DD), adapting a state-of-the-art mathematical coarse-graining scheme. The stress-strain response of mesoscopic samples at realistic, slow, loading rates up to appreciable values of strain is obtained, with significant speed-up in compute time compared to conventional DD. Effects of crystal orientation, loading rate, and the ratio of the initial mobile to sessile dislocation density on the macroscopic response, for both load and displacement controlled simulations are demonstrated. These results are obtained without using any phenomenological constitutive assumption, except for thermal activation which is not a part of microscopic DD. The results also demonstrate the effect of the internal stresses on the collective behavior of dislocations, manifesting, in a set of examples, as a Stage I to Stage II hardening transition.

preprint2019arXiv

The $Λ$-invariant and topological pathways to influence sub-micron strength and crystal plasticity

In small volumes, sample dimensions are known to strongly influence mechanical behavior, especially strength and crystal plasticity. This correlation fades away at the so-called mesoscale, loosely defined at several micrometers in both experiments and simulations. However, this picture depends on the entanglement of the initial defect configuration. In this paper, we study the effect of sample dimensions with a full control on dislocation topology, through the use of a novel observable for dislocation ensembles, the $Λ$-invariant, that depends only on mutual dislocation linking: It is built on the natural vortex character of dislocations and it has a continuum/discrete correspondence that may assist multiscale modeling descriptions. We investigate arbitrarily complex initial dislocation microstructures in sub-micron-sized pillars, using three-dimensional discrete dislocation dynamics simulations for finite volumes. We demonstrate how to engineer nanoscale dislocation ensembles that appear virtually independent from sample dimensions, either by biased-random dislocation loop deposition or by sequential mechanical loads of compression and torsion.

preprint2016arXiv

A phenomenological dislocation mobility law for bcc metals

Dislocation motion in body centered cubic (bcc) metals displays a number of specific features that result in a strong temperature dependence of the flow stress, and in shear deformation asymmetries relative to the loading direction as well as crystal orientation. Here we develop a generalized dislocation mobility law in bcc metals, and demonstrate its use in discrete Dislocation Dynamics (DD) simulations of plastic flow in tungsten (W) micro pillars. We present the theoretical background for dislocation mobility as a motivating basis for the developed law. Analytical theory, molecular dynamics (MD) simulations, and experimental data are used to construct a general phenomenological description. The usefulness of the mobility law is demonstrated through its application to modeling the plastic deformation of W micro pillars. The model is consistent with experimental observations of temperature and orientation dependence of the flow stress and the corresponding dislocation microstructure.

preprint2016arXiv

Comparison of dislocation density tensor fields derived from discrete dislocation dynamics and crystal plasticity simulations of torsion

The importance of accurate simulation of the plastic deformation of ductile metals to the design of structures and components is well-known. Many techniques exist that address the length scales relevant to deformation pro- cesses, including dislocation dynamics (DD), which models the interaction and evolution of discrete dislocation line segments, and crystal plasticity (CP), which incorporates the crystalline nature and restricted motion of dis- locations into a higher scale continuous field framework. While these two methods are conceptually related, there have been only nominal efforts focused on the system-level material response that use DD-generated information to enhance the fidelity of plasticity models. To ascertain to what degree the predictions of CP are consistent with those of DD, we compare their global and microstructural response in a number of deformation modes. After using nominally homogeneous compression and shear deformation dislocation dynamics simulations to calibrate crystal plasticity flow rule parameters, we compare not only the system-level stress-strain response of prismatic wires in torsion but also the resulting geometrically necessary dislocation density fields. To establish a connection between explicit description of dislocations and the continuum assumed with crystal plasticity simulations, we ascertain the minimum length-scale at which meaningful dislocation density fields appear. Our results show that, for the case of torsion, the two material models can produce comparable spatial dislocation density distributions.

preprint2016arXiv

Controlling Strain Bursts and Avalanches at the Nano-to-Micro Scale

We demonstrate, through 3-dimensional discrete dislocation dynamics simulations, that the com- plex dynamical response of nano and micro crystals to external constraints can be tuned. Under load rate control, strain bursts are shown to exhibit scale-free avalanche statistics, similar to critical phenomena in many physical systems. For the other extreme of displacement rate control, strain burst response transitions to quasi-periodic oscillations, similar to stick-slip earthquakes. External load mode control is shown to enable a qualitative transition in the complex collective dynamics of dislocations from self-organized criticality to quasi-periodic oscillations.

preprint2016arXiv

Formation and interaction mechanisms of dipolar dislocation loops in fcc metals

Dipolar dislocation loops, prevalent in fcc metals, are widely recognized as controlling many physical aspects of plastic deformation. We present results of 3D dislocation dynamics simulations that shed light on the mechanisms of their formation, motion, interactions, and large-scale patterning. We identify two main formation mechanisms, enabled by cross-slip, and show that arrays of dipoles can be easily formed as a result of the interaction between glide screw dislocations. We present a systematic analysis of the spectrum of possible junctions that can form as a result of mutual interaction between dipoles, and between dipoles and glide dislocations. We show that fully immobile dislocation segments arise in particular cases of these interactions, leading to hardening and Frank-Read type sources. We reveal that the collective motion of dipolar loop arrays can be induced by glide dislocations in the channels of Persistent Slip Bands (PSB), and result in their clustering within PSB channel walls. An efficient tripolar drag mechanism is found to contribute to the clustering of dipolar loops near channel walls.

preprint2016arXiv

Microstructural Comparison of the Kinematics of Discrete and Continuum Dislocations Models

The Continuum Dislocation Dynamics (CDD) theory and the Discrete Dislocation Dynamics (DDD) method are compared based on concise mathematical formulations of the coarse graining of discrete data. A numerical tool for converting from a discrete to a continuum representation of a given dislocation configuration is developed, which allows to directly compare both simulation approaches based on continuum quantities (e.g. scalar density, geometrically necessary densities, mean curvature). Investigating the evolution of selected dislocation configurations within analytically given velocity fields for both DDD and CDD reveals that CDD contains a surprising number of important microstructural details.

preprint2016arXiv

The atomistic representation of first strain-gradient elastic tensors

We derive the atomistic representations of the elastic tensors appearing in the linearized theory of first strain-gradient elasticity for an arbitrary multi-lattice. In addition to the classical (2nd-Piola) stress and elastic moduli tensors, these include the rank-three double-stress tensor, the rank-five tensor of mixed elastic moduli, and the rank-six tensor of strain-gradient elastic moduli. The atomistic representations are closed-form analytical expressions in terms of the first and second derivatives of the interatomic potential with respect to interatomic distances, and dyadic products of relative atomic positions. Moreover, all expressions are local, in the sense that they depend only on the atomic neighborhood of a lattice site. Our results emanate from the condition of energetic equivalence between continuum and atomistic representations of a crystal, when the kinematics of the latter is governed by the Cauchy-Born rule. Using the derived expressions, we prove that the odd-order tensors vanish if the lattice basis admits central-symmetry. The analytical expressions are implemented as a KIM compliant algorithm to compute the strain gradient elastic tensors for various materials. Numerical results are presented to compare representative interatomic potentials used in the literature for cubic crystals, including simple lattices and multi-lattices. We observe that central potentials exhibit generalized Cauchy relations for the rank-six tensor of strain-gradient elastic moduli. In addition, this tensor is found to be indefinite for many potentials. We discuss the relationship between indefiniteness and material stability. Finally, the atomistic representations are specialized to central potentials in simple lattices. These expressions are used with analytical potentials to study the sensitivity of the elastic tensors to the choice of the cutoff radius.

preprint2015arXiv

Atomistically enabled nonsingular anisotropic elastic representation of near-core dislocation stress fields in $α$-iron

The stress fields of dislocations predicted by classical elasticity are known to be unrealistically large approaching the dislocation core, due to the singular nature of the theory. While in many cases this is remedied with the approximation of an effective core radius, inside which ad hoc regularizations are implemented, such approximations lead to a compromise in the accuracy of the calculations. In this work, an anisotropic non-singular elastic representation of dislocation fields is developed to accurately represent the near-core stresses of dislocations in $α$-iron. The regularized stress field is enabled through the use of a non-singular Green's tensor function of Helmholtz-type gradient anisotropic elasticity, which requires only a single characteristic length parameter in addition to the material's elastic constants. Using a novel magnetic bond-order potential to model atomic interactions in iron, molecular statics calculations are performed, and an optimization procedure is developed to extract the required length parameter. Results show the method can accurately replicate the magnitude and decay of the near-core dislocation stresses even for atoms belonging to the core itself. Comparisons with the singular isotropic and anisotropic theories show the non-singular anisotropic theory leads to a substantially more accurate representation of the stresses of both screw and edge dislocations near the core, in some cases showing improvements in accuracy of up to an order of magnitude. The spatial extent of the region in which the singular and non-singular stress differ substantially is also discussed. The general procedure we describe may in principle be applied to accurately model the near-core dislocation stresses of any arbitrarily shaped dislocation in anisotropic cubic media.

preprint2013arXiv

The solid angle and the Burgers formula in the theory of gradient elasticity: line integral representation

A representation of the solid angle and the Burgers formula as line integral is derived in the framework of the theory of gradient elasticity of Helmholtz type. The gradient version of the Eshelby-deWit representation of the Burgers formula of a closed dislocation loop is given. Such a form is suitable for the numerical implementation in 3D dislocation dynamics (DD).