Certified Rapid Solution of Parametrized Linear Elliptic Equations: Application to Parameter Estimation

We present a technique for the rapid and reliable evaluation of linear-functional output of elliptic partial differential equations with affine parameter dependence. The essential components are (i) rapidly uniformly convergent reduced-basis approximations — Galerkin projection onto a space WN spanned by solutions of the governing partial differential equation at N (optimally) selected points in parameter space; (ii) a posteriori error estimation — relaxations of the residual equation that provide inexpensive yet sharp and rigorous bounds for the error in the outputs; and (iii) offline/online computational procedures — stratagems that exploit affine parameter dependence to de-couple the generation and projection stages of the approximation process. The operation count for the online stage — in which, given a new parameter value, we calculate the output and associated error bound — depends only on N (typically small) and the parametric complexity of the problem. The method is thus ideally suited to the many-query and real-time contexts. In this paper, based on the technique we develop a robust inverse computational method for very fast solution of inverse problems characterized by parametrized partial differential equations. The essential ideas are in three-fold: first, we apply the technique to the forward problem for the rapid certified evaluation of PDE input-output relations and associated rigorous error bounds; second, we incorporate the reduced-basis approximation and error bounds into the inverse problem formulation; and third, rather than regularize the goodness-of-fit objective, we may instead identify all (or almost all, in the probabilistic sense) system configurations consistent with the available experimental data — well-posedness is reflected in a bounded "possibility region" that furthermore shrinks as the experimental error is decreased.

Optimization of inverse reflector design

Aquesta tesi presenta un nou mètode pel disseny invers de reflectors. Ens hem centrat en tres temes principals: l’ús de fonts de llum reals i complexes, la definició d’un algoritme ràpid pel càlcul de la il•luminació del reflector, i la definició d’un algoritme d’optimització per trobar més eficientment el reflector desitjat. Les fonts de llum estan representades per models near-field, que es comprimeixen amb un error molt petit, fins i tot per fonts de llum amb milions de raigs i objectes a il•luminar molt propers. Llavors proposem un mètode ràpid per obtenir la distribució de la il•luminació d’un reflector i la seva comparació amb la il•luminació desitjada, i que treballa completament en la GPU. Finalment, proposem un nou mètode d’optimització global que permet trobar la solució en menys passos que molts altres mètodes d’optimització clàssics, i alhora evitant mínims locals.

A survey on sampling and probe methods for inverse problems

The goal of the review is to provide a state-of-the-art survey on sampling and probe methods for the solution of inverse problems. Further, a configuration approach to some of the problems will be presented. We study the concepts and analytical results for several recent sampling and probe methods. We will give an introduction to the basic idea behind each method using a simple model problem and then provide some general formulation in terms of particular configurations to study the range of the arguments which are used to set up the method. This provides a novel way to present the algorithms and the analytic arguments for their investigation in a variety of different settings. In detail we investigate the probe method (Ikehata), linear sampling method (Colton-Kirsch) and the factorization method (Kirsch), singular sources Method (Potthast), no response test (Luke-Potthast), range test (Kusiak, Potthast and Sylvester) and the enclosure method (Ikehata) for the solution of inverse acoustic and electromagnetic scattering problems. The main ideas, approaches and convergence results of the methods are presented. For each method, we provide a historical survey about applications to different situations.

Unification of the probe and singular sources methods for the inverse boundary value problem by the no-response test

In this article, we use the no-response test idea, introduced in Luke and Potthast (2003) and Potthast (Preprint) and the inverse obstacle problem, to identify the interface of the discontinuity of the coefficient gamma of the equation del (.) gamma(x)del + c(x) with piecewise regular gamma and bounded function c(x). We use infinitely many Cauchy data as measurement and give a reconstructive method to localize the interface. We will base this multiwave version of the no-response test on two different proofs. The first one contains a pointwise estimate as used by the singular sources method. The second one is built on an energy (or an integral) estimate which is the basis of the probe method. As a conclusion of this, the probe and the singular sources methods are equivalent regarding their convergence and the no-response test can be seen as a unified framework for these methods. As a further contribution, we provide a formula to reconstruct the values of the jump of gamma(x), x is an element of partial derivative D at the boundary. A second consequence of this formula is that the blow-up rate of the indicator functions of the probe and singular sources methods at the interface is given by the order of the singularity of the fundamental solution.

