We consider a system governed by the wave equation with index of refraction @math , taken to be variable within a bounded region @math and constant in @math . The solution of the time-dependent wave equation with initial data, which is localized in @math , spreads and decays with advancing time. This rate of decay can be measured (for @math , and more generally, @math odd) in terms of the eigenvalues of the scattering resonance problem, a non--self-adjoint eigenvalue problem governing the time-harmonic solutions of the wave (Helmholtz) equation which are outgoing at @math . Specifically, the rate of energy escape from @math is governed by the complex scattering eigenfrequency, which is closest to the real axis. We study the structural design problem: Find a refractive index profile @math within an admissible class which has a scattering frequency with minimal imaginary part. The admissible class is defined in terms of the compact suppo...
Genetic algorithms have also been employed to minimize energy flux through the boundary @cite_3 @cite_24 . In @cite_1 @cite_32 @cite_43 the inverse method'' is employed, where a desired mode shape is chosen and then the material properties which produce that mode are found algebraically. In @cite_4 , the time-dependent problem is solved to steady state using a finite-difference method with perfectly matched layers to approximate the outgoing boundary conditions. The design problem is solved using a Nelder-Mead method.
@cite_4 @cite_1 @cite_32 @cite_3 @cite_24 @cite_43
2105969869 2004797988 2048738839 2098579833
We consider the inverse electrocardiographic problem of computing epicardial potentials from a body-surface potential map. We study how to improve numerical approximation of the inverse problem when the finite-element method is used. Being ill-posed, the inverse problem requires different discretization strategies from its corresponding forward problem. We propose refinement guidelines that specifically address the ill-posedness of the problem. The resulting guidelines necessitate the use of hybrid finite elements composed of tetrahedra and prism elements. Also, in order to maintain consistent numerical quality when the inverse problem is discretized into different scales, we propose a new family of regularizers using the variational principle underlying finite-element methods. These variational-formed regularizers serve as an alternative to the traditional Tikhonov regularizers, but preserves the L2 norm and thereby achieves consistent regularization in multiscale simulations. The variational formulation also enables a simple construction of the discrete gradient operator over irregular meshes, which is difficult to define in traditional discretization schemes. We validated our hybrid element technique and the variational regularizers by simulations on a realistic 3-D torso heart model with empirical heart data. Results show that discretization based on our proposed strategies mitigates the ill-conditioning and improves the inverse solution, and that the variational formulation may benefit a broader range of potential-based bioelectric problems. This paper describes the development of a novel metaheuristic that combines an electromagnetic-like mechanism (EM) and the great deluge algorithm (GD) for the University course timetabling problem. This well-known timetabling problem assigns lectures to specific numbers of timeslots and rooms maximizing the overall quality of the timetable while taking various constraints into account. EM is a population-based stochastic global optimization algorithm that is based on the theory of physics, simulating attraction and repulsion of sample points in moving toward optimality. GD is a local search procedure that allows worse solutions to be accepted based on some given upper boundary or level'. In this paper, the dynamic force calculated from the attraction-repulsion mechanism is used as a decreasing rate to update the level' within the search process. The proposed method has been applied to a range of benchmark university course timetabling test problems from the literature. Moreover, the viability of the method has been tested by comparing its results with other reported results from the literature, demonstrating that the method is able to produce improved solutions to those currently published. We believe this is due to the combination of both approaches and the ability of the resultant algorithm to converge all solutions at every search process. Consider a system governed by the time-dependent Schrodinger equation initialized in its ground state. When subjected to weak (size e) parametric forcing by an “ionizing field” (time-varying), the state decays with advancing time due to coupling of the bound state to radiation modes. The decay rate of this metastable state is governed by Fermi’s golden rule, Γ[V], which depends on the potential V and the details of the forcing. We pose the following potential design problem: find Vopt which minimizes Γ[V] (maximizes the lifetime of the state) over an admissible class of potentials with fixed spatial support. We formulate this problem as a constrained optimization problem and prove that an admissible optimal solution exists. Then, using quasi-Newton methods, we compute locally optimal potentials. These have the structure of a truncated periodic potential with a localized defect. In contrast to optimal structures for other spectral optimization problems, our optimizing potentials appear to be interior point... The aim of the present work is to validate a full vectorial electromagnetic inverse scattering algorithm against experimental data. Data were provided courtesy of Institut Fresnel, Marseille, France. These data were carried out in an anechoic chamber and correspond to different canonical targets as well as one mysterious object which is known only by experimentalists who measured the associated scattered field. The inverse algorithm was first developed in the optical domain and is adapted herein to the microwave domain. It is an iterative approach where the parameter of interest, namely the relative permittivity distribution, is updated gradually by minimizing a cost function describing the discrepancy between data and those that would be obtained via a forward solver for the best available estimate of the relative permittivity. The forward solver is based on the coupled dipole method which was introduced in the seventies to study the scattering of light by non-spherical dielectric grains. The forward and inverse schemes are briefly described and various examples are presented that demonstrate the efficiency of the inverse algorithm.