Search

The manipulation of mechanical waves is a long-standing challenge for scientists and engineers, as numerous devices require their control. The current forefront of research in the control of classical waves has emerged from a seemingly unrelated field, namely, non-Hermitian quantum mechanics. By drawing analogies between this theory and those of classical systems, researchers have discovered phenomena that defy conventional intuition and have exploited them to control light, sound, and elastic waves. Here, we provide a brief perspective on recent developments, challenges and intricacies that distinguish non-Hermitian elastodynamics from optics and acoustics. We close this perspective with an outlook on potential directions such as topological phases in non-Hermitian elastodynamics and broken Hermitian symmetry in materials with electromomentum couplings.

Algorithmic materials discovery is a multi-disciplinary domain that integrates insights from specialists in alloy design, synthesis, characterization, experimental methodologies, computational modeling, and optimization. Central to this effort is a robust data management system paired with an interactive work platform. This platform should empower users to not only access others data but also integrate their analyses, paving the way for sophisticated data pipelines. To realize this vision, there is a need for an integrative collaboration platform, streamlined data sharing and analysis tools, and efficient communication channels. Such a collaborative mechanism should transcend geographical barriers, facilitating remote interaction and fostering a challenge-response dynamic. In this paper, we present our ongoing efforts in addressing the critical challenges related to an accelerated Materials Discovery Framework as a part of the High-Throughput Materials Discovery for Extreme Conditions Initiative. Our BIRDSHOT Center has successfully harnessed various tools and strategies, including the utilization of cloud-based storage, a standardized sample naming convention, a structured file system, the implementation of sample travelers, a robust sample tracking method, and the incorporation of knowledge graphs for efficient data management. Additionally, we present the development of a data collection platform, reinforcing seamless collaboration among our team members. In summary, this paper provides an illustration and insight into the various elements of an efficient and effective workflow within an accelerated materials discovery framework while highlighting the dynamic and adaptable nature of the data management tools and sharing platforms., Comment: 28 pages, 9 figures, 2 tables, with appendix that has 8 pages, accepted for publication at IMMI

In this study, we introduce a groundbreaking framework for materials discovery, we efficiently navigate a vast phase space of material compositions by leveraging Batch Bayesian statistics in order to achieve specific performance objectives. This approach addresses the challenge of identifying optimal materials from an untenably large array of possibilities in a reasonable timeframe with high confidence. Crucially, our batchwise methods align seamlessly with existing material processing infrastructure for synthesizing and characterizing materials. By applying this framework to a specific high entropy alloy system, we demonstrate its versatility and robustness in optimizing properties like strain hardening, hardness, and strain rate sensitivity. The fact that the Bayesian model is adept in refining and expanding the property Pareto front highlights its broad applicability across various materials, including steels, shape memory alloys, ceramics, and composites. This study advances the field of materials science and sets a new benchmark for material discovery methodologies. By proving the effectiveness of Bayesian optimization, we showcase its potential to redefine the landscape of materials discovery., Comment: 12 pages, 6 figures, with Supplementary Appendix that has 17 pages, 9 figures

This research focuses on the identification and causality analysis of coherent structures that arise in turbulent flows in square and rectangular ducts. Coherent structures are first identified from direct numerical simulation data via proper orthogonal decomposition (POD), both by using all velocity components, and after separating the streamwise and secondary components of the flow. The causal relations between the mode coefficients are analysed using pairwise-conditional Granger causality analysis. We also formulate a nonlinear Granger causality analysis that can account for nonlinear interactions between modes. Focusing on streamwise-constant structures within a duct of short streamwise extent, we show that the causal relationships are highly sensitive to whether the mode coefficients or their squared values are considered, whether nonlinear effects are explicitly accounted for, and whether streamwise and secondary flow structures are separated prior to causality analyses. We leverage these sensitivities to determine that linear mechanisms underpin causal relationships between modes that share the same symmetry or anti-symmetry properties about the corner bisector, while nonlinear effects govern the causal interactions between symmetric and antisymmetric modes. In all cases, we find that the secondary flow fluctuations (manifesting as streamwise vorticial structures) are the primary cause of both the presence and movement of near-wall streaks towards and away from the duct corners.