Uniqueness results for direct and inverse scattering by infinite surfaces in a lossy medium

We consider the Dirichlet boundary-value problem for the Helmholtz equation, Au + x2u = 0, with Imx > 0. in an hrbitrary bounded or unbounded open set C c W. Assuming continuity of the solution up to the boundary and a bound on growth a infinity, that lu(x)l < Cexp (Slxl), for some C > 0 and S~< Imx, we prove that the homogeneous problem has only the trivial salution. With this resnlt we prove uniqueness results for direct and inverse problems of scattering by a bounded or infinite obstacle.

Experimental confirmation of transformation pathways between inverse double diamond and gyroid cubic phases

A macroscopically oriented double diamond inverse bicontinuous cubic phase (QIID) of the lipid glycerol monooleate is reversibly converted into a gyroid phase (QIIG). The initial QIID phase is prepared in the form of a film coating the inside of a capillary, deposited under flow, which produces a sample uniaxially oriented with a ⟨110⟩ axis parallel to the symmetry axis of the sample. A transformation is induced by replacing the water within the capillary tube with a solution of poly(ethylene glycol), which draws water out of the QIID sample by osmotic stress. This converts the QIID phase into a QIIG phase with two coexisting orientations, with the ⟨100⟩ and ⟨111⟩ axes parallel to the symmetry axis, as demonstrated by small-angle X-ray scattering. The process can then be reversed, to recover the initial orientation of QIID phase. The epitaxial relation between the two oriented mesophases is consistent with topologypreserving geometric pathways that have previously been hypothesized for the transformation. Furthermore, this has implications for the production of macroscopically oriented QIIG phases, in particular with applications as nanomaterial templates.

Data Assimilation and Inverse Methods in Terms of a Probabilistic Formulation

The weak-constraint inverse for nonlinear dynamical models is discussed and derived in terms of a probabilistic formulation. The well-known result that for Gaussian error statistics the minimum of the weak-constraint inverse is equal to the maximum-likelihood estimate is rederived. Then several methods based on ensemble statistics that can be used to find the smoother (as opposed to the filter) solution are introduced and compared to traditional methods. A strong point of the new methods is that they avoid the integration of adjoint equations, which is a complex task for real oceanographic or atmospheric applications. they also avoid iterative searches in a Hilbert space, and error estimates can be obtained without much additional computational effort. the feasibility of the new methods is illustrated in a two-layer quasigeostrophic model.

Inverse opal structure of SnO2 and SnO2 : Zn for gas sensing

In the present work, we propose a low cost synthetic sol-gel route that allows to produce high quality oxide nanostructures with inverse opal architecture which, transferred on alumina substrates provided with Pt interdigitated contacts and heater, are tested as gas sensing devices. An opal template of sintered monodisperse polystyrene spheres was filled with alcoholic solutions of metal oxide precursors and transferred on the alumina substrate. The polystyrene template was removed by thermal treatment, leading to the simultaneous sintering of the oxide nanoparticles. Beside SnO₂, a binary oxide well known for gas sensing application, a Zn containing ternary solid solution (SnO₂:Zn, with Zn 10% molar content) was taken into account for sensor preparation. The obtained high quality macro and meso-porous structures, characterized by different techniques, were tested for pollutant (CO, NO₂) and interfering (methanol) gases, showing that very good detection can be reached through the increase of surface area offered by the inverse opal structure and the tailoring of the chemical composition. The electrical characterization performed on the tin dioxide based sensors shows an enhancement of the relative response towards NO₂ at low temperatures in comparison with conventional SnO₂ sensors obtained with sputtering technique. The addition of Zn increases the separation between the operating temperatures for reducing and oxidizing gases and results in a further enhancement of the selectivity to NO₂ detection.

Inverse kinematics analysis of 6-RRCRR parallel manipulators

This paper presents a solution to the inverse kinematics of 6-RRCRR parallel manipulators with orthogonal non-intersecting RR-joint configuration. The inverse kinematics solution of such parallel robots compared with that of parallel robots with orthogonal intersecting RR-joint or universal joint configuration is more complex due to the existence of RR-joint variables. A novel methodology is established to define 6 independent variables of the actuators and 12 dependent RR-joint variables using the pose of the mobile platform with respect to the base frame. The constraint of RR-joints are analysed and the numerical algorithm to obtain joint variables is assessed. The forward kinematics of a 6- RRCRR parallel manipulator is modelled and computational analysis is performed in order to numerically verify the accuracy and effectiveness of the proposed methodology for the inverse kinematics analysis. Numerical results of a trajectory tracking simulation are provided. The results verify high accuracy for the proposed inverse kinematics solution of this special family of parallel micromanipulators.

Analysis of the inverse kinematics problem for 3-DOF axis-symmetric parallel manipulators with parasitic motion

 Determining an analytical solution to the inverse kinematics problem for a parallel manipulator is typically a straightforward problem. However, lower mobility parallel manipulators with 2-5 degrees of freedom (DOFs) often suffer from an unwanted parasitic motion in one or more DOFs. For such manipulators, the inverse kinematics problem can be significantly more difficult. This paper contains an analysis of the inverse kinematics problem for a class of 3-DOF parallel manipulators with axis-symmetric arm systems. All manipulators in the studied class exhibit parasitic motion in one DOF. For manipulators in the studied class, the general solution to the inverse kinematics problem is reduced to solving a univariate equation, while analytical solutions are presented for several important special cases.

Proposal of a Transmission Line Model Based on Lumped Elements: An Analytic Solution

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Laserlike emission from silica inverse opals infiltrated with rhodamine 6G

Monodisperse latex spheres were obtained by a surfactant free styrene polymerization method and used to obtain colloidal crystals by controlled centrifugation settling. Silica inverse opals were then prepared by using the colloidal crystals as templates and TEOS/ethanol solution. The inverse opals were infiltrated with Rhodamine 6G and laserlike emission was observed at 590 nm under 532 nm pump wavelength. The data show line narrowing of the dye fluorescence and a laser threshold of similar to 0.1 mJ/pulse. Local-field effects and light scattering due to structural defects are the main mechanisms contributing to generation of the laser-action observed. (c) 2005 Elsevier B.V. All rights reserved.

Rare earth doped synthetic opals and inverse opals

Monodisperse spheres of silica and latex were obtained by a surfactant free styrene polimerization and the Stober method respectively. Controlling settling either by centrifugation or by dip-coating colloidal crystals could be obtained. Silica inverse opals were prepared by using the latex colloidal crystals as templates and TEOS/ethanol solution. Eu3+ containing silica spheres were obtained dispersing silica spheres in Eu(NO3)(3) isopropanol solutions. Emission spectra suggest the formation of an amorphous Eu3+ containing phase well adhered at the spheres surface. The utilization of solutions of trifluoroacetates salts of Pb2+ and Eu3+ was observed to destroy the silica spherical pattern when samples are treated at 1000degreesC. In that case nanocrystals of PbF2 and amorphous silica were obtained after heat treatment.

Substitution-newton-Raphson method for the solution of electric network equations

The power flow problem, in transmission networks, has been well solved, for most cases, using Newton-Raphson method (NR) and its decoupled versions. Generally speaking, the solution of a non-linear system of equations refers to two methods: NR and Successive Substitution. The proposal of this paper is to evaluate the potential of the Substitution-Newton-Raphson Method (SNR), which combines both methods, on the solution of the power flow problem. Simulations were performed using a two-bus test network in order to observe the characteristics of these methods. It was verified that the NR is faster than SNR, in terms of convergence, considering non-stressed scenarios. For those cases where the power flow in the network is closed to the limits (stressed system), the SNR converges faster. This paper presents the power flow formulation of the SNR and describes its potential for its application in special cases such as stressed scenarios. © 2006 IEEE.

A hybrid method for the solution of a sparse power system matrices on vector computers

This paper describes a methodology for solving a linear system of equations on vector computer. The methodology combines direct and inverse factors. The decomposition and implementation of the direct solution in a CRAY Y-MPZE/232, and the performance results are discussed.