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

Kaushik Dayal

Kaushik Dayal 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-scicond-mat.mes-hallmath.APcond-mat.softphysics.comp-phCell BehaviorComputational Engineering, Finance, and Sciencecond-mat.stat-mechmath-phmath.DSmath.MPphysics.app-phphysics.atom-ph

Catalog footprint

What is connected

15works
13topics
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

15 published item(s)

preprint2023arXiv

Distal and non-symmetrical crack nucleation in delamination of plates via dimensionally-reduced peridynamics

Exploiting the framework of peridynamics, a dimensionally-reduced plate formulation is developed that allows for the through-thickness nucleation and growth of fracture surfaces, enabling the treatment of delamination in a lower-dimensional model. Delamination fracture nucleation and propagation are treated by choosing the kinematics to be composed of an absolutely continuous part and a zone where jumps in the displacements are allowed. This assumption allows the explicit derivation of the dimensionally-reduced elastic energy, which shows a hierarchy of terms characterising the stored energy in the plane element. An interpretation of the various terms of the reduced energy is shown by means of the most straightforward paradigm of bond-based peridynamics. A striking feature of the reduced energy is that, despite the small-displacement assumption, there is a coupling between the membrane and bending terms. Semi-analytical solutions for simplified settings are obtained through a minimization procedure, and a range of nonstandard behaviours such as \textit{distal} crack nucleation and curved crack path are captured by the model. Finally, convergence of the proposed peridynamic reduced model to a local elastic theory for vanishing nonlocal lengthscale is determined, giving a local cohesive model for fracture.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Comparison of simulated and measured grain volume changes during grain growth

The three-dimensional microstructure of Ni, observed after five annealing intervals, was compared to simulations of grain growth using the threshold dynamics method with the assumption of capillarity as the only driving force. A grain-by-grain comparison made it possible to identify the sources of differences between the simulation and experiment. The most significant difference was for grains of the smallest sizes, which the simulation predicted would lose volume and disappear at a greater rate than observed in the experiment. The loss of grains created errors in the numbers of neighbors of the remaining grains, and it was found that errors in the simulated grain volume were correlated to errors in the number of near neighbors. While anisotropic grain boundary properties likely play a role in the differences, the size dependence of the errors suggest that it might be necessary to include a size dependence in the model for grain boundary migration kinetics.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Multiband homogenization of metamaterials in real-space: Higher-order nonlocal models and scattering at external surfaces

This work develops a dynamic homogenization approach for metamaterials. It finds an approximate macroscopic homogenized equation with constant coefficients posed in space and time; however, the resulting homogenized equation is higher order in space and time. The homogenized equation can be used to solve initial-boundary-value problems posed on arbitrary non-periodic macroscale geometries with macroscopic heterogeneity, such as bodies composed of several different metamaterials or with external boundaries. First, considering a single band, the higher-order space derivatives lead to additional continuity conditions at the boundary between a homogeneous material and a metamaterial. These provide predictions of wave scattering in 1-d and 2-d that match well with the exact fine-scale solution; compared to alternative approaches, they provide a single equation that is valid over a broad range of frequencies, are easy to apply, and are much faster to compute. Next, the setting of two bands with a bandgap is considered. The homogenized equation has also higher-order time derivatives. Notably, the homogenized model provides a single equation that is valid over both bands and the bandgap. The continuity conditions are applied to wave scattering at a boundary, and show good agreement with the exact fine-scale solution. Using that the order of the highest time derivative is proportional to the number of bands considered, a nonlocal-in-time structure is conjectured for the homogenized equation in the limit of infinite bands. This suggests that homogenizing over finer length and time scales -- with the temporal homogenization being carried out through the consideration of higher bands in the dispersion relation -- is a mechanism for the emergence of macroscopic spatial and temporal nonlocality, with the extent of temporal nonlocality being related to the number of bands considered.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Phase-Field Finite Deformation Fracture with an Effective Energy for Regularized Crack Face Contact

Phase-field models are a leading approach for realistic fracture problems. They treat the crack as a second phase and use gradient terms to smear out the crack faces, enabling the use of standard numerical methods for simulations. This regularization causes cracks to occupy a finite volume in the reference, and leads to the inability to appropriately model the closing or contacting -- without healing -- of crack faces. Specifically, the classical idealized crack face tractions are that the shear component is zero, and that the normal component is zero when the crack opens and identical to the intact material when the crack closes. Phase-field fracture models do not replicate this behavior. This work addresses this shortcoming by introducing an effective crack energy density that endows the regularized (finite volume) phase-field crack with the effective properties of an idealized sharp crack. The approach is based on applying the QR (upper triangular) decomposition of the deformation gradient tensor in the basis of the crack, enabling a transparent identification of the crack deformation modes. By then relaxing over those modes that do not cost energy, an effective energy is obtained that has the intact response when the crack faces close and zero energy when the crack faces are open. A highlight of this approach is that it lies completely in the setting of finite deformation, enabling potential application to soft materials and other settings with large deformation or rotations. The model is applied to numerically study representative complex loadings, including (1) cyclic loading on a cavity in a soft solid that shows the growth and closing of cracks in complex stress states; and (2) cyclic shear that shows a complex pattern of crack branching driven by the closure of cracks.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Energetic Formulation of Large-Deformation Poroelasticity

The modeling of coupled fluid transport and deformation in a porous medium is essential to predict the various geomechanical process such as CO2 sequestration, hydraulic fracturing, and so on. Current applications of interest, for instance, that include fracturing or damage of the solid phase, require a nonlinear description of the large deformations that can occur. This paper presents a variational energy-based continuum mechanics framework to model large-deformation poroelasticity. The approach begins from the total free energy density that is additively composed of the free energy of the components. A variational procedure then provides the balance of momentum, fluid transport balance, and pressure relations. A numerical approach based on finite elements is applied to analyze the behavior of saturated and unsaturated porous media using a nonlinear constitutive model for the solid skeleton. Examples studied include the Terzaghi and Mandel problems; a gas-liquid phase-changing fluid; multiple immiscible gases; and unsaturated systems where we model injection of fluid into soil. The proposed variational approach can potentially have advantages for numerical methods as well as for combining with data-driven models in a Bayesian framework.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Effective Response of Heterogeneous Materials using the Recursive Projection Method

This paper applies the Recursive Projection Method (RPM) to the problem of finding the effective mechanical response of a periodic heterogeneous solid. Previous works apply the Fast Fourier Transform (FFT) in combination with various fixed-point methods to solve the problem on the periodic unit cell. These have proven extremely powerful in a range of problems ranging from image-based modeling to dislocation plasticity. However, the fixed-point iterations can converge very slowly, or not at all, if the elastic properties have high contrast, such as in the case of voids. The paper examines the reasons for slow, or lack of convergence, in terms of a variational perspective. In particular, when the material contains regions with zero or very small stiffness, there is lack of uniqueness, and the energy landscape has flat or shallow directions. Therefore, in this work, the fixed-point iteration is replaced by the RPM iteration. The RPM uses the fixed-point iteration to adaptively identify the subspace on which fixed-point iterations are unstable, and performs Newton iterations only on the unstable subspace, while fixed-point iterations are performed on the complementary stable subspace. This combination of efficient fixed-point iterations where possible, and expensive but well-convergent Newton iterations where required, is shown to lead to robust and efficient convergence of the method. In particular, RPM-FFT converges well for a wide range of choices of the reference medium, while usual fixed-point iterations are usually sensitive to this choice.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Statistical Mechanical Analysis of the Electromechanical Coupling in an Electrically-Responsive Polymer Chain

Polymeric materials that couple deformation and electrostatics have the potential for use in soft sensors and actuators with potential applications ranging from robotic, biomedical, energy, aerospace and automotive technologies. In contrast to the mechanics of polymers that has been studied using statistical mechanics approaches for decades, the coupled response under deformation and electrical field has largely been modeled only phenomenologically at the continuum scale. In this work, we examine the physics of the coupled deformation and electrical response of an electrically-responsive polymer chain using statistical mechanics. We begin with a simple anisotropic model for the electrostatic dipole response to electric field of a single monomer, and use a separation of energy scales between the electrostatic field energy and the induced dipole field energy to reduce the nonlocal and infinite-dimensional statistical averaging to a simpler local finite-dimensional averaging. In this simplified setting, we derive the equations of the most likely monomer orientation density using the maximum term approximation, and a chain free energy is derived using this approximation. These equations are investigated numerically and the results provide insight into the physics of electro-mechanically coupled elastomer chains. Closed-form approximations are also developed in the limit of small electrical energy with respect to thermal energy; in the limit of small mechanical tension force acting on the chain; and using asymptotic matching for general chain conditions.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Symmetry-Adapted Tight-Binding Electronic Structure Analysis of Carbon Nanotubes with Defects, Kinks, Twist, and Stretch

This paper applies a symmetry-adapted method to examine the influence of deformation and defects on the electronic structure and band structure in carbon nanotubes. First, the symmetry-adapted approach is used to develop the analog of Bloch waves. Building on this, the technique of perfectly-matched layers is applied to develop a method to truncate the computational domain of electronic structure calculations without spurious size effects. This provides an efficient and accurate numerical approach to compute the electronic structure and electromechanics of defects in nanotubes. The computational method is applied to study the effect of twist, stretch, and bending, with and without various types of defects, on the band structure of nanotubes. Specifically, the effect of stretch and twist on band structure in defect-free conducting and semiconducting nanotubes is examined, and the interaction with vacancy defects is elucidated. Next, the effect of localized bending or kinking on the electronic structure is studied. Finally, the paper examines the effect of 5-8-5 Stone-Wales defects. In all of these settings, the perfectly-matched layer method enables the calculation of localized non-propagating defect modes with energies in the bandgap of the defect-free nanotube.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Coherent motion of monolayer sheets under confinement and its pathological implications

Coherent angular rotation of epithelial cells is thought to contribute to many vital physiological processes including tissue morphogenesis and glandular formation. However, factors regulating this motion, and the implications of this motion if perturbed remain incompletely understood. In the current study, we address these questions using a cell-center based model in which cells are polarized, motile, and interact with the neighboring cells via harmonic forces. We demonstrate that, a simple evolution rule in which the polarization of any cell tends to orient with its velocity vector can induce coherent motion in geometrically confined environments. In addition to recapitulating coherent rotational motion observed in experiments, our results also show the presence of radial movements and tissue behavior that can vary between solid-like and fluid-like. We show that the pattern of coherent motion is dictated by the combination of different physical parameters including number density, cell motility, system size, bulk cell stiffness and stiffness of cell-cell adhesions. We further observe that, perturbations in the form of cell division can induce a reversal in the direction of motion when cell division occurs synchronously. Moreover, when the confinement is removed, we see that the existing coherent motion leads to cell scattering, with bulk cell stiffness and stiffness of cell-cell contacts dictating the invasion pattern. In summary, our study provides an in-depth understanding of the origin of coherent rotation in confined tissues, and extracts useful insights into the influence of various physical parameters on the pattern of such movements

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A Dynamic Multiscale Phase-field Model for Structural Transformations and Twinning: Regularized Interfaces with Transparent Prescription of Complex Kinetics and Nucleation

The motion of microstructural interfaces is important in modeling materials that undergo twinning and structural phase transformations. Continuum models fall into two classes: sharp-interface models, where interfaces are singular surfaces; and regularized-interface models, such as phase-field models, where interfaces are smeared out. The former are challenging for numerical solutions because the interfaces need to be explicitly tracked, but have the advantage that the kinetics of existing interfaces and the nucleation of new interfaces can be transparently and precisely prescribed. In contrast, phase-field models do not require explicit tracking of interfaces, thereby enabling relatively simple numerical calculations, but the specification of kinetics and nucleation is both restrictive and extremely opaque. This prevents straightforward calibration of phase-field models to experiment and/or molecular simulations, and breaks the multiscale hierarchy of passing information from atomic to continuum. We present the formulation of a phase-field model -- i.e., a model with regularized interfaces that do not require explicit numerical tracking -- that allows for easy and transparent prescription of complex interface kinetics and nucleation. The key ingredients are a re-parametrization of the energy density to clearly separate nucleation from kinetics; and an evolution law that comes from a conservation statement for interfaces. This enables clear prescription of nucleation through the source term of the conservation law and of kinetics through an interfacial velocity field. A formal limit of the kinetic driving force recovers the classical continuum sharp-interface driving force, providing confidence in both the re-parametrized energy and the evolution statement. We present a number of numerical calculations in one and two dimensions to characterize and demonstrate the formulation.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Anomalous phonon behavior of carbon nanotubes: First-order influence of external load

External loads typically have a indirect influence on phonon curves, i.e., they influence the phonon curves by changing the state about which linearization is performed. In this paper, we show that in nanotubes, the axial load has a direct first-order influence on the long-wavelength behavior of the transverse acoustic (TA) mode. In particular, when the tube is force-free the TA mode frequencies vary quadratically with wave number and have curvature (second derivative) proportional to the square-root of the nanotube's bending stiffness. When the tube has non-zero external force, the TA mode frequencies vary linearly with wave number and have slope proportional to the square-root of the axial force. Therefore, the TA phonon curves -- and associated transport properties -- are not material properties but rather can be directly tuned by external loads. In addition, we show that the out-of-plane shear deformation does {\em not} contribute to this mode and the unusual properties of the TA mode are exclusively due to bending. Our calculations consist of 3 parts: first, we use a linear chain of atoms as an illustrative example that can be solved in close-form; second, we use our recently-developed symmetry-adapted phonon analysis method to present direct numerical evidence; and finally, we present a simple mechanical model that captures the essential physics of the geometric nonlinearity in slender nanotubes that couples the axial load directly to the phonon curves. We also compute the Density of States and show the significant effect of the external load.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Atomistic-to-Continuum Multiscale Modeling with Long-Range Electrostatic Interactions in Ionic Solids

We present a multiscale atomistic-to-continuum method for ionic crystals with defects. Defects often play a central role in ionic and electronic solids, not only to limit reliability, but more importantly to enable the functionalities that make these materials of critical importance. Examples include solid electrolytes that conduct current through the motion of charged point defects, and complex oxide ferroelectrics that display multifunctionality through the motion of domain wall defects. Therefore, it is important to understand the structure of defects and their response to electrical and mechanical fields. A central hurdle, however, is that interactions in ionic solids include both short-range atomic interactions as well as long-range electrostatic interactions. Existing atomistic-to-continuum multi-scale methods, such as the Quasicontinuum method, are applicable only when the atomic interactions are short-range. In addition, empirical reductions of quantum mechanics to density functional models are unable to capture key phenomena of interest in these materials. To address this open problem, we develop a multiscale atomistic method to coarse-grain the long-range electrical interactions in ionic crystals with defects. In these settings, the charge density is rapidly varying, but in an almost-periodic manner. The key idea is to use the polarization density field as a multiscale mediator that enables efficient coarse-graining by exploiting the almost-periodic nature of the variation. In regions far from the defect, where the crystal is close-to-perfect, the polarization field serves as a proxy that enables us to avoid accounting for the details of the charge variation. We combine this approach for long-range electrostatics with the standard Quasicontinuum method for short-range interactions to achieve an efficient multiscale atomistic-to-continuum method.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A fundamental improvement to Ericksen-Leslie kinematics

We demonstrate theory and computations for finite-energy line defect solutions in an improvement of Ericksen-Leslie liquid crystal theory. Planar director fields are considered in two and three space dimensions, and we demonstrate straight as well as loop disclination solutions. The possibility of static balance of forces in the presence of a disclination and in the absence of flow and body forces is discussed. The work exploits an implicit conceptual connection between the Weingarten-Volterra characterization of possible jumps in certain potential fields and the Stokes-Helmholtz resolution of vector fields. The theoretical basis of our work is compared and contrasted with the theory of Volterra disclinations in elasticity. Physical reasoning precluding a gauge-invariant structure for the model is also presented.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Coarse Variables of Autonomous ODE Systems and Their Evolution

Given an autonomous system of ordinary differential equations (ODE), we consider developing practical models for the deterministic, slow/coarse behavior of the ODE system. Two types of coarse variables are considered. The first type consists of running finite time averages of phase functions. Approaches to construct the coarse evolution equation for this type are discussed and implemented on a 'Forced' Lorenz system and a singularly perturbed system whose fast flow does not necessarily converge to an equilibrium. We explore two strategies. In one, we compute (locally) invariant manifolds of the fast dynamics, parameterized by the slow variables. In the other, the choice of our coarse variables automatically guarantees them to be 'slow' in a precise sense. This allows their evolution to be phrased in terms of averaging utilizing limit measures (probability distributions) of the fast flow. Coarse evolution equations are constructed based on these approaches and tested against coarse response of the 'microscopic' models. The second type of coarse variables are defined as (non-trivial) scalar state functions that are required by design to evolve autonomously, to the extent possible, with the goal of being candidate state functions for unambiguously initializable coarse dynamics. The question motivates a mathematical restatement in terms of a first-order PDE. A computational approximation is developed and tested on the Lorenz system and the Hald Hamiltonian system.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Symmetry-Adapted Phonon Analysis of Nanotubes

The characteristics of phonons, i.e. linearized normal modes of vibration, provide important insights into many aspects of crystals, e.g. stability and thermodynamics. In this paper, we use the Objective Structures framework to make concrete analogies between crystalline phonons and normal modes of vibration in non-crystalline but highly symmetric nanostructures. Our strategy is to use an intermediate linear transformation from real-space to an intermediate space in which the Hessian matrix of second derivatives is block-circulant. The block-circulant nature of the Hessian enables us to then follow the procedure to obtain phonons in crystals: namely, we use the Discrete Fourier Transform from this intermediate space to obtain a block-diagonal matrix that is readily diagonalizable. We formulate this for general Objective Structures and then apply it to study carbon nanotubes of various chiralities that are subjected to axial elongation and torsional deformation. We compare the phonon spectra computed in the Objective Framework with spectra computed for armchair and zigzag nanotubes. We also demonstrate the approach by computing the Density of States. In addition to the computational efficiency afforded by Objective Structures in providing the transformations to almost-diagonalize the Hessian, the framework provides an important conceptual simplification to interpret the phonon curves.

0 citations0 reviewsNo code linkedOpen access