In this paper, we introduce a new approach for soft robot shape formation and morphing using approximate distance fields. The method uses concepts from constructive solid geometry, R-functions, to construct an approximate distance function to the boundary of a domain in $\Re^d$. The gradients of the R-functions can then be used to generate control algorithms for shape formation tasks for soft robots. By construction, R-functions are smooth and convex everywhere, possess precise differential properties, and easily extend from $\Re^2$ to $\Re^3$ if needed. Furthermore, R-function theory provides a straightforward method to creating composite distance functions for any desired shape by combining subsets of distance functions. The process is highly efficient since the shape description is an analytical expression, and in this sense, it is better than competing control algorithms such as those based on potential fields. Although the method could also apply to swarm robots, in this paper it is applied to soft robots to demonstrate shape formation and morphing in 2-D (simulation and experimentation) and 3-D (simulation).

In this paper, we present a powerful method (Atomistic Green's Function, AGF) for calculating the effective Hamiltonian of acoustic and elastic wave-scatterers. The ability to calculate the effective Hamiltonian allows for the study of scattering problems in infinite systems without the introduction of any artificial truncating boundaries such as perfectly matched layers or Dirichlet to Neumann (DtN) maps. Furthermore, the AGF formalism also allows for the efficient calculation of the Green's function of the scatterer as well as all relevant scattering metrics including reflection and transmission ratios. The formalism presented here is especially suited to scattering problems involving waveguides, phononic crystals, metamaterials, and metasurfaces. We show the application of the method to three scattering problems: scattering from a slab (1D), scattering from a finite phononic crystal (1D), and scattering from defects in a waveguide (2D).

Physics-Informed Neural Networks (PINNs) are a class of deep learning neural networks that learn the response of a physical system without any simulation data, and only by incorporating the governing partial differential equations (PDEs) in their loss function. While PINNs are successfully used for solving forward and inverse problems, their accuracy decreases significantly for parameterized systems. PINNs also have a soft implementation of boundary conditions resulting in boundary conditions not being exactly imposed everywhere on the boundary. With these challenges at hand, we present first-order physics-informed neural networks (FO-PINNs). These are PINNs that are trained using a first-order formulation of the PDE loss function. We show that, compared to standard PINNs, FO-PINNs offer significantly higher accuracy in solving parameterized systems, and reduce time-per-iteration by removing the extra backpropagations needed to compute the second or higher-order derivatives. Additionally, FO-PINNs can enable exact imposition of boundary conditions using approximate distance functions, which pose challenges when applied on high-order PDEs. Through three examples, we demonstrate the advantages of FO-PINNs over standard PINNs in terms of accuracy and training speedup., Comment: 16 pages, 11 figures, Selected for ML4PS workshop at NeurIPS 2022

The aim of this work is to analyse the formation mechanisms of large-scale coherent structures in the flow around a wall-mounted square cylinder, due to their impact on pollutant transport within cities. To this end, we assess causal relations between the modes of a reduced-order model obtained by applying proper-orthogonal decomposition to high-fidelity-simulation data of the flow case under study. The causal relations are identified using conditional transfer entropy, which is an information-theoretical quantity that estimates the amount of information contained in the past of one variable about another. This allows for an understanding of the origins and evolution of different phenomena in the flow, with the aim of identifying the modes responsible for the formation of the main vortical structures. Our approach unveils that vortex-breaker modes are the most causal modes, in particular, over higher-order modes, and no significant causal relationships were found for vortex-generator modes. We validate this technique by determining the causal relations present in the nine-equation model of near-wall turbulence developed by Moehlis et al. (New J. Phys, vol. 6, 2004, p. 56), which are in good agreement with literature results for turbulent channel flows., Comment: 19 pages, 10 figures

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.

In this work, a parameterized eigenvalue problem is analyzed for a phononic array in a 2D stress wave scattering setup, and a corresponding sensing application of this system is proposed to achieve source angle localization. The phononic domain consists of a periodic micro-structured medium, of which the eigen-wavevector band structure and the eigen-modes are exploited. The eigen-modes are naturally angle dependent due to changes in phases and periodic mode shapes determined by the incident angle. Intriguingly, the band exhibits angle-dependent transitions at the exceptional points (EPs) and critical angles (CAs), where the eigenvalues coincide or vanish. Coupled with these transitions, it is found that the eigen-modes switch their energy characteristics and symmetry patterns at these branch points, leading to enhanced angle dependence. Moreover, these eigen-modes also serve as the basis functions of the scattered waves. Therefore, the scattering response of the medium inherently possesses the angle-dependent properties, making this system naturally suitable for sensing applications. An artificial neural network (ANN) is trained with randomly weighted eigen-modes to achieve deep learning of the eigen features and angle dependence. The training data is derived only based on the eigen-modes of the unit cells. Nevertheless, the trained ANN can accurately identify the incident angle of an unknown scattering signal, with minimal side lobe levels and suppressed main lobe width. The ANN shows superior performance in comparison with standard delay-and-sum technique of estimating angle of arrival. The proposed application of ANN and micro-structured media highlights the physical importance of band structure topology and eigen-modes to a technological application, adds extra strength to the existing localization methods, and can be easily enhanced with the fast-growing data-driven techniques.

In this paper, we take an in-depth look at the analytical structure of the material transfer functions which govern acoustic and elastic response. These include wavenumber ($\kappa$) in such media and refractive index ($n$), density ($\boldsymbol{\rho}$) and its inverse, stiffness ($\boldsymbol{C}$) and compliance ($\boldsymbol{D}$) tensors as well as the Bulk modulus ($B$), and finally the broader generalization of these properties which is now known as the Willis tensor ($\boldsymbol{L}$). Our goal is to clarify the appropriate dispersion relations applicable to these properties from the perspective of passivity. Under some mild assumptions, causality ensures that these properties are analytical in the upper half but deriving dispersion relations for them requires one to know how they behave in the limit $|\omega|\rightarrow\infty$. Unlike electromagnetism, such a determination cannot be made on physical grounds since in that limit the continuum approximation breaks down. Instead, we can exploit the properties of the Herglotz-Nevanlinna (H-N) functions along with their tensorial counterparts which characterize the transfer functions of certain passive systems and for which the appropriate dispersion relation is known. Our aim, therefore, is to clarify the relationship that these transfer functions have with Herglotz functions, which in turn determines the appropriate dispersion relation for them. Our analysis shows that based upon passivity alone, dispersion relations of \emph{minimum} order 1 apply to the Fourier transforms of $\boldsymbol{D},\boldsymbol{\rho}, n'$, and the inverse of $B$, order 3 apply to $\boldsymbol{C},B$, and the inverse of $\boldsymbol{\rho}$, and order 2 applies to $\kappa$.

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 $\phi$, an approximate distance function to the boundary of a domain. To exactly impose homogeneous Dirichlet boundary conditions, the trial function is taken as $\phi$ 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., Comment: 50 pages, 46 figures

We clarify the claims and errors in the paper "Criterion for Negative Refraction with Low Optical Losses from a Fundamental Principle of Causality" by Mark Stockman. Contrary to the central assertion in that paper, simple examples consistent with the basic inequality which Stockman discovered show that it is possible to have negative refraction and low loss in an arbitrarily large frequency window. Further examination of the paper reveals additional errors that invalidate his argument that active materials cannot have low loss and negative refraction in a frequency window. Also, we point out that for active materials non-analyticity of the electrical permittivity in the upper half complex frequency plane does not necessarily imply noncausality, as Stockman infers., Comment: 6 pages, 1 figure

Computation of Mutual Information (MI) helps understand the amount of information shared between a pair of random variables. Automated feature selection techniques based on MI ranking are regularly used to extract information from sensitive datasets exceeding petabytes in size, over millions of features and classes. Series of one-vs-all MI computations can be cascaded to produce n-fold MI results, rapidly pinpointing informative relationships. This ability to quickly pinpoint the most informative relationships from datasets of billions of users creates privacy concerns. In this paper, we present Distributed Differentially Private Mutual Information (DDP-MI), a privacy-safe fast batch MI, across various scenarios such as feature selection, segmentation, ranking, and query expansion. This distributed implementation is protected with global model differential privacy to provide strong assurances against a wide range of privacy attacks. We also show that our DDP-MI can substantially improve the efficiency of MI calculations compared to standard implementations on a large-scale public dataset.

In this paper, we compute the band structure of one- and two-dimensional phononic composites using the extended finite element method (X-FEM) on structured higher-order (spectral) finite element meshes. On using partition-of-unity enrichment in finite element analysis, the X-FEM permits use of structured finite element meshes that do not conform to the geometry of holes and inclusions. This eliminates the need for remeshing in phononic shape optimization and topology optimization studies. In two dimensions, we adopt rational B{\'e}zier representation of curved (circular) geometries, and construct suitable material enrichment functions to model two-phase composites. A Bloch-formulation of the elastodynamic phononic eigenproblem is adopted. Efficient computation of weak form integrals with polynomial integrands is realized via the homogeneous numerical integration scheme -- a method that uses Euler's homogeneous function theorem and Stokes's theorem to reduce integration to the boundary of the domain. Ghost penalty stabilization is used on finite elements that are cut by a hole. Band structure calculations on perforated (circular holes, elliptical holes, and holes defined as a level set) materials as well as on two-phase phononic crystals are presented that affirm the sound accuracy and optimal convergence of the method on structured, higher-order spectral finite element meshes. Several numerical examples demonstrate the advantages of $p$-refinement made possible by the spectral extended finite element method. In these examples, fourth-order spectral extended finite elements deliver $\mathcal{O}(10^{-8})$ accuracy in frequency calculations with more than thirty-fold fewer degrees-of-freedom when compared to quadratic finite elements.

In this review, we take an extensive look at the role that the principles of causality and passivity have played in various areas of physics and engineering, including in the modern field of metamaterials. The aim is not to provide a comprehensive list of references as that number would be in the thousands, but to review the major results and contributions which have animated these areas and to provide a unified framework from which to understand the developments in different fields. Towards these goals, we chart the early history of the field through its dual beginnings in the analysis of the Sellmeier equation and in Hilbert transforms, giving rise to the far reaching dispersion relations in the early works of Kramers, Kronig, and Titchmarsh. However, these early relations constitute a limited result as they only apply to a restricted class of transfer functions. To understand how this restriction can be lifted, we take a quick detour into the distributional analysis of Schwartz. This approach expands the reach of the dispersion analysis to distributional transfer functions and also to those functions which exhibit polynomial growth properties. To generalize the results even further to tensorial transfer functions, we consider the concept of passivity - originally studied in the theory of electrical networks. We clarify why passivity implies causality. Subsequently, as special cases, we present examples of dispersion relations from several areas of physics including electromagnetism, acoustics, seismology, reflectance measurements, and scattering theory. We discuss sum rules, derivative analyticity relations, and nearly local approximations. Finally we review the clever applications of ideas from causality and passivity to the recent field of metamaterials. These ideas have provided limits to what can be achieved in metamaterial property design and metamaterial device performance.

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).

Deep generative models, such as Generative Adversarial Networks (GANs) or Variational Autoencoders (VAs) have been demonstrated to produce images of high visual quality. However, the existing hardware severely limits the size of the images that can be generated. The rapid growth of high dimensional data in many fields of science therefore poses a significant challenge for generative models. In cosmology, the large-scale, three-dimensional matter distribution, modeled with N-body simulations, plays a crucial role in understanding the evolution of the universe. As these simulations are computationally very expensive, GANs have recently generated interest as a possible method to emulate these datasets, but they have been, so far, mostly limited to two dimensional data. In this work, we introduce a new benchmark for the generation of three dimensional N-body simulations, in order to stimulate new ideas in the machine learning community and move closer to the practical use of generative models in cosmology. As a first benchmark result, we propose a scalable GAN approach for training a generator of N-body three-dimensional cubes. Our technique relies on two key building blocks, (i) splitting the generation of the high-dimensional data into smaller parts, and (ii) using a multi-scale approach that efficiently captures global image features that might otherwise be lost in the splitting process. We evaluate the performance of our model for the generation of N-body samples using various statistical measures commonly used in cosmology. Our results show that the proposed model produces samples of high visual quality, although the statistical analysis reveals that capturing rare features in the data poses significant problems for the generative models. We make the data, quality evaluation routines, and the proposed GAN architecture publicly available at https://github.com/nperraud/3DcosmoGAN

In this paper, we consider the operator properties of various phononic eigenvalue problems. We aim to answer some fundamental questions about the eigenvalues and eigenvectors of phononic operators. These include questions about the potential real and complex nature of the eigenvalues, whether the eigenvectors form a complete basis, what are the right orthogonality relationships, and how to create a complete basis when none may exist at the outset. In doing so we present a unified understanding of the properties of the phononic eigenvalues and eigenvectors which would emerge from any numerical method employed to compute such quantities. We show that the phononic problem can be cast into linear eigenvalue forms from which such quantities as frequencies, wavenumbers, and desired components of wavevectors can be directly ascertained without resorting to searches or quadratic eigenvalue problems and that the relevant properties of such quantities can be determined apriori through the analysis of the associated operators. We further show how the Plane Wave Expansion (PWE) method may be extended to solve each of these eigenvalue forms, thus extending the applicability of the PWE method to cases beyond those which have been considered till now. The theoretical discussions are supplemented with supporting numerical calculations. The techniques and results presented here directly apply to wave propagation in other periodic systems such as photonics.

In this paper we report metamaterial properties including negative and singular effective properties for what would traditionally be considered non locally resonant 2-phase phononic unit cells. The negative effective material properties reported here occur well below the homogenization limit and are, therefore, acceptable descriptions of overall behavior. The material property combinations which make this possible were first revealed by a novel level set based topology optimization process which we describe. The optimization process revealed that a 2-phase unit cell in which one of the phases is simultaneously lighter and stiffer than the other results in dynamic behavior which has all the attendant characteristics of a locally resonant composite including negative effective properties far below the homogenization limit. We investigate this further using the Craig-Bampton decomposition and clarify that these properties emerge through an interplay between the fundamental internal modeshape of the unit cell and a rigid body mode. Through explicit numerical calculations on 1-D, 2-phase unit cells, we show that negative effective properties only appear for the specific material property combination mentioned above. Furthermore, we provide a proof which supports this conclusion. The concept is also shown to hold for 2-D unit cells where we show that an appropriately designed hexagonal unit cell made of 2 material phases exhibits negative effective shear modulus and density in an appropriate frequency regime in which it also exhibits negative refraction. An important conclusion of this paper is that the class of unit cells expected to result in negative properties can be expanded beyond the classic unit cell (three-phase unit cells with an explicit locally resonant phase) to include topologically simpler 2-phase unit cells as well.

Integrated Computational Materials Engineering (ICME) calls for the integration of computational tools into the materials and parts development cycle, while the Materials Genome Initiative (MGI) calls for the acceleration of the materials development cycle through the combination of experiments, simulation, and data. As they stand, both ICME and MGI do not prescribe how to achieve the necessary tool integration or how to efficiently exploit the computational tools, in combination with experiments, to accelerate the development of new materials and materials systems. This paper addresses the first issue by putting forward a framework for the fusion of information that exploits correlations among sources/models and between the sources and `ground truth'. The second issue is addressed through a multi-information source optimization framework that identifies, given current knowledge, the next best information source to query and where in the input space to query it via a novel value-gradient policy. The querying decision takes into account the ability to learn correlations between information sources, the resource cost of querying an information source, and what a query is expected to provide in terms of improvement over the current state. The framework is demonstrated on the optimization of a dual-phase steel to maximize its strength-normalized strain hardening rate. The ground truth is represented by a microstructure-based finite element model while three low fidelity information sources---i.e. reduced order models---based on different homogenization assumptions---isostrain, isostress and isowork---are used to efficiently and optimally query the materials design space., Comment: 19 pages, 11 figures, 5 tables

We show that deep convolutional neural networks (CNN) can massively outperform traditional densely-connected neural networks (both deep or shallow) in predicting eigenvalue problems in mechanics. In this sense, we strike out in a new direction in mechanics computations with strongly predictive NNs whose success depends not only on architectures being deep, but also being fundamentally different from the widely-used to date. We consider a model problem: predicting the eigenvalues of 1-D and 2-D phononic crystals. For the 1-D case, the optimal CNN architecture reaches $98\%$ accuracy level on unseen data when trained with just 20,000 samples, compared to $85\%$ accuracy even with $100,000$ samples for the typical network of choice in mechanics research. We show that, with relatively high data-efficiency, CNNs have the capability to generalize well and automatically learn deep symmetry operations, easily extending to higher dimensions and our 2D case. Most importantly, we show how CNNs can naturally represent mechanical material tensors, with its convolution kernels serving as local receptive fields, which is a natural representation of mechanical response. Strategies proposed are applicable to other mechanics' problems and may, in the future, be used to sidestep cumbersome algorithms with purely data-driven approaches based upon modern deep architectures.

In this paper we consider the problem of avoided crossings (level repulsion) in phononic crystals and suggest a computationally efficient strategy to distinguish them from normal cross points. This process is essential for the correct sorting of the phononic bands and, subsequently, for the accurate determination of mode continuation, group velocities, and emergent properties which depend on them such as thermal conductivity. Through explicit phononic calculations using generalized Rayleigh quotient, we identify exact locations of exceptional points in the complex wavenumber domain which results in level repulsion in the real domain. We show that in the vicinity of the exceptional point the relevant phononic eigenvalue surfaces resemble the surfaces of a 2 by 2 parameter-dependent matrix. Along a closed loop encircling the exceptional point we show that the phononic eigenvalues are exchanged, just as they are for the 2 by 2 matrix case. However, the behavior of the associated eigenvectors is shown to be more complex in the phononic case. Along a closed loop around an exceptional point, we show that the eigenvectors can flip signs multiple times unlike a 2 by 2 matrix where the flip of sign occurs only once. Finally, we exploit these eigenvector sign flips around exceptional points to propose a simple and efficient method of distinguishing them from normal crosses and of correctly sorting the band-structure. Our proposed method is roughly an order-of magnitude faster than the zoom-in method and correctly identifies > 97% of the cases considered. Both its speed and accuracy can be further improved and we suggest some ways of achieving this. Our method is general and, as such, would be directly applicable to other eigenvalue problems where the eigenspectrum needs to be correctly sorted.

All metamaterial applications are based upon the idea that extreme material properties can be achieved through appropriate dynamic homogenization of composites. This homogenization is almost always done for infinite domains and the results are then applied to finite samples. This process ignores the evanescent waves which appear at the boundaries of such finite samples. In this paper we first clarify the emergence and purpose of these evanescent waves in a model problem consisting of an interface between a layered composite and a homogeneous medium. We show that these evanescent waves form boundary layers on either side of the interface beyond which the composite can be represented by appropriate infinite domain homogenized relations. We show that if one ignores the boundary layers then the displacement and stress fields are discontinuous across the interface. Therefore, the scattering coefficients at such an interface cannot be determined through the conventional continuity conditions involving only propagating modes. Here we propose an approximate variational approach for sidestepping these boundary layers. The aim is to determine the scattering coefficients without the knowledge of evanescent modes. Through various numerical examples we show that our technique gives very good estimates of the actual scattering coefficients beyond the long wavelength limit.

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.

In this paper the salient features of the Plane Wave Expansion (PWE) method and the mixed variational technique are combined for the fast eigenvalue computations of arbitrarily complex phononic unit cells. This is done by expanding the material properties in a Fourier expansion, as is the case with PWE. The required matrix elements in the variational scheme are identified as the discrete Fourier transform coefficients of material properties, thus obviating the need for any explicit integration. The process allows us to provide succinct and closed form expressions for all the matrices involved in the mixed variational method. The scheme proposed here preserves both the simplicity of expression which is inherent in the PWE method and the superior convergence properties of the mixed variational scheme. We present numerical results and comment upon the convergence and stability of the current method. We show that the current representation renders the results of the method stable over the entire range of the expansion terms as allowed by the spatial discretization. When compared with a zero order numerical integration scheme, the present method results in greater computational accuracy of all eigenvalues. A higher order numerical integration scheme comes close to the accuracy of the present method but only with significantly more computational expense., Comment: Submitted to Journal of Engineering Mechanics. arXiv admin note: text overlap with arXiv:1411.2996, arXiv:1310.6380

In this paper gradient based topology optimization (TO) is used to discover 3-D phononic structures that exhibit ultra-wide normalized all-angle all-mode band gaps. The challenging computational task of repeated 3-D phononic band-structure evaluations is accomplished by a combination of a fast mixed variational eigenvalue solver and distributed Graphic Processing Unit (GPU) parallel computations. The TO algorithm utilizes the material distribution-based approach and a gradient-based optimizer. The design sensitivity for the mixed variational eigenvalue problem is derived using the adjoint method and is implemented through highly efficient vectorization techniques. We present optimized results for two-material simple cubic (SC), body centered cubic (BCC), and face centered cubic (FCC) crystal structures and show that in each of these cases different initial designs converge to single inclusion network topologies within their corresponding primitive cells. The optimized results show that large phononic stop bands for bulk wave propagation can be achieved at lower than close packed spherical configurations leading to lighter unit cells. For tungsten carbide - epoxy crystals we identify all angle all mode normalized stop bands exceeding 100%, which is larger than what is possible with only spherical inclusions.

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., Comment: arXiv admin note: text overlap with arXiv:1411.2996

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., Comment: Accepted for publication in Proc. R. Soc. A

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., Comment: Paper under review

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, Comment: arXiv admin note: substantial text overlap with arXiv:1310.6380

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.

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.

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.

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.

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.

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., Comment: 18 pages, In Press (Mechanics of Materials)

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., Comment: 26 pages, 2 figures, submitted to Proceedings of the Royal Society A

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]., Comment: 23 pages, 6 figures, submitted to JMPS