Extremal basic frequency of non-homogeneous plates

In this paper we propose two numerical algorithms to derive the extremal principal eigenvalue of the bi-Laplacian operator under Navier boundary conditions or Dirichlet boundary conditions. Consider a non-homogeneous hinged or clamped plate $\Omega$, the algorithms converge to the density functions on $\Omega$ which they yield the maximum or minimum basic frequency of the plate

Robust topology optimization of three-dimensional photonic-crystal band-gap structures

We perform full 3D topology optimization (in which "every voxel" of the unit cell is a degree of freedom) of photonic-crystal structures in order to find optimal omnidirectional band gaps for various symmetry groups, including fcc (including diamond), bcc, and simple-cubic lattices. Even without imposing the constraints of any fabrication process, the resulting optimal gaps are only slightly larger than previous hand designs, suggesting that current photonic crystals are nearly optimal in this respect. However, optimization can discover new structures, e.g. a new fcc structure with the same symmetry but slightly larger gap than the well known inverse opal, which may offer new degrees of freedom to future fabrication technologies. Furthermore, our band-gap optimization is an illustration of a computational approach to 3D dispersion engineering which is applicable to many other problems in optics, based on a novel semidefinite-program formulation for nonconvex eigenvalue optimization combined with other techniques such as a simple approach to impose symmetry constraints. We also demonstrate a technique for \emph{robust} topology optimization, in which some uncertainty is included in each voxel and we optimize the worst-case gap, and we show that the resulting band gaps have increased robustness to systematic fabrication errors.Comment: 17 pages, 9 figures, submitted to Optics Expres

Topological Derivative for Imaging of Thin Electromagnetic Inhomogeneity: Least Condition of Incident Directions

It is well-known that using topological derivative is an effective noniterative technique for imaging of crack-like electromagnetic inhomogeneity with small thickness when small number of incident directions are applied. However, there is no theoretical investigation about the configuration of the range of incident directions. In this paper, we carefully explore the mathematical structure of topological derivative imaging functional by establishing a relationship with an infinite series of Bessel functions of integer order of the first kind. Based on this, we identify the condition of the range of incident directions and it is highly depending on the shape of unknown defect. Results of numerical simulations with noisy data support our identification

On a continuation approach in Tikhonov regularization and its application in piecewise-constant parameter identification

We present a new approach to convexification of the Tikhonov regularization using a continuation method strategy. We embed the original minimization problem into a one-parameter family of minimization problems. Both the penalty term and the minimizer of the Tikhonov functional become dependent on a continuation parameter. In this way we can independently treat two main roles of the regularization term, which are stabilization of the ill-posed problem and introduction of the a priori knowledge. For zero continuation parameter we solve a relaxed regularization problem, which stabilizes the ill-posed problem in a weaker sense. The problem is recast to the original minimization by the continuation method and so the a priori knowledge is enforced. We apply this approach in the context of topology-to-shape geometry identification, where it allows to avoid the convergence of gradient-based methods to a local minima. We present illustrative results for magnetic induction tomography which is an example of PDE constrained inverse problem

Extremal Spectral Gaps for Periodic Schr\"odinger Operators

The spectrum of a Schr\"odinger operator with periodic potential generally consists of bands and gaps. In this paper, for fixed m, we consider the problem of maximizing the gap-to-midgap ratio for the m-th spectral gap over the class of potentials which have fixed periodicity and are pointwise bounded above and below. We prove that the potential maximizing the m-th gap-to-midgap ratio exists. In one dimension, we prove that the optimal potential attains the pointwise bounds almost everywhere in the domain and is a step-function attaining the imposed minimum and maximum values on exactly m intervals. Optimal potentials are computed numerically using a rearrangement algorithm and are observed to be periodic. In two dimensions, we develop an efficient rearrangement method for this problem based on a semi-definite formulation and apply it to study properties of extremal potentials. We show that, provided a geometric assumption about the maximizer holds, a lattice of disks maximizes the first gap-to-midgap ratio in the infinite contrast limit. Using an explicit parametrization of two-dimensional Bravais lattices, we also consider how the optimal value varies over all equal-volume lattices.Comment: 34 pages, 14 figure

Optimization of the electromagnetic scattering problem based on the topological derivative method

A new optimization method based on the topological derivative concept is developed for the electromagnetic design problem. Essentially, the purpose of the topological derivative method is to measure the sensitivity of a given shape functional with respect to a singular domain perturbation, so that it has applications in many relevant fields such as shape and topology optimization for imaging processing, inverse problems, and design of metamaterials. The topological derivative is rigorously derived for the electromagnetic scattering problem and used as gradient descent direction to find local optima for the design of electromagnetic devices. We demonstrate that the resulting topology design algorithm is remarkably simple and efficient and naturally leads to binary designs, while depending only on the solution of the conventional finite element formulation fir electrodynamics. Finally, several numerical experiments in two and three spatial dimensions are presented to illustrate the performance of the proposed formulation27233358633605CONSELHO NACIONAL DE DESENVOLVIMENTO CIENTÍFICO E TECNOLÓGICO - CNPQCOORDENAÇÃO DE APERFEIÇOAMENTO DE PESSOAL DE NÍVEL SUPERIOR - CAPESFUNDAÇÃO CARLOS CHAGAS FILHO DE AMPARO À PESQUISA DO ESTADO DO RIO DE JANEIRO - FAPERJFUNDAÇÃO DE AMPARO À PESQUISA DO ESTADO DE SÃO PAULO - FAPESP302036/2018-0; 310512/2017-4; 408274/2018-2; 438272/2018-888881.311020/2018E-26/203.041/20172015/24517-8; 2016/19270-6; 2018/25339-