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

Graeme W. Milton

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

ResearcherUnclaimed source record
math-phmath.MPmath.APcond-mat.mtrl-sciphysics.class-phphysics.opticsgr-qcmath.AGmath.CVmath.NAmath.OCNumerical Analysisphysics.gen-phphysics.med-ph

Catalog footprint

What is connected

52works
14topics
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

52 published item(s)

preprint2026arXiv

Broadband quasistatic passive cloaking: bounds and limitations in the near-field regime

We consider here several aspects of the following challenging question: is it possible to use a passive cloak to make invisible a dielectric inclusion on a finite frequency interval in the quasistatic regime of Maxwell's equations for an observer close to the object? In this work, by considering the Dirichlet-to-Neumann (DtN) map, we not only answer negatively this question, but we go further and provide some quantitative bounds on this map that provide fundamental limits to both cloaking as well as approximate cloaking. These bounds involve the following physical parameters: the length and center of the frequency interval, the volume of the cloaking device, the volume of the obstacle, and the relative permittivity of the object. Our approach is based on two key tools: i) variational principles from the abstract theory of composites and ii) the analytic approach to deriving bounds from sum rules for passive systems. To use i), we prove a new representation theorem for the DtN map which allows us to interpret this map as an effective operator in the abstract theory of composites. One important consequence of this representation is that it allows one to incorporate the broad and deep results from the theory of composites, such as variational principles, and to apply the bounds derived from them to the DtN map. These results could be useful in other contexts other than cloaking. Next, to use ii), we show that the passivity assumption allows us to connect the DtN map (as function of the frequency) with two important classes of analytic functions, namely, Herglotz and Stieltjes functions. The sum rules for these functions, combined with the variational approach, allows us to derive new inequalities on the DtN map which impose fundamental limitations on passive cloaking, both exact and approximate, over a frequency interval. We consider both cases of lossy and lossless cloaks.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

A possible explanation of dark matter and dark energy involving a vector torsion field

A simple gravitational model with torsion is studied, and it is suggested that it could explain the dark matter and dark energy in the universe. It can be reinterpreted as a model using the Einstein gravitational equations where spacetime has regions filled with a perfect fluid with negative energy (pressure) and positive mass density, other regions containing an anisotropic substance that in the rest frame (where the momentum is zero) has negative mass density and a uniaxial stress tensor, and possibly other "luminal" regions where there is no rest frame. The torsion vector field is inhomogeneous throughout spacetime, and possibly turbulent. Numerical simulations should reveal whether or not the equations are consistent with cosmological observations of dark matter and dark energy.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

A unifying perspective on linear continuum equations prevalent in science. Part II: Canonical forms for time-harmonic equations

Following some past advances, we reformulate a large class of linear continuum science equations in the format of the extended abstract theory of composites so that we can apply this theory to better understand and efficiently solve those equations. Here in part II we elucidate the form for many time-harmonic equations that do not involve higher order gradients.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

An Energy Conserving Mechanism for Temporal Metasurfaces

Changing the microstructure properties of a space-time metamaterial while a wave is propagating through it, in general requires addition or removal of energy, which can be of exponential form depending on the type of modulation. This limits the realization and application of space-time metamaterials. We resolve this issue by introducing a novel mechanism of conserving energy at temporal metasurfaces in a non-linear setting. The idea is first demonstrated by considering a wave-packet propagating in a discrete medium of 1-d chain of springs and masses, where using our energy conserving mechanism we show that the spring stiffness can be incremented at several time interfaces and the energy will still be conserved. We then consider an interesting application of time-reversed imaging in 1-d and 2-d spring-mass systems with a wave packet traveling in the homogenized regime. Our numerical simulations show that, in 1-d, when the wave packet hits the time-interface two sets of waves are generated, one traveling forward in time and the other traveling backward. The time-reversed waves re-converge at the location of the source and we observe its regeneration. In 2-d, we use more complicated initial shapes and, even then, we observe regeneration of the original image or source. Thus, we achieve time-reversed imaging with conservation of energy in a non-linear system. The energy conserving mechanism can be easily extended to continuum media.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Limit analysis of strut nets

Truss structures composed of members that work exclusively in tension or in compression appear in several problems of science and engineering, e.g., in the study of the resisting mechanisms of masonry structures, as well as in the design of spider web-inspired web structures. This work generalizes previous results on the existence of cable webs that are able to support assigned sets of nodal forces under tension. We extend such a problem to the limit analysis of compression-only 'strut nets' subjected to fixed and variable nodal loads. These systems provide discrete element models of masonry bodies, which lie inside the polygon/polyhedron with vertices at the points of application of the given forces ('underlying masonry structures'). It is assumed that fixed nodal forces are combined with variable forces growing proportionally to a scalar multiplier (load multiplier), and that the supporting strut net is subjected to kinematic constraints at given nodal positions.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

An extremal problem arising in the dynamics of two-phase materials that directly reveals information about the internal geometry

In two phase materials, each phase having a non-local response in time, it has been found that for some driving fields the response somehow untangles at specific times, and allows one to directly infer useful information about the geometry of the material, such as the volume fractions of the phases. Motivated by this, and to obtain an algorithm for designing appropriate driving fields, we find approximate, measure independent, linear relations between the values that Markov functions take at a given set of possibly complex points, not belonging to the interval [-1,1] where the measure is supported. The problem is reduced to simply one of polynomial approximation of a given function on the interval [-1,1] and to simplify the analysis Chebyshev approximation is used. This allows one to obtain explicit estimates of the error of the approximation, in terms of the number of points and the minimum distance of the points to the interval [-1,1]. Assuming this minimum distance is bounded below by a number greater than 1/2, the error converges exponentially to zero as the number of points is increased. Approximate linear relations are also obtained that incorporate a set of moments of the measure. In the context of the motivating problem, the analysis also yields bounds on the response at any particular time for any driving field, and allows one to estimate the response at a given frequency using an appropriately designed driving field that effectively is turned on only for a fixed interval of time. The approximation extends directly to Markov-type functions with a positive semidefinite operator valued measure, and this has applications to determining the shape of an inclusion in a body from boundary flux measurements at a specific time, when the time-dependent boundary potentials are suitably tailored.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A unifying perspective on linear continuum equations prevalent in physics. Part V: resolvents; bounds on their spectrum; and their Stieltjes integral representations when the operator is not selfadjoint

We consider resolvents of operators taking the form ${\bf A}=Γ_1{\bf B}Γ_1$ where $Γ_1({\bf k})$ is a projection that acts locally in Fourier space and ${\bf B}({\bf x})$ is an operator that acts locally in real space. Such resolvents arise naturally when one wants to solve any of the large class of linear physical equations surveyed in Parts I, II, III, and IV that can be reformulated as problems in the extended abstract theory of composites. We review how $Q^*$-convex operators can be used to bound the spectrum of ${\bf A}$. Then, based on the Cherkaev-Gibiansky transformation and subsequent developments, that we reformulate, we obtain for non-Hermitian ${\bf B}$ a Stieltjes type integral representation for the resolvent $(z_0{\bf I}-{\bf A})^{-1}$. The representation holds in the half plane $\Re(e^{i\vartheta}z_0)>c$, where $\vartheta$ and $c$ are such that $c{\bf I}-[e^{i\vartheta}{\bf B}+e^{-i\vartheta}{\bf B}^\dagger]$ is positive definite (and coercive).

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A unifying perspective on linear continuum equations prevalent in science. Part I: Canonical forms for static, steady, and quasistatic equations

Following some past advances, we reformulate a large class of linear continuum science equations in the format of the extended abstract theory of composites so that we can apply this theory to better understand and efficiently solve those equations. Here in part I we elucidate the form for many static, steady, and quasistatic equations.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A unifying perspective on linear continuum equations prevalent in science. Part III: Canonical forms for dynamic equations with moduli that may, or may not, vary with time

Enlarging on Parts I and II we write more equations in the desired format of the extended abstract theory of composites. We focus on a multitude of full dynamic equations, including equations where the medium is moving or otherwise changing in time. The motivation is that results and methods in the theory of composites then extend to these equations.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A unifying perspective on linear continuum equations prevalent in science. Part IV: Canonical forms for equations involving higher order gradients

Enlarging on Parts I, II, and III we write more equations in the desired format of the extended abstract theory of composites. We focus on a multitude of equations involving higher order derivatives. The motivation is that results and methods in the theory of composites then extend to these equations.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A unifying perspective on linear continuum equations prevalent in science. Part VI: rapidly converging series expansions for their solution

We obtain rapidly convergent series expansions of resolvents of operators taking the form ${\bf A}=Γ_1{\bf B}Γ_1$ where $Γ_1({\bf k})$ is a projection that acts locally in Fourier space and ${\bf B}({\bf x})$ is an operator that acts locally in real space. Such resolvents arise naturally when one wants to solve any of the large class of linear physical equations surveyed in Parts I, II, III, and IV that can be reformulated as problems in the extended abstract theory of composites. We show how the information about the spectrum of ${\bf A}$ can be used to greatly improve the convergence rate.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Inclusions of general shapes having constant field inside the core and non-elliptical neutral coated inclusions with anisotropic conductivity

For certain shapes of inclusions embedded in a body, the field inside the inclusion is uniform for some boundary condition. We provide a construction scheme for inclusions of general shapes having such a uniformity property in two dimensions based on the conformal mapping technique for the potential problem. Using this complex analysis method, we also design non-elliptical neutral coated inclusions with anisotropic conductivity. Neutral coated inclusions do not perturb a background uniform field when they are inserted into a homogeneous matrix. Although coated inclusions of various shapes are neutral to a single field, only concentric ellipses or confocal ellipsoids can be neutral to all uniform fields. This paper presents our work relating to the construction of non-elliptical coated inclusions with anisotropic conductivity in two dimensions that are neutral to all uniform fields, where the assignment of the flux condition on the boundary of the core depends on the applied background field. Using these neutral inclusions, we obtain cylindrical neutral inclusions in three dimensions, with no flux applied to the boundary of the core and with the anisotropic conductivity function of the shell given in accordance with the background uniform field.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Substitution of subspace collections with nonorthogonal subspaces to accelerate Fast Fourier Transform methods applied to conducting composites

We show the power of the algebra of subspace collections developed in Chapter 7 of the book "Extending the Theory of Composites to Other Areas of Science (edited by Milton, 2016). Specifically we accelerate the Fast Fourier Transform schemes of Moulinec and Suquet and Eyre and Milton (1994, 1998) for computing the fields and effective tensor in a conducting periodic medium by substituting a subspace collection with nonorthogonal subspaces inside one with orthogonal subspaces. This can be done when the effective conductivity as a function of the conductivity $σ_1$ of the inclusion phase (with the matrix phase conductivity set to $1$) has its singularities confined to an interval $[-β,-α]$ of the negative real $σ_1$ axis. Numerical results of Moulinec and Suquet show accelerated convergence for the model example of a square array of squares at $25\%$ volume fraction. For other problems we show how $Q^*_C$-convex functions can be used to restrict the region where singularities of the effective tensor as a function of the component tensors might be found.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

A rigorous approach to the field recursion method for two-component composites with isotropic phases

In this chapter of the book entitled, "Extending the Theory of Composites to Other Areas of Science" [edited by Graeme W. Milton, 2016] we give a rigorous derivation of the field equation recursion method in the abstract theory of composites to two-component composites with isotropic phases. This method is of great interest since it has proven to be a powerful tool in developing sharp bounds for the effective tensor of a composite material. The reason is that the effective tensor $\bf L_*$ can be interpreted in the general framework of the abstract theory of composites as the $Z$-operator on a certain orthogonal $Z(2)$ subspace collection. The base case of the recursion starts with an orthogonal $Z(2)$ subspace collection on a Hilbert space $\cal H$, the $Z$-problem, and the associated $Y$-problem. We provide some new conditions for the solvability of both the $Z$-problem and the associated $Y$-problem. We also give explicit representations of the associated $Z$-operator and $Y$-operator and study their analytical properties. An iteration method is then developed from a hierarchy of subspace collections and their associated operators which leads to a continued fraction representation of the initial effective tensor $\bf L_*$.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Analytic Materials

The theory of inhomogeneous analytic materials is developed. These are materials where the coefficients entering the equations involve analytic functions. Three types of analytic materials are identified. The first two types involve an integer $p$. If $p$ takes its maximum value then we have a complete analytic material. Otherwise it is incomplete analytic material of rank $p$. For two-dimensional materials further progress can be made in the identification of analytic materials by using the well-known fact that a $90^\circ$ rotation applied to a divergence free field in a simply connected domain yields a curl-free field, and this can then be expressed as the gradient of a potential. Other exact results for the fields in inhomogeneous media are reviewed. Also reviewed is the subject of metamaterials, as these materials provide a way of realizing desirable coefficients in the equations.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Analyticity of the Dirichlet-to-Neumann map for the time-harmonic Maxwell's equations

In this chapter of the book entitled, "Extending the Theory of Composites to Other Areas of Science" [edited by Graeme W. Milton, 2016] we derive the analyticity properties of the electromagnetic Dirichlet-to-Neumann map for the time-harmonic Maxwell's equations for passive linear multicomponent media. Moreover, we discuss the connection of this map to Herglotz functions for isotropic and anisotropic multicomponent composites.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

High frequency homogenization for travelling waves in periodic media

We consider high frequency homogenization in periodic media for travelling waves of several different equations: the wave equation for scalar-valued waves such as acoustics; the wave equation for vector-valued waves such as electromagnetism and elasticity; and a system that encompasses the Schr{ö}dinger equation. This homogenization applies when the wavelength is of the order of the size of the medium periodicity cell. The travelling wave is assumed to be the sum of two waves: a modulated Bloch carrier wave having crystal wave vector $\Bk$ and frequency $ω_1$ plus a modulated Bloch carrier wave having crystal wave vector $\Bm$ and frequency $ω_2$. We derive effective equations for the modulating functions, and then prove that there is no coupling in the effective equations between the two different waves both in the scalar and the system cases. To be precise, we prove that there is no coupling unless $ω_1=ω_2$ and $(\Bk-\Bm)\odotΛ\in 2π\mathbb Z^d,$ where $Λ=(λ_1λ_2\dotsλ_d)$ is the periodicity cell of the medium and for any two vectors $a=(a_1,a_2,\dots,a_d), b=(b_1,b_2,\dots,b_d)\in\mathbb R^d,$ the product $a\odot b$ is defined to be the vector $(a_1b_1,a_2b_2,\dots,a_db_d).$ This last condition forces the carrier waves to be equivalent Bloch waves meaning that the coupling constants in the system of effective equations vanish. We use two-scale analysis and some new weak-convergence type lemmas. The analysis is not at the same level of rigor as that of Allaire and coworkers who use two-scale convergence theory to treat the problem, but has the advantage of simplicity which will allow it to be easily extended to the case where there is degeneracy of the Bloch eigenvalue.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

On ideal dynamic climbing ropes

We consider the rope climber fall problem in two different settings. The simplest formulation of the problem is when the climber falls from a given altitude and is attached to one end of the rope while the other end of the rope is attached to the rock at a given height. The problem is then finding the properties of the rope for which the peak force felt by the climber during the fall is minimal. The second problem of our consideration is again minimizing the same quantity in the presence of a carabiner. We will call such ropes \textit{mathematically ideal.} Given the height of the carabiner, the initial height and the mass of the climber, the length of the unstretched rope, and the distance between the belayer and the carabineer, we find the optimal (in the sense of minimized the peak force to a given elongation) dynamic rope in the framework of nonlinear elasticity. Wires of shape memory materials have some of the desired features of the tension-strain relation of a mathematically ideal dynamic rope, namely a plateau in the tension over a range of strains. With a suitable hysteresis loop, they also absorb essentially all the energy from the fall, thus making them an ideal rope in this sense too.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Superfunctions and the algebra of subspace collections and their association with rational functions of several complex variables

A natural connection between rational functions of several real or complex variables, and subspace collections is explored. A new class of function, superfunctions, are introduced which are the counterpart to functions at the level of subspace collections. Operations on subspace collections are found to correspond to various operations on rational functions, such as addition, multiplication and substitution. It is established that every rational matrix valued function which is homogeneous of degree 1 can be generated from an appropriate, but not necessarily unique, subspace collection: the mapping from subspace collections to rational functions is onto, but not one to one. For some applications superfunctions may be more important than functions, as they incorporate more information about the physical problem, yet can be manipulated in much the same way as functions. Previously subspace collections had been introduced when there was an inner product on the vector (or Hilbert) space, and appropriate subspaces were mutually orthogonal. In that setting certain normalization and reduction operations on subspace collections led to a continued fraction expansion of the associated function, which allowed one to bound the function in terms of a set of weight matrices and normalization matrices that are derived from series expansions. Here we also initiate the theory of normalization and reduction operations, appropriate when there is no inner product on the space.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The response of linear inhomogeneous systems to coupled fields: Bounds and perturbation expansions

We consider the response of a multicomponent body to $n$ fields, such as electric fields, magnetic fields, temperature gradients, concentration gradients, etc., where each component, which is possibly anisotropic, may cross couple the various fields with different fluxes, such as electrical currents, electrical displacement currents, magnetic induction fields, energy fluxes, particle fluxes, etc. We obtain the form of the perturbation expansions of the fields and response tensor in powers of matrices which measure the difference between each component tensor and a homogeneous reference tensor ${\bf L}_0$. For the case of a statistically homogeneous or periodic composite the expansion coefficients can be expressed in terms of positive semidefinite normalization matrices alternating with positive semidefinite weight matrices, which at each given level sum to the identity matrix. In an appropriate basis the projection operators onto the relevant subspaces can be expressed in block tridiagonal form, where the blocks are functions of these weight and normalization matrices. This leads to continued fraction expansions for the effective tensor, and by truncating the continued fraction at successive levels one obtains a nested sequence of bounds on the effective tensor incorporating successively more weight and normalization matrices. The weight matrices and normalization matrices can be calculated from the series expansions of the fields which solve the conductivity problem alone, without any couplings to other fields, and then they can be used to obtain the solution for the fields and effective tensor in coupled field problems in composites.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Bounds for the response of viscoelastic composites under antiplane loadings in the time domain

To derive bounds on the strain and stress response of a two-component composite material with viscoelastic phases, we revisit the so-called analytic method (Bergman 1978), which allows one to approximate the complex effective tensor, function of the ratio of the component shear moduli, as the sum of poles weighted by positive semidefinite residue matrices. The novelty of this investigation lies in the application of such a method, previously applied (Milton 1980; Bergman 1980) to problems involving cyclic loadings in the frequency domain, to derive bounds in the time domain for the antiplane viscoelasticity case. The position of the poles and the residues matrices are the variational parameters of the problem: the aim is to determine such parameters in order to have the minimum (or maximum) response at any given moment in time. All the information about the composite is translated for each fixed pole configuration into constraints on the residues, the so-called sum rules. The linearity of the constraints, along with the observation that the response at a fixed time is linear in the residues, enables one to use the theory of linear programming to reduce the problem to one involving relatively few non-zero residues. Finally, bounds on the response are obtained by numerically optimizing over the pole positions. In the examples studied, the results turn out to be very accurate estimates: if sufficient information about the composite is available, the bounds can be quite tight over the entire range of time, allowing one to predict the transient behavior of the composite. Furthermore, the bounds incorporating the volume fractions (and possibly transverse isotropy) can be extremely tight at certain specific times.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Exact determination of the volume of an inclusion in a body having constant shear modulus

We derive an exact formula for the volume fraction of an inclusion in a body when the inclusion and the body are linearly elastic materials with the same shear modulus. Our formula depends on an appropriate measurement of the displacement and traction around the boundary of the body. In particular, the boundary conditions around the boundary of the body must be such that they mimic the body being placed in an infinite medium with an appropriate displacement applied at infinity.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Isotropic realizability of current fields in R^3

This paper deals with the isotropic realizability of a given regular divergence free field j in R^3 as a current field, namely to know when j can be written as sigma Du for some isotropic conductivity sigma, and some gradient field Du. The local isotropic realizability in R^3 is obtained by Frobenius' theorem provided that j and curl j are orthogonal in R^3. A counter-example shows that Frobenius' condition is not sufficient to derive the global isotropic realizability in R^3. However, assuming that (j, curl j, j x curl j) is an orthogonal basis of R^3, an admissible conductivity sigma is constructed from a combination of the three dynamical flows along the directions j/|j|, curl j/|curl j| and (j/|j|^2) x curl j. When the field j is periodic, the isotropic realizability in the torus needs in addition a boundedness assumption satisfied by the flow along the third direction (j/|j|^2) x \curl j. Several examples illustrate the sharpness of the realizability conditions.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Rigorous bounds on the effective moduli of composites and inhomogeneous bodies with negative-stiffness phases

We review the theoretical bounds on the effective properties of linear elastic inhomogeneous solids (including composite materials) in the presence of constituents having non-positive-definite elastic moduli (so-called negative-stiffness phases). We show that for statically stable bodies the classical displacement-based variational principles for Dirichlet and Neumann boundary problems hold but that the dual variational principle for traction boundary problems does not apply. We illustrate our findings by the example of a coated spherical inclusion whose stability conditions are obtained from the variational principles. We further show that the classical Voigt upper bound on the linear elastic moduli in multi-phase inhomogeneous bodies and composites applies and that it imposes a stability condition: overall stability requires that the effective moduli do not surpass the Voigt upper bound. This particularly implies that, while the geometric constraints among constituents in a composite can stabilize negative-stiffness phases, the stabilization is insufficient to allow for extreme overall static elastic moduli (exceeding those of the constituents). Stronger bounds on the effective elastic moduli of isotropic composites can be obtained from the Hashin-Shtrikman variational inequalities, which are also shown to hold in the presence of negative stiffness.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Sensitivity of anomalous localized resonance phenomena with respect to dissipation

We analyze cloaking due to anomalous localized resonance in the quasistatic regime in the case when a general charge density distribution is brought near a slab superlens. If the charge density distribution is within a critical distance of the slab, then the power dissipation within the slab blows up as certain electrical dissipation parameters go to zero. The potential remains bounded far away from the slab in this limit, which leads to cloaking due to anomalous localized resonance. On the other hand, if the charge density distribution is further than this critical distance from the slab, then the power dissipation within the slab remains bounded and cloaking due to anomalous localized resonance does not occur. The critical distance is shown to strongly depend on the the rate at which the dissipation outside of the slab goes to zero.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Sharp inequalities which generalize the divergence theorem--an extension of the notion of quasiconvexity, with an addendum

Subject to suitable boundary conditions being imposed, sharp inequalities are obtained on integrals over a region $Ω$ of certain special quadratic functions $f(\bf{E})$ where $\bf{E}(\bf{x})$ derives from a potential $\bf{U}(\bf{x})$. With $\bf{E}=\nabla\bf{U}$ it is known that such sharp inequalities can be obtained when $f(\bf{E})$ is a quasiconvex function and when $\bf{U}$ satisfies affine boundary conditions (i.e., for some matrix $\bf{D}$, $\bf{U}=\bf{D}\bf{x}$ on $\partialΩ$). Here we allow for other boundary conditions and for fields $\bf{E}$ that involve derivatives of a variety orders of $\bf{U}$. We define a notion of convexity that generalizes quasiconvexity. $Q^*$-convex quadratic functions are introduced, characterized and an algorithm is given for generating sharply $Q^*$-convex functions. We emphasize that this also solves the outstanding problem of finding an algorithm for generating extremal quasiconvex quadratic functions. We also treat integrals over $Ω$ of special quadratic functions $g(\bf{J})$ where $\bf{J}(\bf{x})$ satisfies a differential constraint involving derivatives with, possibly, a variety of orders. The results generalize an example of Kang and the author in three spatial dimensions where $\bf{J}(\bf{x})$ is a $3\times 3$ matrix valued field satisfying $\nabla\cdot\bf{J}=0$. In the addendum the paper is clarified. Notably, much more general boundary conditions are given under which sharp lower bounds on the integrals of certain quadratic functions of the fields can be obtained. As a consequence, strict $Q^*$-convexity is found to be an appropriate condition to ensure uniqueness of the solutions of a wide class of partial differential equations with constant coefficients, in a given domain $Ω$ with appropriate boundary conditions.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Anomalous localized resonance using a folded geometry in three dimensions

If a body of dielectric material is coated by a plasmonic structure of negative dielectric material with nonzero loss parameter, then cloaking by anomalous localized resonance (CALR) may occur as the loss parameter tends to zero. It was proved in other papers by authors that if the coated structure is circular (2D) and dielectric constant of the shell is a negative constant (with loss parameter), then CALR occurs, and if the coated structure is spherical (3D), then CALR does not occur. The aim of this paper is to show that the CALR takes place if the spherical coated structure has a specially designed anisotropic dielectric tensor. The anisotropic dielectric tensor is designed by unfolding a folded geometry.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Bounds on the size of an inclusion using the translation method for two-dimensional complex conductivity

The size estimation problem in electrical impedance tomography is considered when the conductivity is a complex number and the body is two-dimensional. Upper and lower bounds on the volume fraction of the unknown inclusion embedded in the body are derived in terms of two pairs of voltage and current data measured on the boundary of the body. These bounds are derived using the translation method. We also provide numerical examples to show that these bounds are quite tight and stable under measurement noise.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Bounds on the volume of an inclusion in a body from a complex conductivity measurement

We derive bounds on the volume of an inclusion in a body in two or three dimensions when the conductivities of the inclusion and the surrounding body are complex and assumed to be known. The bounds are derived in terms of average values of the electric field, current, and certain products of the electric field and current. All of these average values are computed from a single electrical impedance tomography measurement of the voltage and current on the boundary of the body. Additionally, the bounds are tight in the sense that at least one of the bounds gives the exact volume of the inclusion for certain geometries and boundary conditions.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Effective Conductivities of Thin-Interphase Composites

A method is presented for approximating the effective conductivity of composite media with thin interphase regions, which is exact to first order in the interphase thickness. The approximations are computationally efficient in the sense the fields need to be computed only in a reference composite in which the interphases have been replaced by perfect interfaces. The results apply whether any two phases of the composite are separated by a single interphase or multiple interphases, whether the conductivities of the composite phases are isotropic or anisotropic, and whether the thickness of an interphase is uniform or varies as a function of position. It is assumed that the conductivities of the interphase materials have intermediate values as opposed to very high or very low conductivities.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

On three-dimensional dilational elastic metamaterials

Dilational materials are stable three-dimensional isotropic auxetics with an ultimate Poisson's ratio of -1. We design, evaluate, fabricate, and characterize crystalline metamaterials approaching this ideal. To reveal all modes, we calculate the phonon band structures. On this basis, using cubic symmetry, we can unambiguously retrieve all different non-zero elements of the rank-4 effective metamaterial elasticity tensor, from which all effective elastic metamaterial properties follow. While the elastic properties and the phase velocity remain anisotropic, the effective Poisson's ratio indeed becomes isotropic and approaches -1 in the limit of small internal connections. This finding is also supported by independent static continuum-mechanics calculations. In static experiments on macroscopic polymer structures fabricated by three-dimensional printing, we measure Poisson's ratios as low as -0.8 in good agreement with theory. Microscopic samples are also presented.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The searchlight effect in hyperbolic materials

The quasistatic field around a circular hole in a two-dimensional hyperbolic medium is studied. As the loss parameter goes to zero, it is found that the electric field diverges along four lines each tangent to the hole. In this limit, the power dissipated by the field in the vicinity of these lines, per unit length of the line, goes to zero but extends further and further out so that the net power dissipated remains finite. Additionally the interaction between polarizable dipoles in a hyperbolic medium is studied. It is shown that a dipole with small polarizability can dramatically influence the dipole moment of a distant polarizable dipole, if it is appropriately placed. We call this the searchlight effect, as the enhancement depends on the orientation of the line joining the polarizable dipoles and can be varied by changing the frequency. For some particular polarizabilities the enhancement can actually increase the further the polarizable dipoles are apart.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Bounds on the volume fraction of the two-phase shallow shell using one measurement

We study the size estimate problem for the two-phase shallow shell equations in this paper. Our aim is to derive bounds on the volume fraction of each phase assuming that the material properties of the two phases are given. The approach in this paper is based on the translation method. One of the key steps is to connect the shallow shell equations to the thin plate equation.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Bounds on the volume fractions of two materials in a three dimensional body from boundary measurements by the translation method

Using the translation method of Tartar, Murat, Lurie, and Cherkaev bounds are derived on the volume occupied by an inclusion in a three-dimensional conducting body. They assume electrical impedance tomography measurements have been made for three sets of pairs of current flux and voltage measurements around the boundary. Additionally the conductivity of the inclusion and surrounding medium are assumed to be known. If the boundary data (Dirichlet or Neumann) is special, i.e. such that the fields inside the body would be uniform were the body homogeneous, then the bounds reduce to those of Milton and thus when the volume fraction is small to those of Capdeboscq and Vogelius.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Spectral theory of a Neumann-Poincaré-type operator and analysis of anomalous localized resonance II

If a body of dielectric material is coated by a plasmonic structure of negative dielectric constant with nonzero loss parameter, then cloaking by anomalous localized resonance (CALR) may occur as the loss parameter tends to zero. The aim of this paper is to investigate this phenomenon in two and three dimensions when the coated structure is radial, and the core, shell and matrix are isotropic materials. In two dimensions, we show that if the real part of the permittivity of the shell is -1 (under the assumption that the permittivity of the background is 1), then CALR takes place. If it is different from -1, then CALR does not occur. In three dimensions, we show that CALR does not occur. The analysis of this paper reveals that occurrence of CALR is determined by the eigenvalue distribution of the Neumann-Poincaré-type operator associated with the structure.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Spectral theory of a Neumann-Poincare-type operator and analysis of cloaking due to anomalous localized resonance

The aim of this paper is to give a mathematical justification of cloaking due to anomalous localized resonance (CALR). We consider the dielectric problem with a source term in a structure with a layer of plasmonic material. Using layer potentials and symmetrization techniques, we give a necessary and sufficient condition on the fixed source term for electromagnetic power dissipation to blow up as the loss parameter of the plasmonic material goes to zero. This condition is written in terms of the Newtonian potential of the source term. In the case of concentric disks, we make the condition even more explicit. Using the condition, we are able to show that for any source supported outside a critical radius CALR does not take place, and for sources located inside the critical radius satisfying certain conditions CALR does take place as the loss parameter goes to zero.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Bounds on strong field magneto-transport in three-dimensional composites

This paper deals with bounds satisfied by the effective non-symmetric conductivity of three-dimensional composites in the presence of a strong magnetic field. On the one hand, it is shown that for general composites the antisymmetric part of the effective conductivity cannot be bounded solely in terms of the antisymmetric part of the local conductivity, contrary to the columnar case. So, a suitable rank-two laminate the conductivity of which has a bounded antisymmetric part together with a high-contrast symmetric part, may generate an arbitrarily large antisymmetric part of the effective conductivity. On the other hand, bounds are provided which show that the antisymmetric part of the effective conductivity must go to zero if the upper bound on the antisymmetric part of the local conductivity goes to zero, and the symmetric part of the local conductivity remains bounded below and above. Elementary bounds on the effective moduli are derived assuming the local conductivity and effective conductivity have transverse isotropy in the plane orthogonal to the magnetic field. New Hashin-Shtrikman type bounds for two-phase three-dimensional composites with a non-symmetric conductivity are provided under geometric isotropy of the microstructure. The derivation of the bounds is based on a particular variational principle symmetrizing the problem, and the use of Y-tensors involving the averages of the fields in each phase.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Mathematical analysis of the two dimensional active exterior cloaking in the quasistatic regime

We design a device that generates fields canceling out a known probing field inside a region to be cloaked while generating very small fields far away from the device. The fields we consider satisfy the Laplace equation, but the approach remains valid in the quasistatic regime in a homogeneous medium. We start by relating the problem of designing an exterior cloak in the quasistatic regime to the classic problem of approximating a harmonic function with harmonic polynomials. An explicit polynomial solution to the problem was given earlier in [Phys. Rev. Lett. 103 (2009), 073901]. Here we show convergence of the device field to the field needed to perfectly cloak an object. The convergence region limits the size of the cloaked region, and the size and position of the device.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Minimization variational principles for acoustics, elastodynamics, and electromagnetism in lossy inhomogeneous bodies at fixed frequency

The classical energy minimization principles of Dirichlet and Thompson are extended as minimization principles to acoustics, elastodynamics and electromagnetism in lossy inhomogeneous bodies at fixed frequency. This is done by building upon ideas of Cherkaev and Gibiansky, who derived minimization variational principles for quasistatics. In the absence of free current the primary electromagnetic minimization variational principles have a minimum which is the time-averaged electrical power dissipated in the body. The variational principles provide constraints on the boundary values of the fields when the moduli are known. Conversely, when the boundary values of the fields have been measured, then they provide information about the values of the moduli within the body. This should have application to electromagnetic tomography. We also derive saddle point variational principles which correspond to variational principles of Gurtin, Willis, and Borcea.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Sharp bounds on the volume fractions of two materials in a two-dimensional body from electrical boundary measurements: the translation method

We deal with the problem of estimating the volume of inclusions using a finite number of boundary measurements in electrical impedance tomography. We derive upper and lower bounds on the volume fractions of inclusions, or more generally two phase mixtures, using two boundary measurements in two dimensions. These bounds are optimal in the sense that they are attained by certain configurations with some boundary data. We derive the bounds using the translation method which uses classical variational principles with a null Lagrangian. We then obtain necessary conditions for the bounds to be attained and prove that these bounds are attained by inclusions inside which the field is uniform. When special boundary conditions are imposed the bounds reduce to those obtained by Milton and these in turn are shown here to reduce to those of Capdeboscq-Vogelius in the limit when the volume fraction tends to zero. The bounds of this paper, and those of Milton, work for inclusions of arbitrary volume fractions. We then perform some numerical experiments to demonstrate how good these bounds are.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Spectral super-resolution in metamaterial composites

We investigate the optical properties of periodic composites containing metamaterial inclusions in a normal material matrix. We consider the case where these inclusions have sharp corners, and following Hetherington and Thorpe, use analytic results to argue that it is then possible to deduce the shape of the corner (its included angle) by measurements of the absorptance of such composites when the scale size of the inclusions and period cell is much finer than the wavelength. These analytic arguments are supported by highly accurate numerical results for the effective permittivity function of such composites as a function of the permittivity ratio of inclusions to matrix. The results show that this function has a continuous spectral component with limits independent of the area fraction of inclusions, and with the same limits for both square and staggered square arrays.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Transformation elastodynamics and active exterior acoustic cloaking

This chapter consists of three parts. In the first part we recall the elastodynamic equations under coordinate transformations. The idea is to use coordinate transformations to manipulate waves propagating in an elastic material. Then we study the effect of transformations on a mass-spring network model. The transformed networks can be realized with "torque springs", which are introduced here and are springs with a force proportional to the displacement in a direction other than the direction of the spring terminals. Possible homogenizations of the transformed networks are presented, with potential applications to cloaking. In the second and third parts we present cloaking methods that are based on cancelling an incident field using active devices which are exterior to the cloaked region and that do not generate significant fields far away from the devices. In the second part, the exterior cloaking problem for the Laplace equation is reformulated as the problem of polynomial approximation of analytic functions. An explicit solution is given that allows to cloak larger objects at a fixed distance from the cloaking device, compared to previous explicit solutions. In the third part we consider the active exterior cloaking problem for the Helmholtz equation in 3D. Our method uses the Green's formula and an addition theorem for spherical outgoing waves to design devices that mimic the effect of the single and double layer potentials in Green's formula.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Universal bounds on the electrical and elastic response of two-phase bodies and their application to bounding the volume fraction from boundary measurements

Universal bounds on the electrical and elastic response of two-phase (and multiphase) ellipsoidal or parallelopipedic bodies have been obtained by Nemat-Nasser and Hori. Here we show how their bounds can be improved and extended to bodies of arbitrary shape. Although our analysis is for two-phase bodies with isotropic phases it can easily be extended to multiphase bodies with anisotropic constituents. Our two-phase bounds can be used in an inverse fashion to bound the volume fractions occupied by the phases, and for electrical conductivity reduce to those of Capdeboscq and Vogelius when the volume fraction is asymptotically small. Other volume fraction bounds derived here utilize information obtained from thermal, magnetic, dielectric or elastic responses. One bound on the volume fraction can be obtained by simply immersing the body in a water filled cylinder with a piston at one end and measuring the change in water pressure when the piston is displaced by a known small amount. This bound may be particularly effective for estimating the volume of cavities in a body. We also obtain new bounds utilizing just one pair of (voltage, flux) electrical measurements at the boundary of the body.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Exterior cloaking with active sources in two dimensional acoustics

We cloak a region from a known incident wave by surrounding the region with three or more devices that cancel out the field in the cloaked region without significantly radiating waves. Since very little waves reach scatterers within the cloaked region, the scattered field is small and the scatterers are for all practical purposes undetectable. The devices are multipolar point sources that can be determined from Green's formula and an addition theorem for Hankel functions. The cloaking devices are exterior to the cloaked region.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Minimum variational principles for time-harmonic waves in a dissipative medium and associated variational principles of Hashin--Shtrikman type

Minimization variational principles for linear elastodynamic, acoustic, or electromagnetic time-harmonic waves in dissipative media were obtained by Milton, Seppecher and Bouchitté generalizing the quasistatic variational principles of Cherkaev and Gibiansky. Here a further generalization is made to allow for a much wider variety of boundary conditions, and in particular Dirichlet and Neumann boundary conditions. In addition minimization or maximization principles of the Hashin-Shtrikman type, incorporating "polarization fields", are developed.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Realizability of metamaterials with prescribed electric permittivity and magnetic permeability tensors

We show that any pair of real symmetric tensors $\BGve$ and $\BGm$ can be realized as the effective electric permittivity and effective magnetic permeability of a metamaterial at a given fixed frequency. The construction starts with two extremely low loss metamaterials, with arbitrarily small microstructure, whose existence is ensured by the work of Bouchitt{é} and Bourel and Bouchitté and Schweizer, one having at the given frequency a permittivity tensor with exactly one negative eigenvalue, and a positive permeability tensor, and the other having a positive permittivity tensor, and a permeability tensor having exactly one negative eigenvalue. To achieve the desired effective properties these materials are laminated together in a hierarchical multiple rank laminate structure, with widely separated length scales, and varying directions of lamination, but with the largest length scale still much shorter than the wavelengths and attenuation lengths in the macroscopic effective medium.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Active Exterior Cloaking

A new method of cloaking is presented. For two-dimensional quasistatics it is proven how a single active exterior cloaking device can be used to shield an object from surrounding fields, yet produce very small scattered fields. The problem is reduced to finding a polynomial which is approximately one within one disk and zero within a second disk, and such a polynomial is constructed. For the two-dimensional Helmholtz equation, it is numerically shown that three active exterior devices placed around the object suffice to produce very good cloaking.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Broadband Exterior Cloaking

It is shown how a recently proposed method of cloaking is effective over a broad range of frequencies. The method is based on three or more active devices. The devices, while not radiating significantly, create a ``quiet zone'' between the devices where the wave amplitude is small. Objects placed within this region are virtually invisible. The cloaking is demonstrated by simulations with a broadband incident pulse.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Complete characterization and synthesis of the response function of elastodynamic networks

The response function of a network of springs and masses, an elastodynamic network, is the matrix valued function $W(ω)$, depending on the frequency $ω$, mapping the displacements of some accessible or terminal nodes to the net forces at the terminals. We give necessary and sufficient conditions for a given function $W(ω)$ to be the response function of an elastodynamic network, assuming there is no damping. In particular we construct an elastodynamic network that can mimic a suitable response in the frequency or time domain. Our characterization is valid for networks in three dimensions and also for planar networks, which are networks where all the elements, displacements and forces are in a plane. The network we design can fit within an arbitrarily small neighborhood of the convex hull of the terminal nodes, provided the springs and masses occupy an arbitrarily small volume. Additionally, we prove stability of the network response to small changes in the spring constants and/or addition of springs with small spring constants.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Hybrid electromagnetic circuits

Electromagnetic circuits are the electromagnetic analog at fixed frequency of mass-spring networks in elastodynamics. By interchanging the roles of epsilon and mu in electromagnetic circuits one obtains magnetoelectric circuits. Here we show that by introducing tetrahedral connectors having epsilon=mu=0 one can join electromagnetic and magnetoelectric circuits to obtain hybrid circuits. Their response is governed by a symmetric matrix with negative semidefinite imaginary part. Conversely given any such matrix a recipe is given for constructing a hybrid circuit which has that matrix as its response matrix.

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

Optimizing the Superlens: manipulating geometry to enhance the resolution

We analyze the performance of a planar lens based on realistic negative index material in a generalized geometry. We demonstrate that the conventional superlens design (where the lens is centered between the object and the image) is not optimal from the resolution point-of-view, develop an analytical expression for the resolution limit of a generalized lens, use it to find the optimum lens configuration, and calculate the maximum absorption practical nearfield superlenses may have. We demonstrate that in contrast to the conventional superlens picture, planar imaging is typically accompanied by excitation of surface waves at both interfaces of the lens.

0 citations0 reviewsNo code linkedOpen access