Source author record

Robert Lipton

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

19works
11topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

19 published item(s)

preprint2023arXiv

Quasistatic fracture evolution

Nonlocal quasistatic fracture evolution for interacting cracks is developed and supporting numerical examples are presented. The approach is implicit and is based on local stationarity and fixed point methods. It is proved that the fracture evolution decreases stored elastic energy with each load step as the cracks advance; provided the load increments are chosen sufficiently small. This is also seen in the numerical examples. The numerical examples include evolution of a straight crack, a crack propagating inside an L-shaped domain, and two offset inward propagating cracks.

preprint2022arXiv

Bloch Spectra for High Contrast Elastic Media

Analytic representation formulas and power series are developed to describe the band structure inside periodic elastic crystals made from high contrast inclusions. We use source free modes associated with structural spectra to represent the solution operator of the Lamé system inside phononic crystals. Convergent power series for the Bloch wave spectrum are obtained using the representation formulas. An explicit bound on the convergence radius is given through the structural spectra of the inclusion array and the Dirichlet spectra of the inclusions. Sufficient conditions for the separation of spectral branches of the dispersion relation for any fixed quasi-momentum are identified.

preprint2021arXiv

Peridynamics-based discrete element method (PeriDEM) model of granular systems involving breakage of arbitrarily shaped particles

Usage, manipulation, transport, delivery, and mixing of granular or particulate media, comprised of spherical or polyhedral particles, is commonly encountered in industrial sectors of construction (cement and rock fragments), pharmaceutics (tablets), and transportation (ballast). Elucidating particulate media's behavior in concert with particle attrition (i.e., particle wear and subsequent particle fragmentation) is essential for predicting the performance and increasing the efficiency of engineering systems using such media. Discrete element method (DEM) based techniques can describe the interaction between particles but cannot model intra-particle deformation, especially intra-particle fracture. On the other hand, peridynamics provides the means to account for intra-particle deformation and fracture due to contact forces between particles. The present study proposes a hybrid model referred to as \textit{PeriDEM} that combines the advantages of peridynamics and DEM. The model parameters can be tuned to achieve desired DEM contact forces, damping effects, and intra-particle stiffness. Two particle impacts and compressive behavior of multi-particle systems are thoroughly investigated. The model can account for any arbitrarily shaped particle in general. Spherical, hexagonal, and non-convex particle shapes are simulated in the present study. The effect of mesh resolution on intra-particle peridynamics is explicitly studied. The proposed hybrid model opens a new avenue to explore the complicated interactions encountered in discrete particle dynamics that involve the formation of force chains, particle interlocking, particle attrition, wear, and the eventual breakage.

preprint2020arXiv

Classic dynamic fracture recovered as the limit of a nonlocal peridynamic model: The single edge notch in tension

A simple nonlocal field theory of peridynamic type is applied to model brittle fracture. The fracture evolution is shown to converge in the limit of vanishing nonlocality to classic plane elastodynamics with a running crack. The kinetic relation for the crack is recovered directly from the nonlocal model in the limit of vanishing nonlocality. We carry out our analysis for a single crack in a plate subject to mode one loading. The convergence is corroborated by numerical experiments.

preprint2019arXiv

Multiscale-Spectral GFEM and Optimal Oversampling

In this work we address the Multiscale Spectral Generalized Finite Element Method (MS-GFEM) developed in [I. Babuška and R. Lipton, Multiscale Modeling and Simulation 9 (2011), pp. 373--406]. We outline the numerical implementation of this method and present simulations that demonstrate contrast independent exponential convergence of MS-GFEM solutions. We introduce strategies to reduce the computational cost of generating the optimal oversampled local approximating spaces used here. These strategies retain accuracy while reducing the computational work necessary to generate local bases. Motivated by oversampling we develop a nearly optimal local basis based on a partition of unity on the boundary and the associated A-harmonic extensions.

preprint2016arXiv

Bloch Waves in Crystals and Periodic High Contrast Media

Analytic representation formulas and power series are developed describing the band structure inside periodic photonic and acoustic crystals made from high contrast inclusions. Central to this approach is the identification and utilization of a resonance spectrum for quasi-periodic source free modes. These modes are used to represent solution operators associated with electromagnetic and acoustic waves inside periodic high contrast media. Convergent power series for the Bloch wave spectrum is recovered from the representation formulas. Explicit conditions on the contrast are found that provide lower bounds on the convergence radius. These conditions are sufficient for the separation of spectral branches of the dispersion relation.

preprint2016arXiv

Complex fracture nucleation and evolution with nonlocal elastodynamics

A mechanical model is introduced for predicting the initiation and evolution of complex fracture patterns without the need for a damage variable or law. The model, a continuum variant of Newton's second law, uses integral rather than partial differential operators where the region of integration is over finite domain. The force interaction is derived from a novel nonconvex strain energy density function, resulting in a nonmonotonic material model. The resulting equation of motion is proved to be mathematically well-posed. The model has the capacity to simulate nucleation and growth of multiple, mutually interacting dynamic fractures. In the limit of zero region of integration, the model reproduces the classic Griffith model of brittle fracture. The simplicity of the formulation avoids the need for supplemental kinetic relations that dictate crack growth or the need for an explicit damage evolution law.

preprint2016arXiv

Creating band gaps in periodic media

We identify explicit conditions on geometry and material contrast for creating band gaps in 2-d photonic and 3-d acoustic crystals. This approach is novel and makes use of the electro-static and quasi-periodic source free resonances of the crystal. The source free modes deliver a spectral representation for solution operators associated with propagation of electromagnetic and acoustic waves inside periodic high contrast media. An accurate characterization of the quasi-periodic and electrostatic resonance spectrum in terms of the shape and geometry of the scatters is possible. This information together with the Dirichlet and a Neumann like spectra associated with the inclusions deliver conditions sufficient for opening band gaps at finite contrast. The theory provides a systematic means for the identification of photonic and phononic band gaps within a specified frequency range.

preprint2016arXiv

Uncertain Loading and Quantifying Maximum Energy Concentration within Composite Structures

We introduce a systematic method for identifying the worst case load among all boundary loads of fixed energy. Here the worst case load is defined to be the one that delivers the largest fraction of input energy to a prescribed subdomain of interest. The worst case load is identified with the first eigenfunction of a suitably defined eigenvalue problem. The first eigenvalue for this problem is the maximum fraction of boundary energy that can be delivered to the subdomain. We compute worst case boundary loads and associated energy contained inside a prescribed subdomain through the numerical solution of the eigenvalue problem. We apply this computational method to bound the worst case load associated with an ensemble of random boundary loads given by a second order random process. Several examples are carried out on heterogeneous structures to illustrate the method.

preprint2015arXiv

Cohesive Dynamics and Brittle Fracture

We formulate a nonlocal cohesive model for calculating the deformation state inside a cracking body. In this model a more complete set of physical properties including elastic and softening behavior are assigned to each point in the medium. We work within the small deformation setting and use the peridynamic formulation. Here strains are calculated as difference quotients. The constitutive relation is given by a nonlocal cohesive law relating force to strain. At each instant of the evolution we identify a process zone where strains lie above a threshold value. Perturbation analysis shows that jump discontinuities within the process zone can become unstable and grow. We derive an explicit inequality that shows that the size of the process zone is controlled by the ratio given by the length scale of nonlocal interaction divided by the characteristic dimension of the sample. The process zone is shown to concentrate on a set of zero volume in the limit where the length scale of nonlocal interaction vanishes with respect to the size of the domain. In this limit the dynamic evolution is seen to have bounded linear elastic energy and Griffith surface energy. The limit dynamics corresponds to the simultaneous evolution of linear elastic displacement and the fracture set across which the displacement is discontinuous. We conclude illustrating how the approach developed here can be applied to limits of dynamics associated with other energies that $Γ$- converge to the Griffith fracture energy.

preprint2014arXiv

Tuning gain and bandwidth of traveling wave tubes using metamaterial beam-wave interaction structures

We employ metamaterial beam-wave interaction structures for tuning the gain and bandwidth of short traveling wave tubes. The interaction structures are made from metal rings of uniform cross section, which are periodically deployed along the length of the traveling wave tube. The aspect ratio of the ring cross sections are adjusted to control both gain and bandwidth. The frequency of operation is controlled by the filling fraction of the ring cross section with respect to the period.

preprint2013arXiv

Dynamic Brittle Fracture as a Small Horizon Limit of Peridynamics

We consider the nonlocal formulation of continuum mechanics described by peridynamics. We provide a link between peridynamic evolution and brittle fracture evolution for a broad class of peridynamic potentials associated with unstable peridynamic constitutive laws. Distinguished limits of peridynamic evolutions are identified that correspond to vanishing peridynamic horizon. The limit evolution is associated with dynamic brittle fracture and satisfies a dynamic energy inequality expressed in terms of the kinetic energy of the motion together with a bulk elastic energy and a Griffith surface energy. It corresponds to the simultaneous evolution of elastic displacement and brittle fracture with displacement fields satisfying the wave equation inside the cracking domain. The wave equation provides the dynamic coupling between elastic waves and the evolving fracture path inside the media. The elastic moduli, wave speed and energy release rate for the evolution are explicitly determined by moments of the peridynamic influence function and the peridynamic potential energy.

preprint2012arXiv

Double Negative Dispersion Relations from Coated Plasmonic Rods

A metamaterial with frequency dependent double negative effective properties is constructed from a sub-wavelength periodic array of coated rods. Explicit power series are developed for the dispersion relation and associated Bloch wave solutions. The expansion parameter is the ratio of the length scale of the periodic lattice to the wavelength. Direct numerical simulations for finite size period cells show that the leading order term in the power series for the dispersion relation is a good predictor of the dispersive behavior of the metamaterial.

preprint2010arXiv

Multiscale Analysis of Heterogeneous Media in the Peridynamic Formulation

A methodology is presented for investigating the dynamics of heterogeneous media using the nonlocal continuum model given by the peridynamic formulation. The approach presented here provides the ability to model the macroscopic dynamics while at the same time resolving the dynamics at the length scales of the microstructure. Central to the methodology is a novel two-scale evolution equation. The rescaled solution of this equation is shown to provide a strong approximation to the actual deformation inside the peridynamic material. The two scale evolution can be split into a microscopic component tracking the dynamics at the length scale of the heterogeneities and a macroscopic component tracking the volume averaged (homogenized) dynamics. The interplay between the microscopic and macroscopic dynamics is given by a coupled system of evolution equations. The equations show that the forces generated by the homogenized deformation inside the medium are related to the homogenized deformation through a history dependent constitutive relation

preprint2010arXiv

Optimal Local Approximation Spaces for Generalized Finite Element Methods with Application to Multiscale Problems

The paper addresses a numerical method for solving second order elliptic partial differential equations that describe fields inside heterogeneous media. The scope is general and treats the case of rough coefficients, i.e. coefficients with values in $L^\infty(Ω)$. This class of coefficients includes as examples media with micro-structure as well as media with multiple non-separated length scales. The approach taken here is based on the the generalized finite element method (GFEM) introduced in \cite{107}, and elaborated in \cite{102}, \cite{103} and \cite{104}. The GFEM is constructed by partitioning the computational domain $Ω$ into to a collection of preselected subsets $ω_{i},i=1,2,..m$ and constructing finite dimensional approximation spaces $Ψ_{i}$ over each subset using local information. The notion of the Kolmogorov $n$-width is used to identify the optimal local approximation spaces. These spaces deliver local approximations with errors that decay almost exponentially with the degrees of freedom $N_{i}$ in the energy norm over $ω_i$. The local spaces $% Ψ_{i}$ are used within the GFEM scheme to produce a finite dimensional subspace $S^N$ of $H^{1}(Ω)$ which is then employed in the Galerkin method. It is shown that the error in the Galerkin approximation decays in the energy norm almost exponentially (i.e., super-algebraicly) with respect to the degrees of freedom $N$. When length scales "`separate" and the microstructure is sufficiently fine with respect to the length scale of the domain $ω_i$ it is shown that homogenization theory can be used to construct local approximation spaces with exponentially decreasing error in the pre-asymtotic regime.

preprint2010arXiv

Representation formulas for $L^\infty$ norms of weakly convergent sequences of gradient fields in homogenization

We examine the composition of the $L^{\infty}$ norm with weakly convergent sequences of gradient fields associated with the homogenization of second order divergence form partial differential equations with measurable coefficients. Here the sequences of coefficients are chosen to model heterogeneous media and are piecewise constant and highly oscillatory. We identify local representation formulas that in the fine phase limit provide upper bounds on the limit superior of the $L^{\infty}$ norms of gradient fields. The local representation formulas are expressed in terms of the weak limit of the gradient fields and local corrector problems. The upper bounds may diverge according to the presence of rough interfaces. We also consider the fine phase limits for layered microstructures and for sufficiently smooth periodic microsturctures. For these cases we are able to provide explicit local formulas for the limit of the $L^\infty$ norms of the associated sequence of gradient fields. Local representation formulas for lower bounds are obtained for fields corresponding to continuously graded periodic microstructures as well as for general sequences of oscillatory coefficients. The representation formulas are applied to problems of optimal material design.

preprint2010arXiv

Tunable Double Negative Band Structure from Non-Magnetic Coated Rods

A system of periodic poly-disperse coated nano-rods is considered. Both the coated nano-rods and host material are non-magnetic. The exterior nano-coating has a frequency dependent dielectric constant and the rod has a high dielectric constant. A negative effective magnetic permeability is generated near the Mie resonances of the rods while the coating generates a negative permittivity through a field resonance controlled by the plasma frequency of the coating and the geometry of the crystal. The explicit band structure for the system is calculated in the sub-wavelength limit. Tunable pass bands exhibiting negative group velocity are generated and correspond to simultaneously negative effective dielectric permittivity and magnetic permeability. These can be explicitly controlled by adjusting the distance between rods, the coating thickness, and rod diameters.

preprint2009arXiv

Sub-Wavelength Plasmonic Crystals: Dispersion Relations and Effective Properties

We obtain a convergent power series expansion for the first branch of the dispersion relation for subwavelength plasmonic crystals consisting of plasmonic rods with frequency-dependent dielectric permittivity embedded in a host medium with unit permittivity. The expansion parameter is $η=kd=2πd/λ$, where $k$ is the norm of a fixed wavevector, $d$ is the period of the crystal and $λ$ is the wavelength, and the plasma frequency scales inversely to $d$, making the dielectric permittivity in the rods large and negative. The expressions for the series coefficients (a.k.a., dynamic correctors) and the radius of convergence in $η$ are explicitly related to the solutions of higher-order cell problems and the geometry of the rods. Within the radius of convergence, we are able to compute the dispersion relation and the fields and define dynamic effective properties in a mathematically rigorous manner. Explicit error estimates show that a good approximation to the true dispersion relation is obtained using only a few terms of the expansion. The convergence proof requires the use of properties of the Catalan numbers to show that the series coefficients are exponentially bounded in the $H^1$ Sobolev norm.