BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Graeme W. Milton

Graeme W. Milton contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
math-phmath.MPphysics.class-phmath.APphysics.opticscond-mat.mtrl-scigr-qcmath.CVmath.NANumerical Analysis

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
14works
0followers
10topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

14 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