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

Ankit Srivastava

Ankit Srivastava 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.MPcond-mat.softgr-qcmath.NANeural and Evolutionary ComputingNumerical Analysisphysics.app-phphysics.class-phphysics.comp-ph

Catalog footprint

What is connected

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

16 published item(s)

preprint2022arXiv

Scattering of mechanical waves from the perspective of open systems

In this paper, we consider the problem of mechanical wave scattering from a spatially finite system into an infinite surrounding environment. The goal is to illuminate why the scattering spectrum undergoes peaks and dips (resonances) at specific locations and how these locations connect to the vibrational properties of the scatterer. The resonance locations are connected to the eigenvalues of a finite dimensional effective operator, $H_{eff}$, corresponding to the scatterer. The developments are presented from the perspective of open systems, which seeks to convert the infinite dimensional scattering problem (scatterer+environment) into a finite dimensional effective problem involving only the finite scatterer. This is achieved through a projection operator formalism which allows us to formally calculate $H_{eff}$. An interesting corollary of our analysis is the deep connection between resonance locations in the scattering spectrum and the eigenfrequencies of the scatterer under Neumann boundary condition. We bring out this point further by considering 3D scattering from an elastic shell, connecting our results to classical results in acousto-elastic scattering theory.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Exact imposition of boundary conditions with distance functions in physics-informed deep neural networks

In this paper, we introduce a new approach based on distance fields to exactly impose boundary conditions in physics-informed deep neural networks. The challenges in satisfying Dirichlet boundary conditions in meshfree and particle methods are well-known. This issue is also pertinent in the development of physics informed neural networks (PINN) for the solution of partial differential equations. We introduce geometry-aware trial functions in artifical neural networks to improve the training in deep learning for partial differential equations. To this end, we use concepts from constructive solid geometry (R-functions) and generalized barycentric coordinates (mean value potential fields) to construct $ϕ$, an approximate distance function to the boundary of a domain. To exactly impose homogeneous Dirichlet boundary conditions, the trial function is taken as $ϕ$ multiplied by the PINN approximation, and its generalization via transfinite interpolation is used to a priori satisfy inhomogeneous Dirichlet (essential), Neumann (natural), and Robin boundary conditions on complex geometries. In doing so, we eliminate modeling error associated with the satisfaction of boundary conditions in a collocation method and ensure that kinematic admissibility is met pointwise in a Ritz method. We present numerical solutions for linear and nonlinear boundary-value problems over domains with affine and curved boundaries. Benchmark problems in 1D for linear elasticity, advection-diffusion, and beam bending; and in 2D for the Poisson equation, biharmonic equation, and the nonlinear Eikonal equation are considered. The approach extends to higher dimensions, and we showcase its use by solving a Poisson problem with homogeneous Dirichlet boundary conditions over the 4D hypercube. This study provides a pathway for meshfree analysis to be conducted on the exact geometry without domain discretization.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Scattering under Linear Non Self-Adjoint Operators: Case of in-Plane Elastic Waves

In this paper, we consider the problem of the scattering of in-plane waves at an interface between a homogeneous medium and a metamaterial. The relevant eigenmodes in the two regions are calculated by solving a recently described non self-adjoint eigenvalue problem particularly suited to scattering studies. The method efficiently produces all propagating and evanescent modes consistent with the application of Snell's law and is applicable to very general scattering problems. In a model composite, we elucidate the emergence of a rich spectrum of eigenvalue degeneracies. These degeneracies appear in both the complex and real domains of the wave-vector. However, since this problem is non self-adjoint, these degeneracies generally represent a coalescing of both the eigenvalues and eigenvectors (exceptional points). Through explicit calculations of Poynting vector, we point out an intriguing phenomenon: there always appears to be an abrupt change in the sign of the refraction angle of the wave on two sides of an exceptional point. Furthermore, the presence of these degeneracies, in some cases, hints at fast changes in the scattered field as the incident angle is changed by small amounts. We calculate these scattered fields through a novel application of the Betti-Rayleigh reciprocity theorem. We present several numerical examples showing a rich scattering spectrum. In one particularly intriguing example, we point out wave behavior which may be related to the phenomenon of resonance trapping. We also show that there exists a deep connection between energy flux conservation and the biorthogonality relationship of the non self-adjoint problem. The proof applies to the general class of scattering problems involving elastic waves (under self-adjoint or non self-adjoint operators).

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Negative Refraction, Beam Steering, Mode Switching, and High-pass Filtering in a 1-D Periodic Laminate

In this paper we show that a 1-D phononic crystal (laminate) can exhibit metamaterial wave phenomenon which is traditionally associated with 2-, and 3-D crystals. Moreover, due to the absence of a length scale in 2 of its dimensions, it can outperform higher dimensional crystals on some measures. This includes allowing only negative refraction over large frequency ranges and serving as a near-omnidirectional high-pass filter up to a large frequency value. First we provide a theoretical discussion on the salient characteristics of the dispersion relation of a laminate and formulate the solution of an interface problem by the application of the normal mode decomposition technique. We present a methodology with which to induce a pure negative refraction in the laminate. As a corollary to our approach of negative refraction, we show how the laminate can be used to steer beams over large angles for small changes in the incident angles (beam steering). Furthermore, we clarify how the transmitted modes in the laminate can be switched on and off by varying the angle of the incident wave by a small amount. Finally, we show that the laminate can be used as a remarkably efficient high-pass frequency filter. An appropriately designed laminate will reflect all plane waves from quasi-static to a large frequency, incident at it from all angles except for a small set of near-normal incidences. This will be true even if the homogeneous medium is impedance matched with the laminate. Due to the similarities between SH waves and EM waves it is expected that some or all of these results may also apply to EM waves in a layered periodic dielectric.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Variational Methods For Phononic Calculations

Three fundamental variational principles used for solving elastodynamic eigenvalue problems are studied within the context of elastic wave propagation in periodic composites (phononics). We study the convergence of the eigenvalue problems resulting from the displacement Rayleigh quotient, the stress Rayleigh quotient and the mixed quotient. The convergence rates of the three quotients are found to be related to the continuity and differentiability of the density and compliance variation over the unit cell. In general, the mixed quotient converges faster than both the displacement Rayleigh and the stress Rayleigh quotients, however, there exist special cases where either the displacement Rayleigh or the stress Rayleigh quotient shows the exact same convergence as the mixed-method. We show that all methods converge faster for smoother material property variations, but when density variation is rough, the difference between the mixed quotient and stress Rayleigh quotient is higher and similarly, when compliance variation is rough, the difference between the mixed quotient and displacement Rayleigh quotient is higher. Since eigenvalue problems such as those considered in this paper tend to be highly computationally intensive, it is expected that these results will lead to fast and efficient algorithms in the areas of phononics and photonics.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Causality and Passivity in Elastodynamics

What are the constraints placed on the frequency dependent constitutive tensors of elastodynamics by the requirements that the linear elastodynamic system under consideration be both causal (effects succeed causes) and passive (system doesn't produce energy)? Compared to electromagnetism, elastodynamics is complicated by its generally non diagonalizable constitutive tensors. In this paper we clarify the constraints that causality and passivity place on very general forms of elastodynamic constitutive relations. Specifically we show that the satisfaction of passivity (and causality) directly requires that the hermitian parts, as defined later, of the Fourier transforms of the time derivatives of the elastodynamic constitutive tensors be positive semi-definite at all frequencies. Additionally, we show that the conditions subsequently require that the non-hermitian parts of the Fourier transforms of the constitutive tensors, when divided by the imaginary number $i$, be positive semi-definite for positive values of frequency and negative semi-definite for negative values of frequency. We show that when major symmetries are assumed these definiteness relations apply simply to the real and imaginary parts of the relevant tensors. For diagonal and one-dimensional problems, these positive semi-definiteness relationships reduce to simple inequality relations. Finally we extend the results to highly general constitutive relations which include the Willis inhomogeneous relations as a special case.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Elastic Metamaterials and Dynamic Homogenization: A Review

In this paper we review the recent advances which have taken place in the understanding and applications of acoustic/elastic metamaterials. Metamaterials are artificially created composite materials which exhibit unusual properties which are not found in nature. We begin with presenting arguments from discrete systems which support the case for the existence of unusual material properties such as tensorial and/or negative density. The arguments are then extended to elastic continuums through coherent averaging principles. The resulting coupled and nonlocal homogenized relations, called the Willis relations, are presented as the natural description of inhomogeneous elastodynamics. They are specialized to Bloch waves propagating in periodic composites and we show that the Willis properties display the unusual behavior which is often required in metamaterial applications such as the Veselago lens. We finally present the recent advances in the area of transformation elastodynamics, charting its inspirations from transformation optics, clarifying its particular challenges, and identifying its connection with the constitutive relations of the Willis and the Cosserat types.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Mixed-variational formulation for phononic band-structure calculation of arbitrary unit cells

This paper presents phononic band-structure calculation results obtained using a mixed variational formulation for 1-, and 2-dimensional unit cells. The formulation itself is presented in a form which is equally applicable to 3-dimensiomal cases. It has been established that the mixed-variational formulation presented in this paper shows faster convergence with considerably greater accuracy than variational principles based purely on the displacement field, especially for problems involving unit cells with discontinuous constituent properties. However, the application of this formulation has been limited to fairly simple unit cells. In this paper we have extended the scope of the formulation by employing numerical integration techniques making it applicable for the evaluation of the phononic band-structure of unit cells displaying arbitrary complexity in their Bravais structure and in the shape, size, number, and anisotropicity of their micro-constituents. The approach is demonstrated through specific examples

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

On the Limit and Applicability of Dynamic Homogenization

Recent years have seen considerable research success in the field of dynamic homogenization which seeks to define frequency dependent effective properties for heterogeneous composites for the purpose of studying wave propagation. There is an approximation involved in replacing a heterogeneous composite with its homogenized equivalent. In this paper we propose a quantification to this approximation. We study the problem of reflection at the interface of a layered periodic composite and its dynamic homogenized equivalent. It is shown that if the homogenized parameters are to appropriately represent the layered composite in a finite setting and at a given frequency, then reflection at this special interface must be close to zero at that frequency. We show that a comprehensive homogenization scheme proposed in an earlier paper results in negligible reflection in the low frequency regime, thereby suggesting its applicability in a finite composite setting. In this paper we explicitly study a 2-phase composite and a 3-phase composite which exhibits negative effective properties over its second branch. We show that based upon the reflected energy profile of the two cases, there exist good arguments for considering the second branch of a 3-phase composite a true negative branch with negative group velocity. The results open intriguing questions regarding the effects of replacing a semi-infinite periodic composite with its Bloch-wave (infinite domain) dynamic properties on such phenomenon as negative refraction.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Correlating toughness and roughness in ductile fracture

Three dimensional calculations of ductile crack growth under mode I plane strain, small scale yielding conditions are carried out using an elastic-viscoplastic constitutive relation for a progres- sively cavitating plastic solid with two populations of void nucleating second phase particles. Full field solutions are obtained for three dimensional material microstructures characterized by ran- dom distributions of void nucleating particles. Crack growth resistance curves and fracture surface roughness statistics are calculated using standard procedures. The range of void nucleating particle volume fractions considered give rise to values of toughness, JIC, that vary by a factor of four. For all volume fractions considered, the computed fracture surfaces are self-affine over a size range of about two orders of magnitude with a roughness exponent of 0.54 $\pm$ 0.03. For small void nucleating particle volume fractions, the mean large particle spacing serves as a single dominant length scale. In this regime, the correlation length of the fracture surface corresponding to the cut-off of the self-affine behavior is found to be linearly related to JIC thus quantitatively correlating toughness and fracture surface roughness.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Fast distributed phononic band-structure calculations through a GPU accelerated mixed-variational formulation

In this paper we present a Graphical Processing Unit accelerated mixed variational formulation for fast phononic band-structure calculation of arbitrarily complex unit cells and report speed gains of a hundred fold over unoptimized serial cpu computations. To the author's knowledge this is the first application of gpu computing to a non-FE/FDTD bandstructure algorithm. The formulation is presented in a form which is applicable to 1-, 2-, and 3-D cases. However, in this paper we concentrate on optimizing the formulation within the paradigm of gpu computing, presenting results for 2-D unit cells. We describe the application of the formulation with a long term view towards highly efficient and massively distributed band-structure algorithms suitable for tackling optimization and inverse problems. We report that the accurate band-structure evaluation over the boundary of the Irreducible Brillouin Zone (IBZ) for the first 18 phononic branches of a complex 2-D unit cell (with 1132 different phases) can be achieved in less than 20 seconds on a regular desktop. For a simpler unit cell, the first few phononic branches are calculated in less than 2.5 seconds on the same system. The scheme presented in this paper, therefore, represents considerable savings in the computational costs of phononic band-structures.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

From Newton's Universal Gravitation to Einstein's Geometric Theory of Gravity

Starting with Newton's law of universal gravitation, we generalize it step-by-step to obtain Einstein's geometric theory of gravity. Newton's gravitational potential satisfies the Poisson equation. We relate the potential to a component of the metric tensor by equating the nonrelativistic result of the principle of stationary proper time to the Lagrangian for a classical gravitational field. In the Poisson equation the Laplacian of the component of the metric tensor is generalized to a cyclic linear combination of the second derivatives of the metric tensor. In local coordinates it is a single component of the Ricci tensor. This component of the Ricci tensor is proportional to the mass density that is related to a single component of the energy-momentum stress tensor and its trace. We thus obtain a single component of Einstein's gravitational field equation in local coordinates. From the principle of general covariance applied to a single component, we obtain all tensor components of Einstein's gravitational field equation.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Bounds on Effective Dynamic Properties of Elastic Composites

We present general, computable, improvable, and rigorous bounds for the total energy of a finite heterogeneous volume element or a periodically distributed unit cell of an elastic composite of any known distribution of inhomogeneities of any geometry and elasticity, undergoing a harmonic motion at a fixed frequency or supporting a single-frequency Bloch-form elastic wave of a given wave-vector. These bounds are rigorously valid for \emph{any consistent boundary conditions} that produce in the finite sample or in the unit cell, either a common average strain or a common average momentum. No other restrictions are imposed. We do not assume statistical homogeneity or isotropy. Our approach is based on the Hashin-Shtrikman (1962) bounds in elastostatics, which have been shown to provide strict bounds for the overall elastic moduli commonly defined (or actually measured) using uniform boundary tractions and/or linear boundary displacements; i.e., boundary data corresponding to the overall uniform stress and/or uniform strain conditions. Here we present strict bounds for the dynamic frequency-dependent constitutive parameters of the composite and give explicit expressions for a direct calculation of these bounds.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Overall Dynamic Constitutive Relations of Micro-structured Elastic Composites

A method for homogenization of a heterogeneous (finite or periodic) elastic composite is presented. It allows direct, consistent, and accurate evaluation of the averaged overall frequency-dependent dynamic material constitutive relations. It is shown that when the spatial variation of the field variables is restricted by a Bloch-form (Floquet-form) periodicity, then these relations together with the overall conservation and kinematical equations accurately yield the displacement or stress modeshapes and, necessarily, the dispersion relations. It also gives as a matter of course point-wise solution of the elasto-dynamic field equations, to any desired degree of accuracy. The resulting overall dynamic constitutive relations however, are general and need not be restricted by the Bloch-form periodicity. The formulation is based on micro-mechanical modeling of a representative unit cell of the composite proposed by Nemat-Nasser and coworkers; see, e.g., [1] and [2].

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Overall Dynamic Properties of 3-D periodic elastic composites

A method for the homogenization of 3-D periodic elastic composites is presented. It allows for the evaluation of the averaged overall frequency dependent dynamic material constitutive tensors relating the averaged dynamic field variable tensors of velocity, strain, stress, and linear momentum. The formulation is based on micromechanical modeling of a representative unit cell of a composite proposed by Nemat-Nasser & Hori (1993), Nemat-Nasser et. al. (1982) and Mura (1987) and is the 3-D generalization of the 1-D elastodynamic homogenization scheme presented by Nemat-Nasser & Srivastava (2011). We show that for 3-D periodic composites the overall compliance (stiffness) tensor is hermitian, irrespective of whether the corresponding unit cell is geometrically or materially symmetric.Overall mass density is shown to be a tensor and, like the overall compliance tensor, always hermitian. The average strain and linear momentum tensors are, however, coupled and the coupling tensors are shown to be each others' hermitian transpose. Finally we present a numerical example of a 3-D periodic composite composed of elastic cubes periodically distributed in an elastic matrix. The presented results corroborate the predictions of the theoretical treatment.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Universal Theorems for Total Energy of the Dynamics of Linearly Elastic Heterogeneous Solids

In this paper we consider a sample of a linearly elastic heterogeneous composite in elastodynamic equilibrium and present universal theorems which provide lower bounds for the total elastic strain energy plus the kinetic energy, and the total complementary elastic energy plus the kinetic energy. For a general heterogeneous sample which undergoes harmonic motion at a single frequency, we show that, among all consistent boundary data which produce the same average strain, the uniform-stress boundary data render the total elastic strain energy plus the kinetic energy an absolute minimum. We also show that, among all consistent boundary data which produce the same average momentum in the sample, the uniform velocity boundary data render the total complementary elastic energy plus the kinetic energy an absolute minimum. We do not assume statistical homogeneity or material isotropy in our treatment, although they are not excluded. These universal theorems are the dynamic equivalent of the universal theorems already known for the static case (Nemat-Nasser and Hori 1995). It is envisaged that the bounds on the total energy presented in this paper will be used to formulate computable bounds on the overall dynamic properties of linearly elastic heterogeneous composites with arbitrary microstructures.

0 citations0 reviewsNo code linkedOpen